Abstract

The incidence and mortality of colorectal cancer (CRC) are ranked in the top three worldwide in 2020. Abundant studies have reported that circular RNAs (circRNAs) act critical roles in the genesis and development of tumors, including CRC. Nevertheless, the roles and detailed regulation mechanisms of circRNAs that are related to the initiation and development of CRC have not been fully found and clarified. This research primarily revealed that circTMEM59 was greatly downregulated in CRC tissues and cell lines via qRT-PCR. In addition, the decreased expression of circTMEM59 was closely related to adverse clinicopathological characteristics and the shorter survival time of CRC patients. Then, a further study found that the overexpression of circTMEM59 suppressed cell growth and accelerated the cell death of CRC via a series of experiments in vitro and in vivo. Furthermore, circTMEM59 also repressed the metastatic behaviors of CRC cells. Further study revealed that circTMEM59 played the role of competing endogenous RNAs (ceRNAs) by binding to miR-668-3p to increase the expression of inhibitor of DNA binding 4 (ID4) in CRC. In summary, the results of this study clarified the antitumor effects of circTMEM59/miR-668-3p/ID4 axis in CRC progression and provided potential therapeutic targets and clinical prognostic markers for CRC.

1. Introduction

According to the report of Global Cancer Statistics 2020, colorectal cancer (CRC) is the most common gastrointestinal malignant tumor worldwide. Its incidence (about 188 million new cases) and mortality (about 91 million cases) are among the top three of all malignant tumors [1]. Approximately 555,000 new CRC patients and 286,000 deaths in China in 2020 are a burden to the whole society [2]. Although the survival rate of CRC patients greatly improved with the advances of comprehensive treatment programs, such as neoadjuvant therapy, radical surgery, postoperative radiotherapy, postoperative adjuvant chemotherapy, and targeting immunotherapy, the five-year survival rate of CRC patients with distant metastasis is only 14% [35]. Therefore, to improving the detection rate of early CRC patients [6], it is of great clinical significance and value to explore the detailed molecular mechanisms in the pathogenesis and development of CRC to strengthen the therapeutic effects on CRC patients.

Circular RNAs (circRNAs), which are noncoding RNAs with a loop structure that are more stable than linear transcripts, affect tumor growth, invasion and metastasis, angiogenesis, immune escape, and stemness by regulating the expression of target genes in epigenetic, transcriptional, and posttranscriptional processes [79]. Abundant research has investigated the key effects of circRNAs on the initiation and development of malignant tumors, especially the latent clinical significance of circRNAs in the therapeutic targets and diagnostic biomarkers of CRC [1013]. Accumulative studies have shown that circRNAs acted as competing endogenous RNAs (ceRNAs) via sponging microRNAs to significantly influence tumor progression, including cell growth, metastasis, and chemosensitivity [10, 14, 15]. Hence, exploring and clarifying the mechanisms of ceRNAs in emerging circRNAs can provide a new target for CRC therapeutics.

In general, this study first validated that circTMEM59 was downregulated in CRC tissues and cell lines. In addition, reduced circTMEM59 was correlated with adverse clinical characteristics and worse prognosis of CRC patients. Through multiple functional experiments in vitro and in vivo, we demonstrated that elevated circTMEM59 suppressed the growth and metastasis of CRC cells and accelerated cell apoptosis. With further investigation, we found that circTMEM59 acted as ceRNAs by competitively binding to miR-668-3p to modulate the expression of the inhibitor of DNA binding 4 (ID4), thereby influencing CRC progression. To conclude, this study first clarified the antitumor effects of circTMEM59/miR-668-3p/ID4 axis on CRC and provided potential therapeutic targets and clinical prognostic markers for CRC.

2. Materials and Methods

2.1. Tissue Specimen and Cell Culture

In this study, 100 paired CRC tissues and corresponding adjacent normal tissues were collected from the Department of Colorectal Surgery, the First Affiliated Hospital of Nanjing Medical University (NMU). None of the enrolled CRC patients received chemotherapy, radiotherapy, or immunotherapy before the surgery. All samples were retained at -80°C. All patients signed informed consent before surgery. The study was approved by the Ethics Committee of the First Affiliated Hospital of NMU (Ethics No.: 2020-SR-431).

HT29, LoVo, HCT116, SW480, DLD1, and NCM460 were purchased from the Cell Bank of Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China). HT29 was cultured in McCoy’s 5A medium (Thermo Fisher, Waltham, USA) supplemented with 10% fetal bovine serum (FBS) (Winsent, Quebec, Canada). RPMI 1640 (Thermo Fisher, Waltham, USA) supplemented with 10% FBS was used to culture the DLD1 cells. The cell culture mediums were supplemented with penicillin (100 U/ml) and streptomycin (100 μg/ml). All cells were placed at 37°C in an incubator containing a 5% concentration of CO2.

2.2. RNA Isolation, qRT-PCR, and RT-PCR for microRNAs

The nuclear and cytoplasmic RNAs were isolated with a PARIS™ kit (Thermo Fisher, Waltham, USA). The specific procedures of qRT-PCR and RT-PCR for microRNAs were consistent with the previous publication [16]. The primers for qRT-PCR of circRNA, mRNAs, microRNA, and internal control (U6 and β-actin) and primers for RT-PCR of microRNA are listed in Supplementary Table 1. The relative expression of genes was presented by using the 2ΔΔCT method.

2.3. Cell Transfection

Lentivirus of sh-circTMEM59, circTMEM59, miR-668-3p mimics, miR-668-3p inhibitor, and corresponding negative control (NC) were designed and synthesized by GenePharma (Shanghai, China). GenePharma also constructed pcDNA3.1 vectors containing ID4, sh-ID4, and the matched NC. The cells transfections were performed at 30%–40% confluence. Puromycin and neomycin were used to sift the stable transfected cells.

2.4. Cell Proliferation Assays

The specific procedures of Cell Counting Kit-8 (CCK-8), colony formation assay, and 5-ethynyl-2-deoxyuridine (EdU) were conducted as previously described [17]. The appropriate wavelength (450 nm) was chosen to detect the cell absorbance after adding CCK-8 reagents into cells for 2 hours in the CCK-8 assay. For colony formation assay, crystal violet was used to stain cell colonies after 600 CRC cells were cultivated in six-well plates for 10 days. The images were photographed by Nikon fluorescence microscope after all the experimental procedures were completed according to the manufacturers’ instructions in the EdU assay.

2.5. Transwell Assay

Specific procedures were consistent with the previous description [16]. Transwell inserts without Matrigel coating were used to measure the cell migratory ability. After adding Matrigel coating into inserts, cell invasion ability can be investigated. In each group, five random fields were selected and counted using a Nikon microscope (Nikon, Japan).

2.6. Cell Cycle and Apoptosis Analysis

The cells digested with trypsin were washed twice with PBS. Then, cells were fixed in 75% ethanol and stored at -20°C for at least 12 hours. Afterward, 75% ethanol was removed, and the cells were washed twice with PBS. An appropriate amount of DNA staining dye (propidium iodide, PI) (Vanzyme, Nanjing, China) was added into cells for half an hour in the dark at room temperature.

H2O2 (4 μmol) was added into cells for 2 hours to accelerate cell apoptosis before adding the apoptosis reagents (Annexin V-FITC and PI) (Vanzyme, Nanjing, China). BD FACSCanto II (BD Biosciences, San Jose, CA, USA) was used to analyze the cell cycle distribution and apoptotic rate of CRC cells.

2.7. Western Blot Analysis

The detailed procedures of western blot analysis were performed as previous description [16]. ECL (Millipore, Burlington, MA, USA) solution was used to display the level of protein expression by using the Bio-Imaging System. The details of primary antibodies used in the study are as follows: ki-67 (Abcam, ab15580), Bcl-2 (Abcam, ab32124), E-cadherin (Cell Signaling Technology, #14472), Bax (Proteintech, 50599-2-Ig), ID4 (Santa Cruz Biotechnology, sc-365656; Abcam, ab49261), caspase-3 (Cell Signaling Technology, #9662), CDK4 (Abcam, ab108357), Vimentin (Cell Signaling Technology, #5741), cyclin D1 (Abcam, ab 134175), and GAPDH (Abcam, ab ab9485).

2.8. Animal Models

BALB/c nude mice (male, 4 weeks old) were purchased from the Animal Core Facility of NMU (Nanjing, China). cells were subcutaneously injected on both armpits of each mouse. Then, the width and length of the tumor were measured and recorded every week. The mice were euthanized at 3rd week, and the tumors were striped, photographed, weighed, and fixed in formaldehyde in turn. Subsequently, these tumors were subjected to IHC and TUNEL detection.

Briefly, 106 cells suspended in 20 μL phosphate buffer saline were injected into the distal tip of the spleen of mice to construct the metastatic model of CRC cells. After injection, the spleen was removed. The liver tissues were dissected and fixed in formaldehyde for hematoxylin and eosin (H&E) staining 4 weeks later. The animal experiment was ratified by the Animal Ethics Committee of NMU.

2.9. Immunohistochemistry (IHC)

The specific methods of IHC were consistent with the previous description [16]. IHC scores were composed of two parts (staining intensity and positive percentage). The staining intensity was divided into four levels, namely, 0, 1, 2, and 3, corresponding to undyed, weak, medium, and strong staining, respectively. Similarly, positive percentage were divided into five levels, namely, 0, 1, 2, 3, and 4, corresponding to 0%–5%, 6%–25%, 26%–50%, 51%–75%, and 76%–100% positive staining cells, respectively. Five random fields were chosen to obtain the average IHC score.

2.10. Fluorescent In Situ Hybridization (FISH)

The detailed steps were in accordance with those performed in a previous publication [18]. CircTMEM59-specific Cy3-labeled probe was used to detect the subcellular localization of circTMEM59 in CRC cells by using a FISH Kit (RiboBio, Guangzhou, China). The subcellular localization of miR-668-3p in CRC cells was detected by miR-668-3p detection probe (GenePharma, Shanghai, China).

2.11. Luciferase Reporter Assay

The sequences matched the 3-UTR of circTMEM59 and ID4 mRNA, and those containing the wild-type or mutant miR-668-3p binding sites were designed and synthesized by GenScript (Nanjing, China). Procedures were conducted as the previous description [16].

2.12. RNA Immunoprecipitation (RIP)

The RIP Kit (Millipore, Burlington, MA, USA) was used to conduct the RIP assays according to the manufacturer’s instructions. We extracted the coimmunoprecipitated RNA from the precipitates for further qRT-PCR. RIP assay was performed as previous description [17].

2.13. TUNEL Assay

TUNEL assay was performed as previous description [19]. The nucleus of the cell was stained by DAPI. The fluorescence microscope (Leica, Germany) was used to obtain images, which were then randomly filed.

2.14. Statistical Analysis

Each experiment was conducted independently thrice in this study. SPSS 22.0 (Chicago, USA) and GraphPad Prism 8.0 (CA, USA) software were used for the processing and analyzing data. Chi-square test, -test, two-way ANOVA, Pearson correlation analysis, and Kaplan-Meier analysis were performed based on the experimental requirements. In each analysis, the significance threshold was set to 0.05. The data were shown as .

3. Results

3.1. CircTMEM59 Is Downregulated in CRC Tissues and Cell Lines and Correlates with Adverse Clinical Characteristics

To investigate the key circRNAs involved in CRC advancement, the GEO database (GSE142837) was employed to explore the differential expression of circRNAs in CRC tissues compared with adjacent tissues. The expression of circTMEM59 (circBase ID: hsa_circ_0012634) was remarkably decreased in CRC tissues than in adjacent samples (Figure 1(a)). Furthermore, downregulated cireTMEM59 was confirmed in the matched 100 CRC tissues and 100 adjacent normal tissues (Figure 1(b)). We divided patients into two groups (larger and smaller tumor groups) according to the length of the CRC tumors. As shown in Figure 1(c), circTMEM59 was higher in smaller tumors than in larger tumors. Moreover, the results of circTMEM59 in CRC cell lines and NCM460 indicated that circTMEM59 expression in normal colorectal mucosa epithelial cells was higher than in the five CRC cell lines (Figure 1(d)). In addition, qRT-PCR and FISH assays confirmed that circTMEM59 was largely located in the cytoplasm of CRC cells (Figures 1(e) and 1(f)). RNase R was added into the CRC cells to authenticate the circular structure of circTMEM59, and qRT-PCR showed that circTMEM59 was more resistant to RNase R than linearTMEM59 mRNA in CRC cells (Figure 1(g)). Besides, the results of gel electrophoresis presented that the convergent primers of circTMEM59 and GAPDH could amplify the expected size products from cDNA and genomic DNA (gDNA). Only the divergent primers of circTMEM59 could amplify the PCR products from cDNA but not from gDNA. The divergent primers of GAPDH could not amplify the PCR products from cDNA or gDNA (Figure 1(h)). The abovementioned results showed that circTMEM59’s structure was circular rather than linear. In addition, the relationship between the abundance of circTMEM59 and clinical characteristics and prognosis was explored via the collection and analysis of the data from 100 CRC patients. The data showed that the patients with lower expression of circTMEM59 had a worse prognosis (OS and DFS) than patients with higher expression (Figure 1(i)). Furthermore, circTMEM59 expression was significantly and negatively correlated with tumor size, TNM staging system, tumor stage, and the metastasis of surrounding lymph node (Table 1). It revealed that aberrant circTMEM59 was closely related to the proliferation and metastasis of CRC. Apart from that, the results of IHC also indicated that tumors with low expression of circTMEM59 grew faster and expressed lower levels of ID4 in comparison with tumors with high expression of circTMEM59 or the adjacent tissues (Figure 1(j)). The abovementioned facts suggested that circTMEM59 was downregulated in CRC; it can serve as clinical biomarkers for CRC.

3.2. CircTMEM59 Represses Cell Cycle Transition and Facilitates Cell Apoptosis and Metastasis of CRC In Vitro

HT29 and DLD1 cell lines were chosen to conduct the functional experiments of circTMEM59 via the different abundance of circTMEM59 in CRC cell lines. As shown in Figure 2(a), qRT-PCR affirmed that circTMEM59 was downregulated in HT29 cells and overexpressed in DLD1 cells. Downregulated circTMEM59 promoted the cell growth of CRC via CCK-8 assay, whereas elevated circTMEM59 significantly repressed the cell proliferation of CRC (Figure 2(b)). Similarly, the results of colony formation assay and EdU assay also verified that the reduction of circTMEM59 facilitated cell proliferation, and the overexpression of circTMEM59 restrained the cell growth of CRC (Figures 2(c) and 2(d)). In addition, the overexpression of circTMEM59 arrested the cell cycle transition in the G0-G1 phrase, whereas reduction of circTMEM59 propelled the CRC cells transition to the S phrase (Figure 2(e)). Furthermore, increased circTMEM59 accelerated the cell apoptosis of CRC. However, apoptotic rates of HT29 were significantly decreased with the reduction of circTMEM59 (Figure 2(f)). Transwell assay demonstrated that migration and invasion abilities were greatly improved when circTMEM59 were downregulated in HT29 cells. On the contrary, overexpressed circTMEM59 dramatically suppressed migrating cell number and invasion of DLD1 (Figure 2(g)). The photographs of IHC showed that CRC tissues with a low level of circTMEM59 had higher expressions of Bcl-2, CDK4, Vimentin, and cyclin D1 and lower expressions of ID4, caspase-3, and E-cadherin compared with the high level of circTMEM59 expression in tumors of CRC (Figure 2(h) and Figure S1A). These data demonstrated that the overexpression of circTMEM59 prohibited CRC cells progression by preventing the cell cycle from entering the S phase, promoting cell apoptosis, and restraining cell metastasis in vitro.

3.3. CircTMEM59 Inhibits CRC Progression by Acting as a Molecular Sponge of miR-668-3p

Previous research has reported that circRNAs located in the cytoplasm mainly exerted its effects on tumor progression by endogenously competing with microRNAs [10]. FISH and qRT-PCR had confirmed the cytoplasm location of circTMEM59. miRDB (http://mirdb.org/) and circBank (http://www.circbank.cn/index.html) were analyzed to seek the potential microRNAs that circTMEM59 bound to. Hsa-miR-668-3p was chosen from candidates that probably interacted with circTMEM59 through integrating and analyzing the screening data (Figure 3(a)). Then, a series of experiments were conducted to verify the hypothesis. First, overexpressed miR-668-3p significantly suppressed the luciferase activity of wild-type circTMEM59 and reduced miR-668-3p elevated the luciferase activity of wild-type circTMEM59 in CRC cells, whereas no significance was detected in the cells of mutant circTMEM59, as shown by luciferase reporter assay (Figures 3(b) and 3(c)). Second, the FISH assay validated the colocation of circTMEM59 and miR-668-3p in the cytoplasm of CRC cells (Figure 3(d)). Third, the expression of miR-668-3p was remarkably increased in CRC tissues than in adjacent samples (Figure 3(e)). Moreover, miR-668-3p increased in bigger tumors than that in smaller tumors (Figure 3(e)). In addition, miR-668-3p expression was significantly and positively correlated with tumor size, TNM staging system, tumor stage, and metastasis of the surrounding lymph nodes in CRC patients (Table 1). The expression of miR-668-3p was significantly and negatively correlated with the level of circTMEM59 expression in 100 CRC specimens (Figure 3(f)). The results of IHC showed that the CRC tissues with high miR-668-3p presented more proliferation and metastasis-related proteins (upregulated CDK4, cyclin D1, ki-67, and Vimentin and downregulated ID4 and E-cadherin) and less apoptosis-related proteins (upregulated Bcl-2 and downregulated caspase-3) than CRC tissues with low miR-668-3p (Figure 3(g)). Collectively, these data revealed that a competing endogenous RNA (ceRNA) relationship existed between circTMEM59 and miR-668-3p in CRC.

3.4. circTMEM59 Exerts Its Effects on CRC Cells by Regulating miR-668-3p

A series of rescue experiments were conducted to investigate the biological role of miR-668-3p interacting with circTMEM59 in CRC. In HT29 cells, cotransfected miR-668-3p-inhibitor repressed the sh-circTMEM59-induced miR-668-3p elevation. In DLD1, the expression of miR-668-3p was reversed in the cotransfected miR-668-3p-mimics cells compared with the upregulated circTMEM59 group (Figure 4(a)). Downregulated miR-668-3p restrained the promotion of cell growth by circTMEM59 reduction, whereas overexpressed miR-668-3p abolished the repressed effects of elevated circTMEM59 on cell growth of CRC, as determined by CCK-8, EdU, and colony formation (Figures 4(b)4(e)). As shown in Figure 4(f), decreased miR-668-3p expression abolished the effects of circTMEM59 reduction on cell cycle promotion to the S phrase, and increased miR-668-3p expression reverted the effects of circTMEM59 overexpression on cell cycle arrest in the G0–G1 phrase. Similar results were observed in the apoptosis assay in the cotransfected cells, decreased miR-668-3p expression abolished the effects of circTMEM59 reduction on cell apoptosis inhibition, and increased miR-668-3p expression reverted the effects of circTMEM59 overexpression on cell apoptosis acceleration (Figure 4(g)). In addition, the migration and invasion cell number decreased in the cotransfected sh-circTMEM59 + miR-668-3p-inhibitor HT29 cells compared with the cotransfected sh-circTMEM59 + Ctl HT29 cells, whereas the metastatic cell number increased in the cotransfected circTMEM59 + miR-668-3p-mimics DLD1 cells compared with the cotransfected circTMEM59 + NC DLD1 cells (Figure 4(h)). These data suggested that the suppression of CRC cell proliferation and metastasis by circTMEM59 is regulated by miR-668-3p.

3.5. CircTMEM59 Represses Cell Proliferation and Metastasis of CRC In Vivo

To investigate the biological functions of circTMEM59 in vivo, we chose a subcutaneous mouse model that is more easily and accurately monitored, although the orthotopic mouse model can better provide the anatomical site and microenvironment of CRC [20]. Decreased circTMEM59 HT29 cells (HT29-sh-circTMEM59) and negative control cells (HT29-NC) or elevated circTMEM59 DLD1 cells (DLD1-circTMEM59) and negative control cells (DLD1-Ctl) were subcutaneously injected into the nude mice. As shown in Figure 5(a), downregulated circTMEM59 facilitated the tumor proliferation of CRC cells. Furthermore, the tumors of the sh-circTMEM59 group grew much larger and heavier compared with those in the NC group (Figures 5(b) and 5(c)). Through the TUNEL and IHC assays, downregulated circTMEM59 promoted cell proliferation and inhibited cell apoptosis by increasing Bcl-2, CDK4, cyclin D1, and ki-67 expressions and decreasing ID4, Bax expressions, and TUNEL-positive cells (Figure 5(d) and Figure S1B). Conversely, upregulated circTMEM59 significantly suppressed the tumor growth of CRC cells (Figure 5(e)). In addition, the tumors in the DLD1-circTMEM59 group were much smaller and lighter than those in the DLD1-Ctl group (Figures 5(f) and 5(g)). Upregulated circTMEM59 resulted in tumors with less proliferation-related proteins (decreased CDK4, cyclin D1, ki-67, and increased ID4) and more apoptosis-related proteins (increased Bax and TUNEL-positive cells and decreased Bcl-2) (Figure 5(h) and Figure S1C). The circTMEM59 knockdown group showed the most number of and largest foci in the liver of nude mice (Figure 5(i)). Inversely, overexpressed circTMEM59 in DLD1 cells led to a massive decrease in size of foci and number of liver metastasis nodules (Figures 5(j)). The abovementioned data presented that circTMEM59 suppressed the proliferative and metastatic activities of CRC and regulated the level of proliferation-related proteins, apoptosis-related proteins, and ID4 expression in vivo.

3.6. ID4, a Target of miR-668-3p, Is Indirectly Modulated by circTMEM59

To further investigate the ceRNA network among circTMEM59, miR-668-3p, and downstream targets in CRC, miRWalk (http://mirwalk.umm.uni-heidelberg.de/), miRDB (http://mirdb.org/), and TargetScan (http://www.targetscan.org/vert_72/) databases were employed to explore the direct targets of miR-668-3p (Figure 6(a)). As presented in Figure 6(b), a conserved site of 3-UTR of ID4 mRNA was identified to be directly bound to miR-668-3p. Then, luciferase reporter assay (Figure 6(c)) and RIP assay (Figure 6(d)) were conducted to validate the hypothesis that miR-668-3p directly binds to 3-UTR of ID4 mRNA. The expression of ID4 was measured in the 100 paired CRC tissues and adjacent normal tissues by qRT-PCR. As expected, ID4 was simultaneously decreased in CRC tissues and larger tissues compared with that in adjacent tissues and smaller tissues (Figure 6(e)). Pearson’s regression analysis revealed that the ID4 expression was a positive correlation with circTMEM59 expression, whereas it was a negative correlation with miR-668-3p expression in 100 CRC tissues (Figure 6(f)). Thereafter, immunoblot assays validated that elevated miR-668-3p decreased the expressions of ID4 protein in the CRC cells, whereas reduced miR-668-3p increased them (Figure 6(g)). Reduced circTMEM59 decreased the expressions of ID4 protein in the CRC cells, and elevated circTMEM59 increased them (Figure 6(h)). Furthermore, downregulated miR-668-3p reversed the reduction of ID4 protein resulting from reduced circTMEM59 (Figure 6(i)). Overexpressed miR-668-3p inverted the elevation of ID4 expression caused by elevated circTMEM59 (Figure 6(j)). IHC results showed that the CRC tissues with low ID4 expression presented more proliferation and metastasis-related proteins (upregulated CDK4, cyclin D1, ki-67, and Vimentin and downregulated E-cadherin) and less apoptosis-related proteins (upregulated Bcl-2 and downregulated Caspase-3) than CRC tissues with high ID4 expression (Figure 6(k)). These data demonstrated that circTMEM59 repressed cell growth and metastasis and facilitated cell death through competing binding to miR-668-3p to accumulate ID4 expression of CRC.

3.7. circTMEM59 Exerts Its Effects on CRC Cells by Indirectly Regulating ID4

Previous research has reported that ID4 acts as a tumor suppressor in CRC progression [21, 22]. We conducted the rescue experiments to confirm the interaction between circTMEM59 and ID4 of CRC. Immunoblot assays validated that elevated ID4 reversed the reduction of ID4 protein caused by reduced circTMEM59 (Figure 7(a)). In addition, downregulated ID4 inverted the ID4 increment, which resulted from elevated circTMEM59 (Figure 7(a)). Moreover, upregulated ID4 restrained the promotion of cell growth by circTMEM59 reduction, whereas reduced ID4 abolished the repressive effects of elevated circTMEM59 on CRC cell growth, as shown in CCK-8, EdU, and colony formation assays (Figures 7(b)7(d), Figures S1D and S1E). Similar results were observed in the apoptosis assay, in which the cotransfected cells that upregulated ID4 abolished the effects of circTMEM59 reduction on cell apoptosis inhibition, and the reduced ID4 reverted the effects of circTMEM59 overexpression on cell apoptosis acceleration (Figure 7(e), Figure S1F). As presented in Figure 7(f) and Figure S1G, upregulated ID4 abolished the effects of circTMEM59 reduction on cell cycle promotion, and the reduced ID4 reverted the effects of circTMEM59 overexpression on cell cycle arrest. In addition, the migration and invasion cell number decreased in the cotransfected sh-circTMEM59 + ID4 HT29 cells compared with the cotransfected sh-circTMEM59 + Ctl HT29 cells, whereas the metastatic cell number increased in the cotransfected circTMEM59 + sh-ID4 DLD1 cells compared with the cotransfected circTMEM59 + NC DLD1 cells (Figure 7(g), Figure S1H). Western blot results demonstrated that downregulated miR-668-3p abolished the effects of circTMEM59 reduction on the elevation of Bcl-2, CDK4, cyclin D1, and Vimentin and the depletion of cleaved caspase-3 and E-cadherin. Upregulated miR-668-3p reverted the effects of circTMEM59 overexpression on the depletion of Bcl-2, CDK4, cyclin D1, and Vimentin and the elevation of cleaved caspase-3 and E-cadherin (Figure 7(h)). Furthermore, upregulated ID4 abolished the effects of circTMEM59 reduction on the elevation of Bcl-2, CDK4, cyclin D1, and Vimentin and the depletion of cleaved caspase-3 and E-cadherin. Downregulated ID4 reverted the effects of circTMEM59 overexpression on the depletion of Bcl-2, CDK4, cyclin D1, and Vimentin and the elevation of cleaved caspase-3 and E-cadherin (Figure 7(i)). To conclude, the abovementioned results confirmed that upregulated circTMEM59 inhibits cell growth and metastasis by miR-668-3p/ID4 axis in CRC (Figure 8).

4. Discussion

CircRNAs, which are emerging noncoding RNAs with loop structure, are widely reported in connection with the malignant progression of carcinoma, including sustaining growth, avoiding immune destruction, activating metastasis, and reprograming metabolism [11, 14, 2325]. Reactive oxygen species (ROS) are natural byproducts of oxidative phosphorylation on the electron transport chain in mitochondria. The imbalance between ROS production and degradation can cause oxidative stress, resulting in energy depletion, DNA damage, protein carbonylation, and lipid oxidation, leading to various diseases and disorders [26, 27]. The close relationship between circular RNAs and ROS suggests that circular RNAs play critical roles as a target in the treatment of tumor metabolism [2830]. Focusing on circRNAs as emerging molecular propellants of CRC progression, we intended to determine the mechanisms underlying circular RNA-mediated development of CRC. In this study, circTMEM59 was first explored in CRC progression as a tumor suppressor. CircTMEM59’s downregulated expression in CRC tissues and cell lines was verified (Figure 1). The aberrant profiles of noncoding RNAs in the tumor tissues or the circulation can also be used to predict the clinicopathological characteristics and long-term survival of tumor patients [16, 31, 32]. In this study, CRC patients with high expression of circTMEM59 presented better clinicopathological characteristics (smaller tumor size, shallower tumor infiltration depth, and less lymph node metastasis) and longer survival time (Table 1 and Figure 1). These integral results showed that circTMEM59 acted as a tumor inhibitor in CRC and can be used as a clinical biomarker for diagnosis and prognosis.

The aberrant regulations of the cell cycle, cell apoptosis, and metastatic behaviors have played critical roles in the proliferation and progression of tumors [19, 31, 33, 34]. Some studies revealed that circRNAs play key roles in the regulation of the cell cycle, cell apoptosis, and distant metastases of CRC [11, 14, 15]. A series of functional assays validated that the upregulated circTMEM59 inhibited cell growth and cell metastasis of CRC, whereas downregulated circTMEM59 accelerated the cell cycle shift and metastatic behaviors and repressed apoptosis of CRC cells in vitro and in vivo (Figures 2 and 5). Combined with clinical data, IHC results showed that CRC tissues with a low level of circTMEM59 had higher expressions of Bcl-2, CDK4, cyclin D1, and Vimentin and lower expressions of ID4, caspase-3, and E-cadherin compared with the high level of circTMEM59 in CRC tissues (Figure 2). The abovementioned results explore the potentiality and orientation of manipulating circTMEM59 for CRC therapy, combined with the effective application of noncoding RNA-delivering systems [35].

Recent studies reported that circRNA functions as ceRNA of microRNAs in the cytoplasm to affect the tumor progression [10, 15]. qRT-PCR and FISH assays confirmed that circTMEM59 was largely located in the cytoplasm of CRC cells. miR-668-3p was identified as a target of circTMEM59 by bioinformatic tools and luciferase reporter assay, which is reported to be an oncogenic gene in hepatocellular carcinoma [36]. The miR-668-3p was elevated in 100 CRC tissues and positively correlated with adverse clinical features, according to the analysis of clinical specimens and data (Figure 3). Moreover, miR-668-3p was negative relation to circTMEM59 expression in 100 CRC tissues. IHC images presented that CRC tissues with high miR-668-3p presented more proliferation and metastasis-related proteins (upregulated CDK4, cyclin D1, ki-67, and Vimentin and downregulated ID4 and E-cadherin) and less apoptosis-related proteins (upregulated Bcl-2 and downregulated caspase-3) than CRC tissues with low miR-668-3p expression (Figure 3). In addition, miR-668-3p mediated the effects of circTMEM59 on the cell cycle, cell apoptosis, and metastasis of CRC. Decreased miR-668-3p expression abolished the effects of circTMEM59 reduction on cell proliferative and metastatic promotion and cell apoptosis inhibition. Increased miR-668-3p expression reverted the effects of circTMEM59 overexpression on cell proliferative suppression, metastatic inhibition, and cell apoptosis acceleration (Figure 4). The intrinsic ceRNA mechanisms between circTMEM59 and miR-668-3p were gradually revealed in CRC.

To further determine the detailed ceRNAs mechanisms of circTMEM59 in CRC, three datasets were employed to seek the key targets of miR-668-3p. ID4 was identified as a target of miR-668-3p through luciferase reporter and RIP assays, which was abnormal in different types of cancer [3739]. Importantly, ID4 is reportedly an inhibitor of CRC progression [21, 22]. As expected, ID4 was simultaneously decreased in CRC tissues and larger tissues compared with that in adjacent tissues and smaller tissues. Moreover, ID4 was positively correlated with circTMEM59 and was negatively correlated with miR-668-3p in 100 CRC tissues (Figure 6). In addition, miR-668-3p restrained and circTMEM59 facilitated the abundance of ID4 proteins in CRC cells. Through the results of IHC, the CRC tissues with low ID4 expression presented more proliferation and metastasis-related proteins (upregulated CDK4, cyclin D1, ki-67, and Vimentin and downregulated E-cadherin) and less apoptosis-related proteins (upregulated Bcl-2 and downregulated caspase-3) than CRC tissues with high ID4 expression (Figure 6). Rescue experiments revealed that upregulated ID4 restrained the promotion in cell growth and metastasis of circTMEM59 reduction, whereas reduced ID4 abolished the repressed effects of increased circTMEM59 on cell growth, migration, and invasion of CRC (Figure 7). Taken together, all data suggested that elevated circTMEM59 prohibited cell cycle transition (CDK4 and cyclin D1), inhibited cell metastasis (E-cadherin and Vimentin) and promoted cell apoptosis (cleaved caspase-3 and Bcl-2) by accumulating ID4 by sponging for miR-668-3p in CRC (Figure 8).

5. Conclusion

In general, this study validated for the first time that circTMEM59 was downregulated in CRC tissues and cell lines. In addition, low circTMEM59 expression in CRC patients was closely related to shorter survival time and advanced clinicopathological characteristics. Furthermore, circTMEM59 repressed the cell growth and metastasis and accelerated cell apoptosis of CRC via accumulating ID4 by sponging miR-668-3p. In conclusion, the aberrant circTMEM59 had profound effects on tumor growth, migration, and invasion and was an effective biomarker and therapeutic target for CRC.

Data Availability

The datasets used in the study are available from the corresponding author on reasonable request.

Ethical Approval

The clinical data and specimens used in this study were approved by the Ethics Committee of the First Affiliated Hospital with NMU. The animal experiment was ratified by the Animal Ethics Committee of NMU.

Conflicts of Interest

The authors declare that they have no potential competing interests.

Authors’ Contributions

Yueming Sun and Yifei Feng generated and constructed this research. Yifei Feng, Xiaoyuan Wang, and Changzhi Huang performed a series of functional experiments. Yifei Feng, Dongsheng Zhang, and Chuan Zhang established the animal model. Yifei Feng, Yue Zhang, and Dongjian Ji conducted and repeated the western blot assays. Yifei Feng and Junwei Tang analyzed the data. Yueming Sun supervised the experiments, secured funding, and interpreted results. The authors read and approved the final manuscript. Yifei Feng, Xiaoyuan Wang, Changzhi Huang, and Dongsheng Zhang contributed equally to this work.

Acknowledgments

This research got funding from the Natural Science Foundation of Jiangsu Province (Grant No. BK20211378), National Basic Research Program of China (973 Program) (2017YFC0908200), and Jiangsu Key Medical Discipline (Grant No. ZDxKA2016005).

Supplementary Materials

Supplementary Material Supplementary Figure 1: the regulation of circTMEM59 on CRC cells is mediated by ID4. A The IHC scores of ID4, Bcl-2, Caspase-3, CDK4, cyclin D1, E-cadherin, and Vimentin in the low expression circTMEM59 of tumors and high expression circTMEM59 of tumors. B The IHC scores of ID4, Bcl-2, Bax, CDK4, cyclin D1, and ki-67 and TUNEL data in NC group and sh-circTMEM59 group. C The IHC scores of ID4, Bcl-2, Bax, CDK4, cyclin D1, and ki-67 and TUNEL data in Ctl group and circTMEM59 group. D Effects of circTMEM59 and ID4 on proliferation in cotransfected CRC cell lines were detected by EdU assays (scale bar: 100 µm). E Effects of circTMEM59 and ID4 on proliferation in cotransfected CRC cell lines were detected by colony formation. F ID4 overexpression facilitated cell apoptosis of circTMEM59-downregulating of HT29 cells, while ID4 reduction inhibited cell apoptosis of circTMEM59-overexpressing of DLD1 cells. G Effects of circTMEM59 and ID4 on regulating cell cycle in cotransfected CRC cell lines. H The migration and invasion abilities were detected by Transwell assays in cotransfected CRC cell lines (scale bar: 100 µm). Supplementary Table 1: the primers for qRT-PCR and RT-PCR. (Supplementary Materials)