Oncogenic miR-20b-5p contributes to malignant behaviors of breast cancer stem cells by bidirectionally regulating CCND1 and E2F1

Background Breast cancer is the leading cause of cancer mortality in women worldwide. Therefore, it is of great significance to identify the biological mechanism of tumorigenesis and explore the development of breast cancer to achieve a better prognosis for individuals suffering from breast cancer. MicroRNAs (miRNAs) have become a hot topic in cancer research, but the underlying mechanism of its involvement in cancer remains unclear. Methods The miRNA profile between breast cancer stem cells (BCSCs, CD44+CD24−/low) and control MCF-7 breast cancer cells was obtained in a previous study. Based on biological analysis, miR-20b-5p was hypothesized to be a key factor due to the malignant behavior of BCSCs. Then, agomir-20b-5p and antagomir-20b-5p were transfected into MCF-7 and T47D breast cancer cells to detect cell migration, wound healing and proliferation, and lentivirus vectors silencing or overexpressing miR-20b-5p were transfected into T47D-CSCs to detect proliferation and apoptosis. The effect of miR-20b-5p on xenograft growth was investigated in vivo by transfection of a lentivirus-overexpression vector into T47D cells. The target genes were predicted by the online programs picTar, miRanda and TargetScan and verified by dual luciferase assay, and changes in protein expression were detected by western blot. Results MiR-20b-5p had the highest degree in both the miRNA-gene network and miRNA-GO network to regulate BCSCs. Overexpression of miR-20b-5p significantly promoted the migration and wound healing ability of MCF-7 cells and T47D cells compared with the control (P < 0.05). In addition, miR-20b-5p facilitated the proliferation of MCF-7 cells and T47D-CSCs (P < 0.05) and inhibited the apoptosis of T47D-CSCs (P < 0.05). Moreover, miR-20b-5p promoted xenograft growth compared with the control group (P < 0.05). Accordingly, potential targets of both CCND1 and E2F1 were predicted by bioinformatics analysis. MiR-20b-5p directly targeted both CCND1 and E2F1 in a dual luciferase assay, while antagomir-20b-5p downregulated the protein levels of CCND1 and E2F1. Conclusions Oncogenic miR-20b-5p was confirmed to promote the malignant behaviors of breast cancer cells and BCSCs. The underlying mechanism lies in that miR-20b-5p overall enhanced both CCND1 and E2F1 targets via bidirectional regulation probably involving direct downregulation and indirect upregulation.


Background
Breast cancer is one of the most common cancers in women worldwide [1]. In recent years, although significant progress has been made in the diagnosis and treatment of malignant tumors, breast cancer remains the leading cause of cancer-related deaths in women, largely due to its high rate of relapse, metastasis, and drug resistance [2][3][4][5]. Therefore, it is of great significance to explore the mechanism of tumorigenesis in breast cancer and further identify key regulatory factors. Recent studies, including our previous findings, have demonstrated that miRNAs play an important role in regulating biological functions [6][7][8]. Because of their unique expression profile and special functions, miRNAs have been used as molecular markers to define the direction and processes of cell differentiation by regulating one or several target genes. MiRNAs can also control the selfrenewal or proliferation of stem cells [9]. Abnormal changes in miRNAs have been discovered in breast cancer, and miRNAs have been demonstrated to be associated with the tumorigenesis and progression of breast cancer [8,10]. Previous studies have confirmed that miRNAs play an important role in gene regulation in CD44 + CD24 −/low breast cancer stem cells (BCSCs). For example, we found that miR-200C was significantly downregulated in BCSCs and played a critical role in their biological features [8]. In addition, let-7a plays an important role in BCSCs' self-renewal by inhibiting the expression of H-Ras [11,12].
In our previous research, we performed miRNA profiling between sorted CD44 + CD24 −/low BCSCs and the control MCF-7 breast cancer cells [8]. Finally, miR-20b-5p was chosen for further study due to its highest degree of regulation. We found that miR-20b-5p belongs to the miR-106a-363 cluster, which together with the miR-17-92 cluster and miR-106b-25 cluster forms a large family of highly similar miRNAs called the miR-17 family [13]. In many human malignancies, members of the miR-17 family have been reported to accumulate in tumor cells and are speculated to exert oncogenic effects [14,15], and miR-20b serves as a potential oncogene in gastric cancer [16], breast cancer [17,18] and hepatocellular carcinoma [19]. However, there has been no studies exploring the regulatory role of miR-20b in BCSCs, and the underlying mechanism of miR-20b-5p remains unclear. This study aimed to investigate the role of the miR-20b-5p, a subtype of miR-20b family, in the regulation of the malignant behavior of BCSCs and to identify its target genes.

Methods
Induction culture of BCSCs MCF-7 and T47D cells were cultured in DMEM-H (HyClone, USA) containing 10% fetal bovine serum (FBS) at 37°C in a 5% CO 2 incubator. Cells in the logarithmic growth phase, approximately 70-80% confluent, were harvested and digested into a single cell suspension by trypsin digestion.
The complete MammoCult™ medium containing 5% MammoCult™ proliferation supplement, 4 μg/mL heparin, and 0.48 μg/mL hydrocortisone was used. MCF-7 cells were resuspended in complete MammoCult™ medium, and T47D cells resuspended in complete DMEM/F-12 medium. The cell suspension (4 × 10 3 ) was seeded into a 6-well plate and cultured at 37°C and 5% CO 2 . After 7 days of culture, the spheres were collected and centrifuged at 350 g for 5 min, and the supernatant was discarded. Then, 1 mL of Accutase cell dispersion solution was added to digest the spheres, followed by adding 9 mL of sterile PBS solution and centrifugation at 350 g for 5 min. Finally, the supernatant was discarded, and the cells were collected. CD44 + CD24 −/low BCSCs were isolated from MCF-7 and T47D cells by staining with CD44-APC, CD24-PE and ESA-FITC (BD Pharmingen, USA) antibodies via FACS as described in our previous research [8,20].

miRNA profile and miRNA network
Both human miRNA microarray fabrication and hybridization were performed as described previously [8]. The miRNA profiles of both BCSCs and control MCF-7 breast cancer cells were obtained from Capital-Bio Corporation (Beijing, China). All microarray data were uploaded and submitted to the public repository Gene Expression Omnibus (GEO) database (http://www. ncbi.nlm.nih.gov/geo/query/acc.cgi? acc = GSE68271). The differential miRNA profiles between BCSCs and MCF-7 cells were obtained [20].
To identify the miRNAs that regulate BCSCs, we established a miRNA-gene network and miRNA-GO network through analysis of the significant target genes and Gene Ontology (GO) terms performed by Shanghai Qiming Corporation (Shanghai, China).

Target prediction by bioinformatics
Chromosome localization, sequence analysis and target prediction of the miRNAs were carried out by the online programs picTar (http://pictar.mdc-berlin.de/), miRanda (http://microrna.sanger.ac.uk), and TargetScan (http:// www.targetscan.org). The mRNAs associated with cell proliferation and the cell cycle predicted by at least three algorithms were selected as putative targets. Then, the binding free energy (△G) of the hybridization between miRNAs and their 3′ UTR complementary sites was analyzed with the aid of Mfold software. The mRNAs able to bind with lower free energy at both 5′-70 bp and 3′-70 bp than their average random free energy were deemed to be accessible to specific miRNAs [20,21].

Transwell migration assay
Transient transfection of agomir-20b-5p and antagomir-20b-5p in MCF-7 and T47D breast cancer cells was performed as follows. After digestion, breast cancer cells were mixed with 1 mL DMEM-H containing 10% FBS and subjected to resuspension followed by cell counting. Then, the 24-well plates were inoculated at 4 × 10 4 cells per well, containing 500 μL of medium, and placed in a 37°C, 5% CO 2 incubator overnight. The medium was removed the next day, and the cells were washed with PBS three times, and then 500 μL of Opti-MEM medium was added to each well. Next, Lipofectamine 2000 (Invitrogen, USA) was mixed with 1 μL of agomir of miR-20b-5p (agomir-20b-5p), antagomir of miR-20b-5p (antagomir-20b-5p) (Dharmacon, USA), or miR-control at a final concentration of 30 nM, and the cells were incubated for 24-48 h to achieve successful transfection.
The transwell chamber was placed in a 24-well plate. A total of 800 μL of DMEM-H containing 15% FBS was added to the lower chamber, and 200 μL transfected MCF-7 or T47D breast cancer cells (4 × 10 4 cells resuspended in 3% DMEM-H) was added to the upper chamber and cultured in a 37°C, 5% CO 2 incubator for 24 h. After removal of the medium, the cells were fixed with 1 mL of 95% ethanol solution for 10 min, washed 3 times with 1 mL of PBS for 5 min each, and then stained with 1% crystal violet dye solution for 10 min. The crystal violet dye in the upper chamber was washed off under a small stream of water and gently wiped with a cotton swab. Finally, the cells were observed under an inverted microscope followed by imaging and counting (Leica, Germany). The test was repeated three times.

Wound healing assay
Transient transfection of agomir-20b-5p and antagomir-20b-5p into MCF-7 and T47D breast cancer cells was performed as described above. First, 3-4 straight lines were drawn vertically on the back of the 6-well plate. Then, the cells were seeded according to the grouping and cultured in DMEM-H containing 10% FBS. When the transfected MCF-7 and T47D breast cancer cells reached 60-70% confluence, mitomycin was added for 2 h (final concentration was 10 ng/mL). After washing 3 times with PBS, the cells were scratched gently with a 100 μL pipette tip according to the marked lines, and then the detached cells were removed by washing with PBS. Under an inverted microscope (Leica, Germany), the imaging position was recorded at the 24 h and 48 h time points. The test was repeated three times.

Proliferation assay in MCF-7 cells and T47D-CSCs
Transient transfection of agomir-20b-5p and antagomir-20b-5p into MCF-7 breast cancer cells was performed as described above. MCF-7 cells were incubated with 30 μM EdU, and cell proliferation was evaluated by flow cytometry using an EdU assay kit (GeneCopoeia, USA) according to the manufacturer's instructions. The test was repeated five times. Stable transfection of a lentivirus overexpression system was used in the T47D-CSC proliferation assays. The lentiviral vector (GV369) overexpressing miR-20b-5p was purchased from GeneChem (China). The miR-20b-5p overexpression vector was labeled T47D-CSCs/LV-miR-20b-5p, and the control vector was labeled T47D-CSCs/LV-NC. The induced T47D-CSCs were suspended in EpiCult-B serum-free medium, and Lipofectamine 2000 (Invitrogen, USA) was added together with T47D-CSCs/LV-miR-20b-5p or T47D-CSCs/LV-NC. Then, the viruses were harvested. Cultured cells were infected with 5 μg/mL polybrene and lentivirus (MOI = 50) for 24 h, and then the cells were incubated with fresh medium for another 48 h to establish stable cell lines. The proliferation assay was conducted as described above. The test was repeated five times.

Apoptosis assay in T47D-CSCs
The stable transfection of lentivirus overexpression and lentivirus-miR-sponge constructs was used to detect apoptosis in T47D-CSCs as described above. The induced T47D-CSCs were stained with anti-annexin V/7-AAD antibodies in binding buffer for 15 min and then analyzed for apoptosis by flow cytometry. The test was repeated five times.

Animal experiments
The stable transfection of T47D cells overexpressing lentivirus was conducted as described above, and the cells were used in animal experiments. The vector overexpressing miR-20b-5p was labeled as T47D/LV-miR-20b-5p, and the control vector was labeled as T47D/LV-NC. All animal experiments were approved by the Institutional Animal Care and Use Committee of Xinqiao Hospital, Army Medical University. A total of 15 nude mice (4-6 weeks old) from the SPF Laboratory Animal Center of Xinqiao Hospital were randomly divided into the following three groups with 5 mice each: T47D control, T47D/LV-NC and T47D/LV-miR-20b-5p. A total of 1 × 10 7 cells/100 μL of PBS were subcutaneously injected to establish a xenograft model. The tumor size was measured every 3 days and calculated according to the formula V = (length × width 2 )/2. After 4 weeks, the mice were sacrificed by cutting off the head after injection of Chloral hydrate, and the tumor tissues were collected for miR-20b-5p detection by qRT-PCR. All protocols involving mice were conducted in accordance with the animal care guidelines of Xinqiao Hospital.

Dual luciferase assay
Two hundred ninety-three T cells were cultured with DMEM (supplemented with 10% FBS) to a cell density of 70-80%. We followed the protocol of our previous work for experimental design [22]. With Renilla luciferase as the internal reference, the RLU value determined for firefly luciferase was divided by the Renilla luciferase value. The vectors (GV306) for CCND1 and E2F1 were purchased from GeneChem (China). The inserted sequences were as follows: wild type CCND1 (CCND1-WT) and mutated CCND1 (CCND1-MT), wild type E2F1 (E2F1-WT) and mutated E2F1 (E2F1-MT) (Supplementary Fig. 1). A random sequence and miR-20b-5p mimics were used for negative control (NC) and miRNA overexpression, respectively. The experiment was divided into 4 groups: NC + CCND1-WT, MIMICS+ CCND1-WT, NC + CCND1-MT, and MIMICS+ CCND1-MT. The same experiment was performed with E2F1 reagents.

Western blot analysis of target genes
Transfected cells were washed twice with PBS, lysed in RIPA buffer (Beyotime, China) containing a protease inhibitor cocktail (Roche, Switzerland) for 30 min and then centrifuged at 13000×g for 10 min. After boiling in loading buffer for 10 min, the supernatants were loaded onto SDS-PAGE gels and then transferred onto PVDF membranes (Millipore, USA). After blocking, the membranes were incubated separately with specific primary antibodies overnight at 4°C, followed by incubation with goat anti-rabbit IgG or goat anti-mouse IgG antibodies (ZSGB-Bio, China) for 1 h at room temperature. Finally, immunoreactive proteins were visualized using a chemiluminescence detection system (FluorChem HD2, USA).

Statistical analysis
Data were analyzed by one way ANOVA or t-test, and statistical analyses were performed using SPSS 21.0 and GraphPad Prism 5.0 software. The comparisons with a P < 0.05 were considered statistically significant.

Bioinformatics analysis of target genes
By bioinformatics analysis, we explored the target genes of miR-20b-5p related to cell proliferation and cell cycle and found 10 potential targets of miR-20b-5p, namely, CyclinD1 (CCND1), E2F1, MAPK1, STAT3, R2b23b, RAB5BR, RABEP1, TAOK3, PPARDR and XIAP. To roughly screen the potential target genes of miR-20b-5p, agomir-20b-5p was transfected into MCF-7 cells and the gene expression was detected by qRT-PCR. We found that miR-20b-5p increased the expression of some of the Fig. 1 Screening of breast cancer stem cells (BCSCs)-associated miRNAs via miRNA networks. a In the miRNA-gene network, the rectangle in the figure represents upregulated miRNAs, the circle represents genes, and the line represents the regulatory relationship between miRNAs and genes; b The regulatory degrees of miRNAs in the miRNA-gene network; c In the miRNA-GO network, the rectangle in the figure represents miRNAs, the circle represents GO terms, and the straight line represents the regulatory relationship between miRNAs and GO terms. The more miRNAs regulating a GO term, the larger its area; d The regulatory degrees of miRNAs in the miRNAs-GO network potential target genes while inhibited that of the others (Fig. 2a). The dramatically increased expression of CCND1 and E2F1 was intriguing and unexpected, since miRNAs usually negatively regulate their target genes. Therefore, the contradictory preliminary results of qRT-PCR compelled us to choose CCND1 and E2F1 for subsequent validation and functional studies. The binding free energy between CCND1 and miR-20b-5p at the 5′ UTR and 3′ UTR was − 14.50 kcal/mol and − 12.10 kcal/ mol, respectively, lower than the average random free energy of CCND1 (0 kcal/mol). Likewise, the binding free energy between E2F1 and miR-20b-5p at the 5′ UTR and 3′ UTR was − 18.10 kcal/mol and − 19.00 kcal/ mol, respectively, lower than the average random free energy of E2F1 (− 3.60 kcal/mol) (Fig. 2b) (Supplementary Fig. 2).

Verification of target genes
A crucial step to clarifying the regulatory mechanism of miR-20b-5p was to determine its target genes. Three experiments were conducted as follows.
First, we compared the correlation of miR-20b-5p expression with CCND1 and E2F1 mRNA levels between BCSCs and MCF-7 cells or T47D cells. The CD44 + CD24 −/low subpopulation accounted for 2.58% of MCF-7 cells and 64.3% of induced MCF-7-CSCs (Fig. 4a). The expression of miR-20b-5p was 3.98 times and 2.21 times higher in BCSCs than in the control MCF-7 and T47D cells, respectively. Consistently, the mRNA expression levels of E2F1 and CCND1 were higher in BCSCs than in MCF-7 and T47D cells, although the difference was not statistically significant (Fig. 6a), thus confirming the regulatory relationship between miR-20b-5p and the two potential targets CCND1 and E2F1.
Based on all the three results above, it seemed that miR-20b-5p directly downregulated the targets of CCND1 and E2F1, but indirectly upregulated these two proteins.

Discussion
We made researches on the reports on the miR-20b family in different tumors. The expression level of miR-20b was reported to be increased in gastric cancer [16], breast cancer [17,18] and hepatocellular carcinoma [19]. Wang B et al. [25] found that miR-20b promoted proliferation, migration, invasion and tumorigenesis in esophageal cancer cells. It was also reported that the aberrant expression of miR-20b contributed to tumorigenesis and progression of breast cancer [26]. In the above literatures, miR-20b is shown to play an oncogenic role. However, miR-20b includes both miR-20b-3p and miR-20b-5p, which might have different biological functions. Only one report showed that overexpressed miR-20b-5p in cancer tissue and patient serum had a regulatory effect on the proliferation and migration of breast cancer cells [27]. In addition, there has been no exploration of the regulatory role of miR-20b in BCSCs. In the current study, we focused on the function of miR-20b-5p and found that miR-20b-5p acted as an oncogene in breast cancer and malignant BCSCs.
Cancer stem cells (CSCs) are capable of self-renewal and multidirectional differentiation [28]. Studies have shown that the molecular regulation of BCSCs is different from that of breast cancer cells [29]. BCSCs are associated with epithelial-mesenchymal transformation (EMT), miRNAs, the tumor microenvironment and other factors [30]. Our experimental results indicated that the upregulation of miR-20b-5p in MCF-7 and T47D cells significantly promoted the migration ability of breast cancer cells. Furthermore, overexpression and knockdown systems of miR-20b-5p in BCSCs were established. It was confirmed that miR-20b-5p promoted the proliferation of BCSCs and inhibited the apoptosis of BCSCs. As an integrated outcome of miR-20b-5p on breast cancer cells and stem cells, the xenografts of breast cancer cells were significantly enlarged by miR-20b-5p overexpression, which indicated that miR-20b-5p acted as an oncogene in vivo. Similarly, one study has demonstrated that miR-128-3p inhibits the stem-like cell features of BCSCs via inhibition of the Wnt signaling pathway by downregulating NEK2, which provides a new target for breast cancer treatment [31]. Then, another study indicates that CD24 and CD44 are cancer stem cells which can promote the development of breast cancer [32]. Many miRNAs have been shown to regulate the self-renewal and differentiation of CSCs, such as the let-7 miRNA family [33,34]. The let-7 miRNA family appears to play a substantial role in the CSC phenotype. In breast cancer, let-7 is found to be downregulated. In normal tissues, it regulates self-renewal, acting as a prodifferentiation miRNA, whereas in breast cancer it is repressed by the Wnt/β-catenin pathway. Therefore, we also compared its expression difference between normal tissues and tumor tissues (http://mirnamap.mbc.nctu. edu.tw/) in our study, and found that the expression of miR-20b-5p in normal tissues was also different. The expression level was high in the prostate, kidney, lung and other normal tissues, and was low in the liver, adipose, muscle. The highest expression level was in the thymus. No expression was observed in breast. ( Supplementary  Fig. 3).
Although several researches showed that miR-20b-5p acts as an oncogenic molecule in various tumors, we successfully screened BCSCs-associated miRNA profiles and identified miR-20b-5p as a key molecule in mediating BCSCs function. To our best knowledge, it's the first time that miR-20b-5p has been proved to have important effects on BCSCs. It's a novel regulation mechanism for the malignant progression of BCSCs, such as proliferation and apoptosis.
To further clarify the underlying mechanism of miR-20b-5p in the regulation of breast cancer cells and stem cells, we chose CCND1 and E2F1 as potential targets based on bioinformatics analysis and preliminary screening by qRT-PCR. CCND1, as an oncogene, is overexpressed in tumors and plays an important regulatory role in normal breast development, damage repair, maintenance of breast epithelial stem cell proliferation and self-renewal [35][36][37]. E2F1, a member of the E2F transcription factor family, regulates gene expression related to cell proliferation, differentiation and apoptosis and controls the cell cycle via a two-way regulatory mechanism. As a tumor suppressor gene or oncogene, E2F1 is closely related to tumor progression and drug resistance [38][39][40][41]. The literature shows that CCND1 and E2F1 are regulated by the miR-17-92 cluster, which is a miRNA cluster located on chromosome 13 that is composed of miR-17, miR-20a and miR-92a-1 [8,9]. This suggests that E2F1 is targeted by certain miRNAs and is involved in the regulation of tumor cells [42]. To confirm that miR-20b-5p regulated CCND1 and E2F1 expression, we performed three experiments, namely, a comparison between BCSCs and the control breast cancer cells, a dual-luciferase assay and western blot analysis after overexpression or knockdown of miR-20b-5p. The dual-luciferase reporting system showed that CCND1 and E2F1 were the real targets of miR-20b-5p. Our western blot analysis showed that antagomir-20b-5p significantly decreased the protein levels of both CCND1 and E2F1. The luciferase activity was inhibited by miRNA mimics, indicating that miR-20b-5p downregulated the targets. On the contrary, if the luciferase activity was enhanced by miRNA mimics, it may suggest that miR-20b-5p upregulated the targets. Additionally, western blot assay may not demonstrate a direct regulation, as it is a comprehensive reflection of direct and indirect regulation by the miRNA.
This seems to be a contradictory but interesting phenomenon. One speculation is that the observed upregulation of CCND1 and E2F1 might well be the side effect of the upregulation of miR-20b-5p. Theoretically a high level of miR-20b-5p may cause the occupation of the RNA-induced silencing complex (RISC) by this miRNA and thus RISCs are less available for other miR-NAs. Namely, the mature miRNA (the guide strand) is incorporated into RISCs, whereas the passage miRNA* strand can be loaded in the RISC as well or degraded [43]. The mature miRNA guides the Argonaute protein of the RISC to the complementary mRNA sequence on the target to repress its expression [44]. In this case, some genes targeted by other miRNAs might get released from their negative regulation and might get upregulated [45]. Another speculation is that the antagomirs or mimic will produce off-target effects [46]. Off-target phenomenon usually occurs when miR is overexpressed more than hundreds of folds. In our experiment, we used agomir to make miR overexpression less than 80 folds, which is a proper range and effectively avoids off-target effects. Therefore, the second speculation is out of our consideration.
Our results revealed an indirect upregulation of CCND1 and E2F1 by miR-20b-5p. There are likely to be unknown targets (X and Y) negatively regulated by miR-20b-5p. The unknown targets (X and Y) could inhibit CCND1 and E2F1 protein expression. Overall, miR-20b-5p regulated both CCND1 and E2F1 via bidirectional regulation, namely, direct downregulation and indirect upregulation. As a tug-of-war mechanism, the indirect promoting effect of miR-20b-5p on CCND1 and E2F1 may outweigh its direct downregulation of CCND1 and E2F1. However, it is still necessary to reveal the underlying mechanism, which will be explored in further studies. The hypothesis involving the effect of miR-20b-5p on CCND1 and E2F1 is shown in a regulatory network ( Supplementary Fig. 4).

Conclusion
In summary, this study shows that miR-20b-5p can promote tumorigenesis in both breast cancer cells and stem cells. We speculate that the underlying mechanism of miR-20b-5p contributing to the malignant progression of breast cancer is that miR-20b-5p overall upregulates both CCND1 and E2F1 via bidirectional regulation.