Skip to main content

Systematic analysis of histone acetylation regulators across human cancers

Abstract

Background

Histone acetylation (HA) is an important and common epigenetic pathway, which could be hijacked by tumor cells during carcinogenesis and cancer progression. However, the important role of HA across human cancers remains elusive.

Methods

In this study, we performed a comprehensive analysis at multiple levels, aiming to systematically describe the molecular characteristics and clinical relevance of HA regulators in more than 10000 tumor samples representing 33 cancer types.

Results

We found a highly heterogeneous genetic alteration landscape of HA regulators across different human cancer types. CNV alteration may be one of the major mechanisms leading to the expression perturbations in HA regulators. Furthermore, expression perturbations of HA regulators correlated with the activity of multiple hallmark oncogenic pathways. HA regulators were found to be potentially useful for the prognostic stratification of kidney renal clear cell carcinoma (KIRC). Additionally, we identified HDAC3 as a potential oncogene in lung adenocarcinoma (LUAD).

Conclusion

Overall, our results highlights the importance of HA regulators in cancer development, which may contribute to the development of clinical strategies for cancer treatment.

Peer Review reports

Introduction

Cancer is a major public health problem in the whole world, with high global morbidity and mortality, causing enormous economic burden and distracting human development. Beyond genome alternations were initially thought as the driver of carcinogenesis, cancer mechanisms have expended to a complex roster including genomic and non-genomic derangements [1]. Epigenetic modifications are functioned as heritable and reversible alternations in genes or genomes without any changes in primitive DNA sequence, and thus act as a vital role in modulation of gene transcription. Emerging evidence highlights that epigenomic dynamics are implicated in regulation of cancers initiation, progression, metastasis and even impacting status of immune cells in the context of cancer immunotherapy[2]. Regard this as cutting point, exploration of epigenetic adaptation in cancer transformation and epidrugs discovery is a burgeoning field.

Epigenetic regulation, including DNA and histone modifications, non-coding RNA modulation and chromatin remodeling, exerts epigenomic imprinting on genomic DNA and endows general flexibility and plasticity on gene expression. As the decisive advances in decrypting mechanistic of cancerization processes, the cellular memory of epigenetic modifications is lost or perturbated in the circumstances of tumor occurrence and development. Histone acetylation (HA), an important and common epigenetic pathway, regulates gene transcription at relevant genomic sites by controlling the tension of nucleosome structure. Thus, there is no doubt that this epigenetic mechanism could be hijacked by tumor cells during carcinogenesis and cancer progression. As the actual fact is, growing evidence shows that the perturbations on histone acetylation are commonly observed in the cases of tumor occurrence [3,4,5].The “eraser” of HA HDAC1 has been reported to be highly expressed in non-small cell lung cancer (NSCLC) cell lines, and down-regulation of HDAC1 inhibits lung cancer cell proliferation, migration, invasion, tumor angiogenesis, and induces cell apoptosis [6]. And the up-regulation of HDAC2 in breast cancer cells and tissues can also regulate the malignant biological behavior of breast cancer cells. In addition, HA dynamics also regulates immune-related processes under tumor microenvironment and could be harnessed stochastically for evasion of immunosurveillance [7]. In addition, HA dynamics also regulates immune-related processes under tumor microenvironment and could be harnessed stochastically for evasion of immunosurveillance [8]. These findings suggest that HA regulators are involved in important biological processes in cancer.

Studies on histone modifications have never stopped. However, the important role of HA across human cancers remains elusive. Decrypting the landscape and gaining insights of the underlying mechanisms of HA could help delineate alternative and preventive strategies for clinical cancer therapy. Thus, in this study, we performed a comprehensive analysis at multiple levels, aiming to systematically describe the molecular characteristics and clinical relevance of HA regulators in 33 cancer types. We found that there were extensive genetic alterations in HA regulators in human cancers. We also evaluated the correlation of expression perturbations of HA regulators with the activity of cancer pathways, and explored the clinical prognostic value of HA regulators. Overall, our comprehensive analysis of HA regulators provides important resources for understanding HA biology.

Materials and methods

Data source and preprocessing

Histone acetylation (HA) regulators (Table S1) were identified from a previous report [9]. The FPKM-based mRNA expression data, TCGA thresholded SCNA scores (ISAR_GISTIC.all_data_by_genes.txt), somatic mutation data, and clinical data were download from the Genomic Data Commons (https://portal.gdc.cancer.gov/), Broad GDAC Firehose (https://gdac.broadinstitute.org/) and the Xena Browser (https://xenabrowser.net/datapages/). The GDC.h38 GENCODE v22 GTF (gencode.v22.annotation.gtf.gz) (https://gdc.cancer.gov/download-gdc-reference-files) was used as the genome annotation file.

Differential expression analysis and Pearson correlation analysis

From 33 cancer types, 19 cancer types with at least 3 matched tumors and normal samples were selected for differential expression analysis. Differentially expressed genes (DEGs) were defined as described by Li et al.[10]. The Wilcox’s test was used to identify DEGs with adjusted p-values by BH method. The threshold was set as the adjusted p-value less than 0.05. The Pearson correlation between somatic copy number alterations and the expression of HA regulators were investigated as described by Liu et al.[11]. The correlation with p-value less than 0.05 was considered as significant or considerable.

Hallmark pathway activity across cancer types

Gene Set Variation Analysis (GSVA) is a non-parametric unsupervised analysis method, which is mainly used to evaluate whether different metabolic pathways are enriched between different samples by converting the expression matrix of genes between different samples into the expression matrix of gene sets between samples. Unlike gene-set enrichment analysis (GSEA), GSVA does not require pre-grouping of samples and can calculate enrichment scores for specific gene sets in each sample. Here, GSVA was performed on the normalized gene expression to calculate the hallmark pathway activity with GSVA scores. To identify the HA regulators associated with pathway activation or inhibition, we calculated Pearson correlation coefficient (PCC) between the HA regulators expression and pathway activity. The HA regulators with significant correlations were defined as those with absolute value of regulator-pathway pair correlation coefficients greater than 0.2 and adjusted p-value less than 0.001.

Cross talks among histone acetylation (HA) regulators

Based on the gene expression of 33 cancer types, we calculated the PCC among HA regulators and visualized them with “corrplot” R package. Furthermore, we also investigated the protein-protein interaction (PPI) networks among these HA regulators through the String database (https://cn.string-db.org/) (Browsing date: October 7, 2022). This interaction was subsequently further visualized in the Cytoscape_v3.9.0. software.

Survival analysis and CellMiner analysis for HA regulators

To investigate the correlation between the expression of HA regulators and patient survival, we divided all patients into two groups according to the median expression of each HA regulator. Cox regression was used to examine the differences in survival between the two groups. In this process, the “survival” package in R was adopted. In addition, in order to screen potential therapeutic agents targeting HA regulators, we downloaded drug information and sequencing data from the CellMiner database (https://discover.nci.nih.gov/cellminer/home.do). Here, we calculated the PCC between the expression level of each HA regulator and drug half inhibitory concentration (IC50).

Unsupervised consensus clustering

Based on the mRNA expression of HA regulators in KIRC tumor samples, we performed unsupervised consensus clustering to identify different patient clusters. The “ConsensuClusterPlus” package in R was adopted. In this process, 1000-time repetitions were performed to guarantee the classification stability. The specific setting parameters were shown as follows: maxK = 5, reps = 1000, pItem = 0.8, pFeature = 1, seed = 73, clusterAlg = “km”, distance = “euclidean”.

Cell culture

A549, H1299 and Beas-2B cells were all derived from the ATCC. A549 cells were cultured in F-12k medium, H1299 in RPMI 1640 medium, and human normal bronchial epithelial cells Beas-2B in DMEM medium. In these media, we supplemented them with 10% fetal bovine serum (FBS) and 1% streptomycin and 1% penicillin. The cells were cultured in a cell incubator at 37℃ and 5%CO2.

Cell transfection

The cell density of 2.5 × 105 cells per well was inoculated in 6-well plate, and siRNA transfection was performed after A549 cells were attached to the wall at 70–90% confluent. Transfection experiment was divided into four groups: negative control group (si-NC), si-HDAC3-001, si-HDAC3-002 and si-HDAC3-003. The si-HDAC3-001 primer sequence was 5’-GAGCAACCCAGCTGAACAA-3’, si-HDAC3-002 primer sequence was 5’-GTCCTGCATTACGGTCTCT-3’, and si-HDAC3-003 primer sequence was 5’-GTGGTTATACTGTCCGAAA-3’. The transfection procedure was carried out according to LipofectamineTM2000 instructions.

Real time quantitative PCR (RT-qPCR)

Total RNA extraction of cells was performed by TRIpure Total RNA Extraction Reagent (Biosharp, China). RNA was subsequently reverse transcribed to cDNA according to the cDNA Synthesis Kit (Servicebio, Wuhan, China). HDAC3 primer sequence: S: 5’-TACGGAGCTGGACACCCTATG-3’, A: 5’-ATGTAGTCCTCGGAGTGGAAGC-3’. GAPDH primer sequence: S: 5’-GGAAGCTTGTCATCAATGGAAATC-3’, A: 5’-TGATGACCCTTTTGGCTCCC-3’. RT-qPCR was performed using validated specific primers and SYBR Green PCR SuperMix (Servicebio, China). Meanwhile, GAPDH was used as an internal reference to calculate the mRNA content of target gene by 2−ΔΔCt method.

Western blots

Cells were treated with RIPA cell lysis buffer containing PMSF. The protein concentration was determined by BCA method, and the protein solution and 5×Loading Buffer were subsequently mixed according to the volume ratio of 41, and boiled at 100 ℃ for 10 min. After SDS-PAGE gel, 50 µg of total protein samples were added to each sample well for electrophoresis for 1.5 h and a voltage of 100 V. The membrane was transferred at 4 ℃ and 250 mA for 2 h. NC membranes were immersed in 5% skim milk for 1.5 h at room temperature. The NC membranes were then placed within the diluted primary antibody at 4℃ overnight. After this, the NC membranes were placed within the diluted secondary antibody (1:10000) for 2 h. HDAC3, E-cadherin, Vimentin, N-cadherin and GAPDH were diluted at 1: 600, 1 : 1000, 1: 1000, 1: 1000 and 1: 1000 respectively. ECL chemiluminescence solution was used for protein exposure.

CCK-8 assay

A549 cells in the logarithmic growth phase were digested with trypsin (Servicebio, Wuhan, China). Then cells were collected, centrifuged and resuspended with F-12k culture medium. The cell suspension were seeded in 96-well plates at 8000 cells/well, with each group repeated three wells and cultured for 24 h. The RNA transfection was carried out according to LipofectamineTM2000 instructions. Cells were cultured for 24 h, the CCK-8 solution was added at 10 µl/well, and the cells were subsequently incubated in an incubator for 2 h at 37℃ and 5%CO2. The cell proliferation rate was calculated by measuring the OD of each well at 450 nm.

EDU assay

At 48 h after transfection of A549 cells, cell suspension was configured after trypsin digestion and seeded in 6-well plates at a density of 1*105 cells/well and cultured overnight. They were supplemented with 10 µmol/L EdU for further 2 h. Cells were fixed with 4% paraformaldehyde for 30 min and permeabilized with 0.5%Triton X-100 permeabilized solution for 20 min at room temperature. Cells were stained with Apollo staining solution for 30 min and DAPI for 5 min. Photo was taken under a fluorescence microscope. DAPI positive cells (blue) are the total cell number and EdU positive cells (red) are the proliferating cells.

Wound healing assay

Three repeated wells were set for the experimental and control groups. A549 cells from the logarithmic growth phase were seeded at 500 µL per well (2*105/mL) in 6-well plates. Then each well was added with 2 mL of F-12k complete medium and cells were placed in an incubator at 37℃ and 5%CO2. After the cells covered the bottom of the plate, we used a 1 ml pipette tip to make scratches on the well plate. The scratch widths at 0 and 24 h were subsequently observed and photographed.

Transwell assay

The BD matrigel and F-12k were diluted in a 1:3 ratio, then 50ul was absorbed into the transwell upper chamber and placed in an incubator for about 4 h. A549 cells in the experimental and control groups were seeded in a small Transwell chamber with 0.2 mL per well (2*105/mL), with three repeated wells for each group. After adding 500 µl of complete medium to the Transwell plate, the chamber with A549 cells was put into the plate and cultured in a 37℃ incubator for 24 h. The medium in the chamber was then removed and rinsed with PBS and stained with crystal violet for 10 min. The redundant crystal violet on the surface of the chamber was cleaned with running water, and the cells in the upper chamber were wiped clean with cotton swabs. The non-cell inoculant side was photographed, and three fields were randomly selected under an inverted microscope (×100) to calculate the number of cells.

Bioinformatic analysis of HDAC 3 expression differences in LUAD

Two LUAD chips (GSE10072 [12], and GSE32863 [13]) were downloaded from Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) to obtain gene expression data for HDAC3, and we divided all patients into two groups according to the median expression of HDAC3. Statistical differences between the two groups were calculated using Wilcoxon rank sum test. The data was visualized with the “ggplot2” R package [14]. GSE10072 dataset included transcription profile data of 49 normal lung tissues and 58 lung tumor tissues. And GSE32863 dataset included transcription profiling data from 58 normal lung tissues and 58 lung tumor tissues. The difference in protein expression of HDAC3 between LUAD tumors and normal tissues was explored using the UALCAN database (http://ualcan.path.uab.edu/analysis-prot.html). To further verify the differential expression of HDAC3 protein in LUAD, immunohistochemical sections from LUAD tumor and normal tissues were also obtained from the Human Protein Atlas (https://www.proteinatlas.org).

Statistical analysis

Gene differential expression analysis was performed by the “limma” R package. The Wilcoxon rank-sum test was used to calculate the statistical differences between the two groups. Differences in survival between two or more groups was compared using the Kaplan-Meier method. A P-value of < 0.05 was considered statistically significant.

Results

Genetic alterations in HA regulators in human different cancer types

We identified 36 acknowledged histone acetylation (HA) regulators, including 9 (HAT1, KAT2A, KAT2B, KAT8, KAT6A, KAT6B, KAT7, EP300 and CREBBP) in “writers”, 12 (HDAC1, HDAC2, HDAC3, HDAC8, SIRT2, HDAC4, HDAC5, HDAC6, HDAC7, HDAC9, HDAC10, and HDAC11) in “erasers” and 15 (BRD2, BRD3, BRD4, BRDT, BPTF, ATAD2B, BAZ2B, TAF1, YEATS4, DPF3, SMARCA2, SMARCA4, PBRM1, DPF1, and DPF2) in “readers” (Fig. 1A, and Table S1). In Fig. 1B, we initially characterized the key players of these three classes of HA regulators. Generally, histone acetylation is a post-translational modification dynamically regulated by three categories of regulators, namely “writers”, “readers” and “erasers”, to maintain this trimmed process to functional execution. The “writers” transfer acetyl to N-terminal lysine residue of histones in the chromatin to neutralize the electrostatic attraction between DNA chain and histones, thus loose the nucleosome structure to promote chromatin accessibility and specific genes transcription, whereas the “erasers” catalyze the converse process. The “readers”, regarded as effector proteins, refer to recognized site-specific histones acetylation modifications and make it interpretable as well as to perform diverse functions.

Fig. 1
figure 1

Histone acetylation (HA) regulators and their genetic alternations in human cancers. (A) The list of readers, erasers and writers among histone acetylation regulators. (B) The diagram of histone acetylation and deacetylation processes and associated regulators. (C) The mutation frequency of HA regulators across human cancers. (D) The CNV alternation frequency of HA regulators across human cancers

Subsequently, by deeply exploring the patterns of genomic alterations across 33 cancer types, we found low overall levels of mutation frequencies in HA regulators, with the vast majority of HA regulators having less than 5% in the vast majority of cancer types (Fig. 1C). The mutation frequency of HA regulators varies with cancer contexts and regulators. In CHOL, KICH, MESO, PCPG, TGCT, THYM, and UVM, about three quarters of the regulators do not have any mutations. The “reader” PBRM1 has a high mutation frequency in KIRC and CHOL, at 42.42% and 22.22%, respectively, without any mutations in PCPG, THCA and UVM. We found that the overall mutation frequency of HA regulators was higher in UCEC and SKCM than in other cancers, which is consistent with a previous report [10].

We noticed that HA regulators exhibits a cancer-type-dependent pattern of CNV amplification or deep deletion (Fig. 1D). For example, in KICH, the vast majority of HA regulators did not show CNV amplification or deep deletion, but the opposite occurred in OV. Furthermore, HDAC3 had relatively high CNV amplification (22.36%) and lacked deep deletions in KIRC. However, KAT2B, HDAC11, and PBRM1 had a relatively high CNV deep deletion and lacked amplification in KIRC. The genomic alterations of HA regulators in UCEC are dominated by somatic mutations, while in OV, CNV amplification dominated them. Collectively, these results (Fig. 1B-D and Table S2) reveal a highly heterogeneous genetic alteration landscape of HA regulators across different human cancer types.

Aberrant expression of HA regulators among cancer types

We previously investigated the genetic alteration of HA regulators in human different cancer types. To make as clear as possible the potential associations between these alterations and HA regulators expression, we therefore performed differential expression analysis to explore the expression perturbations of HA regulators in 19 cancer types (with at least three normal control samples), and Pearson correlation analysis to investigate the associations between somatic CNV and the expression of HA regulators. We found that the vast majority of HA regulators are differentially expressed in most cancer types. Several HA regulators showed consistent expression patterns in cross-cancer analyses. HDAC3, HDAC10, HDAC2, HDAC7, ATAD2B were significantly upregulated in 9, 14, 13, 9, and 10 cancer types, respectively, whereas KAT2B was downregulated in 16 cancer types (Fig. 2A and Table S3). In addition, the expression patterns of some HA regulators are cancer-type-dependent. For example, HDAC3 was upregulated in most cancer types, including LUAD (log2FC = 0.18, adj.P = 2.94E-09) and LIHC (log2FC = 0.26, adj.P = 1.36E-07), but was down-regulated in KICH (log2FC = -0.21, adj.P = 2.53E-04) and THCA (log2FC = -0.07, adj.P = 1.15E-04). This suggests that HA regulators may play different roles in different cancer types. Our results suggest that CNV alteration may be one of the major mechanisms leading to the expression perturbations in HA regulators. Compared with normal tissues, the expression of CNV-amplified HA regulators was significantly increased in tumor tissues (e.g. HDAC3 in KIRC), while the expression of CNV-deficient HA regulators was significantly decreased (e.g. KAT2B, HDAC11 and PBRM1 in KIRC). Further, we investigated the Pearson correlation between HA regulators gene expression and copy number. We found that the expression of most HA regulators was significantly positively correlated with copy number, except for BRDT in HNSC, HDAC8 in THCA and UVM (Fig. 2B and Table S4). These results suggest that the copy number of HA regulator is abnormal in most cancer types and can affect gene expression. This highlights the importance of the dysregulation of HA regulators in different cancers.

Fig. 2
figure 2

The dysregulated gene expression of HA regulators across human cancers. (A) The heatmap showing different expression alterations of HA regulators in human cancers. (B) The Pearson correlation between somatic copy number alterations and the gene expression of HA regulators

Association between HA regulators and Hallmark pathways in Cancer

To further elucidate the molecular mechanisms underlying the involvement of HA regulators in cancer, we calculated the correlation between the expression of individual HA regulators and the activity of 50 Hallmark pathways. We observed that the expression of each HA regulator was associated with the activation or inhibition of multiple pathways (Fig. 3A and Table S5). For example, HDAC3 was associated with the activity of 46 pathways, of which 23 were positively related (including E2F Targets, G2M Checkpoint, MYC Targets V1, Epithelial Mesenchymal Transition, PI3K/AKT/mTOR signaling), and 23 were negatively related (including Bile acid metabolism, Protein Secretion). Further reinforcing this correlation (|PCC|>0.2 and adj.P < 0.001), we found that HDAC11 was associated with the most pathway activity, which was associated with inhibition of 9 pathways and activation of 15 pathways (Fig. 3B and Table S5). HA regulators of the same functional category are associated with different pathway activities, and HA regulators of different functional categories (“readers”, “writers” or “erasers”) are associated with the same pathway activity, indicating that there is functional crosstalk among HA regulators, that is, HA regulators in the same functional category have different functional roles.

Fig. 3
figure 3

Association of the HA regulators with Hallmark oncogenic pathways. (A) Network diagram showing the correlation between the expression of HA regulators and the activity of cancer pathways. (B) The number of pathways involved by each HA regulator. (C) Pearson correlation among HA regulators. (D) The protein-protein interaction (PPI) networks among these HA regulators

Previous studies [15, 16]have shown close collaboration among “readers”, “writers” and “erasers” belonging to different functional categories in cancer. Our results also preliminarily demonstrate this point. We found that not only HA regulators within the same functional class showed significantly correlated expression patterns (e.g. CREBBP-KAT6A, CREBBP-EP300), but also high correlations between different functional classes (e.g. CREBBP-PBRM1, CREBBP-BRD3, EP300-BRD3) (Fig. 3C). Furthermore, we also noted extensive interactions of these regulators in the protein-protein interaction (PPI) network (Fig. 3D). Taken together, these results suggest that cross-talk among “readers”, “writers” and “erasers” of HA plays key roles in different cancer types.

Clinical relevance of HA regulators

Given the highly heterogeneous genetic and expression alterations of HA regulators in cancer, it is warranted to further investigate the clinical relevance of the HA regulators. First, we investigated the relationship between HA regulators and patient survival (Fig. 4A). We found a significant cancer-specific relationship between HA regulators and patients’ prognosis. For example, HDAC3 is a prognostic unfavorable factor in PAAD and READ, but the opposite factor occurs in LIHC and LGG. Regulators in the same functional category in the same cancer type may also present radically different prognostic values. For example, HDAC2 and HDAC9, both “readers”, had opposite effects on survival in ACC patients. Subsequently, to screen for potential drugs targeting HA regulators, we downloaded drug information and RNA-seq data from the CellMiner database for drug sensitivity analysis. The results showed that HDAC11, HDAC7, BRD 3, YEATS4, and HDAC3 correspond to 6, 6, 3, 3, and 1 drugs, respectively (Fig. 4B and Table S6). For example, the expression of HDAC7 was negatively correlated with the IC50 of selmetinib and cobimetinib, suggesting that its high expression may enhance the drug sensitivity of selmetinib and cobimetinib.

Fig. 4
figure 4

Clinical relevance of HA regulators. (A) Prognostic predictive value of individual HA regulators for patients with different cancer types. (B) The correlation between the expression of HA regulators and drug IC50. Red lines represent positive correlation and blue lines represent negative correlation

Among these 33 cancer types, we observed that KIRC had the largest number of prognosis-related regulators (20/36). In order to explore whether the expression of HA regulators contributes to the prognostic stratification of KIRC, we performed unsupervised clustering based on the overall expression pattern of HA regulators and finally identified three patient clinical clusters. We called them HAsCluster A, HAsCluster B, and HAsCluster C. They contained 175, 177, and 183 patients, respectively. Figure 5 A and Fig. 5B show the overall expression landscape of HA regulators among three clinical clusters in heatmap and bar-graphs, respectively. We found that the vast majority of the HA regulators had the lowest expression in HAsCluster B. We subsequently explored survival differences among the three clusters (Fig. 5C). The results showed that HAsCluster C had significantly better clinical outcomes, no matter in OS, DSS or PFI. Overall, our findings demonstrate the great potential of HA regulators both in the generation of novel therapeutic strategies and in the prognostic stratification of specific cancer types.

Fig. 5
figure 5

The unsupervised clustering identified three patient clinical clusters for KIRC. (A) The heatmap showing the clustering for KIRC patients based on the overall expression pattern of HA regulators. (B) Differential expression of HA regulators among the three clusters. (C) The Kaplan-Meier survival plots showing the survival differences among the three clusters

Validation of the oncogenic role of HDAC3 in LUAD

Histone deacetylases (HDACs) are cellular enzymes that act a crucial role in epigenetic regulation of gene expression by remodeling chromatin [17, 18]. HDAC3, an important member of class I HDACs functioning in many key cellular processes, has recently been found to be overexpressed in multiple malignancies and is considered to be one of the most frequently upregulated genes in malignancies[19, 20]. Considering that the critical functions of HDAC3 in human cancers including LUAD are still not fully interpreted and the important findings were observed in this study, we further revealed the important role of HDAC3 in LUAD through cell function experiments. HDAC3 was highly expressed in multiple cancer types, including LUAD. We also observed a similar situation in two LUAD GEO chip datasets (GSE10072 and GSE32826) (Fig. 6A). This expression pattern was also validated at the protein level (Fig. 6B-C). Subsequently, we also explored at the cellular level and found that HDAC3 was more highly expressed in two LUAD cell lines (A549 and H1299) compared to normal bronchial epithelial cells (Beas-2B) (Fig. 6D). We interfered the expression of HDAC3 with siRNA to further investigate the effect of HDAC3 on the malignant biological behavior of cells. We found that RNA interference significantly reduced HDAC3 expression (Fig. 6E-F). The absorbance at 450 nm represents the number of cells, reflecting the capacity of cell proliferation. According to the results of CCK-8 assay, the cell proliferation rate of si-HDAC3 group was significantly lower than that of si-NC group (Fig. 6G). Also in the EdU assay, the number of EdU-positive cells in the si-HDAC3 group was significantly lower than that in the si-NC group after RNA interference treatment of A549 for 24 h (Fig. 6H). In the subsequent wound healing assay, we noted that after RNA interference treatment of A549 for 24 h, the scratch area of si-HDAC3 group was significantly higher than that of si-NC group (Fig. 6I), indicating that knockdown of HDAC3 could inhibit cell migration of A549. According to the results of Transwell invasion assay, we also found that the invasion ability of A549 was significantly decreased after RNA interference treatment (Fig. 6J). In the investigation results mentioned above, in addition to the observation of high HDAC3 expression in LUAD, we also found that HDAC3 was associated with the activation or inhibition of multiple Hallmark pathways. Our bioinformatic results initially suggested that HDAC3 might participate in the epithelial mesenchymal transition process, and this conclusion was also confirmed by in in vitro experiments. In western blot, 24 h after RNA interference treatment of A549, it was observed that the expression levels of E-cadherin in si-HDAC3 group were significantly higher than those in si-NC group, while the expression levels of N-cadherin and Vimentin were significantly lower than those in si-NC group (Fig. 6K). This suggests that knockdown of HDAC3 inhibits epithelial mesenchymal transformation in human LUAD cells. Collectively, these results suggested that HDAC3 played an oncogenic role in LUAD.

Fig. 6
figure 6

Validation of the oncogenic role of HDAC3 in LUAD. (A) Differential mRNA expression of HDAC3 between LUAD tumor and normal tissues: results from two GEO chips (GSE10072 and GSE32826). (B) Differential protein expression of HDAC3 between LUAD tumor and normal tissues: results from CPTAC samples. (C) Differential protein expression of HDAC3 between LUAD tumor and normal tissues: results from HPA. (D) Differential expression of HDAC3 between LUAD cells (A549 and H1299) and normal bronchial epithelial cells (Beas-2B). (E) qPCR validation of RNA interference effects. (F) Western blot validation of RNA interference effects. (G) CCK-8 assay reveals the cell proliferation rate of si-HDAC3 group and si-NC group. (H) Edu assay reveals the cell proliferation of si-HDAC3 group and si-NC group. (I) Woundhealing assay reaveals the cell migration of si-HDAC3 group and si-NC group. (J) Transwell assay reveals the cell invasion of si-HDAC3 group and si-NC group. (K) Western blot reveals the expression levels of epithelial mesenchymal transition related proteins (E-cadherin, N-cadherin, Vimentin) in different groups of cells

Discussion

HA is a post-translational modification that controls the transcript level of genes and plays a key role in the structural modification of chromosomes and the regulation of gene expression [21, 22]. Increasing evidence suggests that HA is closely related to tumor occurrence and development [23, 24]. For example, histone deacetylases can affect gene expression by direct association with transcription factors, which are deacetylated in malignancy [25]. Bianco-Miotto T et al. found that in prostate cancer, specific histone modifications could predict the progression of cancer [26]. Similarly, in breast cancer, Elsheikh S et al. also observed the close relationship between global histone modification and tumor phenotype and patient prognosis [27]. Molecular drugs targeting HA regulators (e.g. HDAC inhibitors) have been shown to have effective anti-tumor effects in haematological and solid malignancies [28, 29]. With the deepening of research, HA is also considered to be closely related to the tumor immune microenvironment [30]. HDAC inhibitors can reshape the tumor microenvironment through a variety of mechanisms, thereby enhancing the ability of immune surveillance and killing tumor cells. These studies have highlighted the importance of HA in tumor development as well as in anti-cancer drug development. However, existing analyses are limited to a single HA regulator, and the role of HA modification is characterized by highly integrated interactions among many regulators. Therefore, a comprehensive understanding of the regulatory network of HA regulators and their important roles in cancer will facilitate the development of anti-cancer therapeutic strategies.

In this study, we retrospectively investigated genomic data from more than 10,000 tumor samples from TCGA project. We found that the 36 HA regulators had low levels of overall mutation frequency in human cancers, with mutation frequencies varying with cancer context and regulators. Furthermore, we also observed a higher overall mutation frequency of HA regulators in UCEC and SKCM than in other cancers, which was consistent with a previous report [10]. We investigated the correlation between gene expression of the HA regulators and copy number, and found that the expression levels of almost all the regulators were positively correlated with copy number. Compared with normal tissues, the expression of CNV-amplified HA regulators was significantly increased in tumor tissues, while the expression of CNV-deficient HA regulators was significantly decreased, suggesting that the CNV alteration may be one of the main mechanisms leading to expression perturbation of HA regulators.

In the analysis of HA regulator related oncogenic pathways, we observed that the expression of each HA regulator was associated with the activation or inhibition of multiple pathways. Among them, HDAC3 was shown to be associated with the activity of several pathways, including E2F target, G2M checkpoint, MYC target V1, epithelial mesenchymal transition, and PI3K/AKT/mTOR signaling pathway and so on. HDAC3 has been shown to regulate the biological activities of colon cancer cells, including proliferation, differentiation and apoptosis [31, 32]. Moreover, HDAC3 was also confirmed to be overexpressed in triple-negative breast cancer [33], histone deacetylase inhibitors (I-7ab) specifically reduced HDAC3 expression and promoted acetylation of p53 to induce expression of p21, resulting in cell cycle arrest in G1 phase [34]. HDAC3 was reported to be an important player in the development of acute promyelocytic leukemia, and knockdown of HDAC3 can inhibit the PI3k/Akt-mediated signaling pathways and the induction of Caspase activity, thus leading to cell death and apoptosis [35, 36]. HDAC3 has also been found to be involved in the development of many other malignancies, such as melanoma[37], gastric cancer [38], ovarian cancer[39], and children glioma [40]. The important role of HDAC3 in LUAD is still undercharacterized. This study clarified its oncogenic role in LUAD based on bioinformatics and in vitro experiments. In A549, knockdown of HDAC3 could significantly inhibit malignant biological behaviors such as cell proliferation, migration, and invasion, indicating an oncogenic function of HDAC3 in LUAD, which was consistent with it in other cancer types. Epithelial mesenchymal transition is an important biological process for epithelial cell-derived malignant tumor cells to acquire the ability of migration and invasion. In fact, the relationship between HDAC3 and epithelial mesenchymal transition has been reported in malignant tumors and other diseases [41, 42]. For example, targeting HDAC3 can block epithelial mesenchymal transition plasticity in gastric cancer[43], HDAC3 can enhance the migration and invasion properties of fibroblasts by positively affecting the epithelial mesenchymal transition process [44]. In this study, we found that the knockdown of HDAC3 could significantly inhibit the epithelial mesenchymal transition of A549, further revealing an important link between HDAC3 and epithelial mesenchymal transition. Although the oncogenic role of HDAC3 in LUAD has been identified, the specific mechanism remains to be explored. Considering that HDAC3 is a histone deacetylase encoding gene, future exploration of its carcinogenic mechanism can focus on the bridge relationship between the deacetylation modification of downstream regulatory factors by HDAC3 and epithelial mesenchymal transformation.

This study identified genomic alterations and expression perturbations of HA regulators in a variety of cancers, and we are therefore interested in further investigating their clinical translational value. We found a significant cancer-specific relationship between HA regulators and patient outcomes. Of these 33 cancer types, KIRC had the most prognosis-related regulators. Thus, we further explored the significance of the expression of HA regulators in the prognostic stratification of KIRC. We found that KIRC patients could be classified into three clinical clusters with different prognostic characteristics (HAsCluster A, HAsCluster B, and HAsCluster C). This suggested that the transcription profile signature based on HA regulators had prognostic stratification significance. The vast majority of HA regulators were least expressed in HAsCluster B, suggesting that the activity of HA might be relatively low in HAsCluster B. Surprisingly, HAsCluster B had significantly better clinical outcomes, either in terms of OS, DSS, or PFI. This suggested that there might be no one-to-one correspondence between HA and patient outcomes. The development of genomics drives the discovery of tumor-related genes, which will greatly promote the development of anti-cancer drugs. We combined the CellMiner database to explore potential targeted drugs for HA regulators. Using stringent screening criteria, we identified 36 gene-drug pairs involving 13 genes and 27 drugs. These agents was in clinical trials or approved by the FDA, suggesting that they might be potential anti-cancer agents targeting HA regulators, which would provide reference for the development of anti-cancer clinical drugs for follow-up studies. Selumetinib was one of the drugs identified by our CellMiner-based drug sensitivity analysis. Selumetinib was an inhibitor of mitogen-activated protein kinase 1 and 2 (MEK1/2), which were an upstream regulator of the ERK pathway and were frequently activated in a variety of cancers[45,46,47]. Our results revealed correlation of selumetinib and HDAC7, which has been rarely been highlighted in previous studies. This result may provide some direction for future clinical anti-cancer drug development, but still needs further experimental exploration. Taken all together, our results demonstrate the great potential of HA regulators in developing novel therapeutic strategies and in prognostic stratification in specific cancer types.

Overall, this study presents a comprehensive analysis of multiple levels of data on HA regulators in more than 10,000 patients, including somatic mutations, copy number variations, mRNA expression, prognostic value, potential target drugs, etc. This study systematically revealed the molecular features and clinical significance of HA regulators in human cancers, but several key questions still need to be addressed. First, the description of the molecular characteristics of HA regulators in this study is mainly based on the TCGA platform, and there is still a lack of further verification of other platforms and appropriate molecular biology experiments. Although we identified the oncogenic role of HDAC3 in LUAD by cell function experiments, the specific molecular mechanisms remain unelucidated. Secondly, this study only preliminarily identified 36 gene-drug pairs, and a total of 27 drugs may have anti-HA regulator effects, but whether they have anti-cancer activity remains to be investigated. In conclusion, our study for the first time revealed the important role of HA regulators in human cancer, which provides a rich resource for understanding the biology of HA regulators and a new perspective for developing cancer therapeutic strategies based on HA modification.

Data Availability

The datasets generated and analysed during the current study are available in the TCGA GDC repository (https://portal.gdc.cancer.gov), GEO repository (GSE10072 and GSE32863), (https://www.ncbi.nlm.nih.gov/geo/), and so on.

Abbreviations

HA:

Histone acetylation

TGGA:

The Cancer Genome Atlas

ACC:

Adrenocortical Carcinoma

BLCA:

Bladder Urothelial Carcinoma

BRCA:

Breast Cancer

CESC:

Cervical Aquamous Cell Carcinoma And Endocervical Adenocarcinoma

CHOL:

Cholangiocarcinoma

COAD:

Colon Adenocarcinoma

DLBC:

Lymphageal neoplasm diffuse large b-cell lymphoma

ESCA:

Esophageal Caecinoma

GBM:

Glioblastoma Multiforme

HNSC:

Head And Neck Squamous Carcinoma

READ:

Rectum Adenocarcinoma

KICH:

Kidney Chromophobe

KIRC:

Kidney renal clear cell carcinoma

KIRP:

Kidney Renal Clearn Cell Carcinoma

LAML:

Acute Myeloid Leukemia

LGG:

Brain Lower Grade Glioma

LIHC:

Liver Hepatocellular Carcinoma

LUAD:

Lung adenocarcinoma

LUSC:

Lung Squamous cell carcinoma

MESO:

Mesothelioma

OV:

Ovarian serous cystadenocarcinoma

PADD:

Pancereatic Adenocarcinoma

PCPG:

Pheochromocytoma And Paraganglima

PRAD:

Prostate Adenocarcinoma

SARC:

Sarcoma

SKCM:

Skin Cutaneous Melanoma

STAD:

Stomach Adenocarcinoma

TGCT:

Testicular Germ Cell Tumors

THCA:

Thyroid Carcinoma

THYM:

Thymoma

UCEC:

Uterine Corpus Endometrial Carcinoma

UCS:

Uterine Carcinosarcoma

UVM:

Uveal Melanoma

GSEA:

Gene-set enrichment analysis

CNV:

Copy number variation

GSVA:

Gene set variation analysis

SCNA:

Somatic copy number alterations

References

  1. Hanahan D. Hallmarks of Cancer: New Dimensions. CANCER DISCOV. 2022;12(1):31–46.

    CAS  PubMed  Google Scholar 

  2. Dai E, Zhu Z, Wahed S, Qu Z, Storkus WJ, Guo ZS. Epigenetic modulation of antitumor immunity for improved cancer immunotherapy. MOL CANCER. 2021;20(1):171.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. Liu S, Chang W, Jin Y, Feng C, Wu S, He J, Xu T. The function of histone acetylation in cervical cancer development. Biosci Rep 2019, 39(4).

  4. Calcagno DQ, Wisnieski F, Mota E, Maia DSS, Costa DSJ, Leal MF, Gigek CO, Santos LC, Rasmussen LT, Assumpcao PP, et al. Role of histone acetylation in gastric cancer: implications of dietetic compounds and clinical perspectives. EPIGENOMICS-UK. 2019;11(3):349–62.

    CAS  Google Scholar 

  5. Guo P, Chen W, Li H, Li M, Li L. The histone acetylation modifications of breast Cancer and their therapeutic implications. PATHOL ONCOL RES. 2018;24(4):807–13.

    CAS  PubMed  Google Scholar 

  6. Zhang L, Bu L, Hu J, Xu Z, Ruan L, Fang Y, Wang P. HDAC1 knockdown inhibits invasion and induces apoptosis in non-small cell lung cancer cells. BIOL CHEM. 2018;399(6):603–10.

    CAS  PubMed  Google Scholar 

  7. Darvishi N, Rahimi K, Mansouri K, Fathi F, Menbari MN, Mohammadi G, Abdi M. MiR-646 prevents proliferation and progression of human breast cancer cell lines by suppressing HDAC2 expression. Mol Cell Probes. 2020;53:101649.

    CAS  PubMed  Google Scholar 

  8. Cao J, Yan Q. Cancer Epigenetics, Tumor Immunity, and Immunotherapy. Trends Cancer. 2020;6(7):580–92.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Xu Y, Liao W, Luo Q, Yang D, Pan M. Histone Acetylation Regulator-Mediated acetylation patterns define Tumor Malignant Pathways and Tumor Microenvironment in Hepatocellular Carcinoma. FRONT IMMUNOL. 2022;13:761046.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Li Y, Xiao J, Bai J, Tian Y, Qu Y, Chen X, Wang Q, Li X, Zhang Y, Xu J. Molecular characterization and clinical relevance of m(6)a regulators across 33 cancer types. MOL CANCER. 2019;18(1):137.

    PubMed  PubMed Central  Google Scholar 

  11. Liu Z, Zhao Q, Zuo ZX, Yuan SQ, Yu K, Zhang Q, Zhang X, Sheng H, Ju HQ, Cheng H, et al. Systematic analysis of the Aberrances and functional implications of ferroptosis in Cancer. iScience. 2020;23(7):101302.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Landi MT, Dracheva T, Rotunno M, Figueroa JD, Liu H, Dasgupta A, Mann FE, Fukuoka J, Hames M, Bergen AW, et al. Gene expression signature of cigarette smoking and its role in lung adenocarcinoma development and survival. PLoS ONE. 2008;3(2):e1651.

    PubMed  PubMed Central  Google Scholar 

  13. Selamat SA, Chung BS, Girard L, Zhang W, Zhang Y, Campan M, Siegmund KD, Koss MN, Hagen JA, Lam WL, et al. Genome-scale analysis of DNA methylation in lung adenocarcinoma and integration with mRNA expression. GENOME RES. 2012;22(7):1197–211.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Gustavsson EK, Zhang D, Reynolds RH, Garcia-Ruiz S, Ryten M. Ggtranscript: an R package for the visualization and interpretation of transcript isoforms using ggplot2. Bioinformatics. 2022;38(15):3844–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Panneerdoss S, Eedunuri VK, Yadav P, Timilsina S, Rajamanickam S, Viswanadhapalli S, Abdelfattah N, Onyeagucha BC, Cui X, Lai Z, et al. Cross-talk among writers, readers, and erasers of m(6)a regulates cancer growth and progression. SCI ADV. 2018;4(10):r8263.

    Google Scholar 

  16. Barabasi AL, Gulbahce N, Loscalzo J. Network medicine: a network-based approach to human disease. NAT REV GENET. 2011;12(1):56–68.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Kim JH, Yoon CY, Jun Y, Lee BB, Lee JE, Ha SD, Woo H, Choi A, Lee S, Jeong W, et al. NuA3 HAT antagonizes the Rpd3S and Rpd3L HDACs to optimize mRNA and lncRNA expression dynamics. NUCLEIC ACIDS RES. 2020;48(19):10753–67.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Werbeck ND, Shukla VK, Kunze M, Yalinca H, Pritchard RB, Siemons L, Mondal S, Greenwood S, Kirkpatrick J, Marson CM, et al. A distal regulatory region of a class I human histone deacetylase. NAT COMMUN. 2020;11(1):3841.

    PubMed  PubMed Central  Google Scholar 

  19. Adhikari N, Amin SA, Trivedi P, Jha T, Ghosh B. HDAC3 is a potential validated target for cancer: an overview on the benzamide-based selective HDAC3 inhibitors through comparative SAR/QSAR/QAAR approaches. EUR J MED CHEM. 2018;157:1127–42.

    CAS  PubMed  Google Scholar 

  20. Wang L, Zou X, Berger AD, Twiss C, Peng Y, Li Y, Chiu J, Guo H, Satagopan J, Wilton A, et al. Increased expression of histone deacetylaces (HDACs) and inhibition of prostate cancer growth and invasion by HDAC inhibitor SAHA. AM J TRANSL RES. 2009;1(1):62–71.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Urdinguio RG, Lopez V, Bayon GF, Diaz DLGR, Sierra MI, Garcia-Torano E, Perez RF, Garcia MG, Carella A, Pruneda PC, et al. Chromatin regulation by histone H4 acetylation at lysine 16 during cell death and differentiation in the myeloid compartment. NUCLEIC ACIDS RES. 2019;47(10):5016–37.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Vaseghi H, Liu J, Qian L. Molecular barriers to direct cardiac reprogramming. PROTEIN CELL. 2017;8(10):724–34.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Zhou HX, Pang X. Electrostatic interactions in protein structure, folding, binding, and condensation. CHEM REV. 2018;118(4):1691–741.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Munster PN, Marchion D, Thomas S, Egorin M, Minton S, Springett G, Lee JH, Simon G, Chiappori A, Sullivan D, et al. Phase I trial of vorinostat and doxorubicin in solid tumours: histone deacetylase 2 expression as a predictive marker. Br J Cancer. 2009;101(7):1044–50.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Hess-Stumpp H. Histone deacetylase inhibitors and cancer: from cell biology to the clinic. EUR J CELL BIOL. 2005;84(2–3):109–21.

    CAS  PubMed  Google Scholar 

  26. Bianco-Miotto T, Chiam K, Buchanan G, Jindal S, Day TK, Thomas M, Pickering MA, O’Loughlin MA, Ryan NK, Raymond WA, et al. Global levels of specific histone modifications and an epigenetic gene signature predict prostate cancer progression and development. Cancer Epidemiol Biomarkers Prev. 2010;19(10):2611–22.

    CAS  PubMed  Google Scholar 

  27. Elsheikh SE, Green AR, Rakha EA, Powe DG, Ahmed RA, Collins HM, Soria D, Garibaldi JM, Paish CE, Ammar AA, et al. Global histone modifications in breast cancer correlate with tumor phenotypes, prognostic factors, and patient outcome. CANCER RES. 2009;69(9):3802–9.

    CAS  PubMed  Google Scholar 

  28. Yao Y, Liu Q, Adrianto I, Wu X, Glassbrook J, Khalasawi N, Yin C, Yi Q, Dong Z, Geissmann F, et al. Histone deacetylase 3 controls lung alveolar macrophage development and homeostasis. NAT COMMUN. 2020;11(1):3822.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Topper MJ, Vaz M, Marrone KA, Brahmer JR, Baylin SB. The emerging role of epigenetic therapeutics in immuno-oncology. NAT REV CLIN ONCOL. 2020;17(2):75–90.

    PubMed  Google Scholar 

  30. Li X, Yu W, Qian X, Xia Y, Zheng Y, Lee JH, Li W, Lyu J, Rao G, Zhang X, et al. Nucleus-translocated ACSS2 promotes gene transcription for lysosomal Biogenesis and Autophagy. MOL CELL. 2017;66(5):684–97.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Wilson AJ, Byun DS, Popova N, Murray LB, L’Italien K, Sowa Y, Arango D, Velcich A, Augenlicht LH, Mariadason JM. Histone deacetylase 3 (HDAC3) and other class I HDACs regulate colon cell maturation and p21 expression and are deregulated in human colon cancer. J BIOL CHEM. 2006;281(19):13548–58.

    CAS  PubMed  Google Scholar 

  32. Spurling CC, Godman CA, Noonan EJ, Rasmussen TP, Rosenberg DW, Giardina C. HDAC3 overexpression and colon cancer cell proliferation and differentiation. Mol Carcinog. 2008;47(2):137–47.

    CAS  PubMed  Google Scholar 

  33. Hanigan TW, Aboukhatwa SM, Taha TY, Frasor J, Petukhov PA. Divergent JNK phosphorylation of HDAC3 in Triple-Negative breast Cancer cells determines HDAC inhibitor binding and selectivity. CELL CHEM BIOL. 2017;24(11):1356–67.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Yang M, Dang X, Tan Y, Wang M, Li X, Li G. I-7ab inhibited the growth of TNBC cells via targeting HDAC3 and promoting the acetylation of p53. BIOMED PHARMACOTHER. 2018;99:220–6.

    CAS  PubMed  Google Scholar 

  35. Zhang L, Hong Z, Zhang RR, Sun XZ, Yuan YF, Hu J, Wang X. Bakkenolide A inhibits leukemia by regulation of HDAC3 and PI3K/Akt-related signaling pathways. BIOMED PHARMACOTHER. 2016;83:958–66.

    CAS  PubMed  Google Scholar 

  36. Mehdipour P, Santoro F, Botrugno OA, Romanenghi M, Pagliuca C, Matthews GM, Johnstone RW, Minucci S. HDAC3 activity is required for initiation of leukemogenesis in acute promyelocytic leukemia. Leukemia. 2017;31(4):995–7.

    CAS  PubMed  Google Scholar 

  37. Shan X, Fu YS, Aziz F, Wang XQ, Yan Q, Liu JW. Ginsenoside Rg3 inhibits melanoma cell proliferation through down-regulation of histone deacetylase 3 (HDAC3) and increase of p53 acetylation. PLoS ONE. 2014;9(12):e115401.

    PubMed  PubMed Central  Google Scholar 

  38. Xu G, Zhu H, Zhang M, Xu J. Histone deacetylase 3 is associated with gastric cancer cell growth via the mir-454-mediated targeting of CHD5. INT J MOL MED. 2018;41(1):155–63.

    CAS  PubMed  Google Scholar 

  39. Hayashi A, Horiuchi A, Kikuchi N, Hayashi T, Fuseya C, Suzuki A, Konishi I, Shiozawa T. Type-specific roles of histone deacetylase (HDAC) overexpression in ovarian carcinoma: HDAC1 enhances cell proliferation and HDAC3 stimulates cell migration with downregulation of E-cadherin. INT J CANCER. 2010;127(6):1332–46.

    CAS  PubMed  Google Scholar 

  40. Zhu J, Wan H, Xue C, Jiang T, Qian C, Zhang Y. Histone deacetylase 3 implicated in the pathogenesis of children glioma by promoting glioma cell proliferation and migration. BRAIN RES. 2013;1520:15–22.

    CAS  PubMed  Google Scholar 

  41. Mani SA, Barton MC. HDAC3 at the fulcrum of an epithelial-mesenchymal balance. MOL CELL. 2011;43(5):697–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Yang MH, Wu MZ, Chiou SH, Chen PM, Chang SY, Liu CJ, Teng SC, Wu KJ. Direct regulation of TWIST by HIF-1alpha promotes metastasis. NAT CELL BIOL. 2008;10(3):295–305.

    CAS  PubMed  Google Scholar 

  43. Wu SM, Jan YJ, Tsai SC, Pan HC, Shen CC, Yang CN, Lee SH, Liu SH, Shen LW, Chiu CS et al. Targeting histone deacetylase-3 blocked epithelial-mesenchymal plasticity and metastatic dissemination in gastric cancer. CELL BIOL TOXICOL 2022.

  44. Jeong SH, Son ES, Lee YE, Kyung SY, Park JW, Kim SH. Histone deacetylase 3 promotes alveolar epithelial-mesenchymal transition and fibroblast migration under hypoxic conditions. EXP MOL MED. 2022;54(7):922–31.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. Fangusaro J, Onar-Thomas A, Young PT, Wu S, Ligon AH, Lindeman N, Banerjee A, Packer RJ, Kilburn LB, Goldman S, et al. Selumetinib in paediatric patients with BRAF-aberrant or neurofibromatosis type 1-associated recurrent, refractory, or progressive low-grade glioma: a multicentre, phase 2 trial. LANCET ONCOL. 2019;20(7):1011–22.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. Valle JW, Lamarca A, Goyal L, Barriuso J, Zhu AX. New Horizons for Precision Medicine in biliary tract cancers. CANCER DISCOV. 2017;7(9):943–62.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Carter CA, Rajan A, Keen C, Szabo E, Khozin S, Thomas A, Brzezniak C, Guha U, Doyle LA, Steinberg SM, et al. Selumetinib with and without erlotinib in KRAS mutant and KRAS wild-type advanced nonsmall-cell lung cancer. ANN ONCOL. 2016;27(4):693–9.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This work mainly benefited from open databases. We are grateful for the efforts made by the resources and staff to expand and improve the databases.

Funding

This work was supported by grants from the National Natural Science Foundation of China (No81770095) and Science Fund for Creative Research Groups of the Natural Science Foundation of Hubei Province (No. 2020CFA027).

Author information

Authors and Affiliations

Authors

Contributions

CKS, XFL, KL and WCL: made significant contributions to the study design, data collection and analysis, and article writing. CKS, and KL: literature search, data collection, data analysis, data interpretation, Database construction,etc. SZP, and ZLL: participated in drafting and revising the article critically. CKS, DHL and KL: conceptualisation, data curation, formal analysis, methodology. QG and NL: funding acquisition, investigation, project administration, resources. QG and other co-authors read and approved the final manuscript.

Corresponding authors

Correspondence to Ning Li or Qing Geng.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Patient consent for publication

Not applicable.

Competing interests

The authors declared that they had no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Song, C., Liu, X., Lin, W. et al. Systematic analysis of histone acetylation regulators across human cancers. BMC Cancer 23, 733 (2023). https://doi.org/10.1186/s12885-023-11220-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-11220-7

Keywords