- Open Access
ZNF326 promotes malignant phenotype of glioma by up-regulating HDAC7 expression and activating Wnt pathway
Journal of Experimental & Clinical Cancer Researchvolume 38, Article number: 40 (2019)
Zinc-finger protein-326 (ZNF326) was initially found in the NIH3T3 cell line to regulate cell growth, however, the expression and underlying role of ZNF326 in human tumours, especially in glioma, is not fully understood.
Immunohistochemistry was applied to detect the expression of ZNF326 in glioma tissues, and statistical analysis was used to analyse the relationship between ZNF326 expression and clinicopathological factors. The effect of ZNF326 on glioma cells proliferation and invasion was conducted by functional experiments both in vivo and in vitro. Chromatin immunoprecipitation and dual-luciferase assays were performed to demonstrate that histone deacetylase enzyme-7 (HDAC7) is the target gene of ZNF326. Immunoblotting, real-time PCR, GST-pulldown and co-immunoprecipitation assays were used to clarify the underlying role of ZNF326 on Wnt pathway activation.
High nuclear expression of ZNF326 was observed in glioma cell lines and tissues, and closely related with advanced tumour grade in the patients. Moreover, ectopic ZNF326 expression promoted the proliferation and invasiveness of glioma cells. Mechanistically, ZNF326 could activate HDAC7 transcription by binding to a specific promoter region via its transcriptional activation domain and zinc-finger structures. The interaction of the up-regulated HDAC7 with β-catenin led to a decrease in β-catenin acetylation level at Lys-49, followed by a decrease in β-catenin phosphorylation level at Ser-45. These changes in β-catenin posttranscriptional modification levels promoted its redistribution and import into the nucleus. Additionally, ZNF326 directly associated with β-catenin in the nucleus, and enhanced the binding of β-catenin to TCF-4, serving as a co-activator in stimulating Wnt pathway.
Our findings elucidated ZNF326 promotes the malignant phenotype of human glioma via ZNF326-HDAC7-β-catenin signalling. This study reveals the vital role and mechanism of ZNF326 in the malignant progression of glioma, and provides the reference for finding biomarkers and therapeutic targets for glioma.
Glioma is the most common intracranial tumour originating from the epithelial cells of the central nervous system, accounting for more than 60% of primary brain tumours [1,2,3]. Current therapeutic strategies for glioma consist of neurosurgical resection, chemotherapy, and radiotherapy. However, all these strategies have failed to yield an expected good prognosis of malignant glioma. This could be due to the highly aggressive nature of glioma cells that are capable of infiltrating into the adjacent normal brain tissue [4, 5]. Therefore, there is always an urgency to develop novel strategies for timely diagnosis and therapeutic agents for patients with glioma.
The Wnt pathway elicits an important regulatory signal that is able to influence embryonic development of different tissues and organs, including the nervous system. Dysregulation of the Wnt signal is implicated in the development and progression of glioma [6, 7]. Indeed, several studies have revealed that the expression and nuclear localisation of β-catenin and its transcription factor TCF4 are significantly higher in glioma compared to that in normal brain tissue, and these characteristics positively correlate with the glioma grade . Similarly elevated expression of two Wnt pathway activators, TCF4 and SOX, were also reported [9, 10]. Moreover, high levels of canonical Wnt factors such as DKK1, FZD1, and LEF1 were found to be associated with very poor clinical outcome . Oncogenic activities such as proliferation, inhibition of apoptosis, and invasion have also been coupled to abnormal Wnt/β-catenin signalling in glioma cell lines [12, 13]. Altogether, these findings indicate that the Wnt pathway plays a fundamental role in gliomagenesis.
In the canonical Wnt pathway, β-catenin, which is a central player of this signalling cascade, is entrapped into a protein complex formed by Axin1, glycogen synthase kinase-3β (GSK-3β), casein kinase 1α (CK1α), and adenomatous polyposis coli (APC). This complex favours β-catenin degradation by proteasomal ubiquitination. However, when a Wnt protein binds to the membrane receptors of the frizzled (FZD) and low-density lipoprotein receptor-related protein (LRP5/6) families, the degradation complex appears to be inhibited. Thus, unphosphorylated β-catenin accumulates, enters the nucleus and binds to TCF-4, and activates expression of Wnt signalling pathway target genes such as AXIN2, C-MYC, CCND1, and MMP7 [14,15,16,17,18].
Zinc-finger protein-326 (ZNF326) was first identified in NIH3T3 cells and is believed to play an important role in neuronal differentiation . Although the molecular mechanism of ZNF326 is not yet completely understood, it is essentially a protein molecule of 582 amino acids, with C2H2 zinc-finger domain, and acts as a potential transcription factor. The main functional domains include: a transcriptional activation domain (TAD) near the N-terminus (1-124aa), an intra-nuclear localisation sequence between 242-260aa (NLS), and a central region containing two C2H2 zinc-finger domains (313-336aa and 407-430aa) . Until date, the expression of ZNF326 in human glioma, its effect on the malignant phenotype of glioma cells, and the possible signal transduction pathway involved have not been reported.
In this study, we report the clinical relevance of ZNF326 in glioma and its regulatory effect on the Wnt/β-catenin signalling pathway. Initially we measured the expression level of ZNF326 in human resected glioma specimens, and analysed the relationship between ZNF326 expression and clinicopathological factors of glioma. We also investigated the effects of ZNF326 on the proliferation and invasiveness of glioma cells. At the molecular level, a series of ZNF326 mutants were constructed and studied. ChIP-seq, luciferase reporter assay, and co-immunoprecipitation were applied to detect the effects of ZNF326 on the transcriptional activity of target genes. The results here provide a theoretical and experimental basis for glioma treatment strategies.
Human glioma tissue collection and ethics statements
A total of 133 tumour specimens (GradeI-IV) were obtained from patients (average age: 50 years) during surgery at the First Affiliated Hospital of China Medical University from 2008 to 2017. All patients were chemotherapy- and radiotherapy-nave prior to this surgical resection. The clinical data such as age, gender, tumour location and WHO grade were recorded for statistical analysis. All clinical investigations were conducted according to the principles expressed in the Declaration of Helsinki. Written informed consent was obtained from all patients, and all procedures were approved by the Institute Research Ethics Committee of China Medical University.
Assays were performed as described previously . Briefly, the tissue sections were incubated with anti-ZNF326 and HDAC7 mouse monoclonal antibody (sc-390,606 and sc-74,563, respectively 1:50 and 1:20, Santa Cruz Biotechnology Inc. CA, USA). The intensity of ZNF326 and HDAC7 staining was scored as follows: 0 (no staining and weak), 1 (moderate), or 2 (strong). Percentage scores were assigned as follows: 1 (1–25%), 2 (26–50%), 3 (51–75%), and 4 (76–100%). The scores of each tumour sample were multiplied to give a final score of 0 to 8. According to the staining intensity and the staining extent scores, the immunohistochemistry (IHC) result was classified as follows: 0–3, negative (−); 4–5, weakly positive (+); 6–7, moderately positive (++) and ≥ 7, strongly positive (+++). PBS and goat serum were used as negative controls.
Cell lines and cell culture
Human glioma cell lines, including U251, U87, SHG44 and U118 were purchased from the American Type Culture Collection (Manassas, VA, USA). All culture media were supplemented with 10% foetal bovine serum (FBS, Hyclone, Logan, UT, USA). Normal human astrocyte (NHA) is a gift supplied by Professor Anhua Wu (China Medical University). U87 and NHA was cultured in Eagle’s Minimum Essential Medium (MEM, Hyclone), while U251, U118, and SHG44 cells were all cultured in Dulbecco’s modified Eagle’s medium (DMEM, Hyclone). Cells were maintained at 37 °C in a humidified atmosphere of 5% CO2. All cell lines were authenticated by short tandem repeat (STR) DNA profiling.
Plasmids and reagents
pCMV6 empty vector (#PS-100001) and Myc/DDK-pCMV6-ZNF326 (#RC-210339) were purchased from Origene (Rockville, MD, USA). pcDNA3.1 empty vector (#52535), pcDNA3.1-FLAG-HDAC7 (#13824), TCF4 plasmid (#16512), pEGFP-N1 empty vector (#86776), pEGFP-N1-β-catenin (#71367) and Super8 × TOPflash (#12456) were purchased from Addgene (Cambridge, MA, USA). pRL-TK (#E2241) was purchased from Promega (Madison, Wisconsin, USA). Control siRNA (sc-37,007), siRNA-ZNF326 (sc-88,338) and siRNA-HDAC7 (sc-35,546) were purchased from Santa Cruz Technology. The nucleotide sequence shRNA-ZNF326 was provided by Dr. Roberto Rangel and Professor Nancy A. Jenkins at the Anderson Cancer Center of the US. ShRNA-ZNF326, shRNA-HDAC7 plasmid and the lentivirus envelope shZNF326/ZNF326 were constructed by Genechem company (Shanghai, China). The mutants ZNF326-△TAD, ZNF326-△Zn1, ZNF326-△Zn2 and ZNF326-△TAD&△Zn1 + 2 were also constructed by Genechem. HA-CBP plasmid is a gift from Prof. Liu Cao (Department of Translational Medicine, China Medical University). Lipofectamine 3000 (Invitrogen, Carlsbad, CA, USA) transfection reagent was used for plasmid transfection. Puromycin (Sigma-Aldrich, St. Louis, MO, USA) was used to select stably transfected cells.
Cell extraction and immunoblotting
Assays were performed as described previously . Trichostatin A (TSA) was purchased from Sigma-Aldrich. Expression was quantified using densitometry and ImageJ software. Detailed information about primary antibody is provided in the Additional file 1.
Cells were fixed with 4% paraformaldehyde, blocked with 3% BSA, and incubated with anti-β-catenin antibody (1:100, BD Biosciences, #610153) overnight at 4 °C, followed by incubation with fluorescein isothiocyanate (FITC)-conjugated secondary antibodies, at room temperature for 1 h. Cells were counterstained with 4′,6-diamidino-2-phenylindole (DAPI). Confocal microscopy was performed using a Radiance 2000 laser scanning confocal microscope (Carl Zeiss, Thornwood, NY, USA).
GST (glutathione S transferase)-conjugated β-catenin protein was expressed in E. coli BL21 cells and was purified using standard protocols. Glutathione-Sepharose beads (GE Healthcare, Waukesha, WI, USA) coupled with either GST or with the GST-β-catenin purified protein were incubated with the HEK293 cells lysates transfected with Myc/DDK-pCMV6-ZNF326 plasmid overnight at 4 °C. Complexes were washed and subjected to immunoblotting and coomassie blue staining.
Chromatin immunoprecipitation and sequencing
ZNF326-overexpressing H1299 cells were crosslinked and lysed by ultrasound treatment. The pyrolysis liquid was divided into four groups and treated with different antibodies. Protein G agarose was added at 4 °C, followed by rotation for 1 h. The corresponding antibody was then added, followed by incubation overnight at 4 °C. After protein/DNA complex elution, the complex was disintegrated. After recovery of DNA samples, high-throughput sequencing was performed (#17–371; Millipore, Burlington, MA, USA).
ChIP experiment was carried out according to the procedure described by manufacturer (Millipore, MA, USA). The immunoprecipitated DNAs were amplified by qPCR. The primers used are listed in Additional file 1.
Assays were performed as described previously . Briefly, recombinant human Wnt3a (#5036-WN, R&D Systems, France) was reconstituted in PBS containing 0.2% BSA to a concentration of 10 μg/mL and used in experiments at a final concentration of 100 ng/mL. Luciferase-reporter plasmids are described in detail in Additional file 1.
RNA extraction and real-time RT-PCR
Colony formation, Matrigel invasion, and MTT assays
Glioma cells were seeded in 6-cm cell culture dish (1000 per dish) and incubated for 14 days. The plates were then washed with PBS and stained with Giemsa before counting the number of colonies consisting of > 50 cells.
Matrigel invasion and MTT assays were used in this study. All experiments were performed in triplicate. The detailed protocol is provided in Additional file 1.
Transplantation of tumour cells into nude mice
The nude mice used in this study were treated by following the experimental animal ethics guidelines issued at China Medical University. The study was approved by Institutional Animal Research Committee of China Medical University. The nude mice (BALB/c, SPF grade, 16–18 g, 4 weeks old, and female) were purchased from Charles River (Beijing, China), and the axilla of each mouse was subcutaneously inoculated with 5 × 106 glioma cells in 0.2 mL of sterile PBS. Six weeks after inoculation, the mice were euthanised and autopsied to examine tumour growth. A portion of explant tumour was fixed in 4% formaldehyde and embedded in paraffin. Serial 4-μm-thick sections were cut and stained by IHC, and the stained sections were examined under a microscope. The proliferation rate was evaluated by counting Ki-67-positive nuclei in more than 30 high-power fields (HPFs) per group.
The statistical software SPSS 22.0 (SPSS, Chicago, IL, USA) was used for all analyses. The chi-squared test was used to assess correlations between ZNF326 expression and clinicopathological factors. Differences between the groups were tested by Student’s t-test. The correlation between ZNF326 and HDAC7 expression in glioma specimens was tested by Pearson-correlation test. A p value of < 0.05 was considered to indicate statistically significant differences.
ZNF326 is highly expressed in glioma and positively correlated with tumour grade
To explore a potential role of ZNF326 in glioma tumorigenesis, we performed IHC in a cohort of 133 human patients with glioma samples to examine the expression profiles of ZNF326. We found that nearly 60.9% (81/133) of patients with glioma had high level of nuclear ZNF326 (+, ++ and +++) in the glioma samples. In addition, ZNF326 was negative in 5 cases of glioma tissues with grade-I and the ZNF326 staining was significantly associated with tumour grade (Fig. 1a, and c, P = 0.000) and age (P = 0.012), but not with gender and tumour location (Table 1). Consistently, the TCGA database (http://gepia.cancer-pku.cn/) indicated that ZNF326 mRNA levels in glioma are significantly higher than those in normal brain (Fig. 1b, P < 0.05). Similar results were obtained in a panel of four glioma cell lines in vitro, compared with that in normal human astrocyte cell line NHA (Fig. 1d). Immunofluorescence assay indicated ZND326 located in the nucleus of glioma cells (Fig. 1e). Altogether, these results suggest that ZNF326 could be used as a potential predictor of malignancy in gliomas.
Interestingly, we found that the expression of HDAC7 in NHA cell lines was also significantly lower than that in other glioma cell lines, similar to ZNF326 expression (Fig. 1d). Consistently, compared with the ZNF326 IHC staining results, the expression of HDAC7 was positively correlated with the expression of ZNF326 in gliomas (Additional file 2: Figure S1, Additional file 3: Table S1).
ZNF326 promotes proliferation and invasion of glioma cells in vitro
To verify whether ZNF326 has a causal role in regulating glioma cell phenotypes, we stably overexpressed ZNF326, using a lentivirus vector–based ZNF326 plasmid, in U87 and U251 cell lines. Compared with that in the control group, ZNF326 overexpression enhanced the following in U87 and U251 cells: colony formation (U87: CTL vs ZNF326, 100 ± 11 vs 190 ± 10, P < 0.01; U251: CTL vs ZNF326, 26 ± 3 vs 48 ± 2, P < 0.01; Fig. 2a-b), invasiveness (U87: CTL vs ZNF326, 75 ± 5 vs 132 ± 6, P < 0.01; U251: CTL vs ZNF326, 51 ± 2 vs 87 ± 5, P < 0.001; Fig. 2e-f). and proliferation (MTT assay; U87: CTL vs ZNF326, 0.904 ± 0.035 vs 1.254 ± 0.062, P < 0.01; U251: CTL vs ZNF326, 0.804 ± 0.049 vs 1.194 ± 0.032, P < 0.001; Fig. 2i-j), Conversely, ZNF326 knockdown in U87 and U251 cells by a lentivirus vector–based ZNF326 shRNA technique significantly weakened the following in both cell lines: colony formation (U87: CTL vs siZNF326, 89 ± 7 vs 57 ± 6, P < 0.05; U251: CTL vs siZNF326, 43 ± 2 vs 25 ± 2, P < 0.01, Fig. 2c-d), invasiveness (U87: CTL vs siZNF326, 63 ± 5 vs 34 ± 6, P < 0.05; U251: CTL vs siZNF326, 41 ± 2 vs 23 ± 3, P < 0.01; Fig. 2g-h) and proliferation (U87: CTL vs siZNF326, 0.808 ± 0.026 vs 0.528 ± 0.028, p < 0.001; U251: CTL vs siZNF326, 1.028 ± 0.050 vs 0.786 ± 0.061, p < 0.05; Fig. 2k-l).
ZNF326 positively regulates the Wnt signalling pathway
The close association that has been reported between Wnt/β-catenin signalling and glioma tumorigenesis, combined with our KEGG analysis, predicted that ZNF326 is closely related to the Wnt/β-catenin pathway (Additional file 2: Figure S2-A). To verify this prediction, we first assessed the effect of ZNF326 on the activities of the Wnt pathway in glioma cell lines, using luciferase reporter assays. To readily detect differences, we pre-stimulated the Wnt signalling pathway using Wnt-3a . In U87 and U251 cells, ectopic ZNF326 expression significantly increased the TOPflash activity of Wnt signalling induced by Wnt3a (Fig. 3a), and the Wnt pathway activity was gradually increased in a dose-dependent manner upon ZNF326 transfection in HEK293 cells (Fig. 3b).
ZNF326 overexpression in the U87 cell line significantly upregulated MMP7, C-myc and AXIN2 at both the protein and mRNA levels. Conversely, ZNF326 knockdown in U87 cells downregulated the expression of Wnt-related genes (Fig. 3c, d). This is consistent with the positive correlation between ZNF326 and Wnt target genes obtained by online analysis at the GEPIA website (http://gepia.cancer-pku.cn/, Additional file 2: Figure S2-B). In particular, we have noticed that the changes of ZNF326 expression have no significant effect on total β-catenin expression level. Notably, endogenous interaction between ZNF326 and β-catenin was detected in the U87 cell line (Fig. 3e) and in vitro glutathione S-transferase pull-down assays (Fig. 3f) confirmed their direct interaction. We next tested whether ZNF326 can promote TCF4–β-catenin interaction. Immunoprecipitation indicated that ZNF326 overexpression increased the interaction between β-catenin and TCF4 (Fig. 3g), whereas ZNF326 silencing decreased this interaction (Fig. 3).
ZNF326 regulates target-gene HDAC7 expression through its transcriptional activation domain and zinc-finger structure
Interestingly, ZNF326 overexpression in HEK293 cells promotes β-catenin nuclear import, as shown using cytosolic fractionation assay and immunofluorescence (Fig. 4a-b). Next we investigated the underlying mechanism in this process. We transfected ZNF326 plasmid with MYC tag into the H1299 cell line with high transfection efficiency and searched for the potential target genes of ZNF326, using ChIP-seq assay. From the Broad-Peak model analysis, eight potential target genes were found to be associated with these promoter regions , including HDAC7, which was directly related to Wnt pathway activation [23, 24].
To study the phenomenon further, we designed eight pairs of qPCR primers for the HDAC7 promoter region (0 to 2000 bp), and then transfected U87 cells with the Myc-tagged ZNF326 plasmid. ChIP-qPCR assay indicated that ZNF326, but not the control protein IgG, was bound to the promoter region of HDAC7, corresponding to primers no. 3 (− 1552 bp to − 1301 bp) and primers no. 5 (− 1073 bp to − 788 bp) (Fig. 4c). The HDAC7 promoter was examined to identify the region critical to its activity and responsiveness to ZNF326 expression. Truncated reporters were constructed as shown in Fig. 4d. Luciferase reporter assay was performed to show that the ZNF326 regulated HDAC7 transcription activity by associating with the − 1552 to − 1301-bp and − 1073 to − 788-bp regions of the HDAC7 promoter (Fig. 4e). Additionally, we explored the domains of ZNF326 that contribute to the binding with the HDAC7 promoter. We designed a series of ZNF326 mutant plasmids, including the deletion of TAD region, single deletion of a zinc-finger structure, and complete deletion of both zinc-finger structures (Fig. 4f). HEK239 cells were then co-transfected with HDAC7 promoter-reporter construct and these ZNF326 mutant plasmids. The results indicated that upon deletion of TAD and complete deletion of both zinc-finger structures, the mutant ZNF326 could not activate the HDAC7 promoter-luciferase reporter activity (Fig. 4g). Consistently, RT-qPCR (Fig. 4h) and western blotting (Fig. 4j) also demonstrated that the deletion of the TAD region and complete deletion of the two zinc-finger structures in ZNF326 prevented the increase in mRNA and protein levels of HDAC7. Similarly, HDAC7 expression was downregulated by ZNF326 knockdown (Fig. 4i, k). Furthermore, the online analysis at the GEPIA website (http://gepia.cancer-pku.cn/) revealed that ZNF326 was positively correlated with HDAC7 expression (Additional file 2: Figure S3).
HDAC7 deacetylates β-catenin at Lys49 and promotes β-catenin nuclear import
We further examined the role of HDAC7 in the Wnt pathway. We found that HDAC7 knockdown in the U87 cell line, using siRNA-HDAC7, significantly decreased Wnt target genes, C-myc and MMP7 expression, which was in accordance with the positive correlation between HDAC7 and some other Wnt target genes (such as CD44, PTGS2, and MMP2) revealed by GEPIA online analysis (Additional file 2: Figure S4). Meanwhile, the total acetylation level of β-catenin was also upregulated (Fig. 5a). Conversely, HDAC7 ectopic expression in HEK293 cells increased the expression of Wnt-related genes and decreased β-catenin acetylation level (Fig. 5b). Notably, HDAC7 overexpression promoted the nuclear import of β-catenin (Fig. 5c-d).
Based on the experimental data above, we hypothesised that HDAC7 regulates the level of β-catenin nuclear import by affecting its acetylation level. To test this, immunoprecipitation assay in U87 and U251 cells was conducted to demonstrate that HDAC7 could interact with β-catenin (Fig. 5e-f). According to the literature that states that HDAC6 can deacetylate β-catenin at Lys49 , we tested whether HDAC7 also has a similar effect. We thus co-transfected HDAC7 and HA-CBP (CREB binding protein, a protein that can increase β-catenin acetylation at Lys-49) plasmids in the HEK293 cell line . The results showed that HDAC7 could also reverse the up-regulation of β-catenin acetylation at Lys49 induced by CBP (Fig. 5g, h), whereas HDAC7 silencing in U87 cells significantly increased β-catenin acetylation at Lys49 (Fig. 5i). CK1α-mediated phosphorylation of β-catenin at the serine 45 residue is a key step for β-catenin degradation . Our results showed that HDAC7 knockdown significantly upregulated the level of phosphorylation of β-catenin Ser45 (Fig. 5i). Accordingly, the levels of β-catenin acetylation at Lys49 and β-catenin phosphorylation at Ser45 were significantly decreased in a dose-dependent manner after HDAC7 overexpression in HEK293 cells (Fig. 5j).
ZNF 326 promotes β-catenin nuclear import by increasing HDAC7 expression
To further understand HDAC7 regulation by ZNF326, siRNA-HDAC7 or TSA were used to knock out HDAC7 expression or inhibit HDAC7 activity, respectively. The proliferation and invasiveness of glioma cells were reversed by HDAC7 silencing or HDAC7 inhibition, as confirmed by MTT assay and Transwell assay, respectively (Fig. 6a-d, Additional file 2: Figure S5). On the other hand, immunoprecipitation assay demonstrated that ZNF326 overexpression in HEK293 clearly decreased the total β-catenin acetylation level in a dose-dependent manner (Fig. 6e). Overexpression of ZNF326 in U87 cells also deacetylated β-catenin at Lys49 (Fig. 6f). Furthermore, we transfected ZNF326 wild-type and ZNF326 mutant plasmid lacking both the zinc-finger structures and the TAD domain in HEK293. The results showed that wild-type ZNF326 could increase the level of HDAC7 protein, and reduce the level of β-catenin acetylation at Lys49 and phosphorylation at Ser45, and that the mutant abrogated this effect (Fig. 6g). Finally, after overexpression of mutant ZNF326 or knocking down HDAC7 or adding TSA in U87 and U251 cells, the effect of ZNF326 on the reduction of β-catenin acetylation at Lys49 and phosphorylation at Ser45 disappeared (Fig. 6h, i).
ZNF326 promotes tumour growth in xenograft model of nude mice in vivo
To verify the impact of ZNF326 on the tumour growth of glioma cells in vivo, we evaluated the role of ZNF326 in tumour formation of U87 and U251 cells, using a xenograft model of nude mice. As shown in Fig. 7a-c, stable expression of ZNF326 using lentivirus in U87 cells (selected with puromycin, 5 μg/mL) significantly promoted tumour growth in vivo when compared with that in the control group (CTL vs ZNF326, volume:0.168 ± 0.035 vs 0.400 ± 0.046, P < 0.01; weight:0.136 ± 0.024 vs 0.378 ± 0.056, P < 0.01). Consistent with this observation, Ki-67 staining, and transcription of Wnt target genes (AXIN2, CCND1, and MMP7) and HDAC7 were markedly enhanced in tumours in which ZNF326 was overexpressed (Fig. 7d-e). In addition, we used the lentivirus-shZNF326 to knock down ZNF326 (selected with puromycin, 10 μg/mL) in the U251 cell line. The results showed that, the volume and weight of xenografts in nude mice after ZNF326 knockdown significantly decreased (CTL vs shZNF326, volume: 0.522 + 0.126 vs 0.013 + 0.004, P < 0.01; weight: 0.620 + 0.097 vs 0.101 + 0.012, P < 0.001, Fig. 7f-h); and the Ki-67 staining and mRNA levels of Wnt pathway target genes and HDAC7 decreased (Fig. 7i-j). These results suggested that ZNF326 promotes tumor growth by up-regulating the expression of HDAC7 and activating Wnt pathway.
To confirm the hypothesis above, we attempted to compare the effects of ZNF326 and ZNF326 + shRNA-HDAC7 on tumour growth. We found the volume and weight of xenografts in ZNF326 + shRNA-HDAC7 group markedly reduced compared to the ZNF326 group (ZNF326 vs ZNF326 + shRNA-HDAC7, volume: 1.180 + 0.141 vs 0.736 + 0.052, P < 0.05; weight: 0.763 + 0.060 vs 0.366 + 0.041, P < 0.05), indicating that ZNF326 promotes tumour growth, at least partially through HDAC7 (Fig. 8a-c). In particular, we noticed transfection of ZNF326 or shRNA-HDAC7 had no significant effect on the total amount of β-catenin in xenografts tissue (Fig. 8d), which was consistent with the experimental results in vitro.
The zinc-finger protein ZNF326 was initially found by Lee et al. in the NIH3T3 cell line, and detected to be highly expressed in the brain and neural tube of the E11.5 embryo, which indicated that ZNF326 might play an important role in the process of nerve development [19, 20]. Recently, it was reported that ZNF326 may form a DBIRD complex with DBC1in human embryonic kidney cell line HEK293 through its S1RNA domain. This complex acts on mRNA particles and RNAII (polymerase), which regulates gene transcription and alternative splicing .These led us to explore the underlying role of ZNF326 in the development and progression of glioma. In this study, we uncovered that ZNF326 expression is upregulated in glioma specimens, which is also consistent with the result obtained from the TCGA database. Statistical analysis indicated high ZNF326 expression positively correlated with tumour grade, which hinted that ZNF326 might function as an oncogene in glioma. In vitro and in vivo experiments verified that ZNF326 can promote glioma cell proliferation and invasiveness and tumour xenograft formation in nude mice. These results were sufficient to conclude that ZNF326 is a tumour-promoting factor in glioma.
So far, the underlying mechanism of ZNF326 on gliomagenesis is almost unknown. Mounting evidence shows that overstimulation of Wnt signalling can lead to malignant progression of glioma. We found ZNF326 transfection in glioma cells significantly increased the activity of the Wnt signalling pathway and its target gene expression, and there was an opposite effect upon ZNF326 knockdown. This is consistent with the positive correlation between ZNF326 and Wnt target genes as obtained from our GEPIA website (http://gepia.cancer-pku.cn/) and gene enrichment analysis. After fully determining the relationship between ZNF326 and the Wnt signalling pathway, we focused on the role of ZNF326, which is multifaceted. On one hand, ZNF326 can directly bind with β-catenin in the nucleus. This combination can further promote the binding of β-catenin and transcription factor TCF4, thus activating the transcription of target genes of the Wnt pathway. ZNF326 functions as a co-activator in this process. Unfortunately, we are not clear about the molecular structural basis of the combination of ZNF326 and β-catenin, and will explore it in our future research.
On the other hand, ZNF326 can also play a vital role as a transcription factor. We used ChIP-seq and ChIP-qPCR to screen and prove that HDAC7 is one of the target genes downstream to ZNF326. Moreover, ZNF326 transfection can significantly up-regulate the transcriptional level of HDAC7. Subsequently, we further designed and constructed a series of ZNF326 mutant plasmids, and proved that wild-type ZNF326 can bind to the HDAC7 promoter (− 1552 bp to − 1301 bp and − 1073 to − 788 bp) through its TAD and zinc-finger structures, whereas the mutant ZNF326, which cannot bind to the HDAC7 promoter region, abrogates this effect. Thus, it was further confirmed that ZNF326 acts as a transcription factor to positively regulate the transcriptional activity of target gene HDAC7.
HDAC7 is a member of the HDAC family, and plays an important role in many biological processes and is closely related to the occurrence of cancer. Li et al reported that inactivation of HDAC6 increases the level of β-catenin acetylation at Lys49, which leads to increased EGF-induced β-catenin nucleation, and promotes the malignant phenotype of colon cancer cells . We found that after transfection of ZNF326, the total β-catenin expression level was not significantly altered; so we focused our attention on the effect of HDAC7 on β–catenin function. Immunoprecipitation assay showed that the interaction between HDAC7 and β-catenin decreased the level of β-catenin acetylation at Lys49 and then downregulated phosphorylation at Ser45. Transfection of ZNF326 mutants that could not bind to the HDAC7 promoter abrogated this effect. In addition, after overexpressing mutant ZNF326 and knocking down HDAC7 or adding TSA to inhibit the function of HDACs, the down-regulation effect of ZNF326 on β-catenin acetylation at Lys49 and phosphorylation at Ser45 disappeared. Thus, it is clear that the change in β-catenin acetylation and phosphorylation is achieved by ZNF326 through HDAC7. Our results support the fact that many non-histone proteins found in recent studies can also serve as substrates of HDACs to participate in many biological processes [28,29,30]. The high expression of HDACs is usually associated with progression of the tumour and poor prognosis of the patient [31,32,33]; the current HDAC inhibitors are also becoming a new class of anti-tumour drugs . Notably, we found the glioma cells with co-transfecion of ZNF326 and shRNA-HDAC7 still exerted the higher tumor formation ability compared to the control group, or in vitro, transfection siRNA-HADC7 just partially abolished the promotion effect of ZNF326 in glioma proliferation, which indicates ZNF326 also could promotes glioma progression via HDAC7-independent manner. In our previous study, eg. ERCC1, a gene closely related to cell proliferation, is also the target gene of ZNF326 , If ZNF326 can promote the proliferation of glioma through both of them, HDAC7 knockdown alone cannot completely inhibit the growth of glioma cells induced by ZNF326. Therefore, which gene, such as HDAC7, ERCC1, LTBP4 and ZNF383, could play a major role in glioma proliferation need further study and investigation.
So, is there any correlation between the effect of HDAC7 on the acetylation and phosphorylation of β-catenin? Why does HDAC7 inhibit the phosphorylation of β-catenin while the total level of β-catenin remains unchanged both in vitro and in vivo? HDAC7 interacts with β-catenin in the cytoplasm and decreases the acetylation level of β-catenin at Lys-49. Thereafter, the steric hindrance was changed and phosphorylation level with the adjacent 45th serine was inhibited. Therefore, up-regulation of HDAC7 expression can decrease the level of p-β-catenin, which is consistent with previous studies on HDAC6 . In order to explain the decrease in p-β-catenin, the total level of β-catenin did not yet change significantly, we examined whether CK1α expression changes after dual-regulation of ZNF326 and HDAC7, because CK1α-mediated phosphorylation of β-catenin at the serine 45 residue is a key step for β-catenin degradation . The results showed that CK1α did not change while ZNF326 and HDAC7 were changed (Additional file 2: Figure S6). HDAC7 inhibits the phosphorylation of β-catenin by inhibiting the acetylation of β-catenin, while non-phosphorylated β-catenin imports the nucleus and activates the Wnt pathway. To conclude, the inhibition of phosphorylation of β-catenin by HDAC7 only changed its cytosolic-nuclear redistribution.
In the initial stage of our experiment, we noticed that the opposite results were reported in three negative breast cancers (TNBC) by Rangel et al, who found that the expression of ZNF326 in TNBC specimens was decreased and ZNF326 overexpression was able to weaken the ability of the transplanted tumour formation , indicating ZNF326 acts as a tumor suppressor gene. In this regard, we analyzed the possible reasons: the previous experiments used rabbit polyclonal antibody against ZNF326, and there may be some differences between the polyclonal antibody and the monoclonal antibody we used. We also attempted to compare the possible differences between them, but the polyclonal antibody was not available. In addition, Madhumitha et al. also reported that PRMT5/WDR77 can regulate the expression of ZNF326 in the TNBC cell line MDA-MB-231 by regulating ZNF326 promotor methylation, which then regulates the expression of the downstream genes such as REPIN1/AP4 and ST3GAL6 . These data suggest that ZNF326 may play different biological roles in tumours of different tissue types and different genetic backgrounds. It is not clear why ZNF326 is overexpressed in gliomas, whether it is due to gene amplification, regulation of upstream genes or impaired metabolism, and whether the high expression of ZNF326 is tissue-specific, etc. need further study and confirmation.
To conclude, ZNF326 expression is deregulated in human glioma, and its up-regulation is associated with tumour grade in patients with glioma. Both in vivo and in vitro experiments demonstrated that ZNF326 plays a role as a tumour-promoting factor through activation of the Wnt pathway. ZNF326, as a transcription activator, binds to the HDAC7 promoter region and activates the transcription of HDAC7. HDAC7 interacts with β-catenin and reduces the level of β-catenin acetylation at Lys49, leading to the reduction of phosphorylation levels at Ser45. This leads to the nuclear accumulation of β-catenin and activates the Wnt signalling pathway. The β-catenin in the nucleus combines with ZNF326 and acts as a transcriptional co-activator of Wnt target genes (Fig. 9). Therefore, ZNF326-HDAC7-β-catenin forms a regulatory loop that activates and strengthens Wnt pathway activity, thereby promoting the malignant phenotype of glioma cells. These results not only reveal the role and mechanism of ZNF326 in carcinogenesis and glioma progression, but also suggest new targets for drug development and drug resistance research.
In summary, our findings demonstrated the expression level of ZNF326 in glioma tissue was positively correlated with its grades. ZNF326 could activate HDAC7 transcription via its transcriptional activation domain and zinc-finger structures. The interaction of the up-regulated HDAC7 with β-catenin led to β-catenin posttranscriptional modification and promoted its import into the nucleus, then activates the Wnt signalling pathway. On the other hand, ZNF326 directly associated with β-catenin in the nucleus, and enhanced the binding of β-catenin to TCF-4, serving as a co-activator in stimulating Wnt pathway. Therefore, ZNF326 promotes the malignant phenotype of human glioma via ZNF326-HDAC7-β-catenin signalling, which is one of its biological mechanisms.
Reverse transcription quantitative polymerase chain reaction
Clarke J, Butowski N, Chang S. Recent advances in therapy for glioblastoma. Arch Neurol. 2010;67:279–83.
Jemal A, Siegel R, Ward E, Murray T, Xu J, Smigal C, et al. Cancer statistics. CA Cancer J Clin. 2006;56:106–30.
Liu HL, Fan CH, Ting CY, Yeh CK. Combining microbubbles and ultrasound for drug delivery to brain tumours: current progress and overview. Theranostics. 2014;4:432–44.
Galli R, Binda E, Orfanelli U, Cipelletti B, Gritti A, De Vitis S, et al. Isolation and characterization of tumorigenic, stem-like neural precursors from human glioblastoma. Cancer Res. 2004;64:7011–21.
Singh SK, Clarke ID, Terasaki M, Bonn VE, Hawkins C, Squire J, et al. Identification of a cancer stem cell in human brain tumours. Cancer Res. 2003;63:5821–8.
Kahlert UD, Maciaczyk D, Doostkam S, Orr BA, Simons B, Bogiel T, et al. Activation of canonical Wnt/beta-catenin signalling enhances in vitro motility of glioblastoma cells by activation of ZEB1 and other activators of epithelial-to-mesenchymal transition. Cancer Lett. 2012;325:42–53.
Cilibrasi C, Riva G, Romano G, Cadamuro M, Bazzoni R, Butta V, et al. Resveratrol impairs glioma stem cells proliferation and motility by modulating the Wnt signalling pathway. PLoS One. 2017;12:e0169854.
Denysenko T, Annovazzi J, Cassoni P, Melcarne A, Mellai M, Schiffer D. WNT/beta-catenin signalling pathway and downstream modulators in low- and high-grade glioma. Cancer Genom Proteom. 2016;13:31–45.
Phillips HS, Kharbanda S, Chen R, Forrest WF, Soriano RH, Wu TD, et al. Molecular subclasses of high grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. Cancer Cell. 2006;9:157–73.
Chen L, Huang K, Han L, Shi Z, Zhang K, Pu P, et al. Beta-catenin/TCF-4 complex transcriptionally regulates AKT1 in glioma. Int J Oncol. 2011;39:883–90.
Gong A, Huang S. FoxM1 and Wnt/beta-catenin signalling in glioma stem cells. Cancer Res. 2012;72:5658–62.
Zheng H, Ying H, Wiedemeyer R, Yan H, Quayle SN, Ivanova EV, et al. PLAGL2 regulates Wnt signalling to impede differentiation in neural stem cells and gliomas. Cancer Cell. 2010;17:497–509.
Jin X, Jeon HY, Joo KM, Kim JK, Jin J, Kim SH, et al. Frizzled 4 regulates stemness and invasiveness of migrating glioma cells established by serial intracranial transplantation. Cancer Res. 2011;71:3066–75.
de Sousa E, Melo F, Vermeulen L. Wnt signalling in cancer stem cell biology. Cancers (Basel). 2016;8(7). https://doi.org/10.3390/cancers8070060.
Clevers H. Wnt/β-catenin signalling in development and disease. Cell. 2006;127:469–80.
Clevers H, Nusse R. Wnt/β-catenin signalling and disease. Cell. 2012;149:1192–205.
Cruciat CM. Casein kinase 1 and Wnt/β-catenin signalling. Cell Biol. 2014;31:46–55.
Aberle H, Bauer A, Stappert J, Kispert A, Kemler R. β-Catenin is a target for the ubiquitin-proteasome pathway. EMBO J. 1997;16:3797–804.
Lee JY, Nakane Y, Koshikawa N, Nakayama K, Hayashi M, Takenaga K. Characterization of a zinc finger protein ZAN75: nuclear localization signal, transcriptional activator activity, and expression during neuronal differentiation of P19 cells. DNA Cell Biol. 2000;19:227–34.
Lee JY, Kambe M, Hayashi M, Takenaga K. Cloning and characterization of a novel zinc finger protein that associates with nuclear matrix. DNA Cell Biol. 1998;17:849–58.
Han Q, Lin X, Zhang X, Jiang G, Zhang Y, Miao Y, et al. WWC3 regulates the Wnt and hippo pathways via Dishevelled proteins and large tumour suppressor 1, to suppress lung cancer invasion and metastasis. J Pathol. 2017;242:435–47.
Wu J, Zhang X, Han Q, Han X, Rong X, Wang M, et al. ZNF326 promotes proliferation of non-small cell lung cancer cells by regulating ERCC1 expression. Lab Invest. 2018. https://doi.org/10.1038/s41374-018-0148-y.
Bradley EW, Carpio LR, Olson EN, Westendorf JJ. Histone deacetylase 7 (Hdac7) suppresses chondrocyte proliferation and β-catenin activity during endochondral ossification. J Biol Chem. 2015;290:118–26.
Margariti A, Zampetaki A, Xiao Q, Zhou B, Karamariti E, Martin D, et al. Histone deacetylase 7 controls endothelial cell growth through modulation of beta-catenin. Circ Res. 2010;106:1202–11.
Li Y, Zhang X, Polakiewicz RD, Yao TP, Comb MJ. HDAC6 is required for epidermal growth factor-induced beta-catenin nuclear localization. J Biol Chem. 2008;283:12686–90.
Mann M, Iaconelli J, Huang JH, Berkovitch SS, Chattopadhyay S, Mazitschek R, Schreiber SL. Axin-mediated CKI phosphorylation of beta-catenin at Ser 45: a molecular switch for the Wnt pathway. Genes Dev. 2002;16:1066–76.
Close P, East P, Dirac-Svejstrup AB, Hartmann H, Heron M, Maslen S, et al. DBIRD complex integrates alternative mRNA splicing with RNA polymerase II transcript elongation. Nature. 2012;484:386–9.
Tao R, Zhang B, Li Y, King JL, Tian R, Xia S, et al. HDAC-mediated deacetylation of KLF5 associates with its proteasomal degradation. Biochem Biophys Res Commun. 2018;500:777–82.
Dai J, Bercury KK, Jin W, Macklin WB. Olig1 acetylation and nuclear export mediate oligodendrocyte development. J Neurosci. 2015;35:15875–93.
Vishwakarma BA, Nguyen N, Makishima H, Hosono N, Gudmundsson KO, Negi V, et al. Runx1 repression by histone deacetylation is critical for Setbp1-induced mouse myeloid leukemia development. Leukemia. 2016;30:200–8.
Lei Y, Liu L, Zhang S, Guo S, Li X, Wang J, et al. Hdac7 promotes lung tumorigenesis by inhibiting Stat3 activation. Mol Cancer. 2017;16:170.
Witt AE, Lee CW, Lee TI, Azzam DJ, Wang B, Caslini C, et al. Identification of a cancer stem cell-specific function for the histone deacetylases, HDAC1 and HDAC7, in breast and ovarian cancer. Oncogene. 2017;36:1707–20.
Peixoto P, Blomme A, Costanza B, Ronca R, Rezzola S, Palacios AP, et al. HDAC7 inhibition resets STAT3 tumorigenic activity in human glioblastoma independently of EGFR and PTEN: new opportunities for selected targeted therapies. Oncogene. 2016;35:4481–94.
Lasko LM, Jakob CG, Edalji RP, Qiu W, Montgomery D, Digiammarino EL, et al. Discovery of a selective catalytic p300/CBP inhibitor that targets lineage-specific tumours. Nature. 2017;550:128–32.
Rangel R, Guzman-Rojas L, Kodama T, Kodama M, Justin Y, Neal G, et al. Identification of new tumour suppressor genes in triple-negative breast cancer. Cancer Res. 2017. https://doi.org/10.1158/0008-5472.CAN-17-0785.
Rengasamy M, Zhang F, Vashisht A, Song WM, Aguilo F, Sun Y, et al. The PRMT5/WDR77 complex regulates alternative splicing through ZNF326 in breast cancer. Nucleic Acids Res. 2017;45:11106–20.
The shRNA-ZNF326 sequence was provided by Dr. Roberto Rangel and Professor Nancy A. Jenkins at the Anderson Cancer Center in the United States, and HA-CBP plasmid was provided by Professor Cao Liu at the Department of Transformation Medicine in China Medical University.
This study was supported by the National Natural Science Foundation of China (No. 81602209 to Minghao Wang; No.81702881 to Xinmiao Yu).
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files. Further details were available from the corresponding author upon request.
Ethics approval and consent to participate
The study was approved by the Medical Research Ethics Committee of China Medical University and inform consent was obtained from all patients. The nude mice used in this study were treated following the experimental animal ethics guidelines issued at China Medical University and approved by the Institutional Animal Research Committee of China Medical University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Materials and Methods (DOCX 22 kb)
Figure S1. The positive correlation between ZNF326 and HDAC7 expressions in gliomas with different grade. (A, D, G, J): Representative images of glioma specimens using H&E staining were shown (Grade I-IV, Magnification 400×). IHC staining of ZNF326 and HDAC7 expressions in gliomas with different grades outlined below: (B, C): ZNF326 and HDAC7 were both negative expressed in pilocytic astrocytoma (grade I, Magnification 400×). (E, F): ZNF326 (<25%, +) and HDAC7 (<25% positive cells) were detected in diffuse astrocytoma (grade II, Magnification 400×). (H, I): ZNF326 was nuclear positive (nearly 50%, ++), while HDAC7 (>50% positive cells, ++) was perinuclear cytoplasm positive in anaplastic astrocytoma (grade III, Magnification 400×), (K, L): ZNF326 and HDAC7 were strongly expressed (nucleus and cytoplasm, respectively, >75% positive cells, +++) in glioblastoma with grade IV (Magnification 400×). Figure S2. ZNF326 expression is associated with expression of Wnt target genes. A: Bioinformatics KEGG test was used to analyse the correlation between ZNF326 and the Wnt pathway. B: Positive correlation between ZNF326 and the four common Wnt signalling pathway target genes in glioma, analysed at the GEPIA website. Figure S3. Positive correlation between ZNF326 and HDAC7 in glioma, analysed at the GEPIA website. Figure S4. Positive correlation between HDAC7 and Wnt signalling pathway target genes in glioma, analysed at the GEPIA website. Figure S5. ZNF326 and siRNA-HDAC7 were co-transfected, or TSA (10nM) was added in U87 cells, and Transwell assays were performed to detect the changes in the invasiveness of the glioma cells. Figure S6. (A-D): ZNF326, siRNA-ZNF326, HDAC7 and siRNA-HDAC7 were transfected in U87 cells, respectively, and immunoblotting assay was performed to detect the changes in the expression of β-catenin and CK1α. GAPDH was used as a loading control. (ZIP 12193 kb)
Table S1. The correlation between the expression of ZNF326 and HDAC7 in glioma. (DOCX 14 kb)