Identification of differentially expressed circRNAs and a novel hsa_circ_0000144 that promote tumor growth in gastric cancer

Background 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. Methods 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. Results 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. Conclusions 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.


Background
Gastric cancer (GC) is the fifth most common cancer and the third leading cause of cancer-related deaths [1]. 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% [2]. The mechanisms underlying GC pathogenesis need to be clarified to improve the clinical outcome.
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.

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=GSE78 092), 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 ://bioco nduct or.org/biocL ite.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 ://circi ntera ctome .nia. nih.gov/bin/mirna searc h) and built an miRNA txt file. We then, downloaded four files, including miRDB.tsv, miR-TarBase.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 [11].

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% CO 2 . 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).

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.

Statistical analyses
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.

Constructed interaction of the circRNA-miRNA-mRNA ceRNA network
We designed the type and network txt files. The cir-cRNA-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 Fig. 1 Differential expression of circRNAs from GEO database GSE78092 samples. a Hierarchical clustering heat map of differential expression of circRNAs. b Volcano plot. Red represents upregulated circRNAs and green represents downregulated circRNAs 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).

Discussion
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 [14]. 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 [19]. 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 ceR-NAs 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 [18].
Presently, we analyzed the ceRNA network of cir-cRNA-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 [28], 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 [29]. 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 circR-NAs 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).

Conclusion
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.
Additional file 1. KEGG pathway analysis of genes differentially expressed circRNAs targeting.