Introduction
A large number of people in the world suffer from gastric cancer (GC), which is characterized by high mortality [1]. As reported, approximately 723,000 deaths occurred in 2012 and this figure accounted for 8.8% of all cancer-related deaths in 2012 [2]. The incidence of GC varies from region to region. Eastern Asia exhibits the highest incidence of GC whereas Northern Africa exhibits the lowest incidence of GC [1]. Several exogenous risk factors for GC have been confirmed: dietary habits, smoking, alcohol consumption, obesity, aging, viruses and bacteria such as Helicobacter pylori [3]. Besides environmental and behavioral factors, a variety of tumor-related genes have also been reported in GC, including both oncogenes and tumor suppressive genes [4]. Surgery is so far the most effective approach for GC and it is usually carried out with chemotherapy or radiotherapy [5]. The 5-year survival rate decreases from 60% to 2% when tumor recurrence, invasion or metastasis occurs [5]. Since at least 15% of GC patients experience tumor metastasis and surgical procedures may fail to improve the survival status of GC patients, studies are conducted to discover alternative therapies that are able to ameliorate this issue [6, 7].
MicroRNAs (miRNAs) are endogenous small non-coding RNA molecules with a length of 18–24 nucleotides and they modulate gene expression as posttranscriptional regulators [8, 9]. MicroRNAs regulate gene expression by controlling mRNA stability, degradation, and gene translation. The discovery of miRNAs was traced back to the early 1990s [10]; however, they were not recognized as biological regulators until the early 2000s [11]. Since many miRNAs are deregulated in various types of cancers, researchers suspected that miRNAs may act as tumor suppressors or oncogenes [12–14]. Several miRNAs are aberrantly expressed in GC tissues [15, 16] and they are likely to be used as biomarkers for GC screening and diagnosis [17]. For example, microRNA-429 functions as a tumor suppressor of GC cells by targeting FSCN1 [18] while microRNA-16a functions as an oncogene in GC, stimulating proliferation and metastasis in vitro and in vivo [19]. Therefore, ongoing research in this area is worthwhile.
MiR-135b-5p is a microRNA which has been suggested to be associated with cancer cell proliferation and migration since overexpression of miR-135b-5p was often discovered in cancer cells [20]. It has also been reported that miR-135b-5p can directly target several genes, including MEF2C [21] and DISC1 [22]. Krüppel-like factor 4 (KLF4) is predicted as an another target of miR-135b-5p and it is characterized as a pluripotency factor [23]. For example, KLF4 is able to regulate cell cycle progression and differentiation, which are important for the corresponding life cycles of some cancer viruses such as HPV [24]. KLF4 is suppressed by miR-10b and downregulation of KLF4 facilitates the proliferation and metastases of osteosarcoma cells [25]. However, in a subset of melanoma, ectopic expression of KLF4 promotes melanoma cell growth and reduces cell apoptosis, suggesting that depletion of KLF4 could help to control melanoma growth [26]. In our research, we confirmed the low expression of KLF4 in GC tissues compared with adjacent tissues. Hence, we attempted to uncover the regulatory relationship between miR-135b-5p and KLF4, which could provide a profound understanding of the corresponding mechanism of GC carcinogenesis.
In this study, sufficient evidence indicated that ectopic miR-135b-5p negatively regulated KLF4 expression, which also contributed to the proliferation, migration and invasion of GC cells. Therefore, downregulation of miR-135b-5p could be a novel strategy for GC clinical treatment.
Material and methods
Patients and tissue specimens
Sixty (33 male and 27 female patients, average age 45) gastric cancer specimens were obtained from patients who were hospitalized in the Department of Gastrointestinal Colorectal and Anal Surgery of China-Japan Union Hospital of Jilin University. Written consent forms were signed by all the patients and this study was approved by the hospital ethics committee (IRB number: CJH5413-16-05-15). All patients were first diagnosed with gastric cancer and treated without radiotherapy or chemotherapy. Pathology was graded in accordance with World Health Organization criteria. All specimens consisted of gastric cancer tissues and adjacent non-tumorous tissues (distance to GC tissues > 5 cm).
Cell culture
Human gastric carcinoma cell lines BGC-823, MKN-28, SGC-7901 and the human normal gastric epithelial cell line (GES1) (BNBio, Beijing, China) were cultured in RPMI 1640 medium with 10% FBS, 100 U/ml penicillin and 100 U/ml streptomycin under standard culture conditions, 2–3 passages.
Cell transfection
MiR-135b-5p mimics and KLF4 cDNA were synthesized by Molbase biological reagent company. After amplification by PCR, KLF4 cDNA was inserted into the pcDNA3.1 plasmid. BGC-823 cells were first seeded in six-well plates and when cell growth reached 70–80% confluence, transfection was performed with the help of Lipofectamine 2000 reagents (Invitrogen, Carlsbad, CA, USA). Cells were harvested at 48 h after transfection. Transfection efficiency was confirmed using RT-qPCR and Western blot.
QRT-PCR
Total RNA was isolated with Trizol reagent (Thermo Fisher Scientific, MA, USA) and was reverse-transcribed into cDNA before RT-qPCR analysis. The miRNA and mRNA expression levels were normalized against U6 and GAPDH respectively and quantified using the 2–ΔΔCt method. The primers are listed in Table I.
Western blots
Proteins were extracted using RIPA lysis (Beyotime, Shanghai, China) and the protein concentrations were quantified with BCA Kit (Sigma-Aldrich, St. Louis, MO, USA). Proteins were transferred from SDS-PAGE gel to PVDF membranes, which were then blocked with 5% nonfat milk and subsequently incubated with KLF4 and GAPDH antibodies (Abcam, Cambridge, MA, USA) at 4°C overnight. Then the membranes were incubated with secondary antibodies (Abcam) for 1 h at room temperature, followed by analysis with enhanced chemiluminescence (ECL) Western blotting detection reagents (GE Healthcare, Little Chalfont, Buckinghamshire, UK).
Dual-luciferase reporter assay
According to TargetScan (http://www.targetscan.org/), KLF4 is one of miR-135b-5p’s targets. Sequences containing the miR-124-3p target region of EDNRB 3′UTR and its mutation (Sangon, Shanghai, China) were inserted into pmirGLO reporting vector (Promega, USA). Then, pmirGLO recombinant vector was co-transfected with miR-124-3p mimics, inhibitor and control oligo in HEK293T cells. Firefly and Renilla luciferase activities were detected consecutively at 48 h after co-transfection.
MTS assay
BGC-823 cells of different groups were seeded in a 96-well plate at a density of 1 × 103 per well. The absorptions of cells was measured at various time point (0, 12, 24, 36, 48, 60, 72 h) using a Cell Titer 96 AQueous One Solution Cell Proliferation Assay (MTS) (Promega, Beijing, China). The cell absorbance at OD 570 nm was measured after staining.
Colony formation assay
Transfected BGC-823 cells were seeded in six-well plates (500 cells per well) and incubated in a standard environment for 2 weeks. Then colonies were fixed with 4% paraformaldehyde and stained with 0.1% crystal violet. Five random fields were photographed and the number of colonies was counted.
Transwell assay
Transwell assays were performed using 8 μm-pore transwell chambers (Corning company) coated with Matrigel (BD Bioscience). 500 μl of medium containing 10% FBS (the attractant) was added into lower chambers, and 100 μl serum-free medium into upper chambers. The transfected BGC-823 cells were plated in the upper chambers. After 48-hour incubation, the cells on the bottom surface of the upper chambers were cleaned with a cotton swab. Invading cells on another surface were fixed with 4% paraformaldehyde and stained with 0.1% crystal violet for 30 min. Five random fields under a microscope were photographed and then the cell numbers were counted.
Wound healing assay
At 48 h after transfection, 8 × 105 cells were cultured into 35 mm2 culture dishes and 200 μl sterile pipette tips were used to scratch the culture surface. The wound closure was measured to determine the migration of cells. Photographs were taken at 0 and 24 h. Each experiment was repeated 3 times.
Flow cytometry
For cell cycle analysis, transfected cells were stained with PI using the Cycletest PLUS DNA Reagent Kit (BD Biosciences). The FITC Annexin V Apoptosis Detection Kit (BD Biosciences, Bedford, MA, USA) was used for apoptosis analysis according to the protocol of the manufacturer. Afterwards, a flow cytometer (BD FACS Aria III, San Jose, CA, USA) was applied to count cells in different status.
Statistical analysis
The results were presented as mean ± SD. Normally distributed measurement data were compared with the t test or analysis of variance (ANOVA); otherwise data were compared with the nonparametric rank-sum test. A p-value less than 0.01 or 0.05 was considered to indicate a statistically significant difference.
Results
Expression levels of miR-135b-5p and KLF4 in GC tissues and para-carcinoma normal tissues, GC cells and normal cells
We analyzed the expression level of miR-135b-5p and KLF4 in GC tissues (n = 60) and corresponding para-carcinoma normal tissues using RT-qPCR. The results indicated that miR-135b-5p expression levels in GC tissues were significantly higher than that in para-carcinoma tissues (p < 0.05, Figure 1 A), whereas the KLF4 expression levels in GC tissues were significantly lower than in para-carcinoma tissues (p < 0.05, Figure 1 B). Meanwhile, the results of RT-qPCR showed that miR-135b-5p expression levels in BGC-823, MKN-28 and SGC-7901 GC cell lines were significantly higher compared with that in the GES1 cell line (p < 0.05, Figure 1 C), and miR-135b-5p expression levels in BGC-823 were the most significant. Both RT-qPCR results and western blot results indicated that the KLF4 expression in GC cell lines was remarkably suppressed (Figures 1 D, E).
MiR-135b-5p directly targets KLF4 3′UTR
To further investigate the association between miR-135b-5p and KLF4, we first found the potential miR-135b-5p targeting site in the KLF4 3′UTR region using the TargetScan database (Figure 2 A). Wild-type and mutated 3′UTR were synthesized and harbored in the pmirGLO reporting plasmid. Co-transfection was performed in HEK 293T cells and dual luciferase activity was detected at 48 h after transfection. Compared with the control group, luciferase activities of cells co-transfected with miR-135b-5p mimics and KLF4 wild-type 3′UTR-pmirGLO plasmid were remarkably decreased in HEK293T cells (p < 0.05), whereas luciferase activities of cells co-transfected with miR-135b-5p mimics and mutated KLF4 3′UTR pmirGLO plasmid showed no significant differences compared with the control group (Figure 2 B), which demonstrated that KLF4 3′UTR was the target of miR-135b-5p.
MiR-135b-5p inhibits expression of KLF4
In order to investigate the regulatory relationship between miR-135b-5p and KLF4, we first manipulated the expression of miR-135b-5p and KLF4 in BGC-823 cells via transfection. The cells were divided into the following five groups: the blank control group treated with no reagents, the negative control group transfected with pcDNA3.1 vector and miR mimics control, the pcDNA3.1 + mimics group transfected with pcDNA3.1 vector and miR-135b-5p mimics, the oligo + KLF4 group transfected with miR mimics control and pcDNA3.1-KLF4 cDNA, and the mimics + KLF4 group transfected with miR-135b-5p mimics and pcDNA3.1-KLF4 cDNA. QRT-PCR and western blot were applied to determine the expression of miR-135b-5p and KLF4 and confirm the transfection efficiency. MiR-135b-5p expression significantly increased in the pcDNA3.1 + mimics group and mimics + KLF4 group (p < 0.05, Figure 3 A). There was no significant difference in miR-135b-5p expression among the blank control, negative control and oligo + KLF4 group, indicating that KLF4 has no effect on miR-135b-5p expression. However, both mRNA expression and protein level of KLF4 were repressed in the pcDNA3.1 + mimics group and overexpressed in the oligo + KLF4 group compared to the blank control group (p < 0.05, Figures 3 B, C). Intriguingly, the upregulation of KLF4 by pcDNA3.1-KLF4 cDNA was attenuated by co-transfection of miR-135b-5p in BGC-823 cells. Collective data indicated that miR-135b-5p could directly target KLF4 and suppress its expression.
MiR-135b-5p modulates biological functions of BGC-823 cells via inhibiting KLF4
Since the regulatory relationship between miR-135b-5p and KLF4 has been confirmed, the effects of the miR-135b-5p/KLF4 axis on gastric cancer cell behaviors which may explain the aberrant expression of miR-135b-5p and KLF4 in gastric cancer tissues and cells drew our attention. Hence, after manipulating miR-135b-5p and KLF4 expression, we accessed the proliferation, migration, invasion, cell cycle and cell apoptosis of transfected cells.
MTS assays suggested that there were no significant differences among blank control, negative control and mimics + KLF4 group. In the 24–72 h period, BGC-823 cells viability significantly increased in the miR-135b-5p high-expressed group and decreased in the KLF4 overexpressed group (p < 0.05, Figure 4 A). Expectedly, overexpression of miR-135b-5p and KLF4 in the mimics + KLF4 group showed a result of antagonism. Similar results were obtained via colony formation assay. At 48 h after transfection, the colony number of different groups was counted (Figure 4 B). Overexpression of miR-135b-5p facilitated colony formation compared to the blank control group while upregulation of KLF4 efficiently repressed the colony number (p < 0.05). The repression was impaired by co-transfection of miR-135b-5p and there was no significant difference in colony number among blank control, negative control and mimics + KLF4 group. All those observations suggested that miR-135b-5p promotes the proliferation of BGC-823 cells via repressing KLF4.
Transwell assays and wound healing assays indicated that there were no significant differences between the blank control group and the negative control group (p > 0.05, Figures 4 C, D). Compared with the blank control group, the number of invading cells of the miR-135b-5p high-expressed group increased significantly (p < 0.05), which meant cell invasion capacity was enhanced; the wound closure in pcDNA3.1 + mimics was also larger, suggesting an enhanced cell migratory capacity. However, the number of invading cells of the KLF4 high-expressed group decreased significantly, and the wound closure was smaller (p < 0.05). Still, the weakened cell invading and migratory capacity could be repaired by co-transfection of miR-135b-5p mimics.
To further study the changes of cell cycle and cell apoptosis rate, flow cytometry and specialized staining were applied. Compared with the control group, the proportion of BGC-823 cells in G0/G1 phase decreased significantly, the proportion of cells in S and G2/Mphases increased significantly, and the cell apoptosis rate decreased significantly in the miR-135b-5p mimics group (p < 0.05, Figures 5 A, B). In the oligo + KLF4 group, the cell proportion of G0/G1 phase increased significantly, the cell proportion of S and G2/Mphases decreased significantly, and the apoptosis rate increased (p < 0.05). Hence, miR-135b-5p could efficiently promote cell mitosis and suppress cell apoptosis while KLF4 plays an opposite role. Then in the mimics + KLF4 group, the inhibitory effects of KLF4 were attenuated by miR-135b-5p. Taken together, the evidence proved that miR-135b-5p promotes cell proliferation, migration, invasion, and mitosis, and reduces the cell apoptosis rate of BGC-823 cells via downregulating KLF4.
Discussion
Accumulating evidence has proven that many miRNAs are altered in gastric cancer cells and lead to the deregulation of various genes, indicating their roles as either oncogenes or tumor suppressors. For instance, miR-524-5p suppressed the growth and invasion of GC cells [27]. MiR-216b was found downregulated in human GC cells, which could inhibit proliferation and cell cycle progression by targeting HDACB [28]. Also, miR-429 functioned as a tumor suppressor in GC cells, downregulating the oncogene FSCN1 [29]. On the other hand, there were studies revealing that miR-532-5p could promote cell growth, migration and invasion as an oncogenic miRNA by directly targeting RUNX3 [30]. In addition, miR-935 was markedly upregulated in GC tissues and cells, and facilitated cell proliferation through targeting SOX7 [31]. These studies showed that the dysregulation of miRNAs may play different roles in GC, working together to regulate the metabolism in cancer cells. In the present study, we found that miR-135b-5p was overexpressed in GC cells, compared with the GES1 cell line. Our experiment results first demonstrated that the overexpression of miR-135b-5p could significantly enhance the proliferation, migration, invasion and mitosis of GC cells. Apoptosis was also suppressed by miR-135b-5p, indicating a tumor-promoting role of miR-135b-5p in gastric cancer cells.
MiR-135b-5p has been found to modulate gene expression in various cancers. For instance, by targeting MEF2C, miR-135b-5p could promote cell proliferation and migration in atherosclerosis [32]; Wang et al. reported that miR-135b-5p could regulate the level of DISC1 mRNA, and they also found that DISC1 regulation was allele specific, i.e. miR-135b-5p only bound to the major allele (A), not to the minor allele (G) in a neuropsychiatric disorder [33]. However, a limited number of studies have been carried out to investigate the interaction of miR-135b and its target genes in gastric cancer progression. Only those by Gao et al. and Liu et al. demonstrated that miR-135b and miR-29c were promising tissue-specific biomarkers of gastric carcinogenesis [34, 35]. Tseng et al. found that miR-135b-5p was upregulated and could be modulated by interference with β-catenin in gastric cancer [36]. The interplay between miR-135b-5p and KLF4 possibly played an important role in gastric cancer development, as suggested by our study described herein.
KLF4 is a critical cellular transcription factor which is able to regulate cell cycle progression and differentiation [37]. It is also known as the gut-enriched Kruppel-like factor, engaging in gastrointestinal epithelial cell proliferation and differentiation, and is found downregulated in gastrointestinal cancers [38]. KLF4 could also promote ghrelin expression by binding to the ghrelin promoter [39]. Althubiti et al. have suggested that the downregulation of KLF4 was related to diffuse-type gastric cancer [40]. Furthermore, it was shown that KLF4 negatively regulates epithelial-mesenchymal transition (EMT) in GI cells, indicating its role in cancer cells as a tumor suppressor gene [41]. Dopamine was also shown to hinder the proliferation of GC cells by upregulating KLF4, which is an inhibitor of the cell cycle, by reducing AKT phosphorylation [42]. Thus, the dysregulation of KLF4 is assumed to be closely related to GC pathogenesis. To be specific, the low level of KLF4 is speculated to lead to the progression of gastric cancers. In line with former studies, we observed low-level expression of KLF4 in GC cells, compared with GES1 cell lines, and replenishment of KLF4 could significantly inhibit the proliferation, migration and invasion and induce G0/G1 arrest and cell apoptosis of BGC-823 cells, indicating that KLF4 had a tumor suppressive effect in GC cells. However, the limitations of this study should be taken into consideration. For example, KLF4 is not the only target of miR-135b-5p, and similarly miR-135b-5p is not the only miRNA which regulates KLF4. Thus because of the complicated regulation network, our study just contributed to the comprehensiveness of the gastric cancer mechanism.
In conclusion, we have proved that miR-135b-5p is overexpressed in GC cells while KLF4 expression is inhibited, supplementing the previous studies. Furthermore, the overexpression of miR-135b-5p can promote the proliferation, migration as well as invasion of GC cells, suggesting that miR-135b-5p may serve as an onco-miRNA. More importantly, we have demonstrated that there is a negative relationship between miR-135b-5p and KLF4, which has never been studied before, showing the possibility that the knockout of miR-135b-5p in GC tumor cells may serve as a therapeutic strategy by overexpressing anti-oncogenes such as KLF4. However, further in vivo studies need to be carried out to investigate the influence of the miR-135b-5p/KLF4 axis on gastric cancer.