CircCDKN2B-AS1 interacts with IMP3 to stabilize Hexokinase 2 mRNA and facilitate cervical squamous cellcarcinoma aerobic glycolysis progression

doi:10.21203/rs.3.rs-49280/v2

Download PDF

Research

CircCDKN2B-AS1 interacts with IMP3 to stabilize Hexokinase 2 mRNA and facilitate cervical squamous cellcarcinoma aerobic glycolysis progression

https://doi.org/10.21203/rs.3.rs-49280/v2

This work is licensed under a CC BY 4.0 License

Journal Publication

published 11 Dec, 2020

Read the published version in Journal of Experimental & Clinical Cancer Research →

You are reading this latest preprint version

Background: Circular RNAs (circRNAs) have been reported to play key roles in the development of various cancers. However, the biological functions and clinical significance of most circRNAs are still elusive. The purpose of this study was to explore the function and mechanism of a certain circRNA named circCDKN2B-AS1 in cervical cancer development and its potential value in the clinic.

Methods: qRT-PCR was used to verify the expression level of circCDKN2B-AS1. CCK-8, Transwell, and flow cytometry (FCM) assays were performed to detect cellular proliferation, migration, and apoptosis, respectively. A Seahorse XFe96 Analyzer was used to measure glycolysis metabolism level. RNA pull-down, RNA immunoprecipitation (RIP), actinomycin-D addition assays and Western blotting were used to screen and elucidate the potential mechanisms involved. BALB/c nude mice and zebrafish embryos (AB, WT) were used as animal models to investigate tumorigenesis capability. ¹⁸FDG-microPET/CT imaging and lactic acid (LA) and pyruvic acid (PA) content detection assays were used to detect the level of glucose metabolism in subcutaneous tumors from nude mice.

Results: CircCDKN2B-AS1, a circular isoform of the long noncoding RNA (lncRNA) CDKN2B-AS1, was upregulated in cervical cancer and precancerous tissues. We found that circCDKN2B-AS1 associated with the IMP3 protein depending on a specific binding site and regulated the stability of Hexokinase 2 (HK2) mRNA, the rate-limiting enzyme of the aerobic glycolysis pathway. The expression level of circCDKN2B-AS1 fated the binding of IMP3 to the 3’ untranslated region (UTR) of HK2 mRNA, consequently affecting the malignant cell phenotype and aerobic glycolysis in cervical cancer in vitro and in vivo. Mutant circCDKN2B-AS1, lacking the IMP3 binding site, did not have such effects. Utilization of an inhibitory peptide to block the interaction between circCDKN2B-AS1 and the IMP3 protein impeded the binding of IMP3 to the 3’UTR of HK2 mRNA and suppressed aerobic glycolysis in cervical cancer cells.

Conclusions: Our findings demonstrate that circCDKN2B-AS1 facilitates aerobic glycolysis by sponging the IMP3 protein to stabilize HK2 mRNA, consequently promoting the malignant phenotype in cervical cancer, which may provide a potential approach for cervical cancer therapeutics.

Oncology

CircRNA

CircCDKN2B-AS1

Cervical cancer

Aerobic glycolysis

IMP3

HK2

Cervical cancer is the fourth most common cancer among women globally, leading to an estimated 530,000 new cases and 270,000 deaths each year [1, 2]. Screening and human papillomavirus (HPV) vaccination have remarkably reduced the incidence of cervical cancer in developed countries, but the application of both methods is restricted, and advanced cervical cancer is still very common in low-income countries [3], where approximately 90% of the world’s cervical cancer-related deaths occur [4–6]. The curative effect is not ideal, and the 5-year survival rate is only 16.7% in advanced cervical cancer patients [7]. New therapeutic methods are expected for cervical cancer.

Circular RNAs (circRNAs), formed by a covalently closed loop, are a class of noncoding RNAs [8, 9]. CircRNAs resist RNaseR digestion and are more stable than their linear isoforms [9–11]. Recently, circRNAs have been increasingly recognized as being implicated in gene regulatory networks and influencing diseases [12, 13]. The effect of circRNAs on the development of various cancers, such as breast cancer, oral cancer and hepatocellular carcinoma, has been reported [14–16]. Most of the discoveries have revealed that circRNAs can function by sponging microRNAs and regulating the expression levels of target genes [17-18]. Recent studies have showed that circRNAs can bind to RNA binding proteins (RBPs) [19–21]; however, the biological function and potential value of the circRNA-RBP complex is still elusive.

IMP3, also known as IGF2BP3, is an RBP. A previous study revealed that IMP3 is highly expressed in various cancers and can regulate the expression of various genes at the posttranscriptional level [22–24]. HK2 is one of the rate-limiting enzymes in the glycolysis pathway and has been recognized to regulate the malignant phenotype of cancer cells (i.e., cellular apoptosis and migration capabilities) [25-27]. Aerobic glycolysis is the energy metabolism characteristic of tumor cells and fates the malignant phenotype of tumor cells, including cervical cancer cells [28-31]. Subcutaneous tumor formation in BALB/c nude mice combined with microPET/CT can be used to effectively explore tumorigenesis and the level of glucose metabolism in tumors [32]. Zebrafish embryo injection is a novel model for cancer research that has the advantages of a short experimental period and a simple operation procedure [33].

In our previous circRNA sequencing analysis in normal and cervical cancer tissues, we identified a circular isoform of the long noncoding RNA (lncRNA) CDKN2B-AS1, named circCDKN2B-AS1. We verified that circCDKN2B-AS1 was upregulated in HPV16-positive cervical cancer and precancerous tissues compared with normal tissues. We further identified that circCDKN2B-AS1 facilitates aerobic glycolysis by interacting with the IMP3 protein to stabilize HK2 mRNA, consequently promoting the malignant phenotype in cervical cancer cells. To better explore the potential therapeutic role of circCDKN2B-AS1, we selected the above two animal models and verified the capability of circCDKN2B-AS1 to regulate HK2 expression and glycolysis by binding IMP3 protein in vivo. Our study may provide a potential approach for cervical cancer therapeutics.

Human tissue samples

In total, 46 normal cervical epithelium tissues, 41 high-grade squamous intraepithelial lesions (HSILs), and 75 cervical cancer tissues were collected for qRT-PCR assays. Cervical cancer samples were obtained from patients who underwent radical hysterectomy, HSIL samples were obtained from patients who underwent colposcopy biopsy, and normal cervical tissues were collected from patients who underwent hysterectomy because of benign gynecological diseases. All the samples were collected from September 2015 to September 2020 at Women’s Hospital, Zhejiang University School of Medicine, China. Patients provided informed consent to obtain samples, and the study was subject to approval by the Hospital Ethical Committee. Patient information for all specimens is shown in Additional file 1: Table S1. Tissue samples were stored in RNA storage solution at 4°C and stored at -80°C until use.

Identification of HPV type

The HPV type of tissue samples was identified using the 21 HPV GenoArray Diagnostic Kit (Innovation Technologies of Hybribio). With the use of PCR principles to amplify extracted HPV DNA from cervical samples, amplified DNA amplicons were then hybridized with immobilized specific HPV probes on a HybriMem membrane under the patented flow-through hybridization technique. The enzyme immunoassay method was applied for color development to obtain test results. HPV-negative specimens used for RNA sequencing were amplified with the common primer my09/11 PCR of HPV L1 to further determine the HPV-negative status.

RNA extraction and qRT-PCR

RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. For tissue samples used to validate circRNA expression, the RNA integrity number (RIN) was assessed. qRT-PCR analyses were performed using the PrimeScript RT Reagent Kit and SYBR Premix Ex Taq (TaKaRa, Japan). All the primers are presented in Additional file 2: Table S2.

Cell line culture

The human cervical cancer cell line SiHa and the human embryonic renal cell line HEK293 were obtained from the American Type Culture Collection (ATCC, USA). The human cervical cancer CaSki cell line was obtained from the Cell Resource Center, Shanghai Institute of Life Sciences, Chinese Academy of Sciences (China), where it was tested and authenticated. It was not cultured continuously for more than 3 months. The SiHa cell line was cultured in DMEM (BI, Israel) supplemented with 10% FBS and maintained at 37°C in 5% CO₂, whereas the cervical cancer cell line CaSki and the human embryonic renal cell line HEK293 were cultured in RPMI-1640 (BI, Israel) containing 10% FBS.

RNA-fluorescence in situ hybridization (FISH) assay

CY3-labeled probes targeting the junction point sequence were used to visualize circCDKN2B-AS1 in situ (details of the probes are shown in Additional file 2: Table S2). SiHa and CaSki cells were plated onto slides and incubated overnight. Triton X-100 (0.5%) was added after washing the slides with PBS 3 times and then incubated at room temperature for 15 min. Then, the cells were fixed with 4% paraformaldehyde for 5 min and incubated with 100% ethanol. After discarding ethanol, the denatured probe was added to immerse the slides. Then, the slides were incubated at 37°C overnight. Nuclei were counterstained with DAPI-antifade solution, and images were taken on a laser confocal microscope (TCS SP2 AOBS).

RNaseR resistance assay

RNA samples were incubated with or without 2 U/µg RNaseR at 37°C for 15 min.

Northern blot assay

The probe targeting the junction point of circCDKN2B-AS1 was synthesized and labeled with digoxigenin (DIG, details of probes are shown in Additional file 2: Table S2). Fifteen micrograms of total RNA was separated by 1% denaturing formaldehyde gel electrophoresis at 25 V overnight and transferred to Hybond-N + membranes (Amersham, UK, RPN303B). Prehybridization was carried out at 50°C for 2 h in DIG Easy Hyb solution (Roche). Hybridization was performed at 50°C overnight. The membrane was washed stringently and blocked in blocking solution for 1 h. Blots were detected by anti-DIG antibody staining and recorded using X-ray films with the chemiluminescence substrate CSPD (Roche).

Western blot assay

Cellular proteins were extracted using lysis buffer. In total, 20 µg of total protein was separated using 10% SurePAGE (GenScript, USA, M00665) and transferred to PVDF membranes (Bio-Rad, USA, 1620177). The antibodies used were as follows: IMP3 (EMD Millipore), HK2 (Abcam), GAPDH (Diagbio) E-cadherin (CST) and β-Catenin (CST).

Gene knockdown and overexpression

The human circCDKN2B-AS1 linear sequence was inserted into the plasmid vector pLent-EF1a-circRNA-CMV-RFP-P2A-Puro, and lentiviruses stably overexpressing circCDKN2B-AS1 were constructed (Weizhen, Shandong, China). A nonfluorescent circCDKN2B-AS1 overexpression plasmid was constructed for apoptosis detection using the plasmid vector (Weizhen, Shandong, China). The full-length products and the junction point of the circCDKN2B-AS1-overexpressing products were determined by Sanger sequencing after RT-PCR. Short hairpin RNAs (shRNAs) specific to circCDKN2B-AS1 were inserted into the lentiviral vector GV334 (GeneChem, Shanghai, China). Stable cell lines were obtained by puromycin resistance. The human HK2 overexpression plasmid was purchased from GeneChem. All of the small interfering RNA (siRNA) and primer sequences used in this study are presented in Additional file 2: Table S2.

Cell viability, migration, invasion and apoptosis assays

In vitro cell viability was detected using Cell Counting Kit-8 assays (Dojindo, Japan). Transwell assays were performed at 48 h after transfection: 1 x 10⁴cells were resuspended in serum-free Opti-MEM medium, and cell invasion/migration was examined in Transwell cell culture chamber ﬁlters coated on the upper side with/without Matrigel (Corning Biocoat) as previously described [34].The wound healing assays we performed by using culture insert(ibidi, Germany, 80206) according to the instructions of the manufacturer. The apoptotic rates of SiHa and CaSki cells were determined at 72 h after transfection using an Annexin V-FITC/PI Apoptosis Kit (Mutisiences, China, AP101-100-kit). For cells transfected with circCDKN2B-AS1-overexpressing or empty plasmids, apoptosis was induced by incubation with serum-free medium for 6 h.

Extracellular acidification rate (ECAR)

Assays were performed using a Seahorse XFe96 analyzer (Seahorse Bioscience, Agilent) according to the manufacturer’s instructions. The ECAR was measured using a Seahorse XF Glycolytic Rate Assay Kit (Seahorse Bioscience, Agilent, 103044-100) and a Seahorse XF Glycolysis Stress Test Kit (Seahorse Bioscience, Agilent, 103020-100). The glycolytic capacities of cells were analyzed with the Seahorse XF Glycolysis Rate/Stress Test Report Generator package. The %PER from glycolysis was calculated by subtracting the acidification from CO₂ produced by the mitochondria.

Inhibition of glycolysis in cells

We inhibited the glycolysis level of cervical cancer cells in vitro by adding 2-DG (5mM, Seahorse Bioscience, Agilent,103044-100).

Biotin-labeled RNA pull-down assay and mass spectrometry analysis

Biotin-labeled RNA pull-down assays were performed using a Pierce™ Magnetic RNA-Protein Pull-down Kit (Thermo Scientific, USA, 20164). Briefly, cell lysates were prepared using standard lysis buffers (Thermo Scientific, USA). Biotin-labeled DNA probes (details of probes are shown in Additional file 2: Table S2) were incubated with streptavidin magnetic beads for 30 minutes at room temperature with agitation. Cell lysates were incubated with streptavidin magnetic beads at 4°C overnight. Magnetic beads were thoroughly washed. The proteins bound to the magnetic bead were separated with SurePAGE (GenScript, USA, M00665), and the region of interest was excised and subjected to mass spectrometry analysis (Lumingbio, Shanghai, China).

RNA immunoprecipitation (RIP) assay

RIP assays were performed with a Magna RIP™ Quad RNA-Binding Protein Immunoprecipitation Kit (Sigma-Aldrich, USA, 17-704) according to the instructions of the manufacturer.

Determination of mRNA half-life

To assess the half-life of HK2 mRNA, actinomycin-D (5 µg/ml, Sigma, A4262) was added to block mRNA synthesis. Total RNA was collected at different time points and subjected to RT-qPCR analysis. The expression level of HK2 mRNA was normalized to that of 18S rRNA and plotted as a percentage of the value at time 0 (set at 100%).

Design and synthesis of inhibitory peptides

The inhibitory peptides for blocking the interaction between circCDKN2B-AS1 and IMP3 were designed and synthesized by ChinaPeptides Co. Ltd. (Shanghai, China). The inhibitory peptides were synthesized by linking the biotin-labeled 11-amino acid cell-penetrating peptide YGRKKRRQRRR of the Tat protein transduction domain with the core amino acids of IMP3 at the N-terminus (purity greater than 95%).

Biotin peptide pull-down assay

Biotin-labeled peptides were incubated with streptavidin magnetic beads (Thermo Scientific, USA, 88817) for 4 h at 4°C. Then, the bead-peptide complex was incubated with total RNA at 4°C overnight. Beads were extensively washed, and RNAs that precipitated were measured by qRT-PCR.

Nude mice xenograft experiments

This study was conducted in full accordance with the ARRIVE guidelines for reporting animal research [35]. All animal experiments were approved by the Animal Ethical and Welfare Committee of Zhejiang Chinese Medical University under an Affidavit of Approval of Animal Ethical and Welfare license (No. IACUC-20190128-01) and in accordance with the Animals (Scientific Procedures) Act, 1986 (UK) (amended 2013). BALB/c nude mice (4 weeks old±1week, weighing 18±5 grams, female) were ordered from Shanghai SLAC Laboratory Animal Company, Ltd. (China). Nude mice were randomly divided into two groups and anesthetized by inhalation of isoflurane gas (2-3%) also under the effect of muscle relaxation. In total, 10⁷ stable SiHa cells were resuspended in 100 μL PBS and injected subcutaneously under the left armpit of 4-week-old nude mice (n = 10 per group). The length and width of the subcutaneous tumor were measured once a week for 5 weeks, and the volume of the subcutaneous tumor was calculated according to the following formula: volume (cm³) = (length×width²)/2. Five weeks after injection, 6 mice in each group were randomly selected and sacrificed by intravenous injection of pentobarbital (100 mg/kg), and subcutaneous tumors were harvested. A portion of the tumor was stored in RNA storage solution to extract RNA, and the remainder of the tumor was fixed with paraformaldehyde for immunohistochemistry and hematoxylin and eosin (HE) staining. The other 4 mice in each group were subjected to microPET/CT scans. Another 20 nude mice (4 weeks old) were randomly divided into 4 groups (n=5 per group) and injected subcutaneously under the left armpit with 10⁷ stable SiHa cells. siRNAs targeting a negative control (NC) or HK2 were mixed with FECT and Opti-MEM medium and injected into the subcutaneous tumor every week. The volume of the subcutaneous tumor was calculated according to the following formula: volume (cm³) = (length×width²)/2. Four weeks after the injection of cells, mice were sacrificed, and subcutaneous tumors were harvested, weighed, and stored at -80°C for RNA extraction, lactic acid (LA)/pyruvic acid (PA) detection and RIP assays.

MicroPET/CT scans in nude mice

MicroPET/CT imaging of nude mice was performed using the SuperArgus PET/CT 4R system (Sedecal, Spain) five weeks after cell injection. Briefly, 4 tumor-bearing nude mice in each group were anesthetized by inhalation of isoflurane gas (2-3%). Then, nude mice were injected with 7.4 MBq (200 μCi) of ¹⁸F-FDG via the abdominal cavity and scanned at 120 min after injection. Nude mice were subjected to a 10-min microCT scan and then to a 15-min microPET scan. Scans were performed under isoflurane (1-1.5%) inhalation to maintain anesthesia. Images of nude mice were reconstructed manually. The PET/CT camera Vista explore was used to calculate the average standardized uptake volume (SUV_AVG). Then, the mice were sacrificed by intravenous injection of pentobarbital (100 mg/kg).

Digestion of subcutaneous tumors

Subcutaneous tumors were digested with 10 mg type 1 collagenase (Worthington, 49E19342, USA) at 37°C for 60 minutes. The digestion products were centrifuged at 1200 r for 5 minutes, and the precipitate was subjected further to RIP assays and RNA extraction.

Detection of LA and PA

The levels of LA and PA in 10 mg subcutaneous tumors from each mouse were detected by using an LA detection kit (Solarbio, China, BC2230) and a PA detection kit (Solarbio, China, BC2200).

Zebrafish embryo model

Zebrafish (AB, wild type [WT], 48-h embryo) embryos were obtained from the Core Facilities, Zhejiang University, School of Medicine, cultured in Petri dishes with an appropriate amount of water and placed in a 28°C incubator under constant temperature and light. Then, 0.003% 1-phenyl-2-thiourea (PTU) (Sigma) was added within 24 h after the embryo was produced. In total, 200 SiHa cells were harvested after transfection and resuspended in PBS. Plasma membranes were stained with red-fluorescent Alexa Fluor® 594 wheat germ agglutinin of the Image-iT™ LIVE Plasma Membrane and Nuclear Labeling Kit (Molecular Probes, USA, I34406). Zebrafish embryos were anesthetized with 0.016% tricaine (Sigma). Then, the cells were injected into the yolk sac of zebrafish embryos. Then, we obtained images with a fluorescence microscope at 48 h after injection. The area of red fluorescence in the zebrafish yolk sac was calculated using ImageJ software [33]. Next, the zebrafish embryos were anesthetized with an overdose of tricaine.

Immunohistochemical staining

Immunohistochemical staining with an antibody specific to HK2 (Proteintech, USA) and quantitative evaluation were performed as previously described [36].

Statistical analysis

Statistical significance was calculated by using GraphPad Prism 5 software. All experiments were repeated at least three times.

CircCDKN2B-AS1 is upregulated in cervical cancer and precancerous tissues

To explore circRNA role in cervical cancer development, we analyzed aberrantly expressed circRNAs among 6 HPV16 positive cervical squamous-cell carcinoma, 6 HPV16 positive normal cervical epithelium, and 6 HPV negative normal cervical epithelium tissues in our circRNA sequencing data, available in the Gene Expression Omnibus (GEO) database (GSE147009) (Figure 1A-B). We selected 4 circRNAs (circCDKN2B-AS1, circFNACB, circEPSTI1, and circOBSL1) that were upregulated or downregulated in HPV16 positive cervical squamous-cell carcinoma compared with both HPV16 positive and negative normal cervical epithelium samples for further validation.

The expression levels of these 4 circRNAs were detected in 25 HPV-negative normal cervical tissues, 20 HPV16-positive HSIL, and 21 HPV16-positive cervical cancer by qRT-PCR. Among these circRNAs, a novel circRNA formed by the circularization of exon 2, exon 3 and exon 4 of CDKN2B-AS1 was significantly aberrantly expressed in both HSIL tissues and cancer tissues compared with normal cervical tissues (Figure 1C-F). Products amplified by divergent primers of these four circRNAs were confirmed by Sanger sequencing (Figure 1G, Additional file 3: Figure S1A-C, Additional file 4, Table S3).

CircCDKN2B-AS1 was formed by the circularization of exon 2, exon 3 and exon 4 of CDKN2B-AS1 (Figure 1G). Linear CDKN2B-AS1 is a lncRNA that has been reported to facilitate tumorigenesis in several cancers, such as breast cancer and gastric cancer. To verify the endogenous circularization of exons 2-4 of CDKN2B-AS1, we designed convergent primers and divergent primers to amplify the linear and circular forms of CDKN2B-AS1. PCR using the convergent primers amplified products in both genomic DNA and cDNA formed by reverse transcription, but the divergent primers only amplified products in cDNA (Additional file 3: Figure S1D). Products amplified by divergent and convergent primers were confirmed by Sanger sequencing (Additional file 4, Table S3). RT-PCR products amplified by two pairs of full-length primers designed in two different points of circCDKN2B-AS1 were used to amplify full-length circCDKN2B-AS1 (Figure 1H), and the sequence of the products was confirmed by Sanger sequencing (Additional file 4, Table S3). CircCDKN2B-AS1 was more resistant to RNaseR digestion than the canonical form of CDKN2B-AS1 (Figure 1I).

We noted that backsplicing of the novel 380-nt circCDKN2B-AS1 also formed a 1719 nt isoform, as reported in circBase (hsa_circ_0008796). However, the length of this circCDKN2B-AS1 annotated in CircPrimer software was 380 nt instead of 1719 nt [37]. A previous study reported the alternative splicing of circRNAs and confirmed the length of circRNAs using Northern blot analysis [38]. Northern blotting results showed that the length of the main isoform of circCDKN2B-AS1 in cervical cancer cell lines was 380 nt (Additional file 3: Figure S1E).

In addition, we detected the level of circCDKN2B-AS1 expression by qRT-PCR assay in 46 normal cervical tissues, 41 HSIL, and 75 cervical cancer tissues regardless HPV positive or negative, and found that circCDKN2B-AS1 expression was upregulated in HSIL and cervical cancer, and the level of circCDKN2B-AS1 expression in cervical cancer and HSIL tissues was 15.16 and 5.48-fold, respectively, compared to normal cervical tissues (Figure1J).Then we analyzed the correlation between the level of circCDKN2B-AS1 expression and clinical prognostic parameters in 75 cervical cancer patients, and found that the level of circCDKN2B-AS1 expression was positively correlated with poorer prognostic parameters including tumor size, FIGO stage, deep stromal invasion, and lymph node metastasis (Table 1).

CircCDKN2B-AS1 facilitates the malignant phenotype of cervical cancer cells

To investigate the biological function of circCDKN2B-AS1 in the malignant phenotype of cervical cancer cells, we knocked down circCDKN2B-AS1 by using two specific siRNAs in SiHa and CaSki cells and found that circCDKN2B-AS1 downregulation did not affect the expression of the linear forms of CDKN2B-AS1 (Figure 2A) but inhibited cellular proliferation, migration and invasion, and promoted cellular apoptosis (Figure 2B-F,Additional file 5: FigureS2).

We next observed the role of circCDKN2B-AS1 overexpression in SiHa and CaSki cells using a constructed plasmid to upregulate circCDKN2B-AS1 expression in SiHa and CaSki cells (Additional file 6: Figure S3A). Sanger sequencing followed by RT-PCR verified that the sequence of circCDKN2B-AS1 was correct (Additional file 4, Table S3). We found that circCDKN2B-AS1 overexpression promoted cellular proliferation, migration and invasion and repressed the cellular apoptotic rate (Additional file 6: Figure S3B-E). Collectively, these data are consistent with the notion that circCDKN2B-AS1 facilitates the malignant phenotype of cervical cancer and acts as a tumor promoter in cervical cancer development.

CircCDKN2B-AS1 cooperates with IMP3 to promote glycolysis in cervical cancer

RNA-FISH showed that circCDKN2B-AS1 was mainly localized in the cytoplasm of SiHa and CaSki cells (Figure 3A), implying that circCDKN2B-AS1 probably regulated gene expression at the posttranscriptional level by sponging microRNAs or RBPs.

Then, we performed RNA pull-down assays to determine the RBPs that interact with circCDKN2B-AS1. Mass spectrometry analysis after RNA pull-down assays revealed 35 differential proteins between the oligo probe and the circCDKN2B-AS1-specific probe (Figure 3B, Additional file 7: Figure S4A-B, Additional file 8: Table S4). Then, we utilized the RBPmap database to predict candidate RBPs that interact with circCDKN2B-AS1. The results showed 43 proteins that may interact with circCDKN2B-AS1 (Additional file 9: Table S5). Of these, IMP2 and IMP3 were included in both the mass spectrometry analysis and the RBPmap database (Figure 3B). Moreover, IMP3 showed higher coverage [%] in the mass spectrometry results (Additional file 8: Table S4) and higher Z-scores in the RBPmap results (Additional file 9: Table S5). Thus, IMP3 protein was more possibly to be interacted with circCDKN2B-AS1. RIP and RNA pull-down assays confirmed the interaction between circCDKN2B-AS1 and IMP3 in SiHa and CaSki cells (Figure 3C-D, Additional file 7: Figure S4C). These data suggest that circCDKN2B-AS1 interacts with the IMP3 protein in cervical cancer cells.

A previous study showed that the IMP3 protein was a protein partner of RNAs or a protein that regulates glycolysis levels in cancer cells [39-40]. Thus, we investigated the role of circCDKN2B-AS1 in modulating glycolysis in cervical cancer cells. We found that the ECAR and %PER from glycolysis were decreased with circCDKN2B-AS1 knockdown in SiHa and CaSki cells (Figure 3E-F, Additional file 7: Figure S4D), suggesting that circCDKN2B-AS1 knockdown weakens glycolysis in cervical cancer cells. We further used 2-DG to inhibit glycolysis in cervical cancer cells, and found that the inhibition of glycolysis decreased cellular proliferation, migration and invasion, and increased apoptotic rate in cervical cancer cells (Figure 3G-I).

We also found that IMP3 protein expression was not affected by circCDKN2B-AS1 knockdown in SiHa and CaSki cells; however, the amount of IMP3 protein recruited by circCDKN2B-AS1 was decreased with circCDKN2B-AS1 knockdown (Figure 3J).

Additionally, we designed a mutant circCDKN2B-AS1 overexpression vector with the mutation in the binding site (AAAAACACACA) predicted by RBPmap. We verified the interaction between circCDKN2B-AS1 and IMP3 through the specific binding site in HEK293 cells, and we found that IMP3 was recruited by upregulated WT circCDKN2B-AS1 but not by upregulated mutant circCDKN2B-AS1 in HEK293 cells (Figure 3K).Consistent with these results, WT circCDKN2B-AS1 overexpression promoted glycolysis in SiHa and CaSki cells, but this phenomenon was reversed by mutating the binding site between circCDKN2B-AS1 and IMP3 (Figure 3L, Additional file 7: Figure S4E-F). Together, our results suggest that circCDKN2B-AS1 cooperates with IMP3 to promote glycolysis of cervical cancer cells.

CircCDKN2B-AS1 regulates the stability of HK2 mRNA via recruiting IMP3 protein

Then, we compared the differentially expressed mRNAs obtained with our RNAseq analysis (available in GSE147010) and the CLIP-Seq analysis of IMP3 in the POSTAR2 database (http://lulab.life.tsinghua.edu.cn/postar/) and found that HK2, ENO1, PGK1 and PKM were the potential downstream molecules of circCDKN2B-AS1 and IMP3. Among these molecules, HK2 mRNA was the most significantly downregulated in SiHa cells following circCDKN2B-AS1 knockdown (Additional file 10: Figure S5A). While a previous study showed that the IMP3 protein could be recruited to stabilize HK2 mRNA [38]. We thus suspected that HK2 was the common downstream molecule of circCDKN2B-AS1 and IMP3. We found that the 3’ untranslated region (UTR) of HK2 mRNA was significantly recruited by the IMP3 protein (Figure 4A). Accordingly, knockdown of IMP3 reduced the expression of HK2 mRNA in SiHa cells (Figure 4B). We further observed the significantly downregulated expression of both HK2 mRNA and protein in SiHa and CaSki cells following circCDKN2B-AS1 knockdown (Figure 4C) and that the stability of HK2 mRNA was lower in SiHa and CaSki cells following circCDKN2B-AS1 knockdown than in control cells (Figure 4D).

Furthermore, we constructed a SiHa cell line with stable circCDKN2B-AS1 knockdown using shRNAs and with circCDKN2B-AS1 overexpression using lentivirus, respectively. We found that cells with stable circCDKN2B-AS1 knockdown exhibited suppressed glycolysis compared to control cells, but this suppression of glycolysis was reversed by HK2 overexpression (Figure 4E, Additional file 10: Figure S5B), vice versa (Figure4F).

In addition, the percentage of HK2 mRNA 3’UTR recruited to IMP3 were decreased after circCDKN2B-AS1 knockdown in SiHa cells (Figure 4I) and increased after circCDKN2B-AS1 overexpression in HEK293 cells, but mutant circCDKN2B-AS1 didn’t have such effects (Figure 4J). Taken together, these results suggest that circCDKN2B-AS1 cooperates with IMP3 to regulate the stability of HK2 mRNA.

Inhibitory peptides repress glycolysis by blocking the interaction between the circCDKN2B-AS1 and IMP3 protein in cervical cancer cells

It has been reported that inhibitory peptides can block the interaction between specific circRNAs and RBPs [41, 42]. We used an inhibitory peptide to block the binding between circCDKN2B-AS1 and IMP3. Based on the catRapid and RNABindRPlus prediction results, we designed cell-penetrating inhibitory peptides that include the core amino acids of IMP3 and RRM1, named IIP (Figure 5A), and confirmed that IIP could interact with circCDKN2B-AS1 in SiHa and CaSki cells (Figure 5B). Accordingly, we found that the interaction between circCDKN2B-AS1 and IMP3 became weaker after SiHa and CaSki cells were treated with IIP (Figure 5C). Consequently, the binding between IMP3 and the 3’UTR of HK2 mRNA (Figure 5C) and HK2 expression were reduced (Figure 5D-E). Then, we detected the glycolysis capabilities of IIP-treated cells and found that glycolysis in SiHa and CaSki cells treated with IIP was also decreased (Figure 5F), suggesting that IIP represses glycolysis by blocking the interaction between circCDKN2B-AS1 and the IMP3 protein in cervical cancer cells.

CircCDKN2B-AS1 facilitates aerobic glycolysis and cancer progression via interacting with IMP3 and HK2 mRNA in vivo.

We constructed nude mouse models to observe circCDKN2B-AS1 effect on aerobic glycolysis and progression of cervical cancer, and found that the growth of xenografted tumors was markedly retarded (Figure 6A) and the levels of both HK2 mRNA and protein expression were lower, with a significantly positive correlation between HK2 mRNA and circCDKN2B-AS1 expression (Figure 6B-C), in 6 pairs of nude mice with depletion of circCDKN2B-AS1 compared to controls. Another 4 pairs of tumor-bearing nude mice were subjected to ¹⁸F-FDG microPET/CT imaging to assess the level of glucose metabolism in xenograft tumors. The average SUV in tumor tissues with circCDKN2B-AS1 knockdown was lower than that in control tissues (Figure 6D). Then, we injected SiHa cells with upregulated circCDKN2B-AS1 into the yolk sac of zebrafish embryos and found an increased survival rate of cells with upregulated circCDKN2B-AS1 compared with control cells (Figure 6E).

Further, we constructed additional 4 groups of nude mouse models with subcutaneous tumors that were circCDKN2B-AS1 overexpressed or not, and found that circCDKN2B-AS1 overexpression promoted tumor growth (Figure 6F)), elevated LA and PA contents (Figure 6G), elevated the level of HK2 mRNA expression of in tumor tissues with a positive correlation between them (Figure 6H), and facilitated recruitment rate of HK2 mRNA 3’UTR to IMP3 (Figure 6I, Additional file 11: Figure S6).Then, siRNAs targeting NC or HK2 were injected into the subcutaneous tumor every week from day 7 after SiHa cell injection. We found that HK2 knockdown revised tumor growth and the LA/PA content promoted by circCDKN2B-AS1 overexpression (Figure 6F, Figure 6I). Thus, our results suggest that circCDKN2B-AS1 facilitates aerobic glycolysis and the progression of cervical cancer via interacting with IMP3 and HK2 mRNA in vivo.

CircRNAs are noncoding RNAs that are more stable than their linear isoforms. It has been reported that circRNAs are dysregulated and play a tumor suppressive or oncogenic role in cancers, such as breast cancer, ovarian cancer and bladder cancer [16, 43, 44]. Here, we identified a circRNA, circCDKN2B-AS1, that is upregulated in HPV16-positive cervical cancer and precancerous tissues. Furthermore, we showed that circCDKN2B-AS1 acts as a tumor promoter in cervical cancer development.

As circCDKN2B-AS1 is mainly localized in the cytoplasm of cervical cancer cells, we suspected that circCDKN2B-AS1 binds to microRNAs or RBPs, as most circRNAs do [44–48]. Then, we confirmed that IMP3 was the protein partner of circCDKN2B-AS1 by conducting a series of RNA-protein binding validation experiments. Our results indicated that HK2 was a common downstream molecule of circCDKN2B-AS1 and IMP3.

Thus, we hypothesized that circCDKN2B-AS1 stabilizes HK2 mRNA and facilitates aerobic glycolysis in cervical cancer by recruiting the IMP3 protein to the 3’UTR of HK2 mRNA (Figure 7). Supporting this hypothesis, we showed that knockdown of circCDKN2B-AS1 weakened the binding of IMP3 to the 3’UTR of HK2 mRNA, reduced the stability of HK2 mRNA, decreased the expression of HK2 and suppressed aerobic glycolysis. In contrast, overexpression of circCDKN2B-AS1 had the opposite effects. However, overexpression of circCDKN2B-AS1 lacking the IMP3 binding site did not yield the above effects. Next, we showed that circCDKN2B-AS1 knockdown decreased the level of glucose metabolism and the expression of HK2 in xenograft tumors. Additionally, there was a significant positive correlation between the expression of circCDKN2B-AS1 and HK2 mRNA in xenografted tumors regardless circCDKN2B-AS1 knockdown or overexpression, suggesting that circCDKN2B-AS1 facilitates cellular malignant phenotypes of cervical cancer via interacting with IMP3 and HK2 mRNA in vivo.

Furthermore, we designed an inhibitory peptide, IIP, to block the interaction between circCDKN2B-AS1 and IMP3. We designed IIP based on the evidence that RRM1 of IMP3 is the key domain involved in binding CA-rich RNA elements [22, 23]. Our results showed that the inhibitory peptide IIP blocked the interaction between circCDKN2B-AS1 and IMP3, weakened the binding of IMP3 to the 3’UTR of HK2 mRNA, decreased the expression of HK2 and suppressed aerobic glycolysis.

There were some limitations of our study. For instance, it needs to be determined whether some of the potential circCDKN2B-AS1 binding proteins, such as IMP2, are involved in the gene regulatory network that promotes aerobic glycolysis in cervical cancer cells. In addition, the role of circCDKN2B-AS1 in other tumors and other diseases also needs to be investigated in the future.

In the study, we found, for the first time to our knowledge, that circCDKN2B-AS1 acted as a promoter that facilitated aerobic glycolysis by sponging IMP3 protein to stabilize HK2 mRNA, consequently promoted malignant phenotypes in cervical cancer in vitro and in vivo, and the synthesizing inhibitory peptide IIP could effectively block the interaction between circCDKN2B-AS1 and IMP3 protein. Our findings may provide a new approach for cervical cancer therapeutics

ATCC: American type culture collection; CircRNAs: Circular RNAs; CDKN2B-AS1:CDKN2B antisense RNA 1;CircCDKN2B-AS1:Circular isoform of CDKN2B antisense RNA 1; HPV: human papillomavirus; IMP3/IGF2BP3:insulin like growth factor 2 mRNA binding protein 3; HK2: Hexokinase 2; ENO1:enolase 1;PKG1:Phosphoglycerate kinase 1;PKM:pyruvate kinase M 1/2; qRT-PCR: Quantitative real-time polymerase chain reaction;3'-UTR:3’ end Untranslated Regions; FCM: Flow cytometry analysis; CCK8: Cell counting kit-8; MS: Mass spectrometry; RNA-Seq: Genome-wide transcriptome analysis; FITC: fluorescein isothiocyanate; GAPDH:Glyceraldehyde-3-phosphate dehydrogenase; RIP: RNA Binding Protein Immunoprecipitation; PET: positron emission tomography; CT: Computed Tomography; RNA FISH: RNA Fluorescence in situ hybridization; mRNA: Message RNA; CLIP-seq: crosslinking-immunprecipitation and high-throughput sequencing; NC: Negative control; IIP: Inhibitory IMP3 peptide; CTL: Control peptide; shRNA: short hairpin RNA; siRNA: Small interference RNA; cDNA: Complementary DNA; gDNA: Genomic DNA; LA: Lactic acid; PA: Pyruvic acid; EMT: epithelial-mesenchymal transition.

Ethics approval and consent to participate

All animal experiments were carried out in accordance with Animal Research: Reporting In Vivo Experiments (ARRIVE) guidelines and were approved by the Animal Ethical and Welfare Committee of Zhejiang Chinese Medical University. All the samples were collected at Women’s hospital, Zhejiang University School of Medicine, China. Patients were informed consent for obtaining samples and the study was subject to approval by the Hospital Ethical Committee.

Consent for publication

Not applicable.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article (and its additional files)

Competing interests

The authors declare no conflicts of interest.

Funding

This work was supported by the National Key Research and Development Program of China (No. 2016YFC1302900), the Key research and development plan in Zhejiang

province (No. 2020C03025), the Fundamental Research Funds for the Central Universities (2017XZZX001-01, 2019QNA7035), National Natural Science Foundation of China (81702552,82072855), Medical and Health science and Technology project of Zhejiang Province (2017KY099), Postdoctoral Science Foundation(2017M621954) and the Natural Science Foundation of Zhejiang Province (No. LQ20H160049).

Authors' contributions

XX and XW Conceptualized and designed the study, revised the manuscript, supervised and administrated the study. YZ, LZ and SY performed most of the experiments, analyzed the data, wrote the manuscript. YZ, LZ, LW collected and classified the human tissue samples and contributed to the RT-PCR and real-time PCR experiments, biotin-labeled RNA pull-down assay and RIP assay. SY, LZ and YZ performed primer design, RNA-FISH, CCK8 assays, FCM assay and ECAR assay. JZ, YC and TZ performed the transwell assays and western blots. JX, LX, YL (Yuwan Liu) performed most of the xenograft experiments. ZJ, YL (Yang Li), XC and WL analyzed the data, provided technical and theoretical support. All authors read and approved the final manuscript.

Acknowledgements

We thank Prof Mei Tian, Prof Pengfei Xu, Prof Xiaohang Yang and Zhihong Xu (School of medicine, Zhejiang University, Hangzhou, China) for providing technical and scientific research guidance. Thanks for the technical support by the Core facilities，Zhejiang University, School of medicine.

Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2017. CA Cancer J Clin. 2017; 67:7–30.
Small W, Bacon MA, Bajaj A, Chuang LT, Fisher BJ, Harkenrider MM, et al. Cervical cancer: A global health crisis. Cancer. 2017;123: 2404–12.
Zhao F, Qiao Y. Cervical cancer prevention in China: a key to cancer control. The Lancet. 2019; 393:969–70.
Cohen PA, Jhingran A, Oaknin A, Denny L. Cervical cancer. The Lancet. 2019; 393:169–82.
Tewari KS, Sill MW, Long HJ, Penson RT, Huang H, Ramondetta LM, et al. Improved survival with bevacizumab in advanced cervical cancer. N Engl J Med. 2014; 370:734–43.
Sundström K, Elfström KM. Advances in cervical cancer prevention: Efficacy, effectiveness, elimination? PLoS Med. 2020; 17: e1003035.
Noone AM, Howlader N, Krapcho M, Miller D, Brest A, Yu M, Ruhl J, Tatalovich Z, Mariotto A, Lewis DR, Chen HS, Feuer EJ, Cronin KA (eds). SEER Cancer Statistics Review, 1975-2015, National Cancer Institute. Bethesda, MD, November 2017 SEER data submission, posted to the SEER web site, April 2018.
Della Bella E, Menzel U, Basoli V, Tourbier C, Alini M, Stoddart MJ. Differential Regulation of circRNA, miRNA, and piRNA during Early Osteogenic and Chondrogenic Differentiation of Human Mesenchymal Stromal Cells. Cells 2020.
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013; 495:384–8.
Kristensen LS, Andersen MS, Stagsted LVW, Ebbesen KK, Hansen TB, Kjems J. The biogenesis, biology and characterization of circular RNAs. Nat Rev Genet. 2019; 20:675–91.
Chen L-L. The biogenesis and emerging roles of circular RNAs. Nat Rev Mol Cell Biol. 2016; 17:205–11.
Geng Y, Jiang J, Wu C. Function and clinical significance of circRNAs in solid tumors. J Hematol Oncol. 2018; 11:98.
Zhang M, Xin Y. Circular RNAs: a new frontier for cancer diagnosis and therapy. J Hematol Oncol. 2018; 11:21.
Chen L, Zhang S, Wu J, Cui J, Zhong L, Zeng L, Ge S. circRNA_100290 plays a role in oral cancer by functioning as a sponge of the miR-29 family. Oncogene. 2017; 36:4551–61.
Xiao Y, Liu G, Sun Y, Gao Y, Ouyang X, Chang C, et al. Targeting the estrogen receptor alpha (ERα)-mediated circ-SMG1.72/miR-141-3p/Gelsolin signaling to better suppress the HCC cell invasion. Oncogene 2020.
Wu J, Jiang Z, Chen C, Hu Q, Fu Z, Chen J, et al. CircIRAK3 sponges miR-3607 to facilitate breast cancer metastasis. Cancer Lett. 2018; 430:179–92.
Chen, X., Jiang, J., Zhao, Y. et al. Circular RNA circNTRK2 facilitates the progression of esophageal squamous cell carcinoma through up-regulating NRIP1 expression via miR-140-3p. J Exp Clin Cancer Res 39, 133 (2020).

Chen J, Liu G, Wu Y, Ma J, Wu H, Xie Z, et al. CircMYO10 promotes osteosarcoma progression by regulating miR-370-3p/RUVBL1 axis to enhance the transcriptional activity of β-catenin/LEF1 complex via effects on chromatin remodeling. Mol Cancer. 2019; 18:150.
Pandey PR, Yang J-H, Tsitsipatis D, Panda AC, Noh JH, Kim KM, et al. circSamd4 represses myogenic transcriptional activity of PUR proteins. Nucleic Acids Res 2020.
Li Q, Wang Y, Wu S, Zhou Z, Ding X, Shi R, et al. CircACC1 Regulates Assembly and Activation of AMPK Complex under Metabolic Stress. Cell Metab. 2019; 30:157-173.e7.
Hanniford D, Ulloa-Morales A, Karz A, Berzoti-Coelho MG, Moubarak RS, Sánchez-Sendra B, et al. Epigenetic Silencing of CDR1as Drives IGF2BP3-Mediated Melanoma Invasion and Metastasis. Cancer Cell. 2020; 37:55-70. e15.
Min Jia, Heinz Gut, and Jeffery A. et al. Structural basis of IMP3 RRM12 recognition of RNA. RNA 24:1659–1666; October 18, 2019.
Schneider T, Hung L-H, Aziz M, Wilmen A, Thaum S, Wagner J, et al. Combinatorial recognition of clustered RNA elements by the multidomain RNA-binding protein IMP3. Nat Commun. 2019; 10:2266.
Hentze MW, Castello A, Schwarzl T, Preiss T. A brave new world of RNA-binding proteins. Nature Reviews Molecular Cell Biology. 2018; 19:327–41.
Lis P, Dyląg M, Niedźwiecka K, Ko YH, Pedersen PL, Goffeau A, Ułaszewski S. The HK2 Dependent “Warburg Effect” and Mitochondrial Oxidative Phosphorylation in Cancer: Targets for Effective Therapy with 3-Bromopyruvate. Molecules 2016.
Feng J, Li J, Wu L, Yu Q, Ji J, Wu J, et al. Emerging roles and the regulation of aerobic glycolysis in hepatocellular carcinoma. Journal of Experimental & Clinical Cancer Research. 2020; 39:126.

Chen J, Yu Y, Li H, Hu Q, Chen X, He Y, et al. Long non-coding RNA PVT1 promotes tumor progression by regulating the miR-143/HK2 axis in gallbladder cancer. Mol Cancer. 2019; 18:33.
Ward PS, Thompson CB. Metabolic reprogramming: a cancer hallmark even warburg did not anticipate. Cancer Cell. 2012;21: 297–308.
Mayers JR, Vander Heiden MG. Nature and Nurture: What Determines Tumor Metabolic Phenotypes? Cancer Res. 2017; 77:3131–4.
Frezza, C. Metabolism and cancer: the future is now. Br J Cancer 122, 133–135 (2020).
Jiang, Wei; He, Tiancong; Liu, Shuai; Zheng, Yingying; Xiang, Libing; Pei, Xuan et al: The PIK3CA E542K and E545K mutations promote glycolysis and proliferation via induction of the β-catenin/SIRT3 signaling pathway in cervical cancer. In Journal of hematology & oncology ,2018;11(1), p.139.
Zhou, L., Wu, N., Liang, Y. et al. Monitoring response to gefitinib in nude mouse tumor xenografts by ¹⁸F-FDG microPET-CT: correlation between ¹⁸F-FDG uptake and pathological response. World J Surg Onc 2015,13, 111.

Astone M, Dankert EN, Alam SK, Hoeppner LH. Fishing for cures: The alLURE of using zebrafish to develop precision oncology therapies. NPJ Precis Oncol. 2017
Ramirez-Moya J, Wert-Lamas L, Riesco-Eizaguirre G, Santisteban P. Impaired microRNA processing by DICER1 downregulation endows thyroid cancer with increased aggressiveness. Oncogene. 2019; 38: 5486–99.
Kilkenny C, Browne WJ, Cuthill IC, Emerson M, Altman DG (2010) Improving bioscience research reporting: the ARRIVE guidelines for reporting animal research. PLoS Biol 8: e1000412.
Xu J, Li Y, Wang F, Wang X, Cheng B, Ye F, et al. Suppressed miR-424 expression via upregulation of target gene Chk1 contributes to the progression of cervical cancer. Oncogene. 2013; 32:976–87.
Shanliang Zhong, Jinyan Wang, Qian Zhang, Hanzi Xu and Jifeng Feng. CircPrimer: a software for annotating circRNAs and determining the specificity of circRNA primers. BMC Bioinformatics. 2018;19(1):292.

Yang Y, Gao X, Zhang M, Yan S, Sun C, Xiao F, et al. Novel Role of FBXW7 Circular RNA in Repressing Glioma Tumorigenesis. J Natl Cancer Inst 2018.
An M-X, Li S, Yao H-B, Li C, Wang J-M, Sun J, et al. BAG3 directly stabilizes Hexokinase 2 mRNA and promotes aerobic glycolysis in pancreatic cancer cells. J Cell Biol. 2017; 216:4091–105.
Fei Ma, Xu Liu, Shibo Zhou, Wenjie Li, Chunxiao Liu, Michelle Chadwick, Cheng Qian, Long non-coding RNA FGF13-AS1 inhibits glycolysis and stemness properties of breast cancer cells through FGF13-AS1/IGF2BPs/Myc feedback loop, Cancer Letters, 2019, Volume 450:63-75,

Chen Y, Yang F, Fang E, Xiao W, Mei H, Li H, et al. Circular RNA circAGO2 drives cancer progression through facilitating HuR-repressed functions of AGO2-miRNA complexes. Cell Death Differ. 2019; 26:1346–64.
Li H, Yang F, Hu A, Wang X, Fang E, Chen Y, et al. Therapeutic targeting of circ-CUX1/EWSR1/MAZ axis inhibits glycolysis and neuroblastoma progression. EMBO Mol Med. 2019; 11: e10835.
Yang C, Yuan W, Yang X, Li P, Wang J, Han J, et al. Circular RNA circ-ITCH inhibits bladder cancer progression by sponging miR-17/miR-224 and regulating p21, PTEN expression. Mol Cancer. 2018; 17:19.
Zhang L, Zhou Q, Qiu Q, Hou L, Wu M, Li J, et al. CircPLEKHM3 acts as a tumor suppressor through regulation of the miR-9/BRCA1/DNAJB6/KLF4/AKT1 axis in ovarian cancer. Mol Cancer. 2019; 18:144.
Chen R-X, Chen X, Xia L-P, Zhang J-X, Pan Z-Z, Ma X-D, et al. N6-methyladenosine modification of circNSUN2 facilitates cytoplasmic export and stabilizes HMGA2 to promote colorectal liver metastasis. Nat Commun. 2019; 10:4695.
Zheng Q, Bao C, Guo W, Li S, Chen J, Chen B, et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat Commun. 2016; 7:11215.
Du WW, Fang L, Yang W, Wu N, Awan FM, Yang Z, Yang BB. Induction of tumor apoptosis through a circular RNA enhancing Foxo3 activity. Cell Death Differ. 2017; 24:357–70.
Du WW, Yang W, Chen Y, Wu Z-K, Foster FS, Yang Z, et al. Foxo3 circular RNA promotes cardiac senescence by modulating multiple factors associated with stress and senescence responses. Eur Heart J. 2017; 38:1402–12.
Zhong Q, Huang J, Wei J, Wu R. Circular RNA CDR1as sponges miR-7-5p to enhance E2F3 stability and promote the growth of nasopharyngeal carcinoma. Cancer Cell Int. 2019; 19:252.

Due to technical limitations, full-text HTML conversion of Table 1 could not be completed. However, the table can be downloaded and accessed in the Supplementary Files.

Table1.pdf
Table1.pdf
Additionalfile1TableS1.xlsx
Additional file 1: Table S1. Patient information for all specimens.
Additionalfile1TableS1.xlsx
Additional file 1: Table S1. Patient information for all specimens.
Additionalfile2TableS2.xlsx
Additional file 2: Table S2. Sequences of primers, probes, siRNAs and peptides.
Additionalfile2TableS2.xlsx
Additional file 2: Table S2. Sequences of primers, probes, siRNAs and peptides.
Additionalfile3FigureS1.tif
Additional file 3: Figure S1. Identification of circCDKN2B-AS1, circFNACB, circEPSTI1 and circOBSL1. (A) Sequencing analysis of products amplified by the divergent primers showed the junction in circFNACB. (B) Sequencing analysis of products amplified by the divergent primers showed the junction in circEPSTI1. (C) Sequencing analysis of products amplified by the divergent primers showed the junction in circOBSL1. (D) Linear and circular isoforms of CDKN2B-AS1 were amplified from cDNA or gDNA from SiHa cells with the convergent primers and divergent primers, respectively. (E) The backsplicing point-specific probe was used in Northern blot analysis to detect endogenous circCDKN2B-AS1 in cervical cancer cell lines. Left lane: RNA molecular weight markers (2661, 1821, and 1517); right lane: circCDKN2B-AS1-specific probe.
Additionalfile3FigureS1.tif
Additional file 3: Figure S1. Identification of circCDKN2B-AS1, circFNACB, circEPSTI1 and circOBSL1. (A) Sequencing analysis of products amplified by the divergent primers showed the junction in circFNACB. (B) Sequencing analysis of products amplified by the divergent primers showed the junction in circEPSTI1. (C) Sequencing analysis of products amplified by the divergent primers showed the junction in circOBSL1. (D) Linear and circular isoforms of CDKN2B-AS1 were amplified from cDNA or gDNA from SiHa cells with the convergent primers and divergent primers, respectively. (E) The backsplicing point-specific probe was used in Northern blot analysis to detect endogenous circCDKN2B-AS1 in cervical cancer cell lines. Left lane: RNA molecular weight markers (2661, 1821, and 1517); right lane: circCDKN2B-AS1-specific probe.
Additionalfile4TableS3.docx
Additional file 4: Table S3. Sanger sequencing results of PCR products.
Additionalfile4TableS3.docx
Additional file 4: Table S3. Sanger sequencing results of PCR products.
Additionalfile5FigureS2.tif
Additional file 5: Figure S2. CircCDKN2B-AS1 knockdown suppresses EMT of cervical cancer cells. (A) E-cadherin and β-Catenin protein expression levels in SiHa (left) and CaSki (right) cells were analyzed by Western blotting after transfection with two circCDKN2B-AS1 backsplicing-specific siRNAs or a negative control siRNA. GAPDH served as a loading control.
Additionalfile5FigureS2.tif
Additional file 5: Figure S2. CircCDKN2B-AS1 knockdown suppresses EMT of cervical cancer cells. (A) E-cadherin and β-Catenin protein expression levels in SiHa (left) and CaSki (right) cells were analyzed by Western blotting after transfection with two circCDKN2B-AS1 backsplicing-specific siRNAs or a negative control siRNA. GAPDH served as a loading control.
Additionalfile6FigureS3.tif
Additional file 6: Figure S3. CircCDKN2B-AS1 overexpression facilitates the growth and vitality of cervical cancer cells. (A) The expression levels of circCDKN2B-AS1 in SiHa (left panel) and CaSki (right panel) cells with or without circCDKN2B-AS1 overexpression were determined by qRT-PCR (mean ± SEM, n = 4, unpaired Student’s t-test). (B) The growth curve of SiHa and CaSki cells transfected with circCDKN2B-AS1 overexpression plasmids or negative control plasmids was determined with CCK-8 assays. Upper panel: SiHa cells; lower panel: CaSki cells (mean ± SEM, n = 3, unpaired Student’s t-test). (C) The apoptosis level of SiHa and CaSki cells transfected with circCDKN2B-AS1 overexpression plasmids or negative control plasmids. Upper panel: SiHa cells; lower panel: CaSki cells. Right panel: quantification results of early apoptosis (mean ± SEM, n = 5, unpaired Student’s t-test). (D) Representative images (left panel) and quantification (right panel) of migration assays showing the migration capability of SiHa and CaSki cells transfected with circCDKN2B-AS1 overexpression plasmids or negative control plasmids (mean ± SEM, n = 5, unpaired Student’s t-test). (E) Representative images (left panel) and quantification (right panel) of Matrigel invasion assays showing the invasion capability of SiHa and CaSki cells transfected with circCDKN2B-AS1 overexpression plasmids or negative control plasmids (mean ± SEM, n = 5, unpaired Student’s t-test). *P < 0.05, **P < 0.01, ***P< 0.001.
Additionalfile6FigureS3.tif
Additional file 6: Figure S3. CircCDKN2B-AS1 overexpression facilitates the growth and vitality of cervical cancer cells. (A) The expression levels of circCDKN2B-AS1 in SiHa (left panel) and CaSki (right panel) cells with or without circCDKN2B-AS1 overexpression were determined by qRT-PCR (mean ± SEM, n = 4, unpaired Student’s t-test). (B) The growth curve of SiHa and CaSki cells transfected with circCDKN2B-AS1 overexpression plasmids or negative control plasmids was determined with CCK-8 assays. Upper panel: SiHa cells; lower panel: CaSki cells (mean ± SEM, n = 3, unpaired Student’s t-test). (C) The apoptosis level of SiHa and CaSki cells transfected with circCDKN2B-AS1 overexpression plasmids or negative control plasmids. Upper panel: SiHa cells; lower panel: CaSki cells. Right panel: quantification results of early apoptosis (mean ± SEM, n = 5, unpaired Student’s t-test). (D) Representative images (left panel) and quantification (right panel) of migration assays showing the migration capability of SiHa and CaSki cells transfected with circCDKN2B-AS1 overexpression plasmids or negative control plasmids (mean ± SEM, n = 5, unpaired Student’s t-test). (E) Representative images (left panel) and quantification (right panel) of Matrigel invasion assays showing the invasion capability of SiHa and CaSki cells transfected with circCDKN2B-AS1 overexpression plasmids or negative control plasmids (mean ± SEM, n = 5, unpaired Student’s t-test). *P < 0.05, **P < 0.01, ***P< 0.001.
Additionalfile7FigureS4.tif
Additional file 7: Figure S4. CircCDKN2B-AS1 cooperates with IMP3 to promote glycolysis in cervical cancer. (A) Analysis of differential bands in Figure 3B. Left upper: the partial magnification of the glue drawing in Figure 3B.The red arrow point marks the band we chose(Band1),and the blue marks the other two increased bands(Band2 and Band3); Left lower: the Gray value distribution map analyzed by using image J software; Right: quantification of the ratio of grey value between the two groups in the differential bands. (B) Silver staining indicating the differential proteins pulled down by the circCDKN2B-AS1 junction-specific probe and the oligo probe in CaSki cells. The red box indicating the differential proteins. (C) Western blot assay showing that the IMP3 protein was pulled down by biotin-labeled circCDKN2B-AS1 probes from the lysates of CaSki cells. (D) The ECAR was measured in SiHa (left upper) and CaSki (left lower) cells transfected with two junction-specific siRNAs or a negative control siRNA with Seahorse XFe assays. Right panel: quantification of glycolytic capacity and glycolysis in SiHa (right upper) and CaSki (right lower) cells (mean ± SEM, n = 6, one-way ANOVA). (E) Quantification of the %PER from the glycolysis of SiHa (left panel) and CaSki (right panel) cells transfected with negative control plasmids, circCDKN2B-AS1-overexpressing plasmids and mutant circCDKN2B-AS1-overexpressing plasmids (mean ± SEM, n = 6, one-way ANOVA). (F) The ECAR was measured in SiHa (left upper) and CaSki (left lower) cells transfected with the circCDKN2B-AS1 overexpression vector or the negative control vector with Seahorse XFe assays. Right panel: quantification of glycolytic capacity and glycolysis in SiHa (right upper) and CaSki (right lower) cells (mean ± SEM, n = 6, unpaired Student’s t-test). *P < 0.05, **P < 0.01, ***P< 0.001.
Additionalfile7FigureS4.tif
Additional file 7: Figure S4. CircCDKN2B-AS1 cooperates with IMP3 to promote glycolysis in cervical cancer. (A) Analysis of differential bands in Figure 3B. Left upper: the partial magnification of the glue drawing in Figure 3B.The red arrow point marks the band we chose(Band1),and the blue marks the other two increased bands(Band2 and Band3); Left lower: the Gray value distribution map analyzed by using image J software; Right: quantification of the ratio of grey value between the two groups in the differential bands. (B) Silver staining indicating the differential proteins pulled down by the circCDKN2B-AS1 junction-specific probe and the oligo probe in CaSki cells. The red box indicating the differential proteins. (C) Western blot assay showing that the IMP3 protein was pulled down by biotin-labeled circCDKN2B-AS1 probes from the lysates of CaSki cells. (D) The ECAR was measured in SiHa (left upper) and CaSki (left lower) cells transfected with two junction-specific siRNAs or a negative control siRNA with Seahorse XFe assays. Right panel: quantification of glycolytic capacity and glycolysis in SiHa (right upper) and CaSki (right lower) cells (mean ± SEM, n = 6, one-way ANOVA). (E) Quantification of the %PER from the glycolysis of SiHa (left panel) and CaSki (right panel) cells transfected with negative control plasmids, circCDKN2B-AS1-overexpressing plasmids and mutant circCDKN2B-AS1-overexpressing plasmids (mean ± SEM, n = 6, one-way ANOVA). (F) The ECAR was measured in SiHa (left upper) and CaSki (left lower) cells transfected with the circCDKN2B-AS1 overexpression vector or the negative control vector with Seahorse XFe assays. Right panel: quantification of glycolytic capacity and glycolysis in SiHa (right upper) and CaSki (right lower) cells (mean ± SEM, n = 6, unpaired Student’s t-test). *P < 0.05, **P < 0.01, ***P< 0.001.
Additionalfile8TableS4.xlsx
Additional file 8: Table S4. All proteins verified by mass spectrometry analysis.
Additionalfile8TableS4.xlsx
Additional file 8: Table S4. All proteins verified by mass spectrometry analysis.
Additionalfile9TableS5.xlsx
Additional file 9: Table S5. All proteins predicted by RBPmap.
Additionalfile9TableS5.xlsx
Additional file 9: Table S5. All proteins predicted by RBPmap.
Additionalfile10FigureS5.tif
Additional file 10: Figure S5. CircCDKN2B-AS1 regulates the stability of HK2 mRNA by recruiting IMP3. (A) The expression levels of HK2, ENO1, PGK1 and PKM in SiHa cells with or without circCDKN2B-AS1 knockdown were determined by qRT-PCR (mean ± SEM, n = 5, one-way ANOVA). (B) HK2 mRNA (lower) and HK2 protein (upper) levels in SiHa cells stably transfected with sh-circCDKN2B-AS1 or negative control and HK2 overexpression plasmids (mean ± SEM, n = 3, unpaired Student’s t-test). *P < 0.05, **P < 0.01, ***P< 0.001.
Additionalfile10FigureS5.tif
Additional file 10: Figure S5. CircCDKN2B-AS1 regulates the stability of HK2 mRNA by recruiting IMP3. (A) The expression levels of HK2, ENO1, PGK1 and PKM in SiHa cells with or without circCDKN2B-AS1 knockdown were determined by qRT-PCR (mean ± SEM, n = 5, one-way ANOVA). (B) HK2 mRNA (lower) and HK2 protein (upper) levels in SiHa cells stably transfected with sh-circCDKN2B-AS1 or negative control and HK2 overexpression plasmids (mean ± SEM, n = 3, unpaired Student’s t-test). *P < 0.05, **P < 0.01, ***P< 0.001.
Additionalfile11FigureS6.jpg
Additional file 11: Figure S6. The expression levels of circCDKN2B-AS1 and IMP3 protein in in transplanted tumor tissues. The levels of IMP3 protein expression (left) in transplanted tumor tissues (mean ± SEM, n = 6, unpaired Student’s t-test). Right: the association between circCDKN2B-AS1 and IMP3 protein expressions (n=12, Pearson’s correlation coefficient).
Additionalfile11FigureS6.jpg
Additional file 11: Figure S6. The expression levels of circCDKN2B-AS1 and IMP3 protein in in transplanted tumor tissues. The levels of IMP3 protein expression (left) in transplanted tumor tissues (mean ± SEM, n = 6, unpaired Student’s t-test). Right: the association between circCDKN2B-AS1 and IMP3 protein expressions (n=12, Pearson’s correlation coefficient).

Download PDF

Journal Publication

published 11 Dec, 2020

Read the published version in Journal of Experimental & Clinical Cancer Research →

Review #2 received at journal
30 Nov, 2020
Editorial decision: Accept
30 Nov, 2020
Review #1 received at journal
28 Nov, 2020
Reviewer #2 agreed at journal
22 Nov, 2020
Reviewer #1 agreed at journal
20 Nov, 2020
Reviewers invited by journal
18 Nov, 2020
Editor assigned by journal
17 Nov, 2020
Submission checks completed at journal
17 Nov, 2020
Editor invited by journal
17 Nov, 2020

You are reading this latest preprint version

CircCDKN2B-AS1 interacts with IMP3 to stabilize Hexokinase 2 mRNA and facilitate cervical squamous cellcarcinoma aerobic glycolysis progression

Status:

Journal Publication

Version 2

Abstract

Figures

Background

Methods

Results

Discussion

Conclusions

List Of Abbreviations

Declarations

References

Table

Supplementary Files

Status:

Journal Publication

Version 2