Circular RNA, hsa_circRNA_102049, promotes colorectal cancer cell migration and invasion via binding and suppressing miRNA-455-3p.

Colorectal cancer (CRC) is the second most prevalent malignant gastrointestinal tumor type worldwide, displaying poor prognosis. Accumulating studies have reported the significance of circular RNAs (circRNAs) and microRNAs (miRNAs) in CRC carcinogenesis and development. At present, the functions and mechanisms of action underlying the circular RNA, hsa_circRNA_102049, in CRC are not completely understood. The present study aimed to establish the involvement of hsa_circRNA_102049 in CRC, as well as the associated mechanisms. The expression levels of hsa_circRNA_102049 and miRNA-455-3p were measured in CRC cell lines and tissues via reverse transcription-quantitative PCR. CRC progression was evaluated by performing Cell Counting Kit-8, flow cytometry, wound healing and Transwell invasion assays. The results demonstrated that hsa_circRNA_102049 was highly expressed in both CRC tissues and cell lines, which was associated with enhanced CRC cell proliferation, migration and invasion. Furthermore, miR-455-3p expression was downregulated in CRC cells and served as a target of has_circRNA_102049, which was validated by performing the dual luciferase reporter assay. hsa_circRNA_102049 knockdown significantly increased miR-455-3p expression, which was significantly reversed by co-transfection with the miR-455-3p inhibitor. Notably, miRNA-455-3p overexpression alleviated hsa_circRNA_102049-mediated induction of CRC cell proliferation, migration and invasion. The present study clearly demonstrated that miRNA-455-3p was a target of hsa_circRNA_102049. Moreover, the results indicated that the circular RNA, hsa_circRNA_102049, may function as a tumor promoter in CRC via directly sponging miRNA-455-3p. Copyright: © Zhu et al.

Introduction

Colorectal cancer (CRC) is the third most common malignancy and the fourth leading cause of cancer-related death worldwide (1-3). CRC is responsible for 1 million morbidities and half a million mortalities annually worldwide (4). Owing to its high rate of metastasis and aggressive malignancy, CRC has become both a public health and political issue on a global scale (5). In China, the morbidity and mortality rates of CRC continue to increase every year owing to changes in lifestyle, including dietary structure and living environments, making CRC one of the most prevalent malignancies with the fastest growing incidence. Despite the continuous development of surgical techniques and adjuvant therapy, the survival rate of patients with CRC has not significantly improved, displaying a ~50% death rate due to local recurrence or metastasis after surgery (6). The occurrence and progression of CRC is a multi-factorial, multi-step and multi-stage process closely related to environment, genetics, diet and inflammation (7). Moreover, aberrant expression of multiple signals that regulate cell proliferation and differentiation serves a significant role in tumor occurrence and may be effectively used as a predictor of prognosis (8,9). However, the precise mechanisms underlying the progression of CRC require further investigation.

MicroRNAs (miRNAs/miRs) are short non-coding RNAs (22 nucleotides in length) that are widely distributed in animals, plants, viruses and single-celled organisms. miRNAs induce transcriptional and post-transcriptional inhibition by complementing the 3'-untranslated regions (3'UTRs) of target mRNAs (10,11). miRNAs serve critical roles in diverse biological processes, including the regulation of cell proliferation, differentiation and apoptosis, in multiple types of cancer (12-15). Several studies have demonstrated that miRNAs display specific expression patterns in tumors, serving as tumor suppressors or oncogenes with roles in tumor development via negative regulation of specific genes (16-18). miR-455-3p has been reported to be expressed at low levels and serve a tumor suppressor role in numerous types of cancer, including prostate cancer, breast cancer, pancreatic cancer and melanoma (19-22). Zheng et al (23) reported that miR-455-3p overexpression inhibits tumor cell proliferation and induces apoptosis, whereas miR-455-3p knockdown displays the opposite effects in human HCT116 colon cancer cells in vitro. However, the expression patterns and precise roles of miRNA-455-3p in human CRC are not completely understood.

Circular RNAs (circRNAs), a type of cyclic non-coding RNA, have been extensively identified in various eukaryotic cells and display high expression in tumor tissues. circRNAs have a highly stable covalent closed continuous loop structure and covalently combine with 3' and 5' ends after reverse splicing (24-27). circRNAs participate in tumorigenesis via regulation of several biological processes, including survival, proliferation, invasion and differentiation (28-30). In view of their functions as miRNA sponges, RNA-binding protein sequestering agents and transcriptional regulators, circRNAs have become a focus of competing endogenous RNA (ceRNA) research (31,32). The alias of hsa_circRNA_102049 is hsa_circ_0043278 and the gene symbol is transcriptional adaptor 2A (TADA2A; www.circbase.org/cgi-bin/simplesearch.cgi) (33). Previously, Zhang et al (34) demonstrated the involvement of TADA2A in HCT116 colon cancer cells. However, the specific functions and mechanisms of action underlying hsa_circRNA_102049 in human CRC cells and tissues require further investigation. With the aid of various experimental methods and microRNA Target Prediction Database (miRDB) bioinformatics software, the present study explored the molecular mechanisms underlying the functional interplay between hsa_circRNA_102049 and miR-455-3p, as well as how hsa_circRNA_102049 loss-of-function influenced miR-455-3p expression and colorectal cancer cell migration and invasion.

Materials and methods

Patients and tumor tissues

Human CRC and adjacent non-tumor tissues, located 10 cm from the tumor margin, were acquired from 68 patients (45 male patients and 23 female patients, mean age, 52.26±9.99 years; age range, 32-78 years) with histologically diagnosed CRC who underwent curative surgery with no prior anticancer treatment at Yiwu Central Hospital (Yiwu, China) between January 2017 and December 2017. All tissue samples were pathologically confirmed (WC and YZ) and immediately snap-frozen in liquid nitrogen and stored at -80˚C until further analysis. The exclusion criteria were as follows: i) CRC related to a genetic syndrome; and ii) previous anticancer treatments (radiotherapy, chemotherapy or pharmacotherapy). Patient characteristics, including sex, age, tumor grade, TNM staging (35), lymphatic invasion and histological classification, were extracted from patient medical records (Table I). The present study was approved by the Clinical Ethical Committee of Yiwu Central Hospital of Zhejiang Province (approval no. 20170112). Written informed consent was obtained from all patients prior to initiation of the study.

Table I

Associations between hsa_circRNA_102049 or miR-455-3p and the clinicopathological characteristics of 68 patients with colorectal cancer.

		hsa_circRNA_102049		miR-455-3p
Characteristic	Patient (n)	Expression	P-value	Expression	P-value
Sex			0.60		0.49
Male	45	0.78±0.21		0.54±0.16
Female	23	0.75±0.24		0.51±0.18
Age			0.50		0.23
≥50	38	0.79±0.18		0.58±0.11
<50	30	0.76±0.18		0.55±0.09
Differentiation grade			P<0.05		P<0.05
Low	23	0.95±0.16		0.47±0.12
Middle	20	0.81±0.11		0.55±0.12
High	25	0.76±0.08		0.68±0.24
TNM stage			P<0.05		P<0.05
I-II	32	0.65±0.12		0.64±0.20
III-IV	36	0.89±0.20		0.50±0.10
Lymphatic migration and invasion			P<0.05		P<0.05
Yes	40	0.92±0.24		0.43±0.06
No	28	0.72±0.14		0.53±0.11
Histological classification			P<0.05		P<0.05
Adenocarcinoma	27	0.75±0.19		0.66±0.12
Carcinoma muciparum	23	0.87±0.15		0.59±0.11
Undifferentiated carcinoma	18	1.00±0.19		0.51±0.09

Data are presented as the mean ± SD. circRNA, circular RNA; miR, microRNA.

Cell culture

A normal human intestinal epithelium cell line (NCM460), CRC cell lines (SW480, LOVO, HT29, DLD-1, SW620 and HCT116) and the 293T cell line were purchased from The Cell Bank of Type Culture Collection of The Chinese Academy of Science. All cell lines were tested and authenticated via DNA typing. NCM460, SW480, SW620, HCT116 and 293T cells were cultured in DMEM (Sigma-Aldrich; Merck KGaA) supplemented with 10% FBS (Gibco; Thermo Fisher Scientific, Inc.), 1% penicillin/streptomycin (Gibco; Thermo Fisher Scientific, Inc.) at 37˚C under 5% CO2 in a humidified atmosphere. LOVO, HT29 and DLD-1 cells were cultured in RPMI-1640 (Invitrogen; Thermo Fisher Scientific, Inc.) supplemented with 10% FBS and 1% penicillin/streptomycin at 37˚C with 5% CO2. All cells were incubated in a humidified chamber 37˚C with at 5% CO2. After incubation for 72 h, when cells reached 60-70%, subsequent experiments were conducted.

Reverse transcription-quantitative PCR (RT-qPCR)

Total RNA was extracted from cells or tissues using TRIzol® reagent (Invitrogen; Thermo Fisher Scientific, Inc.) according to the manufacturer's protocol. Total RNA was purified with RNase-free DNase treatment (Promega Corporation) following the manufacturer's protocols. miRNA and mRNA were reverse transcribed into cDNA using the Mir-X™ miRNA First-Strand Synthesis kit (Takara Bio, Inc.) and the PrimeScript™ RT reagent kit (Takara Bio, Inc.), respectively and the temperature protocol was 21˚C for 10 min, 48˚C for 50 min and 90˚C for 2 min. Subsequently, qPCR was performed using an optimal system (total volume, 25 µl) consisting of Mg2+ (1.00 µl), 10X PCR Buffer (2.00 µl), 0.5 µl dNTPs, 0.5 U Taq DNA polymerase (Thermo Fisher Scientific, Inc.), 0.6 µl forward and reverse primers, 2 µl template DNA and ddH2O. For hsa_circRNA_102049 expression, qPCR was performed using a PowerUp™ SYBR™ Green Master Mix (Thermo Fisher Scientific, Inc.). For hsa-miR-455-3p, qPCR was carried out by using the miRNA-specific TaqMan MiRNA Assay Kit (cat. no. TAP02280; Xinhai Gene Testing Co., Ltd.; https://www.haigene.cn/index.php?a=show&c=index&catid=75&id=245&m=content) as recommended by the manufacturers. The fluorophore used in this study was carboxyfluorescein. The following thermocycling conditions were used for qPCR: Pre-denaturation at 94˚C for 5 min; followed by 30 cycles of denaturation at 94˚C for 30 sec, annealing at 54.5˚C for 30 sec and extension at 72˚C for 30 sec; followed by final extension at 72˚C for 10 min. All primers were synthesized by Guangzhou RiboBio Co., Ltd. The following primers were used for qPCR: GAPDH forward, 5'-TATGATGATATCAAGAGGGTAGT-3' and reverse, 5'-TGTATCCAAACTCATTGTCATAC-3'; hsa_circRNA_102049 forward, 5'-AATGTGCACCAAGACCAAGG-3' and reverse, 5'-CCAAAGCCACAGTCCATCAC-3'; U6 forward, 5'-CTCGCTTCGGCAGCACATA-3' and reverse, 5'-AACGATTCACGAATTTGCGT-3'; hsa-miR-455-3p forward, 5'-GCAGUCCACGUGGGCAUAUACAC-3' and reverse, 5'-GCAGUCCAUGGGTGCAUAUA CAC-3'. miRNA and mRNA expression levels were quantified using the 2-∆∆Cq method (36) and normalized to the internal reference genes U6 and GAPDH, respectively. A relative fold change of 2.0 was considered significant. RT-qPCR was performed in triplicate.

RNA interference and transfection assays

The following cell groups were established: i) Control [SW480 cells transfected with the pcDNA3.1 (+) circRNA mini vector (cat. no. 60648; Addgene, Inc.)]; ii) hsa_circRNA_102049 [SW480 cells transfected with pcDNA-hsa_circRNA_102049 cDNA (accession no. NM_001488; designed and synthesized by Shanghai Genepharma Co., Ltd.)]; iii) negative control (NC) small interfering (si)RNA (HCT116 cells transfected with scrambled siRNA as the NC siRNA); and iv) hsa_circRNA_102049 siRNA (HCT116 cells transfected with hsa_circRNA_102049 siRNA). Small interfering RNA (si-RNA) targeting to the junction region of hsa_circ_102049 was induced for the loss-of-function study (Shanghai GenePharma Co., Ltd.). The sequence of the hsa_circRNA_102049 siRNA was 5'-UCUGAAGUAGUGAAAUGGAAU-3' and the sequence of NC siRNA was 5'-AGACUUCUAGUGAAAUGGAAU-3'.

To further investigate the interactions between hsa_circRNA_102049 and miR-455-3p, SW480 and HCT116 cells were transfected with miR-455-3p mimic (5'-GCAGUCCAUGGGCAUAUACAC-3'), inhibitor (5'-GUGUAUAUGCCCAUGGACUGC-3'), NC mimic (sense, 5'-UUUGUACUACACAAAAGUACUG-3' and antisense, 5'-CAGUACUUUUGUGUAGUACAAA-3') or NC inhibitor (5'-CAGUACUUUUGUGUAGUACAAA-3') (all purchased from Shanghai GenePharma Co., Ltd.) to overexpress or knockdown miR-455-3p expression, respectively. The following cell groups were established: i) hsa_circRNA_102049 + NC mimic (SW480 cells co-transfected with pcDNA-hsa_circRNA_102049 and NC mimic); ii) hsa_circRNA_102049 + miR-455-3p mimic (SW480 cells co-transfected with pcDNA-hsa_circRNA_102049 and miR-455-3p mimic); iii) hsa_circRNA_102049 siRNA + NC inhibitor (HCT116 cells co-transfected with hsa_circRNA_102049 siRNA and NC inhibitor); and iv) hsa_circRNA_102049 siRNA + miR-455-3p inhibitor (HCT116 cells co-transfected with hsa_circRNA_102049 siRNA and miR-455-3p inhibitor).

Cells (3x104 cells/well) were placed into 96-well plates and then transfected with 50 nM hsa_circRNA_102049 siRNA, 50 nM NC siRNA, 100 nM miR-455-3p inhibitor, 100 nM NC inhibitor, 100 nM miR-455-3p mimic, 100 nM NC mimic, 50 ng pcDNA-hsa_circRNA_102049 or 50 ng empty vector using Lipofectamine® 2000 Transfection reagent (Invitrogen; Thermo Fisher Scientific, Inc.) for 6 h at 37˚C according to the manufacturer's protocol. Cells were harvested after 48 h transfection for subsequent experimentation.

Dual luciferase reporter assay

The microRNA database miRDB (version 2.0; https://mirdb.org/) and TargetScan (version no. 6.2; https://www.targetscan.org) were used to determine the target gene of miR-455-3p and to verify whether hsa_circRNA_102049 is a direct target gene of miR-455-3p. To further explore whether hsa_circRNA_102049 targeted miR-455-3p, the 3'-UTR of hsa_circRNA_102049 was cloned into the pGL3 luciferase reporter vector (Promega Corporation). After reaching 70-80% confluence, 293T cells were co-transfected with 50 nmol/l miR-455-3p mimic or NC mimic and hsa_circRNA_102049-wild-type (WT) (100 ng/well) or hsa_circRNA_102049-mutant (Mut) (100 ng/well) luciferase reporter plasmid for 6 h at 37˚C using the Lipofectamine® 2000 Transfection reagent (Invitrogen; Thermo Fisher Scientific, Inc.). The luciferase reporter plasmid contained two luciferase reporter genes: i) Trepang kidney luciferase reporter gene; and ii) pRL-TK vector (Promega Corporation) provides constitutive expression of Renilla luciferase, which was co-transfected with the firefly luciferase reporter vector as an internal control. The 3'-UTR of hsa_circRNA_102049 was cloned directly downstream of the trepang kidney luciferase reporter gene. In cases where the 3'-UTR was recognized by miR-455-3p, the expression levels of genes upstream of this region were suppressed, and the ratio between the genes and internal reference was altered. At 24 h post-transfection, luciferase activities were determined using the Dual Luciferase Reporter Assay System (Promega Corporation) according to the manufacturer's protocol. The ratio of firefly luciferase activity to Renilla luciferase activity was calculated as the normalized luciferase activity.

Functional analyses of hsa_circRNA_102049 and miR-455-3p in vitro

To assess the functions of hsa_circRNA_102049 and miR-455-3p in vitro, cell proliferation, cell cycle, apoptosis, migration and invasion assays were conducted as previously described. Cell proliferation was assessed by performing Cell Counting Kit-8 (CCK-8; Beyotime Institute of Biotechnology) and colony formation assays (33,37). Briefly, the cells (5x103 cells/well) were seeded into 24-well plates. Subsequently, 10 µl CCK-8 was added to each well at 48 h after transfection and the cells were further incubated for 4 h at 37˚C. The optical density was measured at 490 nm using an iMARK plate reader (Bio-Rad Laboratories, Inc.).

For colony formation assay, cells (2x103 cells/well) were seeded into 6-well plates and incubated at 37˚C with 5% CO2 for 14 days. To visualize colonies, cells were stained with Giemsa solution (Sigma-Aldrich, Merck KGaA) for 20 min at room temperature. The number of colonies each containing >25 cells was determined for each dish and the number of colonies was counted using a light microscope (magnification, x100).

Flow cytometry was performed to assess the cell cycle distribution and cell apoptosis (38). Transfected cells (1x106 cells/ml) were collected and stained with 50 µg propidium iodide (PI) (Sigma-Aldrich; Merck KGaA) containing RNAase (Sigma-Aldrich; Merck KGaA) in the dark for 30 min at 37˚C and filtered with a 100-mesh nylon filter. Cell cycle distribution was measured using a BD FACSCalibur™ flow cytometry (BD Biosciences) by recording the red fluorescence at 488 nm. For cell apoptosis, Annexin-V-fluorescein isothiocyanate (FITC) apoptosis detection kit (cat. no. ab14085; Abcam) was applied. The cells were cultured at 37˚C with 5% CO2 for 48 h, cells at a density of 1x104 cells/well were stained with 5 µl Annexin V-FITC and 5 µl PI at room temperature in the dark for 15 min. Cell apoptosis was determined using a BD FACSCalibur™ flow cytometer (BD Biosciences) at 488 nm. The flow cytometry data were analyzed with CellQuest Pro software (version 3.3; BD Biosciences).

Cell migration was determined by conducting the wound healing assay (33). Cells (2x103 cells/well) were seeded into 6-well plates and incubated at 37˚C with 5% CO2. After reaching 70-80% confluence, a straight wound was made in the cell layer by a 50 µl sterile plastic pipette, followed by PBS washing twice. During the wound healing assay, cells were cultured in RPMI-1640 supplemented with 1% FBS (Costar; Corning, Inc.) at 37˚C for 24 h. Images were taken under an light microscope (Olympus Corporation) at x100 magnification and the cell migration distance was measured by Image-Pro Plus 4.1 (Media Cybernetics, Inc.). Six to eight horizontal lines were drawn randomly and mean width of the wound was calculated.

Cell invasion was assessed by performing Transwell invasion assays (Costar; Corning, Inc.) (37). Transwell apical chambers (cat. no. 353097; 24-well format; 8-µm pores; Corning Inc.) coated with 1% Matrigel (98 µl/chamber; cat. no. 356234, BD Biosciences) at 37˚C for 30 min were used for cell invasion. Transfected cells (1x105 cells) were maintained in 100 µl serum-free DMEM and seeded into the upper chamber at 37˚C with 5% CO2 for 16 h. By contrast, 0.2% BSA containing fibronectin (10 µg/ml; Bio-Techne) dissolved in 600 µl serum-free DMEM medium was added as a chemoattractant to the lower chamber. The non-invading cells on the apical chambers were removed with a cotton swab, and invading cells were fixed with 70% ethanol for 15 min at 25˚C followed by staining with Hemacolor® Rapid staining solution (Merck KGaA) at 25˚C for 30 min. The invasive cells were quantified by a light microscope (Olympus Corporation) at x200 magnification.

Statistical analysis

Statistical analyses were performed using SPSS software (version 15.0; SPSS, Inc.). The experiments were repeated for at ≥ three times and data are presented as the mean ± SD. The distribution of the data was analyzed using the Kurtosis test. Parametric tests were used for normally distributed data, whereas non-parametric tests were used for non-normally distributed data. For normally distributed quantitative data, comparisons between two groups were analyzed using the paired or unpaired Student's t-test, and comparisons between multiple groups were analyzed using one-way or two-way ANOVA followed by Tukey's post hoc test. For non-normally distributed data, comparisons between two groups were analyzed using the Mann-Whitney U test. To analyze categorical data, the χ2 test or Fisher's exact test was used. P<0.05 was considered to indicate a statistically significant difference.

Results

Expression levels of hsa_circRNA_102049 and miR-455-3p in CRC tissues and cells

The RT-qPCR results demonstrated that the relative hsa_circRNA_102049 expression levels were significantly upregulated (P<0.001; Fig. 1A), whereas miR-455-3p expression levels were significantly downregulated (P<0.001; Fig. 1B) in human CRC tumor tissues compared with adjacent non-tumor tissues.

Figure 1

hsa_circRNA_102049 is upregulated and miR-455-3p is downregulated in primary human CRC tumor tissues and CRC cell lines. Relative expression of (A) hsa_circRNA_102049 and (B) miR-455-3p in 68 paired CRC tumor tissues and adjacent non-tumor tissues. Relative expression of (C) hsa_circRNA_102049 and (D) miR-455-3p in CRC cell lines (SW480, LOVO, HT29, DLD-1, SW620 and HCT116) and a normal human intestinal epithelium cell line (NCM460) in six samples. *P<0.05, **P<0.01 and ***P<0.001 vs. NCM460 or non-tumor tissues. circRNA, circular RNA; miR, microRNA; CRC, colorectal cancer.

The associations between hsa_circRNA_102049 or miR-455-3p expression levels and the clinicopathological characteristics of patients with CRC are presented in Table I. Patients with higher grade malignancy are generally associated with increased hsa_circRNA_102049 expression and decreased miR-455-3p expression. hsa_circRNA_102049 and miR-455-3p expression levels in CRC were significantly associated with differentiation (both P<0.05), TNM stage (both P<0.05), lymphatic migration and invasion (P<0.05, respectively) and histological classification (both P<0.05), but not sex (P=0.60 and P=0.49, respectively) and age (P=0.50 and P=0.23, respectively).

The RT-qPCR results demonstrated a significant increase in hsa_circRNA_102049 expression levels in SW480 (P<0.05), LOVO (P<0.05), HT29 (P<0.01), DLD-1 (P<0.001), SW620 (P<0.001) and HCT116 (P<0.001) cells compared with that in the NCM460 cells (Fig. 1C). Furthermore, miR-455-3p expression levels were significantly lower in LOVO (P<0.05), HT29 (P<0.001), DLD-1 (P<0.001), SW620 (P<0.001) and HCT116 (P<0.001) cells compared with those in NCM460 cells (Fig. 1D). Among these cell lines, hsa_circRNA_102049 was highly expressed in HCT116 cells, but the least expressed in SW480 cells. However, miR-455-3p was highly expressed in SW480 cells, but the least expressed in HCT116 cells. Therefore, SW480 and HCT116 cells were selected for hsa_circRNA_102049 overexpression and knockdown, respectively.

miR-455-3p is a target gene of hsa_circRNA_102049

To determine whether miR-455-3p was a target gene of hsa_circRNA_102049, the effect of hsa_circRNA_102049 mRNA on miR-455-3p expression was assessed. The interaction between miR-455-3p and hsa_circRNA_102049 was verified by performing the dual-luciferase reporter assay. The RT-qPCR results demonstrated that hsa_circRNA_102049 overexpression significantly decreased miR-455-3p expression levels in SW480 cells compared with the control group (P<0.001; Fig. 2A). By contrast, hsa_circRNA_102049 knockdown significantly increased miR-455-3p expression levels in HCT116 cells compared with the NC siRNA group (P<0.001; Fig. 2B).

Figure 2

miR-455-3p is a direct target of hsa_circRNA_102049. Relative expression of hsa_circRNA_102049 and miR-455-3p in (A) hsa_circRNA_102049-overexpression SW480. Control is the pcDNA3.1 (+) circRNA mini vector. ***P<0.001 vs. Control. (B) Relative expression of hsa_circRNA_102049 and miR-455-3p in hsa_circRNA_102049-knockdown HCT116 cells. *P<0.05 and ***P<0.001 vs. NC siRNA. (C) Predicted 2D structure of the binding site between hsa_circRNA_102049 and miR-455-3p. (D) Luciferase activities of hsa_circRNA_102049-WT and hsa_circRNA_102049-Mut in NC mimic- or miR-455-3p mimic-transfected 293T cells. ***P<0.001 vs. NC mimic. miR/miRNA, microRNA; circRNA, circular RNA; WT, wild-type; Mut, mutant; UTR, untranslated region; NC, negative control; siRNA, small interfering RNA.

A binding site between miR-455-3p and hsa_circRNA_102049 was predicted using the microRNA database and the TargetScan online tool (Fig. 2C). The dual-luciferase reporter assay results demonstrated that, compared with the NC mimic group, miR-455-3p mimic significantly suppressed the luciferase activity of hsa_circRNA_102049-WT (P<0.001), but did not significantly affect the luciferase activity of hsa_circRNA_102049-Mut (Fig. 2D).

hsa_circRNA_102049 stimulates CRC cell proliferation by targeting miR-455-3p

miR-455-3p expression was overexpressed and knocked down by transfection with miR-455-3p mimic and inhibitor, respectively (Fig. 3A and B). The RT-qPCR results demonstrated that, compared with the control group, hsa_circRNA_102049 overexpression significantly decreased miR-455-3p expression levels (P<0.001), which were significantly increased by co-transfection with miR-455-3p mimic compared with NC mimic (P<0.001) in SW480 cells (Fig. 3C). By contrast, compared with the NC siRNA group, hsa_circRNA_102049 knockdown significantly increased miR-455-3p expression levels (P<0.001), which were significantly downregulated by co-transfection with miR-455-3p inhibitor compared with NC inhibitor (P<0.001) in HCT116 cells (Fig. 3D). Additionally, compared with the control group, hsa_circRNA_102049 overexpression significantly increased cell proliferation (P<0.05) potentially via downregulating miR-455-3p (Fig. 4A). Co-transfection with miR-455-3p mimic significantly reduced hsa_circRNA_102049-overexpression SW480 cell proliferation compared with the NC mimic group (P<0.05). HCT116 cell proliferation was significantly inhibited by hsa_circRNA_102049 knockdown compared with the NC siRNA group (P<0.05; Fig. 4B). miR-455-3p inhibitor significantly enhanced hsa_circRNA_102049-knockdown HCT116 cell proliferation compared with the NC inhibitor group (P<0.05).

Figure 3

Effect of hsa_circRNA_102049 and miR-455-3p expression on miR-455-3p expression levels. (A) Transfection efficiency of miR-455-3p mimic in SW480 cells. ***P<0.001 vs. miR-455-3p mimic. (B) Transfection efficiency of miR-455-3p inhibitor in HCT116 cells ***P<0.001 vs. miR-455-3p inhibitor. (C) Effect of hsa_circRNA_102049 and miR-455-3p overexpression on miR-455-3p expression. Control is the pcDNA3.1 (+) circRNA mini vector. ***P<0.001 vs. Control. ###P<0.001 vs. hsa_circRNA_102049 + NC mimic. (D) Effect of hsa_circRNA_102049 knockdown and miR-455-3p inhibitor transfection on miR-455-3p expression. ***P<0.001 vs. NC siRNA, ###P<0.001 vs. hsa_circRNA_102049 siRNA + NC inhibitor. MiR, microRNA; circRNA, circular RNA; NC, negative control; siRNA, small interfering RNA.

Figure 4

hsa_circRNA_102049 promotes CRC cell proliferation via regulating miR-455-3p in vitro. (A) Effects of hsa_circRNA_102049 and miR-455-3p on cell viability of SW480. Control is the pcDNA3.1 (+) circRNA mini vector. *P<0.05 vs. Control. #P<0.05 vs. hsa_circRNA_102049 + NC mimic. (B) Effects of hsa_circRNA_102049 knockdown and miR-455-3p inhibitor on HCT116 cell viability. *P<0.05 vs. NC siRNA. #P<0.05 vs. hsa_circRNA_102049 siRNA + NC inhibitor. (C) Colony formation was assessed. Control is the pcDNA3.1 (+) circRNA mini vector. (D) Colony formation was quantified for SW480. Control is the pcDNA3.1 (+) circRNA mini vector. ***P<0.001 vs. Control. ###P<0.001 vs. hsa_circRNA_102049 + NC mimic. (E) Colony formation was quantified for HCT116 cells. ***P<0.001 vs. NC siRNA. ###P<0.001 vs. hsa_circRNA_102049 siRNA + NC inhibitor. Scale bar, 100 µm. circRNA, circular RNA; CRC, colorectal cancer; miR, microRNA; NC, negative control; siRNA, small interfering RNA.

The colony formation assay results were similar to the CCK-8 assay results (Fig. 4C-E). In SW480 cells, hsa_circRNA_102049 overexpression significantly increased clone formation compared with the control group (P<0.001; Fig. 4C and D). Consistently, miR-455-3p mimic significantly decreased hsa_circRNA_102049 overexpression-induced clone formation compared with the NC mimic group (P<0.001). In HCT116 cells, hsa_circRNA_102049 knockdown significantly reduced clone formation compared with the NC siRNA group (P<0.001; Fig. 4C and E), whereas co-transfection with miR-455-3p inhibitor significantly increased clone formation compared with the NC inhibitor group (P<0.001).

hsa_circRNA_102049 promotes cell cycle entry and inhibits apoptosis by targeting miR-455-3p

In particular, hsa_circRNA_102049 overexpression resulted in a significantly decreased number of SW480 cells arrested at the G0/G1 phase and a significantly increased number of cells arrested at the S phase compared with the control group (both P<0.001; Fig. 5A and B). However, co-transfection with miR-455-3p mimic significantly increased the number of cells arrested at the G0/G1 phase and significantly decreased the number of cells arrested at the S phase in hsa_circRNA_102049-overexpression SW480 cells compared with the NC mimic group (both P<0.001). Moreover, a significantly increased number of HCT116 cells were arrested at the G0/G1 phase and a significantly decreased number of HCT116 cells were arrested at the S phase in the hsa_circRNA_102049 siRNA group compared with the NC siRNA group (both P<0.001; Fig. 5A and C). Co-transfection with miR-455-3p inhibitor significantly reduced the number of arrested cells at the G0/G1 phase and significantly increased cell cycle arrest at the S phase in hsa_circRNA_102049-knockdown HCT116 cells compared with the hsa_circRNA_102049 siRNAs + NC inhibitor group (both P<0.001).

Figure 5

hsa_circRNA_102049 overexpression promotes cell cycle entry in SW480 cells, whereas hsa_circRNA_102049 knockdown suppresses cell cycle entry in HCT116 cells via regulation of miR-455-3p in vitro. (A) Flow cytometry was performed to determine the DNA content of PI-stained cells at the G0/G1, S and G2 phases of the cell cycle. Control is the pcDNA3.1 (+) circRNA mini vector. (B) Effects of hsa_circRNA_102049 and miR-455-3p on the percentage of PI-stained cells at the G0/G1, S, and G2 phases of the cell cycle in SW480. Control, pcDNA3.1 (+) circRNA mini vector. ***P<0.001 vs. Control. ###P<0.001 vs. hsa_circRNA_102049 + NC mimic. (C) Effects of hsa_circRNA_102049 knockdown and miR-455-3p inhibitor on cell cycle distribution in HCT116 cells. ***P<0.001 vs. NC siRNA, ###P<0.001 vs. hsa_circRNA_102049 siRNA + NC inhibitor. circRNA, circular RNA; miR, microRNA; NC, negative control; siRNA, small interfering RNA.

The Annexin V-FITC/PI staining results suggested that hsa_circRNA_102049 overexpression significantly inhibited CRC cell apoptosis compared with that in cells transfected with the empty pcDNA vector (Fig. 6). Compared with the control group, the apoptotic rate of hsa_circRNA_102049-overexpression SW80 cells was significantly reduced (P<0.001; Fig. 6A and B). Conversely, co-transfection with miR-455-3p mimic significantly enhanced the apoptotic rate of hsa_circRNA_102049-overexpression SW80 cells compared with the NC mimic group (P<0.001; Fig. 6A and B). Moreover, compared with the NC siRNA group, the apoptotic rate of HCT116 cells was significantly increased in the hsa_circRNA_102049 siRNA group (P<0.001; Fig. 6A and C), whereas co-transfection with miR-455-3p inhibitor significantly decreased the apoptotic rate of hsa_circRNA_102049-knockdown HCT116 cells compared with the hsa_circRNA_102049 siRNAs + NC inhibitor group (P<0.001; Fig. 6A and C).

Figure 6

hsa_circRNA_102049 overexpression inhibits SW480 cell apoptosis, whereas hsa_circRNA_102049 knockdown promotes HCT116 cell apoptosis via regulation of miR-455-3p in vitro. (A) Cell apoptosis was assessed by performing flow cytometry. Control is the pcDNA3.1 (+) circRNA mini vector. The upper left quadrant represents necrotic cells (annexin V-/PI+), the upper right quadrant represents late apoptotic cells (annexin V+/PI+), the lower left quadrant represents live cells (annexin V-/PI-), and the lower right quadrant represents early apoptotic cells (annexin V+/PI-), respectively. (B) Effects of hsa_circRNA_102049 and miR-455-3p on early and late apoptosis in SW480. Control is the pcDNA3.1 (+) circRNA mini vector. ***P<0.001 vs. Control. ###P<0.001 vs. hsa_circRNA_102049 + NC mimic (C) Effects of hsa_circRNA_102049 knockdown and miR-455-3p inhibitor on early and late apoptosis in HCT116 cells. ***P<0.001 vs. NC siRNA, ###P<0.001 vs. hsa_circRNA_102049 siRNA+NC inhibitor. circRNA, circular RNA; miR, microRNA; NC, negative control; siRNA, small interfering RNA.

hsa_circRNA_102049 enhances CRC cell migration and invasion by targeting miR-455-3p

The wound healing assay results demonstrated that hsa_circRNA_102049 overexpression significantly increased SW480 cell migration compared with the control group (P<0.001; Fig. 7A and C). Co-transfection with miR-455-3p mimic significantly decreased the rate of wound healing in hsa_circRNA_102049-overexpression SW480 cells compared with the NC mimic group (P<0.001). Additionally, hsa_circRNA_102049 knockdown significantly inhibited HCT116 cell migration compared with the NC siRNA group (P<0.001), which was significantly reversed by co-transfection with miR-455-3p inhibitor compared with NC inhibitor (P<0.001) (Fig. 7B and D).

Figure 7

hsa_circRNA_102049 facilitates cell migration via regulation of miR-455-3p in vitro. (A) Representative images of the wound healing assay in SW480 cells. Control is the pcDNA3.1 (+) circRNA mini vector. (B) Representative images of the wound healing assay in HCT116 cells (scale bar, 50 µm). (C) Quantification of the wound healing rate in SW480 cells. Control is the pcDNA3.1 (+) circRNA mini vector. ***P<0.001 vs. Control. ###P<0.001 vs. hsa_circRNA_102049 + NC mimic. (D) Quantification of the wound healing rate in HCT116 cells. ***P<0.001 vs. NC siRNA, ###P<0.001 vs. hsa_circRNA_102049 siRNA + NC inhibitor. circRNA, circular RNA; miR, microRNA; NC, negative control; siRNA, small interfering RNA.

The Transwell invasion assay results indicated that hsa_circRNA_102049 overexpression significantly increased the number of invading SW480 cells compared with the control group (P<0.001; Fig. 8A and B). Co-transfection with miR-455-3p mimic significantly decreased cell invasion in hsa_circRNA_102049-overexpression SW480 cells compared with the NC mimic group (P<0.001). A significantly reduced number of invading HCT116 cells was observed in the hsa_circRNA_102049 siRNA group compared with the NC siRNA group (P<0.001), which was significantly increased by co-transfection with miR-455-3p inhibitor compared with NC inhibitor (P<0.001) (Fig. 8A and C).

Figure 8

hsa_circRNA_102049 promotes cell invasion via regulation of miR-455-3p in vitro. (A) Representative images of the Transwell invasion assay in SW480 and HCT116 cells (scale bar, 50 µm). Control is the pcDNA3.1 (+) circRNA mini vector. (B) Quantification of invading cells in four fields of view for SW480 cells. Control is the pcDNA3.1 (+) circRNA mini vector. ***P<0.001 vs. Control group. ###P<0.001 vs. hsa_circRNA_102049 + NC mimic. (C) Quantification of invading cells in four fields of view for HCT116 cells. ***P<0.001 vs. NC siRNA; ###P<0.001 vs. hsa_circRNA_102049 siRNA+NC inhibitor. circRNA, circular RNA; miR, microRNA; NC, negative control; siRNA, small interfering RNA.

Discussion

circRNAs, which are characterized as ceRNAs, have become a topic of interest in ceRNA research due to their ability to bind miRNA targets and block their inhibitory actions on target gene expression transcriptionally and post-transcriptionally (31,32). Several circRNAs functioning as ceRNAs are aberrantly expressed and serve crucial roles in important biological processes, including cell proliferation, apoptosis, migration and invasion, in diverse tumor tissues (32,39,40), including CRC (41-43). To the best of our knowledge, the present study demonstrated for the first time that hsa_circRNA_102049 derived from the TADA2A gene (alias of hsa_circ_0043278) was significantly increased in multiple CRC cell lines, including SW480, LOVO, HT29, DLD-1, SW620 and HCT116, and CRC tissues compared with NCM460 cells and non-tumor tissues. Compared with the NCM460 cell line and non-tumor tissues, miR-455-3p expression was significantly decreased in CRC cell lines and tissues, displaying the lowest and highest expression levels in HCT116 cells, and SW480 and LOVO cells, respectively. hsa_circRNA_102049 and miR-455-3p expression levels were associated with clinicopathological features of patients with CRC, including differentiation, TNM staging, lymphatic migration and invasion and histological classification, but not sex and age. The RT-qPCR and dual-luciferase reporter assay results demonstrated that miR-455-3p was a target gene of hsa_circRNA_102049. Therefore, to the best of our knowledge, the present study demonstrated for the first time that hsa_circRNA_102049 may serve an oncogenic role in CRC cells, whereas its target miR-455-3p may exert an inhibitory effect on CRC development and progression.

TADA2A (the gene symbol of hsa_circRNA_102049) was identified in the transcriptome of the HCT116 colon cancer cell line using RNA-sequencing (34), but its expression levels and functions in human CRC cells and tissues are not completely understood. In a recent study conducted by Wu et al (44), circTADA2A was highly expressed in the cytoplasm of both osteosarcoma tissues and cell lines, where it was characterized as a promoter of malignant tumor behavior, including migration, invasion and proliferation, as well as tumorigenesis and migration and invasion in vivo via the miR-203a-3p/cAMP responsive element binding protein 3 axis. Moreover, hsa_circ_0043278 (alias of hsa_circRNA_102049) was reported to be significantly upregulated in non-small cell lung cancer (NSCLC), where it promoted NSCLC progression (proliferation, invasion, and migration) in vitro and in vivo by increasing Rho associated coiled-coil containing protein kinase 1, cyclin dependent kinase inhibitor 1B and AKT3 expression levels via direct downregulation of miR-520f (33). To the best of our knowledge, the present study demonstrated for the first time that hsa_circRNA_102049 was significantly upregulated in CRC tissues and SW480, LOVO, HT29, DLD-1, SW620 and HCT116 CRC cell lines compared with non-tumor tissues and the NCM460 cell line, respectively. Moreover, aberrant hsa_circRNA_102049 expression was highly related to the clinicopathological features of patients with CRC. The results also demonstrated that hsa_circRNA_102049 functioned as a tumor promoter by enhancing CRC cell proliferation, migration and invasion, as well as promoting cell cycle entry by inducing CRC cell arrest at the S and G2/M phases and inhibiting apoptosis in vitro.

To further explore the mechanisms underlying the tumor-promoting activities of hsa_circRNA_102049 in CRC, miRDB bioinformatics software was used to predict the target genes. miR-455-3p was identified as a target gene of hsa_circRNA_102049, which was verified by assessing the effects of hsa_circRNA_102049 on miR-455-3p expression levels and performing dual-luciferase reporter assays. Compared with the corresponding control groups, hsa_circRNA_102049 overexpression significantly downregulated miR-455-3p expression, whereas hsa_circRNA_102049 knockdown significantly increased miR-455-3p expression levels. The results also indicated a tumor inhibitor role of miR-455-3p CRC cells in vitro. These results suggested that hsa_circRNA_102049 promoted CRC progression and malignant behavior via sponging miR-455-3p, highlighting the potential use of hsa_circRNA_102049 and miR-455-3p as novel therapeutic targets for CRC. In a study conducted by Zheng et al (23), miR-455-3p overexpression significantly suppressed tumor cell proliferation and increased cell apoptosis, which was reversed by miR-455-3p knockdown in human HCT116 colon cancer cells. However, miR-455-3p expression levels in HCT116 colon cancer cells had not been previously reported. The results of the present study filled the gap in knowledge regarding the expression patterns and functions of miR-455-3p in human CRC. The present study was consistent with the finding that miR-455-3p is expressed at low levels in diverse types of cancer and functions as a tumor suppressor (19-22).

The results of the present study suggested that hsa_circRNA_102049 may promote malignant behavior in CRC by sponging miR-455-3p. To the best of our knowledge, the present study was the first to characterize the expression patterns and roles of hsa_circRNA_102049 in CRC, supporting its potential use as a novel biomarker for the diagnosis and treatment of CRC.