- Open Access
Circular RNA circRHOBTB3 acts as a sponge for miR-654-3p inhibiting gastric cancer growth
Journal of Experimental & Clinical Cancer Research volume 39, Article number: 1 (2020)
Circular RNAs (circRNAs) have recently emerged as a new family of noncoding RNAs that are involved in the causation and progression of various cancers. However, the roles of circRNAs in the tumorigenesis of gastric cancer (GC) are still largely unknown.
The expression profiles of circRNAs in GC were identified in open GEO database and were evaluated at the mRNA level in clinical GC samples compared with paired non-tumorous tissues. Kaplan-Meier survival curve was used to analyze the correlation of circRNA and patients’ prognosis. Subsequently, the circular structures of candidate circRNAs were validated by Sanger sequencing, divergent primer PCR, and RNase R treatments. Gain- and loss-of-function analyses were performed to evaluate the functional significance of it in GC initiation and progression. Dual-luciferase reporter and RNA pull-down assays were used to identify the microRNA (miRNA) sponge mechanism of circRNAs.
The expression of circRHOBTB3 was lower in GC tissues and cell lines. Downregulation of circRHOBTB3 was significantly correlated with poor differentiation and unfavorable prognosis in patients with GC. Overexpression of circRHOBTB3 in GC cells led to decreased proliferation and induced G1/S arrest in vitro, accompanied with inhibited xenograft tumor growth in vivo, while the opposite effects were achieved in circRHOBTB3-silenced cells. Furthermore, we demonstrated that circRHOBTB3 acts as a sponge for miR-654-3p and verified that p21 is a novel target of miR-654-3p.
Taken together, this study revealed that circRHOBTB3 might function as competing endogenous RNA (ceRNA) for miR-654-3p, which could contribute to growth inhibition of GC through activating p21 signaling pathway. Our data suggested that circRHOBTB3 would serve as a novel promising diagnosis marker and therapeutic target for GC.
Gastric cancer (GC) is one of the most common cancers and the third leading cause of cancer-related death worldwide  . Although in recent years, development of surgical treatment and adjuvant therapies significantly improves the prognosis of GC patients, the mortality rate remains high, mainly due to the heterogeneity and complicated regulatory relation at molecular level [2,3,4,5]. Hence, better understanding of the underlying molecular mechanism in GC tumorigenesis and progression is urgently needed.
CircularRNAs (circRNAs) are a novel class of endogenous non-coding RNAs mainly formed by the mechanism of “direct back splicing” or “Exon skipping” through producing a covalently closed loop, which usually originates from exon of genes [6, 7]. Although being discovered for about 40 years, the pathological and physiological process of circRNAs remains largely unknown [8,9,10]. With the advent of high-throughput sequencing technique and novel bioinformatic analysis thousands of circRNAs have been successfully identified in multiple cell lines and across various species. These circRNAs contain many salient features including high stability, evolutionary conservation and tissue-specific or cell type-specific expression [7, 11, 12]. Moreover, several studies have found that many circRNAs were aberrant expressed and could involve in gene regulation rather than by-products of splicing or splicing errors. Also, evidences showed that circRNAs may act as sponges for microRNAs (miRNAs) or bind to proteins, and abnormal circRNAs expression could lead to alteration of gene products that may contribute to tumor biology including cell proliferation, apoptosis, angiogenesis and metastasis [13,14,15,16]. In fact, evidence from recent literatures and our previous studies has implied that circRNAs play a pivotal role in the tumorigenesis of GC [17,18,19]. These findings suggested that circRNAs might be a novel biomarker for the diagnosis and treatment of cancer. However, these studies only exhibited the preliminary results on the circRNA-miRNA regulatory network in GC. The overall pathophysiological role of circRNAs in GC needs to be further investigated.
In the present study, we identified one circRNA originating from exon 6 and exon 7 of RHOBTB3 gene and termed it circRHOBTB3. Meanwhile, the expression of circRHOBTB3 was determined in paired GC tissues and the relationship between circRHOBTB3 with patients’ clinicopathologic characteristic was also analyzed. Besides, the functions of circRHOBTB3 in the growth of GC were explored both in vitro and in vivo. Finally, we confirmed that circRHOBTB3 could sponge miR-654-3p and promoted the expression of p21, a major effector molecule of the cell cycle inhibiting protein.
Materials and methods
Patients and tissue specimens
Seventy-five pairs of GC tissue and paired non-cancerous tissue were obtained from patients with GC who had endoscopically proven primary GC and received radical surgical treatment at Nanfang Hospital (Guangzhou, China) between 2017 and 2019. All patients enrolled in this study did not receive prior surgery, radiotherapy or chemotherapy. All specimens were immediately snap-frozen in liquid nitrogen before storage at − 80 °C. Fifty nine samples with clinico-pathological parameters was used to analyze the significance and prognostic value of circRHOBTB3. Written informed consent was obtained from all patients. The procedure of human tissue samples collection was conducted in accordance with the international Ethical guidelines for biomedical research involving in human subjects. This study was approved by Southern Medical University Ethics Committee.
Human GC cell lines (AGS, HGC27, MKN45) were obtained from the Committee of Type Culture Collection of Chinese Academy of Sciences (Shanghai, China). In addition, the gastric mucosal cell GES-1 was provided by the Department of Pathology, Nanfang Hospital, Southern Medical University. The cells were cultured in RPMI Medium 1640 (Gibco, Australiaorigin) supplemented with 10% FBS (BI, China) at 37 °C in a humidified atmosphere of 5% CO2.
RNA preparation, qRT-PCR and genomic DNA purity
Total RNA was extracted from tissues and cells by using TRIzol reagent (Takara, Otsu, Japan). Total RNA was reverse transcribed to cDNA and qRT-PCR were conducted by using a SYBR Green PCR Kit (Takara, Otsu, Japan) as described [20, 21]. For miRNA detection, reverse transcription was performed and expression of miRNA was measured by All-in-One™ miRNA qRT-PCR Detection Kit (Genecopoaie, Lot#QP015) according to the use manual. Genomic DNA (gDNA) was isolated from tissues or cultured cells according to easy pure genomic DNA kit (Transgen Biotech, Lot#L61221). All primers were listed in Additional file 2: Table S1.
Nucleic acid electrophoresis
The cDNA and gDNA PCR products were detected using 4% agarose gel electrophoresis with TAE running buffer at 100 V for 30 min. DNA marker used was DL2000 (Takara, Japan). The bands were examined by UV irradiation.
Treatment with RNase R and Actinomycin D
Total RNA (2 μg) was incubated for 15 min at 37 °C with or without 3 U/mg of RNase R (Epicentre Technologies, Madison, WI, USA). Then, the resulting RNA was directly reverse transcribed using the prime script RT master mix (Takara, Japan) for real-time PCR analyses. To assess the stability of circRHOBTB3 and its linear isoform, The culture medium was added to Actinomycin D or DMSO (sigma Aldrich, St. Louis, MO, USA) as described .
Cell proliferation and cell cycle
Cell proliferation was examined using cell counting kit-8 assay (Dojindo Laboratories, Kumamoto, Japan) and the cell-light™ Edu staining kit (RiboBio, Guangzhou, China) according to manufacturer’s instructions. Colony formation assay were performed to monitor cell cloning capability as followed: 1 × 103 cells was seeded into 6-wells plates and cultured at 37 °C in atmosphere of 5% CO2 for 2 weeks. Then the colonies were washed twice with PBS, fixed with 4% paraformaldehyde for 10 min and dyed with Wright-Giemsa Stain. The number of colonies was photographed and counted using a microscope. Cell cycle was analyzed by flow cytometry. In brief, cells were collected, fixed with cold ethanol for 2 h at 37 °C, washed by PBS, stained with propidium iodide (PI) (Keygentec, Nanjing, China) containing RNase A and then the cell cycle was detected by a flow cytometer (FACS Calibur, Becton Dickinson).
Fluorescence in situ hybridization (FISH)
Cy3-labled probes including circRHOBTB3 and 18 s and FAM-labled locked nucleic acid miRNA probes were designed and synthesized by RiboBio (Guangzhou, China) and GenePharm (Suzhou, China), respectively, and the probe sequences were obtained on request. The signals of the probes were detected by Fluorescence In Situ Hybridization kit (RiboBio, Guangzhou, China) according to the manufacturer instruction. The images were acquired on Nikon AISi Laser Scanning Confocal Microscope (Nikon instruments Inc., japan).
SiRNA or miRNA mimic was designed and synthesized by Suzhou GenePharma (Suzhou, china) and RioBio (Guangzhou, china), respectively. The sequences used were shown in Additional file 2: Table S1. The cells were transfected using lipofectamine2000 (Thermo Fisher, Shanghai, China).
Over-expressing plasmids construction and stable cells transfection
To establish circRHOBTB3 over-expression plasmids, human circRHOBTB3 full-length cDNA was synthesized by vigenebio (Shandong, China) and inserted into PKO-ciR vector including a front circular frame and a back circular frame. Transfection was carried out using lipofectamine2000 (Thermo Fisher, Shanghai, China) according to manufacturer’s instructions. Stable over-expression circRHOBTB3 cells were constructed using over-expressing circRHOBTB3 lentivirus carrier (Vigenebio, Shandong, China). We designed two small interfering RNA (siRNAs) targeting the junction sites of circRHOBTB3 to silence circRHOBTB3 expression in AGS and HGC27 cell lines. These siRNA could significantly knock down circRHOBTB3 levels without affecting its linear isoform, and we selected si-circRHOBTB3–1 to insert into lentivirus carrier to establish stable silencing circRHOBTB3 cell lines due to its higher inhibitory efficacy of circRHOBTB3.
Western blot analysis was performed as illustrated [22,23,24]. In brief, cells were harvested, lysed with radio immunoprecipitation assay buffer (RIPA, Beyotime, China), and quantified by bicinchoninic acid (BCA) analysis (Beyotime, China). Then, protein extractions were separated by 10% SDS-PAGE, transferred onto polyvinylidene fluoride (PVDF) membranes (millipore Corporation, Billercica, MA, USA), washed by TBST, and incubated with a high affinity anti-p21 antibody (1:1000), and anti-Tubulin antibody (1:1000) (Cell Signaling Technology, USA) overnight at 4 °C. Subsequently, the membranes were incubated with a secondary antibody (Cell Signaling Technology, USA) (1:5000). After washes, signals were visualized by the enhanced chemiluminescence (ECL) detection system (Pierce Biotechnology, Rockford, IL, USA) conducting in accordance with the manufacturer’s instructions.
Biotinylated RNA pull-down assay
Pull down assay was carried out as described . In brief, for circRHOBTB3 pull down miRNA, about 1 × 107 cells were harvested, lysed and sonicated. The circRHOBTB3 probe was incubated with streptavidin magnetic beads (Beaver, Suzhou, china) for 2 h to generate probe-coated beads, then incubated with cell lysates, followed by eluted with Trizol (Takara, japan), and qRT-PCR. for miR-654-3p pulled down, about 1 × 107 cells were collected, lysed, sonicated, and incubated with streptavidin magnetic beads (beaver, Suzhou, china) after transfected with biotinylated miR-654-3p mimics or mutant using lipofectamine2000 (Thermo Fisher, Shanghai, China), followed by washed, eluted and qRT-PCR.
Luciferase reporter assays
For circRHOBTB3 and miRNA luciferase assay, the circRHOBTB3 sequences containing wild-type or mutated miRNA binding sites were respectively synthesized and inserted into pEZX-MT06 luciferase vector (Genecopoaie, Guangzhou, China), and then co-transfected with miRNA mimics into GC cell lines using lipofectamine2000 (Thermo Fisher, Shanghai, China). After 48 h transfection, cells were harvested, lysed, and subjected to luciferase activity detection by the Luc-pair™ Duo-Luciferase HS assay kit (Genecopoaie, Guangzhou, China). Relative luciferase activity was normalized to the Renilla luciferase internal control.
All animal experiments were approved by the Animal Care Committee and Use Committee of Southern Medical University. The experiments were in accordance with the guidelines for the ethical treatment of animals. To investigate the circRHOBTB3 growth effect on GC cells in vivo, 4-week-old male BALB/C nude mice were randomly divided into two groups (n = 5 for each group). The stable over-expressed circRHOBTB3 MKN45 cells or control group were subcutaneously injected into flank region of legs (5 × 106 cells per mouse), respectively. Tumor size were measured from perpendicular axes and calculated as followed formula: volume = (length×width2)/2. One month later, mice were sacrificed by cervical vertebra dislocation. The primary tumors were removed and harvested, and weighed. Then, the tumors were fixed, paraffin-embedded and sectioned. Subsequently, the sections were visualized under a microscope following hematoxylin and eosin (H&E) staining.
Immunohistochemical staining was performed as described . Slides from nude mice were incubated overnight with primary antibodies against ki67 (#ab15580, Abcam), p21 (#ab15580, Abcam) at 4 °C. The complex was observed by DAB complex, and the nuclei were counterstained with haematoxylin. The immunoreactivity in each section was assessed by at least two experienced pathologists and scored by semi-quantitative H-score approach .
Statistical analyses were performed using SPSS 20.0 software (IBM, Armonk, NY, USA) and Graphpad prism version 5.0 (Graphpad software, La Jolla, CA, USA). Data was expressed as mean ± SD. The statistical differences between groups were tested by two-tailed student’s t-test. The correlation between the circRHOBTB3 expression and clinicopathologic data in GC tissues and paired normal tissues were examined by χ2 test. The relationship between the circRHOBTB3 and linear RHOBTB3 or miR-654-3p was examined by Pearson’s correlation coefficients. Survival prognosis was analyzed by Kaplan-Meier curves and log-rank test for significance. p < 0.05 was considered statistically significant.
CircRHOBTB3 was significantly downregulated in GC tissues, cell lines and correlated with prognosis of patients with GC
To identify a candidate circRNA that may implicate in the tumorigenesis of GC, we firstly analyzed the high-throughput sequencing results of GC (accession code: GSE77661) from the open GEO database . Then we selected five candidate circRNAs from the RNA-seq data by the criteria: fold change > 2 (Circbank numbers: circN4BPL2L, circRELL1, circFAM120A, circFNDC3B, circRHOBTB3). Moreover, the relative expression of these five circRNAs were detected by qRT-PCR in GC and paired non-tumorous tissues (cohort 1, n = 30). Results showed that the relative expressions of circN4BPL2L displayed no significant changes, but the other four circRNAs were low-expressed. Interestingly, we found that circRHOBTB3 was the most significantly downregulated circRNA among them (P < 0.0001 and fold-change (GC/ANT) < 0.274, Fig. 1a), which suggested that circRHOBTB3 might be a circular molecular correlated with progression of GC.
To further confirm the relationship between the expression level of circRHOBTB3 and the clinical progress of GC patients. qRT-PCR was used to detected the expressions of circRHOBTB3 in newly internalized GC tissues and adjacent tissues (cohort 2, n = 75). Results displayed that circRHOBTB3 was frequently repressed in GC tissues compared with non-tumorous tissues (Fig. 1b). When primary cancers were paralleled to corresponding normal tissues, the downregulation of circRHOBTB3 (less than a 0.5-fold change) was observed in 65.4% (49/75) of patients with GC (Fig. 1c, d). Consistently, we also found that the expression of circRHOBTB3 in GC cell lines was obviously downregulated compared with GES-1 cells (the normal gastric mucosal cells) (Fig. 1e). Subsequently, these patients were stratified into high and low groups based on the median value of circRHOBTB3 expression. The relationship between circRHOBTB3 expression in GC tissues and clinicopathological features. As shown in Table 1, low expression of circRHOBTB3 was dramatically correlated with tumor stage, whereas no significant correlation with other clinicopathologic parameters. In addition, survival analyses of these patients revealed that patients with low expression of circRHOBTB3 had poor prognosis in GC than patients in the high circRHOBTB3 expression group (Fig. 1f). These collective data suggested that circRHOBTB3 might play pivot role in the tumorigenesis of GC.
It has been reported that some circRNAs may modulate the corresponding linear RNA transcripts expression and then execute function [28, 29]. Therefore, the regulatory relationship between circRHOBTB3 and its linear RNA transcript (RHOBTB3) was explored. Firstly, the expression level of RHOBTB3 was examined in the 19 paired GC and adjacent non-tumorous tissues (Fig. 1g). However, no significant changes of RHOBTB3 mRNA was observed. Pearson’s correlation analysis revealed a significant positive correlation between circRHOBTB3 and its linear RHOBTB3 in GC tissues (r = 0.54, P = 0.018, Fig. 1h). Nevertheless, RHOBTB3 did not change the mRNA expression levels when the expression of circRHOBTB3 was artificially changed in GC cells (Additional file 1: Figure S2A-C). These results indicated that RHOBTB3 is not the target gene of circ RHOBTB3.
Characteristics of circRHOBTB3
CircRHOBTB3 was generated from exon 6 and exon7 of RHOBTB3 gene (CircBase ID: hsa_circ_00074444, splicing length: 479 nucleic acid base). To further confirm circular characteristics of circRHOBTB3, the transcripts of both circRHOBTB3 and RHOBTB3 mRNA was tested by qRT-PCR in three tumor tissues, AGS and HGC27 cell lines after treatment with or without RNase R. Results showed that the fragment of linear form of RHOBTB3 gene was digested by RNase R while cirRHOBTB3 was retained after RNase R treatment (Fig. 2a, b), which verified that circRHOBTB3 was resistant to RNase R due to its loop structure. Secondly, to rule out the possibility of head-to-tail sequencing produced by trans-splicing or genomic rearrangement, Divergent primers and convergent primers were designed to amplify circRHOBTB3 and RHOBTB3 mRNA, respectively. cDNA and gDNA (genomic DNA) from three GC tissues and AGS, HGC27 cell lines were used as templates. We found that circRHOBTB3 was only amplified by divergent primers in cDNA, but no amplification product was visualized in gDNA. Meanwhile, the head-to-tail junction sequences were validated by Sanger sequencing (Fig. 2c, d). Then, inhibiting transcription experiment was utilized to reveal the stability of circRHOBTB3, and illustrated that it was more stable than its linear mRNA (Fig. 2e). Additionally, the subcellular localization of circRHOBTB3 was determined in nucleoplasmic separation and FISH experiments. Results indicated that circRHOBTB3 was preferentially localized in cytoplasm (Fig. 2f, g and Additional file 1: Figure S1). Taken together, the above results indicated that circRHOBTB3 is an abundant, circular and stable transcript that mainly localized in cytoplasm of GC cells.
CircRHOBTB3 inhibited GC cell growth and cell cycle progression in vitro
To better understand the role of circRHOBTB3 in GC cells. We selected si-circRHOBTB3–1 to insert into lentivirus carrier to establish stable silencing circRHOBTB3 (SH-circRHOBTB3) in AGS and HGC27 cell lines due to its higher inhibitory efficacy of circRHOBTB3. Data demonstrated that stable SH-circRHOBTB3 AGS and HGC27 cell lines were established successfully (Additional file 1: Figure S2A, B). Moreover, circRHOBTB3 were stably over-expressed by circRHOBTB3-overexpressed lentivirus vector in MKN45, AGS and HGC27 cells lines (Additional file 1: Figure S2C). Subsequently, functional assays were executed to reveal the effects of circRHOBTB3 on GC cell proliferation. The cck8 assay showed that knocked-down of circRHOBTB3 enhanced the growth ability of AGS and HGC27 cell lines significantly (Fig. 3a and Additional file 1: Figure S2D). Whereas, over-expression of circRHOBTB3 suppressed the proliferation of MKN45 cells (Fig. 3b). In addition, the colony formation assay indicated that SH-circRHOBTB3 AGS and HGC27 cell lines produced more clonogenicities compared with control group but elevation of circRHOBTB3 displayed opposite effect in MKN45 cells (Fig. 3c and Additional file 1: Figure S2E). Consistently, by using Edu assay, compared with control group, more proliferative cells were observed in AGS and HGC27 cell lines with down-regulation of circRHOBTB3 while less in MKN45 cells with up-regulation of circRHOBTB3 (Fig. 3d and Additional file 1: Figure S2F). Moreover, Flow cytometry analysis was performed to determine whether circRHOBTB3 affected the cell cycle profile. As presented in Fig. 3e and Additional file 1: Figure S2G, less cells stopped in G1 phase after silencing of circRHOBTB3 in AGS and HGC27 cell lines, which suggested that silencing of circRHOBTB3 induced G1/S cell cycle progression. However, over-expression of circRHOBTB3 arrested cell cycle in G1 phase in MKN45 cells (Fig. 3f). These data collectively demonstrated that circRHOBTB3 inhibited the growth phenotype of GC cells and cell cycle progression.
CircRHOBTB3 acted as a molecule sponge for miR-654-3p
Previous studies suggested that circRNAs mainly function as miRNA sponge to sequester miRNA and then regulate subsequent gene expression [30, 31]. To investigate whether circRHOBTB3 served as “miRNA sponge” in GC cells, we selected the eleven top potential binding miRNA (miR-244, miR-494, miR-548p, miR-570, miR-586, miR-600-3p, miR-619-3p, miR-654-3p, miR-1200-5p, miR-1265-5p, miR-1305) with a context score percentile≥90 predicted by CircInteractome database. The overexpression efficacy of circRHOBTB3 was quantified by qRT-PCR in AGS and HGC27 cell lines transfected with elevated circRHOBTB3 plasmid or empty (Additional file 1: Figure S3A). Subsequently, a 3′ terminal-biotinylated-circRHOBTB3 probe was designed to determine which miRNAs interacted with circRHOBTB3. As illustrated in Fig. 4a, the probe was confirmed to pull down circRHOBTB3 in AGS and HGC27 cell lines and overexpression of circRHOBTB3 increased pull-down efficacy. Then, qRT-PCR analyses revealed that several miRNAs including miR-548p, miR-570, miR-586, miR-654-3p and miR-1200-5p could all be pulled down in AGS and HGC27 cell lines. Among them, miR-654-3p was most abundantly pulled down by circRHOBTB3 in HGC27 and AGS cells (Fig. 4b, c). Meanwhile, miR-654-3p was reported and predicted by Target Scan website to target p21 protein, a well-known proliferative suppressor in various tumors [32,33,34]. Therefore, we focused on miR-654-3p for further study. To adequately consolidate the binding between circRHOBTB3 and miR-654-3p, biotin-coupled miR-654-3p and its mutant were utilized to pull down circRHOBTB3 in MKN45 and HGC27 cell lines with stable circRHOBTB3 overexpression. Results showed that wild-type miR-654-3p captured more circRHOBTB3 compared with its mutant. GAPDH served as negative control and showed no significant changes between biotin-labeled miR-654-3p and mutant group (Fig. 4d). Furthermore, we conducted luciferase reporter assay and verified that up-regulation of miR-654-3p relatively reduced the luciferase activity of vector containing full circRHOBTB3 sequences, but didn’t influence the luciferase activity of vector including mutant binding sites of miR-654-3p in AGS, HGC27, and MKN45 cell lines (Fig. 4e and Additional file 1: Figure S3B, C). Moreover, we conducted FISH assay to assess whether there was co-location between circRHOBTB3 and miR-654-3p. Results indicated that co-location of circRHOBTB3 and miR-654-3p was mainly visualized in cytoplasm (Fig. 4f and Additional file 1: Figure S3D). Based on the above, data demonstrated that circRHOBTB3 could serve as a sponge molecule for miR-654-3p.
CircRHOBTB3 inhibits GC growth through the circRHOBTB3/ miR-654-3p/p21 pathway
To address whether circRHOBTB3 plays tumor-inhibiting role via circRHOBTB3/miR-654-3p/p21 pathway in GC, we performed the following experiments: Firstly, the expression levels of miR-654-3p were detected by qRT-PCR in 30 paired GC tissues and adjacent non-cancerous tissues. Results indicated that miR-654-3p was highly expressed in GC tissues compared with control group (Additional file 1: Figure S4A). In addition, we found that the expression of miR-654-3p was negatively correlated with circRHOBTB3 in GC tissues analyzed by Pearson’s correlation (Additional file 1: Figure S4B), which cued that miR-654-3p might severe as proliferation-promoting role in GC. Consistently, the expression level of miR-654-3p was up-regulated in GC cell lines in comparison to GES-1 cells (Additional file 1: Figure S4C). Subsequently, cck8 assay was conducted, and confirmed that miR-654-3p enhanced the growth of AGS and HGC27 cell lines transfected with miR-654-3p compared with NC-mimic (Additional file 1: Figure S4D, E). Then, we examined the protein levels of p21 in AGS and HGC27 cell lines transfected with miR-654-3p or NC-mimic, and found that p21 expression was significantly repressed in miR-654-3p group compared with NC-mimic, as shown in Additional file 1: Figure S4F. Additionally, we examined the protein expression of p21 using western blot in stable over-expression or down-regulation circRHOBTB3 GC cell lines. As expected, results showed that protein expression of p21 was upregulated in AGS and HGC27 cells with stably elevated circRHOBTB3 and downregulated in SH-circRHOBTB3 AGS and HGC27 cell lines (Additional file 1: Figure S4G, H). However, qRT-PCR analyses for the mRNA levels of p21 indicated that there were no significant changes in AGS and HGC27 cell lines treated with miR-654-3p or NC-mimic or stable overexpressed circRHOBTB3 GC cell lines (Additional file 1: Figure S3I), suggesting that miR-654-3p or circRHOBTB3 didn’t affect p21 mRNA levels. To further explore whether circRHOBTB3 acted as tumor-inhibitor in GC cell by attenuating activity of miR-654-3p to up-regulate protein expression of p21, rescue experiments was performed. cck8 assay displayed that AGS cells with elevated circRHOBTB3 plus miR-654-3p grown slower than miR-654-3p group (Fig. 5a, b), and consistent results was reproduced in HGC27 cells. Likewise, the colony assay showed that overexpressed circRHOBTB3 GC cells together with miR-654-3p exerted less cloned cells than GC cells transfected with miR-654-3p (Fig. 5c). Moreover, cell cycle distribution was analyzed by flow cytometry in GC cells and results suggested that miR-654-3p combined with circRHOBTB3 reversed miR-654-3p-induced cell cycle progression, resulting in more cells ceased in G1 phase (Fig. 5d, e). Then, to see whether circRHOBTB3 could restore the expression levels of p21 inhibited by miR-654-3p, p21 was examined using western blot, and the brightness of the protein band was observed to be stronger in circRHOBTB3 + miR-654-3p group than in miR-654-3p group (Fig. 5f). These results indicated that circRHOBTB3 reversed miR-654-3p-induced enhancement of GC cell growth, and could restore the expression of miR-654-3p target p21 at protein levels, forming the circRHOBTB3/miR-654-3p/p21 regulating axis (Fig. 5h).
Stable overexpression of circRHOBTB3 inhibits GC growth in vivo by targeting p21
To assess whether overexpression of circRHOBTB3 affects tumor growth in vivo, MKN45 transfected with elevated circRHOBTB3 lentivirus or control vector was injected subcutaneously into nude mice. Haematoxylin and eosin (H&E) staining was applied to confirm the histopathological features of tumors harvested. Tumor volume shrink and reduced average tumor weight were observed in overexpressed circRHOBTB3 group compared with control group (Fig. 6a-d). Immunohistochemistry (IHC) demonstrated that the expression of ki67 proliferation antigen was significantly weakened in overexpressed circRHOBTB3 xenografts. Moreover, the staining of p21 was observed stronger than control group (Fig. 6e). Therefore, circRHOBTB3 could present tumor-suppressing property to inhibit cell growth and produce cell cycle arrest.
CircRNAs has been long considered to be molecular flukes or by-products of splicing ever since they were observed in eukaryotic cells by electron microscopy almost 40 years ago [10, 35]. However, the understanding of circRNAs has gradually changed with the application and development of the next-generation sequencing. A large number of circRNAs were detected in various animal samples and multiple cell lines, and many of them were found to be highly stable and abundantly expressed, based on which, subsequent studies demonstrated that circRNAs were dysregulated in diverse cancers including colorectal cancer [13, 30], hepatocellular carcinoma , basal cell carcinoma , laryngeal cancer , bladder cancer . It has been reported that circRNAs could serve as the prognostic biomarkers for many cancers. For example, elevation of circHIPK3 was correlated with poor prognosis in bladder cancer and colorectal cancer [20, 30]. Patients with lower circ-ITCH expression had a poorer prognosis in bladder cancer , and higher levels of circPRCKI displayed as a worse progression marker in lung adenocarcinoma . Whereas, their roles in GC remained largely unknown. In this study, we found that circRHOBTB3 was frequently low-expressed in GC tissues compared with adjacent normal tissues. The low expression of circRHOBTB3 was correlated with tumor stage and provided the poor prognosis of GC patients. Astoundingly, clinical correlation analyses showed that tumor size was not associated with circRHOBTB3 expression, which was not consistent with its functional effects in GC. One possible reason might be that the number of GC specimens was too few. Thus, it’s necessary to analyze the relationship between circRHOBTB3 and clinicopathological parameter using more GC tissue samples in the future study.
CircRNAs could act as oncogenic or tumor suppressive factors in GC. CircFAT1 displayed inhibiting effects in GC , repressing proliferation and invasion of GC cells, but circPVT1/circDLST enhanced the malignancy of GC [19, 39]. Currently, we identified a functional role of circRHOBTB3 in GC cells and validated that re-expression of circRHOBTB3 suppressed the growth of GC in vitro and vivo and arrested cell cycle in G1 phase, while knockdown of circRHOBTB3 produced contrary effects. These results indicated that circRHOBTB3 might serve as a potential tumor suppressive marker in GC.
Recently, increasing evidences have shown that circRNAs functions as sponge for miRNAs to affect tumorous biological process . For instance, circSMARCA5 could act as sponge for miR-17-3p and miR-181b-5p to inhibit the growth and metastasis in hepatocellular carcinoma . CircPVT1 interacted with miR-125b to exert proliferative effect in GC . Herein, we found that circRHOBTB3 contained the binding sites of miR-654-3p predicted by bioinformatic website circinteractome and verified there was direct binding relationship between circRHOBTB3 and miR-654-3p using biotinylated-nucleic probe pull down assay and luciferase reporter experiments, which suggested that circRHOBTB3 might serve as a sponge for miR-654-3p to inhibit the proliferation of GC.
MiRNAs have been shown to be a large family of gene regulator that negatively regulate their target mRNAs in a sequence-specific manner and its role in cancers have been extensively studied [31, 42, 43]. Numerous evidences have confirmed that miRNAs play essential roles in multiple biological process related to cancer, including cell differentiation, proliferation, tumorigenesis, angiogenesis, invasion, and metastasis [44, 45]. As described, miR-654-3p was confirmed to target p21 mRNA and induced its protein downregulation in HEK293 cells . Previous studies have demonstrated that p21 could suppress multiple tumorous proliferation by targeting cyclin-dependent kinase (CDK) complexes [46, 47]. However, whether miR-654-3p could regulate the expression of p21 and biological function in GC cells remained unknown. In the present study, we proved that miR-654-3p was up-regulated in GC tissues and cells, and promoted the proliferation of GC cells. Additionally, we also confirmed that miR-654-3p could inhibit the protein expression of p21 in GC cells using western blot, whereas, circRHOBTB3 enhanced the protein expression of p21. Further reversed experiments were conducted and we found that circRHOBTB3 inhibited the activity of miR-654-3p, and thus, reversed the miR-654-3p induced proliferation of GC cells and up-regulated the protein expression of p21 indirectly. Hence, we validated that circRHOBTB3 could inhibit the proliferation of GC by sponging miR-654-3p to up-regulate the protein expression of p21. One thing we have to pay attention is that not all circRNAs can function as “miRNA sponges” . Intronic circRNAs and exon-intronic circRNAs, which predominantly localize in nucleus with lacking enrichment for miRNA target sites, have been reported to modulate parental genes expression via RNA-RNA interaction [29, 49]. Meanwhile, some circRNAs, such as circMbl, circFmn, circDMD, can strongly bind to cognate linear transcripts to prevent mRNA from translation and finally lead to the reduction in protein levels [48, 50]. This process is also termed as “mRNA trap”. Thus, diverse functions of circRNAs in GC need to be further explored.
Our study confirmed that circRHOBTB3 was frequently down-regulated in GC and low expression of circRHOBTB3 was correlated with tumor stage and patients’ unfavorable prognosis. Mechanistically, circRHOBTB3 could serve as a sponge molecule for miR-654-3p, and then rescued the expression of miR-654-3p-inbiting protein p21, finally suppressed GC cell growth. Our findings might provide new insight into GC development and provided a novel potential strategy for GC treatment.
Availability of data and materials
The datasets used and analysed during the current study are available from the corresponding author on reasonable request.
Competing endogenous RNA
Circular RNA RHOBTB3
Fluorescence in situ hybridization
Hematoxylin and eosin
Quantitative real-time polymerase chain reaction
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.
Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136(5):E359–86.
Tan P, Yeoh KG. Genetics and molecular pathogenesis of gastric adenocarcinoma. Gastroenterology. 2015;149(5):1153–62.
Wang K, Kan J, Yuen ST, Shi ST, Chu KM, Law S, et al. Exome sequencing identifies frequent mutation of ARID1A in molecular subtypes of gastric cancer. Nat Genet. 2011;43(12):1219–23.
Nakamura J, Tanaka T, Kitajima Y, Noshiro H, Miyazaki K. Methylation-mediated gene silencing as biomarkers of gastric cancer: a review. World J Gastroenterol. 2014;20(34):11991–2006.
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.
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.
Capel B, Swain A, Nicolis S, Hacker A, Walter M, Koopman P, et al. Circular transcripts of the testis-determining gene Sry in adult mouse testis. Cell. 1993;73(5):1019–30.
Nigro JM, Cho KR, Fearon ER, Kern SE, Ruppert JM, Oliner JD, et al. Scrambled exons. Cell. 1991;64(3):607–13.
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.
Barrett SP, Salzman J. Circular RNAs: analysis, expression and potential functions. Development. 2016;143(11):1838–47.
Rybak-Wolf A, Stottmeister C, Glazar P, Jens M, Pino N, Giusti S, et al. Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed. Mol Cell. 2015;58(5):870–85.
Bachmayr-Heyda A, Reiner AT, Auer K, Sukhbaatar N, Aust S, Bachleitner-Hofmann T, et al. Correlation of circular RNA abundance with proliferation--exemplified with colorectal and ovarian cancer, idiopathic lung fibrosis, and normal human tissues. Sci Rep. 2015;5:8057.
Salzman J, Chen RE, Olsen MN, Wang PL, Brown PO. Cell-type specific features of circular RNA expression. PLoS Genet. 2013;9(9):e1003777.
Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J, et al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA. 2013;19(2):141–57.
Salzman J, Gawad C, Wang PL, Lacayo N, Brown PO. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PLoS One. 2012;7(2):e30733.
Zhao ZJ, Shen J. Circular RNA participates in the carcinogenesis and the malignant behavior of cancer. RNA Biol. 2017;14(5):514–21.
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.
Chen J, Li Y, Zheng Q, Bao C, He J, Chen B, et al. Circular RNA profile identifies circPVT1 as a proliferative factor and prognostic marker in gastric cancer. Cancer Lett. 2017;388:208–19.
Li Y, Zheng F, Xiao X, Xie F, Tao D, Huang C, et al. CircHIPK3 sponges miR-558 to suppress heparanase expression in bladder cancer cells. EMBO Rep. 2017;18(9):1646–59.
Qiu M, Xu Y, Wang J, Zhang E, Sun M, Zheng Y, et al. A novel lncRNA, LUADT1, promotes lung adenocarcinoma proliferation via the epigenetic suppression of p27. Cell Death Dis. 2015;6:e1858.
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(1):19.
Qiu M, Xia W, Chen R, Wang S, Xu Y, Ma Z, et al. The circular RNA circPRKCI promotes tumor growth in lung adenocarcinoma. Cancer Res. 2018;78(11):2839–51.
Shi J, Li F, Yao X, Mou T, Xu Z, Han Z, et al. The HER4-YAP1 axis promotes trastuzumab resistance in HER2-positive gastric cancer by inducing epithelial and mesenchymal transition. Oncogene. 2018;37(22):3022–38.
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.
Wang H, Deng G, Ai M, Xu Z, Mou T, Yu J, et al. Hsp90ab1 stabilizes LRP5 to promote epithelial-mesenchymal transition via activating of AKT and Wnt/beta-catenin signaling pathways in gastric cancer progression. Oncogene. 2019;38(9):1489–507.
Zeng K, Wang Z, Ohshima K, Liu Y, Zhang W, Wang L, et al. BRAF V600E mutation correlates with suppressive tumor immune microenvironment and reduced disease-free survival in Langerhans cell histiocytosis. Oncoimmunology. 2016;5(7):e1185582.
Li X, Wang J, Zhang C, Lin C, Zhang J, Zhang W, et al. Circular RNA circITGA7 inhibits colorectal cancer growth and metastasis by modulating the Ras pathway and upregulating transcription of its host gene ITGA7. J Pathol. 2018;246(2):166–79.
Zhang Y, Zhang XO, Chen T, Xiang JF, Yin QF, Xing YH, et al. Circular intronic long noncoding RNAs. Mol Cell. 2013;51(6):792–806.
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.
Ye G, Huang K, Yu J, Zhao L, Zhu X, Yang Q, et al. MicroRNA-647 targets SRF-MYH9 Axis to suppress invasion and metastasis of gastric Cancer. Theranostics. 2017;7(13):3338–53.
Luo J, Liu K, Yao Y, Sun Q, Zheng X, Zhu B, et al. DMBX1 promotes tumor proliferation and regulates cell cycle progression via repressing OTX2-mediated transcription of p21 in lung adenocarcinoma cell. Cancer Lett. 2019;453:45–56.
Zhu L, Wang J, Kong W, Huang J, Dong B, Huang Y, et al. LSD1 inhibition suppresses the growth of clear cell renal cell carcinoma via upregulating P21 signaling. Acta Pharm Sin B. 2019;9(2):324–34.
Liu YW, Xia R, Lu K, Xie M, Yang F, Sun M, et al. LincRNAFEZF1-AS1 represses p21 expression to promote gastric cancer proliferation through LSD1-mediated H3K4me2 demethylation. Mol Cancer. 2017;16(1):39.
Cocquerelle C, Mascrez B, Hetuin D, Bailleul B. Mis-splicing yields circular RNA molecules. FASEB J. 1993;7(1):155–60.
Sand M, Bechara FG, Sand D, Gambichler T, Hahn SA, Bromba M, et al. Circular RNA expression in basal cell carcinoma. Epigenomics. 2016;8(5):619–32.
Xuan L, Qu L, Zhou H, Wang P, Yu H, Wu T, et al. Circular RNA: a novel biomarker for progressive laryngeal cancer. Am J Transl Res. 2016;8(2):932–9.
Fang J, Hong H, Xue X, Zhu X, Jiang L, Qin M, et al. A novel circular RNA, circFAT1(e2), inhibits gastric cancer progression by targeting miR-548g in the cytoplasm and interacting with YBX1 in the nucleus. Cancer Lett. 2019;442:222–32.
Zhang J, Hou L, Liang R, Chen X, Zhang R, Chen W, et al. CircDLST promotes the tumorigenesis and metastasis of gastric cancer by sponging miR-502-5p and activating the NRAS/MEK1/ERK1/2 signaling. Mol Cancer. 2019;18(1):80.
Han D, Li J, Wang H, Su X, Hou J, Gu Y, et al. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology. 2017;66(4):1151–64.
Yu J, Xu QG, Wang ZG, Yang Y, Zhang L, Ma JZ, et al. Circular RNA cSMARCA5 inhibits growth and metastasis in hepatocellular carcinoma. J Hepatol. 2018;68(6):1214–27.
Calin GA, Croce CM. MicroRNA-cancer connection: the beginning of a new tale. Cancer Res. 2006;66(15):7390–4.
Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, et al. A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci U S A. 2006;103(7):2257–61.
Calin GA, Croce CM. MicroRNA signatures in human cancers. Nat Rev Cancer. 2006;6(11):857–66.
Gregory RI, Shiekhattar R. MicroRNA biogenesis and cancer. Cancer Res. 2005;65(9):3509–12.
Wu S, Huang S, Ding J, Zhao Y, Liang L, Liu T, et al. Multiple microRNAs modulate p21Cip1/Waf1 expression by directly targeting its 3′ untranslated region. Oncogene. 2010;29(15):2302–8.
Harper JW, Adami GR, Wei N, Keyomarsi K, Elledge SJ. The p21 Cdk-interacting protein Cip1 is a potent inhibitor of G1 cyclin-dependent kinases. Cell. 1993;75(4):805–16.
Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol. 2014;32(5):453–61.
Li Z, Huang C, Bao C, Chen L, Lin M, Wang X, et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 2015;22(3):256–64.
Ashwal-Fluss R, Meyer M, Pamudurti NR, Ivanov A, Bartok O, Hanan M, et al. circRNA biogenesis competes with pre-mRNA splicing. Mol Cell. 2014;56(1):55–66.
This work was supported by the grants from the State’s Key Project of Research and Development Plan (2017YFC0108300), National Natural Science Foundation of China (81672446), Guangdong Provincial Science and Technology Key Project (2014A020215014), Research Fund of Public Welfare in the Health Industry, the National Health and Family Planning Commission of China (201402015), the Southern Medical University Clinical Research Start-Up Project (LC2016ZD003), the Key Clinical Specialty Discipline Construction Program (121), and Presidential Foundation of Nanfang Hospital, Southern Medical University (2016C034) (2018B002).
Ethics approval and consent to participate
The collection of human specimens was approved by the Ethics Committee of Nanfang Hospital and the written informed consent were obtained from each patient.
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.
Additional file 1: Figure S1.
RNA FISH was carried out to detect circRHOBTB3’s subcellular localization in AGS cells. Figure S2. Silencing of circRHOBTB3 promoted proliferation and progression of cell cycle in HGC27 cells. Figure S3. CircRHOBTB3 served as sponge of miR-654-3p. Figure S4. CircRHOBTB3 modulated the expression of endogenous miR-654-3p target p21.
Additional file 2: Table S1.
Primers and RNA sequences used in this study.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
About this article
Cite this article
Deng, G., Mou, T., He, J. et al. Circular RNA circRHOBTB3 acts as a sponge for miR-654-3p inhibiting gastric cancer growth. J Exp Clin Cancer Res 39, 1 (2020). https://doi.org/10.1186/s13046-019-1487-2
- Gastric cancer