- 1Department of Neurosurgery, The Fourth Affliated Hospital of Guangxi Medical University, Liuzhou, China
- 2College of Traditional Chinese Medicine, North China University of Science and Technology, Tangshan, China
Gliomas are the most common and aggressive malignancies of the central nervous system. Histone deacetylases (HDACs) are important targets in cancer treatment. They regulate complex cellular mechanisms that influence tumor biology and immunogenicity. However, little is known about the function of HDACs in glioma. The Oncomine, Human Protein Atlas, Gene Expression Profiling Interactive Analysis, Broad Institute Cancer Cell Line Encyclopedia, Chinese Glioma Genome Atlas, OmicShare, cBioPortal, GeneMANIA, STRING, and TIMER databases were utilized to analyze the differential expression, prognostic value, and genetic alteration of HDAC and immune cell infiltration in patients with glioma. HDAC1/2 were considerable upregulated whereas HDAC11 was significantly downregulated in cancer tissues. HDAC1/2/3/4/5/7/8/11 were significantly correlated with the clinical glioma stage. HDAC1/2/3/10 were strongly upregulated in 11 glioma cell lines. High HDCA1/3/7 and low HDAC4/5/11 mRNA levels were significantly associated with overall survival and disease-free survival in glioma. HDAC1/2/3/4/5/7/9/10/11 are potential useful biomarkers for predicting the survival of patients with glioma. The functions of HDACs and 50 neighboring genes were primarily related to transcriptional dysregulation in cancers and the Notch, cGMP-PKG, and thyroid hormone signaling pathways. HDAC expression was significantly correlated with the infiltration of B cells, CD4+ T cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells in glioma. Our study indicated that HDACs are putative precision therapy targets and prognostic biomarkers of survival in glioma patients.
Introduction
Glioma is the most common tumor of the central nervous system (CNS). It accounts for 30% of all CNS tumors and 80% of all malignant brain tumors (Ostrom et al., 2019). Progress has been made in glioblastoma (GBM) treatment by combining maximal surgical resection with radiotherapy and concurrent and adjuvant temozolomide chemotherapy. Nevertheless, the two- and five-year survival rates are only 25 and 10%, respectively, (Stupp et al., 2005; Stupp et al., 2009). Tumor heterogeneity is a major challenge in precise glioma diagnosis and therapy (Neftel et al., 2019). Hence, molecular signatures for glioma are urgently required (Song et al., 2014).
Growing evidence suggests that epigenetic and genetic changes are pivotal in malignant disease onset and progression (Jones and Baylin, 2007). Histone deacetylases (HDACs) are important regulators of gene expression and control a broad array of physiological processes such as differentiation, apoptosis, survival, proliferation, and autophagy. They are also involved in cancer pathogenesis (Glozak and Seto, 2007; Wilson et al., 2008; Kang et al., 2014). There are four HDAC families and 18 individual HDACs. The zinc-dependent HDAC family comprises 11 isoforms divided into class I (HDAC1–3 and 8), class IIa (HDAC4, 5, 7, and 9), class IIb (HDAC6 and 10), and class IV (HDAC11) (Kashyap and Kakkar, 2020). Class III HDACs are homologous to yeast Sir2 and participate in transcriptional silencing. However, they have a deoxyhypusine synthase-like NAD/FAD-binding domain that is distinct from those of the other HDAC classes (Park et al., 2021). HDAC1 and two act mainly via nucleosome remodeling and deacetylase, switch independent 3, mitotic deacetylase, and corepressor of REST complexes. HDAC3 is recruited only by the nuclear receptor corepressor complex (Millard et al., 2017). HDAC6 regulates Hsp90, tau, and the cytoskeleton by interacting with tubulin and cortactin. It recognizes ubiquitinated proteins and induces aggregate formation (Hubbert et al., 2002; Kovacs et al., 2005; Dompierre et al., 2007; Zhang et al., 2007; Cohen et al., 2011; Simões-Pires et al., 2013). HDAC10 is a polyamine deacetylase (Ho et al., 2020), and HDAC11 has limited homology to class I and II enzymes (Li et al., 2015).
HDACs play important roles in many diseases. HDAC1 promotes medulloblastoma growth and affects cell cycle progression, microtubule dynamics, and DNA damage response (Abdelfattah et al., 2018). Preclinical studies suggest a role for HDAC1 in the epigenetic regulation of sarcoma tumorigenesis (Bennani-Baiti and Idriss, 2011). HDAC2 is highly expressed in sarcomas (Pacheco and Nielsen, 2012). HDAC3 is a target of the HDAC inhibitors approved by the US Food and Drug Administration to treat lymphoid malignancies (Stengel et al., 2019). HDAC4 overexpression has been observed in acute promyelocytic leukemia (Chauchereau et al., 2004), B cell lymphoma (Sandhu et al., 2012), and multiple myeloma (Kikuchi et al., 2015). HDAC6 is protective in tauopathy and suppresses aberrant tau accumulation. Chronic HDAC6 loss results in accelerated tau pathology, cognitive dysfunction, and reduced survival (Trzeciakiewicz et al., 2020). HDAC7 is a bona fide transcriptional repressor vital to B cell development (Azagra et al., 2016). HDAC8 modulates p53 activity and ensures long-term hematopoietic stem cell maintenance and survival under stress (Hua et al., 2017). HDAC5, HDAC7, and HDAC9 play protective roles during cardiac hypertrophy (Mckinsey, 2011). HDAC11 performs a key function in oncogene-induced, non-homeostatic hematopoiesis (11 deficiency disrupt, 2020).
Nevertheless, the mechanisms by which HDACs are inhibited or promoted in glioma development and progression are unclear. No comprehensive bioinformatics analysis has been conducted to identify the putative role of HDACs in glioma. In the present study, we mined numerous large databases to analyze HDAC expression, mutation, and function and immune infiltration and determine their potential oncogenic and prognostic value in glioma.
Materials and Methods
Oncomine Database
Oncomine (https://www.oncomine.org) contains 715 gene expression datasets and 86,733 samples. It is the largest oncogene chip database and incorporated data mining platform (§ et al., 2007). Transcriptional mRNA expression data for 11 HDACs in various cancers and their corresponding normal adjacent tissues were retrieved from Oncomine. The search content and threshold values were as follows: keywords, HDAC1–HDAC11; primary filter, cancer vs. normal; cancer type, absolute value of log2 │Fold Change│ > 1.5 and p < 0.05; and gene rank, 10%. p-values were calculated using the Student’s t-test.
Human Protein Atlas
The Human Protein Atlas (https://www.proteinatlas.org) is an online tool including immunohistochemistry (IHC) expression data for protein distribution and expression in 20 cancer tissues, 48 human healthy tissues, 47 cell lines, and 12 blood cell types (Asplund et al., 2012). IHC images were used to compare HDAC protein expression levels among normal and cancer tissues. Localization of HDAC immunofluorescence expression in glioma cell lines was explored.
Gene Expression Profiling Interactive Analysis
GEPIA (http://gepia.cancer-pku.cn/detail.php) is a new online interactive web server enabling users to examine RNA sequencing expression data for tumors, normal tissues, and samples in the Genotype-Tissue Expression projects and The Cancer Genome Atlas (TCGA). GEPIA is based on a criterion processing pipeline. It provides customizable functions such as tumor/normal differential expression analysis, profiling by cancer type or pathological stage, similar gene detection, and patient survival, correlation, and dimensionality reduction analyses (Tang et al., 2017). In the present study, differential gene expression analysis was used to compare tumors and normal tissues using GEPIA. The Student’s t-test was used to generate p-values for the expression analysis. A Kaplan–Meier curve for the patient survival analysis was also plotted.
Broad Institute Cancer Cell Line Encyclopedia Database
The CCLE (https://www.broadinstitute.org/ccle) project conducts detailed genetic and pharmacologic characterizations of large panels of human cancer models. It develops integrated computational analyses linking distinct pharmacologic vulnerability to genomic patterns and provides public access to genomic data for the analysis and visualization of approximately 1,000 cell lines (Ghandi et al., 2019). The mRNA expression levels of 11 HDACs in multiple glioma cell lines were compared. The RNA expression dataset for 12 glioma cell lines was retrieved from CCLE and the HDAC expression levels were plotted in heat maps.
Chinese Glioma Genome Atlas Database
CGGA (http://www.cgga.org.cn) is the largest glioma genome database in China. It provides multiple omics and matches clinical data for >2,000 primary and recurrent samples excised from Chinese cohorts (Zhao et al., 2020). In the present study, online tools and the mRNAseq_693 dataset were used to analyze the pathological stages associated with HDAC expression in patients with glioma.
GeneMANIA Database and OmicShare
GeneMANIA (http://www.genemania.org) is a well maintained, user-friendly gene list analysis web interface used to derive hypotheses based on gene functions (Max et al., 2018). Here, GeneMANIA was utilized to construct a gene–gene interaction network of HDAC family members based on their physical interactions, predictions, co-expression levels, co-localizations, genetic interactions, drug-interactions-2020, transcriptional-factor-targets-2020, and so on. GeneMANIA was also used to evaluate HDAC functions. OmicShare (http://www.omicshare.com/tools) is a comprehensive platform for data processing and analysis, learning biological information, and sharing scientific research knowledge. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the free online data analysis platform OmicShare Tools (https://www.omicshare.com/tools).
cBioPortal Databases
The cBioPortal (http://cbioportal.org) is a free asset that downloads large-scale cancer genomics datasets incorporating 245 cancer studies. The cBioPortal was used to explore genetic alterations in HDACs in glioma (Cerami et al., 2012). Here, the cBioPortal database was applied to explore 804 mutations in the brain lower-grade glioma (LGG; TCGA, and PanCancer Atlas) and glioblastoma multiforme (TCGA and Firehose Legacy) datasets for nerve center tumors. The distribution of HDAC mutations in glioma was calculated.
STRING
STRING (https://string-db.org/) is a protein interaction website providing a comprehensive, objective global network. It presents its data using a unique set of computational predictions (Szklarczyk et al., 2019). A protein–protein interaction (PPI) network analysis was conducted through STRING to collect and integrate the expression of HDACs and their potential interactions.
Tumor Immune Estimation Resource Analysis
The TIMER (https://cistrome.shinyapps.io/timer/) (Li et al., 2017) was used to evaluate the infiltration of CD8+ T cells, B cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells in tumors and correlate them with HDAC expression in GBM and LGG. The survival module was used to establish prognostic relevance. p < 0.01 was considered statistically significant.
Quantitative Real-Time Polymerase Chain Reaction of Tissues
Total RNA was extracted from tissue specimens using Animal RNA Isolation Kit (Invitrogen, Beyotime, Shanghai, China) according to the manufacturer’s instructions, and RNA was reversely transcribed into cDNA using Transcription First Strand cDNA synthesis kit (Beyotime, Shanghai, China). Quantitative real-time PCR (qRT-PCR) analyses were quantified with BeyoFast™ SYBR Green (Beyotime, Shanghai, China). The relative expression of HDACs were calculated based on the 2-ΔΔCt method with GADPH as an internal reference. qRT-PCR primers used in the present study were as follows:
HDAC1 forward primer, 5′- TCAAGATGGCCTGAGCAAGG-3’; HDAC1 reverse primer, 5′-TGTGCGCTGGTCCCTATCTA-3′; HDAC2 forward primer, 5′- TTCCAAGCCCGACTGTGAGA-3′; HDAC2 reverse primer, 5′-ACCTGTTAGAGCCAGTAAGCAC-3′; HDAC3 forward primer, 5′- GGCCGATGCTGAAGAGAGAG-3′; HDAC3 reverse primer, 5′-GGGGATACCCAGTTCAGACC-3′; HDAC4 forward primer, 5′- CCTGTGGCCACTGCTCTAAA-3′; HDAC4 reverse primer, 5′-AATGCCATTCTCGGTGCTGA-3′; HDAC5 forward primer, 5′- CAGGCTGCTGCCACTCAAGA-3′; HDAC5 reverse primer, 5′-CACAATGATGAAGCCCAGAGGG-3′; HDAC6 forward primer, 5′- AGCTAGTCCTGTGCCGCTA-3′; HDAC6 reverse primer, 5′-TGTAGGTAATGCCGCTGTGG-3′; HDAC7 forward primer, 5′- ACCTGCGAGTGGGCCAAAG-3′; HDAC7 reverse primer, 5′-TACGGCACTTCGCTTGCTC-3′; HDAC8 forward primer, 5′- AATGAGCCCCATCGAATCCAG-3′; HDAC8 reverse primer, 5′-GATATCCTCCCTCTTTCCCCCTA-3′; HDAC9 forward primer, 5′- ACTTGGTTACCCCAAGGAGC-3′; HDAC9 reverse primer, 5′-ATGCAGTGGAGGTCAGATGC-3′; HDAC10 forward primer, 5′- TCACTGGACAAGCCTCCAC-3′; HDAC10 reverse primer, 5′-GGCAAGATCGTCGTCCTGAA-3′; HDAC11 forward primer, 5′- TGCTAAAGAGGCCATCAGGC′; HDAC10 reverse primer, 5′-TGAGGATGGAGTCGGCGATA-3′; GAPDH forward primer, 5′- CCGCATCTTCTTGTGCAGTG′; GAPDH reverse primer, 5′- TCCCGTTGATGACCAGCTTC -3′
Results
Comparison of HDAC Expression Levels in Glioma and Normal Tissue Samples
To explore the potential prognostic and therapeutic values of various HDACs, we used the Oncomine and HPA databases to analyze their mRNA and protein expression levels in patients with glioma. TCGA statistics revealed that HDAC1 expression was 3.131-fold higher (p = 3.08E-8) in ductal brain glioblastoma than in normal tissues (Table 1). Pomeroy (Pomeroy et al., 2002) has reported a 3.133-fold increase in HDAC2 in desmoplastic medulloblastoma than in normal tissue samples (p = 5.47E-5). Gutmann has (Gutmann et al., 2002) observed a 2.483-fold increase in HDAC3 in pilocytic astrocytoma than in normal samples (p = 0.021). Sun Brain Statistics (Sun et al., 2006) disclosed HDAC6 overexpression in glioblastoma tissues compared with normal tissues and the fold change was 3.221 (p = 6.66E-13). The database search returned a 2.608-fold increase in HDAC6 mRNA expression in anaplastic astrocytoma compared with normal tissues (p = 2.34E-6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-t001.jpg)
TABLE 1. Significant changes in histone deacetylase transcription between various types of glioma and normal tissues (Oncomine database).
We also explored protein-level HDAC expression using the Human Protein Atlas database. Low and medium HDAC1/4/8/9/10 protein expression levels were determined for normal tissues, and high HDAC1/4/8/9/10 protein expression levels were found in glioma tissues (Figures 1A,D,G,H). Figures 2C,E,F show that HDAC3/5/6 proteins are not expressed in normal gastric tissues but are detected at low and medium levels in glioma tissues.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g001.jpg)
FIGURE 1. Comprehensive analysis of prognosis and histone deacetylase expression in glioma. HDAC1/2/3/4/5/7/8/11 expression levels were correlated with pathological stage of glioma patients (p < 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g002.jpg)
FIGURE 2. Representative immunohistochemistry images of histone deacetylase family members in glioma and normal tissues (Human Protein Atlas Database). All tissues were observed at ×100 magnification.
Relationships Between HDAC mRNA Levels and Clinicopathological Parameters of Patients with Glioma
To explore distinct HDAC expression in patients with glioma, HDAC mRNA expression was analyzed using the Oncomine database (Figure 3). The Oncomine data revealed that compared with normal tissues, the HDAC1/2/3/6 transcriptional levels were significantly elevated and those of HDAC2/4/5/11 were significantly decreased in CNS cancer tissues. The GEPIA dataset was used to compare HDAC mRNA expression between glioma and normal tissues. HDAC1 and HDAC2 were relatively upregulated and HDAC11 was downregulated in LGG and GBM tumor tissues (Figure 4). We also compared HDAC expression levels across glioma tumor stages. HDAC1/2/3/4/5/7/8/11 expression significantly varied across glioma stages, whereas HDAC6/9/10 expression did not (Figure 1). Therefore, qRT-PCR was performed to verify the expression of HDAC family gene mRNA in glioma tumor tissues. HDAC1/2/3/4/5/8/11 expression levels were correlated with pathological stage of glioma patients (p < 0.05) (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g003.jpg)
FIGURE 3. Histone deacetylase transcription levels in various cancers (Oncomine). Numbers in colored cells show quantities of datasets with statistically significant target gene mRNA upregulation (red) or downregulation (blue). The following criteria were used: p-value, 0.01; fold change, two; gene rank, 10%; data type, mRNA; analysis type, cancer vs. normal tissue. CNS: central nervous system.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g004.jpg)
FIGURE 4. Gene expression profiling interactive analysis data for histone deacetylase (HDAC) mRNA expression in glioma including glioblastoma (GBM) and lower-grade glioma (LGG) compared with NB tissue (A,B) Scatter diagram and Box plot showing that HDAC1/2 expression levels in GBM (n = 163) and LGG (n = 518) tumor tissues were higher than those in normal tissues (n = 207). HDAC11 expression levels were lower in GBM tumor tissues than in normal tissues (p < 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g005.jpg)
FIGURE 5. Comprehensive analysis of prognosis and histone deacetylase expression in glioma. HDAC HDAC1/2/3/4/5/8/11 expression levels were correlated with pathological stage of glioma patients (p < 0.05).
HDACs mRNA Expression and Immunofluorescence in Glioma Cells
We compared HDAC expression levels among the foregoing databases and used CCLE to explore HDAC1/2 expression levels in glioma cell lines. All HDACs were strongly upregulated in glioma cell lines (Figure 6C). We also found that HDAC1/2/3/10 were highly expressed in 11 glioma cell lines from the CCLE database, namely, SW1783, TM31, KNS60, M059K, GOS3, SF126, KALS1, LN229, GI1, KNS81, KG1C, and AM38 (Figure 6B). We searched the cell atlas of the HDAC family in glioma and observed immunofluorescence of HDAC1/2/3/4/5/7/9 (Figure 6A). In the human protein map, most HDAC proteins were localized to and upregulated in human glioma cell nuclei.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g006.jpg)
FIGURE 6. Cell line histone deacetylase (HDCA) mRNA expression levels (Cancer Cell Line Encyclopedia) (A) HDAC1/2/3/4/5/7/9 immunofluorescence was expressed in glioma U-251 MG cell nuclei (B) HDAC1/2/3/10 mRNA was highly expressed in 12 glioma cell lines (C) Box plot was sorted and colored according to average distribution of HDAC1 and HDAC2 expression in various tumor cell lines. Immunofluorescence analysis of HDAC1/2/3/4/5/7/9 in U-251 MG.
Prognostic Significance of HDACs in Patients with Glioma
We evaluated the prognostic significance of HDACs in patients with glioma using the GEPIA and CGGA databases. A log-rank test analysis revealed that high HDCA1/3/7 and low HDAC4/5/11 mRNA levels were significantly associated with overall survival (OS) and disease-free survival (DFS) (p < 0.05) (Figure 7). HDAC9 mRNA downregulation and HDAC10 mRNA upregulation had prognostic value in DFS, whereas HDAC2 mRNA upregulation had prognostic value in OS (p < 0.05). The GEPIA database showed that HDAC1/2/3/4/5/7/9/10/11 mRNA expression levels were significantly associated with glioma patient prognosis and could, therefore, be exploited as biomarkers for the prediction of glioma patient survival.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g007.jpg)
FIGURE 7. Prognostic value of histone deacetylase mRNA expression in patients with glioma (Gene expression profiling interactive analysis). *p < 0.05. HR, hazard ratio; OS, overall survival; DFS, disease-free survival.
Genetic HDAC Alteration Analysis in Glioma Patients
Genetic alterations in the HDACs of glioma patients were examined using TCGA database and the cBioPortal online tool. HDACs were altered in 804 samples from 1,118 glioma patients in TCGA databases of brain LGG and GBM. The alteration rates were 12.05% (62/514) and 5.52% (16/290), respectively. Deep deletion accounted for most of the observed changes (Figure 8A). Figure 8B shows that genetic HDAC alterations occurred in 78 (10%) of queried samples and that the alteration rates in the individual sequence were in the range of 0.1–3%. HDAC4, HDAC6, and HDAC10 were ranked as the top three of the eleven HDAC members and their deep deletion rates were 3, 1.4, and 1.4%, respectively, (Figure 8B). The “Survival” tab, the Kaplan–Meier plot, and the log-rank test plotted survival curves indicated that cases with or without HDAC alterations had no correlation with OS or DFS (Figures 8C,D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g008.jpg)
FIGURE 8. Genetic histone deacetylase (HDAC) alteration analysis in patients with OS (cBioPortal) (A) HDAC alterations (B) OncoPrint tab summary of query of HDAC alterations. Kaplan–Meier plots comparing (C) overall survival (OS) and (D) disease-free survival (DFS) in patients with and without HDAC gene alterations.
Predicted HDAC and HDAC-Related Neighboring Gene Functions and Pathways in Patients with Glioma
We analyzed 50 neighboring genes significantly related to HDACs and used GeneMANIA to construct an integrated network. Figure 8A shows that according to the weight analysis, HDACs were associated with ubiquitinase-related genes, including USP3/5/13/16/20/22/33/39/44/49/45/51. Gene-related drugs acting on HDACs included panobinostat, vorinostat, romidepsin, SB939, valproic acid, and MGCD-0103. The transcription factors HMEF2_Q6, RSRFC4_Q2, MEF2_02, and MEF2_03 were associated with HDAC-related genes. The LGG and GBM mRNA data in TCGA and the Spearman’s correlation analyses were used to calculate correlations among HDAC family members. The STRING database was used to build a PPI network for the HDAC family (Figures 9A,B). There were significant positive correlations between the following HDACs: HDAC1 with HDAC3; HDAC2 with HDAC3 and HDAC8; HDAC4 with HDAC5 and HDAC6; HDAC5 with HDAC6; HDAC6 with HDAC10; HDAC7 with HDAC1, HDAC3, HDAC6, and HDAC10; HDAC8 with HDAC3 and HDAC6; HDAC9 with HDAC2; and HDAC11 with HDAC4 and HDAC5. There were also significant negative correlations between the following HDACs: HDAC4 with HDAC1 and HDAC3; and HDAC5 with HDAC1 and HDAC3 (Figure 9C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g009.jpg)
FIGURE 9. Protein–protein interactions (PPI) among histone deacetylases (HDACs) (A) Network of HDACs and their 50 neighboring genes was constructed using GeneMANIA (B) PPI network of various HDACs (C) Pearson correlation coefficients among HDACs.
GO enrichment and KEGG pathway analyses of the HDACs and 50 neighboring genes were performed using OmicShare. The top ten most highly enriched GO items were identified (Figures 10A–C). The GO term analysis showed that differentially expressed genes correlated with the HDACs were localized mainly with the histone deacetylase complex, nuclear part, nucleus, nucleoplasm, nuclear lumen, and intracellular organelle lumen; HDACs negatively regulate RNA biosynthesis and transcription, DNA templates, RNA metabolism, histone H3 deacetylation, and histone modification at these sites. HDACs have histone deacetylase-binding, H3-K14-specific NAD-dependent histone deacetylase and histone deacetylase activity, protein deacetylase activity, NAD-dependent histone deacetylase activity, and so on. The KEGG pathway analysis disclosed viral carcinogenesis, apelin signaling, arginine and proline metabolism, transcriptional misregulation in cancer, Notch signaling, and cGMP-PKG signaling. The thyroid hormone signaling pathway and arginine biosynthesis were significantly associated with glioma tumorigenesis and progression (Figure 10D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g010.jpg)
FIGURE 10. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of HDACs and their interactions (OmicShare). GO enrichment analysis of target genes based on (A) cellular component, (B) biological process, and (C) molecular function (D) KEGG pathway enrichment analysis of target genes.
HDAC Expression and Immune Infiltration in Glioma
The TIMER database was used to investigate associations among the HDACs family members and immune cell infiltration. The immune cell count is correlated with cancer cell proliferation and progression (Figure 11). HDAC1 expression in LGG was positively correlated with the infiltration of B cells (p < 0.001), CD4+ T cells (p < 0.01), macrophages (p < 0.01), neutrophils (p < 0.001), and dendritic cells (p < 0.01). HDAC2/8 expression levels in LGG and GBM were positively correlated with purity infiltration (p < 0.01).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840759/fcell-10-840759-HTML-r1/image_m/fcell-10-840759-g011.jpg)
FIGURE 11. Correlations among differentially expressed HDACs and immune cell infiltration (TIMER). Correlations between immune cell abundance and HDAC1–11 expression.
HDAC3 expression in LGG was positively correlated with the infiltration of B cells (p < 0.01), neutrophils (p < 0.01), and dendritic cells (p < 0.01). HDAC7 expression in LGG was positively correlated with the infiltration of CD4+ T cells (p < 0.01), macrophages (p < 0.01), and dendritic cells (p < 0.001). HDAC9 expression in LGG was positively correlated with CD8+ T cell infiltration (p < 0.001). HDAC11 expression in LGG was positively correlated with purity infiltration (p < 0.01) and negatively correlated with neutrophil infiltration (p < 0.01). A Cox proportional hazard model of HDACs and clinical factors in glioma was also evaluated. HDAC3/4 expression, patient age, and CD4+ T cell density were significantly associated with clinical outcomes in GBM (Table 2). HDAC1 expression and patient age were significantly associated with the clinical prognosis of LGG (Table 3).
Discussion
Epigenetic modifications are heritable, reversible genetic changes that do not involve DNA mutation. HDACs are epigenetic regulators that remove the acetyl moieties from acetylated lysine residues of post-translationally modified histone proteins that form the core of the chromatin network. This mechanism compresses chromatin, compromises DNA accessibility, and represses transcription and gene silencing (Sengupta and Seto, 2004). However, distinct roles of each HDAC family member in glioma remain to be clarified. In the present study, we analyzed HDAC expression, mutation, prognostic value, functional enrichment, and immune cell infiltration in patients with glioma. To this end, we consulted various public databases including Oncomine, Human Protein Atlas, GEPIA, cBioPortal, CCLE, GeneMANIA, String, and CGGA. GO enrichment and KEGG pathway were analyzed with OmicShare Tools.
Previous studies have evaluated the roles of certain HDAC family members in glioma (Lodrini et al., 2013; Wang et al., 2017; Lo Cascio et al., 2019; Song et al., 2020). A previous study has indicated that the NFAT2-HDAC1 pathway may maintain the malignant phenotype and promote mesenchymal transition in glioma stem-like cells (GSCs). Hence, GSCs are potential molecular targets for GBM therapy (Song et al., 2020). HDAC1 ablation significantly decreased stemness and proliferation in patient-derived GSCs in a p53-dependent manner. By contrast, it had only minimal impact on normal human neural stem cells and astrocytes (Lo Cascio et al., 2019). HDAC1 knockdown inhibited the epithelial–mesenchymal transition transcription factors TWIST1 and SNAIL, downregulated the mesenchymal marker matrix metalloprotein-9, and upregulated the epithelial marker E-cadherin 34. Hence, HDAC1 may contribute to epithelial–mesenchymal transition in glioma cells (Wang et al., 2017). Our results showed that HDAC1 was highly expressed in glioma tumor tissue and significantly correlated with patient survival and immune infiltrate density. Moreover, HDAC1 had the highest mRNA expression level of all HDACs in 12 glioma cell lines. A Cox proportional hazard model of HDACs and clinical factors disclosed that HDAC1 expression was significantly correlated with clinical outcome in LGG. HDAC1 upregulation was significantly correlated with poor OS and DFS in patients with glioma patients. The foregoing associations merit further investigation. We plan to knock out or use the expression of protease inhibitor HDAC1 in glioma cell lines to observe tumor cell invasion and migration in functional experiments. Besides use tumor-bearing nude mice to further verify whether its HDAC1 can inhibit tumor growth.
Recent observations demonstrated that class I HDACs (HDAC1, HDAC2, and HDAC3) modulate DNA damage signaling (Miller et al., 2010; Thurn et al., 2013), maintain genomic stability, and prevent tumorigenesis in vivo (Dovey et al., 2013; Heideman et al., 2013; Santoro et al., 2013). The mammalian class I deacetylase HDAC2 has been extensively studied. HDAC2 downregulation markedly inhibits tumor growth. Thus, HDAC2 may be an oncogene in tumorigenesis (Zhu et al., 2004). HDAC2 protein overexpression was detected in human gastric, prostate, and breast cancers (Nakagawa et al., 2007; Jinwon et al., 2014). HDAC2 represses gene expression by deacetylating H4K16ac (Ma et al., 2013), determines transcription repression, and participates in the nucleosome remodeling deacetylase complex. In vitro and in vivo experiments demonstrated that silencing TRPS1 inhibited tumor growth, whereas HDAC2 overexpression promoted it (Wang et al., 2018). The results of the present study corroborated these findings. HDAC2 expression was higher in glioma tumor than in normal tissues. HDAC2 expression was also significantly correlated with glioma grade and positively correlated with purity infiltration. However, we found that unlike mRNA expression, immune combination expression was elevated in both glioma and normal tissues possibly because of gene regulation during transcription and translation.
During thymocyte development, HDAC3 is required for positive selection and CD4 lineage development (Stengel et al., 2015). HDAC3 represses CD8 lineage genes to keep double-positive thymocytes in a bipotential state (Philips et al., 2019). Here, HDAC3 expression in LGG was positively correlated with the infiltration of B cells, neutrophils, and dendritic cells. HDAC3 mRNA expression was significantly associated with clinical outcome in GBM. Hence, immunity-related cells may be implicated in glioma pathogenesis.
HDAC4 is a key member of class IIa HDACs and performs a wide variety of functions (Seto and Yoshida, 2014). In various cells, HDAC4 is post-transcriptionally regulated by several microRNAs such as miR-1, miR-29, miR-140, miR-155, miR-200a, miR-206, and miR-365 (Wang et al., 2014). Some studies have shown that the molecular mechanism of glioma formation involved HDAC4-mediated SP1 and KLF5 deacetylation in selective MKK7 transcription and oncogenic JNK/c-Jun cascade activation. These findings were consistent with the findings of our research. The results of the immune combination showed that HDAC4 is highly expressed in glioma tumor tissues. KEGG enrichment analysis has revealed that in cancers, HDAC and its neighboring 50 genes are enriched in microRNAs. HDAC4 upregulation was significantly correlated with poor OS and DFS in all patients with glioma. A Cox proportional hazard model of HDAC and its clinical factors disclosed that HDAC4 expression was significantly correlated with clinical GBM outcome.
HDAC5 is a class II HDAC. Aberrant HDAC5 expression has been observed in multiple cancer types. HDAC5 participates in cell proliferation and invasion, immune response, and maintenance of stemness (Yang et al., 2021). In the present study, HDAC5 expression was low in tumor cell lines and in glioma and normal tissues. However, HDAC5 upregulation was significantly correlated with tumor stage, poor OS, and DFS in all patients with glioma.
HDAC6 and 10 are class IIb HDACs based on their subcellular localizations and expression patterns (Haberland et al., 2009). Previous research has shown that HDAC6-selective inhibitors delayed tumor initiation and progression in vivo without causing any significant adverse effects (Auzmendi-Iriarte et al., 2019). HDAC6 inhibition induces an in vivo delay in tumor growth and downregulates PD‐L1 expression. Several key immunological checkpoint modulators are regulated by HDAC6 (Lienlaf et al., 2016). HDAC10 expression is positively associated with PD-L1 expression and may predict the outcome of patients with non-small cell lung carcinoma. HDAC10 inhibition combined with doxorubicin administration kills neuroblastoma but not non-malignant cells by impeding drug efflux and enhancing DNA damage. This novel mechanism targets chemotherapy resistance (Johannes et al., 2018). In this study, HDAC6 expression was higher in glioma tumor than in normal tissues. Low HDAC6 mRNA expression was significantly correlated with OS and DFS in patients with glioma. HDAC10 expression is moderate in LGG and elevated in high-grade gliomas. However, low HDAC10 mRNA expression was not significantly correlated with DFS in patients with glioma.
HDAC7 sustains human mammary epithelial cell proliferation and favors the establishment of stem‐like cell populations by maintaining a proficient microenvironment. High HDAC8 expression levels in human GBM tissues and GBM-R cell lines were correlated with O-methylguanine-DNA methyltransferase levels (Cutano et al., 2019). An inhibitor of HDAC8 combined with temozolomide administration induced WT-p53-mediated apoptosis via WT-p53-mediated O-methylguanine-DNA methyltransferase inhibition in GBM-R cell lines (Tsai et al., 2021). Elevated HDAC9 expression is associated with poor prognosis and promotes malignancy (Huang et al., 2018; Linares et al., 2019; Li et al., 2020). HDAC9 deficiency promoted tumor progression by decreasing CD8+ dendritic cell (DC) infiltration in the tumor microenvironment. HDAC9 expression was significantly positively correlated with CD8+ cell counts in human lung cancer stroma samples (Ning et al., 2020). In the present study, immune cell infiltration was associated with HDACs in glioma. HDAC9 expression in LGG was positively correlated with CD8+ cell infiltration. HDAC7 expression in LGG was positively correlated with the infiltration of CD4+ T cells, macrophages, and dendritic cells. HDAC8 expression in LGG and GBM was positively correlated with purity infiltration. High HDAC7 mRNA expression was significantly correlated with OS and DFS in patients with glioma. However, HDAC8 mRNA expression levels were not meaningful for prognostic survival analyses of patients with glioma. HDAC9 mRNA expression levels were not associated with OS in patients with glioma.
HDAC11 is the only member of the class IV HDACs (Seto and Yoshida, 2014). There is limited information regarding its functions. Using a novel class of highly selective HDAC inhibitors and genetically deficient mouse models, the foregoing study revealed that HDAC11 was necessary for oncogenic JAK2-driven myeloproliferative neoplasm cell and tissue proliferation and survival (Yue et al., 2020). Other data demonstrated a significant role of HDAC11 in mitotic cell cycle progression and survival of MYCN-amplified neuroblastoma cells and suggested that HDAC11 is potentially a valuable drug target (Thole et al., 2017). The present study arrived at the same conclusions. KEGG enrichment analysis disclosed that enrichment of HDAC and its neighboring 50 genes was associated with the cell cycle. HDAC11 expression was lower in glioma tumor than in normal tissues. HDAC11 downregulation was significantly correlated with tumor stage, poor OS, and DFS in patients with glioma.
There are several limitations to our study. Firstly, we only analyzed data retrieved from online databases, we must validate our findings using large cohorts. Secondly, our HDAC IHC and immunofluorescence data were incomplete. Third, we did not use our cell lines or tissues to further verify the results. Although HDAC family expression data are similar for glioma cell lines and pancancer analyses. Perhaps, Basic experiments must be conducted to elucidate the mechanisms of HDAC family members and determine their impact on epigenetics. At present, a large number of HDAC inhibitors have been used in other tumors, and we can use HDAC inhibitors to act on cells or animals to further observe their inhibition of tumors. The conventional histone acetylases show lack of blood-brain barrier permeability in the treatment of brain diseases, but through the development of new small molecule probes and carriers, they can penetrate the blood-brain barrier well. we may use nano-encapsulation to further break through the blood-brain barrier (Dang, 2014). In the actual human clinical safety and effectiveness test, we can observe the immunohistochemical expression of HDAC in clinical work and combine with traditional glioma gene diagnosis to analyze and judge the prognosis of glioma diseases.
Conclusion
Most of these histone Acetylases are usually adjunctive to other chemoradiotherapy drugs such as temozolomide/bevacizumab. However, most of the current clinical trials are phase I/II clinical trials, and further phase III trials are needed to demonstrate the role of histone acetylases in gliomas. In a phase two Study of Concurrent Radiation Therapy, Temozolomide, and the Histone Deacetylase Inhibitor Valproic Acid for Patients With Glioblastoma, histone inhibitors may improve outcomes for patients with glioma compared with historical data (Krave, 2015). In the present study, we systematically analyzed HDAC expression and prognostic value in glioma and revealed the heterogeneity and complexity of the molecular mechanisms associated with these cancers. We also investigated the mRNA expression patterns, prognostic values, genetic alterations, GO enrichment, and PPI networks of HDACs in patients with glioma. According to the immune prognosis analysis, the age and HDAC1 mRNA expression of patients with LGG were related to disease prognosis. Patient age, CD4+ T cell density, and HDAC3 and HDAC4 mRNA expression were associated with high-grade glioma prognosis. HDAC1 and HDAC2 were relatively upregulated whereas HDAC11 was downregulated in LGG and GBM tumor tissues. High HDCA1 and low HDAC11 mRNA expression levels were significantly related to glioma OS and DFS. HDAC2 mRNA upregulation had prognostic value in OS. The foregoing results indicate that HDAC1/2 are potential prognostic biomarkers for patients with glioma. The present study may facilitate the discovery of novel prognostic biomarkers for glioma among other HDAC family members.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding authors.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Liuzhou Workers’ Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JL and HC participated in the manuscript of the article and the conception of the article, ML, XL, ZW, and JD edited the picture and data analysis, and QL participated in the review of the article. LC and LX for experimental verification and article revision.
Funding
This study was supported by the Science and Technology Department of Guangxi Zhuang Autonomous Region, Grant No. 2018GXNSFAA294026 and by Guangxi Natural Science Foundation, Grant No. 2018GXNSFAA294026.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
CNS, central nervous system; GBM, glioblastoma; IHC, immunohistochemistry; HDAC, histone deacetylase; GEPIA, gene expression profiling interactive analysis; TCGA, The Cancer Genome Atlas; CCLE, Cancer Cell Line Encyclopedia; KEGG, Kyoto Encyclopedia of Genes and Genomes; LGG, lower-grade glioma; TIMER, Tumor Immune Estimation Resource.
References
Abdelfattah, N., Rajamanickam, S., Panneerdoss, S., Timilsina, S., Yadav, P., Onyeagucha, B. C., et al. (2018). MiR-584-5p Potentiates Vincristine and Radiation Response by Inducing Spindle Defects and DNA Damage in Medulloblastoma. Nat. Commun. 9, 4541. doi:10.1038/s41467-018-06808-8
Asplund, A., Edqvist, P.-H. D., Schwenk, J. M., and Pontén, F. (2012). Antibodies for Profiling the Human Proteome-The Human Protein Atlas as a Resource for Cancer Research. Proteomics 12, 2067–2077. doi:10.1002/pmic.201100504
Auzmendi-Iriarte, J., Saenz-Antoñanzas, A., Andermatten, J., Elua-Pinin, A., Aldaba, E., Vara, Y., et al. (2019). P11.46 Discovery of a New HDAC6 Inhibitor for the Treatment of Glioblastoma. Neuro Oncol. 21, iii53–iii54. doi:10.1093/neuonc/noz126.192
Azagra, A., Román-González, L., Collazo, O., Rodríguez-Ubreva, J., de Yébenes, V. G., Barneda-Zahonero, B., et al. (2016). In Vivo conditional Deletion of HDAC7 Reveals its Requirement to Establish Proper B Lymphocyte Identity and Development. J. Exp. Med. 213, 2591–2601. doi:10.1084/jem.20150821
Bennani-Baiti, I. M., and Idriss, M. (2011). Epigenetic and Epigenomic Mechanisms Shape Sarcoma and Other Mesenchymal Tumor Pathogenesis. Epigenomics 3, 715–732. doi:10.2217/epi.11.93
Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data: Figure 1. Cancer Discov. 2, 401–404. doi:10.1158/2159-8290.cd-12-0095
Chauchereau, A., Mathieu, M., Saintignon, J. D., Ferreira, R., Pritchard, L. L., Mishal, Z., et al. (2004). HDAC4 Mediates Transcriptional Repression by the Acute Promyelocytic Leukaemia-Associated Protein PLZF. Oncogene 23 (54), 8777. doi:10.1038/sj.onc.120812
Cohen, T. J., Guo, J. L., Hurtado, D. E., Kwong, L. K., Mills, I. P., Trojanowski, J. Q., et al. (2011). The Acetylation of Tau Inhibits its Function and Promotes Pathological Tau Aggregation. Nat. Commun. 2, 252. doi:10.1038/ncomms1255
Cutano, V., Di Giorgio, E., Minisini, M., Picco, R., Dalla, E., and Brancolini, C. (2019). HDAC7-mediated Control of Tumour Microenvironment Maintains Proliferative and Stemness Competence of Human Mammary Epithelial Cells. Mol. Oncol. 13, 1651–1668. doi:10.1002/1878-0261.12503
Dang, Junsun. (2014). Image-Guided Synthesis Reveals Potent Blood-Brain Barrier Permeable Histone Deacetylase Inhibitors. Acs Chem. Neurosci. 5, 588.
Dompierre, J. P., Godin, J. D., Charrin, B. C., Cordelières, F. P., King, S. J., Humbert, S., et al. (2007). Histone Deacetylase 6 Inhibition Compensates for the Transport Deficit in Huntington's Disease by Increasing Tubulin Acetylation. J. Neurosci. 27, 3571–3583. doi:10.1523/jneurosci.0037-07.2007
Dovey, O. M., Foster, C. T., Conte, N., Edwards, S. A., Edwards, J. M., Singh, R., et al. (2013). Histone Deacetylase 1 and 2 Are Essential for normal T-Cell Development and Genomic Stability in Mice. Blood 121, 1335–1344. doi:10.1182/blood-2012-07-441949
Ghandi, M., Huang, F. W., Jané-Valbuena, J., Kryukov, G. V., Lo, C. C., McDonald, E. R., et al. (2019). Next-generation Characterization of the Cancer Cell Line Encyclopedia. Nature 569, 503–508. doi:10.1038/s41586-019-1186-3
Glozak, M. A., and Seto, E. (2007). Histone Deacetylases and Cancer. Oncogene 26, 5420–5432. doi:10.1038/sj.onc.1210610
Gutmann, D. H., Hedrick, N. M., Li, J., Nagarajan, R., Perry, A., and Watson, M. A. (2002). Comparative Gene Expression Profile Analysis of Neurofibromatosis 1-associated and Sporadic Pilocytic Astrocytomas. Cancer Res. 62, 2085–2091.
Haberland, M., Montgomery, R. L., and Olson, E. N. (2009). The many Roles of Histone Deacetylases in Development and Physiology: Implications for Disease and Therapy. Nat. Rev. Genet. 10, 32–42. doi:10.1038/nrg2485
Heideman, M. R., Wilting, R. H., Yanover, E., Velds, A., de Jong, J., Kerkhoven, R. M., et al. (2013). Dosage-dependent Tumor Suppression by Histone Deacetylases 1 and 2 through Regulation of C-Myc Collaborating Genes and P53 Function. Blood 121, 2038–2050. doi:10.1182/blood-2012-08-450916
Ho, T. C. S., Chan, A. H. Y., and Ganesan, A. (2020). Thirty Years of HDAC Inhibitors: 2020 Insight and Hindsight. J. Med. Chem. 63, 12460–12484. doi:10.1021/acs.jmedchem.0c00830
Hua, W. K., Qi, J., Cai, Q., Carnahan, E., Ramirez, M. A., Li, L., et al. (2017). HDAC8 Regulates Long-Term Hematopoietic Stem Cell Maintenance under Stress by Modulating P53 Activity. Blood 130 (24), 2619. doi:10.1182/blood-2017-03-771386
Huang, Y., Jian, W., Zhao, J., and Wang, G. (2018). Overexpression of HDAC9 Is Associated with Poor Prognosis and Tumor Progression of Breast Cancer in Chinese Females. Ott Vol. 11, 2177–2184. doi:10.2147/ott.s164583
Hubbert, C., Guardiola, A., Shao, R., Kawaguchi, Y., Ito, A., Nixon, A., et al. (2002). HDAC6 Is a Microtubule-Associated Deacetylase. Nature 417, 455–458. doi:10.1038/417455a
Jinwon, S., Min, S. K., Hye-Rim, P., Hoon, K. D., Jung, K. M., Su, K. L., et al. (2014). Expression of Histone Deacetylases HDAC1, HDAC2, HDAC3, and HDAC6 in Invasive Ductal Carcinomas of the Breast. J. Breast Cancer 17, 323.
Johannes, R., Emily, K., Kolbinger, F. R., Katharina, K., Siavosh, M., Lars, H., et al. (2018). Dual Role of HDAC10 in Lysosomal Exocytosis and DNA Repair Promotes Neuroblastoma Chemoresistance. Scientific Rep. 8, 10039.
Jones, P. A., and Baylin, S. B. (2007). The Epigenomics of Cancer. Cell 128, 683–692. doi:10.1016/j.cell.2007.01.029
Kang, Z.-H., Wang, C.-Y., Zhang, W.-L., Zhang, J.-T., Yuan, C.-H., Zhao, P.-W., et al. (2014). Histone Deacetylase HDAC4 Promotes Gastric Cancer SGC-7901 Cells Progression via P21 Repression. PLoS ONE 9, e98894. doi:10.1371/journal.pone.0098894
Kashyap, K., and Kakkar, R. (2020). Exploring Structural Requirements of Isoform Selective Histone Deacetylase Inhibitors: a Comparative In Silico Study. J. Biomol. Struct. Dyn. 39, 1–16. doi:10.1080/07391102.2019.1711191
Kikuchi, S., Suzuki, R., Ohguchi, H., Yoshida, Y., Lu, D., Cottini, F., et al. (2015). Class IIa HDAC Inhibition Enhances ER Stress-Mediated Cell. Death in Multiple Myeloma. Leukemia 29 (9), 1918. doi:10.1038/leu.2015.83
Kovacs, J. J., Murphy, P. J. M., Gaillard, S., Zhao, X., Wu, J.-T., Nicchitta, C. V., et al. (2005). HDAC6 Regulates Hsp90 Acetylation and Chaperone-dependent Activation of Glucocorticoid Receptor. Mol. Cel. 18, 601–607. doi:10.1016/j.molcel.2005.04.021
Krave, A. V. (2015). A Phase 2 Study of Concurrent Radiation Therapy, Temozolomide, and the Histone Deacetylase Inhibitor Valproic Acid for Patients with Glioblastoma. Int. J. Radiat. Oncol. Biol. Phys. 92, 986–992.
Li, H., Li, X., Lin, H., and Gong, J. (2020). High HDAC9 Is Associated with Poor Prognosis and Promotes Malignant Progression in Pancreatic Ductal Adenocarcinoma. Mol. Med. Rep. 21, 822–832. doi:10.3892/mmr.2019.10869
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77, e108–e110. doi:10.1158/0008-5472.can-17-0307
Li, Y., Peng, L., and Seto, E. (2015). Histone Deacetylase 10 Regulates the Cell Cycle G2/M Phase Transition via a Novel Let-7-HMGA2-Cyclin A2 Pathway. Mol. Cel. Biol. 35, 3547–3565. doi:10.1128/MCB.00400-15
Lienlaf, M., Perez-Villarroel, P., Knox, T., Pabon, M., Sahakian, E., Powers, J., et al. (2016). Essential Role of HDAC6 in the Regulation of PD-L1 in Melanoma. Mol. Oncol. 10, 735–750. doi:10.1016/j.molonc.2015.12.012
Linares, A., Assou, S., Lapierre, M., Thouennon, E., Duraffourd, C., Fromaget, C., et al. (2019). Increased Expression of the HDAC 9 Gene Is Associated with Antiestrogen Resistance of Breast Cancers. Mol. Oncol. 13, 1534–1547. doi:10.1002/1878-0261.12505
Lo Cascio, C., Mcnamara, J., Luna Melendez, E., and Mehta, S. (2019). STEM-25. HDAC1 IS ESSENTIAL FOR GLIOMA STEM CELL SURVIVAL. Neuro-Oncology 21, vi239. doi:10.1093/neuonc/noz175.998
Lodrini, M., Oehme, I., Schroeder, C., Milde, T., Schier, M. C., Kopp-Schneider, A., et al. (2013). MYCN and HDAC2 Cooperate to Repress miR-183 Signaling in Neuroblastoma. Nucleic Acids Res. 41, 6018–6033. doi:10.1093/nar/gkt346
Ma, P., Schultz, R. M., and Matzuk, M. M. (2013). Histone Deacetylase 2 (HDAC2) Regulates Chromosome Segregation and Kinetochore Function via H4K16 Deacetylation during Oocyte Maturation in Mouse. Plos Genet. 9, e1003377. doi:10.1371/journal.pgen.1003377
Max, F., Harold, R., Christian, L., Khalid, Z., Jason, M., Bader, G. D., et al. (2018). GeneMANIA Update 2018. Nucl. Acids Res. 46, W60–W64. doi:10.1093/nar/gky311
Mckinsey, T. A. (2011). The Biology and Therapeutic Implications of HDACs in the Heart, Histone Deacetylases: The Biology and Clinical Implication.
Millard, C. J., Watson, P. J., Fairall, L., and Schwabe, J. W. R. (2017). Targeting Class I Histone Deacetylases in a “Complex” Environment. Trends Pharmacol. Sci. 38, 363–377. doi:10.1016/j.tips.2016.12.006
Miller, K. M., Tjeertes, J. V., Coates, J., Legube, G., Polo, S. E., Britton, S., et al. (2010). Human HDAC1 and HDAC2 Function in the DNA-Damage Response to Promote DNA Nonhomologous End-Joining. Nat. Struct. Mol. Biol. 17, 1144–1151. doi:10.1038/nsmb.1899
Nakagawa, M., Oda, Y., Eguchi, T., Aishima, S.-I., Yao, T., Hosoi, F., et al. (2007). Expression Profile of Class I Histone Deacetylases in Human Cancer Tissues. Oncol. Rep. 18, 769–774. doi:10.3892/or.18.4.769
Neftel, C., Laffy, J., Filbin, M. G., Hara, T., Shore, M. E., Rahme, G. J., et al. (2019). An Integrative Model of Cellular States, Plasticity, and Genetics for Glioblastoma. Cell 178, 835–e21. doi:10.1016/j.cell.2019.06.024
Ning, Y., Ding, J., Sun, X., Xie, Y., Su, M., Ma, C., et al. (2020). HDAC9 Deficiency Promotes Tumor Progression by Decreasing the CD8+ Dendritic Cell Infiltration of the Tumor Microenvironment. J. Immunother. Cancer 8. doi:10.1136/jitc-2020-000529
Ostrom, Q. T., Gino, C., Haley, G., Nirav, P., Kristin, W., Carol, K., et al. (2019). CBTRUS Statistical Report: Primary Brain and Other Central Nervous System Tumors Diagnosed in the United States in 2012–2016. Neuro Oncol. 21, v1–v100. doi:10.1093/neuonc/noz150
Pacheco, M., and Nielsen, T. O. (2012). Histone Deacetylase 1 and 2 in Mesenchymal Tumors. Mod. Pathol. 25, 222–230. doi:10.1038/modpathol.2011.157
Park, H., Kim, J., Ahn, S., and Ryu, H. (2021). Epigenetic Targeting of Histone Deacetylases in Diagnostics and Treatment of Depression. Int. J. Mol. Sci. 22. doi:10.3390/ijms22105398
Philips, R. L., Lee, J. H., Gaonkar, K., Chanana, P., Chung, J. Y., Romero Arocha, S. R., et al. (2019). HDAC3 Restrains CD8-Lineage Genes to Maintain a Bi-potential State in CD4+CD8+ Thymocytes for CD4-Lineage Commitment. Elife 8. doi:10.7554/eLife.43821
Pomeroy, S. L., Tamayo, P., Gaasenbeek, M., Sturla, L. M., Angelo, M., Mclaughlin, M. E., et al. (2002). Prediction of central Nervous System Embryonal Tumour Outcome Based on Gene Expression. Nature 415, 436–442. doi:10.1038/415436a
Sandhu, S. K., Volinia, S., Costinean, S., Galasso, M., Neinast, R., Santhanam, R., et al. (2012). miR-155 Targets Histone Deacetylase 4 (HDAC4) and Impairs Transcriptional Activity of B-Cell Lymphoma 6 (BCL6) in the Eμ-miR-155 Transgenic Mouse Model. Proc. Natl. Acad. Sci. U S A. 109, 20047–20052. doi:10.1073/pnas.1213764109
Santoro, F., Botrugno, O. A., Dal Zuffo, R., Pallavicini, I., Matthews, G. M., Cluse, L., et al. (2013). A Dual Role for Hdac1: Oncosuppressor in Tumorigenesis, Oncogene in Tumor Maintenance. Blood 121, 3459–3468. doi:10.1182/blood-2012-10-461988
Sengupta, N., and Seto, E. (2004). Regulation of Histone Deacetylase Activities. J. Cel. Biochem. 93, 57–67. doi:10.1002/jcb.20179
Seto, E., and Yoshida, M. (2014). Erasers of Histone Acetylation: the Histone Deacetylase Enzymes. Cold Spring Harbor Perspect. Biol. 6, a018713. doi:10.1101/cshperspect.a018713
Simões-Pires, C., Zwick, V., Nurisso, A., Schenker, E., Carrupt, P.-A., Cuendet, M., et al. (2013). HDAC6 as a Target for Neurodegenerative Diseases: what Makes it Different from the Other HDACs. Mol. Neurodegeneration 8, 7. doi:10.1186/1750-1326-8-7
Song, Y., Luo, Q., Long, H., Hu, Z., Que, T., Zhang, X., et al. (2014). Alpha-enolase as a Potential Cancer Prognostic Marker Promotes Cell Growth, Migration, and Invasion in Glioma. Mol. Cancer 13, 65. doi:10.1186/1476-4598-13-65
Song, Y., Jiang, Y., Tao, D., Wang, Z., Wang, R., Wang, M., et al. (2020). NFAT2-HDAC1 Signaling Contributes to the Malignant Phenotype of Glioblastoma. Neuro-Oncology 22, 46–57. doi:10.1093/neuonc/noz136
Stengel, K. R., Zhao, Y., Klus, N. J., Kaiser, J. F., Gordy, L. E., Joyce, S., et al. (2015). Histone Deacetylase 3 Is Required for Efficient T Cell Development. Mol. Cel. Biol. 35, 3854–3865. doi:10.1128/mcb.00706-15
Stengel, K. R., Bhaskara, S., Jing, W., Qi, L., Ellis, J. D., Sampathi, S., et al. (2019). Histone Deacetylase 3 Controls a Transcriptional Network Required for B Cell. Maturation. Nucleic Acids Res. 47 (20), 10612–10627. doi:10.1093/nar/gkz816
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING V11: Protein-Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-wide Experimental Datasets. Nucleic Acids Res. 47, D607–D613. doi:10.1093/nar/gky1131
Stupp, R., Hegi, M. E., Mason, W. P., van den Bent, M. J., Taphoorn, M. J., Janzer, R. C., et al. (2009). Effects of Radiotherapy with Concomitant and Adjuvant Temozolomide versus Radiotherapy Alone on Survival in Glioblastoma in a Randomised Phase III Study: 5-year Analysis of the EORTC-NCIC Trial. Lancet Oncol. 10, 459–466. doi:10.1016/S1470-2045(09)70025-7
Stupp, R., Mason, W. P., van den Bent, M. J., Weller, M., Fisher, B., and Taphoorn, M. J. B. (2005). Radiotherapy Plus Concomitant and Adjuvant Temozolomide for Glioblastoma. New Engl. J. Med. 352, 987–996. doi:10.1056/NEJMoa043330
Sun, L., Hui, A.-M., Su, Q., Vortmeyer, A., Kotliarov, Y., Pastorino, S., et al. (2006). Neuronal and Glioma-Derived Stem Cell Factor Induces Angiogenesis within the Brain. Cancer Cell 9, 287–300. doi:10.1016/j.ccr.2006.03.003
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: A Web Server for Cancer and Normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res. 45, W98–W102. doi:10.1093/nar/gkx247
Thole, T. M., Lodrini, M., Fabian, J., Wuenschel, J., Pfeil, S., Hielscher, T., et al. (2017). Neuroblastoma Cells Depend on HDAC11 for Mitotic Cell Cycle Progression and Survival. Cell Death Dis. 8, e2635. doi:10.1038/cddis.2017.49
Thurn, K. T., Thomas, S., Raha, P., Qureshi, I., and Munster, P. N. (2013). Histone Deacetylase Regulation of ATM-Mediated DNA Damage Signaling. Mol. Cancer Ther. 12, 2078–2087. doi:10.1158/1535-7163.mct-12-1242
Trzeciakiewicz, H., Ajit, D., Tseng, J.-H., Chen, Y., Ajit, A., Tabassum, Z., et al. (2020). An HDAC6-dependent Surveillance Mechanism Suppresses Tau-Mediated Neurodegeneration and Cognitive Decline. Nat. Commun. 11, 5522. doi:10.1038/s41467-020-19317-4
Tsai, C., Ko, H., Chiou, S., Lai, Y., Hou, C., Javaria, T., et al. (2021). NBM-BMX, an HDAC8 Inhibitor, Overcomes Temozolomide Resistance in Glioblastoma Multiforme by Downregulating the β-Catenin/c-Myc/SOX2 Pathway and Upregulating P53-Mediated MGMT Inhibition. Int. J. Mol. Sci. 22. doi:10.3390/ijms22115907
Wang, X.-Q., Bai, H.-M., Li, S.-T., Sun, H., Min, L.-Z., Tao, B.-B., et al. (2017). Knockdown of HDAC1 Expression Suppresses Invasion and Induces Apoptosis in Glioma Cells. Oncotarget 8, 48027–48040. doi:10.18632/oncotarget.18227
Wang, Y., Zhang, J., Wu, L., Liu, W., Wei, G., Gong, X., et al. (2018). Tricho-rhino-phalangeal Syndrome 1 Protein Functions as a Scaffold Required for Ubiquitin-specific Protease 4-directed Histone Deacetylase 2 De-ubiquitination and Tumor Growth. Breast Cancer Res. 20, 83. doi:10.1186/s13058-018-1018-7
Wang, Z., Qin, G., and Zhao, T. C. (2014). HDAC4: Mechanism of Regulation and Biological Functions. Epigenomics 6, 139–150. doi:10.2217/epi.13.73
Wilson, A. J., Byun, D.-S., Nasser, S., Murray, L. B., Ayyanar, K., Arango, D., et al. (2008). HDAC4 Promotes Growth of Colon Cancer Cells via Repression of P21. MBoC 19, 4062–4075. doi:10.1091/mbc.e08-02-0139
Yang, J., Gong, C., Ke, Q., Fang, Z., Chen, X., Ye, M., et al. (2021). Insights into the Function and Clinical Application of HDAC5 in Cancer Management. Front. Oncol. 11, 661620. doi:10.3389/fonc.2021.661620
Yue, L., Sharma, V., Horvat, N. P., Akuffo, A. A., Beatty, M. S., Murdun, C., et al. (2020). HDAC11 Deficiency Disrupts Oncogene-Induced Hematopoiesis in Myeloproliferative Neoplasms. Blood 135, 191–207. doi:10.1182/blood.2019895326
Zhang, X., Yuan, Z., Zhang, Y., Yong, S., Salas-Burgos, A., Koomen, J., et al. (2007). HDAC6 Modulates Cell Motility by Altering the Acetylation Level of Cortactin. Mol. Cel. 27, 197–213. doi:10.1016/j.molcel.2007.05.033
Zhao, Z., Zhang, K., Wang, Q., Li, G., and Jiang, T. (2020). Chinese Glioma Genome Atlas (CGGA): A Comprehensive Resource with Functional Genomic Data for Chinese Glioma Patients.
Zhu, P., Martin, E., Mengwasser, J., Schlag, P., Janssen, K.-P., and Göttlicher, M. (2004). Induction of HDAC2 Expression upon Loss of APC in Colorectal Tumorigenesis. Cancer Cell 5, 455–463. doi:10.1016/s1535-6108(04)00114-x
Keywords: bioinformatics analysis, HDACs, immune infiltration, CGGA, Chinese Glioma Genome Atlas, biomarker, glioma
Citation: Li J, Yan X, Liang C, Chen H, Liu M, Wu Z, Zheng J, Dang J, La X and Liu Q (2022) Comprehensive Analysis of the Differential Expression and Prognostic Value of Histone Deacetylases in Glioma. Front. Cell Dev. Biol. 10:840759. doi: 10.3389/fcell.2022.840759
Received: 21 December 2021; Accepted: 31 January 2022;
Published: 10 March 2022.
Edited by:
Ângela Sousa, University of Beira Interior, PortugalReviewed by:
Vijay Thakur, University of Miami Hospital, United StatesYinyan Wang, Capital Medical University, China
Copyright © 2022 Li, Yan, Liang, Chen, Liu, Wu, Zheng, Dang, La and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Quan Liu, liuq0928@163.com; Xiaojin La, laxiaojin@ncst.edu.cn