MicroRNAs overexpressed in ovarian ALDH1-positive cells are associated with chemoresistance
© Park et al.; licensee BioMed Central Ltd. 2013
Received: 6 February 2013
Accepted: 13 March 2013
Published: 22 March 2013
Skip to main content
© Park et al.; licensee BioMed Central Ltd. 2013
Received: 6 February 2013
Accepted: 13 March 2013
Published: 22 March 2013
Ovarian carcinoma is the leading cause of cancer death worldwide among gynecological malignancies, and the majority of cases are related with recurrence and chemoresistance. Cancer stem cells (CSCs) are believed to be one of the causes of recurrent or chemoresistant ovarian cancer, and microRNAs are regulatory molecules newly implicated to control a variety of cellular processes, including CSCs. Therefore, we identified ovarian CSC-specific microRNAs and investigated their clinicopathological implication in ovarian carcinomas.
We isolated ALDH1 (+) cell population using the Aldefluor assay, and examined the differential expression pattern of miRNAs between ALDH1 (+) and ALDH1 (−) cells using a high-throughput microRNA microarray. We further investigated the expression patterns of differentially expressed miRNAs in human ovarian cancer samples using the real-time reverse transcription-polymerase chain reaction and analyzed their clinical impact in patients with ovarian cancer.
We found that high ALDH1 expression was associated with chemoresistance in in vitro and ex vivo samples (p = 0.024). We identified six miRNAs, including miR-23b, miR-27a, miR-27b, miR-346, miR-424, and miR-503, overexpressed in ALDH1 (+) cells, and they were significantly upregulated in chemoresistant ovarian cancer cells (1.4 ~ 3.5-fold) and tumor samples (2.8 ~ 5.5-fold) compared with chemosensitive group. Upregulation of ALDH1 (p = 0.019) and miR-503 (p = 0.033) correlated with high clinical stage, and upregulation of miR-27a was related with distant metastasis (p = 0.046) in patients with ovarian cancer.
Our findings indicate that ALDH1 is a useful marker for enriching ovarian CSCs, and high expression of ALDH1 and its related miRNAs, particularly miR-23b, miR-27b, miR-424, and miR-503, are significantly implicated in chemoresistance and tumor progression in ovarian cancer.
Ovarian carcinoma is the leading cause of gynecologic cancer deaths , and > 70% of patients with advanced stage develop a recurrence within 5 years . Despite multimodality treatment, many patients with advanced disease become refractory to conventional chemotherapeutic agents , resulting in recurrence and death.
Emerging evidence suggests that cancer stem cells (CSCs) represent the most tumorigenic and chemotherapy-resistant cells within a heterogeneous tumor mass [4, 5]. CSCs are characterized by their ability to self-renew, modulate, and balance differentiation according to the genetic background and environmental stimuli . CSCs can survive conventional treatments and become recurrent tumors that are more chemoresistant and more aggressive . However, the CSC regulatory mechanisms at the molecular level are poorly understood. The recent discovery of microRNAs (miRNA) have opened a field of gene regulation implicated in tumorigenesis and CSC modulation .
MiRNAs, noncoding regulatory RNAs of 21–25 nucleotides , are critical regulators of post-transcriptional gene expression. MiRNAs are transcribed as long RNA precursors (primary miRNAs) that are processed to yield mature miRNAs of approximately 25 nucleotides in length by Drosha-Pasha/DGCR8 and Dicer. Mature miRNAs are incorporated into the RNA-induced silencing complex and then target the 3′ untranslated region (3′-UTR) of a specific mRNA by base pairing, leading to translational repression or mRNA degradation . MiRNAs have been predicted to regulate the expression of up to 70% of human genes, implying a potential role for miRNAs in the regulation of nearly every genetic pathway [11, 12]. Taken together, understanding the regulatory role of miRNAs in CSCs may lead to a better understanding of the molecular events involved in chemoresistance and could lead to the development of a new therapeutic target.
In this study, we isolated the ALDH1 (+) cell population enriched CSCs using the Aldefluor assay, and examined the differential expression pattern of miRNAs between ALDH1 (+) and ALDH1 (−) cells using a high-throughput microRNA microarray to identify miRNAs regulating ovarian CSCs. We further investigated the expression patterns and their clinical impact of differentially expressed miRNAs in human ovarian cancer samples using the real-time reverse transcription polymerase chain reaction (qRT-PCR).
The human ovarian carcinoma cell line SKOV3 was obtained from the American Type Culture Collection (Manassas, VA, USA). Paclitaxel (PTX)-resistant cell lines (SKpac) were produced from the parent cell line (SKOV3) by continuous exposure of a stepwise, escalating concentration of PTX from an IC50 of 10% to 1000% over a period of 12 months. Seven different sublines (SKpac-8, -10, -11, -12, -13, -16, and −17) were generated. SKpac cells were 365.5-fold more resistant to PTX (IC50 = 10.21 μM, 7.59 μM, 6.77 μM, 6.57 μM, 8.69 μM, 8.38 μM, and 5.32 μM, respectively) than that of the SKOV3 cell line (IC50 = 22 nM). All cell lines were maintained in 5A medium (Gibco/Invitrogen, Carlsbad, CA, USA) with 10% fetal bovine serum (FBS) (Invitrogen), 100 U/mL penicillin, and 50 μg/mL streptomycin in a humidified atmosphere containing 5% CO2 at 37°C. The primary tumor cells were obtained at the time of surgery from patients who had undergone oophorectomy for ovarian serous carcinoma. Tumors were mechanically dissected into single-cell suspensions and incubated at 37°C for 1 hour in Ca2+/Mg2+-free PBS containing 50 U/ml collagenase A (Roche, Pleasanton, CA, USA). The enzymatic reaction was blocked by adding RPMI medium (Gibco/Invitrogen) containing 10% FBS. Cells were reacted with BerEP4-coated magnetic Dynal beads (Invitrogen) for 30 min to select epithelial cells and then cultured with RPMI medium containing 10% FBS, 1% penicillin-streptomycin, and 10 μg/ml of endothelial cell growth factor (Invitrogen).
Thirty-four ovarian carcinomas were obtained at the time of surgery from patients who had undergone oophorectomies for ovarian epithelial tumors at the CHA Bundang Medical Center. The samples were immediately frozen in liquid nitrogen and stored at −80°C. Tumor cells comprised nearly 80% of frozen section tissue.
Clinicopathological characteristics of patients with ovarian serous carcinoma (n = 34)
51.5 ± 18.5
64.5 ± 18.5
The ovarian cancer cell line (SKOV3), chemoresistant sublines (SKpac-12, 16, and 17), and primary ovarian tumor cell suspensions were counted and ALDH enzymatic activity was detected using the ALDEFLUOR assay kit (StemCell Technologies Vancouver, BC, Canada) as described by the manufacturer. Half of the cell/substrate mixture in each sample was treated with 50 mM diethylaminobenzaldehyde (DEAB), and the cells were incubated for 40 minutes. Gating was established using propidium iodide (PI) exclusion for viability, and ALDEFLUOR/DEAB treated cells were used to define negative gates. FACS was performed with 1 × 106 cells using the BD FACS Aria III (Becton Dickinson, Franklin Lakes, NJ, USA) under low pressure in the absence of UV light. FACS data were analyzed using BD FACS Diva software V6.1.3.
Primer sequence of stemness-related genes
F: 5′- TCT TGC TGC TGG TCA TTC TC-3′
R: 5′-TGC CTC ATC CTC TTC AGT TG-3′
F: 5′-GGA AAG GCT TCC CCC TCA GGG AAA GG-3′
R: 5′-AAG AAC ATG TGT AAG CTG CGG CCC-3′
F: 5′-CCA GAC ACT GCA GAT GGA GA-3′
R: 5′-AAG CTG AGG ATC TGC TGG AA-3′
F: 5′-GGC ATC ACG GTG ACT TCA AT-3′
R: 5′-GGT TTG GGG AAT GCT TCA TA-3′
F: 5′-TTA CCT GGA GAC CAG CAA GT-3′
R: 5′-CAT TAG AGC CAT TGG CAG CA′3′
F: 5′-AGA AGG CCT CAG CAC CTA C-3′
R: 5′-GGC CTG ATT GTT CCA GGA TT-3′
F: 5′-GCG CGG GCG TGA ACC AG-3′
R: 5′-CGG CGC CGG GGA GAT ACA-3′
F: 5′-ACC ACA GTC CAT GCC ATC AC-3′
R: 5′-TCC ACC ACC CTG TTG CTG TA-3′
The FACS-sorted ALDH1 (+) or ALDH1(−) cells were plated in wells of an ultra-low attachment surface 6-well culture plate (Corning, Acton, MA, USA) at a density of 3,000 cells/ml in serum-free DMEM/F12 medium (Invitrogen) supplemented with 20 ng/ml epidermal growth factor (Invitrogen), 10 ng/ml basic fibroblast growth factor (Sigma-Aldrich, St Louis, MO, USA), 0.4% bovine serum albumin (Sigma-Aldrich), and 5 μg/mL insulin (Sigma-Aldrich). Spheroid formation (each with 50–100 cells/sphere) was assessed 12 days after seeding.
miRNA expression profiles of the FACS-isolated ALDH1(+) or ALDH1(−) cells from four PTX-resistant SKpac cell lines (SKpac-16-#1,-16 #2, -17#1, and −17#2) were determined with an Affymetrix GeneChip® miRNA array to select candidate miRNAs associated with ALDH1. The samples were prepared according to the manufacturer’s instructions. Total miRNA was isolated using TriZol reagent (Invitrogen), and RNA quality was assessed with an Agilent 2100 Bioanalyzer using the RNA 6000 Nano Chip (Agilent Technologies, Amstelveen, The Netherlands). One μg of total RNA was used as input for the labeling reaction, and labeled miRNA was hybridized to the array for 16 hours at 48°C and 60 rpm as described in the protocol. After hybridization, the chips were stained and washed in a Genechip Fluidics Station 450 (Affymetrix) and scanned with a Genechip Array scanner 3000 7G (Affymetrix).
We performed real-time qRT-PCR for ALDH1 and six selected miRNAs that showed significantly altered expression in ALDH1(+) cells compared with ALDH1(−) cells to validate the ALDEFLUOR assay and the microarray data. The qRT-PCR for ALDH 1 was conducted with ALDH1(+) or ALDH1(−) cells from SKpac-12, 16, and 17 cells. One μg of total RNA was converted to cDNA with the SuperscriptIII First-strand Synthesis System (Invitrogen) and qRT-PCR was conducted with specific primers and a probe for ALDH1. For miRNAs, one μg of total RNA was reverse transcribed to cDNA using stem-loop specific RT-primers of mature miRNAs and a TaqMan MicroRNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA). The qRT-PCR was performed with a Taqman Universal Master Mix (Applied Biosystems) using the Bio-Rad CFX96 Real-Time PCR Detection System (Bio-Rad). All PCR reactions were run in triplicate, and gene expression relative to GAPDH (for ALDH1) or RNU48 (for miRNAs) was calculated using the comparative threshold method (2-ΔΔCt). Results for the ALDH1 (+) cells are expressed as fold-change relative to ALDH1 (−) cells.
We also conducted the qRT-PCR in various ovarian cancer cell lines (SKOV3, SKpac-8, 10, 11, 12, 13, 16, and 17) and 34 ovarian serous carcinoma samples including 15 chemoresistant and 19 chemosensitive carcinomas to assess the association ALDH1 and ALDH1(+)-related miRNA expression with chemoresistance. The results for the cancer cell lines are expressed as fold-changes of chemoresistant SKpac cells relative to the level of parent SKOV3 cells. The results for carcinoma samples are expressed as fold-change relative to the mean level of three normal ovarian surface epithelial cells.
Population of ALDH1 (+) cells as determined by the ALDEFLUOR assay in various human ovarian cancer cell lines and primary cells
% ALDH1 (+) cell
0.9 ± 0.2
4.46 ± 0.71
1.66 ~ 7
We performed qRT-PCR for ALDH1 mRNA to confirm whether FACS-sorted ALDH1 (+) cells showed increased ALDH1 expression. The FACS-sorted ALDH1(+) population showed 2–5 fold higher expression of ALDH1 mRNA compared with that in the ALDH1(−) population (Figure 1C).
We compared the spheroid formation capability between ALDH1 (+) and ALDH 1(−) cells to test the stem cell-like property. ALDH1(+) cells showed increased size and number (106 ± 37) of spheroids compared with ALDH1(−) cells (63 ±15) (Figure 2B) (t-test, p = 0.05).
To assess the role of ALDH1 in ovarian cancer, we analyzed the relationship of ALDH1 mRNA expression with clinicopathological parameters, including clinical stages, lymph node metastasis and distant metastasis. For comparison, patients were classified into high (≥2-fold) and low (<2-fold) ALDH1 expression groups. The high expression of ALDH1 mRNA was significantly (p = 0.019) associated with advanced clinical stage (Figure 3C). The ALDH 1 expression was not correlated with lymph node or distant metastases.
MicroRNAs significantly altered (>1.5 fold) in ALDH1(+) cells compared with ALDH1 (−) cells in ovarian cancer
Absolute fold change
Putative target genes
WEE1, CDCA4, USP15, LATS2, HIPK2
GSPT1, CDS1, TAB3, SFRP1, MDM4, PRKCB, SESN2, FOSB, SMAD2
APAF1, PPP2R5E, PPP1CB, PPIF, REPS2, CHUK
GSPT1, CDS1, TAB3, PRKCB
We compared expression of six differentially expressed miRNAs between ALDH1(+) and ALDH1(−) cells using real-time RT-PCR to validate the microarray results. As a result, miR-424 (1.62-fold), miR-346 (3.25-fold), miR-503 (1.66-fold), miR-27a (2.08-fold), miR-23b (1.98-fold), and miR-27b (3.09-fold) were upregulated in ALDH1 (+) cells relative to ALDH1 (−) cells (Figure 4B).
The expression levels of miRNAs in the human cancer tissue samples relative to those of normal epithelial cells were calculated and compared between the chemosensitive and chemoresistant groups. The expression levels of miR-23b (3.5-fold vs. 0.6-fold p = 0.037), miR-27b (5.5-fold vs. 1.6-fold, p = 0.040) and miR-424 (2.0-fold vs. 0.7-fold, p = 0.047) were significantly higher in the chemoresistant group than those in the chemosensitive group (Figure 5B).
We analyzed the correlation between the expression of these miRNAs and clinicopathological parameters, including clinical stage, lymph node metastasis, and distant metastasis to verify the role of ALDH1(+)-associated miRNAs in ovarian cancer. Patients were classified into high (≥1.5-fold) and low (<1.5-fold) expression groups. High expression of miR-503 (Figure 5C) was significantly associated (p = 0.033) with advanced clinical stage (stage III and IV), and upregulation of miR-27a (Figure 5D) was related to distant metastasis (p = 0.046). No other miRNAs were associated with the clinicopathological parameters.
CSCs comprise a very small proportion (0.01–1.0%) of malignant tumors but are believed to be the source of recurrent tumors and of resistance to chemotherapy and radiotherapy due to their ability to survive conventional treatments, which usually target rapidly dividing cells [7, 13]. Therefore, identifying CSCs and their molecular phenotype might serve to develop a new therapeutic strategy for lethal human cancers including ovarian carcinoma. Indeed, several studies focusing on identifying and characterizing ovarian CSCs have been reported recently [14–16]. These studies suggested several markers including CD44, CD117, CD133, and ALDH1 as possible ovarian CSC markers. Among them, ALDH1 is considered to be a consistent ovarian CSC marker identified in all ovarian cancer cell lines and primary human cancer samples . Therefore, we used ALDH1 as a putative ovarian CSC marker in the present study and confirmed that ALDH1(+) cells were enriched with ovarian CSCs by demonstrating an increased expression of stem cell markers and increased spheroid formation ability compared with those in ALDH1(−) cells.
A growing body of evidence indicates that miRNAs regulate most human genes in crucial biological processes including tumorigenesis, progression, and therapy resistance. MiRNAs also seem to be involved in controlling CSC self-renewal and differentiation, according to more recent studies reporting differentially expressed miRNAs in CSCs of various human cancers [17–19]. Recently, deregulation of some miRNAs, such as miR-200a, miR-199a, and miR-214, were reported in CD133(+) or CD44(+) ovarian CSCs . However, only limited data are available regarding microRNA expression profiles in ovarian CSCs. Therefore, we identified ovarian CSC-specific miRNAs using high-throughput microRNA microarray and assessed the association of expression of differentially expressed miRNAs with chemoresistance in ovarian cancer.
We found that six miRNAs, including miR-23b, miR-27a, miR-27b, miR-346, miR-424, and miR-503, were significantly overexpressed in CSC-enriched ALDH1(+) cells using miRNA microarray and qRT-PCR. We demonstrated that these miRNAs and ALDH1 increased significantly in chemoresistant SKpac sublines compared with those in the chemosensitive SKOV3 cell line, suggesting that these miRNAs and ALDH1 are associated with ovarian cancer chemoresistance. To validate the relationship of ALDH1 and ALDH1-related miRNAs with chemoresistance, we further studied the expression of ALDH1 and these miRNAs in human ovarian cancer samples. Consistent with the results of the ovarian cell lines, ALDH1 mRNA expression in the chemoresistant group was significantly higher (2.5-fold) than that in the chemosensitive group, and the expression of miR-23b, miR-27b, and miR-424 was significantly upregulated in the chemoresistant group. Upregulation of ALDH1 and miR-503 correlated with high clinical stage, and upregulation of miR-27a was correlated with distant metastasis.
Our result showing the correlation between high ALDH1 expression and chemoresistance agrees with some other recent breast and ovarian cancer studies [21–23]. The present study also demonstrated that PTX-resistant subclones (SKpac cells) were more enriched in ALDH1(+) cells compared with their parental SKOV3 cells. We also found higher ALDH1 expression in chemoresistant ovarian cancer tissues compared with that in the chemosensitive group. It is possible that clonal selection of intrinsic ALDH1(+) cells by chemotherapy may increase ALDH1 expression in chemoresistant cancer cells and may allow the acquisition of chemoresistance in residual cancers. Chang et al.  reported a correlation between low ALDH1 expression and advanced cancer stage and response to cisplatin treatment. However, Chang et al. included heterogeneous subtypes of ovarian cancer, including the endometrioid and mucinous types, and used immunohistochemical staining with a tissue microarray, which may not represent the characteristics of the whole tumor. More recent studies [15, 23] have suggested an association between high ALDH1 expression and advanced cancer stage and poor response to chemotherapy, consistent with our results. Taken together, targeting ALDH1 in ovarian cancer may be a new strategy to treat chemoresistant and advanced ovarian carcinoma.
In the present study, miR-23b, miR-27b, and miR-424 were significantly upregulated in the chemoresistant human cancer group among six miRNAs differentially overexpressed in ALDH1(+) cells. MiR-23b has been reported to function as a renal cancer oncogene , similar to ovarian cancer seen in the present study, whereas it functions as a tumor suppressor in hepatocellular carcinoma . The introduction of pre-miR-27b stimulates invasion in breast cancer cells by targeting the ST14 tumor suppressor , and is upregulated in glioma cells and tumor tissues . However, it is downregulated in lung cancer tissue . A recent study for miRNAs expression signatures in tumorigenesis  indicated that miR-424 was associated with angiogenesis, and that miR-23b and miR-27b were upregulated in metastatic stage. They also demonstrated that miR-424 were up-regulated in metastases compared with primary tumors in both mice and human pancreatic neuroendocrine tumors. Another study for non-small cell lung cancers demonstrated that up-regulation of miR-424 was related with poor prognosis . On the other hand, it was reported that low expression of miR-424 was positively correlated with advanced clinical stage, lymph node metastasis and other poor prognostic parameters in cervical cancers . Taken together, the oncogenic function of miR-23b, miR-27b, and miR-424 seems to be cell type- and context-specific.
The upregulation of miR-27a didnot show a statistically significant correlation with chemoresistance in the validation with paclitaxel-resistant SKpac sublines and patient’s samples in this study. However, a previous study reported that miR-27a was upregulated in paclitaxel-resistant ovarian cancer cell line A2780/Taxol as compared with its parental line A2780 . In the present study, high expression of miR-27a was related with distant metastasis suggesting its role on the progression of ovarian cancer.
The function of miR-503 in tumor development and progression remains unresolved except for a report demonstrating its regulation of metastatic function in hepatocellular cancer cells . Therefore, we searched predicted target genes associated with apoptosis and tumor suppression via web-based computational programs such as Target Scan (http://www.targetscan.org) and PicTar-Vert (http://pictar.mdc-berlin.de) to predict the role of miR-503 implicated in chemoresistance and tumor progression. CDCA4, a putative miR-503 target gene with a high target score, participates in the regulation of cell proliferation mainly through the E2F/pRB pathway . Further studies are necessary to determine whether miR-503 induces cancer cell growth or migration/invasion, and whether CDCA4 is a direct target gene of miR-503, as well as to define the oncogenic function of miR-503 in this subset of ovarian cancers.
We demonstrated in this study that ALDH1 is a marker for enriching ovarian CSCs, and that high ALDH1 expression is associated with chemoresistance and high clinical stage. Six miRNAs, including miR-23b, miR-27a, miR-27b, miR-346, miR-424, and miR-503 were overexpressed in ALDH1(+) cells, and significantly implicated in chemoresistance in ovarian cancers. Upregulation of miR-503 correlated with high clinical stage, and upregulation of miR-27a was related with distant metastasis in patients with ovarian cancer.
This work was supported by grants funded by Korean Government, KHIDI-A111804 and NRF-2012-R1A1B3004095.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.