1. CircRNA_0009092 was decreased in CRC tissues and cell lines.
In order to identify the circRNA expression profile in CRC, the circRNA microarray analysis was performed to identify the aberrant circRNA expression in three paired CRC tissues and adjacent normal tissues (ANT). The results demonstrated that 20 circRNAs were differentially expressed, including 10 upregulated and 10 downregulated circRNAs. The heatmap displayed the most upregulated and downregulated circRNAs (Fig. 1A). In the present study, we selected circ_0009092 as our protagonist. Firstly, the expression level of circ_0009092 in 80 paired CRC tissues and ANT was determined. As shown in Fig. 1B, the expression level of circ_0009092 was significantly downregulated in CRC tissues. Notably, circ_0009092 expression was correlated with TNM stage, lymph node involvement and metastasis (Table S1). Next, we investigated the correlation between circ_0009092 expression and the prognosis of CRC patients. Kaplan–Meier survival analysis showed that a lower level of circ_0009092 expression was associated with poor overall survival (OS) and recurrence-free survival (RFS) (Fig. 1C-D). Based on the above results, we postulated that circ_0009092 is a key circRNA involved in CRC progression.
CircRNAs are characterized by a covalently closed-loop structure that is generated by back-splicing, eventually leading to increased structural stability. To investigate these properties of circ_0009092, we first determined the circ_0009092 expression in CRC cells and normal cells. Results showed that circ_0009092 was highly expressed in Lovo cells and expressed at low levels in HCT-116 cells, then we chose Lovo and HCT-116 cell lines for subsequent experiments. Sanger sequencing and RNase R treatment were used in HCT-116 and Lovo cells. According to circbase, circ_0009092 was derived from exons 19 and 20 of OGDH gene with a length of 202 nucleotides (Fig. 1E). The presence of head-to-tail back-splicing junction in the circ_0009092 was confirmed by Sanger sequencing in the qPCR product of circ_0009092 (Fig. 1F). We then treated tumor cells with RNase R (a degrader of the linear RNA) and observed that circ_0009092 was resistant to RNase R digestion compared with the linear OGDH mRNA, suggesting that circ_0009092 was more stable than liner OGDH mRNA (Fig. 1G, Fig.S1B). In order to certify circ_0009092 was a circRNA, we designed two sets of primers in qPCR: divergent primers, amplifying circular form, and convergent primers, amplifying linear form. Results showed that circ_0009092 could be amplified by divergent primers in cDNA but not genomic DNA (Fig. S1C). We further investigated the subcellular location of circ_0009092 by nuclear/cytoplasmic fractionation and FISH assay. The results revealed that circ_0009092 was predominately expressed in the cytoplasm of CRC cells (Fig. 1H-I). Taken together, these data demonstrated that circ_0009092 expression was strikingly downregulated in CRC tissues and was a stable circRNA in CRC cells.
2. Circ_0009092 inhibited proliferation, migration, and invasion of CRC cells.
Given the results that circ_0009092 expression is decreased in CRC tissues, we speculated that circ_0009092 may function as a tumor suppressor in CRC. The cell proliferation assay, the wound healing assay, and the Transwell assay with or without Matrigel were conducted to assess the effect of circ_0009092 on cell proliferation, migration, and invasion. Edu assay showed that overexpression of circ_0009092 significantly inhibited the cell proliferation of CRC cells, while knockdown of circ_0009092 markedly promoted cell proliferation (Fig. 2A-B). In a clonogenic assay, overexpression of circ_0009092 lowered both the number of clones, and knockdown of circ_0009092 enhanced clonal proliferative capacity of CRC cells (Fig. 2C). CircRNAs play important roles in enhancing cancer migration and invasion through regulating EMT. Therefore, the effect of circ_0009092 on the expression of EMT-related protein was determined by performing Western blot (WB) assay. WB results showed that decreased circ_0009092 significantly increased the Vimentin protein level and suppressed the expression of E-cadherin. On the contrary, there was a decreased Vimentin protein level and an increased E-cadherin expression level in the circ_0009092 overexpressed group (Fig. 2D). Next, the wound healing assay showed that overexpression of circ_0009092 significantly impaired the migration of CRC cells. In the circ_0009092 knockdown group, the wound healing assay showed the contrary tendency (Fig. 2E). In addition, we explored the effects of circ_0009092 on cell migration and invasion by Transwell migration and Transwell invasion assays. Results showed overexpression of circ_0009092 inhibited cell migration and invasion, while knockdown of circ_0009092 enhanced cell migration and invasion (Fig. 2F). Collectively, these data demonstrated that circ_0009092 inhibited proliferation, migration, and invasion of CRC cells.
3. Circ_0009092 functioned as a sponge for miR-665
CircRNA exerted diverse biological functions by acting as microRNA “sponges”. Therefore, we assessed whether circ_0009092 could sponge miRNA to suppress CRC progression. Five potential miRNAs (miR-330, miR-549, miR-665, miR-671, miR-885) were predicted by the Circinteractome and starBase database. We analyzed the expression levels of these five miRNAs in CRC cells and found that miR-665 was markedly correlated with circ_0009092 (Fig. 3A). qPCR results indicated that the expression of miR-665 was remarkably downregulated by circ_0009092 knockdown, while circ_0009092 overexpression resulted in upregulation of miR-665 in CRC cells (Fig. 3B). In addition, miR-665 expression levels were markedly overexpressed (Fig. 3C) and negatively correlated with circ_0009092 expression in CRC tissues (Fig.3D, R=-0.6111, P < 0.001). Kaplan-Meier survival analysis showed that high expression of miR-665 was correlated with worse RFS and OS in CRC patients (Fig. 3E, Fig. S2B). In addition, the FISH assay was applied to assess the subcellular co-location of circ_0009092 and miR-665. The results demonstrated that circ_0009092 and miR-665 were colocalized in the cytoplasm of HCT-116 and Lovo cells, suggesting that circ_0009092 could interact with miR-665 and function in the cytoplasm (Fig. 3F). FISH analysis of CRC samples and ANT confirmed the downregulation of circ_0009092 and upregulation of miR-665 in CRC samples compared to ANT (Fig. S2A).
We further explored whether circ_0009092 exerts its regulatory effects by sponging miR-665. Bioinformatic analysis was used to predict the potential binding sites of miR-665 on circ_0009092 (Fig. 3G). To confirm the results, we performed RNA immunoprecipitation (RIP) for AGO2 in tumor cells and investigated the expression of circ_0009092 and miR-665 by qPCR. Circ_0009092 and miR-665 were highly expressed in the AGO2 pellet compared with those in the IgG control. Moreover, circ_0009092 was highly enriched in the miR-665-overexpressed (miR-665 mimic) group compared with that in the negative control (NC) group (Fig. 3H). In addition, dual-luciferase activity analysis was performed to verify that miR-665 binds to circ_0009092 in CRC cells. The sequences of circ_0009092 containing wild or mutant miR-665 binding site were cloned into the luciferase reporter vector pLG3. The results demonstrated that miR-665 mimic significantly reduces the luciferase activity of the WT circ_0009092, but had no effect on that of the MUT circ_0009092 compared with the NC group (Fig. 3I). Taken together, these results indicated that the interaction of circ_0009092 and miR-665 lead to the down expression of miR-665 in CRC cells.
4. Circ_0009092 inhibited CRC cell progression by targeting miR-665.
To explore that circ_0009092 regulates tumor progression by sponging miR-665, miR-665 knockdown and overexpression cell lines were constructed. The results of EdU (Fig. 4A) and colony formation assay (Fig. 4B) revealed that miR-665 significantly promoted and miR-665 inhibitor markedly inhibited the proliferation of CRC cells. Wound healing (Fig. 4C) and Transwell assay (Fig. 4D-E) showed that the migratory and invasive abilities of CRC cells were distinctly promoted with the transfection of miR-665 mimic and suppressed by miR-665 inhibitor. To further investigate the roles of circ_0009092/miR-665 in CRC progression, we performed a rescue assay. EdU (Fig. 4A) and colony formation assay (Fig. 4B) revealed that miR-665 mimic efficiently restored the proliferative capacity of circ_0009092-overexpressed cells. MiR-665 inhibitor abrogated the increased proliferation of circ_0009092-deleted CRC cells (Fig. S3A-B). Similarly, wound healing (Fig. 4C, Fig. S3C) and Transwell assay (Fig. 4D-E, Fig. S3D-E) showed that the promotive effects on migratory and invasive abilities induced by circ_0009092 knockdown were restored by miR-665 mimic. The suppressive effect induced by circ_0009092 overexpression was reversed by miR-665 inhibitor. Together, these data indicated that circ_0009092 regulates CRC progression via sponging miR-665.
5. NLK was the target of miR-665 in CRC cells
Based on TargetScan, miRWalk, and Starbase databases, NLK, PRDM15, PRSS18, and PPP2R2A were identified as the potential targets of miR-665. Overexpression of circ_0009092 induced the expression of NLK mRNA in HCT-116 and Lovo cells, while knockdown of circ_0009092 exerted the opposite effects. However, exogenous expression or knockdown of circ_0009092 had no effect on the expression of PRDM15, PRSS18, and PPP2R2A (Fig. 5A, Fig.S4A). In addition, miR-665 mimic significantly suppressed the expression of NLK mRNA and miR-665 inhibitor markedly increased NLK mRNA expression in HCT-116 and Lovo cells (Fig. 5B). Similarly, the WB results validated that NLK expression was regulated by circ_0009092 and miR-665 in CRC cells (Fig. 5C-D). We further examined the expression of NLK in tumor tissues of CRC. The NLK mRNA levels were significantly lower in 80 tumor samples than that in ANT (Fig. 5E). Furthermore, the expression level of NLK in CRC tissues was positively correlated with circ_0009092 expression and negatively correlated with miR-665 expression (circ_0009092 p < 0.001, miR-665 p < 0.001, Fig. 5F, Fig.S4B). Then, the binding ability of miR-665 to NLK mRNA was detected by RIP assay using the antibody against AGO2, a component of the RNA-induced silencing complex that mediates miRNA-directed gene silencing. Overexpression of miR-665 significantly enriched the association of NLK mRNA with the Ago2-containing complex, suggesting that the NLK mRNA was directly bound by miR-665 (Fig. 5H). We further found that miR-665 mimics could reduce the luciferase activity of the WT NLK- 3’ -UTR, but had no effect on that of the Mut NLK-3’-UTR compared with the NC group in HCT-116 and Lovo cells (Fig. 5I, Fig.S4C).
It has been reported that NLK is an effector of the Wnt/β-catenin signaling pathway, we speculated that circ_0009092/miR-665/NLK regulates the Wnt/β-catenin signaling pathway. Western blot showed that overexpression of circ_0009092 markedly attenuated the expression of EMT-related proteins, β-catenin, c-myc (target of Wnt/β-catenin) and p-GSK-3β (Fig. 5J), while circ_0009092 knockdown achieved an opposite result (Fig. 5K). We also determined the expression of EMT-related proteins, β-catenin, c-myc and p-GSK-3β inhibited by overexpressed circ_0009092 was significantly increased after transfected with miR-665 mimic (Fig. 5J, M). NLK knockdown reversed the miR-665 inhibitor-mediated Wnt/β-catenin signaling pathway (Fig. 5K-L). Overall, these data indicated that circ_0009092 could impair the suppression of miR-665 to NLK, which further regulated the Wnt/β-catenin signaling pathway. In CRC, the down expression of circ_0009092 led to the inhibition of NLK, then the Wnt/β-catenin signaling pathway was suppressed.
6. Circ_0009092/miR-665/NLK inhibited TAMs recruitment by suppressing CCL2 secretion in CRC.
Since the infiltration of immune cells in the TME is highly relevant to the CRC progression, we analyzed the distribution of different immune cells in CRC tissues and normal tissues though the Bioinformatics algorithm. We integrated six latest algorithms, including TIMER, xCell, MCP-counter, CIBERSORT, EPIC and quanTIseq. Results showed that The CRC tumors displayed massive macrophage infiltration, especially M2 macrophage infiltration (Fig. S5A-B). In addition, circ_0009092 was found to be inversely correlated with the density of TAMs marker CD68 in CRC (Fig. 6A). These observations led us to speculate that circ_0009092 was correlated with TAMs infiltration. An in vitro model of TAMs was then applied. THP-1 was treated with PMA and co-cultured with tumor cells for 48h. We observed that THP-1 macrophages were recruited by circ_0009092 knockdown, and overexpression of circ_0009092 inhibited the recruitment of macrophages (Fig. 6B). In addition, miR-665 mimic markedly promoted the recruitment of macrophages, and the miR-665 inhibitor could restrain the recruitment of macrophages induced by circ_0009092 knockdown (Fig. 6C-D). In addition, overexpression of NLK (oe-NLK) failed to recruit macrophages and knock down of NLK (si-NLK) showed the opposite effect. MiR-665 inhibitor substantially inhibited si-NLK-induced the recruitment of macrophages. While miR-665 mimic recovered the effect of oe-NLK (Fig. S6A-B). Thus, our results substantiated that circ_0009092/miR-665/NLK inhibits TAMs recruitment in the TME of CRC.
We then analyzed the potential chemokines associated with TAMs recruitment. RT-PCR results showed that mRNA expression of CCL2 emerged as the most prominently chemokine in circ_0009092 transfected tumor cells (Fig. 6E). ELISA analysis showed that CCL2 protein levels were significantly increased in the media from miR-665 mimic, sh- circ_0009092, or si-NLK group compared to NC group (Fig. 6F-G). Rescue experiments suggested that miR-665 mimic-induced CCL2 could be reversed by circ_0009092 overexpression (Fig. 6F, Fig.S6C). In addition, CCL2 neutralizing antibody suppressed the circ_0009092 knockdown-enhanced the recruitment of macrophages (Fig. 6H). Collectively, these data suggest that circ_0009092/miR-665/NLK inhibits the recruitment of TAMs via CCL2.
7. NLK bound to STAT3 protein and regulated its activity to suppress CCL2 expression.
To explore the molecular mechanism by which circ_0009092/miR-665/NLK regulated CCL2, we searched for ERK1/2, STAT3 and p53 signaling pathways, of which all had been reported to be the major molecule involved in the regulation of CCL2. WB results showed that p-STAT3 significantly increased after circ_0009092/miR-665/NLK treatment, whereas there was no significant change in the phosphorylation level of ERK1/2 (Fig. 7A). Rescue experiments showed that miR-665 inhibitor could not significantly change the expression of p-STAT3 after circ_0009092 knockdown (Fig. 7B). Furthermore, miR-665 inhibitor did not alter si-NLK-induced p-STAT3 expression (Fig. 7C). These results suggested that the circ_0009092/miR-665/NLK signal can regulate the p-STAT3 expression. It has been reported that NLK could bind to a wide range of transcription factors and regulated the transcriptional activity. We assessed whether NLK was capable of binding to STAT3, thus affecting CCL2 expression. We carried out a co-immunoprecipitation (CoIP) assay in tumor cells using the antibody against STAT3 or NLK. The CoIP results demonstrated the binding of NLK with STAT3 (Fig. 7D). In addition, IF results showed the co-localization of NLK and STAT3 in HCT-116 and Lovo cells (Fig. S7A). Besides of phosphorylation, STAT3 protein can also be modified by O-GlcNAcylation. O-GlcNAcylation is a key posttranslational modification, a single monosaccharide, β-O-GlcNAc is added onto serine or threonine residues of proteins. The O-GlcNAcylation occurs in an analogous fashion to protein phosphorylation, and there is an extensive crosstalk between O-GlcNAcylation and phosphorylation[19]. Previous studies showed that O-GlcNAcylation of STAT3 could promote STAT3 phosphorylation[20], so we wondered whether circ_0009092/miR-665/NLK could affect O-GlcNAcylation of STAT3. CoIP assay was performed and revealed that knockdown of circ_0009092 dramatically increased the O-GlcNAcylation modification level of STAT3, and increased binding between O-linked N-acetylglucosamine transferase (OGT) and STAT3. Meanwhile, overexpression of circ_0009092 reduced STAT3 O-GlcNAcylation modification and conjunction with OGT (Fig. 7E). In addition, miR-665 inhibitor inhibited the STAT3 O-GlcNAcylation modification and conjunction with OGT (Fig. S7B), while miR-665 mimic increased the STAT3 O-GlcNAcylation modification and conjunction with OGT (Fig. S7C). These results demonstrated that circ_0009092/miR-665/NLK regulate STAT3 by binding with STAT3 and regulating its phosphorylation and O-GlcNAcylation. Then we determined STAT3 transcriptional activity by a promoter reporter assay in which luciferase reporter was driven by a CCL2 promoter. The results showed that STAT3 induced the promoter activity of CCL2 in tumor cells transfected with the promoter region of CCL2. Luciferase activity of CCL2 promoter was significantly decreased in transfected with knockdown of STAT3 (Fig. 7F). Meanwhile, si-NLK significantly increased the STAT3-mediated transcription activity, and oe-NLK inhibited STAT3-mediated transcription activity of CCL2 promoter. MiR-665 mimic could increase the luciferase activity of CCL2 promoter that was restrained by oe-NLK, and miR-665 inhibitor could not reverse the luciferase activity of CCL2 promoter that was increased by si-NLK (Fig. 7G). Moreover, overexpression of circ_0009092 reversed the luciferase activity of CCL2 promoter that was increased by miR-665 mimic (Fig. S7D), and knockdown of circ_0009092 increased the luciferase activity of CCL2 promoter that was reduced by miR-665 inhibitor (Fig. S7E). Taken together, these results demonstrate that NLK bound to STAT3 and increases its transcriptional activity of CCL2.
8. EIF4A3 suppressed circ_0009092 expression in CRC cells.
RNA-binding proteins (RBPs), such as QKI, ESRP1, and EIFEA3, can drive cyclization of RNA to regulate the biogenesis of circRNAs. Therefore, CircInteractome was used to predict the RNA binding proteins that match the flanking regions of the circ_0009092, of which EIF4A3 has the most binding sites (Fig. 7H). PCR results showed that the expression of circ_0009092 was significantly decreased after the overexpression of EIF4A3. Knockdown of EIF4A3 markedly increased the expression of circ_0009092 in CRC cells (Fig. 7I). GEPIA was used to examine the expression of EIF4A3 in CRC samples, and results showed that EIF4A3 expression was increased in CRC tissues than in the ANT (Fig. 7J). In addition, qPCR results also revealed that the expression of EIF4A3 was elevated in CRC tissues than in the ANT (Fig. 7K). And the expression level of EIF4A3 mRNA was negatively correlated with circ_0009092 expression level, suggesting that EIF4A3 is related to circ_0009092 expression (Fig. 7L). In addition, the RIP assay using anti-EIF4A3 antibody showed that EIF4A3 bound with circ_0009092 through two upstream binding sites a and b but not circ_0009092 (region c) in the RNA-protein complex (Fig. 7M). Taken together, these data demonstrated that circ_0009092 is regulated by EIF4A3 and EIF4A3 suppresses circ_0009092 expression by binding to flanking sequences.
9. Circ_0009092 hampered the progression of colorectal cancer in vivo.
To evaluate the function of circ_0009092 in vivo, HCT-116 cells with stably expressing either negative control (oe-NC) or oe-circ_0009092 were subcutaneously injected into the BALB/c nude mice, and THP-1 cells were injected into the caudal vein of mice. It was found that the circ_0009092 overexpression group showed significantly lower tumor weights and volumes compared to the NC group (Fig. 8A-C). Circ_0009092 overexpression group showed a lower expression of miR-665 and higher expression of NLK compared with the oe-NC group (Fig. 8D-F). WB results showed that circ_0009092 affects the expression of NLK, CCL2, p-STAT3, OGT, and Wnt/β-catenin singaling pathway related protein (Fig. 8G). FISH analysis showed the same trend as the PCR and WB results (Fig. 8H). Additionally, decreased CCL2 level and CD68 (the macrophage marker) expression level were detected in immunohistochemistry (IHC) staining in the circ_0009092 overexpression group (Fig. 8I), which indicated that the circ_0009092 overexpression group presented less macrophage infiltration compared to oe-NC group. Tumors of oe-circ_0009092 group exhibited markedly reduced recruitment of macrophages as compared with oe-NC group (Fig. 8J). Furthermore, the number of M2 macrophages (CD163, CD206) was decreased in the oe-circ_0009092 group compared with oe-NC group, while the number of M1 macrophages (CD86) showed no significant difference between two groups (Fig. 8J). Therefore, these results corroborated that circ_0009092/miR-665/NLK inhibits CRC growth and macrophage recruitment in vivo. Taken together, our data demonstrated that circ_0009092/miR-665/NLK suppresses tumor progression and TAMs recruitment by inhibiting the Wnt/β-catenin signaling pathway in the CRC (Fig. 9).