Abstract
Huihua Cai1,*, Jie Yao2,*, Yong An1, Xuemin Chen1, Weibo Chen1, Di Wu1, Boyang Luo1, Yong Yang1, Yong Jiang1, Donglin Sun1, Xiaozhou He3
1Department of Hepatobiliary Surgery, The First People's Hospital of Changzhou, The Third Hospital Affiliated to Soochow University, Changzhou, Jiangsu, China
2Department of Hepatobiliary and Pancreatic Surgery, Northern Jiangsu People's Hospital, The Clinic Medical College of Yangzhou University, Yangzhou, Jiangsu, China
3Department of Urology, The First People's Hospital of Changzhou, The Third Hospital Affiliated to Soochow University, Changzhou, Jiangsu, China
*These authors contributed equally to this work
Correspondence to:
Xiaozhou He, email: [email protected]
Donglin Sun, email: [email protected]
Keywords: pancreatic cancer, miR-613, HOTAIR, cell proliferation, invasion and migration
Received: February 10, 2017 Accepted: March 11, 2017 Published: March 22, 2017
ABSTRACT
Pancreatic cancer is one of the most deadly cancers with a poor prognosis. Though studies have implicated the roles of microRNAs in pancreatic cancer progression, little is known about the role of miR-613 in pancreatic cancer. In the present study, the expression of miR-613 was down-regulated in pancreatic cancer tissues and cancer cell lines. Down-regulation of miR-613 was positively correlated with tumor differentiation, advanced TNM stage, nodal metastasis and shorter overall survival in patients with pancreatic cancer. Overexpression of miR-613 suppressed cell proliferation, invasion and migration, and induced cell apoptosis and cell cycle arrest at G0/G1 phase in pancreatic cancer cells. Bioinformatics analysis, luciferase reporter assay and rescue experiments showed that notch3 was a direct target of miR-613. MiR-613 was inversely correlated with notch3 expression in pancreatic cancer tissues. The long non-coding RNA, HOX transcript antisense RNA (HOTAIR) was up-regulated in both pancreatic cancer tissues and cancer cell lines, and HOTAIR suppressed the expression of miR-613 via functioning as a competing endogenous RNA. In vivo studies showed that stable overexpression of miR-613 or knock-down of HOTAIR suppressed tumor growth and also reduced the expression of notch3. In conclusion, these results suggest that HOTAIR functions as a competing endogenous RNA to regulate notch3 expression via sponging miR-613 in pancreatic cancer.
INTRODUCTION
Pancreatic cancer is one of the most deadly cancers with a poor prognosis, and the overall 5-year survival rate is less than 5% due to advanced stage disease at initial diagnosis, frequent recurrence and the lack of effective therapies [1]. Surgical resection and chemotherapy were considered to be the main treatment options for pancreatic cancer [2]. Unfortunately, even with the advancement in medicine, pancreatic cancer is still resistant to current treatment regimens [3]. In this regard, it is important for us to understand the molecular mechanisms underlying pancreatic cancer progression.
Recently, the non-protein-coding portion of the genome is becoming more and more important in the regulation of normal physiology and the pathogenesis of diseases including cancer [4–6]. MicroRNA (miRNA) is a class of short non-coding RNAs, and it complements with the 3’ untranslated region (3′UTR) of mRNA and regulates the expressions of specific genes [7]. In the cancer studies, miRNAs can function to be oncogenic or tumor-suppressive. For example, miR-182 promotes pancreatic cancer cell proliferation and migration by targeting beta-TrCP2 [8]; miR-940 promotes cell proliferation by targeting GSK3β and sFRP1 in human pancreatic carcinoma [9]; on the other hand, miR-34a inhibits pancreatic cancer progression via Snail1-mediated epithelial-mesenchymal transition and the notch signaling pathway [10]; miR-377 inhibits the proliferation of pancreatic cancer by targeting Pim-3 [11]. The long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides, and they can not code proteins [12]. LncRNA was found to play important roles in various types of cancer including breast cancer, liver cancer, pancreatic cancer, colorectal cancer, gastric cancer and so on [12]. Various molecular mechanisms of lncRNA underlying cancer development have been proposed [13]. One of important mechanisms is that the lncRNA acts as a miRNA sponge to regulate the miRNA expression, which in turn regulates the expression of specific genes targeted by miRNA. Previously, our study has shown that the lncRNA, HOX transcript antisense RNA (HOTAIR) regulates the expression of miR-663b via the histone modification, which results in the regulation of pancreatic cancer progression [14]. In addition, HOTAIR was found to control the expression of Rab22a by sponging miR-373 in ovarian cancer [15]; Liu et al., showed that HOTAIR functions as a competing endogenous RNA to regulated human epidermal growth factor receptor 2 (HER2) expression by sponging miR-331-3p in gastric cancer [16]. Whether HOTAIR also functions as a competing endogenous RNA to regulate pancreatic cancer progression is largely unknown.
In the present study, we found that miR-613 was a potential target of HOTAIR by using bioinformatics tool. We then explored the role of miR-613 in pancreatic cancer. In vitro mechanistic studies revealed the tumor suppressive role of miR-613 via targeting neurogenic locus notch homolog protein 3 (notch3) in pancreatic cancer, and further study showed that HOTAIR functions as a competing endogenous RNA to regulate notch3 expression by sponging miR-613. In vivo and clinical results further confirmed the roles of miR-613 in pancreatic cancer progression. Therefore, our results may provide new insights into understanding the molecular mechanisms of miR-613 in pancreatic cancer.
RESULTS
The down-regulation of miR-613 in the pancreatic cancer tissues and cells lines
The miR-613 expression level in the clinical sample tissues from patients with pancreatic cancer was determined by the qRT-PCR assay, and the expression of miR-613 in the pancreatic cancer tissues was significantly lower than that in the adjacent normal pancreatic tissues (Figure 1A). The expression level of miR-613 in pancreatic cancer tissues was classified into low expression of miR-613 and high expression of miR-613 based on the median values of miR-613 expression in p. The association between miR-613 expression level and the clinicopathological parameters was analyzed in these 59 patients with pancreatic cancer, and low expression of miR-613 was significantly correlated with tumor differentiation, TNM stage and nodal metastasis, while miR-613 was not significantly associated with age, gender and tumor size (Table 1). We also followed the survival status of these patients, and the Kplan-Meier survival analysis showed that low expression of miR-613 was significantly correlated with shorter survival in patients with pancreatic cancer (Figure 1B) In addition, the expression levels in the pancreatic cancer cell lines including CFPAC-1, BXPC-3, L3.6pl and Panc-1 were also significantly down-regulated when compared to that in adjacent normal pancreatic tissues (Figure 1C).
Figure 1: MiR-613 was down-regulated in pancreatic cancer tissues and cancer cells, and was correlated with poor survival in patients with pancreatic cancer patients. (A) qRT-PCR analysis of miR-613 expression levels in adjacent normal pancreatic tissues and pancreatic cancer tissues from patients with pancreatic cancer. N = 59, significant difference between groups was shown as **P < 0.01 (Paired t-test). (B) Kplan-Meier survival analysis of patients with pancreatic cancer. Low expression level of miR-163 and high expression level of miR-613 was defined by using the median values of miR-613 expression levels in the 59 pancreatic cancer tissues. (C) qRT-PCR analysis of miR-613 expression levels in adjacent normal pancreatic cancer tissues and pancreatic cancer cell lines. N = 3, significant differences compared to adjacent normal pancreatic cancer tissues were shown as **P < 0.01, ***P < 0.001 (One-way ANOVA followed by Dunnett’s test).
Table 1: The association between miR-613 levels and clinicopathological characteristics of pancreatic cancer patients
miR-613 expression | |||
---|---|---|---|
Low expression | High expression | P-value | |
N = 28 | N = 31 | ||
Age | |||
< 56 years | 12 | 17 | 0.4379 |
≥ 56 years | 16 | 14 | |
Gender | |||
Male | 15 | 12 | 0.3015 |
Female | 13 | 19 | |
Tumor differentiation | |||
1–2 | 9 | 20 | 0.0191 |
3 | 19 | 11 | |
TNM stage | |||
I–II | 7 | 17 | 0.033 |
III/IV | 21 | 14 | |
Nodal metastasis | |||
0 | 7 | 19 | 0.0083 |
1 | 21 | 12 | |
Tumor size | |||
< 2 cm | 10 | 18 | 0.1188 |
> 2 cm | 18 | 13 |
Low expression and high expression of miR-613 was determined by the cut-off values for miR-613, which were defined as the cohort median, and median of age was used as the cut-off values to define the subgroup (< 56 years old and ≥ 56 years old group). Statistical significance between groups was analyzed by Chi-square tests.
Effect of miR-613 on pancreatic cell proliferation, cell invasion and migration, cell apoptosis and cell cycle
To further understand the molecular mechanisms of miR-613 in pancreatic cancer progression, we performed the gain-of-function studies. The overexpression of miR-613 in L3.6pl and Panc-1 cells was achieved by transient transfection with miR-613 mimics, and transfection with miR-613 mimics significantly increased the expression levels of miR-613 in L3.6pl and Panc-1 cells when compared to that transfected with scrambled miRNA (NC) (Figure 2A). The cell proliferative ability was measured by CCK-8 assay, and the cell proliferation was significantly suppressed in the pancreatic cells (L3.6pl and Panc-1) transfected with miR-613 mimics when compared to cells transfected with scrambled miRNA (Figure 2B). The cell growth was assessed by colony formation assay, and consistently, transfection with miR-613 mimics significantly suppressed the number of colonies formed by pancreatic cancer cells (L3.6pm and Panc-1) (Figure 2C). The cell invasion and migration were determined by Transwell invasion and Transwell migration assay, respectively, and the number of invaded and migrated cells in the miR-613 mimics-transfected group was significantly reduced when compared to the group transfected with scrambled miRNA (Figure 2D and 2E). The cell apoptosis and cell cycle were assessed by using flow cytometry experiments, and cell apoptosis was significantly increased by miR-613 transfection in the pancreatic cancer cells when compared to cells transfected with scrambled miRNA (Figure 2F). In the aspect of cell cycle, transfection with miR-613 mimics significantly increased the cell population in the G0/G1 phase accompanied by a decrease of cell population in the S phase in the pancreatic cancer cells (Figure 2G).
Figure 2: Up-regulation of miR-613 inhibited cell proliferation, invasion and migration in pancreatic cancer cells. (A) qRT-PCR analysis of miR-613 expression levels in pancreatic cancer cells (L3.6pl and Panc-1) after miR-613 mimics or scrambled miRNA transfection. NC = negative control (scrambled miRNAs), mimics = miR-613 mimics, n = 3, significant differences between groups were shown as ***P < 0.001 (Unpaired t-test). (B) Cell proliferation of pancreatic cells (L3.6pl and Panc-1) after miR-613 mimics transfection or scrambled miRNA was determined by CCK-8 assay. N = 3, significant differences compared to NC group were shown as *P < 0.05 (Two-way ANOVA followed by Bonferroni’s test). (C) Cell growth, (D) cell invasion, and (E) cell migration of pancreatic cells (L3.6pl and Panc-1) after miR-613 mimics or scrambled miRNA transfection were measured by colony formation assay, Transwell invasion assay, and Transwell migration assay, respectively. N = 3, significant differences compared to NC group were shown as **P < 0.01 (Unpaired t-test). (F) Cell apoptosis and (G) cell cycle of pancreatic cells (L3.6pl and Panc-1) after miR-613 mimics or scrambled miRNA transection was analyzed by flow cytometry. N = 3, significant differences compared to NC group were shown as *P < 0.05, **P < 0.01 (Unpaired t-test).
MiR-613 repressed notch3 expression via targeting its 3′TUR
As miRNAs exert their function by targeting the 3′UTR of specific genes, we predicted the potential targets of miR-613 by using the TargetScan tool (www.targetscan.org), and predicted results showed that notch3 was one of the potential targets of miR-613. To confirm that notch3 was a target of miR-613, we performed the luciferase reporter assay. As shown in Figure 3A, we constructed the luciferase reporter plasmids containing the wide type 3′UTR of notch3 or the mutated 3′UTR of notch3, and pancreatic cancer cells transfected with miR-613 mimics and reporter plasmids containing the wide type 3′UTR of notch3 significantly suppressed the luciferase activity when compared to cells transfected with scrambled miRNA and reporter plasmids containing the wild type 3′UTR of notch3 (Figure 3B). Co-transfection with miR-613 mimics and mutated 3′UTR of notch3 had no effect on the luciferase activity when compared to cells co-transfected with scrambled miRNA and mutated 3′UTR of notch3 (Figure 3B). In addition, transfection with miR-613 mimics significantly suppressed the mRNA and protein expression levels of notch3 when compared to cells transfected with scrambled miRNA (Figure 3C). To further confirm the functional role in pancreatic cancer progression, cells co-transfected with scrambled miRNA and notch3 overexpressing plasmids (pcDNA3.1-notch3) significantly suppressed cell proliferation when compared to the control group, while co-transfection miR-613 mimics and notch3 overexpressing plasmids restored the inhibitory effects of miR-613 mimics in pancreatic cancer cells (Figure 3D). The notch3 mRNA expression level in the clinical sample tissues in these patients were also examined, and the mRNA expression level of notch3 in pancreatic cancer tissues was significantly up-regulated when compared to that in adjacent normal pancreatic tissues (Figure 3E). The Pearson’s correlation analysis further showed that miR-613 expression level was inversely correlated with the expression level of notch3 mRNA in pancreatic cancer tissues (Figure 3F).
Figure 3: Notch3 was a downstream target of miR-613 in pancreatic cancer cells. (A) Putative miR-613-binding sequences in the 3′UTR of notch3 and the reporter constructs showing the wild type (WT) notch3 3′UTR sequence and the mutated (MUT) notch3 3′UTR sequence. (B) miR-613 mimics suppressed the luciferase activity of the WT but not the MUT 3′UTR of notch3 reporter in pancreatic cells (L3.6pl and Panc-1). NC = negative control (scramble miRNAs), mimics = miR-613 mimics, n = 3, significant differences compared to NC group were shown as *P < 0.05 (Unpaired t-test). (C) qRT-PCR and western blotting analysis of notch3 mRNA and protein expression levels in pancreatic cancer cells (L3.6pl and Panc-1) after miR-613 mimics or scrambled miRNA transfection. N = 3, significant differences compared to NC group were shown as **P < 0.01 (Unpaired t-test). (D) Cell proliferation of pancreatic cancer cells (P3.6lp and Panc-1) co-transfected with miR-613 mimics (or scrambled miRNAs) and pcDNA3.1-notch3 (or pcDNA3.1 vector) was determined by CCK-8 assay. Mimics = miR-163 mimcis, nothc3 = pcDNA3.1-notch3, n = 3, significant differences among groups were shown as *P < 0.05, **P < 0.01 (One-way ANOVA followed by Dunnett’s test). (E) qRT-PCR analysis of miR-613 levels in adjacent normal pancreatic tissues and pancreatic cancer tissues from patients with pancreatic cancer. N = 59, significant difference between groups was shown as **P < 0.01 (Paired t-test). (F) The inverse correlation between miR-613 and notch3 levels were analyzed by Pearson’s correlation test (R = –0.6918, P < 0.001).
Effects of HOTAIR on the miR-613 expression levels in pancreatic cancer cells
HOTAIR was found to play an important role in pancreatic cancer progression, which has been shown in our previous study [14], and we further examined whether miR-613 had an interaction with HOTAIR. The bioinformatics prediction using the starBase tool [17, 18] found that HOTAIR was also a target of miR-613, and HOTAIR may act as a sponge for miR-613. To confirm the interaction between miR-613 and HOTAIR, we performed the luciferase reporter assay. As shown in Figure 4A, we constructed the luciferase reporter plasmids containing the wide type HOTAIR or the mutated HOTAIR (Figure 4A), and pancreatic cancer cells transfected with miR-613 mimics and reporter plasmids containing the wide type HOTAIR significantly suppressed the luciferase activity when compared to cells transfected with scrambled miRNA and reporter plasmids containing the wild type HOTAIR in L3.6pl cells (Figure 4B). Co-transfection with miR-613 mimics and mutated HOTAIR had no effect on the luciferase activity when compare to cells co-transfected with scrambled miRNA and mutated HOTAIR L3.6pl (Figure 4B). In agreement with previous study, HOTAIR was up-regulated in pancreatic cancer cells when compared to that in adjacent normal pancreatic tissues (Figure 4C). Transfection with HOTAIR siRNA markedly suppressed the expression level of HOTAIR (Figure 4D), and also increased the expression level of miR-613 in pancreatic cancer cells when compared to cells transfected with scrambled miRNA (Figure 4E). The expression level of HOTAIR was further confirmed to be up-regulated in the pancreatic cancer tissues when compared to that in adjacent normal pancreatic tissues (Figure 4F). More importantly, the expression level of miR-613 was inversely correlated with the expression level of HOTAIR in the pancreatic cancer tissues (Figure 4G).
Figure 4: MiR-613 was inversely correlated with HOTAIR. (A) Putative miR-613-binding sequences of HOTAIR and the reporter constructs showing the wild type (WT) HOTAIR sequence and the mutated (MUT) HOTAIR sequence. (B) miR-613 mimics suppressed the luciferase activity of the WT but not the MUT HOTAIR reporter in pancreatic cells (L3.6pl). N = 3, significant difference between groups was shown as *P < 0.05 (Unpaired t-test). (C) qRT-PCR analysis of miR-613 expression levels in adjacent normal pancreatic cancer tissues and pancreatic cancer cell lines. N = 3, significant differences compared to adjacent normal pancreatic cancer tissues were shown as *P < 0.05, ***P < 0.001 (One-way ANOVA followed by Dunnett’s test). (D) qRT-PCR analysis of HOTAIR level in L3.6pl cells transfected with siHOTAIR or scrambled siRNA. N = 3, significant difference compared to NC group was shown as ***P < 0.001 (Unpaired t test). (E) qRT-PCR analysis of miR-613 level in L3.6pl cells transfected with siHOTAIR or scrambled lcnRNA. N = 3, significant difference compared to NC group was shown as ***P < 0.001 (Unpaired t test). (F) qRT-PCR analysis of HOTAIR levels in adjacent normal pancreatic tissues and pancreatic cancer tissues from patients with pancreatic cancer. N = 59, significant difference between groups was shown as **P < 0.01 (Paired t-test). (G) The inverse correlation between HOTAIR and miR663b levels were analyzed by Pearson’s correlation test (R = –0.5455, P < 0.001).
The effect of miR-613 and HOTAIR on the xenograft tumor growth in vivo
The role of miR-613 in pancreatic cancer was further examined in vivo by using the xenograft nude mice model. The nude mice were inoculated with L3.6pl cells transfected with LV-miR-613 or LV-Control, and in the LV-miR-613 transfected group, the tumor volume was significantly lower in LV-miR-613 group at 28, 35, and 42 d post inoculation than that in LV-Control group (Figure 5A). The mRNA and protein expression level of notch3 in the excised rumors were assessed by qRT-PCR and western blotting, respectively, and the mRNA and protein expression level of notch3 in the LV-miR-613 group were significantly lower than that in LV-Control group (Figure 5B), in addition, the protein levels of factors mediated apoptosis were also examined by western blotting, and the protein levels of caspase-3 was up-regulated and Bcl-2 was down-regulated in LV-miR-613 group (Figure 5B). We also inoculated the L3.6pl cells transfected with sh-HOTAIR or sh-Control, in the sh-HOTAIR group, the tumor volume was significantly lower at 28, 35, and 42 d post inoculation than that in the sh-Control group. In addition, the expression level of miR-613 was up-regulated in the excised tumors in the sh-HOTAIR group when compared to that in sh-Control group (Figure 5D). In addition, the mRNA and protein expression levels of notch3 were significantly lower in the tumors of sh-HOTAIR group than that in sh-Control group (Figure 5E).
Figure 5: Overexpression of miR-613 or knockdown of HOTAIR suppressed xenograft tumor growth in vivo. (A) Tumor volume changes in the mice bearing L3.6pl cells with miR-613 or scrambled miRNAs. N = 6, significant differences between groups were shown as *P < 0.05, **P < 0.01 (Two-way ANOVA followed by Bonferroni’s test). (B) qRT-PCR and western blotting analysis of notch3 mRNA and notch3, caspase-3 and Bcl-2 protein levels in the tumor tissues isolated from the mice bearing L36.pl cells with miR-613 or scrambled miRNAs. N = 6, significant difference between groups was shown as ***P < 0.05 (Unpaired t test). (C) Tumor volume changes in the mice bearing L3.6pl cells transfected with sh-HOTAIR or sh-Control. N = 6, significant differences between groups were shown as *P < 0.05, **P < 0.01 (Two-way ANOVA followed by Bonferroni’s test). (D) qRT-PCR analysis of miR-613 levels in the tumor tissues isolated from the mice bearing L3.6pl cells with sh-HOTAIR or sh-Control. N = 6, significant difference between groups was shown as ***P < 0.001 (Unpaired t test). (E) qRT-PCR and western blotting analysis of notch3 mRNA and protein levels in the tumor tissues isolated from the mice bearing L3.6pl cells with sh-HOTAIR or sh-Control. N = 6, significant difference between groups was shown as **P < 0.01 (Unpaired t test).
DISCUSSION
In the present study, our results showed that the expression of miR-613 was down-regulated in pancreatic cancer tissues and cancer cell lines. Down-regulation of miR-613 was positively correlated with tumor differentiation, advanced TNM stage, nodal metastasis and shorter overall survival in patients with pancreatic cancer. Overexpression of miR-613 inhibited cell proliferation, invasion and migration in pancreatic cancer cells. In addition, overexpression of miR-613 also induced apoptosis and cell cycle arrest at G0/G1 phase in the pancreatic cancer cells. Bioinformatics analysis, luciferase reporter assay and rescue experiment showed that notch3 was a novel a target of miR-613. Further bioinformatics analysis showed that HOTAIR functioned as a competing endogenous RNA to suppress the miR-613 expression. The effects of miR-613 on pancreatic cancer progression were further confirmed in the in vivo studies and in clinical samples.
Numerous studies have shown the important roles of miRNAs in the regulation of cancer progression. In pancreatic cancer, numerous miRNAs have been identified as either oncogenic miRNAs such as miR-191, miR-182, miR-940 [8, 9, 19] or tumor-suppressive miRNAs such as miR-34a and miR-454 [10, 20]. The role of miR-613 in pancreatic cancer has not been investigated. Previous studies have shown the tumor-suppressive role of miR-613 in various types of cancer including colorectal cancer, osteosarcoma, non-small cell lung cancer, papillary thyroid carcinoma, breast cancer, ovarian cancer, prostate cancer and esophageal squamous cell carcinoma [21–28]. Among these cancers, miR-613 was found to down-regulated in tumor tissues, and overexpression of miR-613 suppressed cell proliferation and cell invasion and migration. Consistently, our results for the first time showed that miR-613 was down-regulated in pancreatic cancer tissues and cell lines, and overexpression of miR-613 suppressed cell proliferation, cell invasion and migration. Further mechanistic study demonstrated that miR-613 overexpression induced cell apoptosis and caused cell cycle arrest at G0/G1 phase. Taken together, our results suggest that miR-613 plays a tumor-suppressive role in pancreatic cancer.
By using bioinformatics analysis, notch3 was predicted to be a downstream target of miR-613. Luciferase reporter assay and rescue experiment confirmed that miR-613 targets 3′UTR of notch3 to suppress the expression of notch3 and inhibits pancreatic cancer proliferation. Notch 3 is a protein that in encoded by the NOTCH3 gene in human and is one of the important mediators in the Notch signaling pathway, which has been shown to maintain a pool of pancreatic progenitor cells at the early stage of pancreatic development, and governs pancreatic ductal cell differentiation [29]. Notch3 has been shown to be overexpressed in several types of cancers, which makes it a promising therapeutic target for cancer treatment [30]. Previous studies showed that notch3 was significantly elevated in pancreatic cancer tissues compared to normal pancreatic tissues, and the elevation of notch3 was also associated with more advanced tumors and tumor metastasis [31, 32]. In addition, notch3 was significantly associated with reduced overall and disease-free survival in patients with pancreatic cancer [31], and high levels of notch3 was positively correlated with high invasion and overall survival rates in pancreatic ductal adenocarcinoma [32]. The sensitivity of tarextumab treatment in combination with gemcitabine in pancreatic tumors was associated with higher levels of notch3 expression [33]. Further study also showed the notch3 could enhance sensitivity to gemcitabine in pancreatic cancer via inactivation of PI3K/Akt-dependent pathway [34]. Consistently, our results also showed the up-regulation of notch3 in pancreatic cancer tissues compared to the adjacent normal pancreatic tissues, and overexpression of notch3 attenuated the inhibitory effect of miR-613 in pancreatic cancer cells, suggesting the oncogenic role of notch3 in the pancreatic cancer. As miR-613 may have other potential downstream targets, some of which have been demonstrated in other types of cancers. The present study limited the investigation into notch3, and further study may be performed to examine the role of other targets of miR-613 in pancreatic cancer.
Human HOTAIR, is a 2.2 kb lncRNA transcribed from the HOXC locus and has been shown to bind both to the PRC2 and the LSD1 complexes and aberrant expression of HOTAIR has been associated with tumorigenesis in various types of cancers including breast cancer, colorectal cancer, lung cancer, gastric cancer [35–38]. Our previous study showed that HOTAIR regulated the expression of miR-663b via histone modification in pancreatic cancer [14]. Recently, HOTAIR was proposed to function as a competing endogenous RNA by sponging miRNAs to regulate miRNAs levels. Zhang et al, demonstrated that HOTAIR controls the expression of Rab22a by sponging miR-373 in ovarian cancer [15]. HOTAIR modulated c-KIT expression through sponging miR-193a in acute myeloid leukemia [39]. In addition, HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-31-3p in gastric cancer [16]. Consistently, our data showed that HOTAIR was found to be up-regulated in both pancreatic cancer tissues and cell lines, and HOTAIR was inversely correlated with miR-613 level in pancreatic cancer tissues. Bioinformatics analysis and luciferase reporter assay showed that HOTAIR suppressed the expression miR-613 via sponging miR-613 in the pancreatic cancer cells. Knock-down of HOTAIR in L3.6pl cells suppressed the in vivo tumor growth and also suppressed the expression levels of miR-613, suggesting that the in vivo tumor suppressive role miR-613 may be correlated with HOTAIR. In summary, our results suggest the interaction between HOTAIR and miR-613 via “competing endogenous RNA” mechanism may be important for the pancreatic cancer progression.
In conclusion, the present study identified the interaction between HOTAIR and miR613 in pancreatic cancer. The results suggested that HOTAIR positively regulated the notch3 expression via acting as a competing endogenous RNA for miR-613 binding. The HOTAIR-miR-613-notch3 axis may be a promising therapeutic target for pancreatic cancer.
MATERIALS AND METHODS
Cell culture
The human pancreatic cancer cell lines (BXPC3, CFPAC-1, Panc-1 and L3.6pl) and HEK-293T cells were purchased from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). Cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Sigma, St Louis, USA) supplemented with 10% fetal bovine serum (HyClone, GE Healthcare Life Science, Logan, USA) and incubated in a humidified chamber supplemented with 5% CO2 at 37oC.
Tissue samples
Pancreatic cancer tissues and their matched adjacent normal pancreatic cancer tissues (approximately 5 cm from cancerous tissues) were taken from 59 patients undergoing surgery for pancreatic cancer at the First People’s Hospital of Changzhou. All cases were reviewed by pathologist and histologically confirmed as pancreatic cancer based on histopathological evaluation. The characteristics of the patients were shown in Table 1. All tissues were immediately snap-frozen in liquid nitrogen and stored at −80oC until further experimentation. After surgical resection treatments, patients were further followed up every 2–4 months. No local or systemic treatment was conducted in these patients before surgical operation. Informed consents were obtained from all patients, and the study was approved by the Research Ethics Committee of the First People’s Hospital of Changzhou.
Oligonucleotide transfection, plasmid construction and lentiviral infection
MiR-613 mimics and its negative control, scrambled miRNAs were purchased from Ribobio (Guangzhou, China). The notch3 mRNA sequences were synthesized and subcloned into the pCDNA3.1 vector, and the empty pcDNA3.1 vector served as negative control (Genepharma, Shanghai, China). For HOTAIR, the siRNA specially targeting HOTAIR or its non-target negative control siRNA was synthesized by Genepharma. For transfection, L3.6pl cells or Panc-1 cells were grown on six-well plates to 60% confluence, and miRNA, siRNA or plasmid transfection was performed by using Lipofectamine 2000 (Invitrogen) according to the manufacturer’s protocol. For the rescue experiment, cells were co-transfected with miRNAs (miR-613 mimics or scrambled miRNA) and plasmids (pcDNA3.1-notch3 or pcNDA3.1). Total RNA and protein were extracted at 24 h post transfection and used for qRT-PCR and western blot analysis.
MiR-613-overexpressing lentiviral constructs were generated using synthetic oligonucleotides and the Lv-CMV-GPF vector (Genepharma), and sh-HOTAIR-overexpressing lentiviral constructs were generated by subcloning sh-HOTAIR into pGLV3/H1/GFP lentiviral frame plasmids (Genepharma); empty vectors were used as negative controls, respectively. All the constructed plasmids were confirmed by sequencing (Invitrogen, Carlsbad, USA). Lentivirus carrying miR-613 or HOTAIR was packaged in the HEK293T cells and collected from the supernatants following the manufacturer’s protocol. Stable cell lines for xenograft study were established by infecting the lentivirus into L36.pl cells.
RNA extraction and qRT-PCR analysis
Total RNA was extracted from cells or tissues using TRIzol reagent (Invitrogen) according to the manufacturer’s instruction. Notch3 mRNA level and HOTAIR was quantified by qRT-PCR using a SYBR Premix ExTaq Reverse Transcription PCR kit (Takaka, Dalian, China) and GAPDH was used as an internal control for normalization. For miR-613 detection, the miR-613 level was quantified by qRT-PCR using TaqMan assay kits (Applied Biosystems, Forster City, USA) and U6 was used an internal control for normalization. The primers for qRT-PCR were shown in Supplementary Table 1. The reaction was performed by using an ABI PRISM 7500 Sequence Detection System (Applied Biosystems). The relative expression levels of notch3 mRNA, HOTAIR and miR-613 were calculated using 2-∆∆Ct method.
CCK-8 assay
Cell proliferation assay was performed with Cell Counting Kit-8 (Dojindo, Kumamoto, Japan) according to the manufacturer’s instruction. Twenty-four hour after transfection, cells were seed in 96-well plates at 1 × 104 cells per well. The proliferative ability of cells was determined at 0, 24, 48, and 72 h. For the rescue experiment, cells were co-transfected with miRNAs (miR-613 mimics or scramble miRNA) and plasmids (pcDNA3.1-notch3 or pcNDA3.1), and the proliferative ability of cells was determined at 48 h. CCK-8 (10 μl) was added to each well at different time points, and the plate was incubated for 2 h at room temperature. The absorbance was measured at 450 nm using a microplate spectrophotometer (Molecular Devices, Sunnyvale, USA).
Colony formation assay
Twenty-four hours after transfection, cells were seeded for colony formation in 6-cm culture dishes at a density of 1000 cells per dish. After continuous culture for 14 days, cells were fixed in methanol for 10 min and stained with 0.1% crystal violet for 30 min. Visible colonies were manually counted.
Transwell invasion and migration assay
Transwell assay was performed using a chamber of 6.5 mm in diameter and with an 8-mm pore size (Corning, Corning, USA). For the invasion assay, twenty-four hours after transfection, 5 × 104 cells in DMEM supplemented with 0.1% FBS were seeded onto the Matrigel-coated membrane matrix of the upper chamber, and the lower chamber was filled with DMEM supplemented with 10% FBS. After 24 h incubation, cells invading the bottom of the membrane were stained with 0.1% crystal violet in 20% ethanol, and the number of invaded cells was counted by using a DM2500 bright field microscope (LEICA, Wetzlar, Germany). For the migration assay, twenty-four hours after transfection, 5 × 104 cells in DMEM supplemented with 0.1% FBS were seeded in the chamber of 6.5 mm of the upper chamber in diameter and with an 8-mm pore size (Corning), and the lower chamber was filled with DMEM supplemented with 10% FBS. After 24 h incubation, cells invading the bottom of the membrane were stained with 0.1% crystal violet in 20% ethanol, and the number of invaded cells was counted by using a DM2500 bright field microscope (LEICA).
Cell apoptosis and cell cycle analysis
Twenty-four hour after transfection, cells were trypsinized and fixed with 70% ethanol for 30 min on ice. RNA was degraded by incubation with 20 mg/ml RNase (Sigma) for 1 h at 37ºC. For the cell cycle analysis, cells were stained with propidium iodide at room temperature for 30 min, and then were analyzed by Calibur flow cytometry (BD Biosciences, Franklin Lakers, USA). For cell apoptosis analysis, cells were stained with FITC-Annexin V and propidum iodide (Beyotime, Beijing, China) and then were analyzed by Calibur flow cytometry (BD Biosciences) equipped with CellQuest software (BD Biosciences).
Luciferase reporter assay
Cells were seeded in 96-well plates at 1 × 104 cells per well. When the cells reached 60% confluence, they were co-transfected with wide type pGL3-notch3 3′UTR, mutated pGL3-notch3 3′UTR plasmids, wild type pGL3-HOTAIR, or mutated pGL3-HOTAIR and either scramble miRNA or miR-613 mimics using Lipofectamine 2000 (Invitrogen). Forty-eight hours after transfection, luciferase activity was measured with the Dual-Luciferase Reporter Assay System (Promega, Madison, USA) and expressed as the ratio between firefly and Renilla luciferase activities.
Western blotting
Protein were extracted from cells or tissues using RIPA lysis buffer. Proteins were then separated by SDS-PAGE. After electrophoresis, proteins were transferred onto polyvinylidene difluoride membrane. After blocking with 5% non-fat milk for 2 h at room temperature, the membranes were then incubated with rabbit polyclonal notch3 antibody (ab23426; Abcam, Cambridge, UK), rabbit polyclonal caspase-3 (#9662; Cell Signaling Technology, Beverly, USA), mouse monoclonal Bcl-2 (#15071; Cell Signaling Technology) or mouse monoclonal β-actin (sc-47778; Santa Cruz) overnight at 4°C. The membranes were then incubated 2 h at room temperature with horseradish peroxidase-conjugated goat anti-rabbit or goat anti-mouse (Santa Cruz, Dalla, USA) secondary antibody and visualized with a chemiluminescence kit (Pierce, Rockford, USA).
In vivo animal study
Tumor formation was studied by establishing a xenograft nude mice model. Four-week-old BALB/c nude mice were purchased from the Shanghai Experimental Animal Center (Chinese Academy of Sciences). The animal experiments in this study were approved by the Animal Research Committee of the First People’s Hospital of Changzhou. Care and handling of the animals were in accordance with the guidelines for Institutional and Animal Care and Use Committees. A total of 24 animals were randomly divided into 4 groups with 6 animals in each group. Infected L3.6pl cells (1 × 106/mice) were subcutaneously injected into the neck area of the nude mice. Tumor volumes (mm3) were measured every 7 days and calculated using the following formula: volume = width × length × height/2. The animals were killed 42 days after cell inoculation, and tumor tissues were harvested for further analysis.
Statistical analysis
All the data were expressed as mean ± SD. Statistical analysis were performed by using Student’s t-test or ANOVA followed by multiple comparison tests. The relationship between the expression of miR-613, notch3 and HOTAIR was examined by Pearson’s correlation analysis. The correlation between miR-613 and pathological parameters was determined by Chi-square test. Differences were considered to be statistically significant when P < 0.05. All the results were performed in at least three independent experiments.
ACKNOWLEDGMENTS
We thank for the grant support from General Research Program of Jiangsu Provincial Commission of Health and Family Planning (Project No.H201548).
CONFLICTS OF INTEREST
There is no conflicts of interest.
FUNDING
The project was supported by grants from General Research Program of Jiangsu Provincial Commission of Health and Family Planning (Project No.H201548).
REFERENCES
1. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2017. CA Cancer J Clin. 2017; 67:7–30.
2. Fogel EL, Shahda S, Sandrasegaran K, DeWitt J, Easler JJ, Agarwal DM, Eagleson M, Zyromski NJ, House MG, Ellsworth S, El Hajj I, O’Neil BH, Nakeeb A, et al. A Multidisciplinary Approach to Pancreas Cancer in 2016: A Review. Am J Gastroenterol. 2017; 112:537–554.
3. Coveler AL, Herman JM, Simeone DM, Chiorean EG. Localized Pancreatic Cancer: Multidisciplinary Management. Am Soc Clin Oncol Educ Book. 2016; 35:e217–226.
4. An X, Sarmiento C, Tan T, Zhu H. Regulation of multidrug resistance by microRNAs in anti-cancer therapy. Acta Pharm Sin B. 2017; 7:38–51.
5. Cui J, Zhou B, Ross SA, Zempleni J. Nutrition, microRNAs, and Human Health. Adv Nutr. 2017; 8:105–112.
6. Trager MM, Dhayat SA. Epigenetics of Epithelial-to-Mesenchymal-Transition in pancreatic carcinoma. Int J Cancer. 2017 Jan 30. [Epub ahead of print].
7. Farazi TA, Hoell JI, Morozov P, Tuschl T. MicroRNAs in human cancer. Adv Exp Med Biol. 2013; 774:1–20.
8. Wang S, Ji J, Song J, Li X, Han S, Lian W, Cao C, Zhang X, Li M. MicroRNA-182 promotes pancreatic cancer cell proliferation and migration by targeting beta-TrCP2. Acta Biochim Biophys Sin (Shanghai). 2016; 48:1085–1093.
9. Yang HW, Liu GH, Liu YQ, Zhao HC, Yang Z, Zhao CL, Zhang XF, Ye H. Over-expression of microRNA-940 promotes cell proliferation by targeting GSK3beta and sFRP1 in human pancreatic carcinoma. Biomedicine & pharmacotherapy. 2016; 83:593–601.
10. Tang Y, Tang Y, Cheng YS. miR-34a inhibits pancreatic cancer progression through Snail1-mediated epithelial-mesenchymal transition and the Notch signaling pathway. Sci Rep. 2017; 7:38232.
11. Chang W, Liu M, Xu J, Fu H, Zhou B, Yuan T, Chen P. MiR-377 inhibits the proliferation of pancreatic cancer by targeting Pim-3. Tumour Biol. 2016; 37:14813–14824.
12. Jiang C, Li X, Zhao H, Liu H. Long non-coding RNAs: potential new biomarkers for predicting tumor invasion and metastasis. Mol Cancer. 2016; 15:62.
13. Bartonicek N, Maag JL, Dinger ME. Long noncoding RNAs in cancer: mechanisms of action and technological advancements. Mol Cancer. 2016; 15:43.
14. Cai H, An Y, Chen X, Sun D, Chen T, Peng Y, Zhu F, Jiang Y, He X. Epigenetic inhibition of miR-663b by long non-coding RNA HOTAIR promotes pancreatic cancer cell proliferation via up-regulation of insulin-like growth factor 2. Oncotarget. 2016; 7:86857–86870. doi: 10.18632/oncotarget.13490.
15. Zhang Z, Cheng J, Wu Y, Qiu J, Sun Y, Tong X. LncRNA HOTAIR controls the expression of Rab22a by sponging miR-373 in ovarian cancer. Mol Med Rep. 2016; 14:2465–2472.
16. Liu XH, Sun M, Nie FQ, Ge YB, Zhang EB, Yin DD, Kong R, Xia R, Lu KH, Li JH, De W, Wang KM, Wang ZX. Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer. Mol Cancer. 2014; 13:92.
17. Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014; 42:D92–97.
18. Yang JH, Li JH, Shao P, Zhou H, Chen YQ, Qu LH. starBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data. Nucleic Acids Res. 2011; 39:D202–209.
19. Liu H, Xu XF, Zhao Y, Tang MC, Zhou YQ, Lu J, Gao FH. MicroRNA-191 promotes pancreatic cancer progression by targeting USP10. Tumour Biol. 2014; 35:12157–12163.
20. Fan Y, Shi C, Li T, Kuang T. microRNA-454 shows anti-angiogenic and anti-metastatic activity in pancreatic ductal adenocarcinoma by targeting LRP6. Am J Cancer Res. 2017; 7:139–147.
21. Chen X, Xu Y, Liao X, Liao R, Zhang L, Niu K, Li T, Li D, Chen Z, Duan Y, Sun J. Plasma miRNAs in predicting radiosensitivity in non-small cell lung cancer. Tumour Biol. 2016; 37:11927–11936.
22. Guan S, Wang C, Chen X, Liu B, Tan B, Liu F, Wang D, Han L, Wang L, Huang X, Wang J, Yao B, Shi J, et al. MiR-613: a novel diagnostic and prognostic biomarker for patients with esophageal squamous cell carcinoma. Tumour Biol. 2016; 37:4383–4391.
23. Li B, Xie Z, Li Z, Chen S, Li B. MicroRNA-613 targets FMNL2 and suppresses progression of colorectal cancer. Am J Transl Res. 2016; 8:5475–5484.
24. Li X, Sun X, Wu J, Li Z. MicroRNA-613 suppresses proliferation, migration and invasion of osteosarcoma by targeting c-MET. Am J Cancer Res. 2016; 6:2869–2879.
25. Qiu W, Yang Z, Fan Y, Zheng Q. MicroRNA-613 inhibits cell growth, migration and invasion of papillary thyroid carcinoma by regulating SphK2. Oncotarget. 2016; 7:39907–39915. doi: 10.18632/oncotarget.9530.
26. Ren W, Li C, Duan W, Du S, Yang F, Zhou J, Xing J. MicroRNA-613 represses prostate cancer cell proliferation and invasion through targeting Frizzled7. Biochem Biophys Res Commun. 2016; 469:633–638.
27. Wu J, Yuan P, Mao Q, Lu P, Xie T, Yang H, Wang C. miR-613 inhibits proliferation and invasion of breast cancer cell via VEGFA. Biochemical and biophysical research communications. 2016; 478:274–278.
28. Zhang X, Zhang H. Diminished miR-613 expression as a novel prognostic biomarker for human ovarian cancer. European review for medical and pharmacological sciences. 2016; 20:837–841.
29. Kim W, Shin YK, Kim BJ, Egan JM. Notch signaling in pancreatic endocrine cell and diabetes. Biochemical and biophysical research communications. 2010; 392:247–251.
30. Purow B. Notch inhibition as a promising new approach to cancer therapy. Advances in experimental medicine and biology. 2012; 727:305–319.
31. Mann CD, Bastianpillai C, Neal CP, Masood MM, Jones DJ, Teichert F, Singh R, Karpova E, Berry DP, Manson MM. Notch3 and HEY-1 as prognostic biomarkers in pancreatic adenocarcinoma. PloS one. 2012; 7:e51119.
32. Zhou JX, Zhou L, Li QJ, Feng W, Wang PM, Li EF, Gong WJ, Kou MW, Gou WT, Yang YL. Association between high levels of Notch3 expression and high invasion and poor overall survival rates in pancreatic ductal adenocarcinoma. Oncology reports. 2016; 36:2893–2901.
33. Yen WC, Fischer MM, Axelrod F, Bond C, Cain J, Cancilla B, Henner WR, Meisner R, Sato A, Shah J, Tang T, Wallace B, Wang M, et al. Targeting Notch signaling with a Notch2/Notch3 antagonist (tarextumab) inhibits tumor growth and decreases tumor-initiating cell frequency. Clinical cancer research. 2015; 21:2084–2095.
34. Yao J, Qian C. Inhibition of Notch3 enhances sensitivity to gemcitabine in pancreatic cancer through an inactivation of PI3K/Akt-dependent pathway. Medical oncology (Northwood, London, England). 2010; 27:1017–1022.
35. He H, Wei Z, Du F, Meng C, Zheng D, Lai Y, Yao H, Zhou H, Wang N, Luo XG, Ma W, Zhang TC. Transcription of HOTAIR is regulated by RhoC-MRTF-A-SRF signaling pathway in human breast cancer cells. Cellular signalling. 2017; 31:87–95.
36. Kogo R, Shimamura T, Mimori K, Kawahara K, Imoto S, Sudo T, Tanaka F, Shibata K, Suzuki A, Komune S, Miyano S, Mori M. Long noncoding RNA HOTAIR regulates polycomb-dependent chromatin modification and is associated with poor prognosis in colorectal cancers. Cancer research. 2011; 71:6320–6326.
37. Liu MY, Li XQ, Gao TH, Cui Y, Ma N, Zhou Y, Zhang GJ. Elevated HOTAIR expression associated with cisplatin resistance in non-small cell lung cancer patients. Journal of thoracic disease. 2016; 8:3314–3322.
38. Yan J, Dang Y, Liu S, Zhang Y, Zhang G. LncRNA HOTAIR promotes cisplatin resistance in gastric cancer by targeting miR-126 to activate the PI3K/AKT/MRP1 genes. Tumour biology. 2016 Nov 30. [Epub ahead of print].
39. Xing CY, Hu XQ, Xie FY, Yu ZJ, Li HY, Bin Z, Wu JB, Tang LY, Gao SM. Long non-coding RNA HOTAIR modulates c-KIT expression through sponging miR-193a in acute myeloid leukemia. FEBS letters. 2015; 589:1981–1987.