What this is
- This research investigates the role of circular RNA () in epithelial ovarian cancer (EOC).
- CircMUC16 was identified as having increased expression in EOC tissues compared to healthy tissues.
- The study explores how circMUC16 promotes and contributes to cancer progression through interactions with specific proteins and miRNAs.
Essence
- CircMUC16 is overexpressed in epithelial ovarian cancer and promotes , enhancing cancer cell invasion and proliferation through interactions with miR-199a-5p and ATG13.
Key takeaways
- CircMUC16 expression is elevated in EOC tissues compared to healthy ovarian tissues, correlating with tumor progression.
- Silencing circMUC16 reduces in SKOV3 cells, while its overexpression increases in A2780 cells, indicating its role in promoting cancer cell behavior.
- CircMUC16 interacts with miR-199a-5p, regulating Beclin1 and RUNX1, thereby forming a feedback loop that enhances and cancer progression.
Caveats
- The study relies on a limited number of EOC samples, which may affect the generalizability of the findings.
- Further research is needed to fully elucidate the mechanisms by which circMUC16 influences and cancer behavior.
Definitions
- circRNA: A class of noncoding RNA that forms a covalently closed loop, involved in regulating gene expression.
- autophagy: A cellular process that degrades and recycles cellular components, often linked to cancer progression.
Simplified
Background
Ovarian carcinoma is one of the most common gynecological malignant tumors with a high mortality rate [1]. Despite rapid progress in the chemotherapy and radiotherapy of ovarian carcinoma, patients with recurrent ovarian cancer are essentially incurable [2]. Unfortunately, the molecular etiology of ovarian cancer remains elusive [3].
Circular RNA (circRNA) is a class of noncoding RNAs discovered in recent years [4]. CircRNAs have various biological functions, which includes regulation of cell proliferation, invasion, apoptosis, etc. [5]. CircRNA ITCH inhibited the malignant behavior of esophageal squamous cell carcinoma through suppression of Wnt/β-catenin pathway [6]. CircHIPK3 regulates the growth of cells by binding directly to miR-124 and inhibits its activity [7]. CircRNA_001059 and circRNA_000167 are involved in the development of radiation resistance [8]. CircRNA ITCH inhibits proliferation and promotes apoptosis of human epithelial ovarian cancer (EOC) cells by sponging miR-10a [9]. The expression of circLARP4 decreases in ovarian cancer, and circLARP4 serves as a potential biomarker of ovarian cancer prognosis [10]. CircITCH is related to tumor size, FIGO stage and overall survival. It also suppresses cell proliferation and promotes apoptosis in EOC [11].
However, there is no complete and systematic study on the expression and role of circRNAs in ovarian cancer. In this study, we performed genome-wide circRNA analysis of EOC tissues and healthy ovarian tissues via next-generation sequencing (NGS) technology. The aim was to investigate the role and demonstrate the potential mechanism of circRNAs in EOC.
Materials and methods
Ovarian tissues
The ovarian and serum samples were obtained from the Department of Obstetrics and Gynecology, Second Affiliated Hospital, Chongqing Medical University. The Ethics Committee of Chongqing Medical University approved this study. Three pathologists confirmed the pathological diagnosis of ovarian cancer.
RNA sequencing analysis
First, we characterized circRNA transcripts via sequencing analysis of ribosomal RNA and linear RNA. Total RNA was extracted from three EOC tissues and four healthy ovarian tissues. Each sample was sequenced on Illumina HiSeq and yielded an average of 57.74 million reads, which were mapped to the human reference genome (GRCh38/hg38) by TopHat2. Qualified reads were sent to CIRCexplorer v2.2.3 pipeline to identify and quantify circRNAs with default parameters. A circRNA was called with the support of a minimum of two unique back-spliced reads. The whole process of library construction, sequencing and data analysis was performed at Shanghai Lifegenes Technology Co., Ltd.
Cell culture, transfection procedure, and reagents
SKOV3 human ovarian cancer cell line derived from the ascites from a 64 year old caucasian female with an ovarian tumor. The ES-2 cell line was established from a surgical tumor specimen taken from a 47 year old black woman. The A2780 cell line was established from tumor tissue from an untreated patient. CAOV-3 human ovarian cancer cell line derived from a 54 years old caucasian female with ovarian adenocarcinoma. Human ovarian cancer cells were cultured in RPMI 1640 medium (Sigma-Aldrich, R8758), with 10% fetal bovine serum and streptomycin. The cells were incubated under 5% CO2 at 37 °C [12]. Lentiviral shRNA interference vectors targeting circMUC16 (named LV2–1 and LV2–2) and the circMUC16-lentiviral expression vector (named LV6) were purchased from Genepharma. MIR199A mimics and inhibitors were synthesized by Genepharma. We constructed pCMV5-Beclin1 for overexpression of Beclin1. The siRNA was synthesized by Genepharma (Shanghai, China). The sequences are as follows: CircMUC16–1: 5′-TGACCTTGCTCAGGCCCGA-3′; CircMUC16–1: 5′-ACCTTGCTCAGGCCCGAGA-3′; RUNX1: 5′-GGAUUUCUGUUGUGUUUAAAU-3′; BECLIN1: 5′-GGUGUUUGAUACUGUUUGAGA-3′; ATG5: 5′-CAAUCCCAUCCAGAGUUGCUUGUGA-3′; and NC (negative control) siRNA: 5′-UUCUUCGAAGGUGUCACGUTT-3′.
Tandem mRFP-GFP-LC3 fluorescence
We used a tandem mRFP-GFP-tagged LC3 to monitor autophagy flux based on our previous study [12]. SKOV3 cells were infected with mRFP-GFP-LC3 and LV2–1, LV2–2 or LV2-NC. A2780 cells were infected with mRFP-GFP-LC3 and LV6 or LV6-NC. mRFP-GFP-LC3 distribution in cells was analyzed by confocal microscope. Image pro plus 6.0 software was used to quantify the LC3 dots. All experiments were repeated three times [12].
RIP assay
Biotin-labelled circMUC16 probe was purchased from Genepharma. The specific operation protocol was in accordance with a previous study [13]. Finally, the expression of miRNA was detected using qPCR.
CircRNA pulldown
The specific steps of circRNA pulldown are described in our previous study [14]. Briefly, the specific proteins were identified using SDS-PAGE gel electrophoresis, silver stain and MALDI-TOF-MS [14].
qPCR
Total RNA was isolated using a high-purity total RNA rapid extraction kit (Bioteke Corporation, RP1201) following the manufacturer’s instructions. RNA was reverse-transcribed into cDNA using the iSCRIPT cDNA synthesis kit (Bio-Rad Laboratories, 4,106,228). The forward primer for circMUC16 sequence was 5′-CTGCTCAGGCCTGTGTTC-3′, and the reverse primer sequence was 5′-GGGGCCCAGCTCTTCA-3′. The real-time PCR was performed using the All-in-One qPCR mix kit. Each sample was analyzed in triplicate. The primers were purchased from Genepharma. The 2-ΔΔCT method was used to quantify the gene expression [14].
Chromatin immunoprecipitation (ChIP) assay
We used ChIP assay kit (Cell Signaling Technology, 9003) to conduct a ChIP assay according to the manufacturer’s instructions. Briefly, cross-linked chromatin was sonicated into 200–1000-bp fragments. Then, the chromatin was immunoprecipitated using anti-RUNX1 antibody (ab23980, Abcam). qPCR was performed as described earlier [12].
Dual-luciferase reporter gene assay
Luciferase reporter gene assay was performed using the dual-luciferase reporter assay system (Promega Corporation, E1910) according to the manufacturer’s instructions. The wild-type or mutant reporter of circMUC16, Beclin1 and RUNX1 was constructed. Then, the reporter was co-transfected into SKOV3 cells in 24-well plates with 100 nM miR-199a or 100 nM miR-NC and Renilla plasmid using Endofectin-Plus (Gene-Copeia, Z01010A). All experiments were repeated three times [12].
Western blotting
The expressions of Beclin1, RUNX1, ATG14, and GAPDH proteins were detected by western blot. The primary antibodies used included rabbit monoclonal anti-Beclin1 (Abcam, ab210498); rabbit polyclonal RUNX1 (Proteintech, 25,315–1-AP); rabbit polyclonal ATG14 (Abcam, ab227849), and rabbit monoclonal anti-GAPDH (Abcam, ab181602). The band density was analyzed using a gel imaging system and compared with the internal control.
EdU assay
The cell proliferation was detected using Cell-Light Apollo567 in vitro kit (RiboBio, Guangzhou, China). Ovarian cancer cells were cultured in 1640 medium. EdU was applied at 20 μM. The cells were fixed with 4% paraformaldehyde and stained with Apollo 567 and Hoechst 33342 [12].
Matrigel invasion assays
Matrigel invasion assays were used to determine the cellular invasion ability as described in our previous study [12]. A total of 5 × 104 cells were seeded into the top chamber of the transwell filter and incubated. After 48 h, the cells on the lower side of the chamber were fixed with 4% paraformaldehyde, stained with 0.5% crystal violet (Beyotime Institute of Biotechnology, C0121), and counted using a microscope [12].
Fluorescence in situ hybridization (FISH) assay
Cy3-labelled probe against circMUC16 (5’CY3-TCTTTCTCGGGCCTGAGCAAGGTCAGT- 3’CY3) and FITC-labeled probe against miR-199a were obtained from Genepharma. The FISH assay was performed following the manufacturer’s guide and instructions (F03401↗, Genepharma). Cell nuclei were stained with Hoechst 33342. Finally, images were obtained on a confocal microscope [15, 16]. All experiments were repeated three times.
Mouse xenograft model
The xenograft model was established as described in our previous study [12]. All animal experimental procedures were conducted with the approval of the Committee on the Use and Care of Animals (Chongqing Medical University, Chongqing, China), based on the institution’s guidelines. SKOV3 cells were infected with LV2-NC and LV2–1, and intraperitoneally injected into six-week-old BALB/c nude mice (5 × 106 cells). After five weeks, the animals were examined. Then, the number of mice with ascites were counted and weighed [12].
Statistical analysis
All statistical analyses were performed with SPSS version 17.0 software (Chicago, IL). Statistical analysis was conducted by Student’s t-test or analysis of variance (ANOVA). The chi-square test was used to compare the associations between circMUC16 expression and the clinicopathological variables of ovarian cancer samples. Data were presented as mean ± standard deviation. A two-sided p < 0.05 was considered significant.
Results
Identification of differentially expressed circRNAs in epithelial ovarian cancer tissues

Hierarchical clustering of circRNA differential expression profiles between three EOC samples and four normal ovarian samples. The heat maps are generated from differentially expressed circRNAs (top 50). Red and green colors indicate relative expression above or below the average.Volcano map of deferentially expressed circRNAs (top 50).Functional analysis of differentially expressed circRNAs (top 50).The top 50 expressed circRNAs were analyzed by KEGG pathway. (D) CircMUC16 expression was detected in 30 normal ovarian tissues and 70 EOC samples.The expression of circMUC16 in 30 paired serum samples with normal ovarian tissues and 70 serum samples with epithelial ovarian cancer tissues. Error bars represent the standard error. The symbols * and ** show< 0.05 and 0.01, respectively a b c d e p
| No. of patients | CircMUC16 expression | valueP | ||
|---|---|---|---|---|
| (= 100)n | Low no.(%) | High no.(%) | ||
| Characteristics | ||||
| Age (years) | > 0.05 | |||
| < 50 | 48 | 23(47.92%) | 25(52.08%) | |
| ≥ 50 | 52 | 24(46.15%) | 28(53.85%) | |
| Normal ovarian | 30 | 28(93.33%) | 2(6.67%) | < 0.05 |
| Cancer tissues | 70 | 19(27.14%) | 51(72.86%) | |
| FIGO stage | ||||
| I/II | 40 | 15(37.50%) | 25(62.50%) | < 0.05 |
| III/IV | 30 | 4(13.33%) | 26(86.67%) | |
| Grade | ||||
| 1 | 16 | 12(75.00%) | 4(25.00%) | |
| 2 | 21 | 5(23.81%) | 16(76.19%) | |
| 3 | 33 | 2(6.06%) | 31(93.94%) | |
| Grade 2–3 versus 1 | < 0.05 | |||
CircMUC16 mediated autophagy flux of epithelial ovarian cancer

CircMUC16-mediated autophagy flux of epithelial ovarian cancer.CircMUC16 expression was determined in four cell lines by qPCR.CircMUC16 expression was detected in SKOV3 cells.The expression of circMUC16 was detected in A2780 cells.Western blot was used to detect the expression of LC3-II after silencing or ectopic expression of circMUC16.mRFP-GFP-LC3 distributions in SKOV3 and A2780 cells were analyzed by confocal microscopy after silencing or ectopic expression of circMUC16. The symbols * and ** show< 0.05 and 0.01, respectively. Scale bar: 5 μm a b c d e p

CircMUC16-mediated autophagy promoted cellar proliferation and invasion.,SKOV3 cells were transfected with LV2-NC, LV2–1, LV2–1 + PCMV5-ATG5, or LV2–1 + PCMV5-Beclin1. A2780 cells were transfected with LV6-NC, LV6, LV6 + ATG5 siRNA, or LV6 + Beclin1 siRNA. Then, western blot was used to detected the expression of LC3.-SKOV3 cells were transfected with LV2-NC, LV2–1, LV2–1 + PCMV5-ATG5, or LV2–1 + PCMV5-Beclin1. A2780 cells were transfected with LV6-NC, LV6, LV6 + ATG5 siRNA, or LV6 + Beclin1 siRNA. Then, the cellular proliferation and invasion was detected. Error bars represent the standard error. The symbols * and ** indicate< 0.05 and 0.01, respectively. Scale bar: 100 μm a b c f p
CircMUC16-mediated autophagy accelerated proliferation and invasion of epithelial ovarian cancer cells
Since circMUC16 was found to regulate autophagy, we investigated whether circMUC16-mediated autophagy was involved in cellular proliferation and migration. The expression of LC3-II was decreased in SKOV3 cells after knockdown of circMUC16. However, ectopic expression of ATG5 or Beclin1 partly reversed this effect (Fig. 3a). The LC3-II expression was elevated in A2780 cells infected with LV6-circMUC16. Consequently, when ATG5 or Beclin1 was silenced, this effect was partly reversed (Fig. 3b). The cellular proliferation and invasion of A2780 cells infected with LV6 was increased. Also, when autophagy was inhibited via silencing ATG5 or Beclin1, this effect was partly retarded (Fig. 3c-f). Hence, circMUC16-mediated autophagy promoted proliferation and invasion of ovarian cancer.
CircMUC16-miR-199a-5p-Beclin1 axis regulated autophagy

CircMUC16-miR-199a-5p-Beclin1 axis regulated autophagy.The expression of miRNAs of SKOV3 cells was determined by qPCR.The expression of miRNAs of A2780 cells was determined by qPCR.The location of miR-199a and circMUC16 in SKOV3 cells was detected using FISH assay.RNA pulled down assay was performed. Then, the expression of circMUC16 was detected by qPCR.RNA pulled down assay was performed. Then, the expression of miR-199a was detected by qPCR.The expression of Beclin1 and ATG14 of SKOV3 cells was determined using western blot.SKOV3 cells were transfected with LV2-NC, LV2–1, LV2–2, LV2–1 + miR-199a inhibitors, miR-199a mimics, miR-199a mimic inhibitors. Then, the expression of Beclin1 and ATG14 were detected using western blot.SKOV3 cells were co-transfected with miR-199a mimic or control RNA (NC) using luciferase reporter plasmids containing either wild-type (pMIR-CircMUC16) or mutant (pMIR-CircMUC16m). Luciferase expression was measured. Error bars represent the standard error. The symbols * and ** indicate< 0.05 and 0.01, respectively a b c d e f f g p
CircMUC16 regulated RUNX1 via sponging miR-199a-5p

CircMUC16 regulated RUNX1 via sponging miR-199a-5p.The expression of RUNX1 was determined by western blot.SKOV3 cells were co-transfected with miR-199a mimic or control RNA (NC) using luciferase reporter plasmids containing either wild-type (pMIR-RUNX1 3’UTR) or mutant (pMIR-RUNX1 3’UTRm). Luciferase expression was measured.A2780 cells were transfected with LV6, LV6-NC, LV6 + RUNX1 siRNA, or LV6 + Beclin1 siRNA. After 48 h, the expression of LC3-II was determined by western blot.,A2780 cells were transfected with LV6, LV6-NC, LV6 + RUNX1 siRNA, or LV6 + Beclin1 siRNA. Then, mRFP-GFP-LC3 distributions were determined by confocal microscopy. Error bars represent the standard error. The symbols * and ** indicate< 0.05 and 0.01, respectively. Scale bar: 5 μm a b c d e p
RUNX1 promoted the expression of circMUC16

RUNX1 promoted the expression of circMUC16.The putative binding sites between circMUC16 and RUNX1.SKOV3 cells co-transfected with MUC16 reporter gene plasmid and NC (negative control), RUNX1 siRNA, PCMV5-RUNX1 or PCMV5-NC (empty vector). After 48 h, luciferase activity was measured.SKOV3 cells transfected with NC (negative control), RUNX1 siRNA, PCMV5-RUNX1 or PCMV5-NC (empty vector). After 48 h, the linear mRNAs of MUC16 and circMUC16 were detected by qPCR.Overall survival was analyzed in ovarian cancer patients based on TCGA data. All experiments were repeated three times. Error bars represent the standard error. The symbols * and ** indicate p < 0.05 and 0.01, respectively a b c d
CircMUC16 directly binds to ATG13 protein

CircMUC16 directly bound to ATG13 protein.The circMUC16 RNA pulldown was performed. Then, silver stains was obtained.The circMUC16 specific binding protein gels was identified by RNA pulldown.The expression of ATG13 protein was determined by western blot.The potential interaction sites between circMUC16 and ATG13 protein were predicted using catRAPID.Segments of ATG13 DNA were cloned into a PCMV5 plasmid. SKOV3 cells were transfected with PCMV5 vector, PCMV5-ATG13, PCMV5-△51–102, PCMV5-△200–251 or PCMV5-△475–526. After 48 h, the expression of circMUC16 was determined. Error bars represent the standard error. The symbols * and ** indicate p < 0.05 and 0.01, respectively a b c d e
CircMUC16 mediated metastasis of SKOV3 cells in vivo

CircMUC16 mediated metastasis of EOC cells in vivo.Silencing circMUC16 suppressed pelvic peritoneal invasion. The arrow showed the tumor in the abdominal cavity.Mean tumor weight was calculated.Immunohistochemical analysis of Beclin1, RUNX1 and ATG13 proteins.A graphic abstract depicting the circMUC16 pathway (CircMUC16, RUNX1, Beclin1, ATG13 and miR-199a). All experiments were repeated three times. Error bars represent the standard error. The symbols * and ** indicate p < 0.05 and 0.01, respectively. Scale bar: 100 μm a b c d
Discussion
This study showed that circMUC16 was increased expression in EOC samples than healthy ovarian tissues. The expression of circMUC16 was related to tumor grade and stage. Furthermore, circMUC16 regulated Beclin1 or RUNX1 via sponging miR-199a-5p. In turn, RUNX1 promoted the transcription of circMUC16. Therefore, circMUC16 was a potential therapeutic target and diagnostic marker.
CircRNAs have been regarded as markers of many cancer types. In progressive laryngeal cancer, hsa_circRNA_104912 had increased expression and was a novel biomarker [18]. The expression of hsa_circ_0000190 was related to tumor progression, TNM stage and CA19–9 levels [19]. In this research, we found that circMUC16 was increased expression in EOC samples and sera from patients with EOC. The expression of circMUC16 correlated with tumor grade and stage. Therefore, circMUC16 was a potential biomarker for EOC.
CircRNAs typically regulate the target genes by sponging miRNAs. cir-ITCH enhanced ITCH expression via sponging miR-7 and miR-214, and then suppressed Wnt/β-catenin signaling pathway [20]. CircRNA-MYLK directly bound to and sponged miR-29a, and then relieved the suppression of target VEGFA [21]. This study showed that circMUC16 regulated Beclin1 and RUNX1 by sponging miR-199a-5p. miR-199a suppressed autophagy by GSK3β/mTOR complex signaling [22]. Moreover, miR-199a directly targeted Beclin1 to inhibited autophagy and reversed drug resistance induced by cisplatin of osteosarcoma cells [23]. Our results also demonstrated that miR-199a targetly regulated Beclin1, which was consistent with the previous study. We also verified that circMUC16-mediated autophagy accelerated the cellular proliferation and invasion. Therefore, circMUC16 acted as a cancer-promoting gene.
Silencing RUNX1 inhibited cellar proliferation and promoted of apoptosis [24]. miR-141 can retard cell proliferation and migration. miR-141 also promoted apoptosis of prostatic cancer. Over-expression of RUNX1 reversed this phenotype. RUNX1 was the direct target of miR-141. In this research, we observed that circMUC16 elevated the expression of RUNX1 via sponging miR-199a-5p. Knockdown of RUNX1 partly retarded the autophagy induced by circMUC16. We also demonstrated that RUNX1 facilitated the expression of circMUC16 by promoting transcription. Therefore, circMUC16-miR-199a-5p-RUNX1 feedback loop boosted autophagy of EOC.
In summary, to the best of our knowledge, this is the first report on circMUC16 being abnormally expressed in EOC. CircMUC16 was correlated with the tumor stage and grade. CircMUC16-mediated autophagy accelerated the malignant behavior of EOC. Hence, circMUC16 was a potential target for the diagnosis and treatment of ovarian cancer.