- Primary research
- Open Access
Identification of circular RNA hsa_circ_0044556 and its effect on the progression of colorectal cancer
Cancer Cell International volume 20, Article number: 427 (2020)
Circular RNAs (circRNAs) are a novel class of noncoding RNAs. Increasing evidence indicates that circRNAs play an important role in the occurrence and development of tumors. However, the role of circRNA hsa_circ_0044556 in the progression of colorectal cancer (CRC) remains unclear.
First, we searched for differentially expressed circRNAs using a circRNA microarray in paired CRC and adjacent normal tissues. The circRNA hsa_circ_0044556 was screened out from the existing CRC circRNA microarray in the Gene Expression Omnibus database and our microarray. The clinical significance of hsa_circ_0044556 expression level in CRC patients was then investigated. Finally, the functions of the targets of this circRNA were determined in CRC cell lines.
Hsa_circ_0044556 was highly expressed in CRC patients and was positively correlated with tumor stage and lymph node metastasis. In CRC cell lines, the proliferation, migration, and invasion of cancer cells were inhibited by knocking down hsa_circ_0044556 expression.
Hsa_circ_0044556 promoted the progression of CRC. It is possible that hsa_circ_0044556 will become a novel biomarker or therapeutic target for CRC.
Colorectal cancer (CRC) is a malignant cancer that seriously endangers the health of humans. Currently, the global incidence rate of CRC ranks third among cancer-related diseases with up to 1.2 million new cases each year, and more than 0.6 million deaths are expected each year . From the perspective of the global incidence trend, the European and American regions are higher than the Asian and African regions, with the second place mortality rate (9.2%), and the developing countries are higher than the developed countries . In 2015, China's cancer statistics showed that the incidence and the mortality of CRC ranked fifth among cancer-related diseases in China. The incidence of CRC and the death toll are rising yearly while the population with this disease trends younger . Early diagnosis, accurate prognosis, and recurrence monitoring play important roles in cancer diagnosis and treatment. The 5-year survival rate of patients with advanced CRC is only 12%, while the 5-year survival rate of patients at the early stage can reach more than 90% . Therefore, early diagnosis can significantly improve the prognosis as well as the survival and quality of life of patients with CRC. On this basis, the in-depth study of the pathophysiology mechanisms underlying the occurrence and development of CRC will help us more comprehensively understand this cancer and open up new ideas and methods for the diagnosis and treatment of CRC. However, accurate and reliable prognosis and recurrence monitoring methods for CRC patients are still lacking in clinical practice, which makes it impossible to perform detailed posttreatment management and recurrence monitoring for patients [5, 6].
Circular RNAs (circRNAs) are a new class of noncoding RNAs (ncRNAs), which are single-stranded circular RNAs with no free 5′-end cap or 3′-end poly (A) tail. They are produced by alternative splicing of a specific pre-mRNA (pre-mRNA). Most circRNAs consist of exons and may also contain intergenic or noncoding regions [7,8,9,10]. CircRNAs did not attract much attention from researchers at first. For a long time, circRNAs were considered by-products of incorrect alternative splicing. Then in 2012, Salzman et al. discovered a large number of circRNAs using high-throughput sequencing technology . Now it is recognized that circRNAs are not as rare as previously thought. In contrast, circRNAs are highly stable in cells with high expression, and sometimes their expression levels are even 10 times higher than those of their homologous messenger RNAs. It has been found that circRNAs can act as a miRNA sponge bound to RNA-binding protein to exert biological functions.
In recent years, circRNAs have been considered to play an important role in tumor progression . CircRNAs are expected to become biomarkers for tumor diagnosis and prognosis due to their stable circular structure . Xu et al. found that hsa_circ_0001649 was expressed at abnormally low levels in intrahepatic cholangiocarcinoma tissues and could promote cell proliferation and tumor metastasis . Circ-ITCH can act as a sponge of miR-7 and miR-20a to inhibit the negative regulatory effect of the latter on the target gene ITCH, while ITCH has a tumor suppressor function by inhibition of the Wnt/β-catenin signaling pathway . However, the role of circRNA hsa_circ_0044556 in CRC remains unclear and deserves further investigation.
Materials and methods
In this study, tissue samples of 52 patients with CRC from the Third Xiangya Hospital of Central South University, China between May 2018 and December 2018 were collected. No patient received radiotherapy or chemotherapy before surgery, and the postoperative pathological diagnosis was adenocarcinoma (high, medium, and low differentiation). For each specimen, two copies of cancer tissues and adjacent normal tissues were collected. The surgical specimens were cryopreserved in liquid nitrogen immediately after excision. This study was approved by the Ethics Committee of the hospital, and all patients gave informed consent.
Human CRC cell lines (SW480, SW620, HCT116, HT29) and normal colonic epithelial cell lines (NCM460) were purchased from Wuhan Boster Biological Technology, Ltd. (Wuhan, China). SW480 and SW620 cells were cultured in L15 medium (Nanjing KeyGen Biotech Co., Ltd., Nanjing, China) containing 10% fetal bovine serum (FBS: Biological Industries Israel Beit-Haemek, Beit-Haemek, Israel). HCT116 and HT29 were cultured with McCoy's 5A medium (Nanjing KeyGen Biotech Co., Ltd.) containing 10% FBS. NCM460 was cultured in Dulbecco's Modified Eagle Medium (DMEM, Thermo Fisher Scientific™, Beijing China) containing 10% FBS. All cells were incubated at 37 °C in a 5% CO2 incubator.
siRNA-1, siRNA-2, siRNA-3, and negative-control siRNA were all designed by Suzhou Genepharma Co., Ltd., China. Cell transfection was done according to the manual of Lipofectamine 3000 (Invitrogen; Thermo Fisher Scientific, Inc.).
Tablet cloning assay
The tumor cell lines that were successfully transfected with siRNA-3 were harvested, counted, and adjusted to a cell concentration of 500 cells/ml, then were seeded in 6-well plates (1 × 103/well), and 3 duplicate wells were set for each cell line. The 6-well plates were incubated at 37 °C in a 5% CO2 incubator for 2 weeks. The cell culture medium was changed every 2 to 3 days. The 6-well plates were taken out after 2 weeks, and the medium was removed. Cells were then washed with phosphate-buffered saline (PBS) 2 times, dried naturally, and then fixed by adding 1 ml of 4% paraformaldehyde in each well. After cells were fixed at room temperature for 30 min, the paraformaldehyde solution was aspirated. Cells were dried and stained with 0.1% crystal violet for 30 min. After washing with PBS twice, the cells were photographed and counted.
Cell Counting Kit-8 (CCK-8) assay
SW480 and HCT116 (5 × 103/well) were mixed well in 100 µl normal medium containing 10% FBS, seeded into 96-well plates, and incubated at 37 °C in a 5% CO2 incubator. Each well received 10 μl of CCK-8 reagent at the specified time point (0, 24, 48, and 72 h). After incubation for another 3 h, the absorbance at 450 nm of each well was measured by an EnVision microplate reader (PerkinElmer, Inc., Waltham, MA, USA).
A marker was used to draw horizontal lines with the help of a ruler on the back of a 6-well plate. The lines were evenly drawn to 1 cm long, and there were at least 5 lines for each well. A total of 5 × 105 cells were inoculated into each well of a 6-well plate. Cells were transfected when they reached a density of approximately 70% confluence. When the cells grew to just cover the entire well, a 10 μl pipette tip was used to make a scratch, with the help of a ruler, perpendicular to the horizontal lines on the back of the plate. The detached cells were washed with PBS, and then serum-free medium was added to continue the culture for 72 h. Samples were taken at 0, 24, 48, and 72 h, the cells were photographed, and the scratch width was measured.
Migration and invasion assays
Migration and invasion assays were performed using a transwell chamber (8 μm, 24-well insert; Corning Incorporated, Corning, NY, USA). After 3 × 105 cells were mixed in 200 µl of serum-free medium, they were seeded in the upper chamber, 500 µl normal medium containing 20% FBS was added into the lower chamber, and then the plates were incubated for 48 h at 37 °C in a 5% CO2 incubator. The upper chamber and the lower chamber were washed twice with PBS. The cells in the upper chamber were cleaned with a cotton swab to remove cell debris. The chamber membrane was fixed with 4% paraformaldehyde for 30 min. The paraformaldehyde solution was aspirated and cells were dried, followed by 0.1% crystal violet stain for 30 min and two washes with PBS. Cells on the lower side of the chamber membrane were observed under an inverted fluorescence microscope (Olympus Corporation, Tokyo, Japan). The cells were photographed and cell numbers were calculated.
Reverse transcription and real-time quantitative polymerase chain reaction (RT-qPCR)
Total RNAs in tissues and cells were extracted using TRIzol reagent (Invitrogen; Thermo Fisher Scientific, Inc. Waltham, MA, USA). RT-qPCR was performed according to the manual of the Toyobo RT kit (Toyobo Life Science, Osaka, Japan) and Hieff® qPCR SYBR® Green Master Mix (Yason Biotech Co., Ltd.). CircRNA primers were designed by Beijing Tsingke Biological Technology Co., Ltd. (Hunan, China). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as the internal control, and the relative expression levels were calculated by the 2−ΔΔCt method.
CircRNA–miRNA–mRNA coexpression network construction
CircRNA–miRNA interactions were predicted using Arraystar’s homemade miRNA target prediction software (Rockville, MD, USA) based on TargetScan  and miRanda . The miRNAs were scored and sequenced using the miRNA support vector regression (mirSVR) algorithm to focus on target miRNAs . Therefore, for each circRNA, we identified the top 5 miRNAs in the miRVR score system to establish a top-5 circRNA–miRNA network (1 circRNA connecting to 5 miRNAs). To further predict the interactions between miRNA and mRNAs, the miRNA–mRNA overlapping set predicted by databases of miRDB , miRTarBase , and TargetScan were used to plot circRNA–miRNA–mRNA interaction networks in Cytoscape (version 3.4.0) .
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of circRNAs
The Database for Annotation, Visualization and Integrated Discovery (DAVID; https://david.ncifcrf.gov (Version6.7))  is an online bioinformatics database that integrates biological data and analysis tools. It provides a complete set of gene and protein functional annotation information for users to extract biological information. KEGG is a database resource used to categorize high-level functions and biological systems from large-scale molecular datasets generated by high-throughput experimental techniques . GO is an important bioinformatics tool for annotating genes and analyzing their biological processes . To analyze the functions of circRNAs, the online database of DAVID was used for biological analysis. P < 0.05 indicated statistical significance.
Data analysis: SPSS 19.0 software was used for statistical analysis (IBM Corp., Armonk, NY, USA)
The data were imaged using GraphPad Prism 6 software (GraphPad Software, Inc., La Jolla, CA, USA), and the data are expressed as the mean ± standard deviation. The differences between two groups were analyzed using the independent-sample t test, and intragroup differences were analyzed using one-way analysis of variance (ANOVA). P < 0.05 indicated that a difference was statistically significant.
Expression profiles of circRNAs in human CRC tissues
To study the expression profile of circRNAs in human CRC tissues, we used circRNA microarray technology to detect and analyze the expression of circRNAs in 3 pairs of CRC tissues and adjacent normal tissues. We drew a box plot showing the density distributions of all datasets after normalization and found that the distributions of log2 ratios were similar between all test samples (Fig. 1a). The unsupervised hierarchical clustering shows that circRNAs had different expression in CRC vs. adjacent normal tissues (Fig. 1b). The differential expression of these circRNAs was further confirmed in the volcano plot. With cutoffs of log2(fold change) ≥ 1 and P < 0.05 for differentially expressed circRNAs compared with adjacent normal tissues, 66 circRNAs in CRC tissues were upregulated, while 77 circRNAs were downregulated (Fig. 1c). Among the upregulated circRNAs, 52 were composed of exons, 9 were composed of introns, 1 was composed of intergenic regions, 2 were mixed circRNAs, and 2 were antisense circRNAs. Among the downregulated circRNAs, 67 were composed of exons, 4 were composed of introns, 1 was composed of intergenic regions, and 5 were mixed circRNAs (Additional file 1: Table S1).
Identification of the target circRNAs in CRC
The circRNAs with upregulated expression in our microarrays (Additional file 1: Table S1) were compared with the circRNAs with upregulated expression in the CRC circRNA microarray GSE126094  from the GEO database . A total of 20 circRNAs were found to be upregulated in both microarrays, and the top 10 co-upregulated circRNAs are listed in Table 1. The top two circRNAs (hsa_circ_0004104 and hsa_circ_0044556) with upregulated expression in both microarrays were selected for further experiments. To further verify whether hsa_circ_0004104 and hsa_circ_0044556 had consistent expression between tissues and the microarray, we extracted total RNA from 10 pairs of CRC tissues and adjacent normal tissues for RT-qPCR validation. The results showed that compared with the adjacent normal tissue, hsa_circ_0044556 expression was upregulated in CRC tissues, which was consistent with the microarray results (Additional file 2: Figure S1A). The expression of hsa_circ_0004104 in CRC tissues and paracancerous tissues was not significantly different, which was inconsistent with the microarray results (Additional file 2: Figure S1B).
Expression characteristics of hsa_circ_0044556 in CRC
To further validate the clinical significance of hsa_circ_0044556, 42 pairs of CRC tissues and adjacent normal tissues were collected. We confirmed that hsa_circ_0044556 was highly expressed in CRC tissues compared to adjacent normal tissues (Fig. 2a) by RT-qPCR. Hsa_circ_0044556 was significantly upregulated in 69.23% (36/52) of CRC tissues compared to adjacent normal tissues (Fig. 2b). Based on the expression level of hsa_circ_0044556, the diagnostic value of hsa_circ_0044556 in distinguishing CRC from adjacent normal tissues was calculated using the receiver operating characteristic (ROC) curve. The area under the ROC curve (AUC) was 0.7274 (P < 0.0001) (Fig. 2c).
To further investigate whether the high expression level of hsa_circ_0044556 in patients was related to clinicopathological parameters, Table 2 was developed. We can see from Table 2 that the expression level of hsa_circ_0044556 did not have a significant difference in terms of patient sex (P = 0.3548), age (P = 0.562), tumor size (P = 0.9865), or position (P = 0.8677), but a significant difference was observed for tumor stage (P = 0.0121) and lymph node metastasis (P = 0.0045) (Table 2).
Next, we used CRC cell lines for experiments. Compared to in the colonic epithelial cell line NCM460, hsa_circ_0044556 was highly expressed in the CRC cell lines of HCT116, SW480 and SW620 but was lower expressed in the CRC cell line HT29 (Fig. 2d). Because hsa_circ_0044556 was expressed at a relatively higher level in the two cell lines HCT116 and SW480, these two cell lines were chosen for further experiments.
Loop-forming validation and siRNA design of hsa_circ_0044556
By reviewing the human reference genome, we found that hsa_circ_0044556 (chr17: 48271490–48272189) was composed of exons 21–24 of the collagen type Ι alpha Ι (COL1A1) gene and was located on human chromosome 17q21.33. By RT-qPCR, Sanger sequencing, and anti-ribonuclease R digestion in HCT116 and SW480 cells, we verified that hsa_circ_0044556 had a circular structure (Fig. 3a–c). To evaluate the biological effects of hsa_circ_0044556 in CRC cells, we constructed three siRNAs covering the back-splicing region of hsa_circ_0044556. The results showed that the expression level of hsa_circ_0044556 in CRC cells transfected with siRNA-3 significantly decreased, while COL1A1 mRNA expression level did not change significantly (Fig. 3d–g). Therefore, hsa_circ_0044556 was chosen for further experiments.
Silencing hsa_circ_0044556 inhibits the proliferation, migration, and invasion of CRC cells
The results of the tablet cloning assay and CCK-8 assay showed that silencing of hsa_circ_0044556 lowered the proliferation ability of HCT116 and SW480 cells (Fig. 4a, b). The results of the scratch test and transwell experiment showed that knockdown of hsa_circ_0044556 inhibited the migration and invasion capability of HCT116 and SW480 cells (Fig. 4c–e).
CircRNA–miRNA–mRNA coexpression network of hsa_circ_0044556
We hypothesized that hsa_circ_0044556 acts as a miRNA “sponge” to regulate circRNA–miRNA–mRNA networks. Using miRNA target prediction software, five miRNAs with the highest miRVR scores that might be bound by differentially expressed circRNAs were identified by TargetScan and miRanda (Additional file 1: Table S1). The molecular interaction between hsa_circ_0044556 and the 5 miRNA targets is depicted in Fig. 5a. Then, miRDB, miRTarBase, and TargetScan were used to predict the target genes that those five miRNAs might bind to. A total of 107 predicted genes (Additional file 3: Table S2) from the 3 databases were selected as the potential target genes of hsa_circ_0044556. Cytoscape analysis of the circRNA–miRNA–mRNA interaction network of hsa_circ_0044556 revealed that hsa-mir-214-3p and hsa-mir-761 exhibited the most complex interaction network, followed by hsa-mir-194-3p, hsa-mir-412-3p, and hsa-mir-362-5p (Fig. 5b). Since predicted target miRNAs only hsa-mir-214-3p proved to be downregulated in cancer progression in both ENCORI  and UALCAN  (Additional file 4: Figure S2A, B and C). Through further verification, we found that down-regulation of the expression level of hsa_circ_0044556 in HCT116 and SW480 cells increased the expression level of hsa-mir-214-3p (Additional file 4: Figure S2D).
Bioinformatic analysis of the predicted genes in hsa_circ_0044556
Gene ontology (GO) analysis was performed on hsa_circ_0044556, and the functional roles of the top 10 most enriched target genes were investigated from the perspective of biological processes (Fig. 6a). The results showed that hsa_circ_0044556 had a strong correlation with the cell cycle, responses to steroid hormone stimulation, and angiogenesis. KEGG analysis of hsa_circ_0044556 indicated that its top nine enriched pathways included prostate cancer, the cancer pathway, and the ErbB signaling pathway (Fig. 6b). Seven target genes were enriched in the cancer pathway. These data suggested that hsa_circ_0044556 may play important roles in the malignant behavior of cancer by regulating the expression of target genes involved in these pathways.
In recent years, with the rapid development and extensive application of RNA sequencing technology, circRNAs have become a hotspot in the field of RNA research. Researchers have found that many exon transcripts can form circRNAs through nonlinear reverse splicing or gene rearrangement. Moreover, they account for a large percentage of all spliced transcripts . CircRNAs may come from introns or exons . In mammals, exon-constituted circRNAs have two mechanisms of loop formation: lariat-driven circularization and intron-pairing-driven circularization [7, 29, 30]. A covalently closed circular structure without a 5′-to-3′ polarity or poly(A) tail is then formed by reverse splicing of a typical splice.
More and more studies have investigated the potential functions of circRNAs in various diseases, such as nervous system diseases, cardiovascular diseases, and cancers [31,32,33,34]. Some abnormally expressed circRNAs have been associated with the tumor development, invasion, metastasis, or prognosis of patients [35,36,37,38,39,40,41]. In mammalian cells, compared with other ncRNAs, such as miRNAs and long noncoding RNAs, circRNAs have highly conserved sequences and high stability . These features might let circRNAs become ideal biomarkers and potential therapeutic targets for disease diagnosis.
In this study, high-throughput circRNA microarrays were used to study the expression of circRNAs in human CRC. The results showed that the expression of circRNAs in CRC tissues (n = 3) was significantly different from that in adjacent normal tissues (n = 3) (Fig. 1). Compared with the adjacent normal tissues, our microarray data showed that 66 circRNAs were significantly upregulated while 77 circRNAs were significantly downregulated in CRC tissues (Additional file 1: Table S1). We verified the expression levels of two circRNAs (hsa_circ_0004104 and hsa_circ_0044556) in 10 pairs of CRC and adjacent normal tissue samples, and only hsa_circ_0044556 was determined to be upregulated in CRC tissues, and with statistical significance (P = 0.0304) (S Additional file 2: Figure S1A). The expression level of hsa_circ_0004104 was inconsistent with the microarray results, and there was no statistical significance (P = 0.9208) (Additional file 2: Figure S1B). The above results indicate that validation of differently expressed circRNAs in microarray analysis is an important step in such a screening study. In addition, to verify the results above, we expanded the sample size to detect the expression of hsa_circ_0044556 in other tissue samples (n = 52). The results showed that hsa_circ_0044556 was significantly upregulated in 69.23% (36/52) of CRC tissues, with an average increase of 6.65-fold compared to the adjacent normal tissues. ROC analysis showed that hsa_circ_0044556 level had relatively high sensitivity and specificity, with an AUC of 0.7274 (Fig. 2). More importantly, considering the clinical pathological factors, we found that the high expression level of hsa_circ_0044556 in CRC was closely related to tumor stage and lymph node metastasis (Table 2), important factors in evaluating the prognosis of CRC. These results indicate that hsa_circ_0044556 might be involved in the progression and metastasis of CRC and could be used as a potential biomarker and a new therapeutic target for CRC.
One of the most important things about circRNAs is that they act as miRNA sponges. Certain specific circRNAs can bind and negatively regulate miRNAs involved in the competitive endogenous RNA (ceRNA) network, thereby regulating linear RNA transcription and protein production. Thomas et al. found that circRNA ciRS-7 can strongly inhibit miR-7 activity, leading to increased miR-7 target expression levels [40, 43]. Other functions may include gene expression regulation at the transcriptional or posttranscriptional level , and even encoding proteins [45, 46]. In this study, to further understand the biological functions of hsa_circ_0044556, 5 miRNAs with the highest mirSVR scores were identifed for each differentially expressed circRNA using miRNA target-prediction sofware (namely, hsa-mir-214-3p, hsa-mir-761, hsa-mir-194-3p, hsa-mir-412-3p, and hsa-mir-362-5p) (Additional file 1: Table S1). We used TargetScan, miRDB, and miRTarBase to predict the hsa_circ_0044556-miRNA–mRNA network. In essence, this network diagram shows a cellular RNA network with hsa_circ_0044556 interacting with 5 miRNA nodes and 107 target genes (Fig. 5b). Through the retrieval of ENCORI and UALCAN, it was found that hsa-mir-214-3p had low expression in colorectal adenocarcinoma. Furthermore, it have been found that hsa-mir-214-3p was significantly reduced in epithelial ovarian cancer cells and could affect epithelial ovarian cancer cell proliferation, invasion, increasing cisplatin chemosensitivity and inhibiting in vivo tumor growth proliferation by binding X-inactive specific transcript (XIST) . It have also been reported that the downregulated hsa-miR-145-5p and hsa-mir-214-3p may modulate the expression of both EMT and NGAL/MMP-9 pathways . In CRC lines, we also confirmed that silencing the expression level of hsa_circ_0044556 increased the expression level of hsa-mir-214-3p. Therefore, the decreased expression and inhibited function of hsa-mir-214-3p in CRC further support our hypothesis that hsa_circ_0044556 functions as a miRNA sponge to regulate the hsa_circ_0044556- hsa-mir-214-3p-mRNA network.
We found that a large number of mRNAs may participate in the above hsa_circ_0044556-miRNA–mRNA network, such as ARL2, MAPK1, and PTEN. Therefore, GO and KEGG pathway analysis was performed to detect the functions of these potential target genes. The results of GO enrichment analysis (Fig. 6a) showed that the target genes of hsa_circ_0044556 participated in the regulation of the cell cycle and angiogenesis, indicating that the regulation of these genes in the occurrence and development of CRC has importance in cell responses. The cancer pathway and the ErbB signaling pathway, two KEGG pathways correlated with hsa_circ_0044556 expression (Fig. 6b), may be related to the proliferation, migration, and invasion of CRC cells. Therefore, we speculate that the hsa_circ_0044556-miRNA–mRNA axis is a possible mechanism that promotes the development of CRC, and it is worthwhile to further study the overexpression of hsa_circ_0044556 as an inhibitor of miRNA and its possible mechanism of action.
In summary, this study revealed the expression profile of circRNAs in CRC tissues and demonstrated the abnormal expression of circRNAs in CRC. This study confirmed the significance of the upregulation of hsa_circ_0044556 and analyzed the relationship between hsa_circ_0044556 and the clinicopathological features of CRC patients, suggesting its potential role in the development and progression of CRC and its potential application as a CRC diagnostic biomarker. In the future, it will be necessary to explore the molecular mechanism of hsa_circ_0044556 as a miRNA sponge regulating the development and progression of CRC.
Hsa_circ_0044556 promoted the progression of CRC. It is possible that hsa_circ_0044556 will become a new biomarker or therapeutic target for CRC.
Data availability statement
All datasets presented in this study are included in the article/additional files.
ADP ribosylation factor-like protein 2
Analysis of Variance
Collagen type 1 alpha 1
The Database for Annotation Visualization and Integrated Discovery
Dulbecco's Modified Eagle Medium
The Encyclopedia of RNA Interactomes
Fetal bovine serum
Glyceraldehyde phosphate dehydrogenase
Gene expression omnibus
Itchy E3 ubiquitin protein ligase
Kyoto Encyclopedia of Genes and Genomes
Mitogen activated protein kinase 1
Receiver operating characteristic
Phosphate buffered saline
Polymerase chain reaction
Phosphatase and tensin homolog
Reverse transcription and real-time quantitative PCR
Small interfering RNA
Short hairpin RNA
Statistical Product and Service Solutions
The Cancer Genome Atlas
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.
Bray F, Ferlay J, Soerjomataram I, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;680:1–31.
Chen W, Zheng R, Baade RD, Zhang S, Zeng H, Bray F, Jemal A, Yu XQ, He J. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–32.
American Cancer Society. Colorectal Cancer Facts and Figures 2011–2013. Atlanta: American Cancer Society; 2011.
Imperiale TF, Ransohoff DF, Itzkowitz SH, et al. Fecal DNA versus fecal occult blood for colorectal-cancer screening in an average risk population. N Engl J Med. 2004;351:2704–14.
Lieberman DA. Screening, surveillance, and prevention of colorectal cancer. Gastrointest Endosc Clin N Am. 2008;18:595–605.
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):1–12.
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.
Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J, et al. Circular RNAs are abundant, converved, and associated with ALU repeats. RNA. 2013;19(2):141–57.
Danan M, Schwartz S, Edelheit S, Sorek R. Transcriptome-wide discovery of circular RNAs in Archaea. Nucleic Acids Res. 2012;40(7):3131–42.
Du WW, Fang L, Yang W, Wu N, Awan FM, Yang Z, Yang BB. Induction of tumor apoptosis through a circular RNA enhancing Foxo3 activity. Cell Death Differ. 2017;24:357–e370.
Wang J, Li H. CircRNA circ_0067934 silencing inhibits the proliferation, migration and invasion of NSCLC cells and correlates with unfavorable prognosis in NSCLC. Eur Rev Med Pharmacol Sci. 2018;22:3053–60.
Xu Y, Yao Y, Zhong X, Leng K, Qin W, Qu L, Cui Y, Jiang X. Downregulated circular RNA hsa_circ_0001649 regulates proliferation, migration and invasion in cholangiocarcinoma cells. Biochem Biophys Res Commun. 2018;496:455–61.
Huang G, Zhu H, Shi Y, et al. Circ-ITCH plays an inhibitory role in colorectal cancer by regulating the Wnt/β-catenin pathway. PLoS ONE. 2015;10(6):e0131225.
Agarwal V, Bell GW, Nam J, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. eLife. 2015;4:e05005.
Pasquinelli AE. MicroRNAs and their targets: recognition, regulation and an emerging reciprocal relationship. Nat Rev Genet. 2012;13:271–82.
Betel D, Koppal A, Agius P, Sander C, Leslie C. Comprehensive modeling of microRNA targets predicts functional non-conserved and non-canonical sites. Genome Biol. 2010;11:R90.
Nathan Wong and Xiaowei Wang. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015;43(D1):D146–152.
Chou C-H, Shrestha S, Yang C-D, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46(D1):D296–D302.
Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27:431–2.
Huang DW, Sherman BT, Tan Q, Collins JR, Alvord WG, Roayaei J, Stephens R, Baseler MW, Lane HC, Lempicki RA. The DAVID Gene Functional Classification Tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol. 2007;8:R183.
Kanehisa M. The KEGG database. Novartis Found Symp. 2002;247:91–252.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.
Chen Z, Ren R, Wan D, Wang Y, et al. Hsa_circ_101555 functions as a competing endogenous RNA of miR-597-5p to promote colorectal cancer progression. Oncogene. 2019;38(32):6017–34 .
Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.
Li JH, et al. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42:D92–D9797.
Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Rodriguez IP, Chakravarthi BVSK, Varambally S. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19(8):649–58.
Zhang Y, et al. Circular intronic long noncoding RNAs. Mol Cell. 2013;51:792–806.
Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol. 2014;32:453–61.
Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J, Marzluff WF, Sharpless NE. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA. 2013;19:141–57.
Li J, et al. Circular RNAs in cancer: novel insights into origins, properties, functions and implications. Am J Cancer Res. 2015;5:472–80.
Wang Y, et al. Circular RNAs in human cancer. Mol Cancer. 2017;16:25.
Lu D, Xu AD. Mini review: circular RNAs as potential clinical biomarkers for disorders in the central nervous system. Front Genet. 2016;7:53.
Wang K, et al. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur Heart J. 2016;37:2602–11.
Xia W, et al. Circular RNA has_circ_0067934 is upregulated in esophageal squamous cell carcinoma and promoted proliferation. Sci Rep. 2016;6:35576.
Li P, et al. Using circular RNA as a novel type of biomarker in the screening of gastric cancer. Clin Chim Acta Int J Clin Chem. 2015;444:132–6.
Qin M, et al. Hsa_circ_0001649: a circular RNA and potential novel biomarker for hepatocellular carcinoma. Cancer Biomarkers Sect A Dis Mark. 2016;16:161–9.
Shang X, et al. Comprehensive circular RNA profling reveals Tat hsa_circ_0005075, a new circular RNA biomarker, is involved in hepatocellular crcinoma development. Medicine. 2016;95:e3811.
Wang X, et al. Decreased expression of hsa_circ_001988 in colorectal cancer and its clinical signifcances. Int J Clin Exp Pathol. 2015;8:16020–5.
Xuan L, et al. CircularRNA: a novel biomarker for progressive laryngeal cancer. Am J Transl Res. 2016;8:932–9.
Song X, et al. Circular RNA profile in gliomas revealed by identifcation tool UROBORUS. Nucleic Acids Res. 2016;44:e87.
Memczak S, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495:333–8.
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495:384–8.
Hansen TB, Wiklund ED, Bramsen JB, et al. MiRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA. EMBO J. 2011;30:4414–22.
Chen CY, Sarnow P. Initiation of protein synthesis by the eukaryotic translational apparatus on circular RNAs. Science. 1995;268:415–7.
Perriman R, Ares M Jr. Circular mRNA can direct translation of extremely long repeating-sequence proteins in vivo. RNA. 1998;4:1047–54.
Wang C, Qi S, Xie C, Li C, Wang P, Liu D. Upregulation of long non-coding RNA XIST has anticancer effects on epithelial ovarian cancer cells through inverse downregulation of hsa-miR-214-3p. J Gynecol Oncol. 2018;29(6):e99. https://doi.org/10.3802/jgo.2018.29.e99.
Falzone L, Candido S, Salemi R, et al. Computational identification of microRNAs associated to both epithelial to mesenchymal transition and NGAL/MMP-9 pathways in bladder cancer. Oncotarget. 2016;7(45):72758–66. https://doi.org/10.18632/oncotarget.11805.
The authors wish to acknowledge the technical assistance of KangChen Bio-tech, Shanghai, China for Microarray experiments.
This work was supported by the Hunan Province Technological Innovation Guidance Program Foundation (No. S2017SFYLJS0274) and The Hunan Province Natural Science Foundation (No.2016JJ2155).
Ethics approval and consent to participate
The study which including human and animal approval were acquired from the Ethics Committee of The Third Xiangya Hospital of Central South University [Lot 2016-S086].
Consent for publication
We assure that the material is original and it has not been published elsewhere yet.
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Results of circRNAs expression microarray in CRC.
The expression levels of has_circ_0044556 and has_circ_0004104 in CRC were preliminarily verified.
Hsa_circ_0044556 common target genes in miRDB, miRTarBase and TargetScan.
The expression levels of has-mir-214-3p in the CRC tissues and adjacent normal tissues.
About this article
Cite this article
Jing, L., Wu, J., Tang, X. et al. Identification of circular RNA hsa_circ_0044556 and its effect on the progression of colorectal cancer. Cancer Cell Int 20, 427 (2020). https://doi.org/10.1186/s12935-020-01523-1
- Noncoding RNA
- Circular RNA
- Colorectal cancer