- Open Access
miR-148b-3p functions as a tumor suppressor in GISTs by directly targeting KIT
Cell Communication and Signalingvolume 16, Article number: 16 (2018)
Gain-of-function mutations and overexpression of KIT are characteristic features of gastrointestinal stromal tumor (GIST). Dysregulation in miRNA expression may lead to KIT overexpression and tumorigenesis.
miRNA microarray analysis and real-time PCR were used to determine the miRNA expression profiles in a cohort of 69 clinical samples including 50 CD117IHC+/KITmutation GISTs and 19 CD117IHC−/wild-type GISTs. GO enrichment and KEGG pathway analyses were performed to reveal the predicted targets of the dysregulated miRNAs. Of the dysregulated miRNAs whose expression was inversely correlated with that of KIT miRNAs were predicted by bioinformatics analysis and confirmed by luciferase reporter assay. Cell counting kit-8 (CCK-8) and flow cytometry were used to measure the cell proliferation, cycle arrest and apoptosis. Wound healing and transwell assays were used to evaluate migration and invasion. A xenograft BALB/c nude mouse model was applied to investigate the tumorigenesis in vivo. Western blot and qRT-PCR were used to investigate the protein and mRNA levels of KIT and its downstream effectors including ERK, AKT and STAT3.
Of the six miRNAs whose expression was inversely correlated with that of KIT, we found that miR-148b-3p was significantly downregulated in the CD117IHC+/KITmutation GIST cohort. This miRNA was subsequently found to inhibit proliferation, migration and invasion of GIST882 cells. Mechanistically, miR-148b-3p was shown to regulate KIT expression through directly binding to the 3’-UTR of the KIT mRNA. Restoration of miR-148b-3p expression in GIST882 cells led to reduced expression of KIT and the downstream effectors proteins ERK, AKT and STAT3. However, overexpression of KIT reversed the inhibitory effect of miR-148b-3p on cell proliferation, migration and invasion. Furthermore, we found that reduced miR-148b-3p expression correlated with poor overall survival (OS) and disease-free survival (DFS) in GIST patients.
miR-148b-3p functions as an important regulator of KIT expression and a potential prognostic biomarker for GISTs.
GIST is the most common mesenchymal tumors of the digestive tract, representing 1–3% of gastrointestinal malignancies . Activating mutations in the receptor tyrosine kinase KIT or platelet-derived growth factor receptor alpha (PDGFRA) are the principal oncogenic triggers of GIST, which have been identified in ~ 80% and 5–10% of GISTs, respectively . Approximately 10% of GIST cases do not contain KIT and PDGFRA mutations, which are called ‘wild-type’ (wt) GISTs [3, 4]. Recent studies have revealed that wt GISTs represent a heterogeneous group in clinical, pathological and biological features that is profoundly different from KIT/PDGFRA mutant tumors . Moreover, the wt GISTs often exhibit a greater degree of resistance to the tyrosine kinase inhibitor imatinib .
MicroRNAs (miRNAs) comprise a class of 18–25 nucleotides, non-coding RNAs that specifically bind to the 3′-untranslated region (3’-UTR) of target mRNAs [7, 8]. Growing evidence indicates that miRNAs are involved in multiple cancer-relevant processes, including migration, invasion and angiogenesis [9,10,11]. To date, several studies have explored the function of miRNAs in the tumorigenesis of GISTs, demonstrating the involvement of miRNA dysregulation in gene mutation, chromosomal changes, tumor risk and metastasis [12,13,14,15]. Since activating mutation in the KIT gene is the main oncogenic trigger in > 80% of GISTs, the aberrant expression of miRNAs that interrelates with KIT-induced tumorigenesis is of particular interest .
With this in mind, we profiled for miRNA expression in a cohort of 5 CD117IHC+/KITmutation and 5 CD117IHC−/wild-type GISTs and validated our findings in another cohort of 59 samples comprising 45 CD117IHC+/KITmutant and 14 CD117IHC−/wild-type GISTs. By this approach, we identified 6 miRNAs that were dysregulated and whose expression was inversely related to KIT mutation or KIT overexpression. In particular, miR-148b-3p emerged as a significantly downregulated miRNA that we subsequently showed to suppress proliferation, migration and invasion of GIST882 cells by directly targeting the 3′-UTR of KIT at nucleotides 1639–1656. Restoration of miR-148b-3p expression in GIST882 cells led to reduced expression of not only KIT but also its downstream effectors, including ERK, AKT and STAT3. Overexpression of KIT reversed the suppressive actions of miR-148b-3p on cell proliferation, migration and invasion. Moreover, the reduction of miR-148b-3p was correlated with high risk grade, high mitotic rate and the presence of recurrence or metastasis in GIST patients. Kaplan-Meier analysis indicated that GIST patients with low levels of miR-148b-3p correlated with poor overall survival (OS) and disease-free survival (DFS) rates. These results indicate that miR-148b-3p plays an essential role in promoting the progression of GISTs by reducing KIT expression. Furthermore, our work has identified miR-148b-3p as a potential prognostic biomarker for GIST and a candidate therapeutic target.
Patients and tissue samples
Formalin-fixed paraffin-embedded (FFPE) samples of 115 primary GISTs, including 80 cases of CD117IHC+ GISTs and 35 cases of CD117IHC- GISTs, were selected from over 1000 GIST patients who received surgery at Tongji Hospital between January 2007 and December 2012 (Table 1). No patient received receptor tyrosine kinase inhibitor therapy before the surgery. All were assessed for clinical and histomorphological features, and the risk grade was evaluated on the basis of the National Institute of Health (NIH) consensus classification criteria. Ethical approval was obtained from the Ethics committee of Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology. Each participant provided written informed consent to join this study.
All cases had been confirmed by immunohistochemical and KIT and PDGFRA mutational analyses before diagnosis. Immunohistochemical staining was performed on FFPE specimens using EnVision systems (Dako). The QIAamp DNA FFPE Tissue Kit (56,404; Qiagen; Hilden, Germany) was applied to isolate genomic DNA from FFPE tissues. KIT exons 9, 11, 13, 17 and PDGFRA exons 12, 18 were amplified with specific primers (Table 2). Additionally, for the 24 cases of GISTs lacking mutations, immunohistochemical analysis was carried out to measure the expression of succinate dehydrogenase complex subunits B (SDHB) in order to detect SDH-deficiency GIST. The 19 cases of CD117IHC−/wild-type GISTs employed for microarray and qRT-PCR analysis were all positive for DOG1 and SDHB.
The miRNA microarray profiling was conducted by LC Sciences (Houston, TX) according to Sanger miRBase Version 19.0. Hybridization was carried out on a μParaflo microfluidic chip (Atactic Technologies, Houston, TX, USA). Fluorescent signals were scanned by a microarray scanner (GenePix 4000B, Molecular Device, USA) and analyzed using the GenePix Pro 6.0 software. The Locally Weighted Regression by LOWESS filter was used to normalize the data. The P-value was corrected by the false discovery rate (FDR). Only dysregulated miRNAs that fulfilled the benchmarks of p value < 0.05 and fold change > 1.5 were chosen out for further investigation. The microarray data have been submitted to the Gene Expression Omnibus (GEO) and the GEO accession number is GSE73346.
Antibodies and reagents
The antibodies against SDHB (10620–1-AP, Proteintech, Philadelphia, PA, USA), ERK (9102, Cell Signaling, San Jose, CA, USA), phospho-ERK (4376, Cell Signaling), AKT (4691, Cell Signaling), phospho-AKT (4060, Cell Signaling), STAT3 (sc-482, Santa Cruz, Dallas, TX, USA), phospho-STAT3 (sc-8059, Santa Cruz), β-actin (sc-47,778, Santa Cruz), KIT (SAB4501647, Sigma-Aldrich, St Louis, MO, USA), HRP (horseradish peroxidase)-labeled anti-rabbit IgG (A0208, Beyotime Biotechnology, Shanghai, China), and HRP-labeled anti-mouse IgG (A0216, Beyotime Biotechnology) were used at the appropriate concentration.
Oligonucleotides, plasmid construction and cell transfection
miR-148b-3p mimics, miR-148b-3p inhibitor, mimic control and inhibitor control were obtained RiboBio Co. (Guangzhou, Guangdong, China). KIT-coding sequences without the 3′-UTR were cloned into the GV141 vector by GeneChem (Shanghai, China) (called pGV141-KIT). The psiCHECK-2 vector (C8021) was purchased from Promega (Madison, WI, USA). The psiCHECK-2-KIT-wt, psiCHECK-2-KIT-mut-1378, and psiCHECK-2-KIT-mut-1639 were constructed by Bios Biological (Wuhan, Hubei, China). miR-148b-3p agomir and miR agomir NC were generated by GenePharma Co., Ltd. (Shanghai, China). Lipofectamine 3000 (L3000015, Invitrogen, Grand Island, NY, USA) was used to transfect cells with the desired genes or plasmids.
Total RNA was reverse-transcribed by RevertAid cDNA Synthesis Kit (K1622, Thermo Fisher Scientific, Waltham, MA, USA) and amplified on an ABI Prism 7900HT platform (Applied Biosystems, Foster City, CA, USA) with specific miRNA primers (Table 2). U6 snRNA or GAPDH was used as endogenous control. The comparative Ct method was employed to estimate the relative expression levels.
Target prediction, pathway and microRNA gene network analysis
The prediction of miRNA targets and the effector pathways were revealed by Gene Ontology (GO) enrichment  and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses . GO and KEGG terms with corrected P-value and FDR below 0.01 were defined as dramatically enriched in target candidate genes and pathways. By bioinformatics algorithms TargetScan and miRanda, we obtained the relationships between microRNAs and the predicted target genes from the gene databases. Interaction pattern between miRNAs and target genes were created using CytoScape .
The GIST882 cell line with a KIT homozygous missense mutation (exon 13, K642E) was kindly provided by Dr. Jonathan Fletcher, Dana-Farber Cancer Institute, Boston, MA, USA. GIST882 cells were cultivated in RPMI-1640 (GIBCO, Carlsbad, CA, USA) containing 10% FBS (GIBCO) and cultured in a humidified incubator with 5% CO2 at 37 °C.
Cell viability and proliferation assay
The Cell Counting Kit-8 kit (CCK-8; CK04–100, Dojindo, Kumamoto Prefecture, Kyushu, Japan) was used to evaluate cell viability. GIST882 cells (5 × 103 cells/well) with or without imatinib treatment were placed into a 96-well plate and treated with 10 μl CCK8 solution for 3 h. Absorbance at 450 nm was detected by an ELx800 Absorbance Reader (BioTek Instruments, Inc., Winooski, VT, USA).
Cell cycle and apoptosis analyses
The propidium iodide cell cycle analysis kit (C1052, Beyotime Institute of Biotechnology) was used to assess the cell cycle arrest of GIST882 cells on a FACScan flow cytometer (BD Biosciences, USA). The Annexin V-FITC/PI Apoptosis Detection Kit (556,547, BD Biosciences) was used to detect the apoptosis in GIST882 cells by flow cytometry.
Cell migration and invasion assays
GIST882 cells suspended in RPMI-1640 medium containing 2% FBS were plated in the upper chamber to assess the migration ability, whereas cells suspended in serum-free RPMI-1640 were added into transwell inserts (3422, Corning, Bedford, MA, USA) precoated with Matrigel (356,234, Corning) to detect the invasion ability. The number of cells crossed over the membrane were fixed in 4% paraformaldehyde, dyed with 0.1% crystal violet and calculated under an inverted microscope after 24 h of incubation.
Wound healing assay
Linear wounds were produced in the cell monolayer by a sterile pipette tip after cells cultures reached a confluence of 90–95%. The cells migrated into the wounded area or protruded from the border of the wound were quantified under an inverted microscope after 24 h of incubation in serum-free medium.
Dual-luciferase reporter assay
The miR-148b-3p mimics and vectors (psiCHECK-2-KIT-wt, psiCHECK-2-KIT-mut-1378, or psiCHECK-2-KIT-mut-1639) were co-transfected into the GSIST882 cells via Lipofectamine 3000. The dual-luciferase reporter assay system (E1910, Promega) were conducted to assess the firefly and Renilla luciferase activities. Renilla luciferase was considered as an internal control.
Protein samples were prepared using RIPA lysis buffer (R0278, Sigma-Aldrich), phosphatase inhibitor cocktail tablets (04906845001, Roche) and protease inhibitor cocktails tablets (04693132001, Roche). Total protein was resolved on a 12% SDS-PAGE gel and transferred onto a nitrocellulose membrane. After blocking with 5% BSA, the membranes were blotted with primary antibodies overnight and secondary antibodies for 1 h. The SuperSignal West Femto Maximum Sensitivity Substrate (34,095, Thermo Fisher Scientific, Waltham, MA, USA) with HRP Substrate (34,095, Thermo Fisher Scientific) was employed to visualize the blots.
Xenografts in nude mice
Athymic BALB/c nude mice (female, 4–6 weeks old) were provided by Huafukang Biotechnology Co. Ltd. (Beijing, China). All animal studies were approved by the Animal.
Ethics Committee of Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology. Approximately 5 × 106 GIST882 cells suspended in culture medium were subcutaneously injected into the flanks of mice. Tumor growth was assessed biweekly and estimated using the formula: 1/2 (length × width2). When the tumor volume reached 0.5–0.6 cm3, miR-148b-3p agomir or miR agomir NC was syringe-injected into the tumor mass at the dose of 50 μg per mouse every week. After 3 weeks, the mice were euthanized and tumor tissues were made into paraffin sections for further study.
The statistical analysis was implemented by the SPSS version 19.0 software. Student’s t-test was utilized to compare the variables between two groups, while one-way ANOVA was utilized to compare the variables among multiple groups. Fisher’s exact test and Chi-square test were applied to sort the GO and KEGG pathway categories. Kaplan-Meier method was performed to estimate OS and DFS and the log-rank test was carried out to assess the statistical difference between groups. Data were presented as mean ± SD from three independent experiments. P < 0.05 was considered statistically significant.
MicroRNA profiling of GIST tissue
The microarray data revealed distinct expression patterns for 9 mRNAs that were significantly different between the CD117IHC+/KITmutation GIST and CD117IHC−/wild-type GIST patient groups (fold change> 1.5; P < 0.05). Among these miRNAs, miR-122-3p, miR-483-3p, miR-101-3p, miR-598 and miR-4492 were significantly upregulated, whereas miR-140-5p, miR-148b-3p, miR-1587 and miR-4507 were significantly downregulated in the CD117IHC+/KITmutation GISTs. In complementary experiments, we used real-time PCR to determine changes in expression for 7 of these miRNAs between the two patient groups (Fig. 1b). Except for one, six miRNAs exhibited the same trend in both the microarray and RT-PCR analysis (Fig. 1c). These included four downregulated miRNAs (miR-140-5p, miR-148b-3p, miR-1587 and miR-4507) and two upregulated miRNAs (miR-483-3p and miR-598).
Target prediction and miRNA-mRNA interaction network analysis
TargetScan, miRBase and miRanda algorithms were employed to predict the target genes for the above six miRNAs (Table 3). GO enrichment and KEGG pathway analyses were subsequently performed to evaluate the functions of the target genes. The target genes were analyzed by GO enrichment analysis in molecular function, cellular component and biological process (Fig. 2a). KEGG pathway analysis revealed that 16,353 target genes were involved in 108 biological functions (Fig. 2b). A number of oncogenes, including KIT, PDGFRA, PIK3CA, AKT2, FGFR1, SMAD3, PTEN and EGFR emerged from both the GO and KEGG analysis of the miRNA-regulated gene network (Fig. 3). Because these genes have been shown to played essential roles in GIST tumorigenesis [19,20,21], our miRNA-mRNA regulatory network analysis revealed several candidate for future investigation of miRNA-target functions that control the development and progression of GIST.
Identification of KIT as a direct target of miR-148b-3p through integrative bioinformatics analysis
Because of the pivotal role of KIT in GIST pathogenensis, we next focused our study on identifying the miRNA that regulates KIT expression in GIST. To this end, we performed a computational target prediction analysis by the TargetScan and miRanda bioinformatic algorithms in combination to search for miRNAs that may bind the 3′UTR of the KIT transcript. Of the six candidate miRNAs, we found that only miR-148b-3p could directly target KIT mRNA by interacting with the 3′-UTR of KIT at nucleotides 1378–1393 and 1639–1656 (Fig. 7a).
miR-148b-3p suppressed the proliferation of GIST cells
The biological significance of miR-148b-3p was assessed by transfecting the miR-148b-3p mimics, miR-148b-3p inhibitor, or the corresponding negative control into the GIST882 cells. The effects of the transfection were verified by qRT-PCR, which indicated that the mimic and inhibitor oligos affected miR-148b-3p expression in GIST882 cells as predicted (Fig. 4a). As shown by the CCK-8 assay, upregulation of miR-148b-3p by the mimics suppressed the proliferation of the GIST882 cells whereas downregulation of miR-148b-3p by the inhibitor oligo had an opposite effect (Fig. 4b). Intriguingly, neither the mimics nor the inhibitor had any effect on cell cycle arrest and apoptosis (Fig. 4c-d). These results indicate that miR-148b-3p has an anti-proliferative effect on the GIST882 cells.
miR-148b-3p augmented the effect of imatinib in suppressing GIST cells proliferation
To investigate whether miR-148b-3p sensitized GIST882 cells to imatinib, we firstly detected the IC50 concentration of imatinib using cell viability assay and confirmed that the IC50 values for imatinib in GIST882 cells was 1.654 μM for 24 h (Fig. 5a). We then evaluated the effect of miR-148b-3p in combination with a 24 h treatment of imatinib (5 μM) on proliferation, cycle arrest and apoptosis in GIST882 cells. As the data shown, imatinib treatment repressed cell proliferation and induced G1-phase arrest and apoptosis in GIST882 cells compared to untreated cells (Fig. 5b-d). Overexpression of miR-148b-3p in combination with imatinib treatment more notably inhibited cell proliferation in GIST882 cells compared to the cells treated with imatinib alone (Fig. 5b). However, miR-148b-3p overexpression did not affect cell cycle arrest and apoptosis in GIST882 cells either alone or in combination with imatinib treatment (Fig. 5c-d). These results indicated that miR-148b-3p acts synergistically with imatinib to suppress cell proliferation and sensitizes GIST882 cells to imatinib treatment.
miR-148b-3p inhibited the migration and invasion of GIST cells
Microscopic observations on wound healing demonstrated that miR-148b-3p overexpression dramatically reduced the speed of wound closure in GIST882 cells (P < 0.01; Fig. 6c). In contrast, the wound closure speed was improved after miR-148b-3p inhibitor transfection (P < 0.05; Fig. 6c). Similarly, transwell migration and invasion assays showed that the cells penetrating the transwell membrane were markedly reduced in the miR-148b-3p mimics group (P < 0.05 and P < 0.01; Fig. 6a-b). In contrast, the miR-148b-3p inhibitor group exhibited high rates of mobility and invasiveness (P < 0.01; Fig. 6a-b). Collectively, miR-148b-3p upregulation markedly reduced the migration and invasion ability of GIST882 cells. This observation may have relevance to the progression and metastasis of GIST.
KIT is a direct target of miR-148b-3p in GIST
To establish the miRNA-mRNA relations between miR-148b-3p and KIT, a dual-luciferase detection system was performed to measure the transcriptional activity of the KIT constructs after psiCHECK-2-KIT-wt or psiCHECK-2-KIT-mut transfection. As shown in Fig. 7b, miR-148b-3p mimics remarkably attenuated the relative luciferase activity of psiCHECK-2-KIT-wt and psiCHECK-2-KIT-mut1378 compared to the negative control. However, the miR-148b-3p mimics did not affect the relative luciferase activity of the psiCHECK-2-KIT-mut1639, suggesting that the sequence between nucleotides 1639–1656 is the major binding site of miR-148b-3p. Western blotting showed that transfection with miR-148b-3p mimics result in an obvious reduction of KIT protein, whereas that with miR-148b-3p inhibitor led to an increase in KIT in GIST882 cells (Fig. 7e). However, qRT-PCR indicated that miR-148b-3p alterations did not change the expression of endogenous KIT mRNA in GIST882 cells (Fig. 7c). Taken together, miR-148b-3p could directly target KIT mRNA in GIST882 cells by interacting with the 3′-UTR of KIT, and negatively regulated KIT protein expression by inhibiting the translation of KIT at the post-transcriptional level rather than mRNA degradation as they did not change KIT mRNA expression level.
Downregulation of KIT leads to the inhibition of ERK, AKT and STAT3
To test whether the downregulation of KIT can lead to the inhibition of known signaling mediators downstream of KIT, we treated the GIST882 cells with 5 μM imatinib mesylate (Novartis) for 24 h and determined the protein levels of ERK, AKT and STAT3 by Western blot. As expected, the inhibition of KIT led to a rapid inhibition of phosphorylated ERK, AKT and STAT3 protein in GIST882 cells (Fig. 7d). These results implied that ERK, AKT and STAT3 activation is coupled to KIT activation in the GIST882 cells.
KIT signaling is perturbed by miR-148b-3p mimics or inhibitor transfection
Transfection with miR-148b-3p mimics dramatically reduced phosphorylated and total ERK, AKT and STAT3 protein of GIST882 cells, while miR-148b-3p inhibition have the opposite effect (Fig. 7e). These results suggest that the oncogenic KIT signal transduction was modulated by miR-148b-3p in GISTs.
Overexpression of KIT reversed the suppressor effects of miR-148b-3p on proliferation, migration and invasion of GIST cells
To elucidate whether the suppressive effects of miR-148b-3p on GISTs were induced by KIT repression, KIT expression was restored in GIST882 cells by transfecting the cells with the pGV141-KIT vector after miR-148b-3p mimics transfection. The ectopic overexpression of KIT not only rescued the cell growth inhibition induced by miR-148b-3p but also reversed the suppressor effectors of miR-148b-3p on the migration and invasion of GIST882 cells (Fig. 8a-d). Moreover, miR-148b-3p decreased the expression of phosphorylated and total ERK, AKT and STAT3 in GIST882 cells, all of which could be restored by KIT overexpression (Fig. 8e). These results suggest that KIT was directly responsible for modulating the biological effects of GIST cells induced by miR-148b-3p.
miR-148b-3p suppresses tumor growth in vivo
Consistent with its effects in vitro, miR-148b-3p remarkably suppressed the growth of xenografted tumors in mice. The tumors treated with miR-148b-3p agomir had a lower volume and weight than the tumors treated with miR agomir NC (Fig. 9b-c). Morphological similarities between the xenografts and human GISTs were evident based on stained paraffin-embedded sections (Fig. 9d). Immunohistochemical staining of the xenograft tumor sections indicated that the subcutaneous tumors of the miR-148b-3p-agomir-treated nude mice exhibited reduced signal for KIT compared to the control tumors (Fig. 9e). In parallel, Western blotting demonstrated that the tumors from the miR-148b-3p-agomir-treated mice exhibited a lower protein level of KIT and a significant reduction in phosphorylated and total ERK, AKT and STAT3 (Fig. 9f). These results indicated that exogenous miR-148b-3p suppressed tumor growth in vivo through KIT.
Downregulation of miR-148b-3p correlates with high risk grades and poor survival for GIST patients
We identified statistically significant correlations between miR-148b-3p expression and several clinicopathological parameters, including tumor risk grade, mitotic rate and recurrence or metastasis (Table 3). Interestingly, the miR-148b-3p expression of all GIST samples as well as CD117IHC+/KITmutation GISTs were negatively associated with tumor risk grade in a step-wise manner (Fig. 10a). Moreover, there was a notable trend for lower miR-148b-3p expression among GIST patients with recurrence or metastasis (Fig. 10b). We also examined the relationship between miR-148b-3p expression and overall survival (OS) and disease-free survival (DFS) rates in GIST patients by Kaplan-Meier analysis with log-rank tests, and found that GIST patients with lower expression of miR-148b-3p have poorer OS rates and higher tumor recurrence rates (P < 0.05, Fig. 10c-d).
miRNAs comprise a class of small noncoding RNAs that implicated to undertake prominent roles in multiple tumor-related processes, such as oncogenesis, progression, invasion, and metastasis [9,10,11]. miRNA dysregulation has been established as one of the key characteristics contributing to carcinogenesis . Recently, numerous studies have explored the association between miRNA expression profiles and GIST pathogenesis, which implicated that the dysregulation of miRNAs is associated with anatomical localization, mutation status, chromosome 14q loss, tumor risk, and biological behavior of GISTs [12,13,14,15, 22]. As gain-of-function mutations of KIT and KIT overexpression are characteristic features of GISTs , the aberrant expression of miRNAs that interrelates with KIT-induced tumorigenesis of GISTs is of particular interest [13, 22, 23].
Although the presence of KIT mutations usually leads to strong KIT expression in GISTs, the mechanism of KIT overexpression is not yet well understood . Dysregulation of miRNAs expression might relevant with KIT mutations and KIT overexpression and the resulting tumorigenesis in GISTs. Recently, several studies have indicated that miRNAs can serve as critical modulators of the KIT expression, and these mechanisms are essential for the initiation and progression of GISTs. For example, Koelz et al. found the repressed expression of miR-221 and miR-222 was inversely relevant with KIT expression and might be considered as a putative molecular mechanism in the pathogenesis of GISTs [25,26,27]. miR-494 has been identified to downregulate KIT by directly binding to two separate seed sequences on the KIT mRNA, thereby inhibiting proliferation, migration and invasion of GIST882 cells . Another microRNA, miR-218, was recently found to diminish KIT protein and suppress cell proliferation and invasion in GIST-T1 cell lines . However, only a few miRNAs involved in tumorigenesis that are related to mutation of KIT have been identified in GISTs. Haller et al. assessed the relevance between miRNAs expression signatures and the mutation status in twelve GIST samples and discovered that 16 miRNAs were performanced in a mutation-dependent manner . Bachet et al. also found that 13 mRNAs were specifically associated with the KIT homozygous/heterozygous mutation status in GISTs . Thus, investigations on the function of miRNAs that are specifically associated with KIT mutation and KIT overexpression are required to improve the current knowledge of GIST and may offer novel insights into GIST therapy.
In this study, we evaluated the miRNA expression profiles in a cohort of CD117IHC+/KITmutation GISTs and CD117IHC−/wild-type GISTs and confirmed 6 aberrant miRNAs whose dysregulated expression was inversely related with KIT mutations and KIT overexpression. GO enrichment and KEGG pathway analyses were carried out to detect the role of abnormally expressed miRNAs in GIST pathogenesis, and the results showed that these miRNAs were obviously enriched in tumor-associated signaling pathways, including the TGF-beta signaling pathway, the MAPK signaling pathway and the Wnt signaling pathway. Among the validated miRNAs, miR-148b-3p and miR-140-5p were predicted to target KIT, PDGFRA, PIK3CA, AKT2, PTEN or EGFR genes, which are important intermediate molecules of the KIT signaling pathway as shown by the miRNA gene network analysis [19, 31]. Based on the structural analysis of the six candidate miRNA and the 3′-UTR of KIT miRNA, we found that only miR-148b-3p could directly target KIT mRNA by interacting with the 3′-UTR of KIT at two conserved putative binding sites. According to the above results, miR-148b-3p is postulated to be a tumor suppressor involved in GIST tumorigenesis by directly binding to the 3′-UTR of KIT.
Accumulating evidence has demonstrated that miR-148b-3p downregulation is a frequent event in multiple tumors, such as pancreatic cancer, breast cancer, colorectal cancer, gastric cancer and lung cancer [32,33,34,35,36,37,38]. Functional studies also revealed that miR-148b-3p could act as a tumor suppressor via directly targeting distinct oncogenes such as DNMT3b , CCK2R , DNMT1 , CEA  and AMPKα1  in these cancers. Moreover, miR-148b-3p downregulation is responsible for the inhibition of cell proliferation, progression and invasion in gastric cancer, pancreatic cancer and NSCLC [35,36,37,38], as well as for distant metastasis or poor prognosis of breast cancer and pancreatic cancer [32, 38]. Although the underlying molecular mechanism of miR-148b-3p involved tumorigenesis has been identified in several studies, the probable role of miR-148b-3p in GISTs remained to be elucidated. Therefore, we detected the molecular mechanism of miR-148b-3p in GISTs by gain- and loss-of-function analysis. Our data revealed that miR-148b-3p upregulation dramatically suppressed proliferation, migration and invasion of GIST882 cells as well as reduced the tumorigenic ability of GIST882 cells in nude mice. Besides,miR-148b-3p interacted synergistically with imatinib in the inhibition of cell proliferation and sensitizes GIST882 cells to imatinib treatment. Clinical analysis also demonstrated that the low expression of miR-148b-3p in tumor specimens of GISTs was correlated with poor prognostic factors, including high risk grade, high mitotic rate and the presence of recurrence or metastasis. These results indicated that miR-148b-3p downregulation might contribute to the progression of GISTs, and thus miR-148b-3p could be recommended as a promising prognostic biomarker in GISTs.
Although previous bioinformatics analyses have indicated KIT as a direct target of miR-148b-3p, the dual-luciferase reporter assay was further carried out and confirmed that only miR-148b-3p could directly target KIT by binding to the potential 3′-UTR at nucleotides 1639–1656. Transfection of miR-148b-3p mimics or inhibitor was found to alter the KIT protein level without affecting its mRNA level, which suggested that miR-148b-3p negatively regulates KIT by a post-transcriptional mechanism. Furthermore, western blotting confirmed that miR-148b-3p downregulated KIT expression but that this could be restored by transfection of the cells with KIT plasmid. Using CCK-8 and transwell assays, we observed that restoration of KIT expression not only rescued the cell growth inhibition induced by miR-148b-3p but also reversed the migration and invasion suppressor functions of miR-148b-3p in GIST882 cells. Taken together, these results imply that miR-148b-3p inhibits cell proliferation, migration and invasion of GIST via directly binding to KIT.
It is generally known that GISTs with KIT mutations exhibit a high level of phosphorylated KIT and display a constitutive activation of downstream signal transduction pathways that are relevant with malignant transformation, including PI3K/Akt, Ras/ERK, and JAK/STAT [19, 31]. Thus, effects of miR-148b-3p on various signaling cascades of KIT were demonstrated in this study, showing that the downregulation of KIT induced by miR-148b-3p dramatically suppresses the expression of phosphorylated and total ERK, AKT and STAT3 in vitro and in vivo. Moreover, restoration KIT expression notably increased the expression of phosphorylated and total ERK, AKT and STAT3 in GIST cells transfection with miR-148b-3p. These results indicated that the miR-148b-3p-dependent effects on GIST tumorigenesis are mediated by oncogenic KIT-related signal transductions.
Our study provides convincing evidence that miR-148b-3p is frequently downregulated in GISTs with KIT mutations and KIT overexpression and its expression is negatively relevance with poor prognosis of GISTs. Moreover, miR-148b-3p functions as a tumor suppressor in suppressing cell proliferation, migration, and invasion through directly targeting KIT, and it could be exploited as a prognostic biomarker to predict the aggressive behavior as well as a therapeutic target in GISTs.
False discovery rate
Gastrointestinal stromal tumors
enrichment and Kyoto Encyclopedia of Genes and Genomes
National Institute of Health
Platelet-derived growth factor receptor alpha polypeptide
Succinate dehydrogenase complex subunits B
Joensuu H, Hohenberger P, Corless CL. Gastrointestinal stromal tumour. Lancet. 2013;382:973–83.
Barnett CM, Corless CL, Heinrich MC. Gastrointestinal stromal tumors. molecular markers and genetic subtypes Hematol Oncol Clin North Am. 2013;27:871–88.
Nannini M, Biasco G, Astolfi A, Pantaleo MA. An overview on molecular biology of KIT/PDGFRA wild type (WT) gastrointestinal stromal tumours (GIST). J Med Genet. 2013;50:653–61.
Pantaleo MA, Astolfi A, Indio V, Moore R, Thiessen N, Heinrich MC, et al. SDHA loss-of-function mutations in KIT-PDGFRA wild-type gastrointestinal stromal tumors identified by massively parallel sequencing. J Natl Cancer Inst. 2011;103:983–7.
Pantaleo MA, Astolfi A, Nannini M, Ceccarelli C, Formica S, Santini D, et al. Differential expression of neural markers in KIT and PDGFRA wild-type gastrointestinal stromal tumours. Histopathology. 2011;59:1071–80.
Pantaleo MA, Astolfi A, Di Battista M, Heinrich MC, Paterini P, Scotlandi K, et al. Insulin-like growth factor 1 receptor expression in wild-type GISTs: a potential novel therapeutic target. Int J Cancer. 2009;125:2991–4.
MicroRNAs BDP. Genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Alvarez-Garcia I, Miska EA. MicroRNA functions in animal development and human disease. Development. 2005;132:4653–62.
Hwang HW, Mendell JT. MicroRNAs in cell proliferation, cell death, and tumorigenesis. Br J Cancer. 2006;94:776–80.
Ventura A, Jacks T. MicroRNAs and cancer: short RNAs go a long way. Cell. 2009;136:586–91.
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, et al. MicroRNA expression profiles classify human cancers. Nature. 2005;435:834–8.
Kang HJ, Nam SW, Kim H, Rhee H, Kim NG, Kim H, et al. Correlation of KIT and platelet-derived growth factor receptor alpha mutations with gene activation and expression profiles in gastrointestinal stromal tumors. Oncogene. 2005;24:1066–74.
Haller F, von Heydebreck A, Zhang JD, Gunawan B, Langer C, Ramadori G, et al. Localization- and mutation-dependent microRNA (miRNA) expression signatures in gastrointestinal stromal tumours (GISTs), with a cluster of co-expressed miRNAs located at 14q32.31. J Pathol. 2010;220:71–86.
Kim WK, Yang HK, Kim H. MicroRNA involvement in gastrointestinal stromal tumor tumorigenesis. Curr Pharm Des. 2013;19:1227–35.
Choi HJ, Lee H, Kim H, Kwon JE, Kang HJ, You KT, et al. MicroRNA expression profile of gastrointestinal stromal tumors is distinguished by 14q loss and anatomic site. Int J Cancer. 2010;126:1640–50.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25:25–9.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Patel S. Exploring novel therapeutic targets in GIST: focus on the PI3K/Akt/mTOR pathway. Curr Oncol Rep. 2013;15:386–95.
Quattrone A, Wozniak A, Dewaele B, Floris G, Vanspauwen V, Van Looy T, et al. Frequent mono-allelic loss associated with deficient PTEN expression in imatinib-resistant gastrointestinal stromal tumors. Mod Pathol. 2014;27:1510–20.
Ricci R, Maggiano N, Castri F, Rinelli A, Murazio M, Pacelli F, et al. Role of PTEN in gastrointestinal stromal tumor progression. Arch Pathol Lab Med. 2004;128:421–5.
Nannini M, Ravegnini G, Angelini S, Astolfi A, Biasco G, Pantaleo MA. miRNA profiling in gastrointestinal stromal tumors: implication as diagnostic and prognostic markers. Epigenomics. 2015;7:1033–49.
Kelly L, Bryan K, Kim SY, Janeway KA, Killian JK, Schildhaus HU, et al. Post-transcriptional dysregulation by miRNAs is implicated in the pathogenesis of gastrointestinal stromal tumor [GIST]. PLoS One. 2013;8:e64102.
Rubin BP, Singer S, Tsao C, Duensing A, Lux ML, Ruiz R, et al. KIT activation is a ubiquitous feature of gastrointestinal stromal tumors. Cancer Res. 2001;61:8118–21.
Koelz M, Lense J, Wrba F, Scheffler M, Dienes HP, Odenthal M. Down-regulation of miR-221 and miR-222 correlates with pronounced kit expression in gastrointestinal stromal tumors. Int J Oncol. 2011;38:503–11.
Ihle MA, Trautmann M, Kuenstlinger H, Huss S, Heydt C, Fassunke J, et al. miRNA-221 and miRNA-222 induce apoptosis via the KIT/AKT signalling pathway in gastrointestinal stromal tumours. Mol Oncol. 2015;9:1421–33.
Gits CM, van Kuijk PF, Jonkers MB, Boersma AW, van Ijcken WF, Wozniak A, et al. MiR-17-92 and miR-221/222 cluster members target KIT and ETV1 in human gastrointestinal stromal tumours. Br J Cancer. 2013;109:1625–35.
Kim WK, Park M, Kim YK, Tae YK, Yang HK, Lee JM, et al. MicroRNA-494 downregulates KIT and inhibits gastrointestinal stromal tumor cell proliferation. Clin Cancer Res. 2011;17:7584–94.
Fan R, Zhong J, Zheng S, Wang Z, Xu Y, Li S, et al. MicroRNA-218 inhibits gastrointestinal stromal tumor cell and invasion by targeting KIT. Tumour Biol. 2014;35:4209–17.
Bachet JB, Tabone-Eglinger S, Dessaux S, Besse A, Brahimi-Adouane S, Emile JF, et al. Gene expression patterns of hemizygous and heterozygous KIT mutations suggest distinct oncogenic pathways: a study in NIH3T3 cell lines and GIST samples. PLoS One. 2013;8:e61103.
Li F, Huynh H, Li X, Ruddy DA, Wang Y, Ong R, et al. FGFR-mediated reactivation of MAPK signaling attenuates antitumor effects of Imatinib in gastrointestinal stromal tumors. Cancer Discov. 2015;5:438–51.
Shen J, Hu Q, Schrauder M, Yan L, Wang D, Medico L, et al. Circulating miR-148b and miR-133a as biomarkers for breast cancer detection. Oncotarget. 2014;5:5284–94.
Sandhu R, Rivenbark AG, Coleman WB. Loss of post-transcriptional regulation of DNMT3b by microRNAs: a possible molecular mechanism for the hypermethylation defect observed in a subset of breast cancer cell lines. Int J Oncol. 2012;41:721–32.
Song Y, Xu Y, Wang Z, Chen Y, Yue Z, Gao P, et al. MicroRNA-148b suppresses cell growth by targeting cholecystokinin-2 receptor in colorectal cancer. Int J Cancer. 2012;131:1042–51.
Song YX, Yue ZY, Wang ZN, Xu YY, Luo Y, Xu HM, et al. MicroRNA-148b is frequently down-regulated in gastric cancer and acts as a tumor suppressor by inhibiting cell proliferation. Mol Cancer. 2011;10:1.
Sui C, Meng F, Li Y, Jiang Y. miR-148b reverses cisplatin-resistance in non-small cell cancer cells via negatively regulating DNA (cytosine-5)-methyltransferase 1(DNMT1) expression. J Transl Med. 2015;13:132.
Liu GL, Liu X, Lv XB, Wang XP, Fang XS, Sang Y. miR-148b functions as a tumor suppressor in non-small cell lung cancer by targeting carcinoembryonic antigen (CEA). Int J Clin Exp Med. 2014;7:1990–9.
Zhao G, Zhang JG, Liu Y, Qin Q, Wang B, Tian K, et al. miR-148b functions as a tumor suppressor in pancreatic cancer by targeting AMPKalpha1. Mol Cancer Ther. 2013;12:83–93.
This work was supported by the National Natural Science Foundation of China (No. 81302114, 31271040, 81270176, 81600485, 81000048 and 81300099).
Availability of data and materials
All data in our study are available upon request.
Ethics approval and consent to participate
The study was approved by the Ethics Committee of Tongji Hospital of Tongji Medical College, Huazhong University of Science and Technology, P.R. China. All participants provided their written informed consent to participate in this study.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.