Research Paper Volume 11, Issue 24 pp 12412—12427

CircRNA_100367 regulated the radiation sensitivity of esophageal squamous cell carcinomas through miR-217/Wnt3 pathway

Junqi Liu1, , Nannan Xue1, , Yuexin Guo1, , Kerun Niu2, , Liang Gao3, , Song Zhang1, , Hao Gu1, , Xin Wang1, , Di Zhao4, , Ruitai Fan1, ,

  • 1 Department of Radiation Oncology, The First Affiliated Hospital of Zhengzhou University, Erqi, Zhengzhou 450000, China
  • 2 Department of Molecular and Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg 69120, Germany
  • 3 Center of Experimental Orthopaedics, Saarland University Medical Center, Kirrberger Strasse, Homburg 66421, Germany
  • 4 Endocrinology Department, The First Affiliated Hospital of Zhengzhou University, Zhengzhou 450000, China

Received: July 6, 2019       Accepted: November 26, 2019       Published: December 18, 2019      

https://doi.org/10.18632/aging.102580
How to Cite
This article has been corrected. See Correction. Aging (Albany NY). 2021; 13:23868-23870 . https://doi.org/10.18632/aging.203664  PMID: 34738919

Copyright © 2019 Liu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Background: Circular RNAs (circRNAs) play important roles in regulating the radioresistance of esophageal squamous cell carcinoma (ESCC). This study aimed to determine the role of hsa_circRNA_100367 in regulating radioresistance of ESCC.

Results: Higher expression and potency of endothelial to mesenchymal transformation (EMT) was found in radioresistant ESCC cells (KYSE-150R) than in ESCC cells (KYSE-150). Silencing circRNA_100367 inhibited the proliferation and migration of KYSE-150R cells, and decreased the expression of β-catenin (an important molecule in Wnt pathway) in KYSE-150R cells. Additionally, circRNA_100367 bound to miR-217, and miR-217 targeted Wnt3. Low Wnt3 expression was associated with the short survival time in patients with ESCC and Wnt3 knockdown inhibited the proliferation and migration of KYSE-150R cells. CircRNA_100367 enhanced the radioresistance of KYSE-150R cells through miR-217/Wnt3 pathway. In vivo, circRNA_100367 silence reduced the growth of KYSE-150R cells under radiation.

Conclusion: Our results revealed that circRNA_100367 attenuated radioresistance of ESCC through miR-217/Wnt3 pathway.

Methods: CircRNAs related with the radioresistance of ESCC were analyzed by hierarchical cluster analysis. The relationship between circRNA_100367 and miR-217, Wnt3 was detected by RNA immunoprecipitation (RIP), RNA pull-down and luciferase reporte assays. The proliferation and migration ESCC cells were detected by MTT, Transwell and colony formation assays.

Introduction

Esophageal squamous cell carcinoma (ESCC), as the predominant form of esophageal cancer and accounting for about 88% of esophageal cancer, is the sixth leading cause of cancer death in the world and the second most deadly cancer in China [1, 2]. Preoperative radiotherapy is a vital treatment of ESCC in clinic [3, 4], but it is not effective for all patients with ESCC because of the occurrence of radioresistance of ESCC cells [5]. Studies have shown that the occurrence of radioresistance of ESCC cells is related with the abnormal expression of genes [6, 7]. Therefore, it is important to find abnormally expressed genes that related with the radiation sensitivity of ESCC, which will provide potential therapeutic targets for the treatment of ESCC.

Circular RNAs (circRNAs) are universal diverse endogenous noncoding RNAs that have covalently linked ends which can be acted as microRNA (miRNA) sponges to regulate gene expressions in many cancers [8]. Studies have reported that circRNAs are abnormally expressed in various cancer cells, such as cervical cancer cell, hepatic stellate cells and ESCC cells, which play important roles in regulating the radiation sensitivity of these cancers [911]. However, there are few studies focused on a specific circRNA in regulating the radiation sensitivity of ESCC. Recently, the expressions of 3752 circRNAs are detected in ESCC cell line KYSE-150 and radioresistant ESCC cell KYSE-150R, and hsa_circRNA_100367 (gene Symbol DDB1 and CUL4 associated factor 8 (DCAF8), also named hsa_circRNA_0014879) is found to be highly expressed in KYSE-150R cell line [12]. Therefore, circRNA_100367 may play a vital role in regulating the radiation sensitivity of ESCC. However, the underlying mechanism of circRNA_100367 in the radiation sensitivity of ESCC is still not clear.

Studies have reported that circRNA exert its function in ESCC through competitive binding with miRNAs. For example, circRNA VRK1 inhibited the progression of ESCC through regulating miR-624-3p/ phosphatase and tensin homolog (PTEN) signaling pathway [11]. So, finding the downstream miRNAs/proteins that regulated by circRNA_100367 is important to reveal the underlying mechanism of circRNA_100367 mediated radiation sensitivity of ESCC. miR-217 is a widely expressed miRNA and acts as a tumor suppressor gene DCAF or oncogene in various cancers [1315]. In addition, it has been reported that miR-217 can inhibit the progression of ESCC [16], and enhance the radiation sensitivity to doxorubicin in acute myeloid leukemia [17]. Base on the analysis of bioinformatics software (circinteractome), we found there are binding sties between circRNA_100367 and miR-217. Therefore, circRNA_100367 might regulate the radiation sensitivity of ESCC through sponging miR-217.

Previous studies have shown that the 3’-untranslated region (UTR) of target mRNAs can be bound by specific miRNAs to play important roles in the regulation of radiation sensitivity of ESCC [18, 19]. Wnt3 is a member of Wnt family, and is the key molecule of Wnt-β-catenin signaling pathway. The expression of Wnt3 can activate Wnt-β-catenin signaling pathway in liver cancer, gastric cancer, lung cancer, and ESCC [2022], etc. Studies have shown that Wnt3 is up-regulated in cancer tissues, such as gastric cancer and lung cancer, and knockdown of Wnt3 can suppress the progression of cancers [20, 21]. Wnt3 is also reported to be associated with the tumorigenicity and poor prognosis of ESCC [22, 23]. Importantly, Wnt-β-catenin signaling is found to enhance radioresistance of ESCC [24]. Bioinformatics softwares predicted there are binding sites between miR-217 and Wnt3 3’UTR. However, whether circRNA_100367 can regulate the radiation sensitivity of ESCC through miR-217/Wnt3 pathway is still known. In the present study, we aimed to elucidate this question.

Results

hsa_circRNA_100367 was highly expressed in human radioresistant esophageal cancer cell lines

According to the online circRNA information in radioresistant esophageal cancer cells 12 [12], 74 differentially expressed circRNAs (FDR<0.01, fold change>2 or <0.5) were selected for circRNAs chips (Figure 1A). Then, hierarchical cluster analysis was performed to select the differentially expressed circRNAs (fold change>5 or < 0.2) between the KYSE-150 and KYSE-150R cells, and the expressions of those circRNAs were measured in different esophageal cancer cell lines (TE-1/TE-1R, KYSE-150/KYSE-150R, ECA-109/ECA-109R cells). As shown in Figure 1B, the expression trend of circRNA_100367 was consistent in TE-1/TE-1R, KYSE-150/KYSE-150R, and ECA-109/ECA-109R cells. Moreover, circRNA_100367 showed an increased expression trend with a most extent in KYSE-150/KYSE-150R cells. Therefore, circRNA_100367 was selected for further experiments.

Expression of hsa

Figure 1. Expression of hsa_circRNA_100367 in human radioresistant esophageal cancer cell lines. (A) Hierarchical cluster analysis showed distinguishable circRNA expression profiles between KYSE-150 and KYSE-150R cells. (B) qRT-PCR detected circRNA_100385, circRNA_104983, circRNA_001059, circRNA_100984, circRNA_100367, circRNA_103783 and circRNA_100312 expressions in TE-1, TE-1R, KYSE-150, KYSE-150R, ECA-109, ECA-109R cells. *p<0.05, **p<0.01 vs. TE-1; #p<0.05, ##p<0.01 vs. KYSE-150; aap<0.01 vs ECA-109.

CircRNA_100367 exists in KYSE-150R cells with higher potency of endothelial to mesenchymal transformation (EMT)

As shown in Figure 2A, EMT-related molecule E-cadherin was decreased and vimentin and snail expressions were increased in KYSE-150R cells than that of KYSE-150 cells, indicating stronger molecular characteristics of EMT in KYSE-150R cells. Besides, the number of migrated KYSE-150R cells was greater than KYSE-150 cells (Figure 2B). The number of CD133+ cells in KYSE-150R cells was greater than that in KYSE-150 cells (Figure 2C). The role of circRNA_100367 was further determined in KYSE-150/KYSE-150R cells. The circular properties of cricRNA 100367 was identified with divergent primers and convergent primers (Figure 2D). Compared with mock KYSE-150R cells, circRNA_100367 expression was significantly down-regulated in RNase R-treated KYSE-150R cells. DCAF8 mRNA expression was significantly down-regulated in RNase R-treated KYSE-150 cells and KYSE-150R cells (Figure 2E). However, the degradation degree of circRNA_100367 was much lower than DCAF8 mRNA, indicating circRNA_100367 exists in esophageal squamous cell carcinomas.

CircRNA

Figure 2. CircRNA_100367 exists in KYSE-150R cells with higher potency of endothelial to mesenchymal transformation (EMT). (A) The protein levels characteristic molecules (β-catenin, Vimentin, and Snail) of EMT in KYSE-150 and KYSE-150R cells were detected by western blot. (B) The migration of KYSE-150 and KYSE-150R cells were determined by transwell assay. (C) The CD133 positive cells of KYSE-150 and KYSE-150R cells were measured by flow cytometry. (D) Divergent primers and convergent primers were designed, and PCR product of cricRNA 100367 was detected in 1.5% agarose gel electrophoresis. (E) After the treatment of RNase R, circRNA_100367 and DCAF8 mRNA expressions in KYSE-150 and KYSE-150R cells were detected by qRT-PCR. *p<0.05, p<0.01 vs. KYSE-150 or mock.

Silencing circRNA_100367 inhibited the proliferation and migration of KYSE-150R cells

Since circRNA_100367 was highly expressed in KYSE-150R cells, sh- circRNA_100367 was transfected into KYSE-150R cells to block the circRNA_100367 expression. As shown in Figure 3A, the viability of KYSE-150 and KYSE-150R cells was decreased with the increase of radiation dose. The viability of KYSE-150R cells was significantly decreased after the treatment of sh- circRNA_100367, indicating silencing circRNA_100367 enhanced the radiation sensitivity of KYSE-150R cells. Silencing circRNA_100367 reduced the number of clones of KYSE-150R cells under the radiation dose of 0Gy and 6Gy, and decreased the survival fraction of KYSE-150R cells under the radiation dose of 4Gy, 6Gy and 8Gy (Figure 3B). In addition, the migration ability of KYSE-150R cells was stronger than KYSE-150 cells, and silencing circRNA_100367 decreased the migration of KYSE-150R cells (Figure 3C). Moreover, we found circRNA_100367 silence reduced the expression of β-catenin, Vimentin and Snail in KYSE-150R cells (Figure 3D).

Effects of circRNA

Figure 3. Effects of circRNA_100367 on the proliferation and migration of KYSE-150R cells. (A) KYSE-150R cells were transfected with Sh-circRNA_100367 (or NC) and then irradiated with 0, 3, and 6 Gy X-ray. Cell ability was detected by MTT assay. (B) KYSE-150R cells were transfected with Sh-circRNA_100367 (or NC) and then irradiated with different dose X-ray. The proliferation of cells was evaluated by colony formation assay and survival fraction was analyzed by multi-target single-hit model. (C) KYSE-150R cells were transfected with Sh-circRNA_100367 (or NC). The migration of cells was assessed by wound healing assay. (D) KYSE-150R cells were transfected with Sh-circRNA_100367 (or NC). The protein levels of β-catenin, Vimentin, and Snail were detected by western blot. NC: the negative control of Sh-circRNA_100367. **p<0.01 vs. KYSE-150; #p<0.05, ##p<0.01 vs. KYSE-150+NC.

The interaction between circRNA_100367 and miR-217

Given that circRNAs can be acted as miRNA sponge and circRNA_100367 is abundant in the cytoplasm, we wonder if circRNA_100367 can bind to miR-217 to affect the anti-radioresistant phenotypes. We conducted immunoprecipitation and dual luciferase reporter gene assay to investigate whether circRNA_100367 could bind to miR-217. Anti-AgO2 immunoprecipitation showed circRNA_100367 was enriched in KYSE-150 and KYSE-150R cells treated with AgO2 antibody, suggesting the occupancy of AgO2 in the region of the circRNA_100367 (Figure 4A). Besides, we found most of circRNA_100367 are located in the cytoplasm of KYSE-150 and KYSE-150R cells (Figure 4B). Also, the binding sites between hsa_circRNA_100367 and miR-217 were analyzed by bioinformatics software (circinteractome) (Figure 4C). The luciferase activity of miR-217 was significantly increased in circRNA_100367 mut (Figure 4D), indicating miR-217 could bind to the binding sites of circRNA_100367. Moreover, the enrichment of circRNA_100367 in biotin-labeled miR-217 further proved the interaction between circRNA_100367 and miR-217 (Figure 4E).

The interaction between circRNA

Figure 4. The interaction between circRNA_100367 and miR-217. (A) KYSE-150 and KYSE-150R cells were transfected with Flag-tagged AgO2 (Flag-AgO2) or negative control Flag-tagged Green Fluorescent Protein (Flag-GFP). Then cell lysates were incubated with Flag or IgG (negative control of Flag) antibodies. The enrichment of circRNA_100367 in the immunoprecipitation products was detected by qPCR. **p<0.01 vs. Anti-IgG. (B) The expression of circRNA_100367 in the nucleus and cytoplasm of KYSE-150 and KYSE-150R cells was detected by qRT-PCR. (C) The diagram showed the binding sites between hsa_circRNA_100367 and miR-217. (D) KYSE-150R cells were co-transfected with miR-217 mimic and luciferase reporter plasmids containing circRNA_100367 or circRNA_100367 mut. The relative luciferase activity was detected by dual luciferase reporter gene assay. **p<0.01 vs. circRNA_100367 mut. (E) The interaction between miR-217 and circRNA_100367 was assessed by RNA pull-down assay. **p<0.01 vs. NC.

Wnt3 regulated the proliferation and migration of KYSE-150 cells

According to the bioinformatics software microT-CDS, miRDB and TargetScan, Wnt was one of the potential target genes of miR-217 with higher scores (Figure 5A). The other target genes of miR-217 were shown in the Supplementary Table 2. Wnt3 mRNA level was decreased in KYSE-150R cells after the treatment of miR-217 mimic under the radiation dose of 0Gy and 6Gy (Figure 5B). According to the data from The Cancer Genome Atlas (TCGA), there is an association between lower expression of Wnt3 and survival in patients with ESCC (n=184) (Figure 5C). As shown in Figure 5D, overexpression of Wnt3 promoted the migration of KYSE-150 cells, while knockdown of Wnt3 inhibited the migration of KYSE-150R cells. In addition, sh-Wnt3 reduced the number of migrated KYSE-150R cells (Figure 5D). Wnt3 also affected the expression of nucleus β-catenin. Wnt3 overexpression promoted nucleus β-catenin expression in KYSE-150 cells, and sh-Wnt3 down-regulated nucleus β-catenin expression in KYSE-150R cells (Figure 5D). Wnt3 overexpression increased the number of clones of KYSE-150 cells under the radiation dose of 0Gy and 6Gy, suggesting Wnt3 overexpression enhanced the radioresistant of KYSE-150 cells (Figure 5E). And Sh-Wnt3 inhibited the number of clones of KYSE-150R cells under the radiation dose of 0Gy and 6Gy, suggesting Sh-Wnt3 inhibited the radioresistant of KYSE-150R cells (Figure 5E).

Effects of Wnt3 on the proliferation and migration of KYSE-150 and KYSE-150R cells. (A) The potential target genes of miR-217 were predicted by bioinformatics softwares (microT-CDS, miRDB and TargetScan). (B) KYSE-150R cells were transfected with miR-217 mimic (or NC) and then irradiated with 0 and 6 Gy X-ray. The mRNA level of Wnt3 was detected by qRT-PCR. **pC) The relationship between Wnt3 expression and survival rate of patients with ESCS. (D) KYSE-150 cells were transfected with Wnt3 overexpression vector and KYSE-150R cells were transfected with Wnt3 silence vector. The invasion of KYSE-150 and KYSE-150R cells were detected by transwell assay, and the protein level of nuclear β-catenin (n β-catenin) was determined by western blot. **pE) KYSE-150 cells were transfected with Wnt3 overexpression vector and KYSE-150R cells were transfected with Wnt3 silence vector. Then cells were irradiated with 0 and 6 Gy X-ray. The proliferation of cells was evaluated by colony formation assay.

Figure 5. Effects of Wnt3 on the proliferation and migration of KYSE-150 and KYSE-150R cells. (A) The potential target genes of miR-217 were predicted by bioinformatics softwares (microT-CDS, miRDB and TargetScan). (B) KYSE-150R cells were transfected with miR-217 mimic (or NC) and then irradiated with 0 and 6 Gy X-ray. The mRNA level of Wnt3 was detected by qRT-PCR. **p<0.01 vs NC. (C) The relationship between Wnt3 expression and survival rate of patients with ESCS. (D) KYSE-150 cells were transfected with Wnt3 overexpression vector and KYSE-150R cells were transfected with Wnt3 silence vector. The invasion of KYSE-150 and KYSE-150R cells were detected by transwell assay, and the protein level of nuclear β-catenin (n β-catenin) was determined by western blot. **p<0.01 vs. KYSE-150; #p<0.05 vs. KYSE-150R. (E) KYSE-150 cells were transfected with Wnt3 overexpression vector and KYSE-150R cells were transfected with Wnt3 silence vector. Then cells were irradiated with 0 and 6 Gy X-ray. The proliferation of cells was evaluated by colony formation assay.

CircRNA_100367 coordinated with miRNA-217 regulated the radioresistant of KYSE-150R cells through Wnt3

As shown in Figure 6A, miR-217 overexpression decreased the number of clones of KYSE-150R cells, whereas miR-217 mimic coordinated with circRNA_100367 increased the number of clones under the radiation dose of 0Gy and 6Gy. In addition, miR-217 overexpression reduced the number of KYSE-150R cells, whereas miR-217 mimic coordinated with circRNA_100367 increased the number of KYSE-150R cells (Figure 6B). miR-217 overexpression down-regulated the expressions of Wnt3 and nucleus β-catenin, whereas miR-217 mimic coordinated with circRNA_100367 up-regulated the expressions of Wnt3 and nucleus β-catenin (Figure 6C). Moreover, miRNA-217 overexpression decreased the luciferase activity of Wnt3 3’-UTR, whereas miRNA-217 mimic coordinated with circRNA_100367 increased the luciferase activity of Wnt3, indicating miRNA-217 and circRNA_100367 regulated cooperatively the 3′UTR of Wnt3 (Figure 6D).

CircRNA

Figure 6. CircRNA_100367 coordinated with miRNA-217 regulated the radioresistant of KYSE-150R cells through Wnt3. (A) KYSE-150R cells were transfected with miR-217 mimic (or NC) or co-transfected with miR-217 mimic and circRNA_100367. Then cells were irradiated with 0 and 6 Gy X-ray. irradiated with 0 and 6 Gy X-ray. (BC) KYSE-150R cells were transfected with miR-217 mimic (or NC) or co-transfected with miR-217 mimic and circRNA_100367. (B) The invasion of KYSE-150 and KYSE-150R cells were detected by transwell assay. (C) The protein levels of Wnt3 and nuclear β-catenin (n β-catenin) were detected by western blot. (D) KYSE-150R cells were co-transfected with miR-217 mimic, circRNA_100367 and luciferase reporter plasmids containing Wnt3 WT or Wnt3 mut. The relative luciferase activity was detected by dual luciferase reporter gene assay.

Silencing circRNA_100367 reduced the growth of KYSE-150R cells under radiation

Since silencing circRNA_100367 inhibited the proliferation and migration of KYSE-150R cells in vitro, we determined the role of circRNA_100367 in vivo using a xenograft nude mouse model. KYSE-150 cells transfected with circRNA_100367 stably expressed high circRNA_100367 level, and KYSE-150R cells transfected with sh-circRNA_100367 stably expressed low circRNA_100367 level. We found that circRNA_100367 overexpression accelerated the growth of KYSE-150 cells, and silencing circRNA_100367 reduced the growth of KYSE-150R cells under radiation (Figure 7A). Tumor volume was significantly larger in KYSE-150R+Gy group than that of KYSE-150+Gy group. CircRNA_100367 overexpression significantly increased the tumor volume of KYSE-150+circRNA_100367+Gy group, and silencing circRNA_100367 significantly reduced the tumor volume of KYSE-150R+sh-circRNA_100367+Gy group (Figure 7B). The protein level of E-cadherin was decreased and the protein levels of vimentin, snail, Wnt3, andβ-catenin were increased in the KYSE-150+sh-circRNA_100367+Gy group compared with KYSE-150+ Gy group. Also, the protein level of E-cadherin was elevated and the protein levels of vimentin, snail, Wnt3, andβ-catenin were reduced in the KYSE-150R+sh-circRNA_100367+Gy group compared with KYSE-150R+ Gy group (Figure 7C).

Effect of circRNA

Figure 7. Effect of circRNA_100367 on tumor growth of KYSE-150R cells under radiation. KYSE-150R cells were stably transfected with Sh-circRNA_100367 or negative control (Circ) and then were subcutaneously inoculated into nude mice. 10 days after inoculation, mice were irradiated with 6 Gy X-ray. (A, B) Representative pictures and volumes of excised tumors. (C) The protein levels of E-cadherin, vimentin, snail, Wnt3 and β-catenin in excised tumors were measured by western blot. (D) A schematic diagram representing the role and mechanism of circRNA_100367 in radiation sensitivity of ESCC.

Discussion

Increasing evidences have revealed that the abnormal expressions of circRNAs are related to the radiation sensitivity of cancers [9, 10]. However, few studies focus on the abnormally expressed circRNAs in regulating radiation sensitivity of ESCC. In this study, the upregulation of circRNA_100367 was observed in KYSE-150/KYSE-150R cells with a most extent than the other two ESCC cell lines and their radioresistant cells. Also, previous studies showed abnormally expressed circRNAs are related with the phenotypic change of cancer cells [25, 26]. So we further transfected sh-circRNA_100367 into KYSE-150R cells to determine whether circRNA_100367 changed the phenotype of ESCC radioresistant cells. Results showed that silencing circRNA_100367 decreased the viability and survival fraction of KYSE-150R cells, reduced the number of clones of KYSE-150R cells, and inhibited the migration of KYSE-150R cells under radiation. These results indicated that upregulation of circRNA_100367 suppressed the radiation sensitivity of radioresistant ESCC KYSE-150R.

Previous researches have reported that disease-specific miRNAs can be sponged by circRNAs in many cancers [26, 27]. miR-217 is one of these disease-specific miRNAs and exerts its functional role in a variety of cancers [28]. For example, abnormally expressed miR-217 enhanced the chemosensitivity of acute myeloid leukemia and cervical carcinoma [17, 29]. However, whether abnormally expressed miR-217 involved in regulating the radiation sensitivity of ESCC is not clear. Based on the mechanism of circRNAs sponging miRNAs in cancers and miR-217 predicted as a target of circRNA_100367, we conducted RIP and luciferase reporter gene assay, and proved the interaction between circRNA_100367 and miR-217. To investigate the in-depth underlying mechanism of circRNA_100367/miR-217, miR-217 mimic or miR-217 mimic+circRNA_100367 was transfected into KYSE-150R cells to determine their effect on colony formation and migration of KYSE-150R cells. Results showed miR-217 mimic coordinated with circRNA_100367 promoted colony formation and migration of KYSE-150R cells under the radiation dose, which indicated miR-217 mimic+circRNA_100367 attenuated radiation sensitivity of KYSE-150R cells. So far, no other studies have demonstrated the role of circRNA_100367/miR-217 in the regulation of radiation sensitivity of KYSE-150R cells, which will provide directions for the improving the survival rate of ESCC.

Wnt3, as a member of Wnt family, has been proved to promote the stabilization of β-catenin to regulate the radiation sensitivity of cancer cells [30, 31]. In this study, we found silencing Wnt3 down-regulated the expression of nucleus β-catenin, which was consistent with previous report [31]. However, the exact role of Wnt3 in the regulation of radiation sensitivity of ESCC cells is still unclear, although Wnt-β-catenin signaling has been reported as an important pathway in regulating the radioresistance of ESCC [32]. In this study, results showed that silencing Wnt3 enhanced the radiation sensitivity of KYSE-150R cells. In addition, silencing Wnt3 changed the phenotype of KYSE-150R cells, which suppressed the colony formation and the migration of KYSE-150R cells. These findings suggested that Wnt3 suppressed the radiation sensitivity of ESCC cells. Moreover, studies have demonstrated that Wnt3 can be the target of disease-specific miRNAs [33, 34]. In this study, we investigated the relationship between miR-217 and Wnt3 in KYSE-150R cells, and proved Wnt3 was a target of miR-217 according to the bioinformatics analysis and luciferase reporter gene assay. And Wnt3 expression was found to be negatively correlated with miR-217. Furthermore, circRNA_100367 reversed the effect of miR-217 on Wnt3 expression, the number of clones of KYSE-150R cells and the migration of KYSE-150R cells. These results indicated that miR-217/Wnt3 was the downstream targets of circRNA_100367 in regulating the radiation sensitivity of ESCC.

In conclusion, highly expressed circRNA_100367 is related with the radiation sensitivity of ESCC, and silencing circRNA_100367 decreased the proliferation and migration of KYSE-150R cells in vitro, and inhibited tumor growth in vivo. In addition, miR-217/Wnt3 was demonstrated as the downstream targets of circRNA_100367 in regulating the radiation sensitivity of ESCC (Figure 7D), which will provide potential targets for reducing radiotherapy failure of patients with ESCC.

Materials and Methods

Cell culture and transfection

Esophageal cancer cell lines TE-1, KYSE-150 and ECA-109 were purchased from Cell Bank of Chinese Academy of Sciences, and cultured in RPMI-1640 (Gibco, Grand Island, NY, USA) supplemented with 10% fetal bovine serum (FBS; Gibco), 100 unit/ml of penicillin (Gibco) and 100 mg/ml of streptomycin (Gibco) at 37°C in an incubator containing 5% CO2 and 95% air.

Short hairpin RNA that targets human circRNA_100367 (sh-circRNA_100367), pcDNA3.0 vectors overexpressing circRNA_100367 (circRNA_100367), sh-Wnt3, and negative controls were synthesized by Genechem (Shanghai, China). miR-217 mimic (Cat. No. miR10000274-1-5) and its negative control were synthesized by RiboBio (Guangzhou, China). 100 nmol vectors were transfected into cells using Lipofectamine 2000 (Invitrogen, UAS). The sequences of those vector were as follows: miR-217 mimic 5′UACUGCAUCAGGAACUGAUUGGA3′; Sh-Wnt3 5′GCGCTTCTGCCGCAATTACAT3′; Sh-CircRNA-100367 5′GCGTCTCCTTCAGTGAATCTA3′; Sh-NC 5′TTCTCCGAACGTGTCACGTCT3′.

Development of radioresistant esophageal cancer cells

Radioresistant esophageal cancer cell lines TE-1R, KYSE-150R and ECA-109R were established by multiple fractionated irradiation according to previous report [35]. TE-1, KYSE-150 and ECA-109 cells (1.5×106) were seeded into 25 cm2 culture flasks, and irradiated with 1 Gy X-ray. Culture medium was immediately replaced after irradiation. After cells reached 90% confluence, cells were transferred to new culture flasks. Then, cells received second irradiation when they reached 50% confluence. In total, cells received 1Gy for three times, 2Gy for three times, 4Gy for three times with a total dose of 21Gy to establish radioresistant cells.

Quantitative real-time PCR (qRT-PCR)

Total RNA was extracted from TE-1, KYSE-150 and ECA-109 cells using Trizol Reagent (Invitrogen, Carlsbad, CA, USA). The nucleus and cytoplasm were extracted by NE-PER™ Nuclear and Cytoplasmic Extraction Reagents (Thermo Scientific, Waltham, MA, USA). 1μg RNA was used to synthesize cDNA using PrimeScript II 1st Strand cDNA Synthesis Kit (Takara Biomedical Technology, Beijing, China). qRT-PCR was performed with One Step TB Green PrimeScript PLUS RT-PCR Kit (Takara) using 7500 Real-Time PCR System (Applied Biosystems, FosterCity, CA, USA). GAPDH was used as the internal reference for circRNAs and Wnt3. The 2-ΔΔCT [36] method was used to calculate the relative expressions of circRNA_100385, circRNA_104983, circRNA_001059, circRNA_100984, circRNA_100367, circRNA_103783, circRNA_100312 and Wnt3. The primer of U6 was purchased from Genechem (Cat. No. MQPS0000002-1-100; Shanghai, China). The sequences of other primers were shown in the Supplementary Table 1.

Western blotting

KYSE-150 and KYSE-150R cells were collected and lysed using RIPA Lysis and Extraction Buffer (Thermo Scientific), and protein concentration was measured by BCA Protein Assay Kit (Beyotime Biotechnology, Shanghai, China). The protein samples were separated by 10% SDS-PAGE and transferred to PVDF membranes (Millipore, Billerica, MA, USA). The membrane was incubated with 5% skim milk at room temperature for 2 h, then incubated with primary antibodies against E-cadherin (1:10000; ab40772, Abcam), vimentin (1:2000; ab92547, Abcam), snail (1:200; ab82846, Abcam), β-catenin (1:5000; ab227499, Abcam), Wnt3 (1:1000; ab32249, Abcam), β-actin (1:5000; ab227387, Abcam), histone H3 (1:1000; MA5-15150, Abcam) overnight at 4°C. The membrane was subsequently incubated with HRP-conjugated secondary antibodies (1:5000) at room temperature for 1.5 h. The proteins were visualized using enhanced chemiluminescence (Thermo Scientific), and the band intensity was determined by ChemiDoc MP Imaging System (Bio-Rad, Hercules, CA, USA).

Transwell assay

Transwell assay was conducted to detect the migration ability of KYSE-150 and KYSE-150R cells. Cells (2×104 cells) were seeded on the upper chamber of Transwell, and RPMI-1640 medium with 10% FBS was added to the bottom chamber. After 48 h, cells on the bottom chamber were fixed with 4% paraformaldehyde and stained with 0.1% crystal violet solution and counted.

Flow cytometry

KYSE-150 cells and KYSE-150R cells were stained with anti-CD133 antibody (Invitrogen) according to the manufacturer’s instruction. The percentage of CD133+ cells of KYSE-150 cells and KYSE-150R cells was determined by flow cytometry analysis using FACS Calibur (BD, San Jose, CA, USA).

Colony formation assay

The 1.2% and 0.7% low-point agarose solutions were prepared at first. The 1:1 mixture of 1.2% agrose solution and RPMI-1640 medium was added into 6-well plates (1.5mL/well). The 0.7% agarose solution was mixed with RPMI-1640 medium supplemented with FBS, and KYSE-150 cells and KYSE-150R cells were resuspended in this mixture. Then cells (1×103cells/ well) were seeded in 6-well plates and incubated for two weeks at 37°C. Then, colonies were fixed with 4% paraformaldehyde, stained with 0.1% crystal violet solution and counted.

RNase R treatment

Total RNA (2 μg) was incubated with 3 U/μg RNase R (Epicentre Technologies, Madison, WI, USA) for 15 min at 37°C. Then, the RNA expression levels of DCAF8 and circRNA_100367 were detected by qRT-PCR.

Cell proliferation assay

1-(4,5-Dimethylthiazol-2-yl)-3,5-diphenylformazan (MTT) assay kit (Sigma-Aldrich, Billerica, MA, USA) was used to detect the proliferation of KYSE-150 cells and KYSE-150R cells with different treatments under the radiation dose of 0Gy, 3Gy and 6Gy. Cells (2×103 cells/well) were seeded into a 24-well plate, then 10 μl MTT solution was added to the cells in each well and incubated for 4 h at 37°C. The supernatant was discarded, then 200μl dimethyl sulfoxide (DMSO, Sigma-Aldrich) were added to the cells. The absorbance at the 490 nm was analyzed using a microplate reader.

Wound healing assay

KYSE-150 cells and KYSE-150R cells with different treatments were seeded, and a linear wound was generated within the monolayers by scraping the cells using the sterile pipette tip. After 24 h of incubation, migration was imaged with an inverted microscope.

RNA immunoprecipitation

Magna RIP RNA-Binding Protein Immunoprecipitation Kit (Millipore) was used to perform RNA immunoprecipitation (RIP) assay. The lysate of KYSE-150 cells or KYSE-150R cells was prepared by 1.5×107 cells using 0.25 μl protease inhibitor, 0.125 μl RNase inhibitor and 50 μl RIP lysis buffer. After centrifugation, the supernatant was collected and incubated with RIP wash buffer, and the procedures were performed according to the manufacturer’s instruction. Finally, qRT-PCR was conducted to detect circRNA_100367 in RNA-binding protein complex absorbed by the magnetic beads.

Dual luciferase reporter gene assay

According to the bioinformatics software (circinteractome), we predicted the binding sties between circRNA_100367 and miR-217. According to the bioinformatics software (microT-CDS, miRDB and TargetScan), we predicted the potential target genes of miR-217. CircRNA_100367 wild-type (circRNA_100367) and mutant (circRNA_100367 mut) reporter vectors were constructed and inserted into the psiCHECK-2 (Promega, Madison, Wisconsin, USA). Wnt3 wild-type (Wnt3 WT) and Wnt3 mutant (Wnt3 mut) reporter vectors were also constructed and inserted into the psiCHECK-2 (Promega). 100 nmol reporter plasmids (circRNA_100367, circRNA_100367 mut) together with miR-217 or negative control were transfected into KYSE-150R cells. The reporter plasmids (Wnt3 WT, Wnt3 mut) together with miR-217 mimic, circRNA_100367+miR-217 mimic, or negative control were transfected into KYSE-150R cells. Finally, Dual-Luciferase Reporter System Kit (Promega) was used to detect the luciferase activity.

RNA pull-down

KYSE-150R cells (1.5×107) were collected and lysed. The biotinylated miR-217 probe was synthesized by Genepharm (Shanghai, China) and used for incubation with streptavidin agarose beads (Thermo Scientific). Cell lysate with miR-217 probe or oligo probe was incubated for one night at 4°C. RNA complex bound to the beads was eluted by wash buffer, and qRT-PCR was used to measure circRNA_100367 enrichment pulled down by miR-217 probe.

Xenografts experiments

This animal study was approved by the Ethics Committee of the First Affiliated Hospital of Zhengzhou University. KYSE-150 cells (5×106) transfected with circRNA_100367 or KYSE-150R cells (5×106) transfected with sh-circRNA_100367 were subcutaneously injected into the back of BALB/c female nude mice (6–8 weeks). All mice received radiotherapy 10 days after inoculation, with a radiation dose of 6Gy. So, the nude mice were divided into four groups, namely KYSE-150+Gy group (n=5), KYSE-150R+Gy group (n=5), KYSE-150+circRNA_100367+Gy group (n=5) and KYSE-150R+sh-circRNA_100367+Gy group (n=5). Tumor volume was measured 21 days later using the formula: tumor volume (mm3)=(length×width2)/2.

Statistical analysis

The data were presented as mean±standard deviation (SD), and were analyzed using SPSS software (version 20; Chicago, IL, USA). Log-rank test was used to compare the survival rates between high-expression Wnt3 group and low-expression Wnt3 group by Kaplan-Meier analysis The significance of difference was detected using t test or one-way analysis of variance (ANOVA), with p<0.05 considered significant difference.

Conflicts of Interest

All authors declare that they have no conflict of interest.

Funding

This study was supported by National Natural Science Foundation of China (No. 81703158).

References

  • 1. Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, Parkin DM, Forman D, Bray F. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015; 136:E359–86. https://doi.org/10.1002/ijc.29210 [PubMed]
  • 2. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018; 68:7–30. https://doi.org/10.3322/caac.21442 [PubMed]
  • 3. Chen WK, Chen CA, Chi CW, Li LH, Lin CP, Shieh HR, Hsu ML, Ko CC, Hwang JJ, Chen YJ. Moscatilin Inhibits Growth of Human Esophageal Cancer Xenograft and Sensitizes Cancer Cells to Radiotherapy. J Clin Med. 2019; 8:E187. https://doi.org/10.3390/jcm8020187 [PubMed]
  • 4. Jones CM, Spencer K, Hitchen C, Pelly T, Wood B, Hatfield P, Crellin A, Sebag-Montefiore D, Goody R, Crosby T, Radhakrishna G. Hypofractionated Radiotherapy in Oesophageal Cancer for Patients Unfit for Systemic Therapy: A Retrospective Single-Centre Analysis. Clin Oncol (R Coll Radiol). 2019; 31:356–64. https://doi.org/10.1016/j.clon.2019.01.010 [PubMed]
  • 5. Zhu H, Song H, Chen G, Yang X, Liu J, Ge Y, Lu J, Qin Q, Zhang C, Xu L, Di X, Cai J, Ma J, et al. eEF2K promotes progression and radioresistance of esophageal squamous cell carcinoma. Radiother Oncol. 2017; 124:439–47. https://doi.org/10.1016/j.radonc.2017.04.001 [PubMed]
  • 6. Chen M, Liu P, Chen Y, Chen Z, Shen M, Liu X, Li X, Li A, Lin Y, Yang R, Ni W, Zhou X, Zhang L, et al. Long Noncoding RNA FAM201A Mediates the Radiosensitivity of Esophageal Squamous Cell Cancer by Regulating ATM and mTOR Expression via miR-101. Front Genet. 2018; 9:611. https://doi.org/10.3389/fgene.2018.00611 [PubMed]
  • 7. Luo J, Wang W, Tang Y, Zhou D, Gao Y, Zhang Q, Zhou X, Zhu H, Xing L, Yu J. mRNA and methylation profiling of radioresistant esophageal cancer cells: the involvement of Sall2 in acquired aggressive phenotypes. J Cancer. 2017; 8:646–56. https://doi.org/10.7150/jca.15652 [PubMed]
  • 8. Wang F, Nazarali AJ, Ji S. Circular RNAs as potential biomarkers for cancer diagnosis and therapy. Am J Cancer Res. 2016; 6:1167–76. [PubMed]
  • 9. Yu D, Li Y, Ming Z, Wang H, Dong Z, Qiu L, Wang T. Comprehensive circular RNA expression profile in radiation-treated HeLa cells and analysis of radioresistance-related circRNAs. PeerJ. 2018; 6:e5011. https://doi.org/10.7717/peerj.5011 [PubMed]
  • 10. Chen Y, Yuan B, Wu Z, Dong Y, Zhang L, Zeng Z. Microarray profiling of circular RNAs and the potential regulatory role of hsa_circ_0071410 in the activated human hepatic stellate cell induced by irradiation. Gene. 2017; 629:35–42. https://doi.org/10.1016/j.gene.2017.07.078 [PubMed]
  • 11. He Y, Mingyan E, Wang C, Liu G, Shi M, Liu S. CircVRK1 regulates tumor progression and radioresistance in esophageal squamous cell carcinoma by regulating miR-624-3p/PTEN/PI3K/AKT signaling pathway. Int J Biol Macromol. 2019; 125:116–23. https://doi.org/10.1016/j.ijbiomac.2018.11.273 [PubMed]
  • 12. Su H, Lin F, Deng X, Shen L, Fang Y, Fei Z, Zhao L, Zhang X, Pan H, Xie D, Jin X, Xie C. Profiling and bioinformatics analyses reveal differential circular RNA expression in radioresistant esophageal cancer cells. J Transl Med. 2016; 14:225. https://doi.org/10.1186/s12967-016-0977-7 [PubMed]
  • 13. Yan J, Wu G, Chen J, Xiong L, Chen G, Li P. Downregulated miR-217 expression predicts a poor outcome in acute myeloid leukemia. Cancer Biomark. 2018; 22:73–78. https://doi.org/10.3233/CBM-170936 [PubMed]
  • 14. Chen Q, Wang P, Fu Y, Liu X, Xu W, Wei J, Gao W, Jiang K, Wu J, Miao Y. MicroRNA-217 inhibits cell proliferation, invasion and migration by targeting Tpd52l2 in human pancreatic adenocarcinoma. Oncol Rep. 2017; 38:3567–73. https://doi.org/10.3892/or.2017.6036 [PubMed]
  • 15. Bai M, Zhang M, Long F, Yu N, Zeng A, Wang X. MiR-217 promotes cutaneous squamous cell carcinoma progression by targeting PTRF. Am J Transl Res. 2017; 9:647–55. [PubMed]
  • 16. Li RQ, Ren Y, Liu W, Pan W, Xu FJ, Yang M. MicroRNA-mediated silence of onco-lncRNA MALAT1 in different ESCC cells via ligand-functionalized hydroxyl-rich nanovectors. Nanoscale. 2017; 9:2521–2530. https://doi.org/10.1039/c6nr09668a [PubMed]
  • 17. Xiao Y, Deng T, Su C, Shang Z. MicroRNA 217 inhibits cell proliferation and enhances chemosensitivity to doxorubicin in acute myeloid leukemia by targeting KRAS. Oncol Lett. 2017; 13:4986–94. https://doi.org/10.3892/ol.2017.6076 [PubMed]
  • 18. Lin CH, Tsai CH, Yeh CT, Liang JL, Hung WC, Lin FC, Chang WL, Li HY, Yao YC, Hsu TI, Lee YC, Wang YC, Sheu BS, et al. MiR-193a-5p/ERBB2 act as concurrent chemoradiation therapy response indicator of esophageal squamous cell carcinoma. Oncotarget. 2016; 7:39680–93. https://doi.org/10.18632/oncotarget.9444 [PubMed]
  • 19. Park M, Yoon HJ, Kang MC, Kwon J, Lee HW. MiR-338-5p enhances the radiosensitivity of esophageal squamous cell carcinoma by inducing apoptosis through targeting survivin. Sci Rep. 2017; 7:10932. https://doi.org/10.1038/s41598-017-10977-9 [PubMed]
  • 20. Wang HS, Nie X, Wu RB, Yuan HW, Ma YH, Liu XL, Zhang JY, Deng XL, Na Q, Jin HY, Bian YC, Gao YM, Wang YD, Chen WD. Downregulation of human Wnt3 in gastric cancer suppresses cell proliferation and induces apoptosis. Onco Targets Ther. 2016; 9:3849–60. https://doi.org/10.2147/OTT.S101782 [PubMed]
  • 21. Xing Z, Wang HY, Su WY, Liu YF, Wang XX, Zhan P, Lv TF, Song Y. Wnt3 knockdown sensitizes human non-small cell type lung cancer (NSCLC) cells to cisplatin via regulating the cell proliferation and apoptosis. Eur Rev Med Pharmacol Sci. 2018; 22:1323–32. [PubMed]
  • 22. He J, Zhou M, Chen X, Yue D, Yang L, Qin G, Zhang Z, Gao Q, Wang D, Zhang C, Huang L, Wang L, Zhang B, et al. Inhibition of SALL4 reduces tumorigenicity involving epithelial-mesenchymal transition via Wnt/β-catenin pathway in esophageal squamous cell carcinoma. J Exp Clin Cancer Res. 2016; 35:98. https://doi.org/10.1186/s13046-016-0378-z [PubMed]
  • 23. Oguma J, Ozawa S, Kazuno A, Nitta M, Ninomiya Y, Kajiwara H. Wnt3a expression is associated with poor prognosis of esophageal squamous cell carcinoma. Oncol Lett. 2018; 15:3100–08. https://doi.org/10.3892/ol.2017.7666 [PubMed]
  • 24. Zhao Y, Yi J, Tao L, Huang G, Chu X, Song H, Chen L. Wnt signaling induces radioresistance through upregulating HMGB1 in esophageal squamous cell carcinoma. Cell Death Dis. 2018; 9:433. https://doi.org/10.1038/s41419-018-0466-4 [PubMed]
  • 25. Zhen N, Gu S, Ma J, Zhu J, Yin M, Xu M, Wang J, Huang N, Cui Z, Bian Z, Sun F, Pan Q. CircHMGCS1 Promotes Hepatoblastoma Cell Proliferation by Regulating the IGF Signaling Pathway and Glutaminolysis. Theranostics. 2019; 9:900–19. https://doi.org/10.7150/thno.29515 [PubMed]
  • 26. Xu L, Feng X, Hao X, Wang P, Zhang Y, Zheng X, Li L, Ren S, Zhang M, Xu M. CircSETD3 (Hsa_circ_0000567) acts as a sponge for microRNA-421 inhibiting hepatocellular carcinoma growth. J Exp Clin Cancer Res. 2019; 38:98. https://doi.org/10.1186/s13046-019-1041-2 [PubMed]
  • 27. Zhang Y, Liu H, Li W, Yu J, Li J, Shen Z, Ye G, Qi X, Li G. CircRNA_100269 is downregulated in gastric cancer and suppresses tumor cell growth by targeting miR-630. Aging (Albany NY). 2017; 9:1585–94. https://doi.org/10.18632/aging.101254 [PubMed]
  • 28. Wang B, Qu XL, Liu J, Lu J, Zhou ZY. HOTAIR promotes osteosarcoma development by sponging miR-217 and targeting ZEB1. J Cell Physiol. 2019; 234:6173–6181. https://doi.org/10.1002/jcp.27394 [PubMed]
  • 29. Yin Z, Ren W. MicroRNA-217 acts as a tumor suppressor and correlates with the chemoresistance of cervical carcinoma to cisplatin. Onco Targets Ther. 2019; 12:759–71. https://doi.org/10.2147/OTT.S176618 [PubMed]
  • 30. Wu Y, Ginther C, Kim J, Mosher N, Chung S, Slamon D, Vadgama JV. Expression of Wnt3 activates Wnt/β-catenin pathway and promotes EMT-like phenotype in trastuzumab-resistant HER2-overexpressing breast cancer cells. Mol Cancer Res. 2012; 10:1597–606. https://doi.org/10.1158/1541-7786.MCR-12-0155-T [PubMed]
  • 31. Derksen PW, Tjin E, Meijer HP, Klok MD, MacGillavry HD, van Oers MH, Lokhorst HM, Bloem AC, Clevers H, Nusse R, van der Neut R, Spaargaren M, Pals ST. Illegitimate WNT signaling promotes proliferation of multiple myeloma cells. Proc Natl Acad Sci USA. 2004; 101:6122–27. https://doi.org/10.1073/pnas.0305855101 [PubMed]
  • 32. Su H, Wu Y, Fang Y, Shen L, Fei Z, Xie C, Chen M. MicroRNA-301a targets WNT1 to suppress cell proliferation and migration and enhance radiosensitivity in esophageal cancer cells. Oncol Rep. 2019; 41:599–607. https://doi.org/10.3892/or.2018.6799 [PubMed]
  • 33. Kureel J, John AA, Prakash R, Singh D. MiR 376c inhibits osteoblastogenesis by targeting Wnt3 and ARF-GEF-1 -facilitated augmentation of beta-catenin transactivation. J Cell Biochem. 2018; 119:3293–3303. https://doi.org/10.1002/jcb.26490 [PubMed]
  • 34. Chu Y, Fan W, Guo W, Zhang Y, Wang L, Guo L, Duan X, Wei J, Xu G. miR-1247-5p functions as a tumor suppressor in human hepatocellular carcinoma by targeting Wnt3. Oncol Rep. 2017; 38:343–51. https://doi.org/10.3892/or.2017.5702 [PubMed]
  • 35. Jing Z, Gong L, Xie CY, Zhang L, Su HF, Deng X, Wu SX. Reverse resistance to radiation in KYSE-150R esophageal carcinoma cell after epidermal growth factor receptor signal pathway inhibition by cetuximab. Radiother Oncol. 2009; 93:468–73. https://doi.org/10.1016/j.radonc.2009.08.008 [PubMed]
  • 36. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001; 25:402–08. https://doi.org/10.1006/meth.2001.1262 [PubMed]