The targeted regulation of Gli1 by miR-361 to inhibit epithelia-mesenchymal transition and invasion of esophageal carcinoma cells

Abstract

This article has been retracted, and the online PDF replaced with this retraction notice.

Keywords

MicroRNA-361 Gli1 Snail EMT esophageal cancer invasion

1. Introduction

Esophageal cancer represented a type of malignant tumors with high incidence and mortality rate worldwide. Its occurrence rate ranks 8th among all sorts of cancers, and its mortality rate is 6th highest all over the world [1]. Esophageal cancer is also a popular tumor affecting digestive tract in China, as approximately 50% of all cases have been found so far in China. Esophageal squamous cell carcinoma (ESCC) is the most common pathological type, occupying 90% of all cases [2]. Due to the lack of sensitive and early diagnostic marker, most esophageal cancer patients were already at advanced or terminal phase at the first time of diagnosis, leading to ineffective treatment and unfavorable prognosis, with 5-year survival rate about merely 10% $\sim$ 20% [3]. Early diagnosis is thus one of the most effective pathways to improve the overall survival rate of cancer patients. With timely diagnosis and treatment, the 5-year survival rate of esophageal cancer patients will be increased to 90%. In this regard, the study of pathogenesis of esophageal cancer will benefit early diagnosis and screen potential targets for the treatment. Metastasis and recurrence often occur in esophageal patients even after surgical removal, severely compromising the treatment efficacy. Current study confirmed the epithelial mesenchymal transition (EMT) as the initial step of tumor invasion and metastasis, and was closely related with unfavorable prognosis [4]. Hedgehog signal pathway played a critical role in human embryonic development [5, 6], and participates in the modulation of cell proliferation [7], migration, motility [8] and differentiation [9]. Various studies showed the close relationship between the occurrence and progression of aberrant activation of Hedgehog signal pathway [10, 11, 12]. Nuclear transcriptional factor, glioma-associated oncogene homolog 1 (Gli1), exerted the function as a transcription effector, as it can initiate transcription of nuclear genes, thus transducing and activating Hedgehog signal pathway via transforming extracellular Hedgehog signal to intracellular Gli1 signal [13, 14]. Study showed the correlation between Gli1 expression or functional enhancement of Gli1 and diverse types of tumors, including liver cancer [15], breast carcinoma [16], pancreatic cancer [17], ovarian tumors [18]. Its role in invasion or metastasis in esophageal carcinoma, however, remains elusive. Moreover, as in silico bioinformatics analysis revealed the targeted complementary relationship between miR-361 and 3’-UTR of Gli1, this study thus investigated the potential role of miR-361 on Gli1 in the invasion and metastasis of esophageal carcinoma.

2. Materials and methods

2.1 Major reagents and materials

Esophageal carcinoma cell line EC9706 was obtained from ATCC (US). DMEM, RPMI1640, FBS and penicillin-streptomycin were purchased from Gibco (US). Trizol and Lipofectamine 2000 were purchased from Invitrogen (US). Reverse transcription kit was purchased from Takara (China). SYBR Green real-time PCR Master Mixes were purchased from Life Technologies (US). Fragments of miR-361 nucleotide acids (including mimic, inhibitor and negative control, NC) and PCR primers were synthesized by Ruibo Bio (China). Mouse anti-human E-cadherin and rabbit anti-human N-cadherin antibody were purchased from Santa Cruz (US). BAC protein quantification assay kit was purchased from Beyotime (China). Transwell chamber was purchased from Corning (US). Matrigel was purchased from BD biosciences (US). Dual-luciferase reporter assay system and pGL3-promoter plasmid were purchased from Promega (US).

2.2 Clinical information

A total of 58 esophageal carcinoma patients who received pathological diagnosis and treatment in Tianjin Medical University Cancer Hospital from December 2014 to November 2015 were enrolled in this study. All patients had ESCC, including 30 males and 28 females, plus the average age at 56.7 $\pm$ 10.3 years. All patients were for the first time treated with surgical resection of tumors, nor did receive any radio- or chemo-therapy before surgery. The clinical phases were classified based on TNM standard stipulated by Union for International Cancer Control (UICC) and American Joint Committee on Cancer (AJCC). There were 19, 22 and 17 patients at phase 0 $\sim$ I, phase IIa $\sim$ IIb and phase III, respectively. All tissues samples were collected within 10 min of tumor extraction. Adjacent tissues were collected from normal esophageal mucosa over 5 cm from tumor edge. All samples were frozen in liquid nitrogen, and were transferred to $-$ 80 ${}^{\circ}$ C for storage.

This study has been pre-approved by the ethical committee of Tianjin Medical University Cancer Hospital. All subjects have signed the consent forms before recruitment in this study.

2.3 Culture of EC9706 cells

EC9706 cells were cultured with RPMI1640 medium containing 10% FBS, 100 U/mL penicillin and 100 $\mu$ g/mL streptomycin, and were incubated in a chamber at 37 ${}^{\circ}$ C with 5% CO ${}_{2}$ perfusion. Cells were performed in further experiments when reaching a confluence of 70% $\sim$ 80%.

2.4 Prediction of target gene

Under the prediction of microRNA.org, there were satisfactory targeted complementary relationship between miR361, miR-873, miR-875 and 3’-UTR of Gli1 gene.

2.5 Construction of luciferase reporter assay plasmid

Using HEK293 genome as the template, full length fragment of 3’-UTR of Gli1 gene was amplified (Forward primer, 5’-TGTGA TGGAT GAGAT TC CCT ACTCT-3’; Reverse primer, 5’-TCCTG ATAAT AAAGG AACTG CATCA-3’). PCR products were purified from agarose gel, and were ligated into pGL-3M luciferase reporter plasmid after XbaI/NotI dual digestion. Recombinant plasmid was then used to transform DH5 $\alpha$ competent cells. Positive clones with primary screening were selected for sequencing (Sangon, China). Those plasmids with correct sequence were named as pGL3-Gli1-3’UTR for further cell transfection and following experiments.

2.6 Luciferase reporter assay

HEK293 cells were inoculated in DMEM culture medium containing 10% FBS and 1% penicillin-streptomycin. One day before transfection, cells with satisfactory growth condition at log phase were seeded into 24-well plate at 4 $\times$ 104/cm ${}^{2}$ . When cells reached the density of 60%, Lipofectamine 2000 was used to transfect cells with 100 ng pGL3-Gli1-3’UTR plasmid, 50 nmol miR-361 mimic (or miR-361 inhibitor/NC), and 50 ng internal controlled plasmid pRL-TK. After 6 h transfection, Opti-MEM medium without antibiotics or serum were discarded, with the replacement of normal DMEM medium containing serum. For 48 h continuous incubation, dual-luciferase assay was performed. In brief, culture medium was discarded from the plate, which was washed twice in PBS, with the addition of 100 $\mu$ L passive lysis buffer. With vortex at room temperature for 15 min, the mixture was centrifuged at 300 g for 5 min. 50 $\mu$ L supernatant was kept to mix with luciferase substrate pre-diluted in luciferase assay buffer II, followed by measurement in chemiluminescence apparatus immediately. The enzymatic reaction of firefly luciferase was terminated in 50 $\mu$ L quenching buffer, which activated Renilla luciferase simultaneously. The activity of the latter enzyme was measured in a microplate reader. The relative expression level of reporter gene was calculated by the ratio of luciferase activity against Renilla luciferase activity. Nucleotide acid sequences were: mimic NC, 5’-UUCUC CGAAC GUGUC ACGUT T-3’; miR-361 mimic, 5’-ACCCC UGGAG AUUCU GAUAA UU-3’; inhibitor NC, 5’-CAGGU CAUAG GCAUA CCGU-3’; and miR-361 inhibitor, 5’-ACCUU AGCAG GAUUA CUAG-3’.

2.7 Immunohistochemistry

Esophageal cancer tissue and adjacent tissue were collected for the preparation of paraffin sections (thickness 8 m). After being dewaxed and rinsed three times with pH7.4 PBS, each slice was added with 1 drop of 3% H ${}_{2}$ O ${}_{2}$ for the incubation at room temperature for 10 min. Rabbit polyclonal anti-human Gli1 receptor antibody (1:1000 dilution) was incubated for 2 h incubation, and polymer reinforcing agent was added for another 20 min incubation. The enzyme-labeled anti-human/rabbit polymer was added for the incubation of 30 min and freshly prepared DAB solution (two amino benzidine) was used for observation by 5 min microscope. After process of hematoxylin staining, 0.1% HCl differentiation, water rinse, dehydration with gradient alcohol, transparent with xylene, sealed with neutral gum, the slice was under observation (Olympus, IX71, Shinjuku, Tokyo, Japan). Immunohistochemical method (Ultra, Vision, Detection, System, three step methods) was used to determine the level of Gli1 protein with Image-pro plus image analysis system (Media Cybernetics, Inc, Rockville, MD, USA).

2.8 Construction of recombinant plasmid for Gli1 overexpression

pGC-FU vector was used to load the fragment of Gli1 gene into the cleavage site of Age I. Recombinant plasmid pGC-FU-Gli1 was constructed for the over-expression of Gli1 in cells after transfection.

2.9 Transfection of EC9706 cells

Using human Gli1 gene as the template, the small interference sequence targeting Gli1 gene was synthesized by Gimma (China). Sequences were: si-Gli1 sense: 5’-CCAGG AAUUU GACUC CCAAT T-3’; si-Gli1 anti-sense: 5’-UUGGG AGUCA AAUUC CUGGC T-3’. Experiments were performed in seven parallel groups, namely, mimic NC, miR-361 mimic, si-NC, si-Gli1, miR-361 mimic+si-Gli1, pGC-FU and pGC-FU-Gli1 groups. After 48 h, cells were collected for further experiments.

2.10 RT-PCR for gene expression assay

Cells were lysed by Trizol method, and total RNA was extracted by chloroform. cDNA was synthesized by PrimerScript RT reagent kit by reverse transcription in a 20 $\mu$ L system including 2 $\mu$ g total RNA, 1 $\mu$ L dNTP (10 mmol/L), 4 $\mu$ L RT buffer (5 $\times$ ), 1 $\mu$ L RT primer (2 $\mu$ mol/L), 2 $\mu$ L reverse transcriptase, 1 $\mu$ L RNase inhibitor and ddH ${}_{2}$ O. The reaction conditions were: 42 ${}^{\circ}$ C for 15 min, followed by 85 ${}^{\circ}$ C 5 s. cDNA products were kept at $-$ 20 ${}^{\circ}$ C fridge. Using cDNA as the template, PCR amplification was performed under the direction of TaqDNA polymerase using primers (miR-361PF: 5’-ATAAA GTGCT GACAG TGCAG ATAGT G-3’; miR-361PR: 5’-TCAAG TACCC ACAGT GCGGT-3’; U6PF: 5’-TTGGA ACGAT ACAGA GAAGA TT-3’; U6PR: 5’-GGAAC GCTTC ACGA ATTT G-3’; Gli1PF: 5’-AGCGT GAGCC TGAAT CTGTG-3’; Gli1PR: 5’-CAGCA TGTAC TGGGC TTTGA A-3’; SnailPF: 5’-ACCCC ACATC CTTCT CACTG-3’; SnailPR: 5’-TACAA AAACC CACGC AGACA-3’; $\beta$ -actinPF: 5’-GAACC CTAAG GCCAA C-3’; $\beta$ -actinPR: 5’-TGTCA CGCAC GATTT CC-3’; In a PCR system with 10 $\mu$ L total volume, we added 4.5 $\mu$ L 2 $\times$ SYBR Green Mixture, 0.5 $\mu$ L of forward/reverse primer (at 2.5 $\mu$ m/L), 1 $\mu$ L cDNA, and 3.5 $\mu$ L ddH ${}_{2}$ O. PCR conditions were: 95 ${}^{\circ}$ C for 15 s, 60 ${}^{\circ}$ C for 30 s and 74 ${}^{\circ}$ C for 30 s. The reaction was performed on ABI 7500 fluorescent quantitative PCR cycler for 40 cycles. Each sample was tested in triplicates. Data were analyzed by comparative Ct method (2 ${}^{-\Delta\Delta\text{Ct}}$ ).

2.11 Western blotting

Cells were collected and lysed using RIPA buffer for extracting total proteins. BCA method was used to determine the protein concentration. Proteins were diluted using three volumes of 4 $\times$ loading buffer, and were boiled for 5 min. 50 $\mu$ g samples were loaded for separation at 10% SDS-PAGE for 2.5 h $\sim$ 3 h. Proteins were then transferred to PVDF membrane under 60 min electrical field. The membrane was then blocked by 5% defatted milk powder for 60 min, followed by primary antibody incubation at 4 ${}^{\circ}$ C overnight. PBST was used to wash the membrane for three times to remove unbounded antibody. Secondary antibody with HRP labels was then added for 60 min room temperature incubation, followed by three changes of PBST washing. ECL reagent was then added for room temperature incubation (1 $\sim$ 3 min) on the membrane, followed by exposure, development in a dark room. Image analysis was performed using Quantity One software.

2.12 Plate clonal formation assay for malignant growth of cells

Cells at log growth phase were digested in 0.25% trypsin and prepared into single cell suspension using RPMI1640 medium with 10% FBS. Cells were then inoculated into 10 cm plate containing 10 mL culture medium at 100 cells per dish density. The dish was gently turned to spread cells evenly, and were incubated in a humidified chamber perfused with 5% CO ${}_{2}$ at 37 ${}^{\circ}$ C for 14 $\sim$ 21 days. Incubation was stopped when visible clones were formed in the culture dish. The medium was discarded. Cells were rinsed twice in PBS, and fixed by 4% paraformaldehyde for 15 min at room temperature. Giemsa staining buffer was added for 20 min incubation after discarding the fixation buffer. The staining was removed by slow tap water and the plate was dried in air. The number of clones with more than 10 cells was counted under 4 $\times$ objectives. Clonal formation rate was calculated as (clone number/inoculated cell number) $\times$ 100%.

2.13 Transwell assay for cell invasion

100 $\mu$ L Matrigel was paved on the upper surface of Transwell chamber (8 $\mu$ m pore size). The chamber was incubated at 37 ${}^{\circ}$ C for 30 min. 200 $\mu$ L cell suspension (1 $\times$ 10 ${}^{6}$ /mL) in serum-free medium was added into the upper chamber, while the lower chamber was filled with 600 $\mu$ L complete medium with 10% FBS. After 48 h incubation, liquid in the upper chamber was discarded. Matrigel and unpenetrated cells were removed by sterilized cotton swabs, followed by methanol fixation for 20 min. The membrane was washed twice in PBS and air-dried, and was stained by 0.1% crystal violet for 30 min. Under an inverted microscope, five randomly selected fields were counted for the number of cells penetrating the micro-pores and averaged.

2.14 Statistical analysis

SPSS18.0 software was used for data analysis. Measurement data were presented as mean $\pm$ standard deviation (SD), and were compared by student $t$ -test. A statistical significance was defined when $p<$ 0.05.

Figure 1.

Negative correlation between Gli1 and miR-361 expression in esophageal cancer tissues. (A) Gli1 mRNA expression by qRT-PCR; (B) miRNA expression by qRT-PCR; (C) Gli1 protein expression by Western blot. ${}^{***}$ , $p<$ 0.001.

Figure 2.

miR-361 specifically inhibited Gli1 gene expression. (A) Dual luciferase assay; (B) miR-361 expression by qRT-PCR; (C) Intracellular Gli mRNA expression in EC9706 cells; (D) Western blotting for Glia protein expression; (E) binding sites between miR-361 and Gli1 3’-UTR; (E) Immunohistochemistry for Glia protein expression. ${}^{***}$ , $p<$ 0.001.

3. Results

3.1 Gli1 expression in esophageal carcinoma tissues was negatively correlated with miR-361

qRT-PCR results showed that, mRNA and protein levels of Gli1 in tumor tissues of esophageal carcinoma patients were significantly elevated compared to that in adjacent tissues. Such expression levels were gradually increased with advanced TNM stage, indicating the potential involvement of Gli1 up-regulation in pathogenesis of esophageal cancer ( $p<$ 0.001) (Fig. 1A, C and D). Bio-informatics analysis showed satisfactory targeted complementary relationship between miR-361, miR-873, miR-875 and 3’-UTR of Gli1 gene. Our results showed no significant difference of miR0873 or miR-875 expression level between tumor tissue and adjacent tissue, suggesting that these two miRs might not participate in the pathogenesis of esophageal cancer (Fig. 1B). The expression level of miR-361, however, was significantly decreased in tumor tissues, with gradually decreased expression with advanced TNM stage (Fig. 1B). Correlation analysis revealed significantly negative correlation between miR-361 level and Gli1 mRNA (r $=-$ 0.854, $p=$ 0.032), suggesting possible targeted regulation between those two factors. The elevation of Gli1 might be caused by the decrease of miR-361 expression, indicating that the disruption of function and expression between miR-361 and Gli1 might participate in modulating features of esophageal cancer pathogenesis and invasion/metastasis.

3.2 miR-361 specifically inhibited Gli1 expression

As shown in Fig. 2A, the transfection of miR-361 mimic significantly depressed relative luciferase activity of HEK293 cell lysate ( $p<$ 0.001), suggesting that miR-361 could specifically targeted 3’-UTR of Gli1 gene (Fig. 2E), further suppressing expression level of luciferase-Gli1-3’UTR mosaic gene, thus decreasing luciferase activity. Further potentiating or inhibiting miR-361 expression (Fig. 2B) remarkably decreased or increased Gli1 gene and protein expression level (Fig. 2C and D) in EC9706 cells, respectively ( $p<$ 0.001). These results suggested the specific inhibition of miR-361 on Gli1 expression, the target gene of which has been confirmed as Gli1.

Figure 3.

The specific inhibition of miR-361 on the malignant growth and invasion potency of esophageal cancer cells. (A) qRT-PCR for Gli1 mRNA expression in EC9706 cells; (B) Plate clonal formation assay for the ability of malignant growth; (C) Snail mRNA expression of EC9706 cells by qRT-PCR; (D) Protein expression of EC9706 cells by Western blotting; (E) Transwell assay for cell invasion ability. ${}^{*}$ , $p<$ 0.05; ${}^{**}$ , $p<$ 0.01; ${}^{***}$ , $p<$ 0.001.

3.3 miR-361 specifically inhibited Gli1 and decreased malignant growth and invasion potency of esophageal carcinoma cells

As a critical transcriptional factor in Hedgehog signal pathway, Gli1 participated in the modulation of tumor cell EMT process via affecting the regulation of Snail on the expression of E-cadherin and N-cadherin, and was thus correlated with invasion and metastasis of various tumors including hepatocellular carcinoma, gastric cancer and pancreatic carcinoma. Examination of clinical samples revealed higher Gli1 expression level in esophageal cancer tissues with advanced TNM stage, accompanied with lower miR-361. This study aimed to investigate if the imbalance of Gli1 and miR-361 was correlated with malignant growth, invasion and metastasis features of cancer cells. The transfection of miR-361 mimics and/or si-Gli1 significantly depressed mRNA and protein expression of Gli1 in EC9706 cells, whereas the over-expression of Gli1 was also found after transfection with pGC-FU-Gli1 (Fig. 3A and D), demonstrating the satisfactory transfection efficiency for further experiments. Results from clonal formation assay showed weakened ability of EC9706 cells after up-regulating miR-361 and/or silencing Gli1, as the combined transfection had most potent effects. In contrast, the over-expression of Gli1 enhanced the effect of clonal formation (Fig. 3B), indicating the involvement of abnormal decrease (increase) of miR-361 and/or Gli1 expression in the malignant growth of esophageal cancer cells. The antagonizing of such effect could weaken anchor-independent growth of cells. The malignant phenotype of cells was thus inhibited to certain extents. The elevation of miR-361 and/or decrease of Gli1 expression also significantly inhibited the invasiveness of esophageal cancer cells (Fig. 3E), decrease intracellular expression of Snail and N-cadherin, and elevated the expression of E-cadherin. Nevertheless, the over-expression of Gli1 promoted the invasiveness of esophageal cancer cells, elevated the levels of Snail and N-cadherin but reduced the expression of E-cadherin (Fig. 3C and D).

4. Discussion

Esophageal carcinoma is one of the most common malignant tumors in China. It was characterized as the cancer with high incidence, insidious onset, propensity for invasion or metastasis, and significant mortality [19]. The study of its pathogenesis mechanism was of critical importance for the early diagnosis and treatment of esophageal cancer. Recent study found the critical role of various signal pathways related with embryonic and tissue development in the pathogenesis of tumors, such as Hedgehog pathway [5, 6]. Mammalian Hedgehog signal pathway possessed three extracellular Hedgehog signal proteins, namely Sonic hedgehog (Shh), Indian hedgehog (Ihh) and Desert hedgehog (Dhh), among which Shh ligand has the widest distribution and powerful functional roles [20]. Within Hedgehog signal pathway, there were extracellular signal ligand HH, transmembrane protein receptor Ptch, alternative transmembrane protein Smo, intermediate transmitting molecule and transcriptional factor Gli. Transmembrane protein Smoothened (Smo), as the information transducer, can convert extracellular HH signal into intracellular Gli signal for initiating nuclear gene transcription [21]. When Hedgehog signal pathway was abnormally activated, Smo was facilitated to enter the cytoplasm, where it transformed extracellular Shh signal into intracellular Gli signal, thus inducing transcription of multiple oncogenes, leading to cell transformation and malignancy, and eventually tumor occurrence [22, 23]. Gli1 belonged to Gli homolog genes in human (Gli1, Gli2 and Gli3), it can modulate downstream gene transcription via variable self length, and affecting various biological processes including cell proliferation [24], apoptosis [25] and differentiation [26]. Study showed that the abnormal elevation of Gli1 in facilitating invasion and metastasis of liver cancer [15], breast cancer [16], pancreatic carcinoma [17] and ovarian tumors [18]. Its relationship with invasion and metastasis of esophageal cancer, however, is still unknown. This study found significantly elevated mRNA and protein expressions of Gli1 in esophageal carcinoma tissues. Moreover, those tissues with advanced TNM stage had significantly higher Gli1 expression than those with lower TNM stage, indicating the potential tumor facilitating role of Gli in esophageal cancer.

On-line prediction by microRNA.org revealed complementary correlation between miR-361, miR-873, miR-875 and 3’-UTR of Gli1 gene. Its abnormal expression thus may affect Gli1 expression and further accelerate invasion and metastasis of esophageal cancer. Results of this study showed no significant difference between the expression of miR-873 and miR-875 in esophageal cancer and tumor adjacent tissues, indicating its absence in cancer pathogenesis. The expression of miR-361 was gradually decreased with advanced TNM stage, and was significantly correlated with Gli1 expression level. Further luciferase reporter gene assay confirmed the targeted inhibitory role of miR-361 on mRNA expression of Gli1, suggesting that miR-361 might regulate Gli1 and function as a tumor suppressor gene in esophageal cancer. Sun et al. [27] found that miR-361 possibly regulate CXCR6 to inhibit the proliferation of liver cancer cells. Liu et al. [28] found that miR-361 could regulate STAT6 to function as tumor suppressor gene in prostate cancer. Wu et al. [29] found that the function of miR-361 could decrease of recurrent rate of gastric cancer via interplay with 3’-UTR of CD80 gene. Ma et al. revealed that miR-361 could specifically inhibit SND1 gene expression and growth of colorectal tumor body. All these studies indicated the tumor-suppressor role of miR-361 against multiple tumors. This study focused on the biological regulatory role of miR-361 in esophageal cancer. Colony formation assay revealed that, after the up-regulation of miR-361 expression by mimic transfection, the malignant growth of EC9706 cells was significantly decreased. Similar effects can be also obtained by interfering endogenous Gli1 expression. The combine treatment using these two factors had more potent inhibitory effects on malignant growth of cells, suggesting that miR-361 modifies biological effect of esophageal cancer via mediating Gli1 expression.

Invasion and metastasis were major reasons for post-op recurrence of esophageal cancer, and were also important factors affecting treatment efficacy and patient prognosis [30]. EMT occurrence was known to be closely correlated with tumor invasion and metastasis [4]. E-cadherin localized at intra-cellular junction, and consisted of N-terminal extracellular domain, hydrophobic transmembrane domain and C-terminal cytosolic domain. Via terminal repeated sequences of extracellular domain, it can form dimer with another E-cadherin molecule. Whilst intracellular domain can form complex with beta-connectin, which directly connects with cellular cytoskeleton actin and forms zipper like structure to maintain epithelial polarity and cell-to-cell adhesion stability. The decrease or absence of E-cadherin will induce EMT occurrence and facilitate cell migration and invasion [31]. N-cadherin and E-cadherin both belong to cadherin family. The up-regulation of N-cadherin induced mesenchymal-fibrosis transition of epithelial, benefiting cellular motility and migration [32]. The down-regulation of E-cadherin and up-regulation of N-cadherin were critical events during EMT process. Transcription factor Snail was a key modulator for EMT, as it bound onto E-box in E-cadherin gene promoter for inhibiting its transcription and expression, thus facilitating EMT occurrence and cell migration [33]. Moreover, Snail can initiate EMT process via up-regulating N-cadherin expression [34]. Previous study found that Gli1 could modulate expressions of E-cadherin and N-cadherin via affecting Snail [35], and participate in EMT regulation of tumor cells, thus being correlated with invasion and metastasis of various tumors such as liver cancer [15], gastric carcinoma [36] and breast cancer [37]. Moreover, abnormal expression of miR-361 has also been reported to be correlated with metastasis of colon cancer, suggesting its role in modulating invasion and metastasis of tumor cells [38]. This study investigated the relationship between miR-361 and invasiveness of esophageal carcinoma, along with possible role of Gli1. Results showed that the elevation of miR-361 and/or decrease of Gli1 expression significantly decreased Snail and N-cadherin expression, and potentiated E-cadherin protein expression, in addition to the inhibition of invasiveness of esophageal carcinoma cells. In contrast, an opposite result was found when Gli1 was over-expressed. Our data thus showed that miR-361 could impede EMT process and weaken invasiveness of esophageal carcinoma cells via specific inhibition of Gli1 gene and induced Snail expression.

5. Conclusion

In esophageal carcinoma tissues, Gli1 expression was abnormally elevated while miR-361 was down-regulated. MiR-361 played as an oncogene during the pathogenesis of esophageal cancer. It can impede EMT process and weaken invasiveness of esophageal cancer cells via specific inhibition on Gli1 and its induced Snail gene expression.

Footnotes

Acknowledgments

This work was supported by Chengde Science and Technology Bureau (201701A095).

Conflict of interest

None.

References

Matsushima

et al., MicroRNAs and esophageal squamous cell carcinoma, Digestion 82(3) (2010), 138–144.

Siegel

et al., Cancer statistics, 2014, CA Cancer J Clin 64(1) (2014), 9–29.

Wen

and Hu

, Expression of esophageal carcinoma related gene 4 (ECRG4) and its clinical significance in prognosis of esophageal carcinoma, Int J Clin Exp Pathol 8(11) (2015), 14772–14778.

et al., MiR-133a suppresses the migration and invasion of esophageal cancer cells by targeting the EMT regulator SOX4, Am J Transl Res 7(8) (2015), 1390–1403.

Lee

M.Y.

Sun

and Veltmaat

J.M.

, Hedgehog and Gli signaling in embryonic mammary gland development, J Mammary Gland Biol Neoplasia 18(2) (2013), 133–138.

Nagase

et al., Hedgehog signaling: a biophysical or biomechanical modulator in embryonic development? Ann N Y Acad Sci 1101 (2007), 412–438.

Yang

et al., Sonic hedgehog overexpression regulates the neuroepithelial cells proliferation in the spinal cord of dorsal regions during chicken embryo development, Neuro Endocrinol Lett 36(4) (2015), 380–386.

Chang

et al., Activation of sonic hedgehog signaling enhances cell migration and invasion by induction of matrix metalloproteinase-2 and -9 via the phosphoinositide-3 kinase/AKT signaling pathway in glioblastoma, Mol Med Rep 12(5) (2015), 6702–6710.

Villarreal

C.M.

et al., Hedgehog signaling regulates imaginal cell differentiation in a basally branching holometabolous insect, Dev Biol 404(2) (2015), 125–135.

10.

et al., Hedgehog signaling is activated in subsets of esophageal cancers, Int J Cancer 118(1) (2006), 139–148.

11.

Kebenko

et al., ErbB2 signaling activates the Hedgehog pathway via PI3K-Akt in human esophageal adenocarcinoma: identification of novel targets for concerted therapy concepts, Cell Signal 27(2) (2015), 373–381.

12.

Yang

et al., Identification of signature genes for detecting hedgehog pathway activation in esophageal cancer, Pathol Oncol Res 17(2) (2011), 387–391.

13.

Gruber

et al., DYRK1B as therapeutic target in Hedgehog/GLI-dependent cancer cells with Smoothened inhibitor resistance, Oncotarget 7(6) (2016), 7134–7148.

14.

Song

et al., Crosstalk between Wnt/beta-catenin and Hedgehog/Gli signaling pathways in colon cancer and implications for therapy, Cancer Biol Ther 16(1) (2015), 1–7.

15.

Chen

J.S.

et al., Down-regulation of Gli-1 inhibits hepatocellular carcinoma cell migration and invasion, Mol Cell Biochem 393(1–2) (2014), 283–291.

16.

Lei

et al., Gli-1 is crucial for hypoxia-induced epithelial-mesenchymal transition and invasion of breast cancer, Tumour Biol 36(4) (2015), 3119–3126.

17.

Nagai

et al., Gli1 contributes to the invasiveness of pancreatic cancer through matrix metalloproteinase-9 activation, Cancer Sci 99(7) (2008), 1377–1384.

18.

Chen

et al., Down-regulation of Gli transcription factor leads to the inhibition of migration and invasion of ovarian cancer cells via integrin beta4-mediated FAK signaling, PLoS One 9(2) (2014), e88386.

19.

Liu

et al., MMP1 promotes tumor growth and metastasis in esophageal squamous cell carcinoma, Cancer Lett 377(1) (2016), 97–104.

20.

Pyczek

et al., Hedgehog signaling activation induces stem cell proliferation and hormone release in the adult pituitary gland, Sci Rep 6 (2016), 24928.

21.

Philipp

and Caron

M.G.

, Hedgehog signaling: is Smo a G protein-coupled receptor? Curr Biol 19(3) (2009), R125–7.

22.

Bai

et al., Hedgehog Signaling in Pancreatic Fibrosis and Cancer, Medicine (Baltimore) 95(10) (2016), e2996.

23.

Hanna

and Shevde

L.A.

, Hedgehog signaling: modulation of cancer properies and tumor mircroenvironment, Mol Cancer 15 (2016), 24.

24.

Bermudez

et al., Gli1 mediates lung cancer cell proliferation and Sonic Hedgehog-dependent mesenchymal cell activation, PLoS One 8(5) (2013), e63226.

25.

Wang

et al., Gli1 inhibition induces cell-cycle arrest and enhanced apoptosis in brain glioma cell lines, J Neurooncol 98(3) (2010), 319–327.

26.

Fiaschi

et al., Targeted expression of GLI1 in the salivary glands results in an altered differentiation program and hyperplasia, Am J Pathol 179(5) (2011), 2569–2579.

27.

Sun

J.J.

Chen

G.Y.

and Xie

Z.T.

, MicroRNA-361-5p Inhibits Cancer Cell Growth by Targeting CXCR6 in Hepatocellular Carcinoma, Cell Physiol Biochem 38(2) (2016), 777–785.

28.

Liu

et al., MiR-361-5p acts as a tumor suppressor in prostate cancer by targeting signal transducer and activator of transcription-6 (STAT6), Biochem Biophys Res Commun 445(1) (2014), 151–156.

29.

et al., A functional variant at miR-132-3p, miR-212-3p, and miR-361-5p binding site in CD80 gene alters susceptibility to gastric cancer in a Chinese Han population, Med Oncol 31(8) (2014), 60.

30.

Zhu

et al., Activating transcription factor 4 promotes esophageal squamous cell carcinoma invasion and metastasis in mice and is associated with poor prognosis in human patients, PLoS One 9(7) (2014), e103882.

31.

Niessen

C.M.

and Gottardi

C.J.

, Molecular components of the adherens junction, Biochim Biophys Acta 1778(3) (2008), 562–571.

32.

Chen

Munoz-Antonia

and Cress

W.D.

, Trim28 contributes to EMT via regulation of E-cadherin and N-cadherin in lung cancer cell lines, PLoS One 9(7) (2014), e101040.

33.

Cano

et al., The transcription factor snail controls epithelial-mesenchymal transitions by repressing E-cadherin expression, Nat Cell Biol 2(2) (2000), 76–83.

34.

Moreno-Bueno

et al., Genetic profiling of epithelial cells expressing E-cadherin repressors reveals a distinct role for Snail, Slug, and E47 factors in epithelial-mesenchymal transition, Cancer Res 66(19) (2006), 9543–9556.

35.

Min

et al., The glioma-associated oncogene homolog 1 promotes epithelial-mesenchymal transition in human esophageal squamous cell cancer by inhibiting E-cadherin via Snail, Cancer Gene Ther 20(7) (2013), 379–385.

36.

Gao

et al., Resveratrol inhibits the hedgehog signaling pathway and epithelial-mesenchymal transition and suppresses gastric cancer invasion and metastasis, Oncol Lett 9(5) (2015), 2381–2387.

37.

Wang

Z.S.

et al., Significance and prognostic value of Gli-1 and Snail/E-cadherin expression in progressive gastric cancer, Tumour Biol 35(2) (2014), 1357–1363.

38.

et al., MiR-361-5p inhibits colorectal and gastric cancer growth and metastasis by targeting staphylococcal nuclease domain containing-1, Oncotarget 6(19) (2015), 17404–17416.