Abstract
Glioblastoma multiform (GBM) is the most lethal intracranial tumor in adults. Glioblastoma stem-like cells (GSCs) are responsible for tumorigenesis and chemotherapy resistance. BMPs are known to increase temozolomide (TMZ) response in GSCs, however, the intracellular molecular mechanism remains largely unknown. In this study, we built a GSC cell model called U87S, and performed RNA sequencing to identify differentially expressed (DE) miRNA profiles in U87S cells treated with BMP2, TMZ or combined BMP2 and TMZ respectively. Bioinformatics analysis revealed that most DE miRNAs were involved in the cancer pathways, suggesting their crucial roles in gliomagenesis. Eight miRNAs from RNA-seq were validated. Four out of these miRNAs (has-miR-199a-3p, hsa-miR-374b-5p, hsa-miR-320d, and hsa-miR-339-5p) were found significantly up-regulated in GBM tumor tissues. One of them, hsa-miR-199a-3p, was significantly correlated with the survival of GBM patients, and differentially expressed in U87S cells. Expression of hsa-miR-199a-3p was up-regulated by BMP. Overexpression of hsa-miR-199a-3p in U87S cells inhibited cell viability and enhanced the cytotoxicity of TMZ. And activation of BMP boosted the effect of hsa-miR-199a-3p on cell viability and TMZ-mediated cytotoxicity. Besides, expressions of five predicted targets of hsa-miR-199a-3p were evaluated. Four of them were differentially expressed in GBM tumors. And one of them, SLC22A18, was associated with the survival of GBM patients. In the end, a hsa-miR-199a-3p-mediated ceRNA network was constructed for the convenience of future study. Together, our data provided DE miRNA expression profiles associated with BMP2 and TMZ in GSCs, which might lead to finding out miRNA-based target therapies that specially target GSCs.