Skip to main content
  • Research article
  • Open access
  • Published:

Identification of microRNA profile specific to cancer stem-like cells directly isolated from human larynx cancer specimens

Abstract

Background

Emerging evidences proposed that microRNAs are associated with regulation of distinct physio-pathological processes including development of normal stem cells and carcinogenesis. In this study we aimed to investigate microRNA profile of cancer stem-like cells (CSLCs) isolated form freshly resected larynx cancer (LCa) tissue samples.

Methods

CD133 positive (CD133+) stem-like cells were isolated from freshly resected LCa tumor specimens. MicroRNA profile of 12 pair of CD133+ and CD133 cells was determined using microRNA microarray and differential expressions of selvected microRNAs were validated by quantitative real time PCR (qRT-PCR).

Results

MicroRNA profiling of CD133+ and CD133 LCa samples with microarray revealed that miR-26b, miR-203, miR-200c, and miR-363-3p were significantly downregulated and miR-1825 was upregulated in CD133+ larynx CSLCs. qRT-PCR analysis in a total of 25 CD133+/CD133 sample pairs confirmed the altered expressions of these five microRNAs. Expressions of miR-26b, miR-200c, and miR-203 were significantly correlated with miR-363-3p, miR-203, and miR-363-3p expressions, respectively. Furthermore, in silico analysis revealed that these microRNAs target both cancer and stem-cell associated signaling pathways.

Conclusions

Our results showed that certain microRNAs in CD133+ cells could be used as cancer stem cell markers. Based on these results, we propose that this panel of microRNAs might carry crucial roles in LCa pathogenesis through regulating stem cell properties of tumor cells.

Peer Review reports

Background

Larynx Cancer (LCa) is an aggressive neoplasm constituting approximately 1 to 2.5 % of all human cancer cases worldwide [13]. It is known as one of the most common tumor types of the head and neck region [4]. Despite notable enhancements in the therapeutic options; treatment outcome, prognosis, and 5-year survival rates for LCa remained almost unchanged in nearly past two decades [5, 6]. Therefore, more studies exploring the underlying mechanisms of LCa pathogenesis are urgently needed for better understanding of LCa development and providing more effective treatment strategies.

Emerging evidences propose the idea that a highly malignant rare subpopulation of tumor cells exhibits stem cell-like features [7]. This reservoir of stem-like cells within the bulk tumor is considered as tumor-initiating or cancer stem-like cells (CSLCs) with their unique capacity for unlimited self-renewal, multi-lineage differentiation, and ability for initiation, maintenance, and spread of tumor [8]. CSLCs has been proven to be present in a variety of tumors including lung, brain, breast, prostate, colon, ovarian, and head and neck cancers [9, 10] and are considered as the driving force for tumor relapse, metastasis, and chemo-radioresistance [1113].

We recently demonstrated that stem-like cells are highly enriched in CD133 overexpressing LCa cells, which are profoundly positive for stem cell markers including SOX2, OCT4, KLF4 and ABCG2 [14]. Furthermore, several studies have pointed to certain gene expression signatures specific to embryonic stem cells in acquisition and maintenance of the biological features of CSLCs [1517]; however, the underlying mechanisms are not yet completely understood. Therefore, elucidation of genetic and epigenetic circuits regulating the stem cell characteristics of CSLCs might help understanding the molecular basis of carcinogenesis.

There is an increasing body of evidence demonstrating that microRNAs (miRNAs) are associated with regulation of distinct physio-pathological processes including development of normal stem cells and carcinogenesis [18, 19]. MiRNAs are 21–25 nucleotides long, endogenously synthesized, noncoding RNAs that are involved in post-transcriptional gene silencing of target messenger RNAs (mRNAs) through binding 3′-untranslated regions (3′UTR) [20]. Deregulation of miRNAs has been linked to several diseases including cancer, where they can act as oncogenes or tumor suppressors. Recent studies implied miRNAs as crucial molecular players in cancer initiation, progression, and metastasis [2123].

Recently, utilization of Dicer or Dgcr8 knockout mice, lacking global miRNA processing capability, demonstrated that cells failed in self-renewal since stem cell specific markers couldn’t be downregulated. This indicated the significance of miRNAs in establishing stem cell identity [24, 25]. Besides, several miRNAs have been proposed to have direct roles in survival of CSLCs [8, 26, 27]. Therefore, understanding the contribution of miRNAs in acquisition and maintenance of CSLCs will provide the opportunity to develop miRNA-based therapeutic tools [28].

In this study, we investigated genome-wide miRNA expression profile of laryngeal CSLCs enriched for CD133 surface marker to identify a CSLCs specific miRNA signature.

Methods

Patients

This study has been reviewed and approved by an institutional review board of Istanbul University, Cerrahpasa Medical School (IRB No: 35697). 25 LCa tumor tissue specimens were obtained from Department of Otorhinolaryngology, Cerrahpasa Medical School, Istanbul University. None of the patients received radiotherapy, chemotherapy or immunotherapy subsequent to the surgery. The characteristics of the patients including age, gender, T classification and histological grade were summarized in Table 1. Freshly resected tumor tissues were collected immediately after the surgery and processed for CSLCs isolation. Patients were included into the study upon giving their written informed consents. We also obtained consent to publish from the participants.

Table 1 Clinico-pathological information of the patients

Cancer stem cell isolation

CD133 positive (CD133+) cells were isolated from freshly resected and physically/enzymatically dissociated tumor tissue samples using Magnetic-activated Cell Sorting (MACS) technique (Miltenyi Biotech, Bergisch Gladbach, Germany) and “EasySep Positive Selection Human PE Selection Kit (StemCell Technologies, (Vancouver, BC, Canada)” following the manufacturer’s protocol. Shortly, fresh tumor tissue samples were physically minced with a scalpel and exposed to enzymatic dissociation using 400 μg/ml Collagenase enzyme (GIBCO, New York, USA) at 37 °C for 3 h. Dissociated cells were filtered using a 70-μm cell strainer to get a single cell suspension. Cells were labeled with CD133/2-PE (Miltenyi Biotech clone AC133) antibody. After magnetic sorting, CD133 enriched (CD133+) and remaining (CD133) cell populations from the same tissue samples were immediately washed and homogenized in “Lysis/Binding Buffer” of “mirVana miRNA Isolation Kit” (Ambion, Darmstadt, Germany) for further RNA isolation.

RNA isolation

Total RNA was isolated from CD133+ and CD133 cells collected from LCa tumor samples using “mirVana miRNA Isolation Kit” (Ambion, Darmstadt, Germany) following the manufacturer’s instructions. The purities and concentrations of RNA samples were determined spectrophotometrically using NanoDrop ND-2000c (Thermo Fisher Scientific, Inc., Wilmington, DE).

MiRNA microarray and data analysis

Genome wide microRNA profiling of 12 pairs of CD133+and CD133 cell populations collected from 12 tumor samples were performed using Agilent Human miRNA Microarray (V19). 100 ng of total RNA from each sample were labeled with Cy3 by using Agilent miRNA labeling kit following manufacturer’s instructions. Labeled RNAs were heat denatured and hybridized to Agilent 8x15k miRNA microarray V19 comprised of 2006 miRNAs from Sanger miRBase (release 19) at 55 ° C for 20 h. After hybridization, slides were immediately washed and scanned in Agilent Microarray Scanner with Surescan High Resolution Technology (Agilent Technologies, Santa Clara, CA). Feature Extraction v10.7.3.1 (Agilent Technologies, CA) software was used to extract all features of the data obtained from the scanned images. Data were normalized by quantile normalization, using Bioconductor 2.10 with R version 2.15. Tumor samples were profiled on one of two different Agilent grid designs: Agilent-031181 (four pairs of CD133+ and CD133 cell populations collected from four tumor tissue samples) and Agilent-053955 (eight pairs of CD133+ and CD133 cell populations collected from eight tumor tissue samples); to correct for inter-platform differences, values were averaged by probe set, and each patient profile was compared with its corresponding CD133 profile by paired analysis (both pairs being represented on the same platform). P values and fold changes were calculated for each feature, using log-transformed values and paired t-test by patient. Differentially expressed miRNAs with P < 0.01 and 1.4-fold change were selected for further confirmation by RT-PCR. Array data have been deposited into the Gene Expression Omnibus (accession GSE69128).

MiRNA cDNA synthesis and quantitative reverse-transcription PCR

For the miRNA selection after microarray analysis, significantly deregulated miRNA probes were listed according to their fold changes. Then, top 10 upregulated and downregulated probes were selected for further literature search. We investigated the following properties and statuses for every single microRNA; deregulation in cancer, deregulation in larynx cancer, deregulation in head and neck cancers, expression in stem cells, and functional studies in stem cells. For top 10 upregulated microRNAs (hsa-miR-197-3p, hsa-miR-574-3p, hsa-miR-885-5p, hsa-miR-483-3p, hsa-miR-1281, hsa-miR-328, hsa-miR-4254, hsa-miR-4290, hsa-miR-1825, hsa-miR-766-3p), we included those have been shown to be deregulated in cancer, and have either expression data or functional studies in stem cells. Only hsa-miR-574-3p, hsa-miR-328, and hsa-miR-1825 met these criteria. For top 10 downregulated microRNAs (hsa-miR-106b-5p, hsa-miR-26b-5p, hsa-miR-494, hsa-miR-425-5p, hsa-miR-363-3p, hsa-miR-15b-5p, hsa-miR-185-5p, hsa-miR-150-5p, hsa-miR-223-3p, hsa-miR-142-5p), we included those have been shown to be deregulated in cancer (having no controversial expression status; some of these microRNAs have been shown to be upregulated in some cancer types, whereas, downregulated in other cancer types), and have either expression data or functional studies in stem cells. Only hsa-miR-26b-5p, hsa-miR-363-3p, and hsa-miR-223-3p met these criteria. Besides, we included miR-200c and miR-203 since they are strongly associated with stemness and cancer, although these miRNAs are not in the top 10 differentially expressed miRNAs.

To validate the differential expression of miR-26b, miR-200c, miR-203, miR-223, miR-328, miR-363-3p, 574-3p, and miR-1825, a total of 25 pairs of CD133+ and CD133 cell populations collected from 25 tumor samples including those used in microarray experiments were studied. First strand DNA (cDNA) synthesis was carried out with 30 ng of total RNA from each sample using miRNA specific primers purchased from Applied Biosystems and “TaqMan MicroRNA Reverse Transcription Kit” according to the manufacturer’s protocol (Applied Biosystems, Foster City, CA). MiRNA expression analysis by quantitative reverse-transcription PCR was carried out using a Roche LightCycler480-II real-time thermal cycler (Roche, Switzerland). TaqMan Universal Master Mix and TaqMan amplification kits (Applied Biosystems, Foster City, CA) were used. Expression levels of miRNAs in each CD133+ cell population were calculated as compared to CD133cells collected from the same tumor tissue sample. Therefore, expression levels of CD133 cells were fixed to 1 for every sample. RNU43 was used for normalization of miRNA expression analyses. Each experiment was performed in duplicate. The relative quantification analysis was performed by delta-delta-Ct method as described previously [29].

Statistical analysis

Data were plotted as mean ± standard error of the mean. Statistical analysis was carried out using two-sided Student’s t-test. Pearson Correlation test was used to show the correlation of differentially expressed mRNAs. A p-value < 0.05 was considered as statistically significant. MiRWalk 2.0 [30] and miRTarBase [31] in silico tools were used to estimate the predicted miRNA targets and to evaluate the validated miRNA targets, respectively. MiRWalk 2.0 is a freely accessible archive of predicted and experimentally verified miRNA-target interactions [30], whereas miRTarBase is a miRNA-target interactions database, where the collected miRNA-target interactions are validated experimentally by reporter assay, western blot, microarray and next-generation sequencing experiments [31]. In both tools, miRBase IDs were used as inputs. MiRWalk 2.0 and miRTarBase provide the gene list of predicted/validated targets of miRNAs and predicted gene interactors of miRNAs based on functional assays, respectively. The number of tumor suppressor and oncogenic targets of miRNAs were determined using miRWalk 2.0 tool. While predicting the targets of miRNAs, in ‘Step 4: Enriched functional patterns’, oncogene or tumor suppressor was selected as gene class, and the results provided the number of tumor suppressor and oncogenic targets of the specified miRNA. String [32] tool was utilized to prepare schematic representation of miRTarBase results. DIANA-miRPath was used for miRNA pathway analysis web-server [33].

Results

Subject characteristics

Twenty five LCa patients were included in this study to explore the miRNA expression status of CD133+ larynx CSLCs and remaining CD133 non-CSLCs. The diagnosis of patients has been confirmed histopathologically in Istanbul University, Cerrahpasa Medical School. All patients except one were men and their ages ranged from 44 to 84 years (median, 64 years). Histological grades of tumor specimens were determined according to World Health Organization classification, which demonstrated that 9 tumors were grade II and 16 tumors were grade III. Clinical characteristics of the patients are summarized in Table 1.

MiRNA profile of CD133+ larynx CSLCs

To analyze the global miRNA profile of CD133+ cells isolated from freshly resected LCa specimens, we performed microarray analysis using a discovery set comprised of 12 CD133+ and 12 CD133 samples. Microarray profiling revealed that 405 probes were differentially expressed with a p value <0.01 (paired t-test) and with at least 1.4-fold change. A heat map representation of the deregulated miRNAs is shown in Fig. 1 (entire set of differentially expressed miRNAs are provided as a Additional file 1: Data Set). Among those significantly differentially expressed miRNAs, five downregulated (miR-26b-5p, miR-200c-3p, miR-203a, miR-223-3p, miR-363-3p) and three upregulated (miR-328, miR-574-3p, miR-1825) miRNAs were selected as a result of detailed literature search for further confirmation with qRT-PCR.

Fig. 1
figure 1

Heat-map representation of top differentially expressed miRNAs (p < 0.01, paired t-test, fold change > 1.4) in CD133+ larynx CSLCs vs. CD133 larynx tumor cells, across twelve different patients. Yellow, high fold change in CD133+ patient sample as compared to its corresponding CD133 paired sample; blue, high fold change in CD133 sample compared to CD133+ paired sample

The qRT-PCR results confirmed that five of the eight selected miRNAs had a differential expression between groups: miR-26b, miR-200c, miR-203, miR-363-3p, and miR-1825 (Fig. 2, p values and fold changes are provided in Table 2). Among those, miR-26b (Fig. 2a, b), miR-200c (Fig. 2c, d), miR-203 (Fig. 2e, f), and miR-363-3p (Fig. 2g, h) were found to have significantly reduced expression in CD133+ larynx CSLCs, whereas miR-1825 (Fig. 2i, j) were validated to have increased expression in these CD133 enriched LCa cells. However, expression levels of miR-223-3p, miR-328, and miR-574-3p were not significantly different between CD133+ vs. CD133 LCa samples (Additional file 2: Figure S1, p values and fold changes are provided in Table 2). Although there was no statistically significant difference in the expression of miR-328 in CD133+ samples, its expression had a tendency to be elevated in CD133 enriched cell populations (Additional file 2: Figure S1C, D, Table 2). We further analyzed these miRNAs’ expressions with regard to T stage and histological stage of tumor samples. Results showed that miR-203 has lower expression in stage III samples compared to stage II tumor samples. Besides, miR-1825 has a tendency to have increased and miR-363-3p and miR-203 have a tendency to have decreased expression in T4 stage tumors compared to early stage tumors, although not significant. Since, miRNAs work in combination with each other rather than working individually and they operate in overlapping regulatory networks, demonstration of miRNAs’ correlation might be considered as indicative for their collaborative functioning in cells [34]. We, therefore, performed correlation analysis for the miRNAs found significantly deregulated between CD133+ and CD133 cells. To evaluate their correlation, we used Pearson correlation analysis, which demonstrated that miR-26b, miR-200c, and miR-203 expressions were significantly correlated with miR-363-3p, miR-203, and miR-363-3p expressions, respectively, in CD133+ LCa tissue samples (Fig. 3).

Fig. 2
figure 2

a Relative expression levels of miR-26b in each CD133+ and CD133 sample pairs, and (b) mean relative expression levels miR-26b in CD133+ cells with respect to CD133 cells. c Relative expression levels of miR-200c in each CD133+ and CD133 sample pairs, and (d) mean relative expression levels miR-200c in CD133+ cells with respect to CD133 cells. e Relative expression levels of miR-203 in each CD133+ and CD133 sample pairs, and (f) mean relative expression levels miR-203 in CD133+ cells with respect to CD133 cells. g Relative expression levels of miR-363-3p in each CD133+ and CD133 sample pairs, and (h) mean relative expression levels miR-363-3p in CD133+ cells with respect to CD133 cells. i Relative expression levels of miR-1825 in each CD133+ and CD133 sample pairs, and j mean relative expression levels miR-1825 in CD133+ cells with respect to CD133 cells

Table 2 Fold Changes and p values for miRNAs evaluated with qRT-PCR
Fig. 3
figure 3

Scatter plot representation of (a) miR-26b/miR-363-3p, (b) miR-200c/miR-203, (c) miR-203/miR-363-3p expression correlation in CD133+ samples

Relevant biological pathways affected from differentially expressed miRNAs

To explore the relevant biological pathways, which could be affected by the differential expression of miR-26b, miR-200c, miR-203, miR-363-3p, and miR-1825, we utilized DIANA miRPath v2.0, which revealed that several pathways overrepresented with a p-value <0.05, including cancer pathways (Table 3). Furthermore, miRWalk analysis showed that 3’UTR of several oncogenes were predicted to be targeted by miR-26b (202 out of 348 oncogenes), miR-200c (161 out of 348 oncogenes), miR-203 (216 out of 348 oncogenes), and miR-363-3p (175 out of 348 oncogenes). In addition, various tumor suppressors were estimated to be targeted by miR-1825 (29 out of 82 tumor suppressors). As to the analysis of validated targets of these miRNAs, miRTarBase database analysis revealed that miR-26b, miR-200c, miR-203, and miR-363-3p, and miR-1825 cooperatively target stem cell associated signaling pathways like Wnt, Hedgehog, and Notch (Fig. 4). Taken together, these analyses proposed that differential expression of miRNAs reported here might deregulate critical pathways involved in both carcinogenesis and establishment of CSLCs features.

Table 3 Overrepresented pathways, which could be affected by the differential expression of miR-26b, miR-200c, miR-203, miR-363-3p, and miR-1825
Fig. 4
figure 4

Validated targets of miR-26b, miR-200c, miR-203, and miR-363-3p, and miR-1825 are members of stem cell associated signaling pathways

Discussion

MiRNAs have emerged as an abundant class of small RNAs implicated in post-transcriptional gene regulation. Since every single miRNA can potentially target hundreds of genes, their cooperative and additive regulation have been postulated to have profound impacts on multiple pathways simultaneously [35]. In addition to their extensively studied roles in tumor biology, they were also proposed to participate in establishing stem cell associated features [36]. MiRNA-driven pathways were demonstrated to be fundamental for oncogenesis as well as acquisition and maintenance of CSLCs characteristics [37, 38]. However, currently, little is known about the miRNA expression profiles of CSLCs. Therefore, there is a need for comprehensive characterization of the miRNAs that might be involved in the acquisition and maintenance of stemness properties of CSLCs.

In this study, we enriched for CSLCs using CD133 surface marker, which are previously demonstrated to have increased potential for self-renewal and multi-lineage differentiating potency in vivo [39] and display elevated levels of stemness factors in LCa specimens [14]. We further investigated the miRNA profiles in CD133+ CSLCs to find out differences in miRNA expression that could distinguish them from their more differentiated progenies. Our findings identified a set of miRNAs in these cells, which might present valuable information for a better understanding of the molecular basis of carcinogenesis and regulation of cancer stem cell features.

Of the miRNAs investigated here, miR-1825 resides within 20q11.21 chromosomal region, where has been reported to be a recurrent gain of function abnormality in human embryonic stem cells and induced pluripotent stem cells [4042]. Nguyen et al. reported that human embryonic stem cells with 20q11.21 amplification displayed increased colony forming potential and decreased apoptosis [43]. Interestingly, 20q11.21 amplification in human embryonic stem cells resulted in acquisition of a gene-expression signature enriched for cancer-associated genes [44]. Recently, miR-1825 expression was reported to be elevated in majority of prostate cancer samples [45], and its expression was found to be upregulated in pancreatic cancer tissues in comparison to normal pancreatic duct [46]. In this study, we found overexpression of miR-1825 in CD133+ larynx CSLCs and suggest miR-1825 as an important contributor of carcinogenesis as a result of its dysregulation in CD133+ larynx CSLCs.

On the other hand, miR-363-3p, derived from the miR-106a-363 cluster on chromosome X, has been previously reported to be dysregulated in multiple cancers [47]. Although it acts different in distinct tumors, it has been demonstrated to behave like a tumor suppressor in several tumors such as colorectal cancer [47], nasal-type natural killer/T-cell lymphoma [48], head and neck squamous cell carcinoma [49], and hepatocellular carcinoma [50]. Interestingly, MYC, which has been implicated in stem cell self-renewal, maintenance of pluripotency, and control of cell fate decisions as well as carcinogenesis [51], was reported to directly bind to promoter of miR-363-3p and inhibit its expression [50]. MYC was also found to be destabilized by miR-363-3p through directly targeting and inhibiting USP28 [50] in hepatocellular carcinoma, pointing to a putative role for miR-363-3p in contribution to carcinogenesis and establishment of stemness features. Furthermore, miR-363-3p was found to directly target and repress GATA6, which is a transcription factor enhancing the expression of LGR5 in colorectal cancer [47]. LGR5, under the regulation of Wnt pathway, was proposed as a stem cell marker [52, 53] and its expression has been found to be overexpressed in various cancer tissues [54, 55]. These findings strengthen the potential role of miR-363-3p as a CSLCs specific miRNA in larynx pathogenesis. Additionally, miR-363 was reported to be repressed in head and neck squamous cell carcinoma tissues with lymph node metastasis and cell lines with increased invasive potential [49]. Ectopic expression of miR-363-3p decreased in vivo metastatic capacity of human neuroblastoma cells [56] and reversed the resistance of the breast cancer cell to the chemotherapeutic agent cisplatin [57]. Considering these findings and those of our own study, we suggest miR-363-3p as a strong candidate for establishment of stemness of CD133high CSLCs.

MiR-26b expression has been found to be downregulated in tongue [58], nasopharyngeal carcinoma [59], and oral cancers [60]. Exposure to cigarette smoke, as a major risk factor for LCa, was proposed to cause repression of miR-26 family members’ expressions in animal models [61]. Besides, recent findings indicated epigenetic silencing of miR-26a/b in cancer cells particularly through aberrant DNA hypermethylation [62, 63]. MiR-26b also displayed low levels of expression in a human embryonic stem cell line (HUES-17) and in a colorectal cancer cell line with a high metastatic potential (LoVo) [64]. Loss of miR-26b also enhanced migration and invasion in oral squamous cell carcinoma [65]. Additionally, miR-26b expression in neural stem cells was found to be induced during their differentiation into neurons in vivo [66]. MiR-26b was also reported to be overexpressed during osteogenic differentiation of unrestricted somatic stem cells, which comprise a rare subpopulation in human cord blood [67]. We, in the present study, showed that CD133+ LCa cells possess reduced miR-26b expression. Given that miR-26b is downregulated in both stem cells and cancer cells, our findings suggest miR-26b as a CSLC specific miRNA, whose deregulation might participate in oncogenic transformation and maintenance of stem cell state in larynx CSLCs and as well as others.

MiR-200c and miR-203 have previously been extensively studied in various contexts, and both miRNAs are associated with the stemness of normal stem cells and CSLCs. Loss of miR-200 expression was observed during conversion of immortalized human mammary epithelial cells to a stem-like phenotype [68]. CSLCs isolated from metastatic breast tumor tissues also exhibited reduced miR-200 expression [68]. Besides, miR-200c was found to be strongly downregulated in Oct3/4, Sox2, and Nanog overexpressing CD133+ ovarian cells [69]. MiR-200c expression was significantly reduced in ALDH1+/CD44+ cells with cancer stem cell potency in head and neck squamous cell carcinoma [70]. Additionally, miR-200c was reported to repress epithelial-to-mesenchymal transition through directly targeting and inhibiting ZEB1 and ZEB2 transcription factors [71]. As one of the central regulators of epithelial mesenchymal transition (EMT), abnormal expression of miR-200c was found to alter several important biological processes implied in cell–cell contact, cell adhesion and motility [72]. Interestingly, reduced miR-200c expression was significantly correlated with recurrence in LCa [73].

As to miR-203, it has been found to be downregulated in head and neck region cancers including LCa [60, 74, 75]. MiR-203 was found to induce differentiation of normal epidermal stem cells [76, 77] and its expression was reported to be inhibited during EMT in stem cell-enriched cancer cell subpopulation [78]. Lower miR-203 expression was significantly associated with poor differentiation, advanced clinical stages, lymph node metastasis and decreased 5-year overall survival in LCa [79]. Additionally, miR-200c and miR-203 cooperatively inhibit stem cell factors’ expressions in both cancer cells and mouse embryonic stem cells [80]. In this study, we have found miR-200c and miR-203 to be downregulated in CD133+ larynx CSLCs, supporting their potential involvement in carcinogenesis as driving forces for tumor initiation, progression, metastasis, and recurrence. We here also demonstrated that miR-200c, and miR-203 expressions were significantly correlated with miR-203, and miR-363-3p expressions, respectively, in CD133+ LCa tissue samples. Correlation of those microRNAs expressions supports a recent report, which demonstrated that decreased expressions of miR-200c, miR-363, and miR-203 were associated with poor prognosis in human head and neck squamous cell carcinoma [81]. Besides, coordinated loss of miR-200c and miR-203 has been found to result in enhanced translation of the multiple targets and chronic activation of NF-κB, PI3K-Akt, and Ras-Erk pathways, leading to B cell transformation, which suggest that collaborative actions of multiple miRNAs rather than a single miRNA ensure the robustness of biological processes [82].

Conclusion

In addition to miR-200c and miR-203, which have been demonstrated in distinct cancers as having CSLCs specific deregulation pattern, we propose miR-1825, miR-363-3p, and miR-26b as specific miRNAs with potential roles in acquisition and maintenance of stem cell associated features as well as in contributing to tumor initiation, progression, metastasis, chemoresistance, and recurrence. However, further detailed investigations are needed for each of the miRNAs studied here, to elucidate their roles in carcinogenesis and establishment of CSLCs related features.

Abbreviations

3′UTR:

3′-Untranslated regions

cDNA:

First strand DNA

CSLCs:

Cancer stem-like cells

LCa:

Larynx cancer

MACS:

Magnetic-activated cell sorting

miRNAs:

MicroRNAs

mRNAs:

Messenger RNAs

qRT-PCR:

Quantitative reverse transcription PCR

References

  1. Rudolph E, Dyckhoff G, Becher H, Dietz A, Ramroth H. Effects of tumour stage, comorbidity and therapy on survival of laryngeal cancer patients: a systematic review and a meta-analysis. Eur Arch Otorhinolaryngol. 2011;268:165–79.

    Article  Google Scholar 

  2. Farhadieh RD, Rees CG, Yang JL, Salardini A, Russell P, Smee R. Radiotherapy in larynx squamous cell carcinoma is not associated with an increased diagnosis of second primary tumours. Clin Oncol (R Coll Radiol). 2009;21:315–9.

    Article  CAS  Google Scholar 

  3. Yilmaz M, Karatas OF, Yuceturk B, Dag H, Yener M, Ozen M. Alpha-B-crystallin expression in human laryngeal squamous cell carcinoma tissues. Head Neck. 2015;37(9):1344–8.

    Article  Google Scholar 

  4. Chu EA, Kim YJ. Laryngeal cancer: diagnosis and preoperative work-up. Otolaryngol Clin North Am. 2008;41:673–95. v.

    Article  Google Scholar 

  5. Rothman KJ, Cann CI, Flanders D, Fried MP. Epidemiology of laryngeal cancer. Epidemiol Rev. 1980;2:195–209.

    Article  CAS  Google Scholar 

  6. Siegel R, Naishadham D, Jemal A. Cancer statistics, 2013. CA Cancer J Clin. 2013;63:11–30.

    Article  Google Scholar 

  7. Yu F, Yao H, Zhu P, Zhang X, Pan Q, Gong C, Huang Y, Hu X, Su F, Lieberman J, Song E. let-7 regulates self renewal and tumorigenicity of breast cancer cells. Cell. 2007;131:1109–23.

    Article  CAS  Google Scholar 

  8. Qian S, Ding JY, Xie R, An JH, Ao XJ, Zhao ZG, Sun JG, Duan YZ, Chen ZT, Zhu B. MicroRNA expression profile of bronchioalveolar stem cells from mouse lung. Biochem Biophys Res Commun. 2008;377:668–73.

    Article  CAS  Google Scholar 

  9. Hao J, Zhao S, Zhang Y, Zhao Z, Ye R, Wen J, Li J. Emerging role of MicroRNAs in cancer and cancer stem cells. J Cell Biochem. 2014;115:605–10.

    Article  CAS  Google Scholar 

  10. Leal JA, Lleonart ME. MicroRNAs and cancer stem cells: therapeutic approaches and future perspectives. Cancer Lett. 2013;338:174–83.

    Article  CAS  Google Scholar 

  11. Al-Hajj M, Wicha MS, Benito-Hernandez A, Morrison SJ, Clarke MF. Prospective identification of tumorigenic breast cancer cells. Proc Natl Acad Sci U S A. 2003;100:3983–8.

    Article  CAS  Google Scholar 

  12. Huang EH, Wicha MS. Colon cancer stem cells: implications for prevention and therapy. Trends Mol Med. 2008;14:503–9.

    Article  CAS  Google Scholar 

  13. Tan BT, Park CY, Ailles LE, Weissman IL. The cancer stem cell hypothesis: a work in progress. Lab Invest. 2006;86:1203–7.

    Article  CAS  Google Scholar 

  14. Suer I, Karatas OF, Yuceturk B, Yilmaz M, Guven G, Oz B, Cansiz H, Ozen M. Characterization of stem-like cells directly isolated from freshly resected laryngeal squamous cell carcinoma specimens. Curr Stem Cell Res Ther. 2014;9:347–53.

    Article  CAS  Google Scholar 

  15. Chen YC, Hsu HS, Chen YW, Tsai TH, How CK, Wang CY, Hung SC, Chang YL, Tsai ML, Lee YY, et al. Oct-4 expression maintained cancer stem-like properties in lung cancer-derived CD133-positive cells. PLoS One. 2008;3:e2637.

    Article  Google Scholar 

  16. Gangemi RM, Griffero F, Marubbi D, Perera M, Capra MC, Malatesta P, Ravetti GL, Zona GL, Daga A, Corte G. SOX2 silencing in glioblastoma tumor-initiating cells causes stop of proliferation and loss of tumorigenicity. Stem Cells. 2009;27:40–8.

    Article  CAS  Google Scholar 

  17. Ben-Porath I, Thomson MW, Carey VJ, Ge R, Bell GW, Regev A, Weinberg RA. An embryonic stem cell-like gene expression signature in poorly differentiated aggressive human tumors. Nat Genet. 2008;40:499–507.

    Article  CAS  Google Scholar 

  18. Croce CM, Calin GA. miRNAs, cancer, and stem cell division. Cell. 2005;122:6–7.

    Article  CAS  Google Scholar 

  19. Caldas C, Brenton JD. Sizing up miRNAs as cancer genes. Nat Med. 2005;11:712–4.

    Article  CAS  Google Scholar 

  20. Iorio MV, Croce CM. MicroRNA dysregulation in cancer: diagnostics, monitoring and therapeutics. A comprehensive review. EMBO Mol Med. 2012;4:143–59.

    Article  CAS  Google Scholar 

  21. Baranwal S, Alahari SK. miRNA control of tumor cell invasion and metastasis. Int J Cancer. 2010;126:1283–90.

    CAS  Google Scholar 

  22. Karatas OF, Guzel E, Suer I, Ekici ID, Caskurlu T, Creighton CJ, Ittmann M, Ozen M. miR-1 and miR-133b Are differentially expressed in patients with recurrent prostate cancer. Plos One. 2014;9:e98675.

    Article  Google Scholar 

  23. Karatas OF, Yuceturk B, Suer I, Yilmaz M, Cansiz H, Ittmann M, Ozen M. The role of miR-145 in Human Laryngeal Squamous Cell Carcinoma. Head Neck. 2016;38(2):260–6.

    Article  Google Scholar 

  24. Yu Z, Pestell TG, Lisanti MP, Pestell RG. Cancer stem cells. Int J Biochem Cell Biol. 2012;44:2144–51.

    Article  CAS  Google Scholar 

  25. Martinez NJ, Gregory RI. MicroRNA gene regulatory pathways in the establishment and maintenance of ESC identity. Cell Stem Cell. 2010;7:31–5.

    Article  CAS  Google Scholar 

  26. Liu C, Kelnar K, Liu B, Chen X, Calhoun-Davis T, Li H, Patrawala L, Yan H, Jeter C, Honorio S, et al. The microRNA miR-34a inhibits prostate cancer stem cells and metastasis by directly repressing CD44. Nat Med. 2011;17:211–5.

    Article  CAS  Google Scholar 

  27. Wang L, Zhang D, Zhang C, Zhang S, Wang Z, Qu C, Liu S. A microRNA expression signature characterizing the properties of tumor-initiating cells for breast cancer. Oncol Lett. 2012;3:119–24.

    CAS  Google Scholar 

  28. Hao J, Zhao S, Zhang Y, Zhao Z, Ye R, Wen J, Li J: Emerging Role of MicroRNAs in Cancer and Cancer Stem Cells. J Cell Biochem. 2013.

  29. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C (T)) method. Methods. 2001;25:402–8.

    Article  CAS  Google Scholar 

  30. Dweep H, Sticht C, Pandey P, Gretz N. miRWalk--database: prediction of possible miRNA binding sites by “walking” the genes of three genomes. J Biomed Inform. 2011;44:839–47.

    Article  CAS  Google Scholar 

  31. Hsu SD, Tseng YT, Shrestha S, Lin YL, Khaleel A, Chou CH, Chu CF, Huang HY, Lin CM, Ho SY, et al. miRTarBase update 2014: an information resource for experimentally validated miRNA-target interactions. Nucleic Acids Res. 2014;42:D78–85.

    Article  CAS  Google Scholar 

  32. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52.

    Article  Google Scholar 

  33. Vlachos IS, Kostoulas N, Vergoulis T, Georgakilas G, Reczko M, Maragkakis M, Paraskevopoulou MD, Prionidis K, Dalamagas T, Hatzigeorgiou AG. DIANA miRPath v. 2.0: investigating the combinatorial effect of microRNAs in pathways. Nucleic Acids Res. 2012;40:W498–504.

    Article  CAS  Google Scholar 

  34. Mallory AC, Vaucheret H. Functions of microRNAs and related small RNAs in plants. Nat Genet. 2006;38(Suppl):S31–6.

    Article  CAS  Google Scholar 

  35. Liu S, Yin F, Zhang J, Wicha MS, Chang AE, Fan W, Chen L, Fan M, Li Q. Regulatory roles of miRNA in the human neural stem cell transformation to glioma stem cells. J Cell Biochem. 2014;115:1368–80.

    Article  CAS  Google Scholar 

  36. Sanchez-Diaz PC, Hsiao TH, Chang JC, Yue D, Tan MC, Chen HI, Tomlinson GE, Huang Y, Chen Y, Hung JY. De-regulated microRNAs in pediatric cancer stem cells target pathways involved in cell proliferation, cell cycle and development. PLoS One. 2013;8:e61622.

    Article  CAS  Google Scholar 

  37. Nam EJ, Lee M, Yim GW, Kim JH, Kim S, Kim SW, Kim YT. MicroRNA profiling of a CD133(+) spheroid-forming subpopulation of the OVCAR3 human ovarian cancer cell line. BMC Med Genomics. 2012;5:18.

    Article  CAS  Google Scholar 

  38. Mishra L, Banker T, Murray J, Byers S, Thenappan A, He AR, Shetty K, Johnson L, Reddy EP. Liver stem cells and hepatocellular carcinoma. Hepatology. 2009;49:318–29.

    Article  Google Scholar 

  39. 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:5821–8.

    CAS  Google Scholar 

  40. Lund RJ, Närvä E, Lahesmaa R. Genetic and epigenetic stability of human pluripotent stem cells. Nat Rev Genet. 2012;13:732–44.

    Article  CAS  Google Scholar 

  41. Nguyen HT, Geens M, Spits C. Genetic and epigenetic instability in human pluripotent stem cells. Hum Reprod Update. 2013;19:187–205.

    Article  CAS  Google Scholar 

  42. Martins-Taylor K, Nisler BS, Taapken SM, Compton T, Crandall L, Montgomery KD, Lalande M, Xu RH. Recurrent copy number variations in human induced pluripotent stem cells. Nat Biotechnol. 2011;29:488–91.

    Article  CAS  Google Scholar 

  43. Nguyen HT, Geens M, Mertzanidou A, Jacobs K, Heirman C, Breckpot K, Spits C. Gain of 20q11.21 in human embryonic stem cells improves cell survival by increased expression of Bcl-xL. Mol Hum Reprod. 2014;20:168–77.

    Article  CAS  Google Scholar 

  44. Werbowetski-Ogilvie TE, Bossé M, Stewart M, Schnerch A, Ramos-Mejia V, Rouleau A, Wynder T, Smith MJ, Dingwall S, Carter T, et al. Characterization of human embryonic stem cells with features of neoplastic progression. Nat Biotechnol. 2009;27:91–7.

    Article  CAS  Google Scholar 

  45. Haj-Ahmad TA, Abdalla MA, Haj-Ahmad Y. Potential urinary miRNA biomarker candidates for the accurate detection of prostate cancer among benign prostatic hyperplasia patients. J Cancer. 2014;5:182–91.

    Article  CAS  Google Scholar 

  46. Li A, Yu J, Kim H, Wolfgang CL, Canto MI, Hruban RH, Goggins M. MicroRNA array analysis finds elevated serum miR-1290 accurately distinguishes patients with low-stage pancreatic cancer from healthy and disease controls. Clin Cancer Res. 2013;19:3600–10.

    Article  CAS  Google Scholar 

  47. Tsuji S, Kawasaki Y, Furukawa S, Taniue K, Hayashi T, Okuno M, Hiyoshi M, Kitayama J, Akiyama T. The miR-363-GATA6-Lgr5 pathway is critical for colorectal tumourigenesis. Nat Commun. 2014;5:3150.

    Google Scholar 

  48. Ng SB, Yan J, Huang G, Selvarajan V, Tay JL, Lin B, Bi C, Tan J, Kwong YL, Shimizu N, et al. Dysregulated microRNAs affect pathways and targets of biologic relevance in nasal-type natural killer/T-cell lymphoma. Blood. 2011;118:4919–29.

    Article  CAS  Google Scholar 

  49. Sun Q, Zhang J, Cao W, Wang X, Xu Q, Yan M, Wu X, Chen W. Dysregulated miR-363 affects head and neck cancer invasion and metastasis by targeting podoplanin. Int J Biochem Cell Biol. 2013;45:513–20.

    Article  CAS  Google Scholar 

  50. Han H, Sun D, Li W, Shen H, Zhu Y, Li C, Chen Y, Lu L, Zhang J, Tian Y, Li Y. A c-Myc-MicroRNA functional feedback loop affects hepatocarcinogenesis. Hepatology. 2013;57:2378–89.

    Article  CAS  Google Scholar 

  51. Chappell J, Dalton S. Roles for MYC in the establishment and maintenance of pluripotency. Cold Spring Harb Perspect Med. 2013;3:a014381.

    Article  Google Scholar 

  52. Barker N, van Es JH, Kuipers J, Kujala P, van den Born M, Cozijnsen M, Haegebarth A, Korving J, Begthel H, Peters PJ, Clevers H. Identification of stem cells in small intestine and colon by marker gene Lgr5. Nature. 2007;449:1003–7.

    Article  CAS  Google Scholar 

  53. Tian H, Biehs B, Warming S, Leong KG, Rangell L, Klein OD, de Sauvage FJ. A reserve stem cell population in small intestine renders Lgr5-positive cells dispensable. Nature. 2011;478:255–9.

    Article  CAS  Google Scholar 

  54. Tanese K, Fukuma M, Yamada T, Mori T, Yoshikawa T, Watanabe W, Ishiko A, Amagai M, Nishikawa T, Sakamoto M. G-protein-coupled receptor GPR49 is up-regulated in basal cell carcinoma and promotes cell proliferation and tumor formation. Am J Pathol. 2008;173:835–43.

    Article  CAS  Google Scholar 

  55. McClanahan T, Koseoglu S, Smith K, Grein J, Gustafson E, Black S, Kirschmeier P, Samatar AA. Identification of overexpression of orphan G protein-coupled receptor GPR49 in human colon and ovarian primary tumors. Cancer Biol Ther. 2006;5:419–26.

    Article  CAS  Google Scholar 

  56. Qiao J, Lee S, Paul P, Theiss L, Tiao J, Qiao L, Kong A, Chung DH. miR-335 and miR-363 regulation of neuroblastoma tumorigenesis and metastasis. Surgery. 2013;154:226–33.

    Article  Google Scholar 

  57. Zhang R, Li Y, Dong X, Peng L, Nie X. MiR-363 sensitizes cisplatin-induced apoptosis targeting in Mcl-1 in breast cancer. Med Oncol. 2014;31:347.

    Article  CAS  Google Scholar 

  58. Wong TS, Liu XB, Wong BY, Ng RW, Yuen AP, Wei WI. Mature miR-184 as potential oncogenic microRNA of squamous cell carcinoma of tongue. Clin Cancer Res. 2008;14:2588–92.

    Article  CAS  Google Scholar 

  59. Lu J, He ML, Wang L, Chen Y, Liu X, Dong Q, Chen YC, Peng Y, Yao KT, Kung HF, Li XP. MiR-26a inhibits cell growth and tumorigenesis of nasopharyngeal carcinoma through repression of EZH2. Cancer Res. 2011;71:225–33.

    Article  CAS  Google Scholar 

  60. Kozaki K, Imoto I, Mogi S, Omura K, Inazawa J. Exploration of tumor-suppressive microRNAs silenced by DNA hypermethylation in oral cancer. Cancer Res. 2008;68:2094–105.

    Article  CAS  Google Scholar 

  61. Izzotti A, Calin GA, Arrigo P, Steele VE, Croce CM, De Flora S. Downregulation of microRNA expression in the lungs of rats exposed to cigarette smoke. FASEB J. 2009;23:806–12.

    Article  CAS  Google Scholar 

  62. Sandhu R, Rivenbark AG, Coleman WB. Loss of post-transcriptional regulation of DNMT3b by microRNAs: a possible molecular mechanism for the hypermethylation defect observed in a subset of breast cancer cell lines. Int J Oncol. 2012;41:721–32.

    CAS  Google Scholar 

  63. Börno ST, Fischer A, Kerick M, Fälth M, Laible M, Brase JC, Kuner R, Dahl A, Grimm C, Sayanjali B, et al. Genome-wide DNA methylation events in TMPRSS2-ERG fusion-negative prostate cancers implicate an EZH2-dependent mechanism with miR-26a hypermethylation. Cancer Discov. 2012;2:1024–35.

    Article  Google Scholar 

  64. Ma YL, Zhang P, Wang F, Moyer MP, Yang JJ, Liu ZH, Peng JY, Chen HQ, Zhou YK, Liu WJ, Qin HL. Human embryonic stem cells and metastatic colorectal cancer cells shared the common endogenous human microRNA-26b. J Cell Mol Med. 2011;15:1941–54.

    Article  CAS  Google Scholar 

  65. Fukumoto I, Hanazawa T, Kinoshita T, Kikkawa N, Koshizuka K, Goto Y, Nishikawa R, Chiyomaru T, Enokida H, Nakagawa M, et al. MicroRNA expression signature of oral squamous cell carcinoma: functional role of microRNA-26a/b in the modulation of novel cancer pathways. Br J Cancer. 2015;112:891–900.

    Article  CAS  Google Scholar 

  66. Dill H, Linder B, Fehr A, Fischer U. Intronic miR-26b controls neuronal differentiation by repressing its host transcript, ctdsp2. Genes Dev. 2012;26:25–30.

    Article  CAS  Google Scholar 

  67. Trompeter HI, Dreesen J, Hermann E, Iwaniuk KM, Hafner M, Renwick N, Tuschl T, Wernet P. MicroRNAs miR-26a, miR-26b, and miR-29b accelerate osteogenic differentiation of unrestricted somatic stem cells from human cord blood. BMC Genomics. 2013;14:111.

    Article  CAS  Google Scholar 

  68. Lim YY, Wright JA, Attema JL, Gregory PA, Bert AG, Smith E, Thomas D, Lopez AF, Drew PA, Khew-Goodall Y, Goodall GJ. Epigenetic modulation of the miR-200 family is associated with transition to a breast cancer stem-cell-like state. J Cell Sci. 2013;126:2256–66.

    Article  CAS  Google Scholar 

  69. Guo R, Wu Q, Liu F, Wang Y. Description of the CD133+ subpopulation of the human ovarian cancer cell line OVCAR3. Oncol Rep. 2011;25:141–6.

    CAS  Google Scholar 

  70. Tu HF, Lin SC, Chang KW. MicroRNA aberrances in head and neck cancer: pathogenetic and clinical significance. Curr Opin Otolaryngol Head Neck Surg. 2013;21:104–11.

    Article  Google Scholar 

  71. Korpal M, Lee ES, Hu G, Kang Y. The miR-200 family inhibits epithelial-mesenchymal transition and cancer cell migration by direct targeting of E-cadherin transcriptional repressors ZEB1 and ZEB2. J Biol Chem. 2008;283:14910–4.

    Article  CAS  Google Scholar 

  72. Knezevic J, Pfefferle AD, Petrovic I, Greene SB, Perou CM, Rosen JM. Expression of miR-200c in claudin-low breast cancer alters stem cell functionality, enhances chemosensitivity and reduces metastatic potential. Oncogene. 2015;34:5997–6006.

    Article  CAS  Google Scholar 

  73. Cappellesso R, Marioni G, Crescenzi M, Giacomelli L, Guzzardo V, Mussato A, Staffieri A, Martini A, Blandamura S, Fassina A. The prognostic role of the epithelial-mesenchymal transition markers E-cadherin and slug in laryngeal squamous cell carcinoma. Histopathology. 2015;67:491–500.

    Article  Google Scholar 

  74. Mathé EA, Nguyen GH, Bowman ED, Zhao Y, Budhu A, Schetter AJ, Braun R, Reimers M, Kumamoto K, Hughes D, et al. MicroRNA expression in squamous cell carcinoma and adenocarcinoma of the esophagus: associations with survival. Clin Cancer Res. 2009;15:6192–200.

    Article  Google Scholar 

  75. Bian K, Fan J, Zhang X, Yang XW, Zhu HY, Wang L, Sun JY, Meng YL, Cui PC, Cheng SY, et al. MicroRNA-203 leads to G1 phase cell cycle arrest in laryngeal carcinoma cells by directly targeting survivin. FEBS Lett. 2012;586:804–9.

    Article  CAS  Google Scholar 

  76. Jackson SJ, Zhang Z, Feng D, Flagg M, O’Loughlin E, Wang D, Stokes N, Fuchs E, Yi R. Rapid and widespread suppression of self-renewal by microRNA-203 during epidermal differentiation. Development. 2013;140:1882–91.

    Article  CAS  Google Scholar 

  77. Yi R, Poy MN, Stoffel M, Fuchs E. A skin microRNA promotes differentiation by repressing ‘stemness’. Nature. 2008;452:225–9.

    Article  CAS  Google Scholar 

  78. Taube JH, Malouf GG, Lu E, Sphyris N, Vijay V, Ramachandran PP, Ueno KR, Gaur S, Nicoloso MS, Rossi S, et al. Epigenetic silencing of microRNA-203 is required for EMT and cancer stem cell properties. Sci Rep. 2013;3:2687.

    Article  Google Scholar 

  79. Tian L, Li M, Ge J, Guo Y, Sun Y, Liu M, Xiao H. MiR-203 is downregulated in laryngeal squamous cell carcinoma and can suppress proliferation and induce apoptosis of tumours. Tumour Biol. 2014;35:5953–63.

    Article  CAS  Google Scholar 

  80. Wellner U, Schubert J, Burk UC, Schmalhofer O, Zhu F, Sonntag A, Waldvogel B, Vannier C, Darling D, zur Hausen A, et al. The EMT-activator ZEB1 promotes tumorigenicity by repressing stemness-inhibiting microRNAs. Nat Cell Biol. 2009;11:1487–95.

    Article  CAS  Google Scholar 

  81. Jamali Z, Asl Aminabadi N, Attaran R, Pournagiazar F, Ghertasi Oskouei S, Ahmadpour F. MicroRNAs as prognostic molecular signatures in human head and neck squamous cell carcinoma: a systematic review and meta-analysis. Oral Oncol. 2015;51:321–31.

    Article  CAS  Google Scholar 

  82. Yamagishi M, Katano H, Hishima T, Shimoyama T, Ota Y, Nakano K, Ishida T, Okada S, Watanabe T. Coordinated loss of microRNA group causes defenseless signaling in malignant lymphoma. Sci Rep. 2015;5:17868.

    Article  Google Scholar 

Download references

Acknowledgements

We thank Yiqun Zhang for technical assistance in microarray analysis.

Funding

This work was supported by The Scientific and Technological Research Council of Turkey (TUBITAK, grant number: 210 T009) and partially by the United States National Cancer Institute grant P30CA125123.

Availability of data and materials

Array data have been deposited into the Gene Expression Omnibus (accession GSE69128).

Authors’ contributions

MO conceived of and participated in the design and coordination of the study and helped to draft the manuscript. OFK, IS, and BY participated in the data acquisition and drafted the manuscript. MY, HC and BG participated in the clinico-pathological analysis of patients and manuscript editing. GG participated in the participated in the data analysis, interpretation, and manuscript preparation. CJC carried out the statistical analysis. MI helped to design the study and draft the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

This study has been reviewed and approved by an institutional review board of Istanbul University, Cerrahpasa Medical School (IRB No: 35697). Patients were included into the study upon giving their written informed consents.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mustafa Ozen.

Additional files

Additional file 1:

Data Set Entire set of top differentially expressed miRNAs is given in Additional Data Set. For each miRNA probe, fold change and p value information is provided. (XLSX 30 kb)

Additional file 2: Figure S1.

(A) Relative expression levels of miR-223 in each CD133+ and CD133 sample pairs, and (B) mean relative expression levels miR-223 in CD133+ cells with respect to CD133 cells. (C) Relative expression levels of miR-328 in each CD133+ and CD133 sample pairs, and (D) mean relative expression levels miR-328 in CD133+ cells with respect to CD133 cells. (E) Relative expression levels of miR-574-3p in each CD133+ and CD133 sample pairs, and (F) mean relative expression levels miR-574-3p in CD133+ cells with respect to CD133 cells. (JPG 392 kb)

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Karatas, O.F., Suer, I., Yuceturk, B. et al. Identification of microRNA profile specific to cancer stem-like cells directly isolated from human larynx cancer specimens. BMC Cancer 16, 853 (2016). https://doi.org/10.1186/s12885-016-2863-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-016-2863-3

Keywords