Research Paper Volume 12, Issue 10 pp 9151—9172

Downregulation of miRNA-146a-5p promotes malignant transformation of mesenchymal stromal/stem cells by glioma stem-like cells

Xingliang Dai1,2, *, , Yunfei Wang3, *, , Xuchen Dong1, *, , Minfeng Sheng1, , Haiyang Wang1, , Jia Shi1, , Yujing Sheng1, , Liang Liu1, , Qianqian Jiang1, , Yanming Chen1, , Bingshan Wu2, , Xuejun Yang3, , Hongwei Cheng2, , Chunsheng Kang3, , Jun Dong1, ,

  • 1 Brain Tumor Lab, Department of Neurosurgery, The Second Affiliated Hospital of Soochow University, Suzhou 215004, China
  • 2 Department of Neurosurgery, The First Affiliated Hospital of Anhui Medical University, Hefei 230022, China
  • 3 Department of Neurosurgery, Tianjin Medical University General Hospital, Tianjin 300052, China
* Equal contribution

Received: November 3, 2019       Accepted: March 9, 2020       Published: May 25, 2020      

https://doi.org/10.18632/aging.103185
How to Cite

Copyright © 2020 Dai et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Mesenchymal stromal/stem cells (MSCs) are promising carriers in cell-based therapies against central nervous system diseases, and have been evaluated in various clinical trials in recent years. However, bone marrow-derived MSCs (BMSCs) are reportedly involved in tumorigenesis initiated by glioma stem-like cells (GSCs). We therefore established three different orthotopic models of GSC-MSC interactions in vivo using dual-color fluorescence tracing. Cell sorting and micropipetting techniques were used to obtain highly proliferative MSC monoclones from each model, and these cells were identified as transformed MSC lines 1, 2 and 3. Nineteen miRNAs were upregulated and 24 miRNAs were downregulated in all three transformed MSC lines compared to normal BMSCs. Reduced miR-146a-5p expression in the transformed MSCs was associated with their proliferation, malignant transformation and overexpression of heterogeneous nuclear ribonucleoprotein D. These findings suggest that downregulation of miR-146a-5p leads to overexpression of its target gene, heterogeneous nuclear ribonucleoprotein D, thereby promoting malignant transformation of MSCs during interactions with GSCs. Given the risk that MSCs will undergo malignant transformation in the glioma microenvironment, targeted glioma therapies employing MSCs as therapeutic carriers should be considered cautiously.

Introduction

In the glioma microenvironment, there are not only glioma stem-like cell (GSC)-derived glioma cells, but also various kinds of tumor stromal cells, including peritumor astrocytes, microglia, fibroblasts, macrophages, dendritic cells, neutrophils, lymphocytes and endothelial cells [1]. Mesenchymal stromal/stem cells (MSCs) can migrate directly to the glioma microenvironment, and thus have been applied as antitumor delivery vectors for glioma treatment [1]. However, some studies have demonstrated that tumor stromal cells such as MSCs in the glioma microenvironment can not only promote tumor progression, but also enhance tumor resistance to chemotherapeutics [2]. Furthermore, several kinds of tumor stromal cells (e.g., fibroblasts, macrophages and oligodendrocytes) can undergo malignant transformation after interacting with cancer cells, especially cancer stem cells; indeed, our previous report indicated that the overexpression of telomerase reverse transcriptase was associated with GSC-induced MSC transformation [3]. However, the overexpression of telomerase reverse transcriptase in normal MSCs was only mildly tumorigenic in vivo, suggesting that other mechanisms contribute to MSC transformation during GSC-initiated tissue remodeling.

Certain pro-tumor genes can be horizontally transferred from tumor cells to stromal cells, thus malignantly transforming the stromal cells [46]; however, the chromosomes may be lost after fusion [8]. Glioma cells were reported to induce the transformation of MSCs in vitro without directly contacting them [7, 8], and the interleukin-6/signal transducer and activator of transcription 3 pathway was found to be involved in this process [9]. Granulocyte-macrophage colony-stimulating factor/interleukin-4 and soluble interleukin receptor/glycoprotein 130 may also contribute to MSC transformation in vitro [10, 11]. Simple long-term culture may induce the spontaneous malignant transformation of MSCs [12], but this finding has not been widely accepted as fact [13]. Bone marrow stromal cells in the rat brain were found to undergo malignant transformation in vivo in a tumor microenvironment containing tumor stem cell niches formed by orthotopically transplanted C6 glioma cells [14]; however, it is unclear where and how bone marrow stromal cells are transformed. In summary, the mechanisms responsible for the malignant transformation of MSCs in the glioma microenvironment have not been fully elucidated.

The aberrant expression of microRNAs (miRNAs), especially oncogenic or tumor suppressor miRNAs, promotes carcinogenesis, tumor progression, malignant transformation, tumor metastasis and anticancer treatment resistance [1517]. High-throughput miRNA profiling techniques such as RNA sequencing and miRNA microarray analysis have greatly clarified the involvement of miRNAs in malignancies [18, 19]. Dysregulated miRNAs contribute to oncogenic transformation processes such as inflammation and metabolic reprogramming, thus creating a tumorigenic microenvironment that promotes the initiation and progression of neoplasms [20]. Altered miRNA expression profiles have been used to diagnose and stage various human tumors, and to predict their progression, prognosis and treatment response [21, 22]. However, further work is needed to determine the contributions of dysregulated miRNAs to the malignant transformation of MSCs, and to characterize the miRNA profiles of transformed MSCs in the glioma microenvironment.

In the current study, we established three different in vivo GSC-MSC interaction models so that we could observe the morphological and functional changes of MSCs that had interacted with GSCs. We then used RNA sequencing to analyze the miRNA profiles of the transformed MSCs, and examined the involvement of miR-146a-5p in MSC transformation both in vitro and in vivo. This study has provided novel evidence on the molecular mechanisms underlying the malignant transformation of bone marrow-derived MSCs (BMSCs), illustrating the potential risks of MSC-based clinical trials for glioma.

Results

Biomarkers of MSCs were associated with the survival of glioma patients

In this study, we obtained data on glioma patients from the databases of the Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA). We constructed Kaplan-Meier survival curves and performed a log-rank comparison to analyze the correlations between MSC surface marker levels in glioma tissues and overall glioma patient survival. The expression of the positive MSC markers CD44, CD105 and CD29 exhibited a significant inverse correlation with the prognosis of glioma patients (Figure 1A1C). Another positive marker, CD90, was not significantly associated with overall survival (Figure 1D). Additionally, the expression of the negative MSC markers CD11b, CD34 and CD45 did not correlate significantly with the prognosis of glioma patients (Figure 1E1G). These results suggested that higher MSC marker expression was associated with a shorter survival period and a poorer prognosis in glioma patients.

Biomarkers of MSCs in the glioma microenvironment. (A–D) Kaplan-Meier survival curves based on the expression of the MSC biomarkers CD44, CD105, CD29 and CD90 in glioblastoma tissues from the CGGA dataset; (E–G) Kaplan-Meier survival curves based on the expression of the MSC biomarkers CD11b, CD34 and CD45 in glioblastoma tissues from TCGA.

Figure 1. Biomarkers of MSCs in the glioma microenvironment. (AD) Kaplan-Meier survival curves based on the expression of the MSC biomarkers CD44, CD105, CD29 and CD90 in glioblastoma tissues from the CGGA dataset; (EG) Kaplan-Meier survival curves based on the expression of the MSC biomarkers CD11b, CD34 and CD45 in glioblastoma tissues from TCGA.

Preparation of the dual-color fluorescence-tracing intracranial GSC-MSC interaction model

Next, we generated GSC-MSC interaction models, as shown in the schematic diagrams in Figure 2A. Human GSCs (SU3 cells) were stably transfected with the red fluorescent protein (RFP) gene, and the SU3-RFP cells displayed 100% RFP expression (Figure 2B). We also generated transgenic enhanced green fluorescent protein (EGFP)+ BALB/c athymic nude mice, which expressed EGFP ubiquitously, except in enucleated cells (Figure 2C, left). To generate chimeric BALB/c nude mice, we used X-ray irradiation to autologously deplete the bone marrow function of normal BALB/c nude mice, and reconstructed their bone marrow by injecting their tail veins with exogenous EGFP+ bone marrow cells harvested from the EGFP+ BALB/c athymic nude mice; thus, in the chimeric mice, only bone marrow-derived cells expressed EGFP (Figure 2C, right).

Intracerebral xenografting of GSCs. (A) Schematic overview of the establishment of the dual-color fluorescence-tracing intracranial GSC-MSC interaction model; (B) SU3-RFP cells; (C) Transgenic EGFP+ nude mouse, and chimeric nude mice with autologous bone marrow deprivation and exogenous EGFP+ bone marrow cell transplantation; (D) Transplanted tumor observed under natural light, and H&E staining of a transplanted tumor section; (E) Transplanted tumor from an EGFP+ nude mouse under a fluorescence microscope; (F) Monocloned EGFP+ cells with high proliferative capacities, TMEC1, TMEC2. Scale bars: (A) 100 μm; (C) 500 μm; (D) 200 μm; (F) 50 μm.

Figure 2. Intracerebral xenografting of GSCs. (A) Schematic overview of the establishment of the dual-color fluorescence-tracing intracranial GSC-MSC interaction model; (B) SU3-RFP cells; (C) Transgenic EGFP+ nude mouse, and chimeric nude mice with autologous bone marrow deprivation and exogenous EGFP+ bone marrow cell transplantation; (D) Transplanted tumor observed under natural light, and H&E staining of a transplanted tumor section; (E) Transplanted tumor from an EGFP+ nude mouse under a fluorescence microscope; (F) Monocloned EGFP+ cells with high proliferative capacities, TMEC1, TMEC2. Scale bars: (A) 100 μm; (C) 500 μm; (D) 200 μm; (F) 50 μm.

An SU3-RFP cell suspension was injected into the cerebral caudate of the EGFP+ nude mice (Model I) and the chimeric mice (Model II) for approximately 35 days. Then, tumor-bearing mice from both models were sacrificed, and the intracranial tumors were examined under natural light and after hematoxylin and eosin (H&E) staining (Figure 2D). SU3-derived RFP+ tumor cells and EGFP+ host cells in the tumor sections could be distinguished under fluorescence microscopy (Figure 2E). Under confocal microscopy, it was evident that GSCs and MSCs directly contacted each other within the transplanted tumors (white arrow, Supplementary Figure 1). Flow cytometry (MoFlo XDP; Beckman Coulter, USA) was used to sort cells expressing EGFP, which were detected at percentages of 5.46% (Model I) and 8.64% (Model II) (Supplementary Figure 2A, 2B). A micropipetting technique was used to obtain monoclones of EGFP-expressing cells with high proliferative capacities from Model I (tumor microenvironment cells 1, TMEC1) and Model II (TMEC2) (Figure 2F).

Surgical tumor specimens from a 34-year-old male patient diagnosed with glioblastoma were transplanted into the caudate of EGFP+ nude mice via a stereotactic trocar system (1-mm3 specimens, Supplementary Figure 2C). After the mice were sacrificed, their tumors were observed under natural light (Figure 3A). H&E staining of xenograft sections revealed densely arranged tumor cells (Figure 3B), and fluorescence microscopy revealed an irregular round tumor border enclosing some EGFP+ cells (Figure 3C). After primary culture of the transplanted tumors, cells expressing EGFP were sorted, and were found to account for 29.25% of the total cells (Supplementary Figure 2D). From the tumor model based on the transplantation of surgical glioblastoma specimens (Model III), one EGFP+ cell line (TMEC3) was cloned (Figure 3D).

Intracerebral-xenograft-derived glioblastoma specimens. (A) Whole brain of a tumor-bearing mouse; (B) H&E staining; (C) Transplanted tumor viewed under a fluorescence microscope; (D) Monocloned EGFP+ TMEC3 cells. Scale bars: (B) 500 μm (left), 20 μm (right); (C) 10 μm; (D) 100 μm (left), 50 μm (right).

Figure 3. Intracerebral-xenograft-derived glioblastoma specimens. (A) Whole brain of a tumor-bearing mouse; (B) H&E staining; (C) Transplanted tumor viewed under a fluorescence microscope; (D) Monocloned EGFP+ TMEC3 cells. Scale bars: (B) 500 μm (left), 20 μm (right); (C) 10 μm; (D) 100 μm (left), 50 μm (right).

As a control experiment, we performed direct co-culture in vitro to assess whether GSCs directly interacted with MSCs. Using time-lapse photography of a living cell workstation, we did indeed observe interactions, including direct contact, between GSCs and BMSCs. We even detected the exchange of cytoplasmic substances between the cells, both through direct contact points (black arrow, Supplementary Figure 3) and through slender tubular structures (black arrow, Supplementary Figure 4) that turned yellow after the intercellular cytoplasm exchange (white arrow, Supplementary Figure 4). However, when GSCs and MSCs were indirectly co-cultured in a Transwell system in vitro, no significant changes were observed in the growth characteristics of the MSCs (Supplementary Figure 5).

Malignant transformation of BMSCs

All three cloned EGFP+ cell lines (TMEC1, TMEC2 and TMEC3) exhibited strong proliferation abilities in vitro, and their cell expansion activities were much higher than those of normal bone marrow cells, but lower than those of SU3 cells (Figure 4A). TMEC1 and TMEC2 cells exhibited high tumorigenicity 25 days after being subcutaneously transplanted (1×106 cells) into BALB/c nude mice; both cell lines had a tumor formation rate of 100% (5/5). However, subcutaneously transplanted TMEC3 cells did not form tumors in any of the five mice (Supplementary Figure 2E).

Transformation of host stromal cells. (A) Cell proliferation curve; (B) RT-PCR analysis of human and mouse β-actin expression in SU3 cells and three TMEC lines; (C) FISH assay of chromosomes in SU3 cells and transformed cells; (D) Immunofluorescence of the three tMSC lines. Scale bars: (C) 2 μm; (D) 20 μm.

Figure 4. Transformation of host stromal cells. (A) Cell proliferation curve; (B) RT-PCR analysis of human and mouse β-actin expression in SU3 cells and three TMEC lines; (C) FISH assay of chromosomes in SU3 cells and transformed cells; (D) Immunofluorescence of the three tMSC lines. Scale bars: (C) 2 μm; (D) 20 μm.

The three transformed cell lines expressed mouse β-actin, but not human β-actin (Figure 4B). A fluorescence in situ hybridization (FISH) assay of the sex chromosomes revealed that the karyotype of the SU3 cells was XY (X, red fluorescent probe; Y, green fluorescent probe) (Figure 4C), in accordance with clinical data showing that SU3 cells were derived from a male patient [23, 24]. The karyotypes of all three transformed cell lines were XX, consistent with the karyotypes of the female host mice (Figure 4C). Immunofluorescence assays demonstrated that TMEC1, TMEC2 and TMEC3 cells expressed biomarkers of MSCs, as the cells were positive for CD29, CD44, CD105 and CD90 but negative for CD11b, CD31, CD34 and CD45. Normal MSCs were used as the positive control (Figure 4D and Supplementary Figure 2F).

We then performed an osteogenic differentiation assay, and found that mineralized nodules emerged from all three cell lines after seven days of differentiation induction. When we conducted an adipogenic differentiation assay, small fat droplets gradually coalesced into large fat droplets after seven days of differentiation induction, and brown adipose tissue deposition occurred after 14 days of in vitro induction. These multidirectional differentiation assays confirmed that TMEC1, TMEC2 and TMEC3 cells were transformed BMSCs and still had differentiation potential. Specifically, TMEC1 cells exhibited obvious osteogenic and slight adipogenic differentiation potential, TMEC2 cells exhibited notable adipogenic differentiation potential, and TMEC3 cells exhibited strong osteogenic and adipogenic differentiation potential, similar to those of normal MSCs (Figure 5). These data suggested that the three transformed EGFP-expressing TMEC lines from the different models were BMSCs; thus, these cells were named transformed MSCs (tMSC1, tMSC2 and tMSC3 for Models I, II and III, respectively).

Osteogenic and adipogenic differentiation assays in tMSCs. Scale bars: 20 μm.

Figure 5. Osteogenic and adipogenic differentiation assays in tMSCs. Scale bars: 20 μm.

Differentially expressed miRNAs in tMSCs

We then performed small RNA sequencing on the tMSCs and normal BMSCs. We found that 19 miRNAs were significantly upregulated and 24 miRNAs were significantly downregulated in tMSCs (n=3) compared with normal BMSCs (n=3) (Figure 6A). Quantitative real-time PCR (qRT-PCR) was used to further assess nine of these 43 miRNAs (Figure 6B, 6C). Because the tMSC1 and tMSC2 cells were highly proliferative and exhibited strong tumorigenicity, we used these cells to examine the correlations between miRNA levels and biological functions contributing to malignant transformation. The qRT-PCR analysis results for the nine selected miRNAs were consistent with the RNA sequencing results.

MiRNA profiles of tMSCs. (A) Differentially expressed miRNAs; (B) Verification of miRNAs identified in sequencing experiments; (C) Expression of miRNA-146a-5p in tMSCs compared with normal BMSCs. * P

Figure 6. MiRNA profiles of tMSCs. (A) Differentially expressed miRNAs; (B) Verification of miRNAs identified in sequencing experiments; (C) Expression of miRNA-146a-5p in tMSCs compared with normal BMSCs. * P<0.05, *** P<0.001.

We found that miR-146a-5p expression was reduced in all three tMSC lines compared with normal MSCs (927.99-fold in tMSC1, 2849.74-fold in tMSC2 and 1.25-fold in tMSC3). Interestingly, miR-146a-5p was significantly downregulated in the tMSCs with higher tumorigenic abilities (tMSC1 and tMSC2), but only slightly downregulated in those with a lower tumorigenic ability (tMSC3 cells) (Figure 6C). These data suggested that lower miR-146a-5p expression may be associated with MSC transformation and tumorigenicity.

Overexpression of miR-146a-5p partially reversed the malignant phenotype of tMSCs

Next, we used lentiviral vectors to overexpress miR-146a-5p in the three tMSC lines. The transfection induced miR-146a-5p expression 20.06-fold, 19.02-fold and 6.79-fold in tMSC1, tMSC2 and tMSC3 cells, respectively (Figure 7A). Cell proliferation was significantly lower in miR-146a-5p-overexpressing tMSC1 (PC) cells than in the blank control (BC) and negative control (NC) cells (Figure 7B). Overexpression of miR-146a-5p in tMSC1 and tMSC2 cells reduced their colony formation abilities in a colony formation assay (P<0.001); the colony numbers of blank control, negative control and miR-146a-5p-overexpressing tMSC1 cells were 94.75±4.32, 91.50±1.50 and 23.25±2.17, while those of blank control, negative control and miR-146a-5p-overexpressing tMSC2 cells were 88.75±5.63, 83.75±4.15 and 20.50±1.12, respectively (Figure 7C). The invasion abilities of tMSC1 and tMSC2 cells also decreased significantly after miR-146a-5p overexpression (P<0.001) (Figure 7D). The tumors generated from miR-146a-5p-overexpressing tMSC1 cells were appreciably smaller in volume than the tumors generated from control cells 28 days after their subcutaneous transplantation into BALB/c nude mice (Figure 7E).

Overexpression of miR-146a-5p partially reversed the malignant phenotype of tMSCs. (A) Relative expression of miR-146a-5p after lentiviral vector transfection; (B) Cell proliferation activity of miR-146a-5p-overexpressing, negative control (NC) and blank control (BC) tMSC1 cells; (C) Colony formation assays of BC, NC and miR-146a-5p-overexpressing tMSC1 and tMSC2 cells; (D) Invasion assays of NC and miR-146a-5p-overexpressing tMSC1 and tMSC2 cells; (E) Tumorigenicity of NC and miR-146a-5p-overexpressing tMSC1 cells. * P

Figure 7. Overexpression of miR-146a-5p partially reversed the malignant phenotype of tMSCs. (A) Relative expression of miR-146a-5p after lentiviral vector transfection; (B) Cell proliferation activity of miR-146a-5p-overexpressing, negative control (NC) and blank control (BC) tMSC1 cells; (C) Colony formation assays of BC, NC and miR-146a-5p-overexpressing tMSC1 and tMSC2 cells; (D) Invasion assays of NC and miR-146a-5p-overexpressing tMSC1 and tMSC2 cells; (E) Tumorigenicity of NC and miR-146a-5p-overexpressing tMSC1 cells. * P<0.05, ** P<0.01, *** P<0.001.

HNRNPD is a direct target gene of miR-146a-5p

We then performed target gene prediction on the nine differentially expressed miRNAs validated in this study, and conducted Gene Ontology (GO) analyses (top 30 shown in Supplementary Figure 8A) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses (top 20 shown in Supplementary Figure 8B). The TargetScan prediction database identified heterogeneous nuclear ribonucleoprotein D (HNRNPD) as a target gene of miR-146a-5p (Figure 8A). In a luciferase reporter gene assay, miR-146a-5p significantly repressed the expression of HNRNPD by binding to its 3' untranslated region (UTR) (P<0.01). Luciferase activity was not detected in the negative control group (containing a mutated 3'UTR of HNRNPD) or the blank group, suggesting that miR-146a-5p specifically binds to the 3'UTR of HNRNPD (Figure 8B). HNRNPD mRNA was significantly upregulated in tMSC1 and tMSC2 cells compared with BMSCs (Figure 8C), but was not highly expressed in tMSC3 cells. Immunocytochemical staining of tMSC1 and tMSC2 cells further verified the high expression of HNRNPD at the protein level, consistent with the transcriptional results (Figure 8D).

Expression of the target gene HNRNPD in tMSCs. (A) Putative miR-146a-5p target sites in the HNRNPD 3'UTR; (B) Dual-luciferase assay; (C) qRT-PCR; (D) Immunocytochemistry of tMSC1 and tMSC2 cells. ** P

Figure 8. Expression of the target gene HNRNPD in tMSCs. (A) Putative miR-146a-5p target sites in the HNRNPD 3'UTR; (B) Dual-luciferase assay; (C) qRT-PCR; (D) Immunocytochemistry of tMSC1 and tMSC2 cells. ** P<0.01.

HNRNPD contributed to the malignant transformation of MSCs

To further study the involvement of HNRNPD in the transformation of MSCs, we used small interfering RNA (siRNA) to downregulate HNRNPD in tMSC1 cells. Western blotting confirmed that HNRNPD protein expression was downregulated in the siRNA-treated cells (Figure 9A). The siRNA also significantly inhibited cell proliferation, as shown in an MTS assay (Figure 9B).

Correlation of HNRNPD expression with MSC transformation. (A) Protein expression after siRNA interference; (B) Cell proliferation; (C) Apoptosis; (D) Intracranial transplantation and survival time. ** P

Figure 9. Correlation of HNRNPD expression with MSC transformation. (A) Protein expression after siRNA interference; (B) Cell proliferation; (C) Apoptosis; (D) Intracranial transplantation and survival time. ** P<0.01.

Flow cytometry was then used to assess changes in apoptosis and the cell cycle. The average total percentage of apoptotic cells was 63.06±3.48% in the siRNA-treated group, which was significantly higher than that in the control (6.73±0.52%) and blank control (7.76±0.99%) groups (P<0.01, Figure 9C). In a colony formation assay, the colony formation rate of siRNA-treated tMSC1 clones was 51±7.0%, which was significantly lower than the rates of the control (94.33±1.89%) and blank control (95.33±2.05%) groups (P<0.01, Supplementary Figure 6A). Significant cell cycle arrest was observed in G1 phase, and occurred in a higher proportion of siRNA-treated cells (61.01±1.22%) than control (54.58±0.97%) or blank control (53.88±1.04%) cells (P<0.01, Supplementary Figure 6B).

We also intracranially inoculated mice with tMSC1 cells that were stably transfected with HNRNPD siRNA. The tumor volume was significantly smaller and the median survival time was nine days longer in the siRNA-treated group than in the control group (Figure 9D). Thus, the results of HNRNPD downregulation were similar to those of miR-146a-5p overexpression.

To assess whether the downregulation of miR-146-5p promoted the malignant transformation of tMSCs by inducing HNRNPD expression, we overexpressed miR-146a-5p, HNRNPD or both in tMSC1 cells. The proliferation curve and clone formation results demonstrated that tMSC1 cell proliferation was significantly inhibited in the miR-146a-5p overexpression group and significantly elevated in the HNRNPD overexpression group. In the group overexpressing both miR-146a-5p and HNRNPD, the proliferation of tMSC1 cells was higher than that of the control group, indicating that the ability of miR-146a-5p to inhibit tMSC proliferation depended on its ability to downregulate HNRNPD (Supplementary Figure 7).

In conclusion, deficiency of miR-146a-5p led to the overexpression of its target gene HNRNPD, thereby promoting the malignant transformation of MSCs.

Discussion

Our dual-color fluorescence-tracing GSC-MSC interaction models revealed the infiltration of certain host stromal cells into xenograft tumors (Figure 1E, 2C). In frozen tumor sections from Model II, many exogenous bone marrow-derived cells expressing EGFP were observed, suggesting that host BMSCs can migrate from the bone marrow to the tumor parenchyma and undergo malignant transformation in the glioma microenvironment. The malignant transformation of MSCs has been reported previously [25], but the malignant transformation of MSCs recruited from the bone marrow to the glioma microenvironment has never been reported. Our model was an effective visual platform for observing tumor-stromal interactions, and thus was an excellent tool for studying the interactions between GSCs and infiltrating MSCs.

MiRNA sequencing allowed us to identify miRNAs that were differentially expressed between normal MSCs and tMSCs. The abnormal expression of these miRNAs was further verified in a qRT-PCR assay, and miR-146a-5p was found to be significantly downregulated in several kinds of transformed TMECs. However, miR-146a-5p has divergent functions in different types of tumors. For instance, miR-146a-5p functions as an oncogene when it is overexpressed in oral tumors, cervical cancer cells and papillary thyroid carcinomas, but serves as a tumor suppressor in triple-negative breast cancer, colon carcinoma, non-small cell lung cancer cells and prostate cancer [26]. Zhao et al. [27] found that knocking out miR-146a in C57BL/6 mice led to histologically and immunophenotypically defined myeloid sarcomas, and also resulted in some lymphomas with true malignancies. Contreras et al. [28] reported that miR-146a deficiency contributed to B-cell oncogenesis by inducing the early growth response-1 gene. The functions of miR-146a in malignant hematopoietic stem cells have also been studied [29]. We observed that the malignant transformation of BMSCs was associated with the downregulation of miR-146a-5p due to the resulting overexpression of HNRNPD.

Changes in HNRNPD expression are involved in the development and progression of multiple types of solid tumors. Gouble et al. [30] examined transgenic mice that overexpressed p37AUF1 (HNRNPD), and found that they overexpressed the mRNAs for multiple tumor-associated genes, including c-Myc, c-Jun, c-Fos, granulocyte macrophage colony-stimulating factor, tumor necrosis factor-α and cyclin D1, resulting in the spontaneous formation of soft tissue sarcoma. In another study, HNRNPD mRNA levels in tumor tissues and corresponding paracancerous tissues were compared in 154 patients with 19 kinds of cancers. HNRNPD was significantly upregulated in 13% of the tumors, consistent with other reports of its upregulation in esophageal, breast, skin, thyroid and liver cancers [31, 32].

In this study, the binding between miRNA-146a-5p and HNRNPD was predicted in a bioinformatic assay and confirmed in a dual luciferase assay. Suppressing the expression of HNRNPD significantly inhibited cell cycle progression and proliferation in vitro, and also increased apoptosis, inhibited tumor cell proliferation and prolonged survival in tumor-bearing mice. Our results suggested that the downregulation of miRNA-146a-5p directly enhances the malignant characteristics of tMSCs by upregulating HNRNPD.

MSCs are capable of infiltrating the glioma microenvironment, and thus can be applied as carriers of anticancer molecules targeting glioma cells; however, the fate of MSCs after their homing to gliomas remains largely unknown. According to our analyses of the CGGA and TCGA databases, the levels of the MSC biomarkers CD29, CD44 and CD105 were associated with the prognosis of glioma patients. CD44 expression was previously reported to promote BMSC migration toward gliomas [33]. Other studies [34, 35] have noted that, compared with CD90+ MSCs, CD90- MSCs promote tumor vascularization and immunosuppression, which may explain why there was no significant difference in glioma patient survival according to CD90 expression in our study. The in vitro and in vivo experiments in this paper confirmed that BMSCs migrated to the glioma microenvironment and were malignantly transformed after interacting with GSCs. The downregulation of miRNA-146a-5p promoted GSC-initiated MSC transformation by derepressing the HNRNPD gene.

In summary, BMSCs isolated from intracerebral xenograft GSC tumors were found to undergo malignant transformation, and the suppression of miR-146a-5p was found to promote this process by activating the oncogene HNRNPD. These results suggest that targeted glioma therapies using MSCs as cellular carriers should be used cautiously, due to the risk that MSCs will undergo malignant transformation in the glioma microenvironment.

Materials and Methods

Correlation studies between MSC surface marker expression and glioma patients’ prognostic information

We analyzed the association of positive (CD29, CD44, CD105 and CD90) and negative (CD11b, CD31, CD34 and CD45) MSC marker expression with glioma prognosis in both the CGGA and TCGA databases. Kaplan-Meier survival analysis was used to estimate the survival distributions, and the log-rank test was used to assess the statistical significance of survival differences between groups stratified based on the median value of each marker. Pearson correlation was used to determine the significance of differences, and t-tests were used to determine the significance of two-group comparisons. Data are presented as the mean ± standard error. A two-sided P value < 0.05 was considered significant.

Cell lines and animals

We established the RFP-expressing GSC line SU3-RFP by using a lentiviral vector to stably transfect SU3 cells with the RFP gene, as previously described [23, 24]. Cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM)/F12 medium (Gibco, USA) containing 20 ng/mL basic fibroblast growth factor (Gibco) and 20 ng/mL epidermal growth factor (Gibco).

To establish ubiquitous EGFP expression in athymic BALB/c nude mice, we crossed BALB/c nude mice with C57BL/6-EGFP mice [36]. BALB/c nude mice expressing EGFP only in their bone marrow-derived cells (named chimeric BALB/c nude mice) were established through bone marrow damage/reconstruction technology. Briefly, X-ray radiation was used to destroy the bone marrow function of normal BALB/c nude mice, and the mice were then transplanted (via tail vein injection) with 1×107 EGFP+ bone marrow-derived cells harvested from EGFP+ BALB/c nude mice. Female mice were selected for the experiment, and were raised in the specific-pathogen-free experimental animal center of Soochow University. The animal studies and clinical tumor specimen collection were conducted with the approval and permission of the Soochow University Ethics Committee. Informed consent was obtained prior to sample acquisition.

Orthotopic transplantation

We used previously described techniques to perform intracranial orthotopic tumor transplantation [37]. Briefly, a Hamilton syringe was used with stereotactic technique to slowly inject a suspension of SU3-RFP cells (1×105) in 15 μL of DMEM into the cerebral caudate of EGFP+ BALB/c nude mice (n=3, Model I) or chimeric BALB/c nude mice expressing EGFP only in their bone marrow-derived cells (n=3, Model II). For clinical tumor specimens, small fragments (approximately 1 mm3 each) of fresh surgical tumor specimens were transplanted into the cerebral caudate of EGFP+ BALB/c nude mice (n=3, Model III) via manual manipulation with a microinjection trocar system developed in-house (Chinese Patent, CN 101569564A). Tumor specimens were obtained with informed consent from a 34-year-old male patient diagnosed with glioblastoma (World Health Organization grade IV). The tumor exhibited O6-methylguanine-DNA methyltransferase promoter methylation, low O6-methylguanine-DNA methyltransferase expression and wild-type isocitrate dehydrogenase 1/2 expression.

Time-lapse photography in vitro

GSCs and MSCs were co-cultured in vitro in a living cell workstation and were imaged via time-lapse photography. Cell suspensions of GSCs and MSCs were made, and the cells were directly co-cultured at a ratio of 1:8 (GSCs : MSCs). The cell culture parameters were adjusted at the living cell workstation, the appropriate observation field was selected, and pictures were taken every 10 minutes during 96 hours of continuous culture.

Subculturing, sorting and cloning of tumor stromal cells

Tumor-bearing mice were sacrificed approximately 35 days after orthotopic xenotransplantation. After the intracranial xenograft tumors were harvested, pieces of the tumors were dissociated and subcultured in vitro, and other pieces were used for pathological examination. Frozen sections were stained with H&E, and 4',6-diamidino-2-phenylindole (KeyGen, Nanjing, China) was used for nuclear labeling. The sections were then observed under both a fluorescence microscope and a confocal laser scanning microscope (ZEISS, Germany).

For cell subculture, minced tumor tissue was washed with cold phosphate-buffered saline, digested with 0.02% trypsin and cultured in DMEM containing 10% fetal bovine serum (FBS, Gibco). After the cells had been serially subcultivated in vitro for approximately two weeks, fluorescence-activated cell sorting was used to sort suspensions of RFP+ tumor cells and EGFP+ stromal cells, and EGFP-expressing tumor stromal cells were collected. Then, the limiting dilution method and single-cell cloning techniques were used to re-purify EGFP-expressing cells with high proliferation abilities from the different animal models [23, 38]. The EGFP-expressing tumor stromal cells from Models I, II and III were named EGFP+ tumor microenvironment cells 1 (EGFP+TMEC1, EGFP+TMEC2 and EGFP+TMEC3, respectively).

Characterization and analysis of highly proliferative EGFP+ stromal cells

The proliferation of the TMECs (TMEC1, TMEC2 and TMEC3) was analyzed with a Cell Counting Kit-8 (CCK-8; Dojindo, Kumamoto, Japan) according to the manufacturer’s instructions. Primers specific for human and mouse β-actin were used for reverse transcription polymerase chain reaction (RT-PCR) analyses. The primers are listed in Supplementary Table 1. Human and mouse sex-specific FISH assays were performed to analyze whether the EGFP+ TMECs were of mouse or human origin, according to a method described previously [10]. To evaluate the tumor formation ability of the EGFP+ TMEC monoclones, we subcutaneously inoculated the right forelimbs of BALB/c nude mice with the cells (1x106 cells/mouse).

An immunofluorescence assay was used to identify cell surface markers. Antibodies against CD29, CD44, CD105, CD90, CD11b, CD31, CD34 and CD45 (Abcam, UK) were applied to detect the expression of the relevant surface markers, in accordance with the manufacturer’s instructions. Cy3-labeled goat anti-rabbit or anti-rat IgG (Beyotime, Shanghai, China) was applied as a secondary antibody. The expression of the cell surface markers was observed under a fluorescence microscope (Olympus IX51; Tokyo, Japan), and images were acquired and merged in the software in accordance with the manufacturer’s instructions.

Osteogenic and adipogenic differentiation assays were performed to evaluate the potential for pluripotent differentiation. Cells were plated in six-well plates (coated with 0.1% gelatin) at a density of 5×105 cells/well in DMEM medium. After 24 h, the medium was replaced with a conditioned medium designed to induce osteogenic differentiation (Cyagen, Guangzhou, China). The medium was replaced every other day. Specific staining was used to assess differentiation (Alizarin red staining for osteogenic differentiation and Oil-Red-O staining for adipogenic differentiation) on days 7, 14 and 21. D1 ORL UVA mouse MSCs (American Type Culture Collection, P4) were used as the positive control.

Small RNA sequencing and identification of differentially expressed miRNAs

RNA was extracted with a mirVanaTM miRNA Isolation Kit (Ambion, Austin, TX, USA) according to the manufacturer's instructions. The RNA quality was evaluated with an Agilent 2100 Bioanalyzer (Agilent Technologies, USA). Total RNA (10 μg) was used to construct a small RNA sequencing library according to the Illumina miRNA sample preparation protocol. As we have previously reported [39], polyacrylamide gel electrophoresis was used to purify and enrich RNA (16-30 nucleotides). Then, the 5' and 3' termini of the RNA were ligated to proprietary adaptors, and the RNA was transcribed to cDNA. Small RNA sequencing was performed with the high-throughput sequencing technology developed by Illumina.

The sequence annotation, data analysis and screening of differentially expressed miRNAs were performed in accordance with the manufacturer’s manual. Briefly, after the adaptors and low-quality sequences had been filtered out, BLAST (version 2.2.11) was used to map the reads against all the annotated mature miRNA sequences in mice (Sanger miRBase 17.0). The miRNA sequence reads were normalized as reads per million (RPM) values according to the following formula: RPM = (normalized reads/total miRNA matches) x 1,000,000. MiRNAs were deemed to be significantly differentially expressed if the RPM was >2 or <0.5.

qRT-PCR analysis of miRNA and target genes

To verify the sequencing data, we selected eight differentially expressed miRNAs (miR-27b-3p, miR-22-3p, miR-146a-5p, miR-146b-5p, miR-148-3p, miR-181a-5p, miR-181b-5p and miR-142-3p) and one non-differentially expressed miRNA (miR-486-5p) for qRT-PCR validation. The PCR Universal Reagents and primers specific for the mature miRNAs and U6 snRNA were purchased from GenePharma (Suzhou, China). An ABI StepOne system (Applied Biosystems, Foster City, CA, USA) was used to perform the qRT-PCR.

Transcript-level target gene quantification was also performed via qRT-PCR. The primers (Supplementary Table 1) were designed and synthesized by Genewiz Biotech (Suzhou, China). Each sample was analyzed in triplicate, and the experiment was repeated three times. The miRNA and mRNA levels were normalized to those of snRNA U6 and β-actin, respectively, and the 2-ΔΔCt method was used for quantification.

MiRNA-146a-5p overexpression

To establish an EGFP+ TMEC line overexpressing mature miR-146a-5p miRNA, we used pGLV10/U6/RFP/Puro lentiviral vectors containing miR-146a-5p sequences (generated by GenePharma, Suzhou, China). The empty pGLV10/U6/RFP/Puro vector (G05AZ, GenePharma) was used as the control. An infection assay was performed with biohazard safety equipment, in accordance with the manufacturer’s instructions. Stable cell lines were selected with 10 μg/mL puromycin.

Invasion and colony formation assays in vitro

A cell invasion assay was performed with Transwell chambers (Corning, NY, USA), as described previously [38]. Briefly, the upper surface of the chamber was coated with a thin layer of Matrigel matrix (BD Bioscience). Cells (3×104) in 120 μL of DMEM without FBS were seeded into the upper chambers, and the lower chambers were filled with DMEM containing 20% FBS. After incubation at 37.0 °C for 36 hours, the upper chambers were removed and fixed with 4% paraformaldehyde. The cells that had migrated to the lower surface were stained with 0.1% crystal violet and counted under a microscope (Olympus IX51, Japan).

For the colony formation assay, cells were diluted to 200 cells/mL in DMEM containing 10% FBS, and were added to each well of a six-well plate. After two weeks of culture, the cell colonies were fixed with 4% paraformaldehyde and stained with Giemsa (Solarbio, Shanghai, China).

Tumorigenicity assay

To compare the in vivo tumor formation capabilities of transformed cells with or without miR-146a-5p overexpression, we subcutaneously injected 1x106 TMEC2 cells, TMEC2 cells expressing the negative control vector or TMEC2 cells overexpressing miR-146a-5p into the armpit of the right forelimb (n=5) of BALB/c nude mice. The lengths and widths of the solid tumors were measured every seven days. The tumor volume was calculated with the formula: tumor volume = length×width2/2.

Luciferase assays

A luciferase reporter gene assay was performed in 293T cells to detect the binding between miR-146a-5p and the HNRNPD 3'UTR. A 1406-bp segment of the HNRNPD 3'UTR was cloned into the psi-CHECK2 vector (Promega, USA). TMECs (2×104 cells/well) were incubated in 24-well plates for 24 h. Then, Lipofectamine 2000 (Invitrogen) was used to transfect the cells with either the control vector or the miR-146a-5p luciferase reporter gene vector. After 48 hours, the cells were lysed, the substrate was added, and luminescence was measured with a luciferase assay system on a GloMax bioluminescence detector (Promega, USA), in accordance with the manufacturer’s protocol.

SiRNA for silencing HNRNPD

An RNA interference assay was performed to verify the function of HNRNPD. SiRNAs were purchased from Sigma-Aldrich, and the sequences are shown in Supplementary Table 2. TMEC1 and TMEC2 cells were seeded into six-well plates at a density of 5×104 cells per well. The transfection experiment was carried out when the cells were 40% confluent. Lipofectamine™ and RNAiMAX were mixed into Opti-MEM, and the mixture was diluted to three graded concentrations: 25 nM, 50 nM and 100 nM. Six hours later, the transfection reagents were replaced with fresh medium. After 24 hours, qRT-PCR and Western blotting were used to determine the siRNA interference efficiency. An MTS assay was used to assess cell viability; flow cytometry was used to detect the cell cycle distribution and apoptosis; and a clonogenicity assay was performed as described previously [23].

Overexpression of HNRNPD for the functional recovery assay

A pLenO-RIP plasmid vector, which carries a puromycin-resistant gene, was used to overexpress HNRNPD in TMEC1 cells according to the manufacturer’s instructions (Invabio Biotechnology, Shanghai, China). Then, CCK-8 and colony formation assays were performed as above.

Statistical analysis

Each experiment was performed three times, and the data are presented as the mean ± standard deviation. Statistical significance was calculated with Student’s t-test for two-group comparisons. P < 0.05 was considered statistically significant. Data analyses were performed in GraphPad Prism 5.0 software.

Ethics approval

The procedures of this study were approved by the Institutional Review Board of The Second Affiliated Hospital of Soochow University.

Abbreviations

GSCs: glioma stem-like cells; MSCs: mesenchymal stromal/stem cells; RFP: red fluorescent protein; EGFP: enhanced green fluorescent protein; tMSCs: transformed MSCs; FISH: fluorescence in situ hybridization.

Author Contributions

DXL, WYF and DXC performed the in vitro assays; WHY, SJ, SMF and CHW collected clinical samples and analyzed the data; SYJ, LL and JQQ performed the in vivo experiments; DXL, CYM and WBS wrote the manuscript; DJ designed the study; and DJ and YXJ reviewed the manuscript. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Funding

This study was supported by the Clinical Special Diagnosis and Treatment Technology project in Suzhou (No. LCZX201807), the National Natural Science Foundation of China (No. 81702457), the China National Nuclear Corporation Youth Innovation Team Project and the Suzhou Science and Technology Project (No. SYS201723).

References

  • 1. Zhang Q, Xiang W, Yi DY, Xue BZ, Wen WW, Abdelmaksoud A, Xiong NX, Jiang XB, Zhao HY, Fu P. Current status and potential challenges of mesenchymal stem cell-based therapy for malignant gliomas. Stem Cell Res Ther. 2018; 9:228. https://doi.org/10.1186/s13287-018-0977-z [PubMed]
  • 2. Saijo A, Goto H, Nakano M, Mitsuhashi A, Aono Y, Hanibuchi M, Ogawa H, Uehara H, Kondo K, Nishioka Y. Bone marrow-derived fibrocytes promote stem cell-like properties of lung cancer cells. Cancer Lett. 2018; 421:17–27. https://doi.org/10.1016/j.canlet.2018.02.016 [PubMed]
  • 3. Zhao Y, Chen J, Dai X, Cai H, Ji X, Sheng Y, Liu H, Yang L, Chen Y, Xi D, Sheng M, Xue Y, Shi J, et al. Human glioma stem-like cells induce malignant transformation of bone marrow mesenchymal stem cells by activating TERT expression. Oncotarget. 2017; 8:104418–29. https://doi.org/10.18632/oncotarget.22301 [PubMed]
  • 4. Goldenberg DM. Horizontal transmission of malignancy by cell-cell fusion. Expert Opin Biol Ther. 2012; 12:S133–9. https://doi.org/10.1517/14712598.2012.671807 [PubMed]
  • 5. Goldenberg DM, Rooney RJ, Loo M, Liu D, Chang CH. In-vivo fusion of human cancer and hamster stromal cells permanently transduces and transcribes human DNA. PLoS One. 2014; 9:e107927. https://doi.org/10.1371/journal.pone.0107927 [PubMed]
  • 6. Shi S, Zhang M, Guo R, Miao Y, Li B. Bone Marrow-Derived Mesenchymal Stem Cell-Mediated Dual-Gene Therapy for Glioblastoma. Hum Gene Ther. 2019; 30:106–17. https://doi.org/10.1089/hum.2018.092 [PubMed]
  • 7. Liu J, Zhang Y, Bai L, Cui X, Zhu J. Rat bone marrow mesenchymal stem cells undergo malignant transformation via indirect co-cultured with tumour cells. Cell Biochem Funct. 2012; 30:650–56. https://doi.org/10.1002/cbf.2844 [PubMed]
  • 8. Duvvuru MK, Han W, Chowdhury PR, Vahabzadeh S, Sciammarella F, Elsawa SF. Bone marrow stromal cells interaction with titanium; Effects of composition and surface modification. PLoS One. 2019; 14:e0216087. https://doi.org/10.1371/journal.pone.0216087 [PubMed]
  • 9. Cui X, Liu J, Bai L, Tian J, Zhu J. Interleukin-6 induces malignant transformation of rat mesenchymal stem cells in association with enhanced signaling of signal transducer and activator of transcription 3. Cancer Sci. 2014; 105:64–71. https://doi.org/10.1111/cas.12313 [PubMed]
  • 10. Hou Z, Cui N, Cui Y, Xing H, Liu W. Correlation and malignant transformation of MSCs with the overexpression of SIL-R/GP130 under the tumor microenvironment. Med Oncol. 2015; 32:328. https://doi.org/10.1007/s12032-014-0328-6 [PubMed]
  • 11. Zhou XG, Yang Y, Yang JS, Zhou J, Fang TL, Dai WD, Chen ZR. Granulocyte-macrophage colony-stimulating factor and interleukin 4 induce the malignant transformation of the bone marrow-derived human adult mesenchymal stem cells. Chin Med J (Engl). 2011; 124:729–33. [PubMed]
  • 12. Røsland GV, Svendsen A, Torsvik A, Sobala E, McCormack E, Immervoll H, Mysliwietz J, Tonn JC, Goldbrunner R, Lønning PE, Bjerkvig R, Schichor C. Long-term cultures of bone marrow-derived human mesenchymal stem cells frequently undergo spontaneous malignant transformation. Cancer Res. 2009; 69:5331–39. https://doi.org/10.1158/0008-5472.CAN-08-4630 [PubMed]
  • 13. Gou S, Wang C, Liu T, Wu H, Xiong J, Zhou F, Zhao G. Spontaneous differentiation of murine bone marrow-derived mesenchymal stem cells into adipocytes without malignant transformation after long-term culture. Cells Tissues Organs. 2010; 191:185–92. https://doi.org/10.1159/000240246 [PubMed]
  • 14. He Q, Zou X, Duan D, Liu Y, Xu Q. Malignant transformation of bone marrow stromal cells induced by the brain glioma niche in rats. Mol Cell Biochem. 2016; 412:1–10. https://doi.org/10.1007/s11010-015-2602-0 [PubMed]
  • 15. Nguyen DD, Chang S. Development of novel therapeutic agents by inhibition of oncogenic microRNAs. Int J Mol Sci. 2017; 19:65. https://doi.org/10.3390/ijms19010065 [PubMed]
  • 16. Vannini I, Fanini F, Fabbri M. Emerging roles of microRNAs in cancer. Curr Opin Genet Dev. 2018; 48:128–33. https://doi.org/10.1016/j.gde.2018.01.001 [PubMed]
  • 17. Liu L, Shi Y, Shi J, Wang H, Sheng Y, Jiang Q, Chen H, Li X, Dong J. The long non-coding RNA SNHG1 promotes glioma progression by competitively binding to miR-194 to regulate PHLDA1 expression. Cell Death Dis. 2019; 10:463. https://doi.org/10.1038/s41419-019-1698-7 [PubMed]
  • 18. Di Leva G, Croce CM. miRNA profiling of cancer. Curr Opin Genet Dev. 2013; 23:3–11. https://doi.org/10.1016/j.gde.2013.01.004 [PubMed]
  • 19. Wu HY, Wei Y, Pan SL. Down-regulation and clinical significance of miR-7-2-3p in papillary thyroid carcinoma with multiple detecting methods. IET Syst Biol. 2019; 13:225–33. https://doi.org/10.1049/iet-syb.2019.0025 [PubMed]
  • 20. Tomasetti M, Amati M, Santarelli L, Neuzil J. MicroRNA in metabolic re-programming and their role in tumorigenesis. Int J Mol Sci. 2016; 17:754. https://doi.org/10.3390/ijms17050754 [PubMed]
  • 21. Yuan HL, Wang T, Zhang KH. MicroRNAs as potential biomarkers for diagnosis, therapy and prognosis of gastric cancer. Oncotargets Ther. 2018; 11:3891–900. https://doi.org/10.2147/OTT.S156921 [PubMed]
  • 22. Frixa T, Donzelli S, Blandino G. Oncogenic microRNAs: key players in malignant transformation. Cancers (Basel). 2015; 7:2466–85. https://doi.org/10.3390/cancers7040904 [PubMed]
  • 23. Wang A, Dai X, Cui B, Fei X, Chen Y, Zhang J, Zhang Q, Zhao Y, Wang Z, Chen H, Lan Q, Dong J, Huang Q. Experimental research of host macrophage canceration induced by glioma stem progenitor cells. Mol Med Rep. 2015; 11:2435–42. https://doi.org/10.3892/mmr.2014.3032 [PubMed]
  • 24. Wan Y, Fei XF, Wang ZM, Jiang DY, Chen HC, Yang J, Shi L, Huang Q. Expression of miR-125b in the new, highly invasive glioma stem cell and progenitor cell line SU3. Chin J Cancer. 2012; 31:207–14. https://doi.org/10.5732/cjc.011.10336 [PubMed]
  • 25. Vishnubalaji R, Elango R, Al-Toub M, Manikandan M, Al-Rikabi A, Harkness L, Ditzel N, Atteya M, Hamam R, Alfayez M, Aldahmash A, Kassem M, Alajez NM. Neoplastic Transformation of Human Mesenchymal Stromal Cells Mediated via LIN28B. Sci Rep. 2019; 9:8101. https://doi.org/10.1038/s41598-019-44536-1 [PubMed]
  • 26. Kopcalic K, Petrovic N, Stanojkovic TP, Stankovic V, Bukumiric Z, Roganovic J, Malisic E, Nikitovic M. Association between miR-21/146a/155 level changes and acute genitourinary radiotoxicity in prostate cancer patients: A pilot study. Pathol Res Pract. 2019; 215:626–31. https://doi.org/10.1016/j.prp.2018.12.007 [PubMed]
  • 27. Zhao JL, Rao DS, Boldin MP, Taganov KD, O’Connell RM, Baltimore D. NF-kappaB dysregulation in microRNA-146a-deficient mice drives the development of myeloid malignancies. Proc Natl Acad Sci USA. 2011; 108:9184–89. https://doi.org/10.1073/pnas.1105398108 [PubMed]
  • 28. Contreras JR, Palanichamy JK, Tran TM, Fernando TR, Rodriguez-Malave NI, Goswami N, Arboleda VA, Casero D, Rao DS. MicroRNA-146a modulates B-cell oncogenesis by regulating Egr1. Oncotarget. 2015; 6:11023–37. https://doi.org/10.18632/oncotarget.3433 [PubMed]
  • 29. Zhao JL, Starczynowski DT. Role of microRNA-146a in normal and malignant hematopoietic stem cell function. Front Genet. 2014; 5:219. https://doi.org/10.3389/fgene.2014.00219 [PubMed]
  • 30. Gouble A, Grazide S, Meggetto F, Mercier P, Delsol G, Morello D. A new player in oncogenesis: AUF1/hnRNPD overexpression leads to tumorigenesis in transgenic mice. Cancer Res. 2002; 62:1489–95. [PubMed]
  • 31. Malik N, Yan H, Moshkovich N, Palangat M, Yang H, Sanchez V, Cai Z, Peat TJ, Jiang S, Liu C, Lee M, Mock BA, Yuspa SH, et al. The transcription factor CBFB suppresses breast cancer through orchestrating translation and transcription. Nat Commun. 2019; 10:2071. https://doi.org/10.1038/s41467-019-10102-6 [PubMed]
  • 32. Yang WH, Ding MJ, Cui GZ, Yang M, Dai DL. Heterogeneous nuclear ribonucleoprotein M promotes the progression of breast cancer by regulating the axin/β-catenin signaling pathway. Biomed Pharmacother. 2018; 105:848–55. https://doi.org/10.1016/j.biopha.2018.05.014 [PubMed]
  • 33. Yin Q, Zhou YY, Wang P, Ma LI, Li P, Wang XG, She CH, Li WL. CD44 promotes the migration of bone marrow-derived mesenchymal stem cells toward glioma. Oncol Lett. 2016; 11:2353–58. https://doi.org/10.3892/ol.2016.4270 [PubMed]
  • 34. Zhang Q, Yi DY, Xue BZ, Wen WW, Lu YP, Abdelmaksou A, Sun MX, Yuan DT, Zhao HY, Xiong NX, Xiang W, Fu P. CD90 determined two subpopulations of glioma-associated mesenchymal stem cells with different roles in tumour progression. Cell Death Dis. 2018; 9:1101. https://doi.org/10.1038/s41419-018-1140-6 [PubMed]
  • 35. Svensson A, Ramos-Moreno T, Eberstål S, Scheding S, Bengzon J. Identification of two distinct mesenchymal stromal cell populations in human malignant glioma. J Neurooncol. 2017; 131:245–54. https://doi.org/10.1007/s11060-016-2302-y [PubMed]
  • 36. Dong J, Dai XL, Lu ZH, Fei XF, Chen H, Zhang QB, Zhao YD, Wang ZM, Wang AD, Lan Q, Huang Q. Incubation and application of transgenic green fluorescent nude mice in visualization studies on glioma tissue remodeling. Chin Med J (Engl). 2012; 125:4349–54. [PubMed]
  • 37. Fei XF, Zhang QB, Dong J, Diao Y, Wang ZM, Li RJ, Wu ZC, Wang AD, Lan Q, Zhang SM, Huang Q. Development of clinically relevant orthotopic xenograft mouse model of metastatic lung cancer and glioblastoma through surgical tumor tissues injection with trocar. J Exp Clin Cancer Res. 2010; 29:84. https://doi.org/10.1186/1756-9966-29-84 [PubMed]
  • 38. Chen Y, Wang Z, Dai X, Fei X, Shen Y, Zhang M, Wang A, Li X, Wang Z, Huang Q, Dong J. Glioma initiating cells contribute to malignant transformation of host glial cells during tumor tissue remodeling via PDGF signaling. Cancer Lett. 2015; 365:174–81. https://doi.org/10.1016/j.canlet.2015.05.026 [PubMed]
  • 39. Hua D, Mo F, Ding D, Li L, Han X, Zhao N, Foltz G, Lin B, Lan Q, Huang Q. A catalogue of glioblastoma and brain MicroRNAs identified by deep sequencing. OMICS. 2012; 16:690–99. https://doi.org/10.1089/omi.2012.0069 [PubMed]