Home Biomarkers to target and silence stemness of breast cancer stem cell model: silencing MDR1 by siRNA
Article Open Access

Biomarkers to target and silence stemness of breast cancer stem cell model: silencing MDR1 by siRNA

  • Gamze Yıldırım ORCID logo , Meltem D. Kars ORCID logo EMAIL logo , Gökhan Kars ORCID logo and Hamdi Ş. Kılıç ORCID logo
Published/Copyright: March 28, 2022

Abstract

Objectives

Aim of the study was to reveal new biomarker genes to target breast cancer stem-like cells (BCSC-like) and then sensitize BCSC-like cells to chemotherapy by silencing MDR1 gene found to be the most suitable target.

Methods

Drug resistance associated genes were screened by cDNA microarray to unveil biomarker genes in drug resistance breast cancer model cells. Drug resistance was then reversed by silencing MDR1 gene in BCSC-like cells. The effect of silencing was monitored by real-time cell proliferation analysis. Differential expressions of MDR1, ALDH1A3, EGFR and BAG4 genes were identified by real-time PCR. P-glycoprotein (P-gp) expression level and its activity were investigated by Western blot and flow cytometry measurements, respectively.

Results

16 new biomarker genes were identified upon gene expression analysis by cDNA microarray. MDR1 gene was selected as the most potent target gene and silencing of it caused down-regulation of MDR1, ALDH1A3, EGFR, BAG4 expression and P-glycoprotein activity and expression in BCSC-like cells. At the end, silenced BCSC-like cells were found to be more responsive to paclitaxel therapy.

Conclusions

In conclusion, siMDR1 silencing is an effective way to reverse multidrug resistance and malignancy. New biomarker genes revealed in this study require to be investigated to target stemness of BC.

Introduction

Chemotherapy is a common treatment strategy against breast cancer [1]. However, multidrug resistance (MDR) is encountered frequently in cancer patients treated with chemotherapeutic agents [2]. Breast cancer cells develop resistance to chemotherapy agents by several mechanisms. These molecular mechanisms lead to MDR phenotype. The most known mechanism of drug resistance is increased pumping out of therapeutic agents due to highly expressed ATP binding cassette (ABC) transporter proteins. MDR causes the development of malignancy after exposure to chemotherapy as well [2]. There are several genes and proteins associated with multidrug resistance. One of the remarkable genes is ALDH1A3 that encodes an aldehyde dehydrogenase enzyme whose increased activity has been tightly correlated with metastasis. Cancer stem cells with high ALDH1A3 activity bear invasive characters and drug resistance phenotype [3]. Epidermal growth factor receptor (EGFR) is also associated with MDR. EGFR is a type of receptor having tyrosine kinase activity and functions in most of the cancer types by promoting cell proliferation, differentiation and cell survival [4]. EGFR also promotes AKT pathway which then stimulates cell growth mechanisms. Activated or deregulated EGFR signaling results in drug resistance in breast cancer cells to commonly anticancer drugs such as doxorubicin, daunorubicin, paclitaxel, cisplatin and 5-flurouracil as well as ionizing radiation (IR) [5]. The other MDR related protein is Bcl-2 Associate Athanogene 4 (BAG4) which is an anti-apoptotic protein and is called silencer of death domains. It is generally highly expressed in cancer tissues and provides a mechanism to escape from death [6]. BAG4 protein was defined in radiation resistant HeLa cells and BAG4 has been also linked both to cisplatin resistance in various cancer types [7]. In the current study, ALDH1A3, EGFR and BAG4 genes were found to be correlated with paclitaxel resistance in BCSCs.

The parental MCF-7 cells used here do not express MDR1 gene. The parental MCF-7 cells were treated with increasing doses of paclitaxel for approximately two years to develop drug-resistant model cells [2]. As the cells acquired paclitaxel resistance, the cell line (MCF-7/Pac) became unresponsive to paclitaxel therapy and it started to express MDR1 and breast cancer stem cell associated genes differentially with respect to parental drug sensitive cell line. In breast cancer stem cell sub populations, highly expressed CD44, ALDH and CD24 are consensus biomarkers. It is known that about two hundred of CD44 positive, CD24 negative and ALDH overexpressing cancer cells are sufficient to initiate tumors [8]. In line with this, the cells used in this study were CD44+, CD24-, ALDH+, P-gphigh and paclitaxel resistant BCSC-like cells. The method of obtaining BCSC-like cells was patented by our group [9].

Gene silencing has received increased attention in functional genomics research. The method has enabled elucidation of the roles of genes with unknown functions in model organisms [10]. The method of targeting mRNA is named ‘antisense strategy’ and has been widely applied to assess the function of genes. These antisense strategies aim to down-regulate the mRNA encoding the protein to be inhibited [11]. Taking into account the aforementioned arguments, the aim in this study was to reveal target genes for the reversal of stemness of BCSC-like cells.

Materials and methods

Cell lines and culture conditions

Estrogen receptor (ER) positive mammary carcinoma epithelial cells (MCF-7) and ER BCSC-like cells were used in the present study. Paclitaxel resistant breast cancer cell line was developed from drug sensitive MCF-7 by stepwise dose increments and it was named as MCF-7/Pac by our group [2, 12]. Since MCF-7/Pac cell population exhibited breast cancer stem cell phenotype, BCSC-like cell population was sorted from them. The cells have been selected by flow cytometry from a paclitaxel resistant cell population by using ALDH, CD44 and CD24 as BCSC markers. The method for isolation of chemotherapy resistant BCSC-like cell line was patented by Turkish Patent and Trademark Office with the registration number TR 2013 10837 B and was published by our group [9, 13]. MCF-7 and BCSC-like cell lines were cultured in RPMI1640 (Biochrom, Berlin, Germany) solution at 37 °C and in the atmosphere of 5% CO2.

Analysis of cDNA microarray data

High throughput total genome cDNA microarray screening (Affymetrix, Santa Clara, CA, USA) was performed to figure out differentially expressed genes in Paclitaxel resistant MCF-7 cells and drug sensitive MCF-7 cells previously [14]. After scanning of the gene chip arrays, the statistically significant data were selected from high throughput data sets by independent sample t-test (p-Value <0.05) between duplicate data for drug resistant and sensitive cells. In the context of this work, statistically significant microarray results (p-Value <0.05) were re-analyzed in silico. Significantly altered genes between resistant subline and MCF-7/S were listed and gene trees were constructed from these lists by standard correlation (Figure 1). Upregulated and downregulated genes were selected from gene trees by volcano plot filtering (by excluding the genes differentially expressed between −2 and 2 folds) to find out new BCSC related target genes.

Figure 1: 
Gene tree constructed during the analysis of microarray gene expression data.
When we compare the drug resistant cells and drug sensitive cells, the upper part of the tree demonstrates the gene clusters of which expression levels are low (blue) in sensitive cells and high in (yellow-orange-red) in drug resistant cells. The bottom part of the tree demonstrates the gene clusters of which expression levels are low (blue, yellow) in drug resistant cells and high in (orange-red) in drug sensitive cells.
Figure 1:

Gene tree constructed during the analysis of microarray gene expression data.

When we compare the drug resistant cells and drug sensitive cells, the upper part of the tree demonstrates the gene clusters of which expression levels are low (blue) in sensitive cells and high in (yellow-orange-red) in drug resistant cells. The bottom part of the tree demonstrates the gene clusters of which expression levels are low (blue, yellow) in drug resistant cells and high in (orange-red) in drug sensitive cells.

siRNA transfection

P-gp encoding MDR1 gene was targeted to be reversed by siRNA as the most overexpressed gene in BCSC-like population was found to be MDR1 according to the results of microarray analysis. First of all, siRNA transfection conditions were optimized by transfecting the cells with Cy3 labeled siRNA (Silencer™ Cy™3-labeled Control, Ambion, Paisley, UK). Then, the cells were observed under confocal microscope (Nikon, Japan, Cy3 is excited at 550 nm and emits at 570 nm) for visualization of the transfected fluorescent siRNA in the cells and for optimization of siRNA transfection protocol. The cells were transfected with Cy3 labeled siRNA (5 pmol) and then they were cultured in 24-well cell view culture plates (Greiner Bio One, Frickenhausen, Germany) for 24 and 48 h to find the optimum incubation time. MDR1 siRNA (IDT- Iowa, USA) and siRNA MOCK (non-template siRNA, Invitrogen) oligonucleotide sequences are presented in Table 1. BCSC-like cells were seeded (1 × 106 cells each well) in six well culture plates and incubated for 24 h. Then, optimum amount of lipofectamine and siRNA (25 pmol) solutions were applied to cells and they were cultured for 24 h at 37 °C in the atmosphere of 5% CO2. The procedure applied in this study was adopted from the report of Dönmez et al. [15].

Table 1:

Oligonucleotide sequences of siRNA and control (MOCK).

siRNA Oligonucleotide sequence
MDR1 sense 5′GGAUGUGAGUUGGUUUGAUdTdT3′
MDR1 antisense 5′AUCAAACCAACUCACAUCC dTdT3′
MOCK sense 5’GGAAGUUUGUGGUUUGGAUdTdT’3
MOCK antisense 5′AUCCAAACCACAAACUUCC dTdT′3

Analysis of in vitro paclitaxel therapy by a real-time cell analysis system

It was previously demonstrated that the BCSC-like cells used here were resistant to paclitaxel, a microtubule inhibiting anti-cancer drug [13]. The effects of paclitaxel on the proliferation of silenced BCSC-like and siMOCK-treated BCSC-like cells were analyzed using xCELLigence (RTCA system, ACEA Biosciences, San Diego, CA, USA). The method was applied as described previously [13]. In summary, BCSC-like cells were first seeded (10.000 cells/well) and silenced by incubating with siRNA for 24 h on e-plates. After that, paclitaxel was added to the cells at concentrations ranging from 20 to 1.5 µM. Finally, the effect of paclitaxel on siMDR1 silenced breast cancer stem cells was compared to the non-silenced drug resistant BCSC-like cells.

Gene expression analysis

RNA isolation from siMDR1 and siMOCK treated BCSC-like cells and un-transfected control cells were performed by combining two RNA isolation protocols as follows. First, total RNA was extracted from cells (1 × 106) by TRI Reagent (Sigma, MO, USA). Then, the aqueous RNA phase was collected and diluted in 70% ethanol. Finally, the RNA suspension was purified using Qiagen RNA isolation kit (Manchester, UK). After the final elution of RNA samples, the quality of RNA samples was determined by agarose gel electrophoresis (Biorad, California, USA) and nanodrop (Thermo Scientific, Massachusetts, USA) measurements. Following RNA extraction, the effect of MDR1 silencing on the expression level of BCSC associated genes were investigated by reverse transcription polymerase chain reaction (RT-PCR). The specific primers for detecting the genes of interest were designed by utilizing Primer3 and NCBI BLAST programs (Table 2).

Table 2:

Primer sequences for RT-PCR.

Primers Sequence Product, bp Annealing T, °C
MDR1 sense 5′-ACAGAAAGCGAAGCAGTGGT-3′ 62 57
MDR1 antisense 5′-ATGGTGGTCCGACCTTTTC-3′
EGFR sense 5′-GAGCACTTGCTGAATGGTGA-3′ 88 57
EGFR antisense 5′-CCTCCACTCCAAACCAGTGT-3′
ALDH1A3 sense 5’ -CTGGGTCAGTGTTCCCTTGT-3′ 72 59
ALDH1A3 antisense 5′-TCACAGTGCCAGAAATGCTC-3′
BAG4 sense 5′- AGGACCGACTGTACGACCAC-3′ 69 57
BAG4 antisense 5′-CTGATGAAGGCCAGGGATAA-3′
GAPDH sense 5′- CGAGATCCCTCCAAAATCAA-3′ 58 57–59
GAPDH antisense 5′- TGGACTCCACGACGTACTCA-3′

cDNA was synthesized from total RNA by reverse transcription using GAPDH, MDR1, EGFR, ALDH1A3 and BAG4 specific primers (Revert-aid RT enzyme, Fermentas, Ontario, Canada). Gene expression analysis was performed by a real-time thermal cycler (ABI-7500 system, Thermo Scientific, ABD). The reaction mixture consisted of SYBR Green master mix (Biorad, California, USA), cDNA and primers. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) encoding gene was used as internal reference. The 2−ΔΔCt analysis method was applied to calculate the expression levels of the target genes in silenced and non-silenced cells [16].

Evaluation of expression levels of P-gp in BC cells

The P-gp was selected as a representative protein for confirmation of RT-PCR gene expression analysis. For this, total protein was extracted from cells and Western blot analysis was applied for the target protein P-gp. Total lysate of silenced and MOCK control BCSC-like (5 × 106) cells were prepared in RIPA cell lysis buffer (Amresco, USA) as previously mentioned [13]. The protein concentration of the cell lysate was determined as previously described [13]. Protein lysates were analyzed on SDS-PAGE and transferred onto 0.45 µm nitrocellulose membrane for 5 h using BioRad Electro Blot System. The membrane was hybridized with antibodies (Rabbit Anti-Human GAPDH antibody, 1:200 dilution; Rabbit Anti-Human-Pgp antibody, 1:1000; ABCAM). Then membrane was hybridized with secondary antibody (HRP conjugated anti-Rabbit for 1 h). After washing, the detection buffer ECL substrate (Abcam, USA) was added the membrane was analyzed under chemiluminescent gel imaging system and density of each protein band was measured (Vilber-smart imaging, Cedex, France). The band densities of P-gp was normalized by the band densities of GAPDH for drug resistant and silenced BCSC-like cells.

Confirmation of downregulation of P-gp activity by Rhodamine-123 assay

Rhodamine-123 (Sigma Aldrich, MO, USA) is a fluorescent compound that is effluxed out of the cell by an active P-gp. Rhodamine is generally used to determine the activity of P-gp through flow cytometry [2]. Briefly, drug sensitive MCF-7 cells, silenced and MOCK treated BCSC-like cells were incubated with 1 µM rhodamine in serum free media for 10 min at 37 °C. Then, the cells were centrifuged and washed in PBS. Finally, fluorescent cells were counted and fluorescent values were measured by a flow cytometer (FACS AriaIII, Becton Dickinson, NJ, USA.). The accumulation of rhodamine results in high fluorescent value (y axis of flow cytometry histogram) which is the designation of drug accumulation in the cells. The decrease in mean fluorescent value indicates that P-gp exports rhodamine out of the cells. The fluorescent activity ratio (FAR) demonstrates the degree of reversal of drug resistance. FAR is calculated by dividing the mean fluorescent value of inhibited drug resistant cells by the mean fluorescent value of untreated drug resistant cell. When the FAR value is higher than 1.10, it shows that the drug resistant cells experience reversal of drug resistance or decrease in P-gp activity [2].

Statistics

Duplicate set of high throughput microarray results were firstly analyzed by Student’s t-test (p-Value <0.05) and then filtered by Volcano plot (by excluding genes with fold change values between −2 and 2 folds) to select the genes that were differentially expressed in drug resistant cells with respect to drug sensitive cells. The analyses regarding toxicological assessments, gene expression analysis and flow cytometry measurements were repeated at least three times and subjected to Student’s t-test (p-Value <0.05).

Results

Whole genome cDNA analysis

Whole genome cDNA microarray expression analysis results revealed that in addition to MDR1, ALDH, CD44, CD24, EGFR, SNAI2, SMAD3 and WNT-3, 16 cancer stem cell associated genes were differentially expressed in MCF-7/Pac cells (containing BCSC-like cells) with respect to drug sensitive MCF-7 cells (Table 3) [13, 14]. According to the results, MDR1 gene is the most overexpressed gene (95 folds) in BCSC-like populated drug resistant cells. Additionally, MUC18, AKT3, RAB3B, WNT5A, LPL, CAPR, MAP3K12, MAP4K4, MAP4K1, STAT6B, BAG4, MAPKKK4 genes were found to be overexpressed about two to 25 folds more in BCSC-like cell populated cells compared to sensitive BC model. On the other hand, MAP3K3, MAPKAPK2, MAP3K1, VMP1 were down regulated about three to six folds. In the present study, 16 differentially expressed genes were disclosed which will provide insight into studies for targeting stemness of BC.

Table 3:

Whole genome cDNA microarray analysis of the drug resistance related genes in drug resistant BCSCs with respect to drug sensitive BCCs.

Gene Description Fold change ratioa Reference
MDR1 ATP-binding cassette, sub-family B (MDR), member 1 95.18 [13]
EGFR Epidermal growth factor receptor 33.39 [14]
MUC18 CD146, melanoma cell adhesion molecule 24.65 This work
AKT3 v-akt murine thymoma viral oncogene homolog 3 24.37 This work
RAB3B RAB3B, member RAS oncogene family 23.63 This work
SMAD3 SMAD, mothers against DPP homolog 3 17.26 [14]
WNT5A Wingless-type MMTV integration site fam. Member 5A 16.38 This work
ALDH1A3 Aldehyde dehydrogenase 1 family, member A3 16.79 [13]
SNAI2 Snail homolog 2 (Drosophila) 9.249 [14]
LPL Lipoprotein lipase 4.938 This work
CAPR Catenin (cadherin-associated protein), alpha 2 3.868 This work
WNT3 Wingless-type MMTV integration site family, member 3 3.676 [13]
MAP3K12 Mitogen-activated protein kinase kinase kinase 12 2.716 This work
MAP4K4 Mitogen-activated protein kinase kinase kinase kinase 4 2.602 This work
MAP4K1 Mitogen-activated protein kinase kinase kinase kinase 1 2.403 This work
STAT6B Signal transducer and activator of transcription 6 2.176 This work
BAG4 BCL2-associated athanogene 4 2.047 This work
MAPKKK4 Mitogen-activated protein kinase kinase kinase 4 2.043 This work
MAP3K3 Mitogen-activated protein kinase kinase kinase 3 −2.398 This work
CD24 CD24 antigen −2.604 [13]
MAPKAPK2 Mitogen-activated protein kinase-activ. protein kinase 2 −2,825 This work
MAP3K1 Mitogen-activated protein kinase kinase kinase 1 −4.902 This work
VMP1 Transmembrane protein 49 −6.452 This work
  1. aExpression level of drug resistance BCS-like cells/drug sensitive BCCs.

siRNA transfection

siRNA transfection process was optimized by using Cy-3 labelled fluorescent siRNA and by visualizing it inside BCSC-like cells under a confocal microscope. A 24 h incubation period was found to be suitable for transfection (Figure 2). Accordingly, the antisense MDR1 RNA transfection was applied for 24 h for all assays including cell proliferation, gene expression, P-gp activity and expression.

Figure 2: 
Observation of Cy3 labeled control siRNA within the cells under confocal microscope (Nikon, Japan, Cy3 was excited at 550 nm and emission was at 570 nm).
The images obtained by transmission detector (to observe the cells) and by standard detector (to observe fluorescence) were merged by the image analysis program (NIS elements viewer) of confocal microscope. The images were obtained after (A) 24 h incubation, (B) 48 h incubation with siRNA-Cy3 solution and PBS washings.
Figure 2:

Observation of Cy3 labeled control siRNA within the cells under confocal microscope (Nikon, Japan, Cy3 was excited at 550 nm and emission was at 570 nm).

The images obtained by transmission detector (to observe the cells) and by standard detector (to observe fluorescence) were merged by the image analysis program (NIS elements viewer) of confocal microscope. The images were obtained after (A) 24 h incubation, (B) 48 h incubation with siRNA-Cy3 solution and PBS washings.

Cell proliferation assays

Assessment of the effect of chemotherapy on BCSC-like cells was conducted by a Real Time Cell Analysis system. The responses of siMDR1 silenced and siMOCK treated BCSC-like cells to paclitaxel were evaluated in various time periods by treating the cells with different drug concentrations (Table 4 and Figure 3). According to the results, the antiproliferative effect of 1.25 µM paclitaxel application on siRNA silenced BCSC-like cells was much more pronounced than that on control-MOCK cells. On the other hand, 20 µM paclitaxel resulted in similar anti-proliferative effect on both siMDR1 treated and control BCSC-like cells. The anti-proliferative effects of increasing paclitaxel concentrations on siMDR1 and siMOCK treated BCSC-like cells are clearly figured out as time and dose-dependent way (Figure 3). The IC50 values computed at the logarithmic phase of cell proliferation are listed in Table 4. Regarding the results, it can be deduced that BCSC-like cells become sensitive to chemotherapy after MDR1 silencing.

Table 4:

IC50 values (time dependent manner), *Resistance index (RI) was calculated from the formula: resistance index = IC50 mock control/IC50 siMDR1.

Hour (at log phase) IC50, µM (R 2≥0.90)
BCSC-like MOCK BCSC-like siMDR1 Resistance index
104:59:35 1.71 1.01 1.69
109:59:54 1.66 0.85 1.95
115:00:13 1.59 0.73 2.18
120:00:38 1.55 0.64 2.42
125:00:57 1.54 0.71 2.17
130:01:20 1.53 0.70 2.19
135:01:39 1.49 0.67 2.22
140:02:00 1.49 0.83 1.80
Figure 3: 
Real-time cell proliferation analysis of (A) BCSC-like (MOCK) cells and silenced cells treated with 1.25 µM paclitaxel, (B) BCSC-like (MOCK) cells treated with 1.25–20.0 µM paclitaxel, (C) siRNA silenced BCSC-like cells treated with 1.25–20.0 µM paclitaxel.
Figure 3:

Real-time cell proliferation analysis of (A) BCSC-like (MOCK) cells and silenced cells treated with 1.25 µM paclitaxel, (B) BCSC-like (MOCK) cells treated with 1.25–20.0 µM paclitaxel, (C) siRNA silenced BCSC-like cells treated with 1.25–20.0 µM paclitaxel.

Gene expression analysis upon MDR1 silencing

The 2−∆∆CT method was used to examine changes in gene expression levels. The results demonstrated that siRNA silencing of the MDR1 gene dramatically down-regulated the expression levels of MDR1, ALDH1A3, EGFR and BAG4 genes. Indeed, MDR1, ALDH1A3, EGFR and BAG4 gene expression levels decreased about 94, 99, 99 and 50% respectively. There was not any significant alteration in gene expression levels in mock treated BCSC-like cells (Figure 4A). The whole genome cDNA microarray analysis was previously performed by our group to evaluate differential gene expressions between paclitaxel sensitive MCF-7 and paclitaxel resistant MCF-7 cells [17]. These results showed that ALDH1A3, EGFR and BAG4 expression levels were 16.79, 33.39 and 2.05 folds higher respectively in drug resistant cells when compared to drug sensitive cells (Table 3). Here, RT-PCR results evidenced that sensitization of drug resistant BCSC-like cells by siMDR1 reverses the cells to the drug sensitive form by silencing MDR1 gene.

Figure 4: 
(A) results of gene expression analysis for MDR1, ALDH1A3, BAG4 and EGFR genes by 2−∆∆CT method. (B) Western blot analysis for P-gp/GAPDH and densitometry values of the normalized membrane blot bands 1. BCSC-like (resistant) 2. BCSC-like (silenced).
Figure 4:

(A) results of gene expression analysis for MDR1, ALDH1A3, BAG4 and EGFR genes by 2−∆∆CT method. (B) Western blot analysis for P-gp/GAPDH and densitometry values of the normalized membrane blot bands 1. BCSC-like (resistant) 2. BCSC-like (silenced).

Evaluation of P-gp expression level and drug efflux activity

According to the Western blot analysis, the representative protein expression, P-gp, was down-regulated by siRNA silencing of MDR1 (according to Image J densitometry analysis) (Figure 4B). P-gp down-regulation was also examined and revealed by functional analysis through measuring P-gp activity by flow cytometry. When the mean fluorescent values of histograms were analyzed, it was observed that MCF-7 drug sensitive cells accumulated about 5.2 folds more fluorescence than drug resistant BCSC-like/MOCK cells (Figure 5). Additionally, FAR value for the siRNA silenced BCSC-like cells was calculated as 1.26. This clearly demonstrates reversal of MDR1 expression and P-gp activity in BCSC-like cells. The flow cytometry histograms for each cell line were given as supplementary materials (S1A, B and C).

Figure 5: 
Rhodamine exclusion assay results:
Fluorescence accumulation (FL values) in sensitive MCF-7 cells, in drug resistant BCSC-like cells (MOCK) and in silenced BCSC-like cells. As the P-gp activity decreased due to silencing, rhodamine accumulation increased, which meant increased potential for drug accumulation.
Figure 5:

Rhodamine exclusion assay results:

Fluorescence accumulation (FL values) in sensitive MCF-7 cells, in drug resistant BCSC-like cells (MOCK) and in silenced BCSC-like cells. As the P-gp activity decreased due to silencing, rhodamine accumulation increased, which meant increased potential for drug accumulation.

Discussion

Breast cancer is one of the most widespread cancer diseases worldwide. Chemotherapy resistance in breast cancer is still an important difficulty encountered in clinic [18]. Drug resistance against paclitaxel has been reported to develop in nearly 90% of metastatic breast cancer patients [19]. Over-expression of MDR1 gene coding for P-gp plays a significant role in drug resistance by exporting drugs and thereby preventing them reaching their intracellular targets [20]. In this study, the aim was to reveal new marker genes to target stemness of BCSC-like cells based on whole genome cDNA microarray results. And following, the mostly overexpressed MDR1 gene was downregulated by siRNA mediated silencing technique. Breast cancer stem-like cells that were ALDH positive, CD44 positive, CD24 negative, P-glycoprotein positive, ER (estrogen receptor) negative were used in this study. The cell line used in the current work was previously developed and established from drug sensitive MCF-7 cells and was validated by in vitro chemotherapy, toxicological tests and pharmacogenomics analysis [13, 21]. High throughput gene expression analysis results exhibited that BCSC-like cell containing drug resistant cells differentially expresses CSC associated genes significantly when compared to drug sensitive parental MCF-7 cells. For instance, MDR1, EGFR, MUC18, AKT3, RAB3B, SMAD3, WNT5A, ALDH1A3, SNAI2 genes were expressed between nine to 95 folds more in the drug resistant cell population that was crowded by BCSC-like cells with respect to drug sensitive MCF-7 model cell. As a result of this, these genes may be proposed as target genes to reverse drug resistance and to target stemness of BC. In addition, LPL, CAPR, WNT3, MAP3K12, MAP4K4, MAP4K1, STAT6B, BAG4, MAPKKK4 genes were found to be upregulated in drug resistant BC cells about two to five folds. On the other hand, MAP3K3, CD24, MAPKAPK2, MAP3K1 and VMP1 were significantly downregulated in drug resistant BC. In concordance with our results, overexpression of BAG4, MDR1, EGFR, AKT3, SMAD3, ALDH1A3, SNAI2, WNT3, STAT6B, CAPR genes were declared to be closely associated with cancer development, proliferation and drug resistance in many studies [22], [23], [24], [25], [26], [27]. SMAD3 induces stemness in breast cancer (BC) cells and inhibition of SMAD3 results in blockage of the SMAD3 associated TGF-β signal pathway and cancer stem cell (CSC) features [22]. In normal breast tissue, Wnt signaling controls important cellular events such as cell fate determination, proliferation, migration and tissue homeostasis. Noncanonical signaling is activated by Wnt4, 5A, 5B and Wnt 11. Wnt5A is known to be one of the regulators of macrophage mediated malignant invasion in BC [23]. Snail transcription repressor 2 (SNAI2) is highly associated with epithelial to mesenchymal transition (EMT) in BC and was reported as one of the CSC markers [24]. A tumor suppressor, α-catenin (cadherin associated protein-alpha, CAPR) is known as one of the components of the E-cadherin and catenin complex and it takes role in maintaining the stability of the intercellular adherent junctions [25]. STAT6 can act as both transcription factor and signaling molecule. STAT6 is highly associated to IL-4 and IL-13 signal pathway, and also has role in polarization of the T helper cells in immune system. STAT6B is one of the isoforms of STAT6 and has about 110 amino acids deletion at the N- terminal region but the product is intact concerning its biological functions [26]. It was reported that STAT6 might be a sensitive marker for prostatic stromal lesions [27]. As an exception, function of MUC18 (CD146) in cancer development took little attention such that only a group declared that it was a new biomarker that promoted angiogenesis in vitro and in vivo [28]. In line with this, we also found that MUC18 was overexpressed in BC cells with CSC and metastatic features. RAB3B which is an oncogene expressed on cell junctions was found to be upregulated in liver cancer stem cells according to the RNA seq analysis in a report in 2019 [29]. Additionally, RAB3B was reported as a new prognostic marker for extracellular vesicle biopsies of prostate cancer patients [30]. Similarly, we observed 24 folds increase in RAB3B gene expression level in BC for the first time. Based on these results, we also propose RAB3B as a biomarker in BCSC. Manupati et al. asserted that LPL was an effector of CD44 signaling and might be a potential target in BCSCs to inhibit progression [31]. We also found that LPL was overexpressed by dug resistant BC cells five folds more than drug sensitive BC cells. So, here we propound that LPL may be considered as target to combat BC in future studies. Mitogen activated protein kinases (MAPKs) are the members of protein kinases family. They phosphorylate their own serine and threonine residues or that of on their target proteins. MAPKs modulate the genes related with cell growth, metabolism, cell survival, stress response, apoptosis, migration, cellular differentiation and immune defense [32]. MAP kinase cascade has important role in cell proliferation and migration. In the present work, we displayed that while MAP3K12, MAP4K4, MAP4K1, MAPKKK4 were upregulated, MAP3K3, MAPKAPK2 and MAP3K1 were down regulated by development of drug resistance in BC cells. Activation of MAP kinases may exhibit anti apoptotic or pro apoptotic effect and may regulate tumor progression. Based on these discussions, MAP cascade seems to be modulated in favor of drug resistance in BC model. We found VMP1 gene expression to be significantly downregulated (circa 6.5 folds) in drug resistant BC cells with respect to sensitive MCF-7 model. It was declared that VMP1 was associated with regulation of autophagy and inhibition of VMP1 results in blockage of autophagosome formation [33]. In our case, down regulation of VMP1 may have resulted in resistance to cell death through autophagy. Hence, we may propose that there may be a cross regulation between autophagy and apoptosis in BCSC-like model cells.

In recent years, increasing knowledge on the mechanisms of endogenous RNA interference opened a door for the usage of siRNAs as innovative nucleic acid therapeutics for the treatment of cancer and multi-drug resistance. The research and development of siRNA drugs has taken 20 years and three siRNA drugs have been approved by the FDA since 2018 [34]. Furthermore, there has been intensive research on the treatment of cancer with siRNA-based drugs, which is quite promising. In this regard, after revealing potential marker genes, we applied siRNA mediated silencing method to downregulate MDR1 gene. In this context, we figured out the effect of MDR1 silencing on both proliferation of BCSC-like cells and cancer stem cell related genes. Eventually, we intended to propose a mechanism to reverse drug resistance in BCSC-like cells. RTCA results clearly demonstrated that silencing of the MDR1 gene and its consequences together caused two folds sensitization of BCSC model cells to paclitaxel treatment. Similarly, Li et al. successfully down-regulated MDR1 gene expression by about 50% in leukemia cell line models through transfection with 5-triphosphate modified-siMDR1 [35]. They also observed a decrease in cell viability in response to in vitro doxorubicin therapy. They proposed that down-regulation of MDR1 gene expression by siRNA might be a good strategy to combat drug resistance in myeloid leukemia [35]. Dönmez et al. also reported that siRNA silencing might be an effective MDR down-regulation strategy that restores the concentration of the drug inside resistant cells and increases the efficiency of doxorubicin in doxorubicin-resistant MCF-7 cell line [15]. Based on Western blot and flow cytometry analyses, we demonstrated that siMDR1 down regulated P-gp expression and also its function in BCSC-like cells. The decreased drug efflux most probably resulted in the decrease in cell viability after paclitaxel therapy.

Here, we also tested the effect of MDR1 down-regulation on the expression levels of important multiple proliferative genes which showed pronounced expressions in BCSCs. Strikingly, down-regulation of MDR1 gene resulted in down-regulation of ALDH1A, EGFR and BAG4 genes. So, our results indicated that inhibition of MDR1 led to the suppression of genes that play role in proliferation of BCSC-like cells. We can simply conclude that apoptosis was induced by restoring the concentration of paclitaxel within cells and reducing the expression of the anti-apoptotic BAG4 gene. In accordance with our results, Li et al. declared that MDR1 silencing resulted in increased apoptosis by releasing cytochrome C from mitochondria and increasing cleavage of caspase-3 and caspase-9 [35]. As drug resistance develops, cell membrane is also modulated; synthesis of sphingolipids, ceramide and cholesterol are elevated. These lipids are embedded in the membrane leading to decreased permeability (decreased influx). Increased expression of EGFR and CD44 membrane receptors was regarded as the outcomes of these changes in the membrane [36]. In parallel to this finding, it could be concluded that MDR1 gene down-regulation resulted in suppression of membrane changes which led to a decrease in EGFR expression. Changxing et al. showed that silencing of MDR1 gene by nano-carrier siRNA system and paclitaxel treatment exhibited synergistic anti-proliferative and gene suppression effects in colon cancer stem cell models [37], which is in concordance with our results. Li et al. recently published a gene correlation table using the Breast Cancer Gene-Expression Miner v4.5. The bioinformatics data consisted of gene expression, correlation and prognostic modules of 4908 breast cancer patients. The table and the additional results demonstrated that high expression levels of EGFR, CD44 and ALDH1A3 were correlated in triple negative breast cancer cells, known to be resistant to chemotherapy. Silencing of EGFR by inhibiting gamma amino butyric acid type A receptor in triple negative breast cancer patients resulted in down regulation of CD44, ALDH1A1 and hindrance of stemness of breast cancer. In parallel to our discussions, they also argued that inhibition of EGFR enhanced the sensitivity to chemotherapeutic agents [38].

To conclude, siRNA mediated silencing strategy could successfully be applied to target BCSC-like cells and may be proposed as a future chemotherapy strategy. To our knowledge, this is the first report demonstrating that silencing of the MDR1 gene results in down-regulation of survival/death-related genes such as EGFR, ALDH1A3 and BAG4. Considering that chemotherapy drugs targeting the MDR1 gene are in the research and development phase, the outcomes of this research will encourage and inspire those studies. As a future prospect, the CSC marker genes, revealed by microarray analysis, need and deserve to be further investigated for the purpose of targeting BCSCs.


Corresponding author: Meltem Demirel Kars, Faculty of Engineering, Department of Biomedical Engineering, Necmettin Erbakan University, Köyceğiz Campus, Konya, 42140, Turkey, Phone: +90 332 325 2024, E-mail:

Award Identifier / Grant number: 15201105

  1. Research funding: The project was funded by Selcuk University (project number: 15201105).

  2. Author contribution: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Competing interests: Authors state no conflict of interest.

  4. Informed consent: Not applicable.

  5. Ethical approval: Not applicable.

References

1. Krishan, A, Fitz, CM, Andritsch, I. Drug retention, efflux, and resistance in tumor cells. Cytometry 1997;29:279–85. https://doi.org/10.1002/(sici)1097-0320(19971201)29:4<279::aid-cyto3>3.0.co;2-5.10.1002/(SICI)1097-0320(19971201)29:4<279::AID-CYTO3>3.0.CO;2-5Search in Google Scholar

2. Kars, MD, Iseri, OD, Ufuk, G, Ural, AU, Arpacı, F, Molnar, J. Development of rational in vitro models for drug resistance in breast cancer and modulation of MDR by selected compounds. Anticancer Res 2006;26:4559–68.Search in Google Scholar

3. Rodriguez-Torres, M, Allan, AL. Aldehyde dehydrogenase as a marker and functional mediator of metastasis in solid tumors. Clin Exp Metastasis 2016;33:97–113. https://doi.org/10.1007/s10585-015-9755-9.Search in Google Scholar

4. Lo, HW, Hung, MC. Nuclear EGFR signalling network in cancers: linking EGFR pathway to cell cycle progression, nitric oxide pathway and patient survival. Br J Cancer 2006;94:184–8. https://doi.org/10.1038/sj.bjc.6602941.Search in Google Scholar

5. Steelman, LS, Chappell, WH, Akula, SM, Abrams, SL, Cocco, L, Manzoli, L, et al.. Therapeutic resistance in breast cancer cells can result from deregulated EGFR signaling. Adv Biol Regul 2020;78:100758. https://doi.org/10.1016/j.jbior.2020.100758.Search in Google Scholar

6. Takayama, S, Reed, JC. Molecular chaperone targeting and regulation by BAG family proteins. Nat Cell Biol 2001;3:237–41. https://doi.org/10.1038/ncb1001-e237.Search in Google Scholar

7. Mariotto, E, Viola, G, Zanon, C, Aveic, S. A BAG’s life: every connection matters in cancer. Pharmacol Ther 2020;209:107498. https://doi.org/10.1016/j.pharmthera.2020.107498.Search in Google Scholar

8. Enciso-Benavides, J, Alfaro, L, Castañeda-Altamirano, C, Rojas, N, González-Cabeza, J, Enciso, N, et al.. Biological characteristics of a sub-population of cancer stem cells from two triple-negative breast tumour cell lines. Heliyon 2021;7:e07273. https://doi.org/10.1016/j.heliyon.2021.e07273.Search in Google Scholar

9. Kars, MD. Method of isolation of chemotherapy resistant breast cancer stem cells. Turkish Patent Institute; TR: 2013;10837 B: 2020.Search in Google Scholar

10. Kurreck, J. Antisense technologies. Improvement through novel chemical modifications. Eur J Biochem 2003;270:1628–44. https://doi.org/10.1046/j.1432-1033.2003.03555.x.Search in Google Scholar

11. Gewirtz, AM. On future’s doorstep: RNA interference and the pharmacopeia of tomorrow. J Clin Invest 2007;117:3612–4. https://doi.org/10.1172/jci34274.Search in Google Scholar

12. Kars, MD, Iseri, OD, Gündüz, U, Molnar, J. Reversal of multidrug resistance by synthetic and natural compounds in drug-resistant MCF-7 cell lines. Chemotherapy 2008;54:194–200. https://doi.org/10.1159/000140462.Search in Google Scholar

13. Kars, MD, Yildirim, G. Determination of the target proteins in chemotherapy resistant breast cancer stem cell-like cells by protein array. Eur J Pharmacol 2019;848:23–9. https://doi.org/10.1016/j.ejphar.2019.01.052.Search in Google Scholar PubMed

14. Iseri, OD, Kars, M, Arpacı, F, Atalay, C, Pak, I, Gündüz, U. Drug resistant MCF-7 cells exhibit epithelial-mesenchymal transition gene expression pattern. Biomed Pharmacother 2011;65:40–5.10.1016/j.biopha.2010.10.004Search in Google Scholar PubMed

15. Dönmez, Y, Gündüz, U. Reversal of multidrug resistance by small interfering RNA (siRNA) in doxorubicin-resistant MCF-7 breast cancer cells. Biomed Pharmacother 2011;65:85–9. https://doi.org/10.1016/j.biopha.2010.12.007.Search in Google Scholar PubMed

16. 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. https://doi.org/10.1006/meth.2001.1262.Search in Google Scholar PubMed

17. Kars, MD, Iseri, OD, Gündüz, U. A microarray based expression profiling of paclitaxel and vincristine resistant MCF-7 cells. Eur J Pharmacol 2011;657:4–9. https://doi.org/10.1016/j.ejphar.2011.02.001.Search in Google Scholar PubMed

18. Hu, Y, Yague, E, Zhao, J, Wang, L, Bai, J, Yang, Q. Sabutoclax, pan-active BCL-2 protein family antagonist, overcomes drug resistance and eliminates cancer stem cells in breast cancer. Cancer Lett 2018;423:47–59. https://doi.org/10.1016/j.canlet.2018.02.036.Search in Google Scholar PubMed

19. Yang, LT, Leong, Y, Song, WS, Yang, H, Wang, W, Wang, M. Efficient and tumor-specific knockdown of MTDH gene attenuates paclitaxel resistance of breast cancer cells both in vivo and in vitro. Breast Cancer Res 2018;20:113. https://doi.org/10.1186/s13058-018-1042-7.Search in Google Scholar PubMed PubMed Central

20. Li, Y, Gao, X, Yu, Z, Liu, B, Pan, W, Li, N. Reversing multidrug resistance by multiplexed gene silencing for enhanced breast cancer chemotherapy. ACS Appl Mater Interfaces 2018;10:15461–6. https://doi.org/10.1021/acsami.8b02800.Search in Google Scholar PubMed

21. Phi, LTH, Sari, IN, Yang, YG, Lee, SH, Jun, N, Kim, KS. Cancer stem cells (CSCs) in drug resistance and their therapeutic implications in cancer treatment. Stem Cell Int 2018;2018:5416923. https://doi.org/10.1155/2018/5416923.Search in Google Scholar PubMed PubMed Central

22. Wang, N, Weng, J, Xia, J, Zhu, Y, Chen, Q, Hu, D, et al.. SIPA1 enhances SMAD2/3 expression to maintain stem cell features in breast cancer cells. Stem Cell Res 2020;49:102099. https://doi.org/10.1016/j.scr.2020.102099.Search in Google Scholar PubMed

23. Mukherjee, N, Panda, CK. Wnt/β-Catenin signaling pathway as chemotherapeutic target in breast cancer: an update on pros and cons. Clin Breast Cancer 2020;20:361–70. https://doi.org/10.1016/j.clbc.2020.04.004.Search in Google Scholar PubMed

24. Zhang, JY, Luo, Q, Xu, JR, Bai, J, Mu, LM, Yan, Y, et al.. Regulating stem cell-related genes induces the plastic differentiation of cancer stem cells to treat breast cancer. Mol Ther Oncolytics 2020;18:396–408. https://doi.org/10.1016/j.omto.2020.07.009.Search in Google Scholar PubMed PubMed Central

25. Sun, Y, Zhang, J, Ma, L. α-catenin. A tumor suppressor beyond adherens junctions. Cell Cycle 2014;13:2334–9. https://doi.org/10.4161/cc.29765.Search in Google Scholar PubMed PubMed Central

26. Hebenstreit, D, Wirnsberger, G, Horejs-Hoeck, J, Duschl, A. Signaling mechanisms, interaction partners, and target genes of STAT6. Cytokine Growth Factor Rev 2006;17:173–88. https://doi.org/10.1016/j.cytogfr.2006.01.004.Search in Google Scholar PubMed

27. Guner, G, Bishop, JA, Bezerra, SM, Taheri, D, Zahavi, DJ, Mendoza Rodriguez, MA, et al.. The utility of STAT6 and ALDH1 expression in the differential diagnosis of solitary fibrous tumor versus prostate-specific stromal neoplasms. Hum Pathol 2016;54:184–8. https://doi.org/10.1016/j.humpath.2016.03.011.Search in Google Scholar PubMed

28. Jiang, T, Zhuang, J, Duan, H, Luo, Y, Zeng, Q, Fan, K, et al.. CD146 is a coreceptor for VEGFR-2 in tumor angiogenesis. Blood 2012;120:2330–9. https://doi.org/10.1182/blood-2012-01-406108.Search in Google Scholar PubMed

29. Tsunedomi, R, Yoshimura, K, Kimura, Y, Nishiyama, M, Matsukuma, S, Tokumitsu, Y, et al.. Cancer stem-like phenotypes including immune surveillance and its responsible genes in induced liver cancer stem-like cells. Ann Oncol 2019;30:ix59–60. https://doi.org/10.1093/annonc/mdz422.050.Search in Google Scholar

30. Weng, J, Xiang, X, Ding, L, Wong, AL, Zeng, Q, Sethi, G, et al.. Extracellular vesicles, the cornerstone of next-generation cancer diagnosis? Semin Cancer Biol 2021;11:139–5.10.1016/j.semcancer.2021.05.011Search in Google Scholar PubMed

31. Manupati, K, Yeeravalli, R, Kaushik, K, Singh, D, Mehra, B, Gangane, N, et al.. Activation of CD44-Lipoprotein lipase axis in breast cancer stem cells promotes tumorigenesis. Biochim Biophys Acta (BBA) - Mol Basis Dis 2021;1867:166228. https://doi.org/10.1016/j.bbadis.2021.166228.Search in Google Scholar PubMed

32. Mohan, CD, Liew, YY, Jung, YY, Rangappa, S, Preetham, HD, Chinnathambi, A, et al.. Brucein D modulates MAPK signaling cascade to exert multi-faceted anti-neoplastic actions against breast cancer cells. Biochimie 2021;182:140–51. https://doi.org/10.1016/j.biochi.2021.01.009.Search in Google Scholar PubMed

33. Ishaq, M, Ojha, R, Sharma, AP, Singh, SK. Autophagy in cancer: recent advances and future directions. Semin Cancer Biol 2020;66:171–81. https://doi.org/10.1016/j.semcancer.2020.03.010.Search in Google Scholar PubMed

34. Kim, YK. RNA therapy: current status and future potential. Chonnam Med J 2020;56:87–93. https://doi.org/10.4068/cmj.2020.56.2.87.Search in Google Scholar PubMed PubMed Central

35. Li, D, Gale, RP, Liu, Y, Lei, B, Wang, Y, Diao, D. 5’-Triphosphate siRNA targeting MDR1 reverses multi-drug resistance and activates RIG-I-induced immune-stimulatory and apoptotic effects against human myeloid leukaemia cells. Leuk Res 2018;58:23–30. https://doi.org/10.1016/j.leukres.2017.03.010.Search in Google Scholar PubMed

36. Singh, MS, Tammam, SN, Shetab, BMA, Lamprecht, A. MDR in cancer: addressing the underlying cellular alterations with the use of nanocarriers. Pharmacol Res 2017;126:2–30. https://doi.org/10.1016/j.phrs.2017.07.023.Search in Google Scholar PubMed

37. Changxing, L, Zhao, G, Liu, J, Ma, N, Chivukula, P, Perelman, L. Novel biodegradable lipid nano complex for siRNA delivery significantly improving the chemosensitivity of human colon cancer stem cells to paclitaxel. J Contr Release 2009;140:277–83. https://doi.org/10.1016/j.jconrel.2009.08.013.Search in Google Scholar PubMed

38. Li, X, Wang, H, Yang, X, Wang, X, Zhao, L, Zou, L, et al.. GABRP sustains the stemness of triple-negative breast cancer cells through EGFR signaling. Cancer Lett 2021;514:90–102. https://doi.org/10.1016/j.canlet.2021.04.028.Search in Google Scholar PubMed


Supplementary Material

The online version of this article offers supplementary material (https://doi.org/10.1515/tjb-2021-0275).


Received: 2021-09-15
Revised: 2022-01-10
Accepted: 2022-02-24
Published Online: 2022-03-28

© 2022 the author(s), published by De Gruyter, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Frontmatter
  2. Review Article
  3. Planning education for best practices in point-of-care testing
  4. Research Articles
  5. Effects of salt concentration on the production of cytotoxic geodin from marine-derived fungus Aspergillus sp.
  6. Is Vision C interchangeable with the modified Westergren method for the erythrocyte sedimentation rate?
  7. KIM-1 and GADDI-153 gene expression in paracetamol-induced acute kidney injury: effects of N-acetylcysteine, N-acetylmethionine, and N-acetylglucosamine
  8. MAP3K1 SNP rs889312 potential risk and MAP3K9 SNP rs11628333 menopause dependent association for breast cancer
  9. Downregulation of MMP-2 and MMP-9 genes in obesity patients and their relation with obesity-related phenotypes
  10. miR-221, miR-650 and miR-4534 as diagnostic markers in prostate cancer and their relationship with lymphatic invasion
  11. Biomarkers to target and silence stemness of breast cancer stem cell model: silencing MDR1 by siRNA
  12. The correlation between bone biomarkers, glucosylsphingosine levels, and molecular findings in Gaucher type 1 patients under enzyme therapy
  13. Endoplasmic reticulum aminopeptidase-1 polymorphism increases the risk of rheumatoid arthritis
  14. Evaluation of the phenolic compounds and the antioxidant potentials of Vitex agnus-castus L. leaves and fruits
  15. Antioxidant and antimicrobial potential of Ganoderma lucidum and Trametes versicolor
  16. Ginsenoside-Mc1 reduces ischemia/reperfusion-induced cardiac arrhythmias through activating JAK2/STAT3 pathway and attenuating oxidative/endoplasmic reticulum stress in hyperlipidemic rats
  17. Different spacer-arm attached magnetic nanoparticles for covalent immobilization of Jack bean urease
  18. Evaluation of the antifungal activity of essential oils against Alternaria alternata causing fruit rot of Eriobotrya japonica
  19. Emerging insights into the relationship between pre-microRNA146a rs2910164 gene polymorphism and TNF-α in ischemic stroke
  20. Education Section
  21. The effect of virtual laboratory simulations on medical laboratory techniques students’ knowledge and vocational laboratory education
  22. Medical students’ opinions on career planning course: evaluations of the relationship between course and faculty attributes and student characteristics
Downloaded on 25.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/tjb-2021-0275/html
Scroll to top button