- Open Access
Transcriptional factor six2 promotes the competitive endogenous RNA network between CYP4Z1 and pseudogene CYP4Z2P responsible for maintaining the stemness of breast cancer cells
Journal of Hematology & Oncologyvolume 12, Article number: 23 (2019)
The expression of CYP4Z1 and the pseudogene CYP4Z2P has been shown to be specifically increased in breast cancer by our group and others. Additionally, we previously revealed the roles of the competitive endogenous RNA (ceRNA) network mediated by these genes (ceRNET_CC) in breast cancer angiogenesis, apoptosis, and tamoxifen resistance. However, the roles of ceRNET_CC in regulating the stemness of breast cancer cells and the mechanisms through which ceRNET_CC is regulated remain unclear.
Transcriptional factor six2, CYP4Z1-3′UTR, and CYP4Z2P-3′UTR were stably overexpressed or knocked down in breast cancer cells via lentivirus infection. ChIP-sequencing and RNA-sequencing analysis were performed to reveal the mechanism through which ceRNET_CC is regulated and the transcriptome change mediated by ceRNET_CC. Clinical samples were used to validate the correlation between six2 and ceRNET_CC. Finally, the effects of the six2/ceRNET_CC axis on the stemness of breast cancer cells and chemotherapy sensitivity were evaluated by in vitro and in vivo experiments.
We revealed that ceRNET_CC promoted the stemness of breast cancer cells. Mechanistically, six2 activated ceRNET_CC by directly binding to their promoters, thus activating the downstream PI3K/Akt and ERK1/2 pathways. Finally, we demonstrated that the six2/ceRNET_CC axis was involved in chemoresistance.
Our results uncover the mechanism through which ceRNET_CC is regulated, identify novel roles for the six2/ceRNET_CC axis in regulating the stemness of breast cancer cells, and propose the possibility of targeting the six2/ceRNET_CC axis to inhibit breast cancer stem cell (CSC) traits.
Human tumors are composed of heterogeneous cells, in which cancer stem cells (CSCs), here defined as the tumor cells specifically endowed with self-renewal and tumor-seeding potential, have been regarded as the main drivers of tumor progression and chemoresistance ; however, the mechanisms underlying CSC maintenance are not well defined, and currently, no drugs directly kill CSCs. In fact, we have only a scattered understanding of the cellular mechanisms that contribute to CSC attributes. Notably, because tumorigenesis and chemoresistance are led by the deregulation of gene networks, single gene expression analysis cannot completely explain chemotherapy or tumor recurrence, which should be assessed by several synergistic factors.
Competitive endogenous RNAs (ceRNAs) are defined as the transcripts that share common miRNAs and can coregulate each other’s expression . ceRNAs play critical roles in tumor progression; for example, the long non-coding RNA UICLM can facilitate colorectal cancer metastasis by acting as a ceRNA for miR-215  and ROR promotes pancreatic cancer progression by acting as a ceRNA for Nanog . We and others have shown that the expression of CYP4Z1 and the pseudogene CYP4Z2P are significantly and specifically higher in breast cancer tumor tissues [5,6,7]. Furthermore, we revealed that CYP4Z1 and the pseudogene CYP4Z2P form a ceRNA network via competiting with several miRNAs, such as miR-211, miR-125a-3p, miR-197, miR-1226, and miR-204, here called ceRNET_CC, and both we and others have shown that these shared miRNAs exert tumor suppressive effects [8,9,10,11,12]. We further demonstrate that ceRNET_CC promotes breast cancer angiogenesis  and tamoxifen resistance  and suppresses breast cancer apoptosis ; these effects suggest that suppression of ceRNET_CC may allow for inhibition of breast cancer. Due to the roles of ceRNET_CC in promoting breast cancer, here, we focused on understanding the mechanisms underlying the progression of ceRNET_CC and the roles of ceRNET_CC in regulating the stemness of breast cancer. Transcript overexpression occurs through transcriptional and epigenetic control mechanisms in the vast majority of cancers, as demonstrated by Rinath et al., who uncovered a role for the RUNX2-ER (runt-related transcription factor 2-estrogen receptor) complex in stimulating the transcription of a set of genes, including most notably the stem cell factor Sox9, which promotes proliferation and a metastatic phenotype , and TRIM28, which interacts with EZH2 and SWI/SNF to activate genes that promote mammosphere formation . A bioinformatics method was then used to predict the transcriptional factors binding to the promoters of CYP4Z1 and the pseudogene CYP4Z2P, and transcriptional factor six2 attracted our attention based on its critical roles in organ development and promoting roles in cancers [16,17,18].
Transcriptional factor six2 activity is required to maintain the mesenchymal progenitor population in an undifferentiated state , and cells autonomously regulate a multipotent nephron progenitor population throughout mammalian kidney development . Additionally, haploinsufficiency for the six2 gene increases nephron progenitor proliferation, promoting branching and nephron number . These processes are normally maintained by stem cells, and the expression of genes involved in embryonic development is usually reinstated in tumors. Additionally, previous studies have revealed that several cell fate regulators (DeltaNp63, Slug, Sox9, and miR-200c) are molecular links between mammary stem cells and breast tumor-initiating cells that drive renewal activity in both normal and cancerous mammary gland tissues [20,21,22]. Based on these data, we speculated that six2 could drive CSC progression. Previous studies have shown that breast cancer patients with higher six2 levels have shorter time to both relapse and metastasis [18, 23]. Additionally, increased expression and decreased methylation of six2 are correlated with increased tumor size, clinical stage, vascular invasion, and unfavorable histological differentiation in Wilms’ tumor . Moreover, a systematic meta-analysis revealed that higher six2 expression is associated with a greater possibility of tumorigenesis and predicted poor overall survival (OS) in non-small-cell lung cancer (NSCLC) and poor relapse-free survival (RFS) in lung adenocarcinoma (ADC) . However, the molecular mechanisms underlying six2-mediated oncogenic effects, and the downstream effectors of six2 are not fully understood.
In the present study, we showed that the expression of six2, CYP4Z1, and pseudogene CYP4Z2P was significantly increased in breast cancer tumors and that six2 could directly bind to the 5′-TCAG-3′ motif in the promoter of CYP4Z1 and CYP4Z2P and thus promote the progression of ceRNET_CC. Notably, this novel six2/ceRNET_CC regulatory axis was responsible for the stemness and chemoresistance of breast cancer cells. Importantly, the expression of six2, CYP4Z1, and the pseudogene CYP4Z2P was negatively correlated with the OS of breast cancer patients, and the expression of these genes was positively correlated with one another, underscoring the critical roles of this regulatory axis in breast cancer progression. A broader understanding of six2-dependent regulation on ceRNET_CC is needed to effectively target therapy-resistant breast cancer cells with stemness characteristics and address the challenges of tumor heterogeneity.
Materials and methods
Paraffin-embedded breast cancer tissue samples were obtained in our previous work . Eight pairs of fresh breast cancer and normal adjacent tissues were collected from Huai An First People’s Hospital and Jiangsu Province Hospital of TCM between September 2017 and January 2018. Written informed consent from all patients and approval of the hospital ethics review committees were obtained. Online clinical posited data (http://www.firebrowse.org/), Breast Cancer Gene-Expression Miner (version 4.0; http://bcgenex.centregauducheau.fr/), and R2: Genomics Analysis and Visualization Platform (http://hgserver1.amc.nl/cgi-bin/r2/main.cgi) were used for the analyses of gene expression and correlation. The diagnostic values of the genes were analyzed by KM plotter (http://kmplot.com) to obtain KM survival plots in which the number at risk is indicated below the main plot . Hazard ratio, 95% confidence intervals, and log-rank P values were calculated and displayed on the webpage.
Cell culture and chemical reagents
The human breast cancer cell lines MCF-7, MDA-MB-231, and HEK293T were preserved in our laboratory. Adriamycin-resistant MCF-7-Adr cells were purchased from KeyGen BioTECH (Nanjing, China). The cell line was authenticated every year through short tandem repeat (STR) DNA profiling. HEK293T and MCF-7 cells were cultured in DMEM (Gibco, Grand Island, NY, USA), MCF-7-Adr cells were cultured in 1640 medium (Gibco), and MDA-MB-231 cells were cultured in L-15 medium (Gibco) at 37 °C under a humidified atmosphere with 5% CO2. All of the media were supplemented with 10% FBS (Gibco), 80 U/ml penicillin, and 0.08 mg/ml streptomycin. PI3K inhibitor (LY-294002) and ERK1/2 inhibitor (VX-11e) were purchased from APExBIO. Adriamycin was purchased from Zhongda Hospital Southeast University.
Quantitative real-time PCR (qRT-PCR)
Total RNA from the cells was extracted using TransZol Up (Cat. No. ET111-01, TransGen Biotech, Beijing, China) following the manufacturer’s recommendation. Total RNA from paraffin-embedded breast cancer tissues was extracted using a total RNA extraction kit for paraffin-embedded tissues (Cat. No. DP439, TianGen Biotech, Beijing, China) according to standard protocols. Then, complementary DNA (cDNA) was reverse-transcribed using M-MLV (H-) Reverse Transcriptase (Cat. No. R021-01, Vazyme, Nanjing, China) according to the manufacturer’s protocol. qRT-PCR was performed with AceQ Universal SYBR qPCR Master Mix (Cat. No. Q511-02, Vazyme). A melting curve analysis was performed routinely to check the amplification specificity. cDNA templates were analyzed in triplicate, and GAPDH was used as an internal control. The relative expression level of each transcript was calculated by the 2-△△ct method. The qRT-PCR primers are described in Additional file 1: Table S1.
The detailed procedure was described in our previous study . Protein in fresh tissues was extracted using total protein extraction kit (Invent, USA) following the manufacturer’s recommendation. β-actin or GAPDH was used as an internal reference. Detailed information on the antibodies used in this work is given in Additional file 2: Table S2.
Fluorescence-activated cell sorting
CD24 and CD44 expression was analyzed in cells derived from monolayer cultures following dissociation in trypsin-EDTA at 37 °C. At least 1 × 106 cells were pelleted by centrifugation at 300×g and 4 °C for 5 min. Then, cells were washed in PBS, re-suspended with anti-CD24-PE (BD Biosciences, USA) and anti-CD44-APC (BD Biosciences, USA), and then incubated at 4 °C for 30 min in the dark. The labeled cells were washed using PBS and analyzed using a flow cytometer (BD, USA). The negative fraction was determined using appropriate isotype controls.
Chromatin immunoprecipitation assay
A chromatin immunoprecipitation (ChIP) assay was performed using the EZ-Magna ChIP™ A/G Chromatin Immunoprecipitation Kit (Cat. No. 17–10086, Merck) following the manufacturer’s protocols. Primers flanking the six2 binding sites on the promoters of CYP4Z1 and pseudogene CYP4Z2P were used for qRT-PCR. The sequences of the primers for ChIP analysis were denoted in Additional file 3: Table S3.
ChIP-sequencing and data assay
ChIP-sequencing analysis was performed by GENEWIZ (Suzhou, China). ChIP-seq raw reads were aligned to a human reference genome (hg19) using cutadapt (version 1.9.1) to pass filter data and acquire clean data. Up to mismatch per read was allowed. The quality of the sequencing data was assessed using FastQC (v0.10.1), and only uniquely mapped reads were kept for downstream analysis. The data are available in the Gene Expression Omnibus (GEO) database as GSE117145.
RNA sequencing and data analysis
RNA sequencing and data analysis were conducted by Novogene (Beijing, China). The data are available in the Gene Expression Omnibus (GEO) database as GSE116984.
Tissue microarray analysis
A tissue microarray including 30 breast cancer tissues and 30 normal adjacent tissues was purchased from OUTDO IVD (Shanghai, China). Further immunohistochemistry was performed to detect the expression of six2 following the protocols described in our recent work .
Small interfering RNA transfection
When cells confluency reached 50%, cells were transfected with a final concentration of 50 nM small interfering RNA (siRNA) or 50 nM NC (si-NC, inhibitor-NC, mi-NC) which were synthesized by Biomics Biotechnology (Nantong, China) using jetPRIME (Polyplus Transfection, France) following the manufacturer’s protocols. The siRNA sequences are mentioned in Additional file 4: Table S4.
Cell spheroid formation assay
A mammosphere formation assay was performed using MammoCult™ Human Medium Kit (STEMCELL Technologies, Canada). A total of 3 × 103 cells were mixed with 500 μl of Complete MammoCult™ Medium in the presence of 4 μg/ml Heparin (STEMCELL Technologies, Canada) and 0.48 μg/ml hydrocortisone (STEMCELL Technologies, Canada) and seeded in 24-well ultra-low attachment plates (Corning, USA) for 7 days. The spheres were counted and photographed. All images were obtained with a Leica DMI microscope (DE).
Lentivirus package, stable cell lines, and plasmid construction
To construct stable expression cells, the six2 coding area or the CYP4Z1 3′UTR or CYP4Z2P 3′UTR sequences were subcloned into pLVX-ZsGreen. siRNA oligos were purchased from GenePharma (Shanghai, China). After annealing, double-strand oligos were inserted into the lentiviral pLKO.1-puro vector (Addgene). To package lentivirus, HEK293T cells were co-transfected using Lentifectin (ABM, USA) with the lentiviral vector and packaging vectors psPAX2 and pMD2.G. MCF-7 and MDA-MB-231 cells were infected with the virus in the presence of 2 μg/ml polybrene. Cells infected with Plko.1-derived vectors were selected with puromycin (Sigma, 2 μg/ml) for 2 weeks. Cells infected with pLVX-ZsGreen-derived vectors were selected by fluorescent cell sorting. Western blot and qRT-PCR analyses were used to verify expression levels. MCF-7 cells stably overexpressing CYP4Z1 3′UTR, CYP4Z2P 3′UTR, and six2 were designated as MCF-7-Z1-UTR, MCF-7-Z2P-UTR, and MCF-7-six2, respectively. MCF-7 cells stably knocking down CYP4Z1 3′UTR, CYP4Z2P 3′UTR, and six2 stable were denoted as MCF-7-Plko-Z1, MCF-7-Plko-Z2P, and MCF-7-Plko-six2, respectively. MDA-MB-231 cells stably overexpressing six2 were designated as 231-six2, while those stably knocking down CYP4Z1 3′UTR, CYP4Z2P 3′UTR, and six2 were denoted as 231-Plko-Z1, 231-Plko-Z2P, and 231-Plko-six2, respectively. The sequences of the primers used for plasmid constructions were listed in Additional file 5: Table S5.
Luciferase reporter assay
The promoter sequences of CYP4Z1 (1–2403 bp) and the pseudogene CYP4Z2P (1–2998 bp) were inserted into the pGL3-promoter vector, and truncations of these sequences were also inserted into the pGL3-promoter vector. Then, potential six2 binding sites were mutated using the Fast Mutagenesis Kit V2 (Vazyme, China) following the manufacturer’s instructions and inserted into the pGL3-promoter vector as well. All of the abovementioned constructs were verified by DNA sequencing before use. To analyze the activity of the abovementioned constructs, they were individually co-transfected with β-gal and a six2 overexpression vector into MCF-7 cells. Seventy-two hours later, the luciferase activity was measured by a POLARstar Omega multimode microplate reader according to the manufacturer’s protocol and normalized to β-gal. All primer sequences for this experiment are listed in Additional file 6: Table S6.
Chemoresistance or sensitivity assay
Cell viability was assessed using MTT (Amresco, USA) staining. Cells were seeded at 3 × 103 cells/well in 96-well cell plates overnight and then treated with different concentrations of adriamycin for 48 h. During the last 3.5 h, the cells were exposed to MTT (5 mg/ml) and the resulting formazan crystals were dissolved in 150 μl of DMSO and measured using a spectrophotometer(BIO-RAD) at a test wavelength of 490 nM. The experiments were conducted in triplicate.
In vivo tumorigenesis
The procedure was referenced in our previous work . Briefly, six-week-old male athymic BALB/c nude mice were purchased from the Model Animal Research Center of Nanjing University, housed and fed under standard pathogen-free conditions. For the tumor-limiting dilution assay, 1 × 106 and 1 × 105 MCF-7 cells, 1 × 105 and 1 × 104 MDA-MB-231 cells, or 1 × 105 MCF-Adr cells receiving different treatment were orthotopically implanted in the inguinal mammary gland of mice. On day 8, all the mice were euthanized and tumor tissues were collected and weighed. The stem cell frequencies were calculated using an ELDA (http://bioinf.wehi.edu.au/software/elda/) . All the animal experiments were carried out with the approval of the Ethics Committee for Animal Experimentation of China Pharmaceutical University.
GraphPad Prism 5.01 software (GraphPad Software, Inc., La Jolla, CA, USA) was used for the statistical analysis. All the data were obtained from at least three independent experiments and are presented as the means ± standard deviations (SDs). Datasets with only two groups were analyzed using Student’s t test. Differences between multiple groups were analyzed by one-way analysis of variance with the Tukey-Kramer post hoc test. P < 0.05 was considered indicative of a statistically significant difference.
CeRNET_CC promotes the stemness of breast cancer cells in vitro
The expression levels of CYP4Z1 and its pseudogene CYP4Z2P were initially examined in breast tumor and normal adjacent tissues via online clinical deposited data (http://www.firebrowse.org/). As shown in Fig. 1a, b, the expression of CYP4Z1 and the pseudogene CYP4Z2P was significantly increased in breast tumor tissues, which is consistent with our previous work in which CYP4Z1 and pseudogene CYP4Z2P expression was detected in clinical samples [26, 30]. KM-plotter analysis (http://kmplot.com) indicated that CYP4Z1 expression was negatively correlated with the OS of breast cancer patients (Fig. 1c). We defined the CYP4Z1-3′UTR- or CYP4Z2P-3′UTR-regulated transcriptome in MCF-7 cells with CYP4Z1- or CYP4Z2P-3′UTR overexpression and subsequently performed expression profiling, and the results revealed a substantial and highly statistically significant overlap between genes regulated by both CYP4Z1-3′UTR and CYP4Z2P-3′UTR in MCF-7 cells (Fig. 1d). Additionally, CYP4Z1-3′UTR and CYP4Z2P-3′UTR activated (3854 vs. 3788) or repressed (3404 vs. 3934) similar numbers of genes in MCF-7 cells (Fig. 1e, f). Further gene set enrichment analysis (GSEA) of this dataset revealed that a negative enrichment of stem cell-differentiated signatures was observed in MCF-7 cells with CYP4Z2P-3′UTR overexpression (Fig. 1g), and the top positive signatures associated with CYP4Z1-3'UTR overexpression were the embryonic stem cell function and adult tissue stem modules (Fig. 1h, i) , both of which are indicative of a stemness phenotype. Functional annotation analysis revealed that the overexpression of CYP4Z2P-3′UTR or CYP4Z1-3′UTR activated signaling pathways regulating the pluripotency of stem cells (Fig. 1j). Among the differentially expressed genes, we identified a set of gene signatures related to epithelial cancer stem cells in MCF-7 cells with CYP4Z2P-3′UTR and CYP4Z1-3′UTR overexpression, including YY1, KRAS, YAP1, and HMGB2 (Fig. 1k) . Consistent with the activation of stem cell function, we detected an increase in a set of cell cycle-related genes, namely CDK5R1, CDK1, CDK2, and CDK7 (Fig. 1). Notably, the expression of CYP4Z1 and CYP4Z2P displayed a positive correlation in clinical breast cancer samples (P < 0.001, Fig. 1l).
As previous studies have demonstrated that non-adherent spheres are highly enriched for CSCs [32, 33], the expression levels of CYP4Z1 and its pseudogene CYP4Z2P were examined in non-adherent MCF-7 spheres and parental cells. qRT-PCR results showed that non-adherent MCF-7 spheres displayed significantly higher levels of CYP4Z1 and CYP4Z2P relative to monolayer cultures of MCF-7 cells (Additional file 7: Figure S1A). We then determined whether the overexpression of CYP4Z1- or CYP4Z2P-3′UTR conferred stemness upon breast cancer cells in vitro. First, MCF-7 cells with CYP4Z1- or CYP4Z2P-3′UTR stable overexpression or knockdown were subjected to a spheroid formation assay. The infection efficiency was confirmed by qRT-PCR (Additional file 7: Figure S1B and C). Both the sphere size and number were increased in MCF-7 cells with CYP4Z1- or CYP4Z2P-3′UTR overexpression, whereas the knockdown of CYP4Z1- or CYP4Z2P-3′UTR exerted the opposite effects (Additional file 7: Figure S1D). Additionally, the overexpression of CYP4Z1- or CYP4Z2P-3′UTR increased the CD44+/CD24− population, which has been identified as having breast CSC markers  (Additional file 7: Figure S1E). Moreover, the expression of several pluripotent transcription factors, namely Oct3/4, ALDH1, Nanog, and Sox2, was increased in cells with CYP4Z1- or CYP4Z2P-3′UTR overexpression, whereas the knockdown of CYP4Z1- or CYP4Z2P-3′UTR yielded the opposite effects (Additional file 7: Figure S1F and S1G). These results were recapitulated in MDA-MB-231 cells with or without the knockdown of CYP4Z1- or CYP4Z2P-3′UTR. qRT-PCR analysis confirmed the knockdown efficiency of shRNAs against CYP4Z1- or CYP4Z2P-3′UTR (Additional file 8: Figure S2A and B). As expected, the knockdown of CYP4Z1- or CYP4Z2P-3′UTR resulted in fewer primary mammospheres than control cells (Additional file 8: Figure S2C) and decreased the CD44+/CD24− population (Additional file 8: Figure S2D). Additionally, the expression of stemness markers was inhibited in MDA-MB-231 cells with CYP4Z1- or CYP4Z2P-3′UTR knockdown (Additional file 8: Figure S2E and F). Notably, we assessed the expression correlation between Nanog and CYP4Z1 or CYP4Z2P across basal-like breast cancer subtypes using Breast Cancer Gene-Expression Miner (version 4.0; http://bcgenex.centregauducheau.fr/) and found that Nanog expression is positively correlated with CYP4Z1 or CYP4Z2P expression in basal-like breast cancer subtypes (Additional file 8: Figure S2G and S2H). These results indicate that ceRNET_CC is engaged in tumor stemness in breast cancer.
ceRNET_CC promotes the tumor-initiating potential of breast cancer cells in vivo
We further investigated whether ceRNET_CC facilitates the tumor-initiating potential of breast cancer cells in vivo. We compared the capacity of MCF-7 cells with CYP4Z1- or CYP4Z2P-3′UTR overexpression or knockdown to seed tumors at limiting dilutions. Although all of the cell lines could form tumors at a density of 1 × 106 and 1 × 105 cells, CYP4Z1- or CYP4Z2P-3′UTR overexpressed cells showed increased tumor size and weight (Fig. 2a–d) while knockdown of CYP4Z1- or CYP4Z2P-3′UTR decreased the tumor-initiating ability of MCF-7 cells (Fig. 2a–f). Additionally, we performed an in vivo tumorigenic assay with MDA-MB-231 cells after CYP4Z1- or CYP4Z2P-3′UTR knockdown and showed that the knockdown of CYP4Z1- or CYP4Z2P-3′UTR remarkably reduced the tumor-initiating potential of MDA-MB-231 cells (Fig. 2g–k). Therefore, our results demonstrate that ceRNET_CC could promote the stemness of breast cancer cells.
ceRNET_CC promotes the stemness of breast cancer cells partly through the hTERT/PI3K/Akt and ERK1/2 pathways
We next sought to explore the mechanisms through which ceRNET_CC promotes the stemness of breast cancer cells. To address this issue, we characterized pathways regulated by CYP4Z1- or CYP4Z2P-3′UTR based on RNA-sequencing data. As expected, we found that the phosphatidylinositol signaling system and MAPK signaling pathways in cancer were some of the most highly upregulated pathways after CYP4Z1- or CYP4Z2P-3′UTR overexpression (Fig. 3a, b). These results were supported by our previous studies showing that ceRNET_CC could act as a subceRNA network for hTERT and promote hTERT expression, thus activating hTERT/PI3K/Akt and ERK1/2 signaling [13, 26]. We then aimed to detect whether ceRNET_CC indeed facilitated the stemness of breast cancer cells through these two pathways. MCF-7-Z1-UTR or MCF-7-Z2P-UTR cells were transfected with siRNA against hTERT (sihTERT) or treated with a PI3K/Akt inhibitor (LY-294002) or ERK1/2 inhibitor (VX-11e), and we then detected the formation of cell spheroids and the CD44+/CD24− population. As shown in Fig. 3c, d, the increased cell spheroid formation or CD44+/CD24− population induced by CYP4Z1- or CYP4Z2P-3′UTR overexpression was attenuated by sihTERT transfection or by LY-294002 or VX-11e treatment. Additionally, the increased expression of stemness markers induced by CYP4Z1- or CYP4Z2P-3′UTR overexpression was attenuated by sihTERT (Fig. 3e) or by LY-294002 or VX-11e treatment (Fig. 3f) and the knockdown of hTERT or treatment with LY-294002 or VX-11e decreased the p-Akt and p-ERK1/2 levels. Thus, these results demonstrate that ceRNET_CC promotes the stemness of breast cancer cells in a manner dependent on the hTERT/PI3K/Akt and ERK1/2 pathways.
Transcriptional factor six2 induces the progression of ceRNET_CC by directly regulating the transcription of CYP4Z1 and the pseudogene CYP4Z2P
To determine the mechanisms through which ceRNET_CC is regulated, the Genomatix Software Suite v3.10 (https://www.genomatix.de) was used to predict transcriptional factors that could bind to the promoters of CYP4Z1 and CYP4Z2P. We aimed to identify specific transcriptional factors that correlate with tumor progression and tissue development, and six2 attracted our attention because it has been shown to promote breast cancer metastasis  and regulate the expansion of the nephron progenitor pool during nephrogenesis, which involves a process similar to CSC formation . We identified one (− 1458 nt) and two putative binding sites (− 1125 nt and − 2716 nt) on the promoters of CYP4Z1 and CYP4Z2P, respectively (Fig. 4a). To confirm the interaction between six2 and CYP4Z1 or CYP4Z2P, MCF-7 cells showing stable overexpression of six2 were constructed by lentivirus infection. We then performed a genome-wide chromatin immunoprecipitation sequencing (ChIP-seq) analysis to identify the six2-bound chromatin regions using an antibody against six2. As shown in Fig. 4b, binding molecular function and developmental process were enriched in six2-overexpressing cells. Importantly, signaling pathways regulating stem cell pluripotency were enriched in six2-overexpressed cells (Fig. 4c). Peak and motif analysis revealed that the 5′-TCAG-3′ motif was highly enriched (P = 1e−11) (Fig. 4d), and the majority of six2 peaks were located at the promoters/transcription start sites (Fig. 4e). Importantly, six2 occupied the promoters of CYP4Z1 and CYP4Z2P (Fig. 4f). We subsequently performed a ChIP assay for six2 in MCF-7 cells with or without six2 overexpression. Our results indicated that six2 was indeed bound to the promoters of CYP4Z1 and CYP4Z2P and that six2 overexpression increased six2 binding to these promoters (Fig. 4g, h). Additionally, we constructed luciferase reporter vectors containing different regions of the promoters of CYP4Z1 and CYP4Z2P and found that the luciferase activity of the luciferase reporter vectors containing the putative binding sites was enhanced in MCF-7 cells with six2 overexpression (Additional file 9: Figure S3A–C), whereas the activity of luciferase vectors with mutated binding sites or without putative binding sites was unaffected when the binding sites were mutated (Additional file 9: Figure S3B–D). These results indicated that six2 could directly bind to the 5′-TCAG-3′ motif in the promoters of CYP4Z1 and CYP4Z2P. Consistently, the expression levels of CYP4Z2P and CYP4Z1 were significantly increased or decreased by six2 overexpression or knockdown, respectively (Additional file 9: Figure S3E–H). Notably, a qRT-PCR assay of clinical samples showed that both six2 and CYP4Z1 expression and six2 and CYP4Z2P expression were positively correlated in breast cancer tissues (P < 0.001, Fig. 4i, j). These results suggest that six2 directly binds to the promoters of CYP4Z1 and CYP4Z2P, thus increasing their expression and activating ceRNET_CC.
Transcriptional factor six2 promotes the stemness of breast cancer cells
Notably, we indeed observed a large overlap and positive correlation between gene expression profiles regulated by CYP4Z1-3′UTR- or CYP4Z2P-3′UTR and six2 in MCF-7 cells (Fig. 5a, b). We next explored whether six2 exerts similar effects on ceRNET_CC. First, six2 expression was examined in breast tumor and normal adjacent tissues via online clinical deposited data (http://www.firebrowse.org/). The results showed that six2 expression was significantly increased in breast tumor tissues (Fig. 5c), and a further tissue microarray analysis exhibited consistent results (Fig. 5d). Additionally, an analysis of different clinical samples (http://hgserver1.amc.nl/cgi-bin/r2/main.cgi) showed that six2 expression was negatively correlated with OS, disease-free survival, and RFS of breast cancer patients (Additional file 10: Figure S4A–F). Notably, we identified a set of gene signatures regulating the pluripotency of stem cells in MCF-7 cells with overexpressing CYP4Z2P-3′UTR, CYP4Z1-3′UTR, and six2, and these genes included AMOTL2, CCNA2, ID1/2, GSK3β, and FGF2 (Fig. 5e). A GSEA revealed that six2 overexpression resulted in the enrichment of gene sets related to embryonic stem cell and mammary stem cell function (Fig. 5f, g), consistent with the CYP4Z1-3′UTR or CYP4Z2P-3′UTR overexpression-mediated changes. As expected, six2 overexpression potentiated sphere formation in MCF-7 cells, whereas the knockdown of six2 exerted the opposite effects (Additional file 11: Figure S5A). Moreover, six2 overexpression increased the CD44+/CD24− population, whereas the knockdown of six2 decreased this population (Additional file 11: Figure S5B and S5C). Furthermore, the expression of several pluripotent transcription factors was increased in cells with six2 overexpression and decreased in cells with six2 knockdown (Additional file 11: Figure S5D–G). We further investigated whether six2 could promote the tumor-initiating potential of breast cancer cells in vivo. We compared the capacity of six2-overexpressing or six2-knockdown MCF-7 cells to seed tumors at limiting dilutions. Although all cell lines could form tumors at the densities of 1 × 106 and 1 × 105 cells, the six2-overexpressing cells showed an increase of tumor size and weight, whereas six2-knockdown cells exhibited decreased tumor size and weight (Fig. 6a–d). Furthermore, an in vivo tumorigenic assay was performed with MDA-MB-231 cells after six2 overexpression or knockdown and demonstrated that six2 remarkably attenuated the tumor-initiating potential of MDA-MB-231 cells (Fig. 6e–i). Together, our gain- and loss-of-function assays demonstrate that six2 could promote the stemness of breast cancer cells.
ceRNET_CC is sufficient and necessary for six2-induced effects
We continued investigating whether the ability of six2 to promote the stemness of breast cancer cells is dependent on ceRNET_CC. MCF-7-six2 cells were transfected with siRNA against CYP4Z1- (si-Z1) or CYP4Z2P-3′UTR (si-Z2P), and the resulting cells were subjected to cell spheroid formation and CD44+/CD24− population assays. As shown in Fig. 7a–d, the increased cell spheroid formation, CD44+/CD24− population, and stemness marker expression induced by six2 overexpression were nullified by si-Z1 and/or si-Z2P transfection. CYP4Z1- or CYP4Z2P-3′UTR knockdown decreased the hTERT, p-Akt, and p-ERK1/2 levels or even reversed the six2-mediated increase in the hTERT, p-Akt, and p-ERK1/2 levels (Fig. 7e). Notably, co-transfection with si-Z1 and si-Z2P exerted additive effects. Importantly, we constructed MCF-7-six2 cells with CYP4Z1 or CYP4Z2P knockdown, designated as MCF-7-six2-si-Z1 and MCF-7-six2-si-Z2P cells, respectively, and MCF-7-Plko-six2 cells with CYP4Z1- or CYP4Z2P-3′UTR overexpression, denoted as MCF-7-Plko-six2-Z1-UTR and MCF-7-Plko-six2-Z2P-UTR cells, respectively. An in vivo tumorigenic assay showed that six2 overexpression increased the tumor-initiation capacity of MCF-7 cells, an effect that was fully inhibited by CYP4Z1 or CYP4Z2P knockdown, and the decreased tumor-initiation capacity of MCF-7 cells with six2 knockdown was rescued by CYP4Z1- or CYP4Z2P-3′UTR overexpression (Fig. 7f).
Additionally, we detected whether the six2-mediated promotion of breast cancer cell stemness occurs through the downstream effectors of the ceRNA network between CYP4Z1 and the pseudogene CYP4Z2P and the hTERT/PI3K/Akt and ERK1/2 pathways. To address this issue, we characterized pathways regulated by six2 based on RNA-sequencing data. As expected, we found that PI3K-Akt signaling pathway, phosphatidylinositol signaling system, and signaling pathways regulating pluripotency of cells and breast cancer were among the pathways most highly upregulated by six2 overexpression (Additional file 12: Figure S6A). As expected, the increased cell spheroid formation or CD44+/CD24− population induced by six2 overexpression was attenuated by sihTERT transfection or by LY-294002 or VX-11e treatment (Additional file 12: Figure S6B–D). Moreover, the increased expression of stemness markers induced by six2 overexpression was attenuated by sihTERT (Additional file 12: Figure S6E), or by LY-294002 or VX-11e treatment (Additional file 12: Figure S6F). Furthermore, the overexpression or knockdown of six2 increased or decreased the expression of hTERT, p-ERK1/2, and p-AKT in MDA-MB-231 cells, respectively (Additional file 12: Figure S6G). Consistent with these results, the knockdown of CYP4Z1, CYP4Z2P, or hTERT or the treatment with LY-294002 or VX-11e reduced the ability of six2 overexpression to promote hTERT expression, PI3K/Akt signaling, ERK1/2 signaling, or the expression of cell stemness markers in MDA-MB-231 cells (Additional file 12: Figure S6H). Notably, we found that the overexpression of CYP4Z1- or CYP4Z2P-3′UTR rescued the six2 knockdown-mediated inhibition of breast cancer cell stemness, characterized by increasing cell spheroid formation (Additional file 13: Figure S7A and S7B), CD44+/CD24− population (Additional file 13: Figure S7C), and expression of stemness markers and reactivation of the PI3K/Akt and ERK1/2 signaling (Additional file 13: Figure S7D–G). Thus, these results demonstrate that six2 promotes breast cancer cell stemness dependent on ceRNET_CC.
Six2-mediated regulation of ceRNET_CC renders breast cancer cells resistant to adriamycin treatment by promoting cell stemness
We have established the ability of the six2/ceRNET_CC axis to promote CSC traits in breast cancer cells. As conferring CSC traits have been confirmed to endow chemoresistance to tumor cells , we speculated whether this regulatory axis could decrease adriamycin sensitivity in breast cancer cells by promoting cell stemness. First, the expression of six2, CYP4Z1, and the pseudogene CYP4Z2P was examined in MCF-7 and adriamycin-resistant MCF-7 (MCF-7-Adr) cells via qRT-PCR assay, and the expression was shown to be significantly increased in MCF-7-Adr cells (Fig. 8a). Furthermore, IC50 values were determined in MCF-7 and MDA-MB-231 cells with CYP4Z1- or CYP4Z2P-3′UTR overexpression or knockdown, or with six2 overexpression or knockdown, and we found that overexpression of CYP4Z1- or CYP4Z2P-3′UTR or six2 increased the IC50 values of adriamycin, while knockdown decreased the IC50 values in both MCF-7 and MDA-MB-231 cells (Table 1). MCF-7-Adr cells displayed higher expression of P-gp (a multidrug resistance protein) and pluripotent transcription factors and hyperactivation of PI3K/Akt and ERK1/2 signaling, both of which were attenuated by the knockdown of CYP4Z1- or CYP4Z2P-3′UTR or six2, and the overexpression of CYP4Z1- or CYP4Z2P-3′UTR rescued the attenuation mediated by six2 knockdown (Fig. 8b). Additionally, the knockdown of CYP4Z1- or CYP4Z2P-3′UTR or six2 decreased the spheroid formation ability (Fig. 8c) and the CD44+/CD24− population in MCF-7-Adr cells (Fig. 8d), and the decrease induced by six2 knockdown was rescued by the overexpression of CYP4Z1- or CYP4Z2P-3′UTR. Notably, the knockdown of CYP4Z1- or CYP4Z2P-3′UTR or six2 enhanced adriamycin sensitivity, evidenced by the decrease of IC50 values of adriamycin (Table 1). Consistently, the knockdown of CYP4Z1- or CYP4Z2P-3′UTR or six2 impaired the tumor-initiating potential of MCF-7-Adr cells (Fig. 8e, f). Thus, our results indicate that the six2/ceRNET_CC regulatory axis attenuates adriamycin sensitivity by promoting the stemness of breast cancer cells (Additional file 14: Figure S8).
The pseudogene CYP4Z2P was first identified in 2004 by Rieger et al., who showed that the expression of both the pseudogene CYP4Z2P and its functional gene CYP4Z1 were specifically increased in breast cancer . Further study indicated that the expression of CYP4Z1 and the pseudogene CYP4Z2P was associated with PIK3CA mutations in ERalpha-positive breast cancer . These effects suggest that CYP4Z1 and the pseudogene CYP4Z2P might be involved in breast cancer progression. Our previous studies established that ceRNET_CC promotes angiogenesis , confers tamoxifen resistance , and serves as an anti-apoptotic factor in breast cancer cells . In the present study, we further determined its role and underlying molecular mechanism in regulating the stemness of breast cancer, and we obtained the following findings: (1) the expression levels of CYP4Z1 and the pseudogene CYP4Z2P were higher in breast cancer sphere cells and adriamycin-resistant MCF-7-Adr cells; (2) RNA-sequencing and pathway analysis combined with in vitro and in vivo experiments indicated that ceRNET_CC promoted the stemness of breast cancer cells by activating the PI3K/Akt and ERK1/2 signaling pathways, both of which are critically involved in pro-survival and pro-stem-cell maintenance [40,41,42]; (3) ChIP-sequencing and a ChIP assay confirmed that transcriptional factor six2 directly bound to the 5′-TCAG-3′ motif in the promoters of CYP4Z1 and CYP4Z2P; (4) six2 expression was increased in breast cancer tissues and facilitated the progression of ceRNET_CC, thereby enhancing the stemness of breast cancer cells; (5) the six2-mediated regulation of ceRNET_CC contributed to adriamycin resistance by regulating the stemness of breast cancer cells, and this effect was due at least in part to the activation of PI3K/Akt and ERK1/2 signaling pathways; and (6) the expression of both six2 and CYP4Z1 was negatively correlated with the OS of breast cancer patients and positively correlated with the expression of the stemness marker Nanog, and furthermore, six2 and CYP4Z1 expression was positively correlated in breast cancer tissues. These findings provide the first insights into the roles and molecular mechanism of the six2-mediated regulation of ceRNET_CC, which promotes stemness and chemoresistance in breast cancer.
A previous study indicated that tumor cells with higher expression of hTERT displayed stronger stemness characteristics . This work and our previous work have suggested that ceRNET_CC acts as a subceRNA network for hTERT , which fully supports the role of ceRNET_CC in promoting the stemness of breast cancer cells. Additionally, combined with the current work showing that ceRNET_CC attenuates adriamycin sensitivity in estrogen receptor (ER)-positive and ER-negative breast cancer, our previous study showed that ceRNET_CC confers tamoxifen resistance in ER-positive breast cancer . Additionally, KEGG pathway enrichment showed that endocrine resistance and EGFR tyrosine kinase inhibitor resistance were enriched in MCF-7 cells with CYP4Z1- or CYP4Z2P-3′UTR or six2 overexpression (data not shown). These results strongly support ceRNET_CC as a CSC-related marker because CSCs contribute to drug resistance. Nevertheless, we investigated adriamycin sensitivity in this work because it is one of the first-line drugs used for chemotherapy in breast cancer patients, and it remains unclear whether six2/ceRNET_CC contributes to resistance to other drugs. Notably, due to the relatively lower stemness of MCF-7 cells and higher stemness of MDA-MB-231 and MCF-7-Adr cells [7, 44], we mostly chose to overexpress CYP4Z1-3′UTR or CYP4Z2P-3′UTR or six2 in MCF-7 cells and to knock these factors down in MDA-MB-231 and MCF-7-Adr cells.
Importantly, current treatments for triple-negative breast cancer (TNBC), which is the deadliest form of breast cancer, rely mainly on chemotherapy because no targeted therapies are currently approved for TNBC , raising the potential for ceRNET_CC-targeted therapies in TNBC. Furthermore, a previous study demonstrated that six2 could bind to the E-cadherin promoter and enhance its methylation levels in breast cancer cells . However, the exact six2 binding sites remained unclear. Here, a ChIP-sequencing assay indicated that six2 specifically bound to the 5′-TCAG-3′ motif. This is the first study elucidating the binding sites of six2 in breast cancer. However, it remains unclear whether this is a common phenomenon. GSEA analysis indicated that ACCAAAG_MIR9, CATGTAA_MIR496, GTTATAT_MIR410, ACATATC_MIR190, IKEDA_MIR30_TARGETS_UP, and 3′UTR-mediated translational regulation were enriched in MCF-7-six2 cells (data not shown). These five miRNAs were predicted to bind to the 3′UTR of both CYP4Z1 and the pseudogene CYP4Z2P, and KEGG enrichment showed that miRNAs in cancers were also enriched in MCF-7-six2 cells (data not shown), which suggested that six2 might regulate ceRNET_CC through the 3'UTR and shared miRNAs of CYP4Z1 and the pseudogene CYP4Z2P. These findings should be explored in future work. A previous study indicated that six2 defines and regulates a multipotent self-renewing nephron progenitor population throughout mammalian kidney development , and recent work has shown that normal cells and CSCs might share regulatory mechanisms for maintaining self-renewing capacity and resisting differentiation . Moreover, the YAP1/six2 axis is required for DDX3-mediated tumor aggressiveness and cetuximab resistance in KRAS wild-type colorectal cancer . Notably, the WNT signaling pathway, Jak-STAT signaling pathway, and FoxO signaling pathway, which are involved in regulating cell stemness [48,49,50], were enriched in MCF-7 cells with CYP4Z1- or CYP4Z2P-3′UTR or six2 overexpression (data not shown). Although we cannot exclude the possibility that six2 may still function through additional signaling pathways to regulate the stemness of breast cancer cells, our study firmly establishes the critical roles of six2-mediated regulation of ceRNET_CC in these processes.
Notably, six2 protein levels were gradually increased by adriamycin treatment in a concentration-dependent manner (data not shown), suggesting that six2 overexpression and adriamycin resistance might form a positive feedback loop in breast cancer cells. Despite efforts to develop chemotherapies for killing CSCs over the past decades and evidence of early success [43, 51, 52], there have been significant setbacks, presumably due to limited effectiveness in late-stage clinical trials. In addition to toxicity and side effects, the reasons for the setbacks could also be the lack of predictive biomarkers for patients. Our finding that the six2/ceRNET_CC regulatory axis promotes the stemness of breast cancer provides putative targets for the development of new strategies for targeting and compromising the maintenance of breast cancer stemness, and we hypothesize that a gene expression signature comprising all three of these genes will predict chemotherapeutic sensitivity in breast cancer patients.
Competitive endogenous RNA
ceRNA network between CYP4Z1 and pseudogene CYP4Z2P
Cancer stem cells
Fetal bovine serum
Gene Expression Omnibus
Gene set enrichment analysis
Non-small-cell lung cancer
Quantitative real-time PCR
Small interfering RNA
Short tandem repeat
Wang L, Leite de Oliveira R, Huijberts S, Bosdriesz E, Pencheva N, Brunen D, Bosma A, Song JY, Zevenhoven J, Los-de Vries GT, et al. An acquired vulnerability of drug-resistant melanoma with therapeutic potential. Cell. 2018;173(6):1413–25 e1414.
Mehta M, Basalingappa K, Griffith JN, Andrade D, Babu A, Amreddy N, Muralidharan R, Gorospe M, Herman T, Ding WQ, et al. HuR silencing elicits oxidative stress and DNA damage and sensitizes human triple-negative breast cancer cells to radiotherapy. Oncotarget. 2016;7(40):64820–35.
Chen DL, Lu YX, Zhang JX, Wei XL, Wang F, Zeng ZL, Pan ZZ, Yuan YF, Wang FH, Pelicano H, et al. Long non-coding RNA UICLM promotes colorectal cancer liver metastasis by acting as a ceRNA for microRNA-215 to regulate ZEB2 expression. Theranostics. 2017;7(19):4836–49.
Gao S, Wang P, Hua Y, Xi H, Meng Z, Liu T, Chen Z, Liu L. ROR functions as a ceRNA to regulate Nanog expression by sponging miR-145 and predicts poor prognosis in pancreatic cancer. Oncotarget. 2016;7(2):1608–18.
Wilkinson F, Pratt H, Atchison ML. PcG recruitment by the YY1 REPO domain can be mediated by Yaf2. J Cell Biochem. 2010;109(3):478–86.
Srinivasan L, Pan X, Atchison ML. Transient requirements of YY1 expression for PcG transcriptional repression and phenotypic rescue. J Cell Biochem. 2005;96(4):689–99.
Prat A, Parker JS, Karginova O, Fan C, Livasy C, Herschkowitz JI, He X, Perou CM. Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res. 2010;12(5):R68.
Zheng L, Meng X, Li X, Zhang Y, Li C, Xiang C, Xing Y, Xia Y, Xi T. miR-125a-3p inhibits ERalpha transactivation and overrides tamoxifen resistance by targeting CDK3 in estrogen receptor-positive breast cancer. FASEB J. 2018;32(2):588–600.
Tao F, Tian X, Ruan S, Shen M, Zhang Z. miR-211 sponges lncRNA MALAT1 to suppress tumor growth and progression through inhibiting PHF19 in ovarian carcinoma. FASEB J. 2018:fj201800495RR. https://doi.org/10.1096/fj.201800495RR.
Ooi CY, Carter DR, Liu B, Mayoh C, Beckers A, Lalwani A, Nagy Z, De Brouwer S, Decaesteker B, Hung TT, et al. Network modeling of microRNA-mRNA interactions in neuroblastoma tumorigenesis identifies miR-204 as a direct inhibitor of MYCN. Cancer Res. 2018;78(12):3122–34.
Yang Y, Li F, Saha MN, Abdi J, Qiu L, Chang H. miR-137 and miR-197 induce apoptosis and suppress tumorigenicity by targeting MCL-1 in multiple myeloma. Clin Cancer Res. 2015;21(10):2399–411.
Jin C, Rajabi H, Kufe D. miR-1226 targets expression of the mucin 1 oncoprotein and induces cell death. Int J Oncol. 2010;37(1):61–9.
Li C, Zheng L, Xin Y, Tan Z, Zhang Y, Meng X, Wang Z, Xi T. The competing endogenous RNA network of CYP4Z1 and pseudogene CYP4Z2P exerts an anti-apoptotic function in breast cancer. FEBS Lett. 2017;591(7):991–1000.
Jeselsohn R, Cornwell M, Pun M, Buchwalter G, Nguyen M, Bango C, Huang Y, Kuang Y, Paweletz C, Fu X, et al. Embryonic transcription factor SOX9 drives breast cancer endocrine resistance. Proc Natl Acad Sci U S A. 2017;114(22):E4482–91.
Li J, Xi Y, Li W, McCarthy RL, Stratton SA, Zou W, Li W, Dent SY, Jain AK, Barton MC. TRIM28 interacts with EZH2 and SWI/SNF to activate genes that promote mammosphere formation. Oncogene. 2017;36(21):2991–3001.
Self M, Lagutin OV, Bowling B, Hendrix J, Cai Y, Dressler GR, Oliver G. Six2 is required for suppression of nephrogenesis and progenitor renewal in the developing kidney. EMBO J. 2006;25(21):5214–28.
Kobayashi A, Valerius MT, Mugford JW, Carroll TJ, Self M, Oliver G, McMahon AP. Six2 defines and regulates a multipotent self-renewing nephron progenitor population throughout mammalian kidney development. Cell Stem Cell. 2008;3(2):169–81.
Wang CA, Drasin D, Pham C, Jedlicka P, Zaberezhnyy V, Guney M, Li H, Nemenoff R, Costello JC, Tan AC, et al. Homeoprotein Six2 promotes breast cancer metastasis via transcriptional and epigenetic control of E-cadherin expression. Cancer Res. 2014;74(24):7357–70.
Combes AN, Wilson S, Phipson B, Binnie BB, Ju A, Lawlor KT, Cebrian C, Walton SL, Smyth IM, Moritz KM, et al. Haploinsufficiency for the Six2 gene increases nephron progenitor proliferation promoting branching and nephron number. Kidney Int. 2017;93:589–98.
Chakrabarti R, Wei Y, Hwang J, Hang X, Andres Blanco M, Choudhury A, Tiede B, Romano RA, DeCoste C, Mercatali L, et al. DeltaNp63 promotes stem cell activity in mammary gland development and basal-like breast cancer by enhancing Fzd7 expression and Wnt signalling. Nat Cell Biol. 2014;16(10):1004–15 1001-1013.
Guo W, Keckesova Z, Donaher JL, Shibue T, Tischler V, Reinhardt F, Itzkovitz S, Noske A, Zurrer-Hardi U, Bell G, et al. Slug and Sox9 cooperatively determine the mammary stem cell state. Cell. 2012;148(5):1015–28.
Shimono Y, Zabala M, Cho RW, Lobo N, Dalerba P, Qian D, Diehn M, Liu H, Panula SP, Chiao E, et al. Downregulation of miRNA-200c links breast cancer stem cells with normal stem cells. Cell. 2009;138(3):592–603.
Xu HX, Wu KJ, Tian YJ, Liu Q, Han N, He XL, Yuan X, Wu GS, Wu KM. Expression profile of SIX family members correlates with clinic-pathological features and prognosis of breast cancer: a systematic review and meta-analysis. Medicine. 2016;95(27):e4085.
Song D, Yue L, Wu G, Ma S, Guo L, Yang H, Liu Q, Zhang D, Xia Z, Wang L, et al. Assessment of promoter methylation and expression of SIX2 as a diagnostic and prognostic biomarker in Wilms’ tumor. Tumour Biol. 2015;36(10):7591–8.
Liu Q, Li A, Tian Y, Liu Y, Li T, Zhang C, Wu JD, Han X, Wu K. The expression profile and clinic significance of the SIX family in non-small cell lung cancer. J Hematol Oncol. 2016;9(1):119.
Zheng L, Li X, Gu Y, Lv X, Xi T. The 3’UTR of the pseudogene CYP4Z2P promotes tumor angiogenesis in breast cancer by acting as a ceRNA for CYP4Z1. Breast Cancer Res Treat. 2015;150(1):105–18.
Lanczky A, Nagy A, Bottai G, Munkacsy G, Szabo A, Santarpia L, Gyorffy B. miRpower: a web-tool to validate survival-associated miRNAs utilizing expression data from 2178 breast cancer patients. Breast Cancer Res Treat. 2016;160(3):439–46.
Zheng L, Xiang C, Li X, Guo Q, Gao L, Ni H, Xia Y, Xi T. STARD13-correlated ceRNA network-directed inhibition on YAP/TAZ activity suppresses stemness of breast cancer via co-regulating Hippo and Rho-GTPase/F-actin signaling. J Hematol Oncol. 2018;11(1):72.
Hu Y, Smyth GK. ELDA: extreme limiting dilution analysis for comparing depleted and enriched populations in stem cell and other assays. J Immunol Methods. 2009;347(1–2):70–8.
Zheng L, Li X, Gu Y, Ma Y, Xi T. Pseudogene CYP4Z2P 3’UTR promotes angiogenesis in breast cancer. Biochem Biophys Res Commun. 2014;453(3):545–51.
Wong DJ, Liu H, Ridky TW, Cassarino D, Segal E, Chang HY. Module map of stem cell genes guides creation of epithelial cancer stem cells. Cell Stem Cell. 2008;2(4):333–44.
Dontu G, Abdallah WM, Foley JM, Jackson KW, Clarke MF, Kawamura MJ, Wicha MS. In vitro propagation and transcriptional profiling of human mammary stem/progenitor cells. Genes Dev. 2003;17(10):1253–70.
Liu R, Shi P, Nie Z, Liang H, Zhou Z, Chen W, Chen H, Dong C, Yang R, Liu S, et al. Mifepristone suppresses basal triple-negative breast cancer stem cells by down-regulating KLF5 expression. Theranostics. 2016;6(4):533–44.
Cordenonsi M, Zanconato F, Azzolin L, Forcato M, Rosato A, Frasson C, Inui M, Montagner M, Parenti AR, Poletti A, et al. The Hippo transducer TAZ confers cancer stem cell-related traits on breast cancer cells. Cell. 2011;147(4):759–72.
Xu J, Wong EY, Cheng C, Li J, Sharkar MT, Xu CY, Chen B, Sun J, Jing D, Xu PX. Eya1 interacts with Six2 and Myc to regulate expansion of the nephron progenitor pool during nephrogenesis. Dev Cell. 2014;31(4):434–47.
Yang M, Liu P, Huang P. Cancer stem cells, metabolism, and therapeutic significance. Tumour Biol. 2016;37:5735–42.
Rieger MA, Ebner R, Bell DR, Kiessling A, Rohayem J, Schmitz M, Temme A, Rieber EP, Weigle B. Identification of a novel mammary-restricted cytochrome P450, CYP4Z1, with overexpression in breast carcinoma. Cancer Res. 2004;64(7):2357–64.
Cizkova M, Cizeron-Clairac G, Vacher S, Susini A, Andrieu C, Lidereau R, Bieche I. Gene expression profiling reveals new aspects of PIK3CA mutation in ERalpha-positive breast cancer: major implication of the Wnt signaling pathway. PLoS One. 2010;5(12):e15647.
Zheng L, Li X, Meng X, Chou J, Hu J, Zhang F, Zhang Z, Xing Y, Liu Y, Xi T. Competing endogenous RNA networks of CYP4Z1 and pseudogene CYP4Z2P confer tamoxifen resistance in breast cancer. Mol Cell Endocrinol. 2016;427:133–42.
Majumder M, Xin X, Liu L, Tutunea-Fatan E, Rodriguez-Torres M, Vincent K, Postovit LM, Hess D, Lala PK. COX-2 induces breast cancer stem cells via EP4/PI3K/AKT/NOTCH/WNT axis. Stem Cells. 2016;34(9):2290–305.
Kwon SJ, Kwon OS, Kim KT, Go YH, Yu SI, Lee BH, Miyoshi H, Oh E, Cho SJ, Cha HJ. Role of MEK partner-1 in cancer stemness through MEK/ERK pathway in cancerous neural stem cells, expressing EGFRviii. Mol Cancer. 2017;16(1):140.
Huang B, Fu SJ, Fan WZ, Wang ZH, Chen ZB, Guo SJ, Chen JX, Qiu SP. PKCepsilon inhibits isolation and stemness of side population cells via the suppression of ABCB1 transporter and PI3K/Akt, MAPK/ERK signaling in renal cell carcinoma cell line 769P. Cancer Lett. 2016;376(1):148–54.
Bonuccelli G, Peiris-Pages M, Ozsvari B, Martinez-Outschoorn UE, Sotgia F, Lisanti MP. Targeting cancer stem cell propagation with palbociclib, a CDK4/6 inhibitor: telomerase drives tumor cell heterogeneity. Oncotarget. 2017;8(6):9868–84.
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(5):499–507.
Duhachek-Muggy S, Qi Y, Wise R, Alyahya L, Li H, Hodge J, Zolkiewska A. Metalloprotease-disintegrin ADAM12 actively promotes the stem cell-like phenotype in claudin-low breast cancer. Mol Cancer. 2017;16(1):32.
Celia-Terrassa T, Liu DD, Choudhury A, Hang X, Wei Y, Zamalloa J, Alfaro-Aco R, Chakrabarti R, Jiang YZ, Koh BI, et al. Normal and cancerous mammary stem cells evade interferon-induced constraint through the miR-199a-LCOR axis. Nat Cell Biol. 2017;19(6):711–23.
Wu DW, Lin PL, Wang L, Huang CC, Lee H. The YAP1/SIX2 axis is required for DDX3-mediated tumor aggressiveness and cetuximab resistance in KRAS-wild-type colorectal cancer. Theranostics. 2017;7(5):1114–32.
Nabhan AN, Brownfield DG, Harbury PB, Krasnow MA, Desai TJ. Single-cell Wnt signaling niches maintain stemness of alveolar type 2 cells. Science. 2018;359(6380):1118–23.
Ojha R, Singh SK, Bhattacharyya S. JAK-mediated autophagy regulates stemness and cell survival in cisplatin resistant bladder cancer cells. Biochim Biophys Acta. 2016;1860(11 Pt A):2484–97.
Gargini R, Cerliani JP, Escoll M, Anton IM, Wandosell F. Cancer stem cell-like phenotype and survival are coordinately regulated by Akt/FoxO/Bim pathway. Stem Cells. 2015;33(3):646–60.
Doherty MR, Cheon H, Junk DJ, Vinayak S, Varadan V, Telli ML, Ford JM, Stark GR, Jackson MW. Interferon-beta represses cancer stem cell properties in triple-negative breast cancer. Proc Natl Acad Sci U S A. 2017;114:13792–7.
Lee KM, Giltnane JM, Balko JM, Schwarz LJ, Guerrero-Zotano AL, Hutchinson KE, Nixon MJ, Estrada MV, Sanchez V, Sanders ME, et al. MYC and MCL1 cooperatively promote chemotherapy-resistant breast cancer stem cells via regulation of mitochondrial oxidative phosphorylation. Cell Metab. 2017;26(4):633–47 e637.
This work was critically reviewed by Prof. Changying Guo from our university, and we thanked him for his valuable advice.
This work was supported by the National Nature Science Foundation of China, No. 81702957, the China Postdoctoral Science Foundation, No. 2017 M620230, the Postdoctoral Research Funding Scheme of Jiangsu Province (2017), No. 1701197B, and the Priority Academic Program Development (PAPD) of Jiangsu Higher Education Institutions.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Ethics approval and consent to participate
Written informed consent from all patients and approval from the Huai An First People’s Hospital and Jiangsu Province Hospital of TCM ethics review committees were obtained.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Sequences of primers used for qRT-PCR in this study. (DOCX 17 kb)
Table S2. Primary antibodies used in this study. (DOC 35 kb)
Table S3. Sequences of primers used for ChIP qRT-PCR in this study. (DOCX 16 kb)
Table S4. Sequences of siRNA against specific target in this study (DOC 32 kb)
Table S5. Sequences of primers used for plasmid constructions. (DOC 40 kb)
Table S6. Sequences of primers used for Luciferase reporter assay. (DOC 46 kb)
Figure S1. CeRNET_CC promotes the stemness of MCF-7 cells in vitro. (A) The expression of CYP4Z2P and CYP4Z1 in MCF-7 and MCF-7-tumorsphere cells was detected by qRT-PCR. (B and C) The infection efficiency of MCF-7 cells with CYP4Z1- or CYP4Z2P-3′UTR stable overexpression (B) or knockdown (C) was detected by qRT-PCR. (D) Phase contrast images of mammospheres formed by stable expression cells depicted in B and C and quantification of spheres. (E) Representative FACS profile of cells described in B with CD24− and CD44+ markers. (F and G) The mRNA and protein expression of stemness markers (ALDH1, SOX2, OCT4 and Nanog) in cells described in B and C were examined by qRT-PCR and western blot analysis, respectively. The data are presented as the means ± SDs, n = 3, *P < 0.05, **P < 0.01, ***P < 0.001 vs. MCF-7. (PDF 5600 kb)
Figure S2. CeRNET_CC promotes the stemness of MDA-MB-231 cells in vitro. (A and B) The infection efficiency of MDA-MB-231 cells with CYP4Z1- (A) or CYP4Z2P-3′UTR (B) stable knockdown was detected by qRT-PCR. (C) Phase contrast images of mammospheres formed by stable expression cells depicted in A and B and quantification of spheres. (D) Representative FACS profile of cells described in A and B with CD24− and CD44+ markers. (E and F) The mRNA and protein expression of stemness markers (ALDH1, SOX2, OCT4, and Nanog) in cells described in A and B. (G) Pearson correlation analysis of the expression of CYP4Z1 and Nanog in basal-like breast cancer (n = 252, P < 0.01). (H) Pearson correlation analysis of the expression of the pseudogene CYP4Z2P and Nanog in basal-like breast cancer (n = 144, P < 0.01). (PDF 3450 kb)
Figure S3. Transcriptional factor six2 promotes the expression of CYP4Z1 and the pseudogene CYP4Z2P. (A) Computational analysis of the CYP4Z2P and CYP4Z1 promoters showed potential binding sites for six2. (B and C) Fragments of the CYP4Z2P and CYP4Z1 promoters were cloned into the luciferase reporter vector pGL3. MCF-7 cells were co-transfected with six2 and the luciferase constructs or the control construct. 72 h later, luciferase activity was measured. (D) Relative luciferase activity was detected in MCF-7 cells co-transfected with the six2 overexpression vector and CYP4Z1 and CYP4Z2P promoter vectors with six2 binding sites or mutated six2 binding sites. (E and F) The expression of CYP4Z2P, CYP4Z1, and six2 in MCF-7-six2 (E) and 231-six2 (F) cells was examined by qRT-PCR. (G and H) CYP4Z1 protein expression in MCF-7-six2 (G) and MCF7-Plko-six2 (H) cells was detected by western blot. The data were presented as the means ± SDs, n = 3, **P < 0.01 vs. control or MCF-7. (PDF 1270 kb)
Figure S4. The correlation between six2 expression and the survival of breast cancer patients. (A) KM-plotter survival curves showed the disease-free survival probability of patients separated into low and high six2 levels. (B) KM-plotter survival curves showed the OS survival probability of patients separated into low and high six2 level. (C and F) KM-plotter survival curves showed the RFS probability of patients separated into low and high six2 levels (http://hgserver1.amc.nl/cgi-bin/r2/main.cgi). (PDF 576 kb)
Figure S5. Transcriptional factor six2 promotes the stemness of breast cancer cells in vitro. (A) Phase contrast images of mammospheres formed by MCF-7 cells with or without six2 overexpression or knockdown and MDA-MB-231 cells with or without six2 knockdown; spheres were quantified. (B and C) Representative FACS profile of MCF-7 cells with or without six2 overexpression (B) and MDA-MB-231 cells (C) with or without six2 knockdown, with the CD24− and CD44+ markers. (D and E) The mRNA expression of stemness markers (ALDH1, SOX2, OCT4, and Nanog) in MCF-7 (D) or MDA-MB-231(E) cells with six2 stable overexpression or knockdown was detected by qRT-PCR. (F and G) Cells depicted in D and E were subjected to western blot analysis and followed by detecting the expression of six2 and stemness markers (ALDH1, SOX2, and OCT3/4). The data are presented as the mean ± SD, n = 3, *P < 0.05, **P < 0.01, ***P < 0.001 vs. MCF-7 or MDA-MB-231. (PDF 4250 kb)
Figure S6. Transcriptional factor six2 promotes the stemness of breast cancer cells partly through the hTERT/PI3K/Akt and ERK1/2 pathways. (A) Functional annotation analysis of genes coordinately activated by six2 overexpression. (B and C) Phase contrast images of mammospheres formed by MCF-7-six2 cells with LY-294002, VX-11e, or sihTERT treatment (B), and spheres were quantified (C). The data are presented as the means ± SDs, n = 3, *P < 0.05 vs. MCF-7-six2. (D) Representative FACS profile of cells described in B with CD24− and CD44+ markers. (E and F) Cells depicted in B were subjected to western blot analysis and followed by detecting the expression of p-Akt/p-ERK1/2 and stemness markers (ALDH1 and OCT3/4). (G) MDA-MB-231 cells with six2 stable overexpression or knockdown were subjected to western blot analysis and followed by detecting the expression of p-Akt/p-ERK1/2 and hTERT. (H) 231-six2 cells with LY-294002, VX-11e, or sihTERT, or si-Z1, or si-Z2P treatment were subjected to western blot analysis and followed by detecting the expression of p-Akt/p-ERK1/2 and stemness markers (ALDH1 and OCT3/4). (PDF 5040 kb)
Figure S7. CeRNET_CC is sufficient and necessary for six2-induced effects. (A and B) Phase contrast images of mammospheres formed by MCF-7-Plko-six2 cells with Z1-UTR or Z2P-UTR overexpression (A); spheres were quantified (B). The data are presented as the means ± SDs, n = 3, *P < 0.05, **P < 0.01 vs. MCF-7 or MCF-7-six2. (C) Representative FACS profile of cells described in A with CD24− and CD44+ markers. (D–G) MCF-7-Plko-six2 (D and F) and 231-Plko-six2 (E and G) cells with Z1-UTR, Z2P-UTR (D and E), or hTERT (F and G) overexpression were subjected to western blot analysis and followed by detecting the expression of p-Akt/p-ERK1/2 and stemness markers (ALDH1 and OCT3/4). The data were presented as the mean ± SD, n = 3, *P < 0.05, **P < 0.01 vs. MCF-7 or MCF-7-six2. (PDF 3760 kb)
Figure S8. Proposed model in which transcriptional factor six2-mediated regulation of ceRNET_CC is responsible for breast CSC formation and thus drug resistance. Transcriptional factor six2 induces the progression of ceRNET_CC by directly binding to the promoters of CYP4Z1 and the pseudogene CYP4Z2P. This six2/ceRNET_CC regulatory axis results in breast CSC progression and thus enhances drug sensitivity. (TIF 242 kb)