- Primary research
- Open Access
Identification of differentially expressed circRNAs and a novel hsa_circ_0000144 that promote tumor growth in gastric cancer
Cancer Cell International volume 19, Article number: 268 (2019)
Circular RNAs (circRNAs) are involved in regulating tumor pathogenesis. The mechanism of circRNAs in gastric cancer (GC) is still unknown. Our study aimed to identify differentially expressed circRNAs and assess a novel circRNA (hsa_circ_0000144) in the proliferation, migration, and invasion in GC.
Gene ontology (GO) enrichment and analyses of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, pathway network, and the ceRNA regulatory network of hsa_circ_0000144 targeting miRNAs and mRNAs were performed with the help of bioinformatics using R language and Perl software. hsa_circ_0000144 expression and circRNA knockdown in GC cell lines were detected using quantitative PCR (qPCR) in vitro. Cell proliferation, migration, and invasion after circRNA knockdown were measured using the cell counting kit-8 assay and Transwell assay.
The circRNA expression profile GSE78092 downloaded from the Gene Expression Omnibus database included three GC patients and three normal tissues. Thirty-two differentially expressed circRNAs comprised six upregulated circRNAs and 26 downregulated circRNAs. In particular, the ErbB signaling pathway, neurotrophin signaling pathway, cellular senescence, and pathways in bladder cancer and GC played the most important roles in the pathway network. The expression of hsa_circ_0000144 was upregulated in GC cell lines. Hsa_circ_0000144 knockdown suppressed tumor growth in vitro.
Hsa_circ_0000144 promotes GC cell proliferation, migration, and invasion, and the ceRNA regulatory network of hsa_circ_0000144 targeting miRNAs and mRNAs might be biomarkers for GC diagnosis and targeted therapy.
Gastric cancer (GC) is the fifth most common cancer and the third leading cause of cancer-related deaths . In 2018, there were over 1,000,000 new cases and 783,000 deaths. Although many novel therapeutic pathways and diagnosis technologies have been used in GC, the clinical outcome of patients with GC remains very poor, with a 5-year overall survival of GC < 30% . The mechanisms underlying GC pathogenesis need to be clarified to improve the clinical outcome.
Non-coding RNAs (ncRNAs) that regulate proliferation and invasion of GC cells have been evaluated in previous studies. MicroRNAs (miRNAs) and long ncRNAs (lncRNAs) play important roles in tumor biology [3,4,5]. Moreover, circular RNAs (circRNAs) are involved in various human diseases, particularly cancers [6,7,8,9,10].
Hsa_circ_0000144, which was generated from the back splicing of SLAMF6 first intron, is also known as circSLAMF6. It has been previously detected in bladder cancer. However, the underlying mechanism of hsa_circ_0000144 in GC is still unknown. In this study, we demonstrate that circSLAMF6 promotes tumor proliferation and invasion in GC and explore the regulatory network of circRNA targeting miRNA–mRNA. These data provide evidence for targeting circSLAMF6 to explore therapeutic measures and mechanisms underlying the pathogenesis of GC.
Materials and methods
ArrayStar microarray data
We used the keywords “gastric cancer circRNA” to search the National Centre of Biotechnology Information (NCBI) Gene Expression Omnibus database (GEO, https://www.ncbi.nlm.nih.gov/geo/), and the gene chip data sets GSE78092 (https://www.ncbi.nlm.nih.gov/gds/?term=GSE78092), submitted by Huang et al., were downloaded. GSE78092 contained six samples including three GC tissues and three normal tissues based on the GPL21485 platform of the ArrayStar human circular RNA microarray V2.0. The miRNA expression, clinical, and meta- and manifest data on GC were downloaded from The Cancer Genome Atlas (TCGA).
Identification of differential circRNA expression
The circRNA expression data were converted using R language and Perl software. The circRNAs in these platforms were named according to international standard names. Differentially expressed circRNAs were identified by the limma package in the Bioconductor package (source (https://bioconductor.org/biocLite.R, biocLite ("limma")). P-values of 0.05 and log-fold change of 2 were considered evidence of significant difference. Volcano plots and heat maps were constructed using limma and pheatmap packages, respectively.
miRNA target gene prediction and prognosis of miRNAs
We obtained the miRNAs targeting hsa_circ_0000144 in the circinteractome database (https://circinteractome.nia.nih.gov/bin/mirnasearch) and built an miRNA txt file. We then, downloaded four files, including miRDB.tsv, miRTarBase.tsv, and TargetScan.tsv, and identified miRNA targeted genes. Next, we extracted the expression and survival time of these miRNAs. The overall survival curves were analyzed using the R language survival package.
Interaction competitive endogenous RNA (ceRNA) network of circRNA–miRNA–mRNA
Type and network txt files were designed. The interaction ceRNA network of circRNA–miRNA–mRNA was constructed using Cytoscape software 3.6.0 .
GO enrichment and KEGG pathway analyses
GO enrichment and KEGG pathway analyses were performed using R language and Perl software. We installed the packages ("colorspace") ("dose"), biocLite ("DOSE") ("clusterProfiler"), and ("pathview").
Cell culture and transfection
The GC cell line MGC-803 and normal cell line GES-1 were purchased from the Chinese Academy of Sciences (Shanghai, China). MGC-803 cells were cultured in F-12K and DMEM-H medium (Gibco, USA), respectively. All cells were cultured at 37 °C for 18 to 24 h in a humidified incubator containing 5% CO2. MGC-803 cells were transfected with plasmids using Lipofectamine 2000 reagent (Invitrogen, USA) based on the manufacturer's instructions. The expression of hsa_circ_0000144 in the transfected cells was detected by quantitative PCR (qPCR).
CCK-8 cell proliferation assay
After 24 h of incubation, transfected cells were plated onto 96-well plates and cultured for 48 h. Every well contained 3000 cells. 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) solution was added to each well, and cell viability was assessed by measuring the absorbance at 450 nm.
Transwell migration and invasion assay
To examine cell migration and invasion ability, we conducted the Transwell assay according to the manufacturer’s instructions. Cells were incubated for 24 h. Three microscopy fields were randomly selected to acquire images.
Real-time qRT-PCR analysis
Total RNA was extracted using TRIzol reagent (Invitrogen, USA) and synthesized into cDNA using M-MLV reverse transcriptase (TaKaRa Bio, Japan) following the manufacturer’s instructions. qRT-PCR was performed using SYBR Green assay (Roche, Switzerland). Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and U6 were used as controls. The primer sequences were:
circRNA expression data analysis was performed using GraphPad Prism 6.0 software, and survival analysis was performed using the R language survival package. P < 0.05 was considered significant.
Identification of differential expression of circRNAs
We analyzed the data GSE78092 from the GEO database containing three GC tissues and three adjacent normal tissues. A total of 32 circRNAs were differentially expressed in the GC samples compared to the normal samples. Among the differentially expressed circRNAs, 26 were downregulated and six were upregulated. These results are presented in the volcano plot and heat map (Fig. 1). The top 12 circRNAs with the most significant expression were hsa_circ_0009172, hsa_circ_0089153, hsa_circ_0005927, hsa_circ_0040039, hsa_circ_0000993, hsa_circ_0092341, hsa_circ_0005075, hsa_circ_0077527, hsa_circ_0007158, hsa_circ_0042986, hsa_circ_0092342, and hsa_circ_0000144.
miRNAs and target gene prediction
We used the keyword “hsa_circ_0000144” in the circinteractome database (https://circinteractome.nia.nih.gov/bin/mirnasearch) and obtained the miRNAs that could bind to the circRNA hsa_circ_0000144. These included hsa-mir-1178, hsa-mir-1276, hsa-mir-197, hsa-mir-217, hsa-mir-485-3p, hsa-mir-502-5p, hsa-mir-526b, hsa-mir-532-3p, hsa-mir-554, hsa-mir-580, hsa-mir-583, hsa-mir-610, hsa-mir-623, and hsa-mir-942. The relationship of miRNAs targeting genes is listed in Table 1. The relationship of miRNA expression with survival is depicted in Fig. 2; hsa-mir-217 and hsa-mir-942 were significantly associated with survival (P < 0.05).
Constructed interaction of the circRNA–miRNA–mRNA ceRNA network
We designed the type and network txt files. The circRNA–miRNA–mRNA regulatory network was constructed using Cytoscape software 3.6.0. Hsa_circ_0000144 interacted with miRNAs and genes (Fig. 3).
GO enrichment and KEGG pathway analyses
GO enrichment and KEGG pathway analyses were performed using R language and Perl software. GO functional enrichment analysis of genes with P < 0.05 was obtained. The KEGG result was shown in Additional file 1. Nuclear hormone receptor binding was the most significant enrichment (Fig. 4a). The signaling pathways of differentially expressed circRNAs were mainly enriched in the ErbB signaling pathway, neurotrophin signaling pathway, cellular senescence, and pathways involved in bladder cancer and GC (Fig. 4b).
Expression of hsa_circ_0000144 in GC tumor cell lines and normal cell lines
Hsa_circ_0000144 was highly expressed in GSE78902 samples from the GEO database (Fig. 5a). To further verify the results, we chose two cell lines and found that the expression of hsa_circ_0000144 was upregulated in MGC-803 GC cancer cells lines compared to GES-1 normal cells (Fig. 5b).
Expression of hsa_circ_0000144 in different groups
To assess the interference efficiency of different siRNA sequences on hsa_circ_0000144 expression, the expression of hsa_circ_0000144 mRNA was detected by qPCR after transfection of MGC-803 cells with different siRNA sequences. Hsa_circ_0000144 expression with small interfering (si)RNA was lower than with siNC (Fig. 6a).
SiRNA-hsa_circ_0000144 inhibits proliferation in AGS cell lines
The CCK8 assay was used to evaluate the effect of hsa_circ_0000144 on proliferation. SiRNA-hsa_circ_0000144 markedly inhibited cell proliferation (Fig. 6b).
Hsa_circ_0000144 knockdown inhibits migration and invasion of MGC-803 cells
Transwell migration and invasion assays were performed to further examine the effect of hsa_circ_0000144 on GC metastasis. Hsa_circ_0000144 significantly affected cell invasion and migration (Fig. 6c, d).
Globally, GC is a common disease and is the second leading cause of cancer-related deaths [12, 13]. The mechanism underlying the pathogenesis of GC is currently unknown. The regulatory mechanism of GC should be clarified. CircRNAs play different roles in different tissues and cell types . The emerging function of circRNAs in tumorigenesis has garnered recent interest. Herein, we explored the role of hsa_circ_0000144, derived from the SLAF6 gene, in the pathogenesis of GC.
We analyzed data on GC circRNAs from the GEO database using R software and bioinformatics. Thirty-two differentially expressed circRNAs were identified, including six upregulated circRNAs and 26 downregulated circRNAs. Hsa_circ_0000144, hsa_circ_0005529, hsa_circ_0023642, hsa_circ_0061274, hsa_circ_0008035, and hsa_circ_0032821 were upregulated. Hsa_circ_0000144 was the most significantly different circRNA expressed.
From the circinteractome database, we obtained 14 miRNAs targeting hsa_circ_0000144 and tested them using the miRNA expression data of the TCGA database. The expression of hsa-mir-217 and hsa-mir-942 was associated with increased survival in GC (Fig. 2b, j). Previous studies reported that miR-217 has different roles in tumor proliferation, migration, and invasion [15,16,17,18]. Presently, the high expression of miR-942 was associated with better prognosis in GC (Fig. 2b). Conversely, another study showed that miR-942 promoted tumor pathogenesis . This dichotomy underscores the need to further assess the mechanism underlying the role of miR-942 in GC.
The ceRNA of the lncRNA–miRNA–mRNA network affects tumor regulation. CircRNAs can function as ceRNAs in tumor biology [20, 21]. Recent studies have indicated that circRNAs may act as ceRNAs to sequester miRNAs of a particular family, thereby serving as competitive inhibitors that suppress the ability of a miRNA to bind its mRNA targets [22, 23]. A previous study demonstrated that hsa_circ_0000144 serves as an endogenous sponge for miR-217 and validated that hsa_circ_0000144 directly binds to miR-217 using a luciferase reporter assay in bladder cancer cells .
Presently, we analyzed the ceRNA network of circRNA-hsa_circ_0000144 targeting miRNA–mRNA. We found that hsa_circ_0000144 barely interacted with miRNAs and mRNAs (Fig. 3). To detect the gene signal and pathway networks, GO enrichment and KEGG pathway analyses were performed. GO enrichment analysis revealed that genes that interacted with circRNAs by miRNAs were mainly involved in pathways related to nuclear hormone receptor binding, which had the most significant enrichment (Fig. 4a). The signaling pathways of differentially expressed circRNAs were mainly enriched in the ErbB signaling, neurotrophin signaling, and cellular senescence pathways, and pathways in bladder cancer and GC (Fig. 4b). Recent observations indicate that the ErbB signaling pathway contributes to pancreatic tumorigenesis [24, 25]. ErbB signaling is a prerequisite for maintenance of the intestinal epithelium following injury and tumor formation [26, 27]. Neurotrophins are a family of closely structurally related growth factors , which include nerve growth factor, brain-derived neurotrophic factor, neurotrophin 3 (NT-3), and NT-4/5. Neurotrophins exert a range of effects on cell proliferation and migration in non-neuronal cells as well as in cancer cells . Recent findings have revealed that neurotrophin signaling pathways are also a driver of tumor neurogenesis, via the stimulation of NGF receptors on nerve endings.
Recent evidence has also demonstrated that circRNAs affect functions in cancer [30, 31]. Although the mechanism of circRNAs has been reported in many cancers, its mechanism underlying carcinogenesis and cancer progression in GC has not been fully elucidated. In the present study, we identified a novel circRNA, hsa_circ_0000144, from the GEO database, which was increased in GSE78092 (Fig. 1). We detected the expression of circRNA hsa_circ_0000144 in MGC-803 GC cells and GES-1 normal cells and found that hsa_circ_0000144 was significantly upregulated in MGC-803 cells. Further analysis of the GC phenotype in vivo showed that hsa_circ_0000144 knockdown suppressed bladder cancer cell proliferation, invasion, and migration in vitro (Fig. 6b–d).
Hsa_circ_0000144 promotes GC cell proliferation and invasion. The ceRNA regulatory network of hsa_circ_0000144 targeted miRNAs and mRNAs, indicating the potential value as biomarkers for GC diagnosis and targeted therapy.
Availability of data and materials
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;0:23–31.
Allemani C, Weir HK, Carreira H, Harewood R, Spika D, Wang XS, et al. Global surveillance of cancer survival 1995–2009: analysis of individual data for 25,676,887 patients from 279 population-based registries in 67 countries (CONCORD-2). Lancet. 2015;385:977–1010.
Wang QX, et al. Altered MiRNA expression in gastric cancer: a systematic review and meta-analysis. Cell Physiol Biochem. 2015;35:933–44.
Huang YK, Yu JC. Circulating microRNAs and long non-coding RNAs in gastric cancer diagnosis: an update and review. World J Gastroenterol. 2015;21:9863–86.
Virgilio E, et al. Long non-coding RNAs in the gastric juice of gastric cancer patients. Pathol Res Pract. 2018;214:1239–46.
Rong D, Tang W, Li Z, Zhou J, Shi J, Wang H, Cao H. Novel insights into circular RNAs in clinical application of carcinomas. Onco Targets Ther. 2017;10:2183–8.
Chen Y, et al. Circular RNA circAGO2 drives cancer progression through facilitating HuR-repressed functions of AGO2–miRNA complexes. Cell Death Differ. 2019;26(7):1346–64.
Li Y, et al. CircHIPK3 sponges miR-558 to suppress heparanase expression in bladder cancer cells. EMBO Rep. 2017;18:1646–59.
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:1151–64.
Funaki NO, Tanaka J, Kasamatsu T, Ohshio G, Hosotani R, Okino T, Imamura M. Identification of carcinoembryonic antigen mRNA in circulating peripheral blood of pancreatic carcinoma and gastric carcinoma patients. Life Sci. 1996;59:2187–99.
Franz M, Lopes CT, Huck G, Dong Y, Sumer O, Bader GD. Cytoscape.js: a graph theory library for visualisation and analysis. Bioinformatics. 2016;32:309–11.
Guo J, Miao Y, Xiao B, Huan R, Jiang Z, Meng D, Wang Y. Differential expression of microRNA species in human gastric cancer versus non-tumorous tissues. J Gastroenterol Hepatol. 2006;24:652–7.
Li T, Mo X, Fu L, Xiao B, Guo J. Molecular mechanisms of long noncoding RNAs on gastric cancer. Oncotarget. 2016;7:8601–12.
Starke S, Jost I, Rossbach O, et al. Exon circularization requires canonical splice signals. Cell Rep. 2015;10:103–11.
Lin Y, Cheng K, Wang T, Xie Q, Chen M, Chen Q, Wen Q. miR-217 inhibits proliferation, migration, and invasion via targeting AKT3 in thyroid cancer. Biomed Pharmacother. 2017;95:1718–24.
Yu B, Ye X, Du Q, Zhu B, Zhai Q, Li XX. The Long Non-coding RNA CRNDE promotes colorectal carcinoma progression by competitively binding miR-217 with TCF7L2 and enhancing the Wnt/β-catenin signaling pathway. Cell Physiol Biochem. 2017;41:2489–502.
Zhao WG, Yu SN, Lu ZH, Ma YH, Gu YM, Chen J. The miR-217 microRNA functions as a potential tumor suppressor in pancreatic ductal adenocarcinoma by targeting KRAS. Carcinogenesis. 2010;31:1726–33.
Huang W, Lu Y, Wang F, Huang X, Yu Z. Downregulation of circular RNA hsa_circ_0000144 inhibits bladder cancer progression via stimulating miR-217 and suppressing RUNX2 expression. Gene. 2018;678:337–42.
Shan Z, An N, Qin J, Yang J, Sun H, Yang W. Long non-coding RNA Linc00675 suppresses cell proliferation and metastasis in colorectal cancer via acting on miR-942 and Wnt/β-catenin signaling. Biomed Pharmacother. 2018;101:769–76.
Qi X, et al. ceRNA in cancer: possible functions and clinical implications. J Med Genet. 2015;52:710–8.
Ergun S, Oncocers OS. ceRNA-mediated cross-talk by sponging miRNAs in oncogenic pathways. Tumour Biol. 2015;36:3129–36.
Chen Y, Li C, Tan C, Liu X. Circular RNAs: a new frontier in the study of human diseases. J Med Genet. 2016;53:359–65.
Zhao ZJ, Shen J. Circular RNA participates in the carcinogenesis and the malignant behavior of cancer. RNA Biol. 2017;14:514–21.
Ardito CM, Gruner BM, Takeuchi KK, et al. EGF receptor is required for KRAS-induced pancreatic tumorigenesis. Cancer Cell. 2012;22:304–17.
Navas C, Hernandez-Porras I, Schuhmacher AJ, et al. EGF receptor signaling is essential for k-ras oncogene-driven pancreatic ductal adenocarcinoma. Cancer Cell. 2012;22:318–30.
Lee D, et al. Tumor-specific apoptosis caused by deletion of the ERBB3 pseudo-kinase in mouse intestinal epithelium. J Clin Invest. 2009;119:2702–13.
Roberts RB, et al. Importance of epidermal growth factor receptor signaling in establishment of adenomas and maintenance of carcinomas during intestinal tumorigenesis. Proc Natl Acad Sci USA. 2002;99:1521–6.
Skaper SD. The neurotrophin family of neurotrophic factors: an overview. Methods Mol Biol. 2012;846:1–12.
Bradshaw RA, Pundavela J, Biarc J, Chalkley RJ, Burlingame AL, Hondermarck H. NGF and ProNGF: regulation of neuronal and neoplastic responses through receptor signaling. Adv Biol Regul. 2015;58:16–27.
Li F, Zhang L, Li W, et al. Circular RNA ITCH has inhibitory effect on ESCC by suppressing the Wnt/β-catenin pathway. Oncotarget. 2015;6:6001–133.
Chen X, Fan S, Song E. Noncoding RNAs: new players in cancers. Adv Exp Med Biol. 2016;927:1–47.
This work was supported by The National Natural Science Fund of China (81570375), The Scientific and Technological Project of Tianjin Binhai New Area Health Committee (2017BWKY036) and The Incubation Fund Project of Tianjin Medical University General Hospital (ZYYFY2016021).
Ethics approval and consent to participate
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
Wei, J., Wang, J., Gao, X. et al. Identification of differentially expressed circRNAs and a novel hsa_circ_0000144 that promote tumor growth in gastric cancer. Cancer Cell Int 19, 268 (2019). https://doi.org/10.1186/s12935-019-0975-y
- circRNA–miRNA–mRNA regulatory network
- Gastric cancer