- Open Access
A novel protein encoded by circFNDC3B inhibits tumor progression and EMT through regulating Snail in colon cancer
Molecular Cancer volume 19, Article number: 71 (2020)
Colon cancer (CC) is a common malignant cancer. Recently, circFNDC3B was found to exert biological function in multiple cancers. However, it was unclear whether the potential protein encoded by circFNDC3B is involved in carcinogenesis of CC.
We used Sanger sequence and RNase R digestion assay to confirm the existence of circFNDC3B, and quantitative real-time PCR was used to evaluate the circRNA’s expression. Then fluorescence in situ hybridization (FISH) was performed to study location of circFNDC3B. The identification of protein encoded by circFNDC3B was performed using LC-MS/MS. The function of circFNDC3B-218aa on proliferation, invasion and migration were assessed by CCK8 assays, colony formation assays, transwell assays, wound-healing assays and animal experiments. RNA-sequencing and western blot were used to identify the gene regulated by circFNDC3B-218aa. Finally, glucose metabolism-related assays were performed to further investigate function of circFNDC3B-218aa.
CircFNDC3B was localized mostly in the cytoplasm, and was decreased in CC cell lines and tissues. The patients with low circFNDC3B expression had a shorter OS (P = 0.0014) than patients with high expression. Moreover, circFNDC3B inhibited the proliferation, invasion and migration of CC cells. Next, we identified that circFNDC3B could encode a novel protein circFNDC3B-218aa. Furthermore, circFNDC3B-218aa, not circFNDC3B, inhibited the proliferation, invasion and migration of CC. Additionally, the in vivo experiments implied that up-regulated circFNDC3B-218aa exhibited an inhibitory effect on CC progression. By RNA-sequencing, western blot and glucose metabolism-related assays, we found that circFNDC3B-218aa inhibited the expression of Snail, and subsequently promoted the tumor-suppressive effect of FBP1 in CC.
The novel circFNDC3B-218aa may serve as a tumor suppressive factor and potential biomarker which may supply the potential therapeutic target for CC.
Colon cancer (CC) is a malignant cancer with increasingly high prevalence, and is the fourth leading cause of cancer-related death worldwide . Although the clinical outcome of CC patients has met great improvement due to the development of novel therapeutic approaches, CC is still a deadly disease . Hence, it is important to deeply explore the molecular mechanism of CC progression and to identify the potential prognosis biomarkers of CC.
Circular RNAs (circRNAs) are believed to be a group of endogenous non-coding RNAs (ncRNAs) [3, 4]. CircRNAs are circularized by the connection of a 5′ splice site with the 3′ splice site of an upstream exon or intron, which is also called back-splicing reaction . Increasing evidences indicated that circRNAs are not useless byproducts created as splicing errors and that they may have regulatory roles intracellularly . Acting as miRNA sponges is the most commonly reported function of circRNAs [7, 8]. Additionally, circRNAs could bind to related proteins and regulate their biological behaviors [9, 10]. Like most ncRNAs, circRNAs were previously assumed to be untranslatable because of the lack of obvious open reading frames (ORFs). However, several recent studies have reported that circRNAs could actually encode proteins [11, 12].
CircFNDC3B, originated from exons 5 and 6 of FNDC3B gene in chr3, has been recently found to be functional in previous studies [13,14,15]. Among these studies, Liu et al. reported that circFNDC3B is a tumor-suppressive factor in bladder cancer by sponging miR-1178-3p . Another study indicated that circFNDC3B interacts with IGF2BP3 and promotes the migration and invasion of gastric cancer . Although the miRNA sponging and protein binding effects of it have received some attention, there are few reports about the protein translation process of circFNDC3B. In order to explore its potential protein encoding ability, we have matched circFNDC3B with circRNADb (http://126.96.36.199:8000/circrnadb/circRNADb.php) and found a predicted ORF. However, it remains unclear whether the predicted protein encoded by circFNDC3B exists in CC and whether it exerts function in CC.
In this study, we identified that circFNDC3B was low expressed in CC. CircFNDC3B may encode a 218 amino acid novel protein which was termed the circFNDC3B-218aa. The in vitro and in vivo data indicated that circFNDC3B-218aa decreased proliferation, invasion and migration abilities of CC cells. Moreover, we found that circFNDC3B-218aa inhibited tumor progression and epithelial-mesenchymal transition (EMT) via alleviating the repressive effect of Snail on FBP1 in colon cancer.
Materials and methods
Ethical statement and tissue collection
Eighty-seven CC specimens and their paired normal tissues were obtained from patients at the Lanzhou University Second Hospital between 2013 and 2015. The tissues were immediately frozen and stored at − 80 °C after surgical resection. Paraffin-fixation was performed by experienced pathologists. The collection of all clinical data was performed after each surgery. This study was approved by the Ethics Committee of Lanzhou University Second Hospital and the written informed consent of each participant was obtained before surgery.
Cell culture and RNase R treatment
Human CC cell lines (DLD1, HCT116, SW480, LoVo, Caco2 and HT29) and normal colon mucosal epithelial cell line (NCM460) were purchased from ATCC. All cells were cultured in DMEM supplemented with 10% FBS under standard culture conditions (5% CO2, 37 °C). RNase R (Epicentre Technologies, USA) was used to degrade linear mRNA. In brief, we extracted RNAs from CC cells and split RNA to two parts: one for RNase R digestion and another for control with digestion buffer only. The samples were incubated at 37 °C for 30 min. Then the expression levels of circFNDC3B and FNDC3B were detected by quantitative real-time RT-PCR.
Total RNA extraction and quantitative real-time PCR
The extraction of total RNA sample was completed with RNAiso Plus (TaKaRa, Japan) according to manufacturer’s instructions. The concentration and purity of all samples were subsequently measured via NanoDrop 2000 (Thermo Scientific, Wilmington, DE, USA). Corresponding cDNAs were generated using PrimeScript RT Master Mix (TaKaRa, Japan). Quantitative real-time PCRs (qRT-PCR) were performed on a LightCycler® 96 System (Roche, Switzerland). The relative expression was calculated by the ΔCq method. All primer sequences were given in Additional files 1 (Table S1).
RNA fluorescence in situ hybridization
RNA fluorescence in situ hybridization (RNA-FISH) was performed to study the location of circFNDC3B by a Fluorescent in Situ Hybridization Kit (RiboBio, Guangzhou, China) in accordance with the manufacturer’s instructions.
circFNDC3B overexpressing vector and the control plasmid (pCDH-CMV-MCS-EF1-copGFP-Puro) were purchased from General Biosystems (General Biosystems, Anhui, China). The expression of circFNDC3B was knockdown by siRNAs generated by GenePharma (GenePharma Corporation, Shanghai, China). To verify the exist of circFNDC3B-218aa, we constructed four flag labeled vectors for circFNDC3B using lentivirus vectors (pCDH-CMV-MCS-EF1-copGFP-Puro) which were purchased from General Biosystems (General Biosystems, Anhui, China) and transfected them stably into cell lines. The full designs of sequences synthesized and cloned into vector were listed in Additional file 2. Cells were seeded in approximately 60% confluence before transfection. All siRNAs and their vectors were transfected using a lipofectamine 3000 transfection kit (Invitrogen, USA) following the manufacturer’s instructions. The transfection efficiency was assessed via qRT-PCR analysis.
Cell proliferation assay
Cell Counting Kit-8 (CCK-8, ImmunoWay Biotechnology Company Plano, TX, USA) were performed to monitor the growth of cells. In brief, cells transfected with si-circFNDC3B or overexpression vectors were seeded into 96-well plates and cultured for 24, 48, 72 and 96 h, respectively. OD450 values were determined.
Colony formation assay
The cloning capability of CC cells was evaluated by colony formation assays. Cells were seeded into 6-well plates and incubated for 2 weeks. After staining, the visible colonies were stained and counted.
Transwell invasion assay
Transwell assays were used to assess the invasion ability of cells using the transwell chambers (Costar, USA), which were precoated with Matrigel. CC cells in serum-free medium were added to the upper chambers (pore size, 8 μm; Corning Inc., Tewksbury, MA, USA) and DMEM with 10% FBS was added to the lower chambers. The cells migrated into the lower chambers were fixed in 4% paraformaldehyde after incubation for 24 h and stained with crystal violet. Random fields were imaged and number of cells was counted.
Wound-healing assays were used to assess the migration ability of cells. The CC cells were cultured in the six-well plates until the cell confluence reached 95%. After we performed vertical scratched in the six-well plates, serum-free medium was added in it. Afterwards, image collection and migration distance measurement were conducted at 0 h and 24 h.
Western blot analysis
Extraction of proteins was performed with RIPA buffer (CWBIO, China). The samples were electrophoresed by SDS-PAGE and transferred to nitrocellulose membranes, which were then incubated with primary antibodies at 4 °C overnight. After that, the membranes were incubated with secondary antibodies at room temperature for 1 h. Signals detection and images acquisition were performed by Immobilon ECL substrate (Millipore, Germany) and Optimax X-ray Film Processor (Protec, Germany). The antibodies used in this study were listed in Additional file 3.
Analysis of peptide patterns by LC-MS/MS
The identification of protein was performed using LC-MS/MS. In brief, proteins were subjected to 12% SDS-PAGE gel. The protein bands near 25 kDa were excised and chopped into 1mm3 pieces. Peptide mixtures were extracted from the chopped gels for LC-MS/MS analysis finally.
Intracellular ATP production was measured using a CellTiter-Glo Luminescent Cell Viability Assay kit (Promega, Madison, WI, USA). Lactate, pyruvate production and intracellular glucose uptake were measured using lactate, pyruvate and glucose assay kits (BioVision, Milpitas, CA, USA) according to the manufacturers’ instructions, respectively.
Ethical approval was approved by the Ethics Committee of Sun Yat-sen University (Guangzhou, China). All animal care and procedures were performed according to the institutional guidelines. Lentiviral vector (pCDH-CMV-MCS-EF1-copGFP-Puro) mixed with packing plasmid (psPAX2) and envelope plasmid (pMD2G) were transfected into 293 T cells using Lipofectamine 3000. Viral supernatants were collected after 48 h. Lentivirus was infected into HCT116 cells according to the manufacturer’s instructions. The stable lines (transfected respectively with empty vector, circFNDC3B, circFNDC3B-mut and circFNDC3B-218aa) were selected with Puromycin (Life Technologies). The stable transfected HCT116 cells (5 × 106/200 μL PBS) were injected subcutaneously into the right backs of 4-week-old BALB/c nude mice. The size of tumors was measured every 3 days and their volumes were calculated. Thirty days later, the mice were sacrificed, and the tumors were excised for further analysis. To assess the effect of circFNDC3B-218aa on liver metastases, HCT116 cells transfected with four vectors were injected into the spleens of recipient mice, and bioluminescent images were captured. The liver metastatic nodules were excised for immunohistochemistry (IHC).
In brief, 4-μm thick paraffin-embedded sections were prepared and received dewaxing, dehydrating and treating with 3% H2O2 at 37 °C for 10 min. After rinsed twice with PBS, the samples were boiled in ethylene diamine tetraacetic acid (EDTA, pH 8.0) at 95 °C for 20 min and cooled down to 25 °C for antigen repairment. Subsequently, the sections were treated with normal goat serum at room temperature for 10 min to block non-specific antigens and incubated with primary antibodies at 4 °C overnight, followed by incubation with the secondary antibody at 37 °C for 30 min. Finally, the sections were observed under a light microscope (Leica, Germany) after treated with diaminobenzidine. Five different fields (magnification × 400) in each section were selected to evaluate the count of positive cells.
The data are showed as mean ± standard deviation (SD) or as values directly. As appropriate, unpaired two-tailed Student’s t-test was used to evaluate the significance of difference. Statistical analyses were conducted using Stata (version 13.1, StataCorp,College Station, TX) and GraphPad Prism 7 software (GraphPad Software, San Diego, CA). P value < 0.05 was considered to be statistically significant.
The characterization of circFNDC3B in CC
Firstly, we designed two sets of primers. Convergent primers were used to amplify only the linear form of FNDC3B, the mRNA, while divergent primers were to amplify only the circular form, the circFNDC3B. Using cDNA and genomic DNA (gDNA) as templates, the PCR products were validated by electrophoresis. The results revealed that the single and distinct product around the expected size was amplified in the divergent primers for cDNA, but no product for gDNA (Fig. 1a). Furthermore head-to-tail splicing of the amplified circFNDC3B was confirmed by the Sanger sequencing. These results suggested that circFNDC3B was derived from exon 5 and exon 6 of the FNDC3B gene, which was consistent with the circFNDC3B data from circBase (Fig. 1b). Compared with the linear form, circFNDC3B was resistant to RNase R digestion (Fig. 1c), with a longer half-life (Fig. 1d). Further fluorescence in situ hybridization (FISH) showed that circFNDC3B was localized mostly in the cytoplasm of CC cells (Fig. 1e).
CircFNDC3B is downregulated in CC and related to poor survival
Next, we examined the circFNDC3B expression in different CC cell lines (DLD1, HCT116, SW480, LoVo, Caco2 and HT29) using qRT-PCR analysis. CircFNDC3B was significantly lower in CC cell lines in comparison to a normal human colon mucosal epithelial cell line, the NCM460 (Fig. 1f). In addition, we examined the circFNDC3B expression in 87 pairs of CC tissues and paired normal tissues and found that circFNDC3B expression was significantly decreased in cancer tissues compared with their normal counterparts (Fig. 1g). Additionally, circFNDC3B expression was low in 77% of CC patients (67/87) (Fig. 1h). The relationship between circFNDC3B and clinicopathological features was then analyzed. The CC patients were classified into high- and low-level groups according to circFNDC3B expression (Table 1). Briefly, low circFNDC3B was positively correlated with lymphatic metastasis (P = 0.024). No significant correlation was observed concerning other clinicopathological features, including age, T stage, gender, and differentiation. Kaplan–Meier survival curves indicated that patients with lower circFNDC3B expression had a significantly shorter OS than patients with higher expression (P = 0.0014) (Fig. 1i). Furthermore, multivariate COX analysis revealed that the expression of circFNDC3B was an independent prognostic factor for CC prognosis (low expression vs high expression, HR = 2.27; 95% CI 1.19–4.32; P = 0.013) (Table 2).
CircFNDC3B inhibits the proliferation, migration and invasion of CC cells
To investigate the biological function of circFNDC3B in CC, gain-of-function and loss-of-function assays were performed. After transfection of the circFNDC3B overexpression plasmid, there was no significant change of FNDC3B mRNA in CC cells (Fig. 2a). Otherwise, the circFNDC3B was knocked down after the cells being transfected with siRNA targeting the back-splice region, and there was no significant change in its linear counterpart, FNDC3B mRNA (Fig. 2b). CircFNDC3B upregulation significantly inhibited the proliferation of CC cells (Fig. 2c), as confirmed in the CCK8 assays. Colony formation assays also indicated that the clonality of CC cells was suppressed by up-regulating circFNDC3B (Fig. 2e). The result of transwell assays showed that invasion of CC cells was suppressed, again by circFNDC3B upregulation (Fig. 2g). Wound-healing assay showed that circFNDC3B overexpression also inhibited the migration abilities of CC cells (Fig. 2i). In contrast, downregulation of circFNDC3B yielded the reverse effects, which promoted the proliferation, invasion and migration abilities in these CC cells (Fig. 2d, f, h, and j).
A novel protein encoded by circFNDC3B
Upon the annotation of circFNDC3B from circRNA Db, an ORF named the 657-nt ORF with the potential to encode a 218-aa peptide was found. In human circFNDC3B, the tandem “AUG” within the RNA circle could start the translation of a novel protein. The ORF suggested that it took more than one whole circle of circFNDC3B to translate the putative 218 amino-acid (aa) protein. The protein contained a unique sequence and we termed it the ‘circFNDC3B-218aa’ (Fig.3a). An internal ribosomal entrance site (IRES) is required for the 5′-cap independent translation, and we found there was an IRES in the ORF (from + 448 to + 524). Furthermore, we verified the activity of the IRES in circFNDC3B by using a dual-luciferase assay, which was showed in Fig. 3b.
To verify the exist of circFNDC3B-218aa, four kinds of flag-labeled vectors were constructed (Lv-vector; Lv-circFNDC3B-flag: flag-labeled circFNDC3B sequence was cloned into a CMV-induced expression vector; Lv-circFNDC3B-flag-mut: flag-labeled circFNDC3B sequence with start codon mutant (ATG → ACG) was cloned into a CMV-induced expression vector. Lv-circFNDC3B-218aa: the flag-labeled circFNDC3B-218aa sequence was cloned into a CMV-induced expression vector) (Fig. 3c). After transfection with Lv-circFNDC3B-flag vector and Lv-circFNDC3B-flag-mut vector, the circFNDC3B expression was significantly up-regulated. In the CC cells transfected with Lv-vector and Lv-circFNDC3B-218aa, there was no obvious change in the level of circFNDC3B (Fig. 3d). All four vectors transfected did not affect the linear FNDC3B mRNA. Furthermore, the flag-labeled circFNDC3B-218aa protein was detected at approximate 25 kDa after Lv-circFNDC3B-flag or Lv-circFNDC3B-218aa transfection. Besides, the FNDC3B antibodies targeting the middle part of FNDC3B could detect circFNDC3B-218aa in Lv-vectors and Lv-circFNDC3B-flag-mut-transfected cells, while Lv-circFNDC3B-flag or Lv-circFNDC3B-218aa transfected cells showed more obvious expression of circFNDC3B-218aa at 25 kDa. Moreover, we found that these four vectors did not affect the 133Kda FNDC3B expression (Fig. 3e). After separation of the pulldown cell lysates from the cell transfected with Lv-circFNDC3B-flag, protein bands at the size of 25 kDa were excised and submitted for the LC-MS/MS (Fig. 3f). As a result, the circFNDC3B-218aa translated from circFNDC3B was confirmed, and the unique sequences of circFND3CB-218aa (MIQTCKVIEDSTGVRR) were subsequently identified (Fig. 3g). To further prove the circFNDC3B-218aa is endogenous, we used the FNDC3B antibodies which can also recognize the circFNDC3B-218aa with their same amino acids sequence to conduct co-immunoprecipitation in the extracts from the cells without forced expression. The LC-MS/MS successfully identified the specific peptide fragments from circFNDC3B-218aa (MAMGRAIVVEVAEAAAVVVPELRK) (Additional file 4 Fig.S1 a and b). Taken together, circFNDC3B could encode a novel protein in CC.
CircFNDC3B-218aa, not circFNDC3B inhibits the proliferation, invasion and migration abilities of CC in vitro and in vivo
To explore the biological function of circFNDC3B-218aa, the four vectors were transfected into CC cell lines. The results of CCK8 and colony formation assays indicated that the proliferation abilities of CC cells were inhibited by transfection with Lv-circFNDC3B or Lv-circFNDC3B-218aa. However, there was no effect on the proliferation in the CC cells transfected with Lv-circFNDC3B-mut which could not be translated into the circFNDC3B-218aa. The results suggested that only overexpression of circFNDC3B-218aa decreased the proliferation abilities (Fig. 4a and b). In transwell assays and wound-healing assays, overexpression of circFNDC3B-218aa inhibited the invasion and migration abilities of CC cells. However, mutant circFNDC3B did not influence the invasion and migration abilities of CC cells (Fig. 4c and d). Collectively, the above results showed that the protein form (circFNDC3B-218aa), not the circle RNA form (circFNDC3B), inhibits the proliferation, invasion and migration abilities of CC. To further assess whether circFNDC3B-218aa exerts tumor-suppressing effect in vivo, a xenograft mouse model was established by subcutaneous injection of HCT116 cells (n = 5 for each group). The results indicated that the growth rates and tumor weight were significantly lower in the Lv-circFNDC3B and Lv-circFNDC3B-218aa group than in the vector control and the Lv-circFNDC3B-mut group (Fig. 5a). Additionally, in vivo metastatic assay showed that Lv-circFNDC3B and Lv-circFNDC3B-218aa decreased detectable liver metastasis and the incidence of liver metastasis compared with the Lv-vector group and Lv-circFNDC3B-mut group (Fig. 5b).
Effect of CircFND3CB-218aa on CC cells and the underlying mechanism
To explore the potential regulatory mechanisms influenced by circFNDC3B-218aa in CC cell lines, we performed RNA-Seq analyses. Differentially expressed genes (DEGs) in Lv-vector and Lv-circFNDC3B-218aa groups were analyzed (Additional file 5), and Fig. 6a listed the top 10 significant DEGs. We further evaluated the gene expression through qRT-PCR, and found that four of these genes were indeed differentially expressed (Fig. 6b). Among these genes, there was few studies reported that OAS1 and IFI44 were related to tumor progression [16, 17]. Additionally, Snail was downregulated and FBP1 was upregulated in circFNDC3B-218aa group. Further qRT-PCR and western blot results showed that the Snail level was decreased and the FBP1 level was increased when the CC cells were transfected with Lv-circFNDC3B or Lv-circFNDC3B-218aa (Fig. 6c and d). According to reported articles, Snail could induce EMT via inhibiting the expression of FBP1 [18, 19]. Therefore, we speculated that circFNDC3B-218aa exerts a tumor-related role in Snail/FBP1/EMT axis. We assessed the expression of Snail, FBP1, E-cadherin and vimentin in the liver metastases specimens of nude mice. The IHC results showed that the expression levels of FBP1 and E-cadherin were significantly increased, while Snail and vimentin were decreased in the Lv-circFNDC3B and Lv-circFNDC3B-218aa group compared with the control and the Lv-circFNDC3B-mut group (Fig. 6e). Such results further revealed that circFNDC3B-218aa exhibited inhibitory effects on CC EMT and regulated the expression of Snail and FBP1.
Furthermore, we investigated the roles of circFNDC3B-218aa, FBP1 and Snail in EMT. Further Western blot analysis indicated that compared with the control group, the expression of E-cadherin and FBP1 were significantly decreased, whereas that of vimentin was increased in the HCT116 cells with overexpressed Snail. However, the above phenomenon was restored after overexpressing circFNDC3B-218aa in upregulated Snail HCT116 (Fig. 7a and b). Moreover, in transwell invasion assays and wound-healing assays, we found that the invasion and migration abilities of CC cells enhanced by Snail could be alleviated by circFNDC3B-218aa (Fig. 7c and e), suggesting the overexpression of circFNDC3B-218aa restores the oncogenic roles of Snail in CC cells. Given that FBP1 could inhibit EMT and circFNDC3B-218aa could regulate FBP1, it suggested that circFNDC3B-218aa inhibited EMT by enhancing FBP1. To investigate whether circFNDC3B-218aa could inhibit cancer progression by enhancing FBP1, we transfected CC cells with siFBP1 and circFNDC3B-218aa overexpression plasmid. Analogously, the transwell assays and wound-healing assays showed that cell invasion and migration ability were promoted by FBP1 downregulation, but were partially restored by transfection of circFNDC3B-218aa (Fig. 7d and f). These results illustrated that circFNDC3B-218aa reversed the effect of downregulated FBP1.
As a gluconeogenesis regulatory enzyme, FBP1 had been further demonstrated play an important role in impairing aggressive phenotype in various cancers via a metabolic switch from glycolysis to oxidative phosphorylation (OXPHOS) [20, 21]. Therefore, we investigated whether circFNDC3B-218aa exhibited an inhibitory effect on EMT through FBP1 participating in Warburg effect inhibition. The glucose uptake, pyruvate production and lactate production were remarkably increased in HCT116 and DLD1 cell lines with upregulated Snail or downregulated FBP1, which could be reversed by circFNDC3B-218aa overexpression (Fig. 7g-i). In addition, circFNDC3B-218aa could also elevate ATP produced by OXPHOS upon Snail upregulation or FBP1 silencing in CC cells (Fig. 7j). These results indicated that circFNDC3B-218aa enhanced metabolic reprogramming from glycolysis to oxidative phosphorylation via inhibiting Snail-FBP1 signaling axis, subsequently suppressed the EMT progression in CC cells. Taken together, these data suggested that circFNDC3B-218aa could inhibit the cancer progression and EMT via alleviating the repressive effect of Snail on FBP1 (Fig. 8).
CircRNAs are novel RNA molecules, with various biological functions. Among them, acting as a miRNA sponge and interaction with related proteins are the two common function patterns of circRNA. A previous study has reported that artificial circRNA was translatable in eukaryotic cells . Recently, it is also reported that circRNAs translation could be driven by N6-methyladenosine . Meanwhile, the circZNF609 and the circMbl have been found to be protein-encoding in two other studies [23, 24]. These reported evidences all suggested that endogenous circRNAs have the potential to encode protein. Moreover, several protein-encoding circRNAs have been identified to be involved in cancer progression. For example, a functional protein encoded by circFBXW7 could repress glioma tumorigenesis . In the present study, we found that there is an ORF in the circFNDC3B sequence. Further experiments showed that circFNDC3B could encode a novel protein which suppressed colon cancer progression and EMT in vitro and in vivo.
According to previous studies, there are two main reasons why it is difficult to detect the circRNA-translated protein. One is that the protein encoded by circRNAs may overlap with the linear mRNA of the translated products. The other is that the linear mRNA products obscures the protein translated by circRNA because the circRNA level was low compared with their linear counterparts . With the development of testing approaches, growing number of peptides encoded by circRNAs have been detected, and their potential biological functions have also been identified. CircRNA encoded proteins usually show the functional independence with their host gene products. For example, it has been reported that circPPP1R12A-73aa promotes CC progression via Hippo-YAP signaling . In this study, we found that circFNDC3B-218aa, rather than circFNDC3B, regulated the expression of Snail and FBP1. It suggested that circFNDC3B-218aa exerted biological function in CC independently, which is similar to previous studies. In contrast, some products encoded by circRNAs are involved in their host gene functions. For example, circSHPPH-146aa prevents its host gene SHPPR from degradation  and circFBXW7-185aa induced c-Myc degradation . Because circFNDC3B-218aa shares the majority of its amino acids sequence with that of FNDC3B with the exception of the unique amino acids in the ORF, it is possible that circFNDC3B-218aa involved in the functions of FNDC3B. However, the underlying mechanism is needed further validation.
In present study, we found that circFNDC3B-218aa could suppress the Snail expression and subsequently upregulate FBP1 and inhibit EMT in CC. Therefore, we speculated that circFNDC3B-218aa exerted its tumor-suppressive effect in Snail/FBP1/EMT axis. FBP1, one of the related-limiting enzymes in gluconeogenesis, plays an important role in glucose metabolism . It had been widely reported that metabolic shifts control EMT progression and trigger tumor malignancy, especially an enhanced glycolytic phenotype switched from OXPHOS, which was called Warburg effect [28,29,30]. This specific tumor hallmark not only provides cancer cells with nutrients, but also constructs a more acidic tumor environment that leads to extracellular matrix destruction and induces metastasis [31, 32]. In our study, we found that circFNDC3B-218aa suppressed tumor EMT progression was mainly dependent on inhibiting of Snail expression, subsequently enhancing FBP1-induced OXPHOS. Consistently, Dong et al. demonstrated that loss of FBP1 by Snail-mediated repression was required for EMT induction via glycolytic shift . Thus, disruption of glycolysis process might be a potential therapeutic target for circFNDC3B/circFNDC3B-218aa deficient colon cancer.
Taken together, we found that the expression of circFNDC3B was decreased in CC tissues and cells. Additionally, the results showed that circFNDC3B could encode circFNDC3B-218aa. The circFNDC3B-218aa, not circFNDC3B itself, inhibited the proliferation and invasion of CC cells in vitro and in vivo. Furthermore, we speculated that circFNDC3B-218aa could inhibit the cancer progression and EMT via alleviating the repressive effect of Snail on FBP1. Collectively, this study expanded upon the understanding of the coding potential of circRNAs. The novel circFNDC3B-218aa is a potential biomarker which may supply the potential therapeutic target for CC.
Availability of data and materials
All data in our study are availability upon request.
Open reading frame
Fluorescence in situ hybridization
Internal ribosomal entrance site
Cassidy S, Syed BA. Colorectal cancer drugs market. Nat Rev Drug Discov. 2017;16(8):525–6.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69(1):7–34.
Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, et al. Circular RNA: a new star of noncoding RNAs. Cancer Lett. 2015;365(2):141–8.
Kulcheski FR, Christoff AP, Margis R. Circular RNAs are miRNA sponges and can be used as a new class of biomarker. J Biotechnol. 2016;238:42–51.
Wang F, Nazarali AJ, Ji S. Circular RNAs as potential biomarkers for cancer diagnosis and therapy. Am J Cancer Res. 2016;6(6):1167–76.
Glazar P, Papavasileiou P, Rajewsky N. circBase: a database for circular RNAs. RNA. 2014;20(11):1666–70.
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.
Zeng K, Chen X, Xu M, Liu X, Hu X, Xu T, et al. CircHIPK3 promotes colorectal cancer growth and metastasis by sponging miR-7. Cell Death Dis. 2018;9(4):417.
Teng L, Chen Y, Chen H, He X, Wang J, Peng Y, et al. Circular RNA hsa_circ_0021001 in peripheral blood: a potential novel biomarker in the screening of intracranial aneurysm. Oncotarget. 2017;8(63):107125–33.
Du WW, Yang W, Liu E, Yang Z, Dhaliwal P, Yang BB. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 2016;44(6):2846–58.
Schneider T, Bindereif A. Circular RNAs: coding or noncoding? Cell Res. 2017;27(6):724–5.
Yang Y, Fan X, Mao M, Song X, Wu P, Zhang Y, et al. Extensive translation of circular RNAs driven by N(6)-methyladenosine. Cell Res. 2017;27(5):626–41.
Liu H, Bi J, Dong W, Yang M, Shi J, Jiang N, et al. Invasion-related circular RNA circFNDC3B inhibits bladder cancer progression through the miR-1178-3p/G3BP2/SRC/FAK axis. Mol Cancer. 2018;17(1):161.
Chen T, Yu Q, Shao S, Guo L. Circular RNA circFNDC3B protects renal carcinoma by miR-99a downregulation. J Cell Physiol. 2019.
Hong Y, Qin H, Li Y, Zhang Y, Zhuang X, Liu L, et al. FNDC3B circular RNA promotes the migration and invasion of gastric cancer cells via the regulation of E-cadherin and CD44 expression. J Cell Physiol. 2019;234(11):19895–910.
Liu X, Xing H, Gao W, Yu D, Zhao Y, Shi X, et al. A functional variant in the OAS1 gene is associated with Sjogren's syndrome complicated with HBV infection. Sci Rep. 2017;7(1):17571.
Power D, Santoso N, Dieringer M, Yu J, Huang H, Simpson S, et al. IFI44 suppresses HIV-1 LTR promoter activity and facilitates its latency. Virology. 2015;481:142–50.
Dong C, Yuan T, Wu Y, Wang Y, Fan TW, Miriyala S, et al. Loss of FBP1 by snail-mediated repression provides metabolic advantages in basal-like breast cancer. Cancer Cell. 2013;23(3):316–31.
Liu GM, Li Q, Zhang PF, Shen SL, Xie WX, Chen B, et al. Restoration of FBP1 suppressed snail-induced epithelial to mesenchymal transition in hepatocellular carcinoma. Cell Death Dis. 2018;9(11):1132.
Dai Q, Li N, Zhou X. Increased miR-21a provides metabolic advantages through suppression of FBP1 expression in non-small cell lung cancer cells. Am J Cancer Res. 2017;7(11):2121–30.
Qattan AT, Radulovic M, Crawford M, Godovac-Zimmermann J. Spatial distribution of cellular function: the partitioning of proteins between mitochondria and the nucleus in MCF7 breast cancer cells. J Proteome Res. 2012;11(12):6080–101.
Chen CY, Sarnow P. Initiation of protein synthesis by the eukaryotic translational apparatus on circular RNAs. Science. 1995;268(5209):415–7.
Pamudurti NR, Bartok O, Jens M, Ashwal-Fluss R, Stottmeister C, Ruhe L, et al. Translation of CircRNAs. Mol Cell. 2017;66(1):9–21 e7.
Legnini I, Di Timoteo G, Rossi F, Morlando M, Briganti F, Sthandier O, et al. Circ-ZNF609 is a circular RNA that can be translated and functions in Myogenesis. Mol Cell. 2017;66(1):22–37 e9.
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;110(3).
Zhang M, Huang N, Yang X, Luo J, Yan S, Xiao F, et al. A novel protein encoded by the circular form of the SHPRH gene suppresses glioma tumorigenesis. Oncogene. 2018;37(13):1805–14.
Zheng X, Chen L, Zhou Y, Wang Q, Zheng Z, Xu B, et al. A novel protein encoded by a circular RNA circPPP1R12A promotes tumor pathogenesis and metastasis of colon cancer via hippo-YAP signaling. Mol Cancer. 2019;18(1):47.
Halldorsson S, Rohatgi N, Magnusdottir M, Choudhary KS, Gudjonsson T, Knutsen E, et al. Metabolic re-wiring of isogenic breast epithelial cell lines following epithelial to mesenchymal transition. Cancer Lett. 2017;396:117–29.
Lunetti P, Di Giacomo M, Vergara D, De Domenico S, Maffia M, Zara V, et al. Metabolic reprogramming in breast cancer results in distinct mitochondrial bioenergetics between luminal and basal subtypes. FEBS J. 2019;286(4):688–709.
Peppicelli S, Toti A, Giannoni E, Bianchini F, Margheri F, Del Rosso M, et al. Metformin is also effective on lactic acidosis-exposed melanoma cells switched to oxidative phosphorylation. Cell Cycle. 2016;15(14):1908–18.
Liu X, Wang X, Zhang J, Lam EK, Shin VY, Cheng AS, et al. Warburg effect revisited: an epigenetic link between glycolysis and gastric carcinogenesis. Oncogene. 2010;29(3):442–50.
Qiu Z, Guo W, Wang Q, Chen Z, Huang S, Zhao F, et al. MicroRNA-124 reduces the pentose phosphate pathway and proliferation by targeting PRPS1 and RPIA mRNAs in human colorectal cancer cells. Gastroenterology. 2015;149(6):1587–98 e11.
This work was supported by the National Natural Science Foundation of China (81802402, 81900597), the Medical Scientific Research Foundation of Guangdong Province of China (A201742), the Science and Technology Planning Project of Zhuhai city of China (20171009E030080), the PhD Start-up Fund of Natural Science Foundation of Guangdong Province 2018A030310323, the China Postdoctoral Science Foundation (2019 M653197, 2019 M653199).
Ethics approval and consent to participate
All patients consented to an institutional review broad-approved protocol that allows comprehensive analysis of tissues samples (Ethics committee of Lanzhou University Second Hospital).
Consent for publication
Written consents for publication were obtained from all the patients involved in our study.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
The primer sequences used in this study.
Additional file 2.
The full design of vectors was listed.
Additional file 3: Table S2.
Antibodies used in this study
Additional file 4: Figure S1. a
The lysates from indicated cells were separated by SDS-PAGE. Protein at 25 kDa were excised manually and summited for identification by LC-MS/MS. b circFNDC3B-218aa specific peptide (MAMGRAIVVE VAEAAAVVVPELRK) was identified.
Additional file 5.
DEGs were screened.
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.
About this article
Cite this article
Pan, Z., Cai, J., Lin, J. et al. A novel protein encoded by circFNDC3B inhibits tumor progression and EMT through regulating Snail in colon cancer. Mol Cancer 19, 71 (2020). https://doi.org/10.1186/s12943-020-01179-5
- Colon cancer