Hsa_circ_0026416 promotes proliferation and migration in colorectal cancer via miR-346/NFIB axis
Cancer Cell International volume 20, Article number: 494 (2020)
Colorectal cancer (CRC) is one of the most common cancers worldwide. Circular RNAs (circRNAs), a novel class of non-coding RNAs, have been confirmed to be key regulators of many diseases. With many scholars devoted to studying the biological function and mechanism of circRNAs, their mysterious veil is gradually being revealed. In our research, we explored a new circRNA, hsa_circ_0026416, which was identified as upregulated in CRC with the largest fold change (logFC = 3.70) of the evaluated circRNAs via analysing expression profiling data by high throughput sequencing of members of the GEO dataset (GSE77661) to explore the molecular mechanisms of CRC.
qRT-PCR and western blot analysis were utilized to assess the expression of hsa_circ_0026416, miR-346 and Nuclear Factor I/B (NFIB). CCK-8 and transwell assays were utilized to examine cell proliferation, migration and invasion in vitro, respectively. A luciferase reporter assay was used to verify the combination of hsa_circ_0026416, miR-346 and NFIB. A nude mouse xenograft model was also utilized to determine the role of hsa_circ_0026416 in CRC cell growth in vivo.
Hsa_circ_0026416 was markedly upregulated in CRC patient tissues and plasma and was a poor prognosis in CRC patients. In addition, the area under the curve (AUC) of hsa_circ_0026416 (0.767) was greater than the AUC of CEA (0.670), CA19-9 (0.592) and CA72-4 (0.575). Functionally, hsa_circ_0026416 promotes cell proliferation, migration and invasion both in vitro and in vivo. Mechanistically, hsa_circ_0026416 may function as a ceRNA via competitively absorbing miR-346 to upregulate the expression of NFIB.
In summary, our findings demonstrate that hsa_circ_0026416 is an oncogene in CRC. Hsa_circ_0026416 promotes the progression of CRC via the miR-346/NFIB axis and may represent a potential biomarker for diagnosis and therapy in CRC.
Colorectal cancer (CRC) is one of the most common cancers worldwide, with an incidence of 1.8 million cases and 896,000 deaths in 2017 , ranking third in cancer incidence and second in cancer-related death for both sexes combined . The 5-year relative survival rate of CRC patients is 65% but increases to 90% when colorectal cancers are diagnoses at a localized stage . Although excellent progress has been made in elucidating the diagnosis, treatment and molecular mechanisms of CRC and other solid tumours in recent years, with therapies including surgery, chemoradiotherapy [4, 5], molecular therapy [6, 7] and immunotherapy [8,9,10,11], the poor prognosis of colorectal patients is still closely related to recurrence and metastasis [12, 13]. Blood-based markers in current clinical use, such as cancer antigen (CA) 19–9, cancer antigen (CA) 724 and carcinoembryonic antigen (CEA), are used to supervise and monitor response to therapy but have unsatisfactory sensitivity and specificity, ranging from 40 to 70% and 73 to 90%, respectively, making them insufficient as screening or diagnostic markers [14,15,16]. Therefore, there is an urgent need to identify an ideal biomarker to improve the rate of early diagnosis in colorectal cancer.
CircRNAs are novel stars of the non-coding RNA world that are abundant in mammalian cells, featuring relative stability and high tissue and cell specific expression [17,18,19]. In contrast to linear RNAs, circRNAs have covalently closed loop structures lacking 5′–3′ polarity and a polyadenylated tail . Until now, benefiting from advances in RNA deep sequencing technology and bioinformatics, a large number of circRNAs has been detected in various cell lines and species . Recent evidence suggests that circRNAs may play a crucial role in the occurrence and development of a variety of diseases, such as atherosclerosis, neurological disorders, prion diseases and cancers, by working as miRNA sponges, RNA-binding proteins or transcriptional regulators [22, 23]. Some circRNAs that are deregulated in CRC tissues or cell lines, such as hsa_circ_0000069, hsa_circ_0014717, circular BANP and ciRS-7, have been implicated in proliferation, metastasis and poor prognosis of CRC, indicating their crucial roles in this malignancy [24,25,26,27]. Nevertheless, the molecular mechanism of the involvement of circRNAs in CRC remains to be explored.
In our research, we conducted a comprehensive exploration of a new circRNA (hsa_circ_0026416). We found that hsa_circ_0026416 significantly upregulated in CRC tissues and high expression of hsa_circ_0026416 leads to poor prognosis. In addition, hsa_circ_0026416 acts as a ceRNA to regulate NFIB via competition for miR-346, facilitating CRC proliferation and metastasis. In summary, this study extends our knowledge on the biological roles and mechanisms of hsa_circ_0026416 in CRC and provides a novel biomarker for the diagnosis and treatment of CRC.
Materials and methods
Analysing circRNA expression profiling by high throughput sequencing
A circRNA expression profile (GSE77661) was downloaded from GEO. Analysis of colorectal-related differentially expressed circRNAs was performed using the Bioconductor Limma package in R software. The criteria for selection of differentially expressed circRNAs (DEcircRNAs) were P-value < 0.05 and |log2FC|> 1.
Patients and clinical specimens
From January 2017 to November 2018, 169 pairs of CRC tissues and adjacent normal mucosa were collected from Qilu Hospital of Shandong University, China. At the same time, 212 cases of preoperative plasma samples and 64 cases of postoperative plasma samples were collected. All patients met the following inclusion criteria: received enterectomy and were confirmed by pathological diagnosis; complete clinicopathological data, including gender, age, tumour location, tumour diameter, tumour differentiations, pT stage, pN stage, distant metastasis, pTNM, lymphovascular invasion, and perineural invasion; complete follow-up information; and written informed consent. Patients were excluded if they received radiotherapy or chemotherapy before surgery, presented with other malignant disease within the past 5 years, were lost to follow-up, or exhibited incomplete clinicopathological data. After surgery, all patients were followed-up for at least 2 years, and patients with advanced CRC received chemotherapy and/or radiotherapy according to the National Comprehensive Cancer Network (NCCN) Guidelines. In addition, fresh normal plasma samples were obtained from 183 healthy people at Qilu Hospital of Shandong University in February 2017. All tissue and plasma samples were quickly frozen in liquid nitrogen after removal and stored at −80 °C until further experiments. This study was approved by the Ethics Committee of Qilu Hospital. Each study participant provided informed consent.
The normal intestinal epithelial cell line (HCO), HEK297T and four human CRC cell lines (DLD-1, HCT-8, HCT-116, SW480) were purchased from the Culture Collection of Chinese Academy of Sciences (Shanghai, China). DLD-1 and HCT-8 cell lines were maintained in RPMI-1640 media (Hyclone, Logan, Utah, USA) with 10% foetal bovine serum (FBS, Clark Bioscience, California, USA). HCO, HEK293T, HCT-116 and SW480 cell lines were maintained in DMEM (Gibco, Grand Island, NY, USA) supplemented with 10% FBS (Clark Bioscience, California, USA) and cultured in a 37 °C incubator containing 5% CO2. After 2–3 stable generations, cells were used for subsequent experiments and analyses. Short tandem repeat (STR) DNA fingerprinting was utilized to identify all cell lines. Mycoplasma contamination was checked at least once a month using the MycoAlert Mycoplasma Detection Kit (Lonza, Walkersville, MD). No mycoplasma was detected in any of the cell lines.
RNA extraction, reverse transcription PCR and qRT-PCR
Total RNA was extracted using TRIzol Reagent (Thermo Fisher Scientific, MA, USA), and the concentration and purity of total RNA were determined by an ultraviolet spectrophotometer (Eppendorf, Hamburg, Germany). Nuclear and cytoplasmic RNA isolation was performed using a PARIS Kit (Invitrogen, Carlsbad, MA, USA) following the manufacturer’s instructions. PrimeScript™ RT Reagent Kit with gDNA Eraser (Takara, Tokyo, Japan) was used to reverse transcribe RNA to cDNA. Real-time PCR was implemented using SYBR Premix Ex Taq (Takara, Tokyo, Japan) on a LightCycler® 2.0 Real-time PCR System (Roche, Basel, Switzerland). GAPDH was used as an internal reference gene. The results were analysed by the 2-ΔΔCt method. All primer sequences are listed in Additional file 1: Table S1.
Plasmid construction and cell transfection
To overexpress hsa_circ_0026416, hsa_circ_0026416 cDNA (704 bp) was inserted into the PLCDH-ciR vector (BioSune, Jinan, China). Amplified from the cDNA of HCT-8 cells, full length NFIB was cloned into the pcDNA3.1 ( +) vector, full length of hsa_circ_0026416 and 3′UTR of NFIB was inserted into the pmirGLO vector. PmirGLO-mut (hsa_circ_0026416), pmirGLO-mut (NFIB 3′UTR), all miRNA mimics, the miR-346 inhibitor, NFIB si-RNA and 2′-O-Me-Modified-si-circ_0026416 were synthesized by Genepharma (Shanghai, China). Hsa_circ_0026416 si-RNA was obtained from Suzhou Ribo Life Science Co., Ltd (Suzhou, China).
For transient transfection of plasmids, HCT-8, SW480 and HEK293T cells were transfected with Lipofectamine 3000 reagent (Invitrogen, Waltham, MA, USA). Si-RNAs, miRNA mimics and the miRNA inhibitor were transfected into cells using the X-tremeGENE transfection reagent (Roche Applied Science, Indianapolis, IN, USA). All cell transfections were based on the manufacturer’s instructions. The sequences of these nucleic acids are listed in Additional file 1: Table S1.
Matrigel (BD Biosciences, New Jersey, USA) was added into the Transwell chamber, and the 24-well plates containing the chamber were placed into the cell culture incubator to accelerate the solidification of Matrigel. Twenty-four hours after transfection, 1 × 105 (loss of function) HCT-8 cells and 2 × 105 (loss of function) SW480 cells were suspended in 200 µl serum-free medium and placed into the upper chambers of a Transwell (Costar, NY, USA). When we performed gain of function and rescue experiments, 7 × 104 HCT-8 cells and 1.3 × 105 SW480 cells were seeded into the upper chambers of each Transwell. Then, we added 600 µl media containing 20% FBS into the lower well of each chamber. After incubation for 48 h, cells in the upper cavity were removed with cotton swabs, and cells that had invaded the lower surface were fixed, stained using crystal violet and imaged under a microscope at 100 × magnification (Olympus, Tokyo, Japan). When we assessed cell migration, all steps were the same as above, except that we did not add Matrigel.
For CCK-8 assay, 3000 HCT-8 or SW480 cells per well were cultured in 96-well plates, and each treatment featured three repeats. At 0, 24, 48, 72 and 96 h after transfection, 10 μl CCK-8 reagent (Dojindo Molecular Technology Company, Japan) was added to each well with the cell culture medium, and samples were placed back into the incubator for another 2 h. Absorbance was measured at 450 nm by a microplate reader (Eppendorf, Germany).
Flow cytometry analysis
Cell cycle distribution was analysed by flow cytometry. Cells were harvested, washed twice with PBS and fixed in 75% ethanol at −20 °C for one hour. RNA was removed by incubating cells with RNase A (100 μg/ml, Sigma-Aldrich) at 37 °C for 30 min. Cells were then stained with propidium iodide (PI) solution (50 μg/ml, Sigma-Aldrich) at 4 °C for 30 min and analysed by flow cytometry (Beckman Coulter Gallios, California, USA).
Luciferase reporter assay
The HEK293T cells were seeded in 24-well plates overnight and then co-transfected with pmirGLO-circ_0026416 WT or MUT and miR-346 mimics or NC. For HCT-8 and SW480 cells, they were co-transfected with pmirGLO-NFIB 3′UTR WT or MUT and miR-346 mimics or NC. After a 48-h culture, Firefly and Renilla luciferase activities were measured with the Dual-Luciferase Reporter Assay System (Promega, WI, USA) according to the manufacturer’s instructions, and Firefly luciferase activity was standardized by Renilla luciferase activity.
To assess protein expression, both the immunoreactive percentage and intensity were scored. The percentage of positive cells was graded as 0, < 5%; + 1, 5–25%; + 2, 26–50%; + 3, 51–75%; and + 4, 76–100% positive cells, and the intensity of cellular staining was scored as: 0, negative; + 1, weak; + 2, moderate; and + 3, strong. The final staining score was obtained by multiplying the two scores.
Protein extraction and western blot
Protein extraction and western blot were performed as previously described . Antibodies were as follows: NFIB (Abcam, ab186738, Cambridge, UK, 1:1000 dilution) and GAPDH (CST, #5174, MA, USA, 1:3000 dilution).
In vivo tumorigenesis
BALB/c athymic male nude mice (4 weeks old) were obtained from Charles River Biotechnology (Beijing, China). SW480 cell suspensions (2 × 106 cells in 0.1 ml PBS) were subcutaneously inoculated into the right axillary fossa of mice (n = 6 per group). 2′-O-Me-Modified-si-circ_0026416 was injected into the tumours every 3 days after tumours formed. Tumour growth was measured every 3 days for 22 days, at which time the mice were sacrificed. Tumour volume (V) was calculated as follows: V = (length diameter) × (width diameter)2 /2. The longest diameter did not exceed 2.0 cm, and the general condition of all mice was well throughout the experiment. All animal experiments were approved by the Animal Care and Use Committee of Shandong University, as well as the animal welfare settings .
All data were analysed using GraphPad Prism 8.0 (La Jolla, USA). The efficiency of hsa_circ_0026416 in the diagnosis of CRC was evaluated by receiver operating characteristic (ROC) curve. The two groups were compared using two-tailed Student’s t-test. The relationship between the expression of hsa_circ_0026416 and clinicopathological characteristics was assessed by Chi-squared test. Univariate and multivariate Cox proportional hazard regression models were performed to analyse the effects of different clinicopathological factors on overall survival. Kaplan–Meier survival curve and log-rank test were used to determine the OS rate of CRC patients with different hsa_circ_0026416 expression levels. Data are displayed as means ± SD of three independent experiments. P-values < 0.05 were defined as statistically significant.
Characterization of hsa_circ_0026416 in CRC tissues
To explore the roles of circRNAs in the progression of colorectal cancer, we analysed colorectal tissue-related data in the GSE77661 datasets using RNA sequencing analysis. We detected 67 differentially expressed circRNAs (|logFC|> 1 and P < 0.05). Finally, we selected the upregulated circRNA (hsa_circ_0026416) with the largest fold change (logFC = 3.70) for additional research (Fig. 1a). Then, by querying the circBase database (https://circbase.org/), we determined that hsa_circ_0026416 is derived from regions of exons 3, 4, 5, 6, 7 and 8 of the KRT6C gene (GenBank: NM_173086) (Fig. 1b). Next, the splice junction was directly verified by Sanger sequencing (Fig. 1c). Afterwards, we detected hsa_circ_0026416 expression in CRC tissues and their corresponding adjacent normal mucosa (n = 169) by qRT-PCR. By analysing the data, we found that expression of hsa_circ_0026416 was higher in CRC tissues than in adjacent normal mucosal tissues (P < 0.001) (Fig. 1d). We then divided patients into two cohorts based on the median expression value of hsa_circ_0026416 (hsa_circ_0026416 high- and low-expression), and Kaplan–Meier analysis revealed that CRC patients with higher hsa_circ_0026416 expression exhibited reduced overall survival compared to CRC patients with lower hsa_circ_0026416 expression (Fig. 1e). Next, we observed that expression of hsa_circ_0026416 was higher in cancer cells of the colorectum than in normal colorectal cells (Fig. 1f). With the highest expression observed in HCT-8 and SW480, we chose these two cell lines for subsequent study. qRT-PCR analysis of nuclear and cytoplasmic RNAs indicated that hsa_circ_0026416 was primarily localized in the cytoplasm in HCT-8 and SW480 cells (Figs. 1g, h). These results revealed that hsa_circ_0026416 is an upregulated cytoplasmic circRNA and may play a vital role in CRC tumour progression.
Hsa_circ_0026416 is a diagnostic biomarker of CRC
First, we compared hsa_circ_0026416 expression in the plasma of 212 patients with CRC before surgery with that of hsa_circ_0026416 in 63 patients after surgery (Fig. 2a, b). Meanwhile, we compared hsa_circ_0026416 expression in CRC patients’ preoperative plasma (n = 212) with healthy individuals’ plasma (n = 183) (Fig. 2c). These results demonstrated that expression of hsa_circ_0026416 in postoperative plasma was distinctly lower than in preoperative plasma (P < 0.01). Furthermore, compared to hsa_circ_0026416 expression in healthy individuals’ plasma, its expression in CRC patients was significantly increased (P < 0.001). Next, ROC curve analysis was used to assess the diagnostic ability of hsa_circ_0026416 (Fig. 2d). The AUC (area under the curve) of hsa_circ_0026416, CEA, CA19-9 and CA72-4 was 0.767, 0.670 0.592 and 0.575, respectively. To assess the correlation between hsa_circ_0026416 and clinicopathological factors, we divided patients into two cohorts according to their median expression value of hsa_circ_0026416 as a cut-off threshold (hsa_circ_0026416 high- and low-expression). We discovered that levels of hsa_circ_0026416 were significantly associated with tumour differentiation (P = 0.029), distant metastasis (P < 0.001), pTNM stage (P = 0.004), lymphovascular invasion (P = 0.002) and perineural invasion (P = 0.004) (Table 1). Finally, univariate and multivariate Cox regression analysis both revealed that advanced T stage, positive distant metastasis, terminal pTNM stage and high hsa_circ_0026416 expression were all indicative of poor prognosis (Table 2). These results suggested that hsa_circ_0026416 is an independent prognostic factor in CRC patients (P < 0.05) and has value as a biomarker for the diagnosis of CRC.
Hsa_circ_0026416 promotes CRC cell migration, invasion and proliferation in vitro
To assess the effects of hsa_circ_0026416 on colorectal cancer cell behaviour, hsa_circ_0026416 was knocked down and overexpressed by transfection with two small interfering RNAs (si-RNA1 and si-RNA2) targeting hsa_circ_0026416 and an hsa_circ_0026416 overexpression plasmid, respectively (Fig. 3a, b). Transwell assays showed that knockdown of hsa_circ_0026416 markedly inhibited the migratory and invasive capacities of CRC cells (Fig. 3c, d), while overexpression of hsa_circ_0026416 dramatically increased these abilities (Fig. 3e, f). CCK-8 assays revealed that knockdown of hsa_circ_0026416 inhibited proliferation in CRC cells (Fig. 3g, i), whereas overexpression of hsa_circ_0026416 promoted proliferation ability (Fig. 3h, j). Nevertheless, cell cycle experiments indicated that hsa_circ_0026416 had no effect on cell cycle distributions of CRC cells (Fig. 3k–v). These results prove that hsa_circ_0026416 dramatically promotes CRC cell migration, invasion and proliferation but has no effect on cell cycle distributions of CRC cells in vitro.
Hsa_circ_0026416 acts as a ceRNA for miR-346 in CRC
Next, taking into consideration that hsa_circ_0026416 was present primarily in the cytoplasm and many circRNAs have been proven to act as miRNA sponges, we next explored the capability of hsa_circ_0026416 to bind miRNAs. RegRNA2.0 and CircInteractome were used to predict miRNA response elements (MREs) harboured by hsa_circ_0026416. Among all predicted miRNAs, we selected the top 12 as candidate miRNAs (Fig. 4a). Among them, miR-346, miR-1184, and miR-1205 overlapped between the two databases (Fig. 4b).
MiR-346, the most downregulated by the ratio of Firefly and Renilla after transfecting pmirGLO-circ_0026416 and miRNA mimics, was chose for further analysis (Fig. 4c). Then, we utilized The Cancer Genome Atlas (TCGA) sequencing data to verify miR-346 expression in human CRC tissues (n = 601) compared to normal colorectal tissues (n = 11) (Fig. 4d) and performed qRT-PCR in 24 pairs of CRC and normal colorectal tissues to further prove the consequence of TCGA. The results both indicated that the expression of miR-346 in CRC tissues is downregulated compared to normal colorectal tissues (Fig. 4e). We further transfected si-circ_0026416 or a negative control into HCT-8 and SW480 cells using qRT-PCR to detect expression of miR-346. qRT-PCR indicated that knock down of hsa_circ_0026416 obviously increased miR-346 expression in HCT-8 and SW480 cells (Fig. 4f). Next, luciferase reporter assays were employed to ascertain whether miR-346 directly targets hsa_circ_0026416. We constructed two pmirGLO vectors, one containing the full length hsa_circ_0026416 (wild type, WT) and the other one with the miR-346-binding site mutated (MUT) (Fig. 4g). A conspicuous decrease in luciferase reporter activity was observed in HEK293T cells co-transfected with miR-346 mimics and hsa_circ_0026416 (WT), but not with the MUT vector (Fig. 4h). These results demonstrate that hsa_circ_0026416 acts as a sponge for miR-346.
MiRNAs repress transcription or induce target mRNA degradation by binding to the 3′UTR of target mRNAs . To identify potential target genes of miR-346, four miRNA target gene databases (miRDB, miRWalk, starBase v3.0, TargetScan 7.1) were utilized to predict miR-346 direct target genes, and 32 genes appeared in all four databases (Fig. 4i). Among them, NFIB, which was predicted with a high score and is upregulated in human colorectal cancer as a well-known oncogene , was confirmed as the most likely gene directly targeted by miR-346. Next, we used The Human Protein Atlas (THPA, https://www.proteinatlas.org) proteomic data to verify the expression of NFIB in human CRC tissue compared to normal colorectal tissue (Fig. 4j). Immunohistochemistry results of tumour and normal tissues were scored as described in the Materials and Methods, and results indicated that expression of NFIB in CRC tissue was higher than in normal tissue (Fig. 4k). Next, using qRT-PCR to detect the expression of NFIB in 24 pairs of CRC and normal colorectal tissues, we obtained the same consequence (Fig. 4l). Meanwhile, Yangyang Li et al. reported that miR-346 directly targets NFIB to inhibit the growth of glioma cells. In order to verify whether miR-346 could directly target NFIB in CRC cells, we constructed two pmirGLO vectors, one containing the 3′UTR of NFIB (wild type, WT) and the other one with the miR-346-binding site mutated (MUT) (Fig. 4n). A dramatic decrease in luciferase reporter activity was observed in both HCT-8 and SW480 cells co-transfected with miR-346 mimics and NFIB 3′UTR (WT), but not with the MUT vector (Fig. 4o). Next, using qRT-PCR, we discovered that expression of NFIB was decreased in response to hsa_circ_0026416 knockdown and was reduced by miR-346 mimics as well (Fig. 4m). Meanwhile, western blot showed that expression of NFIB was increased when hsa_circ_0026416 was overexpressed, whereas knockdown of hsa_circ_0026416 inhibited NFIB expression (Fig. 4p–s). Consequently, these findings verified the hsa_circ_0026416/miR-346/NFIB ceRNA network in CRC.
Hsa_circ_0026416 facilitates CRC progression through a ceRNA network in vitro.
To assess whether hsa_circ_0026416 exerts biological functions in CRC through a ceRNA network, rescue experiments were performed. First, miR-346 and NFIB were overexpressed by transfection with miR-346 mimics and NFIB overexpression plasmids. Meanwhile, miR-346 and NFIB were knocked down by transfection with a miR-346 inhibitor and a small interfering RNA (si-NFIB) targeting NFIB. qRT-PCR revealed that NFIB and miR-346 were upregulated in HCT-8 and SW480 cells in response to pcDNA3.1-NFIB and miR-346 mimics co-transfection, while si-NFIB and miR-346 inhibitor co-transfection downregulated expression of NFIB and miR-346 in HCT-8 and SW480 cells (Fig. 5a–d). As anticipated, the decreased migratory, invasive and proliferative abilities induced by miR-346 mimics were rescued by co-transfecting pcDNA3.1-NFIB in HCT-8 and SW480 cells (Fig. 5e–h). These results are consistent with the above findings that the increased migratory, invasive and proliferative capacities induced by miR-346 inhibitor were regained through co-transfecting si-NFIB in HCT-8 and SW480 cells (Fig. 5i–l). These results indicated that hsa_circ_0026416 promotes CRC progression in vitro through the hsa_circ_0026416/miR-346/NFIB ceRNA network.
Downregulation of hsa_circ_0026416 inhibits CRC proliferation in vivo
To further explore the tumour-promoting effect of hsa_circ_0026416 in vivo, we purchased stable small interfering RNA (2′-O-Me-Modified-si-circ_0026416) for intratumoral injection. Compared to the NS (normal saline) group, we discovered that mean tumour volume and weight in the 2′-O-Me-Modified-si-circ_0026416 group was markedly reduced (Fig. 6a–c), proving that targeted inhibition of hsa_circ_0026416 effectively inhibits tumour proliferation in vivo.
Although proteins are the basis of life-sustaining activities, more than 80% of the human genome is actively transcribed into RNA, while less than 3% of the human genome encodes translated proteins [33,34,35]. RNA molecules that are not translated into proteins are called non-coding RNAs (ncRNAs). In recent years, massive studies have indicated that ncRNAs play a crucial role in cell biology and cancer. At present, circRNAs, members of the ncRNA family, are considered to be rising stars in the study of tumorigenesis and development.
Circular RNA, which exists in the cytoplasm of eukaryotic cells, was first reported in 1979 . Benefiting from the rapid advancements in high-throughput RNA sequencing, increasing circRNAs are being discovered and identified as playing vital roles in tumour progression. CircRNAs have distinct functions in different cellular locations. Cytoplasmic circRNAs have a common recognition mechanism that occurs through base pairing of the complementary region between circRNAs and their target RNA sequence.
This constitutes an intricate RNA-RNA network. For example, (i) circRNAs can interact with RNA binding proteins, such as circ-Foxo3, of which there are many. By interacting with CDK2 and p21, circRNA inhibits the cell cycle and block the transition from G1 to S phase . (ii) CircRNAs can modulate the stability of mRNAs, such as circ-RasGEF1B, which strengthens the stability of ICAM-1 mRNA in mouse macrophages  and so on. Among them, a classical pattern known as ceRNA drew our attention. In this co-regulation pattern, circRNAs act as miRNA sponges to competitively absorb miRNAs and protect target mRNAs from suppression . The first observation of a circRNA acting as a miRNA sponge was CDR1as, which has more than 70 conserved binding sites for miRNA-7 .
In our research, hsa_circ_0026416 was identified as an upregulated circRNA with the largest fold change (logFC = 3.70) via analysing expression profiling by high throughput sequencing from the GEO dataset (GSE77661). Hsa_circ_0026416 was upregulated in CRC tissues, plasma, and HCT-8 and SW480 cells, conveying a poor prognosis in CRC patients. In addition, our data revealed compared to CEA, CA19-9 and CA724 levels, hsa_circ_0026416 levels in plasma may represent a more promising diagnostic marker for CRC. Hsa_circ_0026416 promotes CRC cell proliferation, migration and invasion both in vitro and in vivo by competitively absorbing miR-346 and upregulating its target gene, NFIB. As previous studies have shown, the majority of circRNAs act as miRNA sponges or competitive endogenous RNAs (ceRNAs) to regulate the expression of target genes . MiRNAs can combine to the matched 3′-UTRs of mRNAs, and circRNAs also include miRNA target sites. By competing for miRNAs, circRNAs indirectly regulate the translation of target mRNAs. However, circRNAs have a variety of functions, including interacting with RNA binding proteins to regulate tumour progression. This implies that hsa_circ_0026416 may also participate in additional regulatory circuits apart from the hsa_circ_0026416/miR-346/NFIB axis. This is a limitation of our research, which needs further investigation.
MiR-346, which functions as a tumour suppressor , was contained in both databases and was the highest-ranked miRNA. Luciferase reporter assay identified that miR-346 binds with hsa_circ_0026416. The Cancer Genome Atlas (TCGA) sequencing data indicated that miR-346 was downregulated in CRC tissues and qRT-PCR experiment in 24 pairs of CRC and normal colorectal tissues further validated the consequence of TCGA. Then, we utilized four online bioinformatics prediction databases to predict potential target genes of miR-346. NFIB appeared in all four databases and has been reported as a target gene of miR-346 in glioma . Moreover, NFIB was identified as an oncogene in CRC [31, 43]. THPA database and qRT-PCR experiment both demonstrated that NFIB was upregulated in CRC. Afterwards, luciferase reporter assay further identified that miR-346 binds with NFIB. In addition, western blot indicated that knocking down hsa_circ_0026416 expression downregulated the expression of NFIB, and increasing the expression of hsa_circ_0026416 upregulated NFIB expression. These results further proved that hsa_circ_0026416 acts as a competing endogenous RNA (ceRNA) for miR-346 and upregulates expression of NFIB, further promoting proliferation and migration in CRC (Fig. 6d).
Our data indicate that hsa_circ_0026416 can be considered a potential biomarker for the diagnosis of CRC. Its expression was tremendously upregulated in CRC tissues and plasma, and its AUC value of 0.767 was higher than that of CEA (0.670), CA19-9 (0.592) and CA72-4 (0.575). In vivo experiments further revealed that targeted inhibition of hsa_circ_0026416 slowed tumour growth. Taken together, these factors all indicate that hsa_circ_0026416 is an excellent biomarker for the diagnosis of CRC and a novel target for the treatment of CRC.
In summary, our study demonstrated that hsa_circ_0026416 is upregulated in CRC tissues and CRC patient plasma. Univariate and multivariate Cox regression analysis both revealed that.
high hsa_circ_0026416 expression was the indicative of poor prognosis. Functionally, hsa_circ_0026416 promotes CRC cell proliferation, migration and invasion both in vitro and in vivo. Mechanistically, hsa_circ_0026416 acts as a ceRNA for miR-346 to upregulate NFIB and promote tumour progression (Fig. 6d). Our results indicate that hsa_circ_0026416 plays a crucial role in CRC progression and may represent an important and promising biomarker for CRC patient diagnosis and targeted therapy.
Our research proves that hsa_circ_0026416 is upregulated in CRC patient tissues and plasma and high expression of hsa_circ_0026416 lead to poor prognosis in CRC patients. Hsa_circ_0026416 promotes CRC cell proliferation, migration and invasion via the miR-346/NFIB axis in vitro and in vivo. Therefore, hsa_circ_0026416 may represent a potential biomarker for diagnosis and therapeutic targeting.
Availability of data and materials
All data generated or analysed during this study are included in the published article and are available in Gene Expression Omnibus (GEO), circBase, The Human Protein Atlas, RegRNA2.0, CircInteractome, miRDB, miRWalk, starBase v3.0 and TargetScan 7.1. (https://www.ncbi.nlm.nih.gov/geo/, https://circbase.org/, https://www.proteinatlas.org, https://regrna2.mbc.nctu.edu.tw/, https://circinteractome.nia.nih.gov/, https://mirdb.org/, https://mirwalk.umm.uni-heidelberg.de/, https://starbase.sysu.edu.cn/, https://www.targetscan.org/vert_71/).
Adjacent normal mucosa
Fetal bovine serum
Small interfering RNA
Cell counting kit-8
Receiver operating characteristic curve
Area under curve
The Cancer Genomic Atlas
The Human Protein Atlas
Fitzmaurice C, Abate D, Abbasi N, Abbastabar H, Abd-Allah F, Abdel-Rahman O, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2017: a systematic analysis for the Global Burden of Disease Study. JAMA Oncol. 2019;5(12):1749–68.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Miller KD, Siegel RL, Lin CC, Mariotto AB, Kramer JL, Rowland JH, et al. Cancer treatment and survivorship statistics, 2016. CA Cancer J Clin. 2016;66(4):271–89.
Mohammad N, Malvi P, Meena AS, Singh SV, Chaube B, Vannuruswamy G, et al. Cholesterol depletion by methyl-β-cyclodextrin augments tamoxifen induced cell death by enhancing its uptake in melanoma. Mol Cancer. 2014;13:204.
Mohammad N, Singh SV, Malvi P, Chaube B, Athavale D, Vanuopadath M, et al. Strategy to enhance efficacy of doxorubicin in solid tumor cells by methyl-β-cyclodextrin: Involvement of p53 and Fas receptor ligand complex. Sci Rep. 2015;5:11853.
Singh S, Chouhan S, Mohammad N, Bhat MK. Resistin causes G1 arrest in colon cancer cells through upregulation of SOCS3. FEBS Lett. 2017;591(10):1371–82.
Muhammad N, Steele R, Isbell TS, Philips N, Ray RB. Bitter melon extract inhibits breast cancer growth in preclinical model by inducing autophagic cell death. Oncotarget. 2017;8(39):66226–36.
Bhattacharya S, Muhammad N, Steele R, Kornbluth J, Ray RB. Bitter melon enhances natural killer-mediated toxicity against head and neck cancer cells. Cancer Prev Res. 2017;10(6):337–44.
Bhattacharya S, Muhammad N, Steele R, Peng G, Ray RB. Immunomodulatory role of bitter melon extract in inhibition of head and neck squamous cell carcinoma growth. Oncotarget. 2016;7(22):33202–9.
Singh SV, Ajay AK, Mohammad N, Malvi P, Chaube B, Meena AS, et al. Proteasomal inhibition sensitizes cervical cancer cells to mitomycin C-induced bystander effect: the role of tumor microenvironment. Cell Death Dis. 2015;6(10):e1934.
Ganesh K, Stadler ZK, Cercek A, Mendelsohn RB, Shia J, Segal NH, et al. Immunotherapy in colorectal cancer: rationale, challenges and potential. Nat Rev Gastroenterol Hepatol. 2019;16(6):361–75.
Lech G, Slotwinski R, Slodkowski M, Krasnodebski IW. Colorectal cancer tumour markers and biomarkers: recent therapeutic advances. World J Gastroenterol. 2016;22(5):1745–55.
Loree JM, Kopetz S. Recent developments in the treatment of metastatic colorectal cancer. Ther Adv Med Oncol. 2017;9(8):551–64.
Das V, Kalita J, Pal M. Predictive and prognostic biomarkers in colorectal cancer: a systematic review of recent advances and challenges. Biomed Pharmacother. 2017;87:8–19.
Sun J, Fei F, Zhang M, Li Y, Zhang X, Zhu S, et al. The role of (m)SEPT9 in screening, diagnosis, and recurrence monitoring of colorectal cancer. BMC Cancer. 2019;19(1):450.
Young GP, Pedersen SK, Mansfield S, Murray DH, Baker RT, Rabbitt P, et al. A cross-sectional study comparing a blood test for methylated BCAT1 and IKZF1 tumor-derived DNA with CEA for detection of recurrent colorectal cancer. Cancer Med. 2016;5(10):2763–72.
Yuan X, Yuan Y, He Z, Li D, Zeng B, Ni Q, et al. The regulatory functions of circular RNAs in digestive system cancers. Cancers. 2020;12(3):770.
Liu J, Li D, Luo H, Zhu X. Circular RNAs: the star molecules in cancer. Mol Aspects Med. 2019;70:141–52.
Barrett SP, Salzman J. Circular RNAs: analysis, expression and potential functions. Development. 2016;143(11):1838–47.
Zhang P, Zhang XO, Jiang T, Cai L, Huang X, Liu Q, et al. Comprehensive identification of alternative back-splicing in human tissue transcriptomes. Nucleic Acids Res. 2020;48(4):1779–89.
Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.
Qian P, Linbo L, Xiaomei Z, Hui P. Circ_0002770, acting as a competitive endogenous RNA, promotes proliferation and invasion by targeting miR-331-3p in melanoma. Cell Death Dis. 2020;11(4):264.
Chen L, Kong R, Wu C, Wang S, Liu Z, Liu S, et al. Circ-MALAT1 functions as both an mRNA translation brake and a microRNA sponge to promote self-renewal of hepatocellular cancer stem cells. Adv Sci. 2020;7(4):1900949.
Guo JN, Li J, Zhu CL, Feng WT, Shao JX, Wan L, et al. Comprehensive profile of differentially expressed circular RNAs reveals that hsa_circ_0000069 is upregulated and promotes cell proliferation, migration, and invasion in colorectal cancer. Onco Targets Ther. 2016;9:7451–8.
Wang F, Wang J, Cao X, Xu L, Chen L. Hsa_circ_0014717 is downregulated in colorectal cancer and inhibits tumor growth by promoting p16 expression. Biomed Pharmacother. 2018;98:775–82.
Weng W, Wei Q, Toden S, Yoshida K, Nagasaka T, Fujiwara T, et al. Circular RNA ciRS-7-A promising prognostic biomarker and a potential therapeutic target in colorectal cancer. Clin Cancer Res. 2017;23(14):3918–28.
Zhu M, Xu Y, Chen Y, Yan F. Circular BANP, an upregulated circular RNA that modulates cell proliferation in colorectal cancer. Biomed Pharmacother. 2017;88:138–44.
Ma F, Li X, Ren J, Guo R, Li Y, Liu J, et al. Downregulation of eukaryotic translation initiation factor 3b inhibited proliferation and metastasis of gastric cancer. Cell Death Dis. 2019;10(9):623.
Hansen BC, Gografe S, Pritt S, Jen KC, McWhirter CA, Barman SM, et al. Ensuring due process in the IACUC and animal welfare setting: considerations in developing noncompliance policies and procedures for institutional animal care and use committees and institutional officials. Faseb j. 2017;31(10):4216–25.
Muhammad N, Bhattacharya S, Steele R, Ray RB. Anti-miR-203 suppresses ER-positive breast cancer growth and stemness by targeting SOCS3. Oncotarget. 2016;7(36):58595–605.
Liu Z, Chen J, Yuan W, Ruan H, Shu Y, Ji J, et al. Nuclear factor I/B promotes colorectal cancer cell proliferation, epithelial-mesenchymal transition and 5-fluorouracil resistance. Cancer Sci. 2019;110(1):86–988.
Li Y, Xu J, Zhang J, Zhang J, Zhang J, Lu X. MicroRNA-346 inhibits the growth of glioma by directly targeting NFIB. Cancer Cell Int. 2019;19:294.
Chen LL. The biogenesis and emerging roles of circular RNAs. Nat Rev Mol Cell Biol. 2016;17(4):205–11.
Consortium, E.P. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489(7414):57–74.
Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, et al. Landscape of transcription in human cells. Nature. 2012;489(7414):101–8.
Hsu MT, Coca-Prados M. Electron microscopic evidence for the circular form of RNA in the cytoplasm of eukaryotic cells. Nature. 1979;280(5720):339–40.
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.
Ng WL, Marinov GK, Liau ES, Lam YL, Lim YY, Ea CK. Inducible RasGEF1B circular RNA is a positive regulator of ICAM-1 in the TLR4/LPS pathway. RNA Biol. 2016;13(9):861–71.
Wang K, Long B, Liu F, Wang JX, Liu CY, Zhao B, et al. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur Heart J. 2016;37(33):2602–11.
Hansen TB, Kjems J, Damgaard CK. Circular RNA and miR-7 in cancer. Cancer Res. 2013;73(18):5609–12.
Han B, Chao J, Yao H. Circular RNA and its mechanisms in disease: From the bench to the clinic. Pharmacol Ther. 2018;187:31–44.
Bai N, Peng E, Qiu X, Lyu N, Zhang Z, Tao Y, et al. circFBLIM1 act as a ceRNA to promote hepatocellular cancer progression by sponging miR-346. J Exp Clin Cancer Res. 2018;37(1):172.
Sun L, Fang Y, Wang X, Han Y, Du F, Li C, et al. miR-302a inhibits metastasis and cetuximab resistance in colorectal cancer by targeting NFIB and CD44. Theranostics. 2019;9(26):8409–25.
This study was supported by grants from the China Postdoctoral Science Foundation [No. 2019M652393] and the Science Foundation of Qilu Hospital of Shandong University [No. 2017QLQN16]. The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing of this manuscript.
Ethics approval and consent to participate
The study was obtained the support of the ethics committee of Qilu Hospital of Shandong University and received informed consent from all patients.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Liang, Y., Shi, J., He, Q. et al. Hsa_circ_0026416 promotes proliferation and migration in colorectal cancer via miR-346/NFIB axis. Cancer Cell Int 20, 494 (2020). https://doi.org/10.1186/s12935-020-01593-1