Abstract
Background: Circular RNAs (circRNAs) are a class of non-coding RNAs broadly expressed in cells of various species. However, the molecular mechanisms that link circRNAs with colorectal cancer (CRC) are not well understood. In the present study, we attempted to provide novel basis for targeted therapy for CRC from the aspect of circRNA–microRNA (miRNA)–mRNA interaction.
Methods: We investigated the expression of circRNAs in five paired CRC tissues and adjacent non-tumor tissues by microarray analysis. Differentially expressed circRNAs were identified between CRC tissues and non-cancerous matched tissues. We focused on hsa_circ_0005100, which is located on chromosome 1 and derived from FMN2, and thus we named it as circFMN2. The expression of circFMN2 was detected in 88 CRC tissues and cell lines by quantitative real-time PCR. Functional assays were performed to evaluate the effects of circFMN2 on proliferation in vitro, and on tumorigenesis in vivo. The relationship between circFMN2 and miR-1182 was confirmed by luciferase reporter assay.
Results: circFMN2 was found to be significantly up-regulated in CRC tissues and cell lines. Moreover, knockdown of circFMN2 significantly inhibited cell proliferation and migration in vitro. Bioinformatics analysis predicted that there is a circFMN2/miR-1182/hTERT axis in CRC progression. Dual-luciferase reporter system validated the direct interaction of circFMN2, miR-1182, hTERT. Western blot verified that inhibition of circFMN2 decreased hTERT expression. Importantly, we demonstrated that circFMN2 was up-regulated in serum exosomes from CRC patients.
Conclusion: In conclusion, circFMN2 is a central component linking circRNAs to progression of CRC via an miR-1182/hTERT axis.
Introduction
Colorectal cancer (CRC) remains one of the most common and deadly cancers [1]. The incidence of CRC has increased in the last 20 years [2]. Despite considerable progress made in the diagnosis and therapy of this disease, the mortality of CRC is still high [3,4]. Therefore, the underlying mechanisms of CRC progression need to be further elucidated.
Circular RNAs (circRNAs) are a novel class of endogenous RNAs that have a covalent closed loop structure [5]. It is highly evolutionarily conserved and stable and particularly resistant to RNases activity [6]. Growing evidence has indicated that circRNAs are widely involved in diverse physiological and pathological processes, especially in the generation and development of tumors [7–10]. It is reported that circRNAs act as a microRNA (miRNA) sponge which absorbs miRNAs and then regulates the expression of miRNA targeted genes [11]. In addition, the inherent stability of circRNAs conferred by the circular structure allows them to be enriched in the exosomes and stably present in plasma, saliva and other peripheral tissues, which renders them potential diagnostic molecular markers for various diseases [12–15].
Exosomes are small membrane-derived vesicles with a diameter of approximately 30–150 nm [16]. It has been well established that exosomes originate from multivesicular endosomes (multivesicular bodies) by inverse budding to form multivesicular bodies, and are released into the extracellular space when the multivesicular body fuses with cell membrane exosomes [17]. Exosomes can regulate the physiological and pathological functions of various cancers as message transmitters in intercellular communication because they contain a variety of proteins, lipids, circRNAs, and microRNAs [18]. Recent studies have demonstrated that circRNAs are abundant and stable in exosomes, however, there is little evidence for the expression profile and potential function of secreted circRNAs in CRC.
In the present study, we compared the expression patterns of circRNAs between CRC patients and controls. We identified 448 significantly dysregulated circRNAs in CRC tissues. We focused on hsa_circ_0005100, which is located on chromosome 1 and derived from FMN2, and thus we named it as circFMN2. We further tested circFMN2 in 58 pairs of CRC samples by qRT-PCR and the results showed that the expression of circFMN2 was markedly elevated both in CRC tissues and exosomes from CRC plasma. Silencing circFMN2 suppressed CRC cell growth both in vitro and in vivo. In mechanism, circFMN2 served as a sponge for miR-1182 promoting CRC progression. Thus, we suggested that circFMN2 might act as an effective therapeutic target for the treatment of CRC. Our findings will provide new insights into the regulatory mechanisms of circFMN2 in CRC progression.
Materials and methods
Tissue specimens
The study was approved by the Clinical Research Ethics Committee of China-Japan Union Hospital of Jilin University and written informed consent was obtained from all the patients. The paired tissue samples were obtained from surgical specimens at the Department of Endoscopy Center, China-Japan Union Hospital of Jilin University. The samples were freshly frozen in liquid nitrogen and stored at −80°C. A total of 62 CRC tissues and adjacent normal colorectal tissues were included in the study. For exosome purification, whole blood samples were collected from these 62 CRC patients and healthy controls. Fresh plasma samples (3 ml) were collected in ethylenediamine tetraacetic acid tubes from each of the subjects. These samples were centrifuged at 3000×g for 10 min at 4°C and then stored at −80°C.
Agilent methods for ceRNA array
Total RNA of exosomes was extracted from plasma using TRIzol LS reagent (Invitrogen, Carlsbad, CA) according to the manufacturer’s instructions and quantified using a NanoDrop ND-2000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE). Total RNA from exosomes were isolated using TRIzol reagent (Invitrogen, Carlsbad, CA) according to the manufacturer’s instructions and purified using an RNeasy Mini Kit (Qiagen GmBH, Hilden, Germany). RNA samples of each group were then used to generate fluorescently labeled cRNA targets for the human ceRNA array V1.0 (4× 180 K; Shanghai Biotechnology Corporation, Shanghai, China). The labeled cRNA targets were then hybridized with the slides. After hybridization, slides were scanned on the Agilent Microarray Scanner (Agilent Technologies, Inc., Santa Clara, CA). Data was extracted using Feature Extraction software 10.7 (Agilent Technologies, Inc.). Raw data was normalized by the quantile algorithm of the Limma package in the R program. Microarray experiments were performed at Shanghai Biotechnology Corporation following the protocol of Agilent Technologies, Inc. Ratios between CRC and normal subjects (NC) were calculated. Genes with a fold change (FC) of at least 2 were selected for further analysis. The selected parent genes of circRNA were grouped into functional categories based on the Gene Ontology database (http://www.geneontology.org/) and functional pathways (Kyoto Encyclopedia of Genes and Genomes) were also analyzed using the online enrichment analysis tool of Shanghai Biotechnology Corporation.
The Cancer Genome Atlas dataset analysis
The data and the corresponding clinical information of patients were collected from the Cancer Genome Atlas (TCGA) database (http://cancergenome.nih.gov/). We used the edgeR package of R packages to perform the difference analysis (http://www.bioconductor.org/packages/release/bioc/html/edgeR.html) and used the pheatmap package of R packages to perform the cluster analysis (https://cran.r-project.org/web/packages/pheatmap/index.html). Sva R package was used to remove the batch effect. Genes with adjusted P-values <0.05 and absolute FCs > 1.5 were considered as differentially expressed genes. Kaplan–Meier survival curves were drawn to analyze the relationships between genes and overall survival in the survival package. The corresponding statistical analysis and graphics were performed in R software (R version 3.3.2).
Cell lines and culture conditions
The human CRC cell lines (HCT116, DLD1, SW480, RKO, HT-29) were purchased from the Institute of Biochemistry and Cell Biology of the Chinese Academy of Sciences (Shanghai, China). The human colonic epithelial cell line HCoEpiC and the 293T cell line were obtained from American Type Culture Collection (Manassas, VA). Cells were cultured in RPMI 1640 or DMEM (Gibco, Grand Island, NY, U.S.A.) supplemented with 10% fetal bovine serum (10% FBS), 100 U/ml penicillin, and 100 mg/ml streptomycin (Gibco) in humidified air at 37°C with 5% CO2.
RNA extraction and qRT-PCR
TRIzol Reagent (Invitrogen) was used to extract total RNA from tissues. For circRNAs, RNase R was used to degrade linear RNA, which have poly (A), and amplified by divergent primer. qRT-PCR analysis on circRNA and mRNA was performed using Prime Script RT reagent Kit (TaKaRa) and SYBR Premix Ex Taq II (TaKaRa). The b-actin was used as an endogenous control. For miR-1182 analysis, miRNA was treated with DNase I to eliminate genomic DNA and cDNA was synthesized by Mir-X miR First-Strand Synthesis Kit (TaKaRa). SYBR Premix Ex Taq II (TaKaRa) was used for qRT-PCR. The expression was normalized to RNU6-2. The 2−ΔΔCt method was adopted to calculate relative expression of circFMN2 and hsa-miR-1182.
Cell transfection
Small interfering RNAs (siRNAs) targeting the back-splice junction sequences of circFMN2 and the respective negative control oligonucleotides were synthesized by GenePharma (Shanghai, China). The siRNA sequences: si-circFMN2#1: AGACTTGAAAGCTGTTGTGAA, si-circFMN2#2: AAGACTTGAAAGCTGTTGTGA, and si-circFMN2#3: AAGAAAGACTTGAAAGCTGTT. A pcDNA 3.1 circRNA mini vector was used to ectopically express circFMN2 level in CRC cells. The miR-1182 mimics, inhibitor and NC were synthesized from GenePharma (Shanghai, China). We used Lipofectamine 3000 for transfection following the protocols of the manufacturer.
Cell counting kit-8 and colony formation assay
Cell proliferation rates were detected with Cell Counting Kit-8 (CCK-8, Dojindo, Japan). A total of 2 × 103 cells were cultured in each 96-well plate and the absorbance at 450 nm was measured after incubation with 10 μl CCK-8 solution at 37°C for 1 h. Cells (500 per well) were seeded in six-well plates for 1–2 weeks. The colonies were stained with 0.1% Crystal Violet solution diluted by anhydrous methanol for 15 min. Colonies were microscopically examined and counted. All detections were carried out in triplicate.
Cell cycle analysis
Cells for cell cycle analysis were stained with Propidium Oxide by the CycleTEST PLUS DNA Reagent Kit (BD Biosciences) following the protocol and were analyzed by FACScan. The percentage of the cells in G0–G1, S, and G2–M phase were counted and compared.
Apoptosis analysis
Cell apoptosis was analyzed using the Annexin V-FITC Apoptosis Detection Kit (KeyGen Biotech, Nanjing, China) according to the manufacturer’s instructions. Forty-eight hours after transfection, cells were stained with FITC and Propidium Iodide. Then flow cytometry was performed using a FACS Canto II (BD Biosciences) and the data were analyzed using BD FCSDiva Software and FCS Express 5 software (De Novo Software, Los Angeles, CA).
Immunohistochemistry
For each patient sample, three paraffin sections of 5 μm were prepared, one for Hematoxylin and Eosin (HE) staining and the other two for immunohistochemical staining. Phosphate-buffered saline (PBS) instead of primary antibodies was used for negative control, and the breast cancer tissue was used for positive control. Sections were dewaxed using xylene, followed by hydration with ethanol solutions, and addition of EDTA for antigen retrieval. Later, sections were blocked with normal goat serum for 30 min to eliminate non-specific binding. Sections were incubated with anti-human FMN2 polyclonal antibody (1:300; ab72052, Abcam Inc., Cambridge, MA, U.S.A.) and and anti-human hTERT polyclonal antibody (1:100; ab183105, Abcam Inc., Cambridge, MA, U.S.A.). Sections were then incubated with biotin-labeled secondary antibodies for 30 min at room temperature, followed by staining with diaminobenzidine (DAB). Finally, the sections were counterstained with Hematoxylin. The result of staining was determined by two doctors who did not know the clinical condition of patients. The proportions of positive cells of 0, 1–5, 6–25, 26–75, and 76–100% were assigned the scores of 0, 1, 2, 3, and 4, respectively. Scores of 0–2 were considered as negative expression, and scores of 3–4 were considered as positive expression.
Xenografts in mice
Approximately 1 × 106 cells were injected subcutaneously into the right neck of male BALB/C nude mice (age, 4–6 weeks; weight, 18–22 g, five mice per group) purchased from Slaccas (Slaccas Laboratory Animal, Shanghai, China). The length and width of tumor xenografts were measured weekly by Vernier calipers, and the tumor volume was calculated using the following formula: Volume (mm3) = 0.5 × width2 × length. Six weeks after injection, the mice were killed by cervical dislocation. The animal studies were conducted in accordance with the ethical guidelines for animal experiment and approved by the Animal Care and Use Committee of the China-Japan Union Hospital of Jilin University. The animal work took place in the animal central of China-Japan Union Hospital of Jilin University.
Luciferase reporter assay
The recombination luciferase plasmids containing circFMN2 full-length sequences with wild-type or mutant miR-1182 binding site were constructed from Genecopoeia company (Rockville, MD, U.S.A.). Subsequently, luciferase plasmids were co-transfected with control or miR-1182 mimics into cells by Lipofectamine 3000 (Invitrogen). After 48 h of treatment, the luciferase activity in each group was measured by using a dual luciferase reporter assay system (Promega, Madison, WI, U.S.A.).
RNA immunoprecipitation assay
RNA immunoprecipitation (RIP) assay was performed using an EZ-Magna RiP Kit (Millipore, Billerica, MA, U.S.A.) in accordance with the manufacturer’s instructions. Cells were lysed at 70–80% confluence in RIP lysis buffer, and then incubated with magnetic beads conjugated with human anti-Ago2 antibody (Millipore) and normal mouse IgG control (Millipore) in RIP buffer. The RNAs in the immunoprecipitates were isolated with TRIzol reagent and analyzed by qRT-PCR.
Plasma exosome isolation
First, the samples were centrifuged twice at 3000×g and 10000×g for 20 min at room temperature to remove cells and other debris in the plasma. The supernatants were then centrifuged at 100000×g for 30 min at 4°C to remove microvesicles that were larger than exosomes, harvested, and again centrifuged at 10000×g for 70 min at 4°C. Subsequently, the supernatants were gently decanted, and the exosome sediments were re-suspended in PBS. Concentration of exosomes was determined using BCA method as recommended by the manufacturer (Thermo Scientific, U.S.A.).
Transmission electron microscopy
The exosome suspension was diluted to 0.5 mg/ml with PBS, and then spotted on to a glow-discharged copper grid placed on a filter paper and dried for 10 min by exposure to infrared light. Next, the exosome samples were stained with one drop of phosphotungstic acid (1% aqueous solution) for 5 min and dried for 20 min by exposure to infrared light. Finally, the exosomes were visualized under a transmission electron microscope (HT7700, Hitachi, Tokyo, Japan) at 100 keV.
Nanoparticle tracking analysis
Briefly, the exosomes were resuspended in PBS and filtered with a syringe filter (Millipore). Then, the samples were diluted until individual nanoparticles could be tracked. The size distribution of the exosomes was evaluated using a NanoSight NS300 instrument (Malvern Instruments Ltd, Worcestershire, U.K.).
Western blot assay
The lysates from cells were collected by RIPA buffers (Beyotime Biotechnology, Shanghai, China) and then boiled for 5 min at 100°C. After that, the proteins were transferred to PVDF membrane and blocked by non-fat dried milk. The membrane was incubated with primary antibodies at 4°C overnight. The next day, the membrane was washed strictly and probed with HRP–conjugated secondary antibody, followed by visualization with ECL Plus chemiluminescence reagent (Beyotime Biotechnology).
Statistical analysis
All statistical analyses were performed by SPSS 19.0 software. The differences between two groups were calculated using Student’s t or Chi-square test. Survival data were evaluated using the Kaplan–Meier method and log-rank test. P<0.05 was considered statistically significant.
Results
Profiles of circRNAs in CRC
A total of 11245 circRNAs were detected in CRC and paired non-tumorous samples by the circRNA microarray analysis. Among them, 1075 circRNAs were significantly aberrantly expressed (P<0.05 and fold-change > 2.0) between CRC tissues and paired non-tumorous tissues. Of these circRNAs, 558 were significantly up-regulated and 517 were significantly down-regulated in CRC tissues compared with paired adjacent normal tissues. Hierarchical clustering was then performed to demonstrate the five most up-regulated circRNAs (hsa_circ_0005100, hsa_circ_0001821, hsa_circ_0084188, hsa_circ_0000372, and hsa_circ_0001178) and five most down-regulated circRNAs (hsa_circ_0000374, hsa_circ_0036592, hsa_circ_0006271, hsa_circ_0009155, and hsa_circ_0002211) expression patterns among the sets (Figure 1A).
Deregulated circRNAs in CRC tumor tissues
The five most up-regulated circRNAs were selected and validated by qRT-PCR using 88 CRC and paired non-tumorous tissue samples. As shown in Figure 1B–F, except for hsa_circ_0000372, the circRNAs displayed a consistent expression level between the microarray and qRT-PCR analyses. There was an increasing trend in hsa_circ_0005100 (chr1:240458121-240497529) levels from non-tumorous tissues to CRC tissues, with more than 10 FC from microarray analysis. By browsing the human reference genome (GRCh37/hg19), we identified that hsa_circ_0005100 is derived from FMN2, which is located on chromosome 1, with a spliced mature sequence length of 612 base pairs (bp), and thus we named it as circFMN2.
To further investigate the role of circFMN2 in CRC, the relationship between circFMN2 expression in CRC tissues and clinicopathological characteristics of CRC patients was analyzed. Using the median expression level of circFMN2 as cut-off value, patients who expressed circFMN2 equal to or greater than the average level were assigned to the ‘circFMN2 high’ group. As shown in Table 1, high expression of circFMN2 in CRC tissues was significantly correlated with larger tumor size (P=0.012), advanced tumor stage (P=0.011), distant metastasis (P=0.001), and TNM stage (P=0.03), but not with other characteristics of CRC that included age (P=0.521), gender (P=0.286), lymph node metastasis (P=0.829), or histological grade (P=0.134). The data provide evidence of the importance of circFMN2 in the growth and tumorigenesis of CRC.
Clinicopathological features . | Ovearall (n=88) . | circFMN2 . | P . | |
---|---|---|---|---|
. | . | High (n=44) . | Low (n=44) . | . |
Age | ||||
<60 | 48 | 26 | 22 | 0.521 |
≥60 | 40 | 18 | 22 | |
Gender | ||||
Male | 46 | 20 | 26 | 0.286 |
Female | 42 | 24 | 18 | |
Tumor size | ||||
>5 cm | 30 | 21 | 9 | 0.012 |
≤5 cm | 58 | 23 | 35 | |
Tumor stage | ||||
T1 and T2 | 60 | 24 | 36 | 0.011 |
T3 and T4 | 28 | 20 | 8 | |
Lymph node metastasis | ||||
N1–2 | 38 | 20 | 18 | 0.829 |
N0 | 50 | 24 | 26 | |
Distant metastasis | ||||
Absent | 66 | 26 | 40 | 0.001 |
Present | 22 | 18 | 4 | |
TNM stage | ||||
I/II | 62 | 26 | 36 | 0.03 |
III/IV | 26 | 18 | 8 | |
Histological grade | ||||
Well/moderately | 48 | 20 | 28 | 0.134 |
Poorly | 40 | 24 | 16 |
Clinicopathological features . | Ovearall (n=88) . | circFMN2 . | P . | |
---|---|---|---|---|
. | . | High (n=44) . | Low (n=44) . | . |
Age | ||||
<60 | 48 | 26 | 22 | 0.521 |
≥60 | 40 | 18 | 22 | |
Gender | ||||
Male | 46 | 20 | 26 | 0.286 |
Female | 42 | 24 | 18 | |
Tumor size | ||||
>5 cm | 30 | 21 | 9 | 0.012 |
≤5 cm | 58 | 23 | 35 | |
Tumor stage | ||||
T1 and T2 | 60 | 24 | 36 | 0.011 |
T3 and T4 | 28 | 20 | 8 | |
Lymph node metastasis | ||||
N1–2 | 38 | 20 | 18 | 0.829 |
N0 | 50 | 24 | 26 | |
Distant metastasis | ||||
Absent | 66 | 26 | 40 | 0.001 |
Present | 22 | 18 | 4 | |
TNM stage | ||||
I/II | 62 | 26 | 36 | 0.03 |
III/IV | 26 | 18 | 8 | |
Histological grade | ||||
Well/moderately | 48 | 20 | 28 | 0.134 |
Poorly | 40 | 24 | 16 |
Knockdown of circFMN2 inhibits CRC growth in vitro
The expression level of circFMN2 in multiple CRC and a normal colonic epithelial cell (HCoEpiC) were measured. The expression level of circFMN2 in CRC cell lines was generally higher than that in the HCoEpiC normal colonic cell line (Figure 2A; P<0.01). The expression of circFMN2 was highest in HCT116 cells and lowest in HT-29 cells. Using the backsplice junction-specific siRNA (respectively si-circFMN2#1, si-circFMN2#2, or si-circFMN2#3), we successfully knocked down circFMN2 expression in HCT116 cells (Figure 2B; P<0.01). The si-circFMN2#3 was used in the following experiments. Using circFMN2-overexpressing lentivirus, the circFMN2 expression in HT-29 cells was shown to be dramatically increased (Figure 2C).
Knockdown of circFMN2 inhibits CRC growth in vitro
Analyses using CCK-8 (Figure 2D,E; P<0.01) and colony formation (Figure 2F) revealed that the knockdown of circFMN2 significantly restrained the proliferation of HCT116 cells, while ectopic expression of circFMN2 significantly promoted the proliferation of HT-29 cells. However, CCK-8 assay revealed that the knockdown of circFMN2 did not inhibit the proliferation of HCoEpiC cells, and ectopic expression of circFMN2 did not promote the proliferation of HCoEpiC cells (Supplementary Figure S1A,B). We further investigated whether circFMN2 affected cell cycle progression or apoptosis of CRC cells using flow cytometry. The knockdown of circFMN2 significantly increased the percent of HCT116 cells in G0/G1 phase and decreased the percent of HCT116 cells in G2/M phase. The overexpression of circFMN2 significantly decreased the percent of HT-29 cells in G0/G1 phase and increased the percent of HT-29 cells in G2/M phase (Figure 3A,B; P<0.01). However, circFMN2 overexpression and knockdown did not affect apoptosis of CRC cells (Figure 3C,D). The collective findings suggested that knockdown of circFMN2 inhibits the proliferation of CRC cells by facilitating DNA synthesis.
Knockdown of circFMN2 inhibits the proliferation of CRC cells by inducing cell cycle arrest
CircFMN2 has no effect on its linear transcript
Some circRNAs regulate the expression and function of the corresponding linear transcripts. Therefore, the regulatory relationship between circFMN2 and its linear transcript (FMN2) was investigated in the present study. We first analyzed the colon adenocarcinoma (COAD) dataset from the TCGA database and found that the level of FMN2 was significantly lower in 275 COAD tissues than 349 normal tissues (P<0.05; Figure 4A). However, Kaplan–Meier survival analysis from TCGA COAD datasets suggested that FMN2 expression is not significantly associated with worse overall survival (OS) (log-rank test, P=0.69, Figure 4B) and disease-free survival (DFS) of COAD patients (log-rank test, P=0.26, Figure 4C). Moreover, we used IHC assay to analyze the expression of FMN2 in our cohort of 88 pairs of CRC tumor tissues. FMN2 was significantly down-regulated in CRC tissues compared with adjacent non-tumorous tissues (Figure 4D). FMN2 did not change in both mRNA and protein levels when the expression of circFMN2 was artificially changed in CRC cells (Figure 4E,F and Supplementary Figure S1C,D). These data indicated that FMN2 is not the target gene of circFMN2.
CircFMN2 has no effect on its linear transcript
CircFMN2 promoted CRC progression by sponging miR-1182
Increasing studies have verified that circRNAs can exert their functions by acting as ceRNAs to sponge miRNAs [16–18]. Considering the cytoplasmic localization of circFMN2 in CRC cells, we investigated whether circFMN2 exerted function in CRC in an miRNA-mRNA-dependent manner. Therefore, we used the StarBase v2.0 target prediction tool to find 28 potential miRNAs that could bind to circFMN2 (Supplementary Data S1). To figure out the miRNA implicated in CRC, all candidate miRNAs were subjected to microarray analysis. Among them, ten miRNAs were significantly down-regulated in response to circFMN2 overexpression (Figure 5A). For further selection, we chose top five down-regulated miRNAs to do circRIP assay. The results demonstrated the interaction of only miR-1182 with circFMN2 (Figure 5B; P<0.01). To validate the interaction between circFMN2 and miR-1182, we obtained the binding sequence between miR-1182 and circFMN2 (Figure 5C). Next, luciferase reporter assays were carried out to verify whether there was a direct interaction between circFMN2 and miR-1182. As shown in Figure 5D, miR-1182 mimics could significantly decrease the luciferase activity of the vector containing the wild-type miR-1182 binding site within circFMN2 in HCT116 cells, while the luciferase activity of the vector containing the mutant binding site was restored. In a further RIP experiment, circFMN2 and miR-1182 simultaneously existed in the production precipitated by anti-AGO2 (Figure 5E; P<0.01), suggesting that miR-1182 is circFMN2-targeting miRNA. These outcomes indicated that the interaction of circFMN2 and miR-1182 was realized by the putative binding site.
CircFMN2 promoted CRC progression by sponging miR-1182
The qRT-PCR analysis indicated that miR-1182 was expressed at low level in CRC tissues (Supplementary Figure S2A; P<0.01). The expression of miR-1182 was obviously decreased in HCT116 cells, indicating the opposite result to circFMN2 expression (Supplementary Figure S2B; P<0.01). The decreased expression level of miR-1182 in circFMN2-overexpressed CRC cells indicated the negative regulatory effect of circFMN2 on miR-1182 expression (Figure 4F and Supplementary Figure S2C; P<0.01). Furthermore, miR-1182 had the similar tumor-suppressing effects as circFMN2 silencing (Figure 5G and Supplementary Figure S2D; P<0.01). These data suggested that circFMN2 might exert its functions by sponging miR-1182 in CRC.
CircFMN2 up-regulates hTERT expression via sponging miR-1182
Previous studies have reported that several genes are directly targeted by miR-1182 in human cancers including hTERT, SKA1, and CDK14 [19–21]. However, we did not detect significant changes in the expression levels of these mRNAs in HCT116 cells transfected with an miR-1182 mimic or in HT-29 cells transfected with an miR-1182 inhibitor except for the expression of the hTERT (Figure 6A,B; P<0.01). TargetScan online software demonstrated that hTERT was a candidate target of miR-1182 (Figure 6C). To verify whether the 3′ UTR of hTERT mRNA was a target of miR-1182 in 293T cells, a luciferase reporter gene assay was used. The wild-type (wt) 3′ UTR sequence or mutant (mu) 3′ UTR sequence of hTERT was cloned into a luciferase reporter vector. The luciferase activity was significantly inhibited by the miR-1182 mimics in wt 3′ UTR sequence-transfected 293T cells. Conversely, the luciferase activity was not inhibited by the miR-1182 mimics in mu 3′ UTR sequence-transfected cells (Figure 6D; P<0.01). These findings indicate a direct interaction between miR-1182 and hTERT mRNA in CRC cells. We determined the expression of hTERT in 88 CRC patient tissues using qRT-PCR and IHC. We also found that hTERT expression was higher in CRC tissues compared with their normal counterparts (Figure 6E,F; P<0.01). Then, we found that levels of hTERT were negatively regulated by circFMN2 knockdown, and the effect of inhibition of circFMN2 was attenuated by miR-1182 inhibitor (Figure 6G; P<0.01). These data further demonstrated the regulatory network of circFMN2/miR-1182/hTERT. According to above data, we confirmed that circFMN2 can exert function in CRC by sponging miR-1182 to up-regulate hTERT expression.
CircFMN2 up-regulates hTERT expression via sponging miR-1182
Knockdown of circFMN2 inhibits CRC growth in vivo
HCT116 cells were used to establish a stable CRC cell line with low circFMN2 expression via lentiviral-circFMN2-RNAi. HCT116 cells with low or normal circFMN2 expression induced by transfection with lentiviral or control vectors were subcutaneously implanted into nude mice. Compared with the control group, smaller tumor size and lower tumor weight were observed in the low circFMN2 expression group group (Figure 7A–C, P<0.01). IHC staining results showed that knockdown circFMN2 significantly inhibited the expression levels of Ki-67 when compared with the control group (Figure 7D). These results demonstrated that circFMN2 can promote the growth of CRC in vivo.
circFMN2 promoted CRC cell growth in vivo
CircFMN2 is secreted by exosomes into serum of CRC patients
Finally, in our current study, we collected abundant serums from 35 CRC patients and 30 normal people. After isolation of serum exosomes by sequential centrifugation, transmission electron microscopy (TEM) analysis showed that isolated CRC-secreted exosomes had similar morphologies (30–150 nm in diameter) and exhibited a round-shaped appearance (Figure 8A). The nanoparticle tracking analysis (NTA) results demonstrated that isolated CRC-secreted exosomes showed a similar size distribution, and the peak size range was 80–135 nm. Moreover, the presence of the exosome markers CD63, TSG101, and HSP70 were confirmed by Western blot (Figure 8B). Our results showed that circFMN2 expression is detectable in extracted serum exosomes derived from CRC patients (Figure 8C; P<0.01). Furthermore, there was a significant inverse correlation between the expression levels of circFMN2 and miR-1182 in serum exosomes derived from CRC patients (Figure 8D, r = −0.821, P<0.0001).
circFMN2 is secreted by exosomes into serum of CRC patients
Discussion
circRNAs are a group of non-coding RNAs that in recent years have emerged as potential regulators of many cell processes [22,23]. Currently, with the development of deep RNA sequencing (RNA-seq) technologies and novel bioinformatics approaches, abundant and diverse circRNAs have been detected and identified, including few circRNA functions, such as regulating transcription in the nucleus, functioning as efficient microRNA sponges, competing with pre-mRNA splicing, and serving as circRNA–protein interactions [24,25]. However, details of the mechanisms of circRNAs in CRC remain ambiguous. Therefore, it is necessary to seek out new biomarkers and explore their functions. In this present study, we analyzed the expression profiles of circRNAs from CRC and matched non-tumor normal tissues by microarray, and focused on the role and underlying mechanism of the increased circFMN2 expression. Our findings report that circFMN2 expression is significantly elevated in CRC tissues and cell lines. Elevated expression of circFMN2 was positively correlated with correlated with larger tumor size, advanced tumor stage, distant metastasis, and TNM stage. We also discovered that circFMN2 exerted oncogenic function in the progression of CRC in vitro and in vivo. Functionally and mechanistically, circFMN2 promoted CRC progression by sponging miR-1182 and up-regulating hTERT expression. Finally, we showed that overexpressed circFMN2 was secreted by exosomes into the serum of CRC patients, suggesting that circFMN2 might be a novel clinical molecular marker for CRC patients. These results revealed that circFMN2 may act as a biomarker and therapeutic target for CRC.
Increasing evidence shows that circRNAs have a range of functions including regulation of RNA binding proteins (RBPs), promotion of mRNA transcription in cis and post-transcriptional regulation through sponging of miRNAs [26–28]. MiRNAs are a class of non-coding small RNAs of approximately 19–23 nucleotides in length, which are widely involved in various biological behaviors of cancer cells [29]. MiRNAs can post-transcriptionally reduce the levels of specific target protein coding gene expression by binding to the 3′UTR of target mRNAs and resulting in translation inhibition or mRNA degradation. Studies showed that miRNAs are closely related to cell biological processes of CRC. Recently, miR-1182 has been reported as a tumor suppressor that is significantly down-regulated in multiple cancer types. We found that miR-1182 was expressed at a significantly lower level in CRC tumor tissues than in paratumor tissues. In our study, hTERT was predicted as the candidate target gene of miR-1182 by miRDB and TargetScan, and was further testified by dual-luciferase reporter assay. Next, we verified that circFMN2 had an endogenous sponge-like effect on miR-1182 in CRC. Furthermore, bioinformatics prediction and a luciferase reporter assay showed that circFMN2 and the hTERT 3ʹUTR share identical miR-1182 response elements and might therefore bind competitively to miR-1182. Third, circFMN2 could bind directly to miR-1182 in an AGO2-dependent manner. Finally, circFMN2 could control the hTERT level by provoking miR-1182. It has recently been reported that circRNAs can act as miRNA sponges to negatively control miRNA. The majority of circRNAs can function as sponges, via a mechanism of back-splicing, as they are enriched in miRNA binding sites. They can also competitively bind to miRNAs and decrease the activity of miRNAs. Our results further showed that circRNAs can serve as competitive endogenous RNAs and play an important role in CRC development.
Exosomes, derived from the endosomal compartment, are released in the extracellular milieu under various physiological and pathological conditions [30]. It has been reported that exosomes are closely related to tumor development, including angiogenesis, metastasis, drug resistance, and immune escape [31], but the role of cancer-secreted exosomal circRNAs in CRC has not been clarified yet. Here, we performed TEM to reveal the shapes and size of exosomes from serum of CRC patients. Notely, we found that the highly expressed circFMN2 could be examined to serum exosomes of CRC patients. We found that circFMN2 levels are negatively correlated with miR-1182 levels in CRC patients’ serum.
Conclusions
In summary, our present study demonstrated that circFMN2 was up-regulated in CRC cell lines and tissues, and its high expression was associated with poor clinicopathological characteristics of CRC patients. Mechanistically, circFMN2 could significantly promote proliferation of CRC through directly binding to miR-1182 and subsequently reduce the suppressing capability of miR-1182 on hTERT. Exosome circFMN2 secreted from CRC patients’ serum play an important role to promote the tumor growth of CRC. Our study suggested that circFMN2 was a novel potential biomarker and therapeutic target in CRC.
Clinical perspectives
The molecular mechanisms that link circRNAs with CRC are not well understood.
CircFMN2 was an oncogenic factor that promotes tumorigenesis by serving as a ceRNA to regulate hTERT expression by sponging miR-1182.
CircFMN2 is a central component linking circRNAs to progression of CRC via an miR-1182/hTERT axis.
Availability of Data and Materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
Funding
This study was approved by National Natural Science Foundation of China [grant number 81872323] and the Science and Technology Department of Jilin Province [grant number 20190303149SF].
Author Contribution
Changfeng Li performed primers design and experiments. Yongchao Li and Ruisi Xu contributed flow cytometry assay and animal experiments. Dandan Li and Bin Zhang collected and classified the human tissue samples. Yongchao Li and Ruisi Xu contributed to qRT-PCR. Changfeng Li and Yun Wang analyzed the data. Changfeng Li wrote the paper. All authors read and approved the final manuscript.
Ethics Approval and Consent to Participate
This study was approved by the Medical Ethics Committee of China-Japan Union Hospital of Jilin University.
Consent for Publication
We have received consent from the individual patients who have participated in the present study. The consent forms will be provided upon request.
Competing Interests
The authors declare that there are no competing interests associated with the manuscript.
Abbreviations
- CCK-8
cell counting kit-8
- ceRNA
competing endogenous RNA
- circRNA
circular RNA
- COAD
colon adenocarcinoma
- CRC
colorectal cancer
- FC
fold change
- hTERT
human telomerase reverse transcriptase
- IHC
immunohistochemistry
- miRNA
microRNA
- qRT-PCR
quantitative real-time PCR
- PBS
phosphate-buffered saline
- RIP
RNA immunoprecipitation
- siRNA
small interfering RNA
- TCGA
the Cancer Genome Atlas
- TEM
transmission electron microscopy