CircZFR functions as a sponge of miR-578 to promote breast cancer progression by regulating HIF1A expression

Background Breast cancer (BC) is the most common malignancy among women. Emerging studies have demonstrated that circular RNA (circRNA) zinc finger RNA binding protein (circZFR) serves as a crucial regulator in many human cancers. However, the role and mechanism of circZFR in BC tumorigenesis remain unclear. Methods The levels of circZFR, miR-578 and hypoxia-inducible factor 1α (HIF1A) were detected by quantitative real-time polymerase chain reaction (qRT-PCR) or western blot. Cell viability, colony formation, apoptosis, migration and invasion capacities in vitro were determined by using the Cell Counting Kit-8 (CCK-8), standard colony formation, flow cytometry and transwell assays, respectively. Glucose uptake, lactate product and adenosine triphosphate (ATP) levels of cells in vitro were measured using the commercial human assay kits. Targeted relationships among circZFR, miR-578 and HIF1A in BC cell lines were verified by dual-luciferase reporter and RNA pulldown assays. Animal studies were performed to assess the effect of circZFR on tumor growth in vivo. Results Our data indicated that circZFR was overexpressed in BC tissues and cells, and the increased circZFR level predicted poor prognosis of BC patients. CircZFR silencing or miR-578 overexpression repressed BC cell viability, colony formation, migration, invasion, and glycolysis and enhanced cell apoptosis in vitro. CircZFR silencing also hampered tumor growth in vivo. Mechanistically, circZFR acted as a sponge of miR-578, and circZFR silencing hindered BC cell malignant behaviors by miR-578. HIF1A was a functional target of miR-578 in regulating BC cell viability, colony formation, migration, invasion, glycolysis and apoptosis in vitro. Furthermore, circZFR modulated HIF1A expression through sponging miR-578. Conclusion Our findings first identified that the silencing of circZFR suppressed BC malignant progression in vitro via the regulation of the miR-578/HIF1A axis, providing evidence for the crucial involvement of circZFR in BC pathogenesis.


Background
Breast cancer (BC) remains the most commonly diagnosed cancer and the leading cause of cancer-associated death among females in 2018 [1]. Although the therapeutic methods have greatly improved over the past two decades, effective treatment against metastatic BC is still limited [2,3]. Therefore, a deeper understanding of what drives this disease is the first step to design innovative interventions.
Circular RNAs (circRNAs) are covalently closed, endogenous RNAs that have crucial non-coding functions in human physiologic and pathologic processes [4]. Work in biological functions has demonstrated the roles of circRNAs as microRNA (miRNA) sponges [5]. Dysregulation of circRNAs has recently implicated in the pathogenesis of BC [6]. For example, Yuan et al. uncovered that hsa_circ_0068033, a down-regulated circRNA, exerted a repressive impact on BC malignant progression via sequestering miR-659 [7]. Cao and colleagues demonstrated that hsa_circ_0087784 functioned as a potential promoter in BC development through sponging miR-487a [8]. Xu et al. reported that circTADA2As attenuated BC progression and metastasis by the regulation of miR-203a-3p [9]. Moreover, hsa_circ_001783 and circABCB10 were reported as oncogenic regulators in BC through functioning as specific miRNA sponges [10,11]. As for circRNA zinc finger RNA binding protein (circZFR, hsa_circ_0072088), it has been identified as an oncogenic modulator in many human cancers, such as renal carcinoma, bladder cancer and non-small cell lung cancer [12][13][14]. Nevertheless, the biological roles of circ-ZFR in BC tumorigenesis remain largely unknown.
MiRNAs modulate gene expression but are frequently dysregulated in human cancers, including BC [15]. Danza et al. reported that miR-578 was underexpressed in BRCA-BC and it regulated tumor angiogenesis [16]. However, the precise function of miR-578 in BC progression is still undefined.
Here, we undertook to examine the biological effect of circZFR in the malignant progression of BC in vitro and in vivo. We identified that circZFR, a prominently up-regulated circRNA in BC, controlled BC progression in vitro via targeting the miR-578/hypoxia-inducible factor 1α (HIF1A) axis.

Clinical tissues and cells
In this study, we enrolled 70 BC patients admitted to The First Affiliated Hospital of Zhengzhou University from April 2013 to June 2014. The clinicopathological features of these patients were provided in Table 1. Seventy pairs of primary tumor tissues and matched healthy breast tissues were collected and stored at − 80 °C to detect the expression levels of circZFR, miR-578 and HIF1A (Accession: NM_181054.3). These patients were followed-up for at least 60 months and the follow-up information was obtained by telephone calls every 3 months. The study was approved by the Ethics Committee of The First Affiliated Hospital of Zhengzhou University, and written informed consent was provided by all participants. BC cell lines MCF7 (ATCC ® HTB-22), BT-549 (ATCC ® HTB-122), MDA-MB-231 (ATCC ® HTB-26) and MDA-MB-453 (ATCC ® HTB-131, American Type Collection Culture, ATCC, Manassas, VA, USA) were cultured in RPMI-1640 medium, supplemented with 10% fetal calf serum (FCS) and 1% antibiotic solution (all from HyClone, Logan, UT, USA). The immortalized MCF10A cell line (ATCC ® CRL-10317, ATCC) was maintained in Dulbecco's modified Eagle's medium/Nutrient Mixture F-12 (DMEM/F-12) with 10% FCS, 20 ng/mL epidermal growth factor, 0.5 μg/mL hydrocortisone and 5 μg/mL insulin (all from HyClone). All cells were cultured in a 5% CO 2 incubator at 37 °C.

Cell viability, colony formation and apoptosis assays
sh-circZFR-infected or sh-NC-transduced cell lines were transfected with or without anti-NC or anti-miR-578, and MCF7 and BT-549 cell lines were carried out the indicated transfections, followed by the incubation for 24 h at 37 °C. Cell viability assay was done using the Cell Counting Kit-8 (CCK-8, Abcam, Cambridge, UK) assay as per the manufacturer's guidance. Cell colony formation was tested using standard colony formation protocols as previously reported [19]. Cell apoptosis measurement was performed by flow cytometry using double-staining with fluorescein isothiocyanate (FITC)-labeled Annexin V and propidium iodide (PI) based on the directions of manufacturers (Invitrogen). 10,000 events were analyzed by a flow cytometer, and the apoptotic cells were determined by calculating the sum of early (Annexin V + /PI − ) and late (Annexin V + /PI + ) apoptotic cells. All experiments were done in triplicate.

Transwell migration and invasion assays
Cell migration and invasion were detected by the transwell assay using modified Boyden chambers in 24-transwell plates (8 μm pores, Corning, Amsterdam, The Netherlands). Chambers of cell invasion assays consisted of Matrigel-precoated membrane inserts (Corning). After transfection, BC cell lines were seeded into the top chamber of a 24-transwell insert, and medium containing 15% FCS was used as a chemoattractant in the lower chamber. 24 h later, the well was stained with 1% crystal violet (Solarbio), and the number of cells that had migrated or invaded to the basal side of the membrane was counted under a microscope (Nikon, Shinagawa, Tokyo, Japan) at 100× magnification. Each experiment was performed in triplicate.

Measurement of glucose uptake, lactate product and adenosine triphosphate (ATP) level
The Colorimetric Glucose Uptake Assay Kit, L-Lactate Assay Kit and ATP Assay Kit (all from Abcam) were used to determine the levels of glucose uptake, lactate product and ATP, according to the recommendations of manufacturers. Briefly, the lysates of transfected cells were prepared using the Assay buffer and incubated with standard protocols for the indicated time point, followed by the measurement of the absorbance with a microplate reader (Invitrogen) at OD 412 nm for glucose uptake, 450 nm for lactate product and 570 nm for ATP level. All assays were done in triplicate.

Animal studies
The xenograft models were constructed to assess the role of circZFR on tumor growth in vivo. Animal studies were implemented in accordance with a protocol approved by the Ethics Committee on Animal Use and Care of the First Affiliated Hospital of Zhengzhou University. Twelve 6-8 weeks BALB/c female mice (Shanghai Animal Laboratory Center, Shanghai, China) were used and randomly divided into two groups (n = 6 per group): sh-NC and sh-circZFR. Approximately 5 × 10 6 sh-NC-infected or sh-circZFR-transduced BT-549 cell line was subcutaneously injected into the dorsal flank of the mice. One week later, tumor size was measured every week with callipers and tumor volume was calculated using the following formula: volume = 0.5 × length × (width) 2 . All mice were sacrificed at 35 days after implantation, and tumor tissues were collected for weight.

Bioinformatics, dual-luciferase reporter and RNA pulldown assays
Analysis for the targeted miRNAs of circZFR was performed using the online web CircInteractome (https :// circi ntera ctome .nia.nih.gov/index .html) and CircBANK (http://www.circb ank.cn/searc hCirc .html). The putative targets of miR-578 were predicted by TargetScan v.7 online software (http://www.targe tscan .org/vert_71/). Targeted relationships among circZFR, miR-578 and HIF1A were confirmed by dual-luciferase reporter and RNA pulldown assays. In dual-luciferase assays, the fragment of circZFR containing the miR-578-binding sites and HIF1A 3′UTR were individually cloned into pmirGLO vector (Promega, Madison, WI, USA) to construct corresponding wild-type reporters (circZFR-WT and HIF1A-3′UTR-WT). The TaKaRa MutanBest Kit was used to construct the corresponding mutations (circZFR-MUT and HIF1A-3′UTR-MUT) as per the insturctions of manufacturers (TaKaRa, Dalian, China). Each reporter construct (200 ng) and 30 nM of miR-578 mimic or miR-NC mimic were cotransfected into BC cell lines. Cell line extracts were prepared with RIPA buffer (TaKaRa) 24 h post-transfection, and the ratio of Renilla to firefly luciferase was detected using the Promega Dualluciferase Assay. In RNA pulldown assays, cell lysates were incubated with the biotin-labeled miR-578 mimic (Bio-miR-578, Ribobio) or nontarget control sequence (Bio-NC, Ribobio) for 4 h at 4 °C before adding to the streptavidin beads (Sigma-Aldrich) for 2 h. The beads were collected, and total RNA was extracted for circ-ZFR quantification. Each experiment was performed in triplicate.

Western blot for HIF1A expression
Western blot was used to determine the expression of HIF1A using standard protocols. The preparation of cell line extracts was done using RIPA buffer with proteinase inhibitors (Roche, Charente, France). Proteins (50 μg) were resolved on a 10% SDS-polyacrylamide gel, electrophoretically blotted onto nitrocellulose membranes (GE Healthcare, Little Chalfont, UK) and probed with antibody against HIF1A (ab51608, 1 μg/mL, Abcam) or GAPDH (MA5-15738, 0.5 μg/mL, Invitrogen). Following the incubation with horseradish peroxidase-coupled IgG secondary antibody (ab97051, 0.1 μg/mL, Abcam), the signals were visualized by Cheniluminescence (GE Healthcare) as recommended by the manufacturers. All experiments were done in triplicate.

Statistical analysis
Data were shown as the mean ± standard deviation from at least three independent assays. Pairwise comparisons were done using a two-tailed Student's t test, Mann-Whitney U test or analysis of variance (ANOVA) with SPSS version 19.0 software (SPSS, Chicago, IL, USA). For survival analysis, the Kaplan-Meier survival curve and log-rank test were used. Correlations among circZFR, miR-578 and HIF1A expression levels in BC tissues were determined by the Spearman correlation test. All tests were considered statistically significant at p-value < 0.05.

Overexpression of circZFR predicted poor prognosis of BC patients
As demonstrated by qRT-PCR, circZFR was significantly up-regulated in BC tissues and cell lines compared with their counterparts (Fig. 1a, b). To determine its clinical relevance, we preliminarily examined the link between circZFR level and the prognosis of BC patients. Kaplan-Meier survival curves showed that the patients in low circZFR level group had a longer survival time than those in high circZFR group (Fig. 1c). Additionally, circZFR expression was remarkably associated with the distant metastasis and TNM stage of these patients ( Table 1).

Silencing of circZFR hindered BC cell viability, colony formation and enhanced apoptosis in vitro and weakened tumor growth in vivo
To study the biological role of circZFR in BC progression, the loss-of-function experiments were carried out using shRNAs against circZFR (sh-circZFR#1, sh-circZFR#2 and sh-circZFR#3). In contrast to the negative control, the transfection of the three shRNAs prominently reduced circZFR expression in both MCF7 and BT-549 cell lines (Fig. 2a, b). Notably, sh-circZFR#1 (also named sh-circZFR) caused the most significant down-regulation in circZFR expression, so we used it for further analyses. Functional experiments data revealed that the silencing of circZFR led to a striking inhibition in cell viability (Fig. 2c), colony formation (Fig. 2d, e), as well as a strong promotion in cell apoptosis (Fig. 2f, g).
To determine whether circZFR regulated BC tumor development in vivo, we performed the xenograft model assays. When we infected BT-549 cell line with sh-circ-ZFR, tumor growth was remarkably blocked compared with the sh-NC control (Fig. 2h, i).

Silencing of circZFR suppressed BC cell migration, invasion and glycolysis
We also asked whether circZFR regulated BC cell migration, invasion and glycolysis in vitro. Transwell assays showed that in comparison to the control group, cell migration (Fig. 3a) and invasion (Fig. 3b) were significantly repressed by circZFR knockdown. Moreover, in both cell lines, circZFR silencing resulted in decreased glucose uptake (Fig. 3c), lactate product (Fig. 3d) and ATP level (Fig. 3e), demonstrating the suppression of circZFR knockdown on cell glycolysis.

CircZFR directly interacted with miR-578 by binding to miR-578
To further understand the role of circZFR in BC pathogenesis, we performed a detailed analysis for its targeted miRNAs. The two online algorithms CircInteractome and CircBANK collectively revealed that circZFR harbored a putative complementary sequence for miR-578, miR-944 and miR-532-3p (Fig. 4a). The data of qRT-PCR showed that the silencing of circZFR led to a striking overexpression in miR-578 and miR-532-3p levels, but miR-944 expression was not affected by circZFR knockdown in the two BC cell lines (Fig. 4b, c). Previous work had reported that circZFR regulated colorectal cancer progression through acting as a sponge of miR-532-3p [20]. So, we aimed to identify whether miR-578 was a molecular mediator of circZFR in BC progression. By contrast, miR-578 was prominently up-regulated in the two BC cell lines transfected with miR-578 mimic (Fig. 4d). We then cloned circZFR fragment containing the miR-578-binding sites into a luciferase plasmid and mutated the miR-578-binding sites (Fig. 4e). The elevated miR-578 expression significantly reduced the activity of circZFR wild-type reporter (circZFR-WT) (Fig. 4f). However, the mutation of the target sites (circZFR-MUT) completely abrogated the effect of miR-578 on reporter gene expression (Fig. 4f), indicating the validity of the target sequence for interaction. Additionally, RNA pulldown assays revealed that the enrichment level of circZFR was remarkably elevated by Bio-miR-578 in both cell lines (Fig. 4g).

Overexpression of miR-578 restrained BC cell viability, colony formation, migration, invasion and glycolysis and promoted apoptosis in vitro
In BC tissues, miR-578 expression was significantly decreased compared to the normal control (Fig. 5a), and it was inversely correlated with circZFR level (Fig. 5b). Moreover, miR-578 level was lower in BC cell lines than that of control (Fig. 5c). Subsequently, we analyzed the biological effect of miR-578 on BC progression. In contrast to the control group, the increased expression of miR-578 induced a distinct repression in cell viability (Fig. 5d) and colony formation (Fig. 5e), and a strong promotion in cell apoptosis (Fig. 5f), as well as a striking reduction in cell migration (Fig. 5G), invasion (Fig. 5H) and glycolysis ( Fig. 5i-k).

CircZFR modulated HIF1A expression via acting as a sponge of miR-578
Using the software TargetScan, a predicted miR-578-binding sequence (ACA AGA A) was identified within the 3′UTR of HIF1A (Fig. 6a). When we cloned the 3′UTR fragment containing the putative miR-578-binding sites downstream of a luciferase coding sequence, the cotransfection of the luciferase reporter (HIF1A-3′UTR-WT) and miR-578 mimic into the two BC cell lines produced lower luciferase activity than in cell lines cotransfected with the miR-NC control (Fig. 7b,  c). However, the mutation of the target sequence (HIF1A-3′UTR-MUT) prominently abolished the suppression of miR-578 (Fig. 7b, c). By contrast, HIF1A mRNA and protein levels were significantly reduced by miR-578 overexpression in the two BC cell lines (Fig. 7d, e). These data together pointed that HIF1A in BC cell lines was directly targeted and inhibited by miR-578.  invasion (a, b), corresponding assay kits for glucose uptake, lactate product and ATP level (c-e) in sh-NC-infected or sh-circZFR-transduced MCF7 and BT-549 cell lines. *P < 0.05 We then examine whether circZFR influenced HIF1A expression in BC cell lines. As expected, in comparison to their counterparts, HIF1A expression was remarkably down-regulated by circZFR silencing at both mRNA and protein levels in the two BC cell lines, and this effect was significantly abrogated by anti-miR-578 introduction (Fig. 7f, g). Additionally, qRT-PCR data showed a striking up-regulation of HIF1A mRNA level in BC tissues (Fig. 7h). Furthermore, in BC tissues, HIF1A mRNA expression was positively correlated with circZFR expression and negatively correlated with miR-578 level (Fig. 7i).

HIF1A was a functional target of miR-578 in regulating BC cell viability, colony formation, migration, invasion, glycolysis and apoptosis in vitro
To provide further insight into the link between miR-578 and HIF1A in BC progression, we elevated HIF1A expression using a recombinant overexpressing plasmid in miR-578 mimic-transfected BC cell lines. As a result, HIF1A protein level was prominently increased by the overexpressing plasmid in the two cell lines (Fig. 8a). Functional experiments results revealed that the up-regulation of HIF1A dramatically abolished miR-578 overexpression-mediated anti-viability (Fig. 8b), pro-apoptosis CircZFR directly interacted with miR-578 by binding to miR-578. a Venn diagrams representing the putative miRNAs that bind to circZFR identified by CircInteractome and CircBANK online algorithms. b, c The levels of miR-578, miR-944 and miR-532-3p by qRT-PCR in sh-NC-infected or sh-circZFR-transduced MCF7 and BT-549 cell lines. d qRT-PCR for miR-578 expression in the two BC cell lines transfected with miR-NC mimic or miR-578 mimic. e Schematic of the miR-578-binding sites within circZFR and the mutation of the seed sequence. f Relative luciferase activity in MCF7 and BT-549 cell lines cotransfected with circZFR-WT or circZFR-MUT and miR-NC mimic or miR-578 mimic. g qRT-PCR for circZFR level in cell lysates incubated with Bio-NC or Bio-miR-578 and streptavidin beads. *P < 0.05 (Fig. 8c), anti-migration (Fig. 8d), anti-invasion (Fig. 8e) and anti-glycolysis ( Fig. 8f-h).

Discussion
To date, the emerging links between circRNAs and BC progression have opened up a novel field for cancer diagnosis and treatment [21,22]. In the meantime, understanding the molecular basis underlying the actions of circRNAs has been still challenging. In this study, we identified the biological role of circZFR in BC progression in vitro and in vivo and investigated the mechanisms governing it.
Here, we firstly demonstrated that circZFR was overexpressed in BC, and the elevated expression of circ-ZFR was associated with the poor prognosis of these patients. By using loss-of-function in vitro and in vivo analyses, we were first to uncover that the silencing of circZFR performed a suppressive effect in BC progression. Previous reports had highlighted the potential oncogenic role of circZFR in several other malignancies, such as hepatocellular carcinoma, papillary thyroid cancer and non-small cell lung cancer [14,23,24]. Conversely, circZFR was reported as a tumor suppressor in colorectal cancer and gastric cancer [20,25]. These contradictory findings may attribute to different type tumor or complex tumor microenvironment.
Using the online algorithms, we first identified that circZFR sequestered miR-578 through sponging miR-578 in BC cell lines. Ji   and BT-549 cell lines were transfected with miR-NC mimic or miR-578 mimic, followed by the determination of cell viability by CCK-8 assay (d), colony formation using a standard colony formation assay (e), cell apoptosis by flow cytometry (f), cell migration (g) and invasion (h) by transwell assay, glucose uptake, lactate product and ATP level using assay kits (i-k). *P < 0.05 associated with the tumor cell apoptosis [27]. As it has been reported [16], our data validated the down-regulation of miR-578 expression in BC tissues and cells. Moreover, we first identified that miR-578 overexpression restrained BC cell malignant behaviors in vitro. More importantly, for the first time, we substantiated that circZFR knockdown hampered BC progression in vitro by miR-578.
Using TargetScan software, we identified that miR-578 in BC cell lines directly targeted and inhibited HIF1A. HIF1A, a transcriptional regulator in response to intratumoral hypoxia [28,29], contributes to BC metastasis and malignant progression [30][31][32]. We also uncovered that the up-regulation of miR-578 hampered BC malignant progression via down-regulating HIF1A in vitro. Previous studies had reported that several Fig. 6 CircZFR knockdown repressed BC malignant progression by miR-578 in vitro. sh-NC-infected or sh-circZFR-transduced MCF7 and BT-549 cell lines were transiently transfected with anti-NC or anti-miR-578. a CCK-8 assay for cell viability. b A standard colony formation assay for cell colony formation. c Flow cytometry for cell apoptosis. d, e Transwell assay for cell migration and invasion. f-h Glucose uptake, lactate product and ATP level using assay kits. *P < 0.05 other miRNAs, such as miR-18a and miR-497, exerted an anti-tumor activity in BC through targeting HIF1A [33,34]. Furthermore, our data first illuminated the role of circZFR as a sponge of miR-578 to mediate HIF1A expression in BC cell lines. The findings by Liang et al. underscored that circRNA circDENND4C contributed to BC cell proliferation under hypoxia via regulating HIF1A [35].
In conclusion, our present study demonstrated that circZFR was overexpressed in BC, and the silencing of CircZFR modulated HIF1A expression via sponging miR-578. a Schematic of the putative miR-578-binding sequence and mutated the target sequence. b, c Relative luciferase activity in MCF7 and BT-549 cell lines cotransfected with HIF1A-3′UTR-WT or HIF1A-3′UTR-MUT and miR-578 mimic or miR-NC mimic. HIF1A mRNA and protein levels by qRT-PCR and western blot in MCF7 and BT-549 cell lines transfected with miR-578 mimic or miR-NC mimic (d, e), sh-NC-infected or sh-circZFR-transduced MCF7 and BT-549 cell lines transfected with anti-NC or anti-miR-578 (f, g). h Relative HIF1A mRNA expression by qRT-PCR in 70 pairs of tumor tissues and adjacent normal tissues. i Correlation between HIF1A expression with circZFR or miR-578 level in BC tissues using Spearman test. *P < 0.05 Fig. 8 The repression of miR-578 up-regulation on BC progression in vitro was mediated by HIF1A. a HIF1A protein level by western blot in MCF7 and BT-549 cell lines transfected with vector or HIF1A. MCF7 and BT-549 cell lines were transfected with miR-NC mimic + vector, miR-578 mimic + vector or miR-578 mimic + HIF1A, followed by the determination of cell viability by CCK-8 assay (b), cell apoptosis by flow cytometry (c), cell migration and invasion by transwell assay (d, e), glucose uptake, lactate product and ATP level using assay kits (f-h). vector: negative control plasmid, HIF1A: recombinant HIF1A overexpressing plasmid. *P < 0.05