Skip to main content
  • Primary research
  • Open access
  • Published:

FOXA1-induced circOSBPL10 potentiates cervical cancer cell proliferation and migration through miR-1179/UBE2Q1 axis

Abstract

Background

Recently, extensive evidence has clarified the crucial role of circular RNAs (circRNAs) as a pro-tumor or anti-cancer participant in human malignancies. A new circRNA derived from oxysterol binding protein like 10 (OSBPL10) (circOSBPL10) has not been researched in cervical cancer (CC) yet.

Methods

The expression of molecules was analyzed by RT-qPCR or western blot. Several functional assays were applied to explore the biological influence of circOSBPL10 on CC. The interaction between RNAs was estimated via luciferase reporter, RNA immunoprecipitation and RNA pull-down assays.

Results

CircOSBPL10 characterized with cyclic structure was revealed to possess elevated expression in CC cells. CircOSBPL10 downregulation elicited suppressive impacts on CC cell proliferation and migration. Interestingly, circOSBPL10 regulated CC progression by interacting with microRNA-1179 (miR-1179). Moreover, ubiquitin conjugating enzyme E2 Q1 (UBE2Q1) targeted by miR-1179 was positively regulated by circOSBPL10 in CC. Furthermore, enhanced UBE2Q1 expression or suppressed miR-1179 level countervailed the repressive effect of circOSBPL10 depletion on the malignant phenotypes of CC cells. Moreover, forkhead box A1 (FOXA1) was confirmed to induce circOSBPL10 expression in CC cells.

Conclusions

FOXA1-induced circOSBPL10 facilitates CC progression through miR-1179/UBE2Q1 axis, highlighting a strong potential for circOSBPL10 to serve as a promising therapeutic target in CC.

Background

As a frequent type of human gynecological malignancies worldwide, cervical cancer (CC) is depicted as one of the dominating causes contributing to cancer-associated death in women [1, 2]. It is estimated that that nearly 500,000 new cases are diagnosed with CC annually [3]. In China, CC is also regarded as one of the most prevalent lethal tumors. Over the past few years, in spite of the application of human papillomavirus (HPV) vaccine in treatment, CC is still a major stumbling block for female health [1, 4]. The majority of patients at an early stage of CC are likely to be cured through surgery [5], whereas no or limited efficient therapeutic approaches for those at the advanced stages [6]. In order to make advances in the treatment of CC, researchers have focused on exploring and developing tumor-particular biomarkers for CC [7]. For example, melatonin was identified as a new adjuvant agent in treating patients with CC [8], so was curcumin [9]. However, more efforts should be made in developing new effective therapeutic strategies for CC. Hence, it is imperative to make in-depth exploration of the underlying molecular mechanisms in CC.

Non-coding RNAs (ncRNAs) have been commonly considered as the potential key modulators in gene regulation to impact on tumor development [10, 11], including genetic and epigenetic manners [12]. In recent years, ncRNAs have been indicated as clinical biomarkers in diverse diseases [13, 14], including cancer [15, 16]. As a member of ncRNAs, circular RNAs (circRNAs) own a covalently closed structure which are tolerant to RNase R-mediated degradation [17, 18]. Increasing analyses have suggested the abnormal expression of circRNAs in various cancers [19, 20], including CC [21, 22]. Emerging researches have testified the implication of circRNAs in tumorigenesis and progression via regulation of different biological processes, which includes cell proliferation, migration and invasion [23,24,25]. Emerged as a new circRNA, circOSBPL10 (circbase ID: hsa_circ_0064669) is derived from back-splicing of OSBPL10 mRNA (messenger RNA). To our knowledge, the critical regulatory mechanism of circOSBPL10 has not been investigated in CC yet.

In this study, the main purpose was to decipher the potential regulatory role of circOSBPL10 in CC. Data from a series of assays uncovered that FOXA1-induced upregulation of circOSBPL10 contributes to CC progression via miR-1179/UBE2Q1 axis, revealing that circOSBPL10 might be a hopeful biomarker for CC.

Methods

Cell culture

Human normal cervical epithelial cells (H8) and human CC cells (C33A, CaSki, HeLa and SiHa) were bought from Chinese Academy of Sciences (Shanghai, China). Cells were cultured with Dulbecco’s Modified Eagle Medium (DMEM; Invitrogen, Carlsbad, CA, USA) adding 10% fetal bovine serum (FBS; Invitrogen), 1% penicillin/streptomycin (Sigma-Aldrich, Milan, Italy), and then incubated in an incubator of 5% CO2 at 37 °C.

Cell transfection

HeLa and SiHa cells were transfected with specific short hairpin RNAs (shRNAs) against circOSBPL10 (sh-circOSBPL10#1#2), FOXA1 (sh-FOXA1#1#2) and their corresponding negative control (NC) sh-NC, as well as pcDNA3.1/circOSBPL10, pcDNA3.1/UBE2Q1, pcDNA3.1/FOXA1 and empty pcDNA3.1 (±) circRNA Mini vector, empty pcDNA3.1 vector, severally. The miR-1179 mimics, miR-1179 inhibitor, NC mimics and NC inhibitor were synthesized by GenePharma (Shanghai, China). Transfection experiments were executed by Lipofectamine 2000 (Invitrogen).

Real-time quantitative polymerase chain reaction (RT-qPCR)

Total RNA of cells was isolated using TRIzol reagent, followed by cDNA (complementary DNA) synthesis with Reverse Transcription Kit (Invitrogen). RT-qPCR was measured by SYBR-Green Real-Time PCR Kit (Takara, Tokyo, Japan) operated on Bio-Rad CFX96 system (Takara). Relative expression level was calculated utilizing 2−ΔΔCt method with normalization to glyceraldehyde-3-phosphate dehydrogenase (GAPDH) or U6. The sequences of primers were presented in Additional file 1: Table S1.

Cell counting kit-8 (CCK-8)

In short, 1 × 103 cells were plated into a 96-well plate. After incubation for specific times (24, 48, 72, 96 h), cells were processed with 10 μL of CCK-8 reagent for additional 4 h. Absorbance at 450 nm was measured via a microplate reader (Olympus, Tokyo, Japan).

Colony formation assay

Transfected cells (1 × 103) were first coated into 6-well plates. After 2 weeks of incubation, cells were rinsed with phosphate buffer saline (PBS; Sigma-Aldrich, San Francisco, USA), fixed in methanol (Sigma-Aldrich) and dyed using crystal violet (Sigma-Aldrich). The visible colonies were counted manually.

ActinomycinD (ActD) and RNase R treatment

To block transcription, 2 mg/ml Actinomycin D (ActD; Sigma-Aldrich) or dimethylsulphoxide (DMSO; Sigma-Aldrich) was added into culture medium. Total RNA was cultivated with or without 3 U/μg of RNase R (Epicentre Technologies, Madison, WI, USA) for 30 min. After treatment with ActD or RNase R, RT-qPCR was applied for determining the expression levels of circOSBPL10 and OSBPL10 mRNA.

Nucleic acid electrophoresis

Convergent primers and divergent primers were designed to amplify OSBPL10 mRNA or circOSBPL10. The level of circOSBPL10 in PCR products from cDNA and genomic DNA was examined by agarose gels with TE (Tris-ethylene diamine tetraacetic acid) buffer from Thermo Scientific (Waltham, MA, USA).

Terminal deoxynucleotidyl transferase-mediated dUTP nick-end labeling (TUNEL) assay

Apoptosis transfected SiHa and HeLa cells were assessed utilizing TUNEL Apoptosis Kit (Invitrogen). 4′,6-diamidino-2-phenylindole (DAPI; Sigma-Aldrich) was employed to dye above cells. The percentage of positive stained cells was observed by fluorescence microscopy (Olympus) and then analyzed.

Flow cytometry analysis

Cell apoptosis analysis was performed via Cell Apoptosis Analysis Kit (Takara). After incubation in 6-well plates, SiHa and HeLa cells were rinsed with PBS and resuspended in binding buffer. Followed by fixation with 70% ice-cold ethanol (Sigma-Aldrich), cells were double-stained by Annexin V-fluorescein isothiocyanate and propidium iodide. Last, cell apoptosis rate was detected by Flow Cytometer (Becton–Dickinson, MA, USA).

Migration assay

Cell migration abilities were examined using transwell chambers (Corning, NY, USA). Transfected cells with serum-free medium were placed into top compartment, while medium with 10% FBS was added into the lower compartment. 48 h later, cells in the lower chamber were immobilized and dyed in methanol and crystal violet, separately. Then migratory cells were then counted in five random chosen fields under a microscope (Olympus).

Wound healing

SiHa and HeLa cells were added in 6-well plates for cultivation. When cell confluence was 80%, scratches were produced in cell layer using sterile pipette tip. Afterward, cells were cleaned using PBS and incubated in a culture medium for 24 h. Images of migrating cells were detected at last.

Western blot

Total protein was reaped in Radio Immunoprecipitation Assay (RIPA) lysis buffer (Beyotime, Shanghai, China). Bicinchoninic acid (BCA) kit (Beyotime) determined protein concentrations. Proteins were separated through sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) and moved onto polyvinylidene fluoride (PVDF) membranes. The membranes were sealed with 5% skimmed milk and cultivated with primary antibodies against UBE2Q1 (orb77618, Biorbyt, San Francisco, California, USA) and GAPDH (ab8245, Abcam, Cambridge, UK). Secondary antibody was added for incubating for 1 h. GAPDH was an internal control. Proteins quantities were evaluated via chemiluminescence detection system.

Luciferase reporter assay

The wild-type (WT) and mutant (Mut) binding sites of miR-1179 in circOSBPL10 or UBE2Q1 3′UTR was sub-cloned into pmirGLO dual-luciferase vector to construct circOSBPL10-WT/Mut or UBE2Q1-WT/Mut. And plasmids were co-transfected with miR-1179 mimics or NC mimics into SiHa and HeLa cells, respectively. The pGL3-OSBPL10 promoter was co-transfected with pcDNA3.1/FOXA1 or pcDNA3.1 into cells. Dual-Luciferase Reporter Assay System (Promega, USA) detected the luciferase activity.

Subcellular fractionation

Fractions of cytoplasmic and nuclear were separated using NE-PERâ„¢ Nuclear and Cytoplasmic Extraction Reagents (Invitrogen) and gathered by RNeasy Midi Kit (Qiagen, Hilden, Germany) to determine the cellular localization of circOSBPL10. RT-qPCR was used to examine the levels of circOSBPL10, U6 (nuclear control) and GAPDH (cytoplasmic control).

RNA pull-down

Briefly, cell lysates were incubated with biotinylated RNA including Bio-miR-1179-WT, Bio-miR-1179-Mut and Bio-NC. Moreover, M-280 streptavidin magnetic beads (Sigma-Aldrich) were added to co-culture for 48 h. The relative enrichment of RNAs pulled down in each group were assayed by RT-qPCR.

RNA immunoprecipitation (RIP)

RIP assays were progressed with Magna RIPTM RNA-Binding Protein Immunoprecipitation Kit (Millipore, Bedford, USA). SiHa and HeLa cells were lysed with RIP lysis buffer, followed by incubation with magnetic beads conjugated with anti-Ago2 (Millipore) or anti-IgG (Millipore). The RT-qPCR was performed to evaluate the expression levels of circOSBPL10, miR-1179 and UBE2Q1 in the precipitates.

Chromatin immunoprecipitation (ChIP)

Via Magna ChIP Kit (Millipore), ChIP experiment was achieved. DNA in cell lysates was interrupted into 200-300-bp chromatin fragments by ultrasound. After that, lysates were subjected to immunoprecipitation with anti-FOXA1 or anti-IgG (negative control group). The precipitated DNA fragments were detected by RT-qPCR.

Statistical analysis

Experimental data from at least three independent experiments were shown as mean ± standard deviation (SD). Statistical analysis was performed using GraphPad Prism 7.0 software (Graph Pad, La Jolla, CA, USA). Significance in differences between 2 or more groups was analyzed via student’s t test or one-way analysis of variance (ANOVA). P < 0.05 had statistical significance in requirements.

Results

CircOSBPL10 is highly expressed in CC and its depletion impedes CC cell proliferation and migration

To study the cellular function of circOSBPL10 in CC, we first applied RT-qPCR analysis and unveiled a marked elevation of circOSBPL10 expression in CC cell lines compared with H8 cells (Fig. 1a). Then, nucleic acid electrophoresis manifested that in SiHa and HeLa cells, divergent primers could produce circOSBPL10 from cDNA but not from genomic DNA (gDNA), while convergent primers amplified linear OSBPL10 from both cDNA and gDNA (Fig. 1b). Besides, OSBPL10 mRNA was greatly degraded by ActD whereas circOSBPL10 exhibited as resistant to ActD (Fig. 1c). Additionally, OSBPL10 expression was dramatically reduced whereas circOSBPL10 expression demonstrated no evident change after SiHa and HeLa cells were treated with RNase R (Fig. 1d). Then, we verified that circOSBPL10 expression was lowered in two CC cells after transfection with sh-circOSBPL10#1/2, while those with sh-circOSBPL10#1 showed higher silencing efficiency (Fig. 1e). Subsequently, cell proliferation assays depicted a notably weakened proliferation ability of SiHa and HeLa cells under circOSBPL10 silence (Fig. 1f, g). Moreover, cell apoptosis capability was proved to be facilitated after silencing circOSBPL10 in SiHa and HeLa cells (Fig. 1h, i). Furthermore, it was uncovered that circOSBPL10 deficiency gave rise to attenuated migration ability of SiHa and HeLa cells (Fig. 1j, k). Taken together, circOSBPL10 is expressed at high levels in CC and knockdown of it impairs malignant behaviors in CC cells.

Fig. 1
figure 1

Circular RNA circOSBPL10 was highly expressed in CC and knockdown of it suppressed CC progression. a CircOSBPL10 expression was detected by RT-qPCR in CC cell lines H8 cells. b It was delineated by nucleic acid electrophoresis analysis that divergent primers amplified circOSBPL10 from cDNA, but not from gDNA. GAPDH was a negative control. c The resistance of circOSBPL10 and OSBPL10 mRNA to ActD in SiHa and HeLa cells was analyzed by RT-qPCR. d RT-qPCR assay was conducted to determine the abundance of circOSBPL10 and linear OSBPL10 mRNA in SiHa and HeLa cells treated with RNase R (normalized to mock treatment). e RT-qPCR was utilized to analyze the efficacy of circOSBPL10 knockdown in SiHa and HeLa cells. f, g The proliferation ability of SiHa and HeLa cells transfected with sh-circOSBPL10#1 or sh-NC was evaluated via CCK-8 and colony formation. h, i Cell apoptosis ability in transfected cells was measured by TUNEL assay and flow cytometry analysis. j, k Transwell and wound healing assays were conducted to analyze the migration of transfected cells. *P < 0.05, **P < 0.01

CircOSBPL10 sponges miR-1179 in CC

For the purpose of investigating the molecular mechanism of circOSBPL10 in regulating CC, we first detected its cellular sublocalization in SiHa and HeLa cells via subcellular fractionation. As illustrated in Fig. 2a, circOSBPL10 was majorly distributed in cytoplasm. Thus, we speculated that circOSBPL10 might affect CC via serving as a sponge of specific miRNA. After searching starBase (http://starbase.sysu.edu.cn/) with certain condition (CLIP Data: strict stringency ≥ 5, Degradome Data: low stringency ≥ 1), three miRNAs (miR-1179, miR-27a-3p and miR-27b-3p) were revealed to have binding potentials with circOSBPL10 (Fig. 2b). Then, we discovered a significant downregulation of miR-1179, whereas no apparent changes on the levels of miR-27a-3p and miR-27b-3p, in CC cell lines compared to normal H8 cells (Fig. 2c). Therefore, miR-1179 was chosen for further analysis. Subsequently, circOSBPL10 and miR-1179 were presented to be conspicuously concentrated in anti-Ago2 group (Fig. 2d). Afterwards, two binding sites between circOSBPL10 and miR-1179 were predicted via starBase (Fig. 2e). Moreover, we validated that miR-1179 bound with circOSBPL10 at site 1 (Fig. 2f). To further test whether circOSBPL10 promoted CC progression via its interaction with miR-1179, we mutated the sequence of circOSBPL10 recognized by miR-1179. As displayed in Fig. 2g, the expression of circOSBPL10 was observably elevated in C33A and CaSki cells after overexpressing circOSBPL10 and circOSBPL10 (mut). Interestingly, it seemed that cell proliferation and migration in C33A and CaSki cells could be fortified by overexpression of full-length circOSBPL10 but not by upregulation of circOSBPL10 with mutated miR-1179 binding sites (Fig. 2h, i), indicating that the function of circOSBPL10 in CC depended on its binding to miR-1179. In conclusion, circOSBPL10 interacts with miR-1179 to drive CC progression.

Fig. 2
figure 2

CircOSBPL10 sponged miR-1179 in CC. a Subcellular fractionation was applied to detect the cellular sublocalization of circOSBPL10 in SiHa and HeLa cells. b MiR-1179, miR-27a-3p and miR-27b-3p were predicted via starBase to have the binding capacity with circOSBPL10. c MiR-1179 expression in CC cell lines and H8 cells was detected via RT-qPCR. d RIP assay unveiled the significant enrichment of circOSBPL10 and miR-1179 in anti-Ago2 group. e Two binding sites between circOSBPL10 and miR-1179 were predicted via starBase. f Luciferase reporter assay validated the interaction between circOSBPL10 and miR-1179. g The expression of circOSBPL10 in SiHa and HeLa cells transfected with pcDNA3.1, pcDNA3.1/circOSBPL10 or pcDNA3.1/circOSBPL10 (mut) was analyzed through RT-qPCR. h, i Cell proliferation and migration in different groups was analyzed via CCK-8 and transwell assays, respectively. *P < 0.05, **P < 0.01

CircOSBPL10 upregulates UBE2Q1 level in CC by sequestering miR-1179

To further investigate the downstream mechanism of circOSBPL10 in CC, starBase was utilized. As predicted under certain circumstances (CLIP Data: strict stringency ≥ 3, Degradome Data: low stringency ≥ 1), four candidates (UBE2Q1, PRPF38B, CACUL1 and LUC7L3) were found as the targets of miR-1179 (Fig. 3a). The following RNA pull-down assay demonstrated the distinct enrichment of UBE2Q1 whereas limited harvest of other three mRNAs in Bio-miR-1179-WT group (Fig. 3b). More importantly, UBE2Q1 level was cut down in SiHa and HeLa cells by circOSBPL10 knockdown as well as by augmented miR-1179 expression (Fig. 3c, and Additional file 2). Besides, RT-qPCR obtained a conspicuous elevation on UBE2Q1 expression in CC cell lines relative to H8 cells, at both mRNA and protein levels (Fig. 3e and Additional file 2). Afterwards, the binding sites between UBE2Q1 and miR-1179 were obtained through starBase prediction (Fig. 3f). Importantly, we observed declined luciferase activity of UBE2Q1-WT owing to miR-1179 upregulation through luciferase reporter assay. However, no overt changes on that of UBE2Q1-Mut were noted between two groups (Fig. 3g). More importantly, a significant enrichment of circOSBPL10, miR-1179 and UBE2Q1 in Anti-Ago2 groups was observed via RIP analysis (Fig. 3h). Further, UBE2Q1 expression was visibly upregulated by circOSBPL10 overexpression, whereas enhanced expression of mutated circOSBPL10 had no such influence on UBE2Q1 expression in these two CC cells (Fig. 3i and Additional file 2). In sum, circOSBPL10 facilitates UBE2Q1 expression by sponging miR-1179 in CC.

Fig. 3
figure 3

CircOSBPL10 upregulated UBE2Q1 expression by competitively binding with miR-1179 in CC. a Venn diagram showed the overlaps of potential miR-1179 targets predicted by PITA, microT, PicTar, miRmap. b The binding capacity between miR-1179 and four mRNAs was verified through RNA pull-down assay. c, d The expression of UBE2Q1 in SiHa and HeLa cells transfected with different plasmids was detected via RT-qPCR and western blot. e UBE2Q1 expression in CC cell lines and H8 cells was detected via RT-qPCR. f The binding sites between UBE2Q1 and miR-1179 obtained from starBase were displayed. g The interaction between UBE2Q1 and miR-1179 was validated by luciferase reporter assay. h RIP analysis revealed that circOSBPL10, miR-1179 and UBE2Q1 co-existed in RISCs (RNA-induced silencing complexes). i RT-qPCR and western blot were utilized to detect the expression of UBE2Q1 in indicated cells. **P < 0.01

CircOSBPL10 accelerates CC progression via miR-1179/UBE2Q1 axis

Given that circOSBPL10 sponged miR-1179 to upregulate UBE2Q1 expression in CC, we intended to detect whether this mechanism contributed to CC progression. Herein, the efficacy of elevating UBE2Q1 expression or lowering miR-1179 level was analyzed at first by RT-qPCR and the outcome appeared to be satisfactory in SiHa cells (Fig. 4a). Afterwards, we testified that UBE2Q1 upregulation or miR-1179 inhibition could offset circOSBPL10 depletion-mediated suppressive effect on cell proliferation (Fig. 4b, c). In addition, cell apoptosis facilitated by circOSBPL10 knockdown was counteracted by UBE2Q1 overexpression or miR-1179 suppression (Fig. 4d, e). Moreover, the attenuated cell migration capability in circOSBPL10-depleted cells was restored by upregulation of UBE2Q1 or inhibition of miR-1179 (Fig. 4f, g). To conclude, circOSBPL10 promotes CC progression via miR-1179/UBE2Q1 axis.

Fig. 4
figure 4

CircOSBPL10 accelerated CC progression via miR-1179/UBE2Q1 axis. a The efficacy of UBE2Q1 overexpression and miR-1179 inhibition in SiHa cells was analyzed by RT-qPCR. b, c The proliferation ability of SiHa cells under diverse conditions were measured by CCK-8 and colony formation assays. d, e The apoptosis ability of transfected cells were measured via TUNEL assay and flow cytometry analysis. f, g Transwell and wound healing assays were carried out to analyze cell migration under different contexts. **P < 0.01

FOXA1 activates circOSBPL10 expression in CC

Since the downstream signaling responsible for the regulation of circOSBPL10 in CC had been studied, we subsequently focused on figuring out its possible upstream mechanism. Through utilizing UCSC (University of California, Santa Cruz: http://genome.ucsc.edu/), FOXA1 seemed to be a probable transcription factor of OSBPL10 (the host gene of circOSBPL10). Prior to testify the influence of FOXA1 on circOSBPL10 expression in CC, we first silenced or overexpressed FOXA1 with satisfactory efficacies in SiHa and HeLa cells (Fig. 5a). Of interest, we discovered that the expression of circOSBPL10 was distinctly declined by FOXA1 depletion whereas overtly augmented by FOXA1 upregulation (Fig. 5b). Later on, we employed JASPAR database (http://jaspar.genereg.net/) and obtained the DNA motif of FOXA1 (Fig. 5c). Seen from Fig. 5d, the sequences of OSBPL10 promoter were fragmented into five parts (P1-5). Interestingly, it was verified by ChIP assay that FOXA1 mainly bound to OSBPL10 promoter at P4 region (Fig. 5e). Moreover, when overexpressing FOXA1 in SiHa and HeLa cells, the luciferase activity of OSBPL10 promoter-WT was observably increased while there were no evident changes on that of OSBPL10 promoter-Mut (with mutated FOXA1 binding sites predicted in P4 region) (Fig. 5f). In a word, FOXA1 activates OSBPL10 transcription and thereby facilitates circOSBPL10 expression in CC.

Fig. 5
figure 5

FOXA1 activated circOSBPL10 expression in CC. a The efficacy of FOXA1 knockdown and overexpression in SiHa and HeLa cells was obtained through RT-qPCR analysis. b The expression of circOSBPL10 in different groups was detected by RT-qPCR. c The DNA motif of FOXA1 was obtained from JASPAR database. d Full-length of OSBPL10 promoter was fragmented into 5 pieces (P1-5), and JASPAR predicted a potential binding site for FOXA1 to OSBPL10 promoter at P4 region. e ChIP assay proved that the P4 region of OSBPL10 promoter was recognized by FOXA1 in both SiHa and HeLa cells. f The interaction between OSBPL10 promoter and FOXA1 at above predicted sites was testified by luciferase reporter assay. **P < 0.01

Discussion

Mounting evidence has manifested that circRNAs serve vital parts in CC initiation and progression. For instance, circRNA SMARCA5 regulates CC progression by sponging miR-620 [26]. Increased expression of circ_0067934 accelerates CC development by targeting miR-545/EIF3C axis [27]. CircRNA hsa_circ_0023404 plays an oncogenic part in CC [28]. CircOSBPL10, derived from back-splicing of OSBPL10 mRNA, is a circular RNA that has not been studied in CC but is worth exploring. In this research, circOSBPL10 was evidenced to have a circular structure and possess high levels in CC. Additionally, silenced circOSBPL10 exerted suppressive impacts on CC cell proliferation and migration.

Increasing researches have elucidated that circRNAs might serve as a molecular sponge of specific miRNAs, which has been suggested to be of significant value in CC [29, 30], so as to regulate the tumorigenesis and development of numerous cancers including CC [22, 23]. In current study, on the basis of bioinformatics prediction and molecular mechanism experiments, miR-1179 that was unveiled as an anti-tumor regulator in some cancers [31, 32], was screened out and validated to be implicated in circOSBPL10- regulated cellular processes in CC.

Previously, UBE2Q1 has been depicted as a critical participator in tumor progression [33, 34]. In present study, UBE2Q1 was verified capable of binding with miR-1179, and its expression was boosted by circOSBPL10 in CC through miR-1179 sequestration. More importantly, the follow-up rescue experiments indicated that UBE2Q1 upregulation or miR-1179 inhibition could rescue circOSBPL10 depletion-mediated suppressive effects on malignant phenotypes in CC.

Increasing researches have indicated that FOXA1 expresses at high levels in many cancers, such as lung cancer [35], glioma [36] and prostate cancer [37]. Besides, a previous study indicated that FOXA1 could directly bind with PLOD2 promoter and activate PLOD2 transcription in NSCLC [38]. Similarly, we revealed that FOXA1 activates OSBPL10 transcription and thereby facilitates circOSBPL10 expression in CC in this study.

Conclusion

In conclusion, FOXA1-induced circOSBPL10 promotes CC progression by targeting miR-1179/UBE2Q1 axis, providing novel insights into exploring more effective treatment of CC. However, the biggest regret is that no clinical data are included in this work, and the findings in our study need to be further testified by clinical samples in the future.

Availability of data and materials

Research data and material are not shared.

Abbreviations

circRNA:

Circular RNA

CC:

Cervical cancer

UBE2Q1:

Ubiquitin conjugating enzyme E2 Q1

FOXA1:

Forkhead box A1

ncRNAs:

Noncoding RNAs

FBS:

Fetal bovine serum

ActD:

Actinomycin D

gDNA:

Genomic DNA

References

  1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA. 2011;61(2):69–90.

    PubMed  Google Scholar 

  2. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA. 2015;65(2):87–108.

    PubMed  Google Scholar 

  3. Tewari KS, Sill MW, Long HJ 3rd, Penson RT, Huang H, Ramondetta LM, Landrum LM, Oaknin A, Reid TJ, Leitao MM, et al. Improved survival with bevacizumab in advanced cervical cancer. New Engl J Med. 2014;370(8):734–43.

    CAS  PubMed  Google Scholar 

  4. Chen W, Zheng R, Zeng H, Zhang S, He J. Annual report on status of cancer in China, 2011. Chinese J Cancer Res. 2015;27(1):2–12.

    Google Scholar 

  5. Testing for cervical cancer: new recommendations from the American Cancer Society, American Society for Colposcopy and Cervical Pathology, and American Society for Clinical Pathology. CA. 2012, 62(3):211–212.

  6. Burki TK. Cervical cancer: screening and risk with age. Lancet Oncol. 2014;15(3):e107.

    PubMed  Google Scholar 

  7. Noordhuis MG, Fehrmann RS, Wisman GB, Nijhuis ER, van Zanden JJ, Moerland PD, Ver Loren van Themaat E, Volders HH, Kok M, ten Hoor KA, et al. Involvement of the TGF-beta and beta-catenin pathways in pelvic lymph node metastasis in early-stage cervical cancer. Clin Cancer Res. 2011;17(6):1317–30.

    CAS  PubMed  Google Scholar 

  8. Shafabakhsh R, Reiter RJ, Mirzaei H, Teymoordash SN, Asemi Z. Melatonin: a new inhibitor agent for cervical cancer treatment. J Cell Physiol. 2019;234(12):21670–82.

    CAS  PubMed  Google Scholar 

  9. Ghasemi F, Shafiee M, Banikazemi Z, Pourhanifeh MH, Khanbabaei H, Shamshirian A, Amiri Moghadam S, ArefNezhad R, Sahebkar A, Avan A, et al. Curcumin inhibits NF-kB and Wnt/β-catenin pathways in cervical cancer cells. Pathol Res Pract. 2019;215(10):152556.

    PubMed  Google Scholar 

  10. Perez DS, Hoage TR, Pritchett JR, Ducharme-Smith AL, Halling ML, Ganapathiraju SC, Streng PS, Smith DI. Long, abundantly expressed non-coding transcripts are altered in cancer. Hum Mol Genet. 2008;17(5):642–55.

    CAS  PubMed  Google Scholar 

  11. Guttman M, Donaghey J, Carey BW, Garber M, Grenier JK, Munson G, Young G, Lucas AB, Ach R, Bruhn L, et al. lincRNAs act in the circuitry controlling pluripotency and differentiation. Nature. 2011;477(7364):295–300.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Khani P, Nasri F, Khani Chamani F, Saeidi F, Sadri Nahand J, Tabibkhooei A, Mirzaei H. Genetic and epigenetic contribution to astrocytic gliomas pathogenesis. J Neurochem. 2019;148(2):188–203.

    CAS  PubMed  Google Scholar 

  13. Mirzaei H. Stroke in women: risk factors and clinical biomarkers. J Cell Biochem. 2017;118(12):4191–202.

    CAS  PubMed  Google Scholar 

  14. Saeedi Borujeni MJ, Esfandiary E, Baradaran A, Valiani A, Ghanadian M, Codoñer-Franch P, Basirat R, Alonso-Iglesias E, Mirzaei H, Yazdani A. Molecular aspects of pancreatic β-cell dysfunction: oxidative stress, microRNA, and long noncoding RNA. J Cell Physiol. 2019;234(6):8411–25.

    CAS  PubMed  Google Scholar 

  15. Asma V, Zahra S, Ahmad M, Soheila M, Sima F, Hamid RM, Afshin N, Amir S, Hamed M. Long non-coding RNAs as epigenetic regulators in cancer. Curr Pharm Des. 2019;25(33):3563–77.

    Google Scholar 

  16. Mirzaei H, Yazdi F, Salehi R, Mirzaei H. SiRNA and epigenetic aberrations in ovarian cancer. J Cancer Res Ther. 2016;12(2):498–508.

    CAS  PubMed  Google Scholar 

  17. Qu S, Zhong Y, Shang R, Zhang X, Song W, Kjems J, Li H. The emerging landscape of circular RNA in life processes. RNA Biol. 2017;14(8):992–9.

    PubMed  Google Scholar 

  18. Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, Sun W, Dou K, Li H. Circular RNA: a new star of noncoding RNAs. Cancer Lett. 2015;365(2):141–8.

    CAS  PubMed  Google Scholar 

  19. Shabaninejad Z, Vafadar A, Movahedpour A, Ghasemi Y, Namdar A, Fathizadeh H, Pourhanifeh MH, Savardashtaki A, Mirzaei H. Circular RNAs in cancer: new insights into functions and implications in ovarian cancer. J Ovarian Res. 2019;12(1):84.

    PubMed  PubMed Central  Google Scholar 

  20. Naeli P, Pourhanifeh MH, Karimzadeh MR, Shabaninejad Z, Movahedpour A, Tarrahimofrad H, Mirzaei HR, Bafrani HH, Savardashtaki A, Mirzaei H, et al. Circular RNAs and gastrointestinal cancers: epigenetic regulators with a prognostic and therapeutic role. Crit Rev Oncol Hematol. 2020;145:102854.

    PubMed  Google Scholar 

  21. Wang L, Peng X, Lu X, Wei Q, Chen M, Liu L. Inhibition of hsa_circ_0001313 (circCCDC66) induction enhances the radio-sensitivity of colon cancer cells via tumor suppressor miR-338-3p: effects of cicr_0001313 on colon cancer radio-sensitivity. Pathol Res Pract. 2019;215(4):689–96.

    CAS  PubMed  Google Scholar 

  22. Cai H, Zhang P, Xu M, Yan L, Liu N, Wu X. Circular RNA hsa_circ_0000263 participates in cervical cancer development by regulating target gene of miR-150-5p. J Cell Physiol. 2019;234(7):11391–400.

    CAS  PubMed  Google Scholar 

  23. Zhang HD, Jiang LH, Hou JC, Zhong SL, Zhou SY, Zhu LP, Li J, Wang DD, Sun DW, Ji ZL, et al. Circular RNA hsa_circ_0052112 promotes cell migration and invasion by acting as sponge for miR-125a-5p in breast cancer. Biomed Pharm. 2018;107:1342–53.

    CAS  Google Scholar 

  24. Yang L, Wang J, Fan Y, Yu K, Jiao B, Su X. Hsa_circ_0046264 up-regulated BRCA2 to suppress lung cancer through targeting hsa-miR-1245. Respir Res. 2018;19(1):115.

    PubMed  PubMed Central  Google Scholar 

  25. Huang W, Lu Y, Wang F, Huang X, Yu Z. Downregulation of circular RNA hsa_circ_0000144 inhibits bladder cancer progression via stimulating miR-217 and suppressing RUNX2 expression. Gene. 2018;678:337–42.

    CAS  PubMed  Google Scholar 

  26. Tian JDC, Liang L. Involvement of circular RNA SMARCA5/microRNA-620 axis in the regulation of cervical cancer cell proliferation, invasion and migration. Eur Rev Med Pharm Sci. 2018;22(24):8589–98.

    Google Scholar 

  27. Hu C, Wang Y, Li A, Zhang J, Xue F, Zhu L. Overexpressed circ_0067934 acts as an oncogene to facilitate cervical cancer progression via the miR-545/EIF3C axis. J Cell Physiol. 2019;234(6):9225–32.

    CAS  PubMed  Google Scholar 

  28. Zhang J, Zhao X, Zhang J, Zheng X, Li F. Circular RNA hsa_circ_0023404 exerts an oncogenic role in cervical cancer through regulating miR-136/TFCP2/YAP pathway. Biochem Biophys Res Commun. 2018;501(2):428–33.

    CAS  PubMed  Google Scholar 

  29. Nahand JS, Taghizadeh-boroujeni S, Karimzadeh M, Borran S, Pourhanifeh MH, Moghoofei M, Bokharaei-Salim F, Karampoor S, Jafari A, Asemi Z, et al. microRNAs: new prognostic, diagnostic, and therapeutic biomarkers in cervical cancer. J Cell Physiol. 2019;234(10):17064–99.

    CAS  PubMed  Google Scholar 

  30. Sadri Nahand J, Moghoofei M, Salmaninejad A, Bahmanpour Z, Karimzadeh M, Nasiri M, Mirzaei HR, Pourhanifeh MH, Bokharaei-Salim F, Mirzaei H, et al. Pathogenic role of exosomes and microRNAs in HPV-mediated inflammation and cervical cancer: a review. Int J Cancer. 2020;146(2):305–20.

    CAS  PubMed  Google Scholar 

  31. Li Y, Qin C. MiR-1179 inhibits the proliferation of gastric cancer cells by targeting HMGB1. Hum Cell. 2019;32(3):352–9.

    CAS  PubMed  Google Scholar 

  32. An J, Shi H, Zhang N, Song S. Elevation of circular RNA circ_0003645 forecasts unfavorable prognosis and facilitates cell progression via miR-1179/TMEM14A pathway in non-small cell lung cancer. Biochem Biophys Res Commun. 2019;511(4):921–5.

    CAS  PubMed  Google Scholar 

  33. Zhang B, Deng C, Wang L, Zhou F, Zhang S, Kang W, Zhan P, Chen J, Shen S, Guo H, et al. Upregulation of UBE2Q1 via gene copy number gain in hepatocellular carcinoma promotes cancer progression through beta-catenin-EGFR-PI3K-Akt-mTOR signaling pathway. Mol Carcinog. 2018;57(2):201–15.

    CAS  PubMed  Google Scholar 

  34. Chang R, Wei L, Lu Y, Cui X, Lu C, Liu L, Jiang D, Xiong Y, Wang G, Wan C, et al. Upregulated expression of ubiquitin-conjugating enzyme E2Q1 (UBE2Q1) is associated with enhanced cell proliferation and poor prognosis in human hapatocellular carcinoma. J Mol Histol. 2015;46(1):45–56.

    CAS  PubMed  Google Scholar 

  35. Li X, Chen S, Sun T, Xu Y, Chen Y, Liu Y, Xiang R, Li N. The transcriptional regulation of SOX2 on FOXA1 gene and its application in diagnosis of human breast and lung cancers. Clin Lab. 2014;60(6):909–18.

    CAS  PubMed  Google Scholar 

  36. Wang L, Qin H, Li L, Feng F, Ji P, Zhang J, Li G, Zhao Z, Gao G. Forkhead-box A1 transcription factor is a novel adverse prognosis marker in human glioma. J Clin Neurosci. 2013;20(5):654–8.

    PubMed  Google Scholar 

  37. Robbins CM, Tembe WA, Baker A, Sinari S, Moses TY, Beckstrom-Sternberg S, Beckstrom-Sternberg J, Barrett M, Long J, Chinnaiyan A, et al. Copy number and targeted mutational analysis reveals novel somatic events in metastatic prostate tumors. Genome Res. 2011;21(1):47–55.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Du H, Chen Y, Hou X, Huang Y, Wei X, Yu X, Feng S, Wu Y, Zhan M, Shi X, et al. PLOD2 regulated by transcription factor FOXA1 promotes metastasis in NSCLC. Cell Death Dis. 2017;8(10):e3143.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are very grateful to all individuals and groups involved in this study.

Funding

This study was supported by The National Natural Science Foundation of China (81872460), The National Cancer Center Climbing Foundation of China(NCC201808B016), Department of Science and Technology of Sichuan Province(20YYJC3815), Department of Science and Technology of Sichuan Province(20GJHZ0088).

Author information

Authors and Affiliations

Authors

Contributions

SY, YJ, XR and DF designed this study; SY, LZ, FZ, YJ and XR devoted to experiment, data curation and interpretation; SY, YJ, XR, DF, DH, SH and LZ were in charge of investigation, experiment record, figures and article writing; SL  edited the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Liaoyun Zhang or Shun Lu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

List of the sequences of primers used in RT-qPCR.

Additional file 2: Figure S1.

The size of marker for western blot gels.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, S., Jiang, Y., Ren, X. et al. FOXA1-induced circOSBPL10 potentiates cervical cancer cell proliferation and migration through miR-1179/UBE2Q1 axis. Cancer Cell Int 20, 389 (2020). https://doi.org/10.1186/s12935-020-01360-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-020-01360-2

Keywords