HMGA1 stimulates MYH9-dependent ubiquitination of GSK-3β via PI3K/Akt/c-Jun signaling to promote malignant progression and chemoresistance in gliomas

Myosin heavy chain 9 (MYH9) plays an essential role in human diseases, including multiple cancers; however, little is known about its role in gliomas. In the present study, we revealed that HMGA1 and MYH9 were upregulated in gliomas and their expression correlated with WHO grade, and HMGA1 promoted the acquisition of malignant phenotypes and chemoresistance of glioma cells by regulating the expression of MYH9 through c-Jun-mediated transcription. Moreover, MYH9 interacted with GSK-3β to inhibit the expression of GSK-3β protein by promoting its ubiquitination; the downregulation of GSK-3β subsequently promoted the nuclear translocation of β-catenin, enhancing growth, invasion, migration, and temozolomide resistance in glioma cells. Expression levels of HMGA1 and MYH9 were significantly correlated with patient survival and should be considered as independent prognostic factors. Our findings provide new insights into the role of HMGA1 and MYH9 in gliomagenesis and suggest the potential application of HMGA1 and MYH9 in cancer therapy in the future.


INTRODUCTION
Glioma is the most common primary brain malignancy. The current main therapy is maximal safe resection, followed by radiotherapy in combination with chemotherapy [1][2][3][4]. Despite consistent progress in treatment, the therapeutic effect is barely ameliorated. Owing to its high malignancy, glioma is hard to cure and easily recurrent. High malignancy indicates rapid and infiltrative growth, formidable invasive and migrant capability of tumor cells. Moreover, the malignant progression of gliomas is hypothesized to be induced by chemotherapy resistance. Therefore, it is an important direction to study high proliferation, invasion, and migration, as well as chemoresistance in gliomas.
A collection of studies recently demonstrated that high mobility group AT-hook 1 (HMGA1), which acts as a tumor promoter, is associated with several types of human cancers, including nonsmall cell lung cancer [5], breast cancer [6], and cervical cancer [7]. Recent studies have suggested the important roles of HMGA1 in gliomas. Its expression correlated significantly with glioma malignancy, proliferation, invasion, and angiogenesis of gliomas [8,9]. Although we have previously identified an oncogenic ZEB2/ miR-637/HMGA1 signaling axis targeting Vimentin which could promote both migration and invasion in glioma [10], little is known about the role of HMGA1 which is worthy of continuous investigation.
Myosin heavy chain 9 (MYH9), which belongs to Myosin family, is involved in cell adhesion and migration [11] and plays essential roles in multiple human diseases [12]. It is overexpressed in multiple cancers, including gastric cancer [13], non-small cell lung cancer [14], colon carcinoma [15], and breast cancer [16]. It promotes the tumorigenesis and development of those cancers and has a positive correlation with the prognosis of patients. Recently, it was confirmed that MYH9 binds to the CTNNB1 promoter to promote CTNNB1 transcription, conferring resistance to anoikis in gastric cancer [17]. Targeting MYH9 blocked HBX-induced GSK-3β ubiquitination to activate the β-catenin destruction complex, suppressing hepatocellular cancer stemness and the epithelial-to-mesenchymal transition [18]. For gliomas, overexpression of MYH9 contributes to cell migration ability [19], suggesting major roles in cell migration and tumor invasion [20]. However, the expression and role of MYH9 in gliomas are still undetermined.
In the present study, we showed that upregulation of HMGA1 and MYH9 in human gliomas is correlated with tumor progression, and knockdown of HMGA1 can inhibit the proliferation, migration, and invasion, as well as chemoresistance in glioma cells by stimulating MYH9-mediated ubiquitinated degradation of GSK-3β through PI3K/Akt/c-Jun pathway in vitro and in vivo. Our data tend to provide new insights into the molecular function of HMGA1 and MYH9 as well as its regulatory mechanisms in gliomas.

Upregulation of HMGA1 and MYH9 is correlated with glioma progression
To determine the roles of HMGA1 and MYH9 in gliomagenesis, the protein and mRNA expression of HMGA1 and MYH9 in clinical sample tissues were detected by quantitative real-time PCR (qPCR) T. Que et al. and immunohistochemistry (IHC) assays. QPCR showed that both HMGA1 and MYH9 mRNA levels were upregulated in glioma tissues compared to normal brain tissues (NB), and overexpressed levels of HMGA1 and MYH9 in glioma patients were positively correlated with pathological classification (WHO 2 vs. WHO 3 vs. WHO 4) (P < 0.05) ( Fig. 1A and B). Further analysis confirmed HMGA1 expression was positive correlated with MYH9 expression (r = 0.8356, P < 0.0001) (Fig. 1C). IHC revealed that the protein expression of HMGA1 and MYH9 also was increased in glioma tissues compared with NB.
In diffuse astrocytoma and glioblastoma, strong cytoplasmic and nuclei staining for HMGA1 protein was observed, whereas only slight nuclear HMGA1 expression was detected in NB ( Fig. 1D). Strong MYH9 staining was primarily observed in the cytoplasm and secondarily in the nucleus in glioma tissues, whereas low expression was observed in NB tissues (Fig. 1E). Subsequently, the data for HMGA1 and MYH9 mRNA expression in gliomas were extracted from the Chinese Glioma Genome Atlas (CCGA) database [21,22]. Both HMGA1 and MYH9 mRNA levels were significantly increased with increasing WHO grade ( Fig. 1F and G), and HMGA1 expression was positively correlated with MYH9 expression (r = 0.28, P < 0.01) (Fig. 1H). These results suggested that HMGA1 and MYH9 were upregulated in glioma tissues and correlated with the progression of gliomas.
HMGA1 reduces glioma sensitivity to temozolomide In our previous work, we demonstrated that overexpressed HMGA1 promoted proliferation, invasion, and migration in glioma cells, and Vimentin could directly bind to HMGA1 [10]. To further determine whether HMGA1 inhibition sensitizes U87 and LN229 cells to TMZ-induced cell death, dose-response TMZ curves were assessed for cells overexpressing HMGA1 and control (Mock) groups (Fig. 1I). The TMZ concentration causing IC 50 in U87 cells increased from 483.00 to 779.39 μM with HMGA1 overexpression (Fig. 1I, left), while HMGA1 overexpression in LN229 cells increased the TMZ IC 50 from 226.92 to 387.45 μM (Fig. 1I right). Importantly, we found that knockdown of HMGA1 increased the expression of γ-H2AX, a biomarker for DNA double-strand breaks (Supplemental Fig. 1A). These data suggested that HMGA1 inhibition sensitizes glioblastoma cells to TMZ-induced apoptosis.
Then, a western blot assay was performed to explore the mechanism of HMGA1 ( Fig. 1J and Supplemental Fig. 1B). Phosphorylation of PI3K and Akt was inactivated by HMGA1 knockdown. Silencing of HMGA1 downregulated the expression of c-Jun, which is a downstream factor in the PI3K/Akt pathway while upregulating GSK-3β. Moreover, expression of epithelialmesenchymal transition (EMT)-related factors, including β-catenin and Vimentin, and cell cycle-related proteins such as Cyclin D1 and CDK4, were reduced after inhibition of HMGA1, as well as stemness-related factors Oct-4 and Sox2. Interestingly, we found that the expression of MYH9 was decreased by HMGA1 silencing, indicating that HMGA1 might directly or indirectly regulate MYH9 expression.

HMGA1 increases MYH9 expression through c-Jun-mediated transcription
To confirm the relationship between HMGA1 and MYH9, we performed qPCR and western blotting. As shown in Fig. 1J and Fig. 2A, HMGA1 silencing led to a decrease in both mRNA and protein expressions of MYH9. However, bioinformatics data suggested that HMGA1, as a transcription factor, could not bind to the promoter region of MYH9 (data was not shown). Using the UCSC Genome Browser and PROMO, we found that the transcription factor c-Jun might bind to the promoter of MYH9, with three c-Jun-binding motifs at −171 to −159 (Site 1), −1547 to −1535 (Site 2), and −1608 to −1596 (Site 3) of the transcription start site present in the promoter region (Fig. 2B).
To validate these findings, we silenced c-Jun expression in glioma cells; and qPCR and western blotting revealed that knockdown of c-Jun markedly reduced MYH9 expression ( Fig. 2C and D). Then, ChIP assay was performed to examine potential protein-DNA interactions between c-Jun and the transcription regulatory region of MYH9. We found that c-Jun bound to all predicted sites inside the transcription regulatory region of MYH9 in U87 cells; moreover, downregulation of c-Jun inhibited the protein-DNA interaction at all three binding sites ( Fig. 2E and Supplemental Fig. 1C). A dual-luciferase reporter assay was performed to confirm whether c-Jun positively mediates the expression of MYH9. We subcloned the promoter region of MYH9, including the predicted c-Jun recognition site or mutated sequences into luciferase reporter plasmids. The reporter plasmid with the promoter region of MYH9 resulted in a significant decrease in luciferase activity after transfection with si-c-Jun, whereas the plasmid without the MYH9 promoter showed no change in luciferase activity (Fig. 2F).
MYH9 promotes proliferation, invasion, migration, and TMZ resistance in gliomas We transfected U87 and LN229 cells with a lentiviral vector carrying the MYH9-shRNA vector (shMYH9) and a control vector (Mock). 3-(4,5-Dimethylthiazol-2-yl)−2,5-diphenyltetrazolium bromide (MTT) assays ( Fig. 3A and B) and Edu incorporation assays (Fig. 3C) showed that MYH9 downregulation inhibited cell growth by blocking the G1/S phase transition. U87-shMYH9 cells and control cells were injected into the cerebrum of nude mice. The mice injected with U87-shMYH9 cells had a smaller mean tumor volume than controls (Fig. 3D). Transwell and Boyden chamber assays revealed that MYH9 significantly increased invasiveness and migration over controls ( Fig. 3E and F). TMZ dose-response curves showed that the IC 50 in U87 cells decreased from 505.20 to 250.05 μM after MYH9 knockdown (Fig. 3G); MYH9 knockdown in LN229 cells reduced the IC 50 from 232.61 to 142.40 μM (Fig. 3H). We used the MYH9-overexpression vector (MYH9 OE) to identify the phenotype of U87 and LN229 cells after functional recovery of MYH9, and the results also suggested that MYH9 promoted proliferation, invasion, migration, and chemotherapy resistance in glioma cells (Supplemental Fig. 2).
Subsequently, western blotting was performed to explore the mechanism of MYH9 and to find the downstream pathway. Fig. 1 Upregulation of HMGA1 and MYH9 was correlated with the progression of gliomas. A, B HMGA1 (A) and MYH9 (B) mRNA expression were upregulated in glioma tissues compared to normal brain tissues, and overexpressed levels of HMGA1 and MYH9 in glioma patients were positively correlated with the status of pathological classification. C The expression of HMGA1 was positively related to MYH9 expression. D Immunohistochemistry assay identified HMGA1 was overexpressed in glioma sample tissues compared to normal brain tissues. E MYH9 was overexpressed in glioma sample tissues. Strong expression of MYH9 was mainly observed in cytoplasmic (B, D, and E) and partly found in the nucleus (C and F) in glioma tissues. F, G Analysis of Chinese Glioma Genome Atlas (CCGA) database suggested HMGA1 (F) and MYH9 (G) mRNA level was significantly increased with the increase of WHO grade. H The expression of HMGA1 was positively related to MYH9 expression in CCGA database. I Dose-response curves of TMZ treatment were examined for HMGA1-overexpression groups and the control groups. J Western blot assay confirmed the change of the downstream pathway of HMGA1.
Inhibition of MYH9 led to the upregulation of GSK-3β expression in both U87 and LN229 cells, similar to the effect of HMGA1 (Fig. 3I). The expression of EMT-related proteins such as β-catenin and vimentin, cell cycle-related factors including Cyclin D1 and CDK4, as well as stemness-related factors including Oct-4 and Sox2 were all downregulated by silencing MYH9. Moreover, the MYH9 knockdown increased the expression of γ-H2AX, suggesting the enhancement of double-strand DNA breaks (Supplemental Fig. 1D).

HMGA1 promotes malignant phenotypes and chemoresistance of glioma cells via MYH9
To confirm whether HMGA1 modulating of MYH9 was responsible for the promotion of the proliferation, migration, and invasion of glioma cells, we utilized an MYH9-overexpression vector in shHMGA1 glioma cells. MTT and Edu incorporation assays showed that MYH9 upregulation reversed the shHMGA1-mediated reduction of proliferation in glioma cells ( Fig. 4A-C). Transwell chamber invasion and Boyden chamber migration assays verified that the invasive and migratory abilities were increased in the shHMGA1 + MYH9 plasmid-treated cells group compared to the shHMGA1treated cells group ( Fig. 4D and E). To determine whether MYH9 reversed the shHMGA1-mediated reduction of TMZ-induced cell death, dose-response curves for TMZ showed the IC 50 was rescued in shHMGA1 + MYH9 plasmid-treated cells relative to shHMGA1-treated cells ( Fig. 4F and G).
Subsequently, we examined the expression of several proteins related to the cell cycle, EMT, stemness, and DNA double-strand breaks in shHMGA1-treated cells overexpressing MYH9, compared with shHMGA1-treated and mock cells separately. As shown in Fig. 4H and Supplemental Fig. 1E, MYH9 overexpression reversed the shHMGA1-mediated increase of GSK-3β and rescued the shHMGA1-mediated reduction of EMT-related proteins (β-catenin and vimentin), cell cycle-related factors (cyclin D1 and CDK4), stemness-related factors (Oct-4 and Sox2), and the biomarker of DNA double-strand breaks (γ-H2AX).

MYH9 interacts with GSK-3β and promotes its ubiquitination
Previously, we showed that MYH9 expression affected the protein expression levels of GSK-3β, which could explain the effect of HMGA1 on GSK-3β expression. Then, we sought to identify the relationship between MYH9 and GSK-3β. We first measured the effect of MYH9 on GSK-3β mRNA levels, finding that silencing MYH9 had no effect on GSK-3β mRNA levels in glioma cells (Fig. 5A). However, overexpression of MYH9 downregulated GSK-3β protein levels and vice versa (Fig. 5B), indicating that GSK-3β protein expression was regulated at the post-transcriptional level. DOMINE data sets (http://domine. utdallas.edu/) [23] predicted that MYH9 might interact with GSK-3β, so a co-IP assay was performed. As expected, GSK-3β and MYH9 proteins coprecipitated (Fig. 5C). Overexpression of MYH9 not only downregulated GSK-3β protein levels but also markedly increased ubiquitin activity, suggesting that MYH9 mediates GSK-3β protein levels via ubiquitination. To test this hypothesis, we constructed three different plasmids, each with a mutation changing the lysine residue at three potential ubiquitination sites (15 K, 27 K, and 36 K), then the plasmids with the GSK-3β ubiquitination sites mutant (mt) and wild-type (wt) were transfected into U87 cells. Co-IP showed that all three mutants reduced the ubiquitination of GSK-3β, but the 27 K and 36 K mutant plasmids led to more significant inhibition of ubiquitination, indicating that they are likely the in vivo ubiquitination sites in GSK-3β (Fig. 5D).
We constructed two plasmids with mutated potential interaction domains of MYH9 (Mut1: 75-777, Mut2: 841-1921). Exogenous co-IP indicated that both mutants upregulated GSK-3β protein expression and inhibited ubiquitin activity, but mutation of domain 1 (75-777) displayed a greater suppression of ubiquitin activity (Fig. 5E). Finally, we detected the colocalization of MYH9 and GSK-3β by immunofluorescence staining. MYH9 and GSK-3β were mainly located in the cytoplasm, and overexpression of MYH9 reduced the cytoplasmic expression of GSK-3β (Fig. 5F). It has been reported that inhibition of GSK-3β promotes the nuclear translocation of non-phosphorylated β-catenin [24]. We also determined whether the ubiquitinated degradation of GSK-3β protein induced by HMGA1/MYH9 signaling promoted the translocation of β-catenin into the nucleus. Immunofluorescent detection of MYH9 and β-catenin was performed on glioma cells overexpressing HMGA1; we observed increased cytoplasmic expression of MYH9 and total expression of β-catenin. Furthermore, we found increased nuclear β-catenin (Supplemental Fig. 3).
Elevated HMGA1 and MYH9 expression as an unfavorable factor in gliomas The oncogenic HMGA1/MYH9 signaling pathway, shown in Fig. 6A, was assessed using the CCGA and The Cancer Genome Atlas (TCGA) databases, as well as using a cohort of 86 prospectively collected primary glioma tissues from our department. Survival analysis based on CCGA and TCGA data suggested that patients with low expression of HMGA1 or low expression MYH9 had a longer survival time compared to patients with high expression (Fig. 6B-E). Moreover, patients with both low expression of HMGA1 and MYH9 had remarkably longer survival times compared to patients with both high expression levels (Fig. 6F). In primary glioma tissues from our department, both HMGA1 and MYH9 were upregulated in glioma samples compared to NB tissues (Table 1). We found that age, histologic type, and WHO grade were significantly correlated with both HMGA1 and MYH9 expression (Table 1).
Using Kaplan-Meier analysis and the log-rank test (  Fig. 6G and H). Moreover, patients with high expression of both HMGA1 and MYH9 had the shortest survival times, but patients with low expression of both HMGA1 and MYH9 had the longest survival times (Fig. 6I). Further multivariate analysis demonstrated Boyden (E) chamber assays verified that the ability of cell invasion and migration was recovered in shHMGA1+ MYH9 plasmid-treated cells group compared to shHMGA1-treated cells group. F, G Dose-response curves for TMZ showed the TMZ concentration causing IC 50 was rescued in shHMGA1+MYH9 plasmid-treated cells groups compared to shHMGA1-treated cells groups. H Western blot assay confirmed MYH9 upregulation rescued shHMGA1-mediated reduction of the EMT-related proteins, the cell cycle-related factors, and the stemness-related factors, and MYH9-overexpression could reverse the shHMGA1-mediated increase of GSK-3β.
that WHO grade, HMGA1, and MYH9 expression significantly correlated with patient survival (P = 0.01, P = 0.002, and P < 0.001, respectively, Table 2), indicating that these three factors should be treated as independent prognostic factors in patients with glioma.

DISCUSSION
The therapeutic challenges presented by gliomas are their high malignancy and infiltrative growth, frequently resulting in incomplete resection and short-term recurrence. These characteristics are closely related to the EMT of tumor cells. The EMT not only promotes the infiltration of tumor cells to the periphery and leads to the residual tumor tissue after surgical resection, but also enables the glioma cells to acquire stem cell-like potentials and chemotherapy resistance [25][26][27]. Therefore, targeting the EMT in glioma cells has always been an effective strategy.
HMGA1 belongs to the HMG protein family, which alters chromatin structure by direct recognition of the A/T-rich sequences in its target promoter [28]. It is upregulated in gliomas and correlates with patient survival [9]; moreover, knockdown of HMGA1 suppresses the stemness of glioma stem cells and increases their sensitivity to temozolomide (TMZ) [29,30]. In accordance with these studies, we have found that HMGA1 promoted glioma cell proliferation, invasion, and migration, suggesting that it is a key regulator of the EMT process [10]. Moreover, a new oncogenic ZEB2/miR-637/HMGA1 signaling axis targeting vimentin was identified. Based on this, we further confirmed that HMGA1 inhibition sensitizes glioma cells to TMZ. A subsequent study on the mechanism of HMGA1 showed that HMGA1 upregulated the expression of the PI3K/Akt downstream factors c-Jun and MYH9, while upregulating GSK-3β. These results confirm the relationships among HMGA1, c-Jun, MYH9, and GSK-3β, but further studies on their precise mechanisms are necessary.
MYH9 is a hexamer that is comprised of two heavy chains, two 17 kDa essential light chains, and two 20 kDa modulated light chains [31]. It contains an IQ domain and a myosin head-like domain, which is involved in several important functions, including cytokinesis, cell motility, and maintenance of cell shape [32,33]. Mutations in this gene have been associated with multiple diseases [12,34,35], including cancers [13,16,[36][37][38][39][40]. MYH9 has been reported to be overexpressed in gastric cancer [41], nonsmall cell lung cancer [14], colon cancer [40], and breast cancer [42], but was downregulated in squamous cell carcinomas [43]. Moreover, MYH9 promoted the invasion of breast cancer cells [16], and the EGF-dependent phosphorylation of MYH9 resulted in increased migration [44]. However, our understanding of the role of MYH9 in gliomas is limited. It has been reported that MYH9 is deregulated in glioblastoma (GBM) serum small extracellular vesicles and is associated with GBM progression and metastasis [45]. But another study reported that overexpression of MYH9 abrogated the migration-inhibiting effects of the endogenous aryl hydrocarbon receptor agonist, indicating that MYH9 is essential for glioma cell migration [19]. Therefore, the actual function of MYH9 in gliomas is worth exploring. In this study, we found that MYH9 was overexpressed in glioma samples, which promoted glioma cell proliferation, invasion, migration, and TMZ resistance. We further found that HMGA1 increased MYH9 expression through c-Jun-mediated transcription, and MYH9 exerted an oncogenic effect by interacting with GSK-3β, followed by its ubiquitination and degradation. These results indicate that MYH9 plays an essential role in gliomas and as such, presents a potential new therapeutic target.
GSK-3β, which belongs to the serine-threonine kinase family, is involved in embryonic development, cell growth, cell differentiation, apoptosis, and tumorigenesis [46]. In gliomas, GSK-3β is downregulated; its overexpression inhibits angiogenesis and tumor growth by upregulating the β-catenin and mTOR/ p70S6K1 pathways [47], and also restores cell differentiation [48]. Moreover, knockdown of GSK-3β induced c-Myc-dependent cell death [49] and promoted glioma stem cell proliferation and colony formation [50]. More interestingly, inhibition of GSK-3β promoted the nuclear translocation of the non-phosphorylated β-catenin, subsequently binding the transcription factors of the TCF/LEF family and activating downstream targets such as MYC, cyclin D1, and BIRC5 [24]. In this study, we found that HMGA1mediated expression of MYH9 inhibited the expression of GSK-3β and induced the nuclear translocation of non-phosphorylated β-catenin. Therefore, we propose that HMGA1/MYH9 signaling promotes the nuclear translocation of β-catenin via ubiquitinated degradation of GSK-3β.
Our study demonstrated that both HMGA1 and MYH9 promote glioma cell proliferation, invasion, migration, and TMZ resistance. HMGA1 regulated PI3K/Akt/c-Jun signaling, and c-Jun promoted the expression of MYH9. Moreover, MYH9 binds to GSK-3β and promotes its ubiquitinated degradation, resulting in the promotion of malignant phenotypes and chemotherapy-resistant of glioma cells. Therefore, HMGA1 and MYH9 play essential roles in gliomas, which should be considered as tumoral enhancers with significant value as unfavorable progression indicators for patients with glioma, and may serve as feasible therapeutic targets in the future.

MATERIALS AND METHODS Patients and tissue samples
A total of 86 paraffin-embedded glioma specimens and 15 paraffinembedded adjacent non-tumor specimens were collected from the Department of Neurosurgery at Nanfang Hospital. Tissues from patients who had received preoperative radiation, chemotherapy, or biotherapy were excluded. All patients in this study underwent maximum safe total resection (partly resected cases were excluded). The subjects were 57 males and 29 females, with ages ranging from 24 to 69 (mean, 48.54) years. Prior consent from patients and approval from the Ethics Committees of Nanfang Hospital for Nationalities were obtained for the use of these clinical materials for research purposes. All specimens had a confirmed pathological diagnosis and were classified according to the 2016 World Health Organization (WHO) criteria. The end date of the follow-up study was December 2020, and the median follow-up for overall survival was 19.64 (range, 9-32) months.

IHC and evaluation of staining
Paraffin sections (4 μm) were deparaffinized and dehydrated, then antigen retrieval was performed in citrate buffer for 3 minutes. Endogenous peroxidase activity and nonspecific antigens were blocked with 3% H 2 O 2 and goat serum, followed by incubation with HMGA1 antibody (1:200, Abcam) or MYH9 antibody (1:250, Abcam) overnight at 4°C. After washing, sections were incubated with HRP-conjugated secondary antibodies and visualized using a DAB substrate (Maixin Biotech. Co., Ltd, Fuzhou, China). Immunostaining results were evaluated and scored by two independent pathologists blinded to the clinical parameters, as described in our previous work [10]. Positive nuclear staining scores were defined as: 0, <20%; 1, 20-49%; 2, 50-79%; and 3, >80%. Sums of the cytoplasmic and nuclear staining scores were used as final staining scores (0-12). For statistical analysis, a final staining score of 0-4 and 5-6 in the cytoplasm or 0-3 and 4-6 in the nucleus were considered to be below or high expression, respectively.

Cell background and transfection
U87 and LN229 were purchased from the Chinese Academy of Sciences (Shanghai, China). U87 and LN229 cells are both MGMT-promoter methylated [51], which is currently the only known biomarker for TMZ response in GBM patients, but these two GBM cell lines are highly resistant to TMZ treatment, indicating another factor rather than MGMT conferred the therapeutic resistance in U87 and LN299 cells. Moreover, the expression of HMGA1 and MYH9 were both high in U87 and LN229 cell lines as shown in Fig. 1J and Fig. 3I. Plasmids and siRNAs were designed and synthesized by RiboBio Co., Ltd. (Guangzhou, China). Lentiviral vectors

RT-PCR and qPCR
Total RNA was isolated from cells or fresh clinical tissues using Trizol (TaKaRa Bio, Inc., Shiga, Japan), and cDNA synthesis was performed using reverse transcription reagents (TaKaRa Bio, Inc., Shiga, Japan). Then cDNA was used as a template for amplification with specific primers (Supplemental Table 1). The specificity of the amplification products was confirmed by melting curve analysis. PCR for each gene was repeated thrice. Independent experiments were performed in triplicates.

Western blotting
Cell-lysate protein concentrations were determined using a BCA protein assay kit (Thermo Scientific, Waltham, MA, USA

Migration and invasion assay
A Transwell assay (BD Biosciences, NJ, USA) was performed to detect cell migration and invasion abilities. Cells were suspended in 100 μl Dulbecco's Modified Eagle Medium (DMEM) without serum and seeded into the top chamber of the Transwell coated with Matrigel (BD Biosciences, NJ, USA) or left uncoated, and the bottom chambers were filled with 500 μl DMEM supplemented with 10% fetal bovine serum. The migrated cells were stained with crystal violet, photographed, and quantified by counting the cell number in five random fields.

MTT and EdU incorporation assay
Cell proliferation and drug sensitivity were tested using the MTT assay. Exponentially growing cells were seeded into 96-well plates and cultured overnight to allow for cell adherence. Cell viability was measured using MTT (5 mg/ml) (Sigma-Aldrich, MO, USA). EdU incorporation was performed using an Apollo567 in vitro imaging kit (RiboBio Co., Ltd, Guangzhou, China) according to the manufacturer's protocol. Cells were incubated with 10 μM EdU for 2 h, then fixed with 4% paraformaldehyde.

TMZ chemosensitivity
Transfected U87 and LN229 cells were seeded at a density of 6000 cells per well, and after 24 h the cells were treated with different concentrations (50, 100, 200, and 400 μM) of TMZ (Merck product) for 48 h. MTT assays were then performed to determine the fraction of cells that survived TMZ exposure. The resistance of the induced cells was measured using the IC50. The TMZ concentration causing 50% inhibition of glioma cell activity was defined as the IC50.

In vivo orthotopic xenograft study
Animal protocols were approved by the Institutional Animal Ethics Committee of the Experimental Animal Center of Nanfang Hospital. Female nude mice (4 weeks old) were purchased from Southern Medical University and divided into two groups of five mice per group. After acclimation for 1 week, the mice were anesthetized with 1% pentobarbital sodium and placed in a stereotactic fixation device (RWD Life Science, Shenzhen, China) prior to xenograft inoculation. A total of 10 6 cells in 5 µl of sterile PBS were stereotactically implanted into the point located 1 mm posterior to the bregma, 2 mm left of the midline, and 2.5 mm deep into the brain. Cells were automatically injected by a stereotactic fixation device over 5 minutes. The mice were sacrificed when the tumors of control groups were grown into the appropriate size confirming by animal MRI. The mice brains were then fixed in 4% paraformaldehyde for 24 hours and embedded in paraffin for further staining.

Dual-luciferase reporter assay
A fragment of the MYH9 promoter was amplified by PCR. Mutagenesis of the c-Jun binding site was performed using the GeneTailor Site-Directed Mutagenesis System (Invitrogen, Guangzhou, China). The wild-type and mutant promoters were separately cloned into the psiCHECK-2 vector. The construct was cotransfected with si-c-Jun and c-Jun plasmids (or their control plasmids) into cells. Luciferase activity was measured 48 h after transfection using the Dual-Luciferase Reporter Assay System (Promega Corporation, Madison, WI, USA).

Chromatin immunoprecipitation
ChIP assays were performed using a commercial assay kit (Thermo Scientific, Waltham, MA, USA) according to the manufacturer's protocol.
Chromatin was cross-linked, isolated, and digested with micrococcal nuclease to obtain DNA fragments. Anti-c-Jun (ChIP Grade, Abcam) or IgG was added for immunoprecipitation. After elution and purification, the recovered DNA fragments were subjected to RT-PCR and qPCR. IgG served as a negative control.

Co-immunoprecipitation assay
Co-IP was carried out using a Pierce Co-Immunoprecipitation kit (Thermo Scientific, Waltham, MA, USA) according to the manufacturer's instructions. Briefly, total protein was extracted from cells and quantified. A total of 5 mg of protein was incubated with 10 μg specific antibodies (Abcam) or IgG overnight at 4°C. After elution, the recovered proteins were subjected to silver staining, mass spectrometry, and western blotting. IgG was used as a negative control.

Statistical analysis
All experiments were replicated at least three times. All data were analyzed using SPSS (version 26.0; SPSS Inc. Chicago, IL, USA) and GraphPad Prism v9.0 (GraphPad Software, San Diego, CA, USA) software. Data are expressed as means ± SD from at least three independent experiments. Statistical significance was determined using the Student's two-tailed t test for two groups, one-way ANOVA for multiple groups, and a general linear model for tumor growth and MTT assays. Correlations between gene expression and clinicopathological characteristics were analyzed using the chi-square test. Log-rank tests were performed on Kaplan-Meier survival curves to elucidate relationships between gene expression and overall survival in patients. Univariate and multivariate survival analyses were performed using the Cox proportional hazards regression model. All statistical tests were two-sided. A P value of <0.05 was considered to be statistically significant.

DATA AVAILABILITY
The data sets used and/or analyzed during the current study are available from the corresponding author on reasonable request.