Histone regulator KAT2A acts as a potential biomarker related to tumor microenvironment and prognosis of diffuse large B cell lymphoma
BMC Cancer volume 23, Article number: 934 (2023)
Recent studies have indicated that epigenetic alterations contribute significantly to lymphoma pathogenesis. A type of epigenetic regulation known as histone acetylation plays a crucial role in transcriptional regulation in eukaryotic cells. Specifically, a significant effect of histone acetylation modifications on the abnormal progression and microenvironment of diffuse large B-cell lymphoma (DLBCL) has been observed.
To provide insight into the significance of histone acetylation-related genes, we developed a HAscore model for analyzing histone acetylation patterns in DLBCL samples. Furthermore, KAT2A, a regulator of histone acetylation, was knocked down in DLBCL cell lines to investigate its role in proliferation, cell cycle, and apoptosis.
The HAscore model has been demonstrated to provide insight into the significance of these patterns, showing that patients with a low HAscore have distinct tumor immune microenvironments and poorer prognoses. Besides, KAT2A was identified as a potential biomarker related to immune infiltration and malignant pathways in DLBCL.
According to these findings, it is evident that the histone acetylation pattern score model is helpful in describing the immune status of DLBCL and that KAT2A may be used as a biomarker for its treatment.
Diffuse large B-cell lymphoma (DLBCL) is a highly aggressive form of non-Hodgkin’s lymphoma (NHL), originating from B-lineage lymphocytes, accounting for 30–40% of cases [1,2,3]. Moreover, this lymphoid neoplasm has highly variable gene expression profiles and genetic alterations [4, 5]. The most common up-front treatment is R-CHOP . However, due to its heterogeneity, 60% of patients are curable with combination therapy and the remainders still succumb to the disease . Unfortunately, those who develop the disease refractory to up-front treatment or relapse after remission have a poor prognosis [8,9,10]. To understand the underlying mechanisms of DLBCL progression, more studies are necessary. As a result, stratifying DLBCL patients and developing predictive models can provide more precise molecular subtyping and, therefore, more customized treatment.
Epigenetic alterations are directly linked to lymphoma pathogenesis. As one of these post-translational modifications, histone acetylation has been extensively studied . Post-translational modifications of histones regulate transcription and DNA repair and are linked to the stable maintenance of repressive chromatin . The N-terminal tail of histones contains lysine residues that are tightly regulated by acetylation regulators. Acetylation reactions are catalyzed by histone acetyltransferases (HATs), while deacetylation reactions are catalyzed by histone deacetylases (HDACs) . Besides, other regulators recognize modified histones and play a role in acetylation. The three types of regulators are regarded as “writers,“ “readers,“ and “erasers.“ Increasing evidence suggests that the tumor microenvironment (TME) plays a critical role in malignancy development and maintenance. The process is accomplished through sustained proliferation and immune evasion [14, 15]. As a result of several signal molecules being activated and inhibited, histone acetylation has recently been linked to cancer and TME. In various types of human malignancies, aberrant expression of HDACs has been reported [16, 17]. It has also been reported that many HDAC inhibitors are effective against several hematologic and solid malignancies [18, 19]. The HAT paralogs p300 and CBP are involved in many vital cellular processes and have critical roles in several pathological conditions, including cancer [20,21,22]. Current studies indicated that somatic mutations affecting CREBBP and EP300 are a hallmark of DLBCL . Multiple studies have also shown that histone deacetylase inhibitors (HDACi) can improve the abilities of the immune system to eradicate tumor cells by changing the TME through various mechanisms [24, 25].
Collectively, histone acetylation is critical in regulating TME. Targeting histone acetylation modulators can disrupt the resistance to cancer immunotherapy. Recent research has concentrated on individual histone acetylation modulators and their impact on cancer treatment and prognosis. We retrospectively collected transcription information from public databases to better understand how histone acetylation regulators impact the immune system. We investigated histone acetylation regulatory variables and infiltration of immune cells, as well as the value of HAscore in targeted DLBCL immunotherapy.
Materials and methods
Dataset acquisition and clinical samples
The Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA) databases provided gene expression data and comprehensive clinical annotations. Transcriptome data was derived from fragments per kilobase million and converted into transcripts per million (TPM) and then transformed into log2 (TPM + 1) values for further analysis. The UCSC Xena database also provided genomic mutations (including somatic mutations and copy number variations (CNV). This study focused on the DLBCL cohorts (GSE10846 and GSE31312) as well as the TCGA-DLBC. The R package limma “normalizeBetweenArrays” package was used for normalization after gene symbol conversion (R version: 4.1.2; Bioconductor version: 3.13). A method in “sva” package called “ComBat” was utilized to adjust for batch effects caused by non-biotechnical bias [26, 27].
Human DLBCL cell lines OCI-LY1, OCI-LY3, OCI-LY8, OCI-LY10, and U2932 were cultured in IMDM (Gibco, MD, USA), supplemented with 10% FBS (Gibco). The cells were maintained under optimal conditions, with a temperature of 37 °C and a humid atmosphere containing 5% CO2. Peripheral blood mononuclear cells and serum were isolated from healthy volunteers.
Lentiviral generation and cell transfection
The stable knockdown of KAT2A was encoded by cloning the shRNAs into lentiviral vectors (Beijing Syngentech Co., Ltd., China). The manufacturer’s instructions were followed for lentivirus infection. To select the stably transfected cells, the medium was supplemented with puromycin (2.0 g/ml, Sigma-Aldrich, USA). A 72-hour timeframe was used for the collection and analysis of cells.
Total RNA was extracted using RNAiso Plus (TaKaRa, Dalian, China) following the manufacturer’s protocol. Reverse transcription reactions were undertaken using reverse transcription reagents (Vazyme). Expression levels of specific genes were then measured using qRT-PCR on LightCycler 480II system (Roche, Basel, Switzerland). Normalized results were determined based on GAPDH expression.
The KAT2A primer sequence was as follows:
KAT2A: forward 5-CCCGCTACGAAACCACTCAT-3, reverse 5-GCATGGACAGGAATTTGGGGA-3.
Cell proliferation assay
The cells were added to 96-well plates at a density of 1 × 104 per well. The proliferation of cells was determined after 24, 48, and 72 h of exposure to CCK-8 (Dojindo Laboratories, Kumamoto, Japan). After incubation for four hours, the OD was measured at 450 nm.
The Western blot analysis was performed as described previously . The primary antibody used in this assay was KAT2A (66575-1-Ig, Proteintech). GAPDH was used as an internal reference.
Cell cycle and apoptosis were determined by Navios flow cytometer. Cells were stained for 15 min and monitored for cell cycle. Annexin V-PE/7-aminoactinomycin (7AAD) apoptosis detection assays (BD Biosciences) were used to detect apoptosis.
Consensus clustering analysis
38 recognized histone acetylation regulators were isolated from the published literature [28,29,30]. An unsupervised clustering algorithm of 884 DLBCL samples was performed based on their expression levels for 38 histone acetylation regulators. The patients were clustered using the R package “consensusClusterplus”, and the classification was performed 1,000 times to ensure stability . Principal Component Analysis (PCA) was employed to illustrate the clustering conditions of given samples with tumors by dimension reduction.
Gene set variation analysis (GSVA)
We used the “GSVA” R package to analyze differences in histone acetylation modification patterns between biological processes. The gene sets “c2.cp.kegg.v7.5.1.symbols” were accessed from MSigDB. Genes related to histone acetylation modifications were functionally annotated with a threshold of FDR < 0.05.
Identification of DEGs in histone acetylation modification characteristics
Patients were divided into three groups by their expression of histone acetylation regulators. An analysis of differentially expressed genes (DEGs) exhibiting various characteristics among histone acetylation patterns was conducted using the limma R package . Accordingly, the DEGs were selected based on |logFC|> 1 and adjusted P values < 0.001. Kyoto Encyclopedia of Genes and Genomes (KEGG) [33,34,35] pathway analysis and Gene Ontology (GO) biological were processed using the R package “clusterProfiler” and “org.Hs.eg.db”.
Estimation of tumor environment cell infiltration
Adapted from Charoentong’s study , the gene set was used to define immune cell types. A fraction of the immune cells infiltrating each sample was determined based on their relative abundance. Based on the single-sample gene set enrichment analysis (ssGSEA) algorithm, ESTIMATE was calculated by the R package “ESTIMATE” based on immune cell and stromal cell-specific gene expression levels to calculate a score reflecting the level of immune cell and stromal cell infiltration. Boxplots of ESTIMATE score, immune score, and stromal score for different risk groups were demonstrated using the R package “ggpubr”.
Evaluation and generation of HAscore
There has been an effort to develop a method for quantifying the pattern of histone acetylation modifications by developing a scoring approach (HAscore). The overlapping DEGs from different HA clusters were first accessed to divide patients into several clusters via unsupervised clustering. Following a univariate Cox regression model, prognostic-related DEGs were selected to construct HA gene signatures. Following determining the prognostic value of gene signature scores, a method similar to the Genome Grading Index  was applied to define the HAscore for each patient: HAscore = ∑(PC1i + PC2i), where i indicates the expression value of each histone acetylation regulator.
Patients’ survival was analyzed by the Kaplan-Meier analysis, and overall survival (OS) between subgroups by the log-rank test was determined. Using the “surv-cutpoint” function in the “survminer” R software package, the optimal cut-off point was determined, resulting in the classification of patients into high- and low-HAscore subgroups, as well as high- and low-KAT2A expression subgroups. The difference between three or more groups was examined using one-way ANOVA and Kruskal-Wallis tests . p values were analyzed using two-sided statistical tests in all statistical analyses, and p < 0.05 was considered statistically significant. At least three independent experiments were conducted, and the mean and standard deviation (SD) were used in describing the experimental data. R and GraphPad Prism (version 8.0) were utilized to conduct all the statistical analyses.
Genetic variants of DLBCL histone acetylation regulators
A systematic review of histone acetylation publications identified 43 histone acetylation regulatory genes, including 10 “writers”, 18 “erasers”, and 15 “readers” (Table S1). The rate of somatic mutations of 43 histone acetylation regulators was assessed to determine the genetic alterations in DLBCL. As evidenced by the TCGA database, 12 of the analyzed samples had mutations in histone acetylation regulators, with a frequency of 32.43% (Fig. 1A). The study results showed that CREBBP had the highest mutation frequency, followed by KAT6A, EP300, HDAC8, HDAC9, BRD4, BRD2, PBRM1, SMARCA2, BAZ2B, and SMARCA4. The remaining genes had no mutations in the DLBCL samples. In terms of the mutation status of CREBBP, DLBCL samples were grouped into two categories: mutation and wild group. Several genes were found to be overexpressed in the mutation group (Fig. 1B). Regarding CNVs, HDAC9, YEATS4, HDAC7, HDAC11, and KAT2B exhibited high amplification rates, whereas KAT6A exhibited copy number loss (Fig. 1C). Chromosome-wide copy number alterations of HA regulators have been identified (Fig. 1D). Subsequently, we investigated the expression levels of histone acetylation regulators in both normal and DLBCL tissues in “writers”, “erasers”, and “readers”, respectively, to investigate the differences in the expression of these regulators in the two groups of sample (Fig. 1E-G). Most genes were upregulated in DLBCL samples compared to normal samples, except for HDAC5, DPF3, and SMARCA2. There is a high degree of heterogeneity in histone acetylation regulators, suggesting that aberrant expression patterns of HA regulators may contribute to carcinogenesis and progression of the disease.
Construction of histone acetylation patterns
A univariate Cox regression model demonstrated the predictive significance of histone acetylation regulators (Fig. 2A and Table S2). The results showed that multiple regulators were risk factors for DLBCL, and some regulators were protective factors. An analysis of spearman correlations was conducted in order to describe the interactions among the regulator connections (Fig. 2B). According to the results, most of the regulators were positively correlated. A significant correlation was found between genes belonging to the same biological group (writers, erasers, and readers), and a positive correlation was also found between genes with opposite biological functions. The results suggest that there is extensive crosstalk among histone acetylation regulators, which collectively regulate histone acetylation modifications and influence DLBCL development. HA regulators were analyzed with unsupervised clustering to identify three patterns of histone acetylation. 884 DLBCL patient samples were grouped according to their gene expression, with 163 cases belonging to cluster A, 455 cases belonging to cluster B, and 266 cases belonging to cluster C (Fig. S1A-D). Significantly, PCA analysis revealed differences between the three different histone acetylation patterns (Fig. 2C). Prognostic analysis of the three subtypes of histone acetylation modifications revealed that HAcluster A had the worst prognosis (Fig. 2D). Additionally, the correlation between clinical factors and gene expression was analyzed. The results presented in Fig. 2E illustrated the correlation among three distinct histone acetylation patterns. Compared to the other two groups, HAcluster A showed a significantly higher incidence of activated B-cell–like (ABC) DLBCL. There is a well-established association between ABC DLBCL and a poor prognosis , indicating decreased survival observed among HAcluster A patients.
Immune infiltration and biological process associated with histone acetylation patterns
In accordance with the KEGG gene set, GSVA enrichment analysis examined the differences in biological behavior among the three modification patterns (Fig. 3A-B, Fig. S1E). The results of GSVA analysis revealed that HAcluster A was enriched in several metabolically relevant pathways. Additionally, HAcluster A was enriched in pathways related to tumors, including the DRUG-METABOLISM/JAK-STAT/MAPK/NOTCH/PARP signaling pathway and apoptosis, in comparison to HAcluster B and HAcluster C (Fig. 3C). Subsequently, TME cell infiltration analysis showed substantial differences in histone acetylation patterns among the three groups, with HAcluster A exhibiting higher levels of activated dendritic cells and MDSCs and lower levels of nature killer cells compared to HAcluster B and HAcluster C (Fig. 3D). Based on RNA expression levels of DLBCL, the ESTIMATE algorithm was employed to calculate the stromal, immune, and estimated scores (Fig. 3E-F). The boxplots showed that HAcluster A had significantly lower StromalScore and ESTIMATEScore and had significantly higher tumor purity than the other two groups. This is consistent with the previous analysis. According to the previous research, we assessed the immunosuppressive, immune cytolytic, and antigen-processing effects (Fig. 3G) [37, 39]. The results showed that HAcluster A had the highest immunosuppression, which indicated that the tumor microenvironment of HAcluster A was significantly different.
Construction of HA gene signature
For a better understanding of the differences among these three histone acetylation patterns, we identified 929 DEGs (Fig. S2A) significantly associated with patient survival in the previous cohort. This was followed by an unsupervised cluster analysis that classified patients into different genotypes (Fig. S2B-E). In the analysis of GO and KEGG enrichment of these DEGs, the main functions were found to be RNA splicing, histone modification, DNA replication, PI3K-AKT, MAPK, and Ras signaling (Fig. 4A-B). A model-based clustering analysis led to the identification of three different modification patterns, including 179 cases in geneCluster A, 450 cases in geneCluster B, and 255 cases in geneCluster C, respectively. In addition, a prognostic crossover could be observed for gene cluster A, which had a good prognosis (Fig. 4C). This suggests that these three different histone acetylation modification patterns are indeed present in DLBCL. As with previous HA modification patterns, significant differences were observed within HA geneClusters (Fig. 4D). Additionally, DEG expression varied significantly among the three HA geneClusters (Fig. 4E). GeneCluster A demonstrated the highest level of DEGs expression, followed by geneCluster B and C, which was also in line with patient survival.
HAscore-related clinical features and TME immune activity
A set of models of histone acetylation modification was developed based on individual differences and complexity of each DLBCL, called the HAscore, which can be used to quantify individual DLBCLs. In Fig. S3, the alluvial chart shows the change in the patient’s individual attributes. In HAcluster A, the majority of samples belonged to geneCluster A, which all had low HAscores. Also of note is that most of the patients in geneCluster B and C were in HAcluster B and C with high HA scores. Using the R package “survminer”, we calculated a cut-off value of -1.18397065. Patients were successfully classified into high and low groups according to their HAscores (Fig. 5A). Additionally, HAscore results indicated that a greater percentage of patients died in the low HAscore group, while surviving patients had significantly higher levels of HAscore than those who died (Fig. 5B-C). Based on these findings, high HAscore patients had a statistically significant survival advantage. Furthermore, Kruskal-Wallis analysis revealed significant differences in HAscores between HAclusters. A significantly lower mean HAscore was observed in HAcluster A and geneCluster A with poor prognoses (Fig. 5D-E). As a result, HAscore levels were strongly correlated with the survival benefits of patients.
For a more comprehensive understanding of the correlation between HAscores and the immune status of patients, we quantified different immune cell subsets and infiltration scores for each patient. Further analysis of the relationship between HAscore and TME revealed a positive correlation between immune cells and high HAscore levels. Based on the CIBERSORT algorithm, we confirmed the association of tumor-initiating cells (TICs) with HAscore (Fig. 5F-G). Furthermore, from the ESTIMATE algorithm, the low HAscore group had low expression in stromalscore, immunescore, and ESTIMATEscore (Fig. 5H). Consequently, the low HAscore group had a higher tumor purity, which is strongly correlated with poor prognosis (Fig. 5I). HAscore groups differed in terms of immune infiltration in this study. This suggests a regulatory role for histone acetylation regulators in the immune microenvironment of patients with DLBCL.
KAT2A was expressed in high levels in DLBCL
In light of the fact that KAT2A is closely related to several histone regulators, HAscore is able to predict the outcomes of DLBCL patients. In addition, we examined whether KAT2A might serve as a potential biomarker in the treatment of DLBCL. Based on the public database, KAT2A was overexpressed in DLBCL patients with poor prognoses (Fig. 6A-B). Similarly, KAT2A was found to be positive in non-Hodgkin’s lymphoma in HPA (Fig. 6C). Furthermore, we observed that KAT2A was mainly overexpressed in many DLBCL cell lines compared to CD19 + B cells in both protein and RNA levels (Fig. 6D-E). KAT2A was found to be associated with T cell activation as well as metabolic processes according to the GO analysis of DEGs (Fig. 6F). Moreover, we analyzed the relationship between KAT2A and TME and found a negative correlation with immune cells. It was found that there is a strong correlation between KAT2A and adaptive immune responses (such as activated CD4 T cells, Type 1 T helper cells, Type 17 T helper cells, Gamma delta T cells) and innate immune responses, such as plasmacytoid dendritic cells and natural killer T cells (Fig. 6G). In particular, low expression of KAT2A is closely associated with cells that perform antitumor functions.
Role of KAT2A in the progression of DLBCL
To elucidate the clinical significance of KAT2A in DLBCL, we knocked down KAT2A in two DLBCL cell lines (OCI-LY1 and U2932), resulting in stable KAT2A KD cell lines (Fig. 7A-B), of which shKAT2A#3 exhibited the highest efficacy. Further function investigations were performed based on shKAT2A#3 cell lines. By staining the cells with PI, we monitored the cell cycle. Compared with shControl cells, there was a significant increase in G2/M phase DLBCL cells after KAT2A was knocked down (Fig. 7C-D). We investigated the effect of KAT2A knockdown on apoptosis (Fig. 7E-F) and found that it significantly increased apoptosis rates. In addition, the stable KAT2A knockdown cells showed growth suppression as opposed to shControl cells (Fig. 7G). Therefore, KAT2A knockdown significantly inhibits the proliferation of the G2/M phase in DLBCL cell lines. As indicated by the findings of the study, KAT2A appeared to be capable of promoting DLBCL progression through its effects on proliferation, apoptosis, and cell cycle progression during the G2/M phase.
Histone acetylation is an epigenetic modification necessary for cancer biology. KAT7, for example, inhibits tumor cell proliferation and invasion by acetylating H4K5 at promoters of FOXO1 and FOXO3a genes , and HDAC5 loss increased H3K27-ac acetylation and circumvented oncogene-related cell cycle repression . Considerable research has been conducted on individual histone acetylation regulators and their role in cancer; however, a comprehensive assessment of these regulators and their interactions remains lacking. Our research aims to develop a more complete understanding of the different histone acetylation patterns and their associated biological properties.
Based on a preliminary analysis of the expression of histone acetylation regulator genes, most genes were significantly different between normal and tumor samples with an association with prognosis. According to unsupervised clustering of 36 modulators, patients were categorized into three histone acetylation phenotypes. Survival analysis revealed three distinct patterns of histone acetylation modification associated with prognosis. To further characterize these different histone acetylation phenotypes, we determined DEGs among them. The HAscore model was developed based on these genes to assess histone acetylation phenotypes in individual patients quantitatively. All patients in HAcluster A have a low HAscore, associated with the poorest survival outcome. Moreover, we found a close relationship between most acetylation-related genes. One of these risk factors, KAT2A, has been found to be closely related to histone acetylation regulators.
In order to shed light on the mechanisms underlying the different prognoses of patients with different phenotypes, the biological characteristics of each pattern were examined. The results of GSVA analysis confirmed that HAcluster A showed enrichment in several metabolically relevant pathways. Tumorigenesis and tumor progression are supported by abnormal metabolic activity, which allows cells to obtain essential nutrients from the environment . Moreover, we found that PARP, JAK-STAT, MAPK, NOTCH, and apoptosis were significantly activated in HAcluster A. Enrichment of these malignancy signaling pathways may be indicative of malignant progression, which may result in a poor prognosis [43,44,45].
A growing number of studies have suggested that TME components may contribute to cancer development [46,47,48]. An investigation into the relationship between histone acetylation modifications and TME cell infiltration was undertaken to to understand the antitumor immune response to DLBCL better. HAcluster A, with the lowest OS, exhibited a significantly higher level of MDSCs, which play a crucial role in immunosuppression . MDSC is characterized by its ability to suppress immune cell function, including inhibition of T cell proliferation . Each pattern is characterized by a different degree of TME infiltration, with immunosuppression characteristic of HAcluster A group. By analyzing histone acetylation patterns, we evaluated the potential therapeutic effects of HAscore based on differences in signaling pathways and tumor microenvironments between the different patterns. In the high-score group, higher levels of immune-activating cells were observed, including activated CD4 T cells and CD56 bright natural killer cells. CD4 + T cells have antitumor activity through the production of effector cytokines that activate CD8 + T cells . By modulating DC and T cell responses, CD56 bright natural killer cells can positively influence the anticancer response . This indicated that the HAscore contributed to further defining the microenvironment and thereby guiding more effective treatment.
Furthermore, our research revealed that KAT2A correlates with DLBCL prognosis and functions independently. An enzyme in the HAT family, lysine acetyltransferase 2 A (KAT2A), is involved in transcriptional activation through histone acetylation, histone succinylation, and recruitment of transcriptional coactivators . Several studies have demonstrated that KAT2A functions as an epigenetic oncogene in several cancers [54, 55]. In previous studies, KAT2A has been shown to be a viable target for reducing the growth of acute myeloid leukemia by significantly promoting myeloid differentiation and apoptosis . However, little was investigated between KAT2A and DLBCL. In this study, we utilized the bioinformatic methods to find that histone regulator KAT2A was a risk factor in DLBCL. The stable KAT2A KD cell lines were constructed for further study. Moreover, our research found that KAT2A deficiency could prohibit cell proliferation, promote cell apoptosis, and arrest cells in the G2/M phase. These findings suggest that inhibiting KAT2A may inhibit tumor growth.
In conclusion, our findings suggest that genetic signatures derived from histone acetylation regulators can be used to provide personalized survival assessments for patients newly diagnosed with DLBCL. It might provide clinicians with valuable information regarding treatment decisions, follow-up, and prognosis for their patients. There are some limitations to the study, despite its strengths. All the data we used were from public databases. It would be beneficial to collect more prospective real-world data to confirm its clinical utility. On the other hand, the specific mechanism of the effect of KAT2A in DLBCL still needs to be further explored.
The datasets presented in this study can be found in online repositories. RNA sequencing data and CNV data were downloaded from the TCGA database (https://www.cancer.gov/) through the UCSC XENA website (https://xena.ucsc.edu/). The RNA data of normal control was downloaded from the Genotype-Tissue Expression (https://www.gtexportal.org/home/). And GSE10846 and GSE31312 were downloaded from the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/).
diffuse large B-cell lymphoma
differentially expressed genes
histone deacetylases inhibition
Gene Expression Omnibus
the Cancer Genome Atlas
transcripts per million
Kyoto Encyclopedia of Genes and Genomes
Gene set variation analysis
principal component analysis
single-sample gene set enrichment analysis
copy number variants
- ABC DLBCL:
activated B-cell–like DLBCL
lysine acetyltransferase 2 A
Zhang Y, Tan H, Daniels JD, Zandkarimi F, Liu H, Brown LM, et al. Imidazole Ketone Erastin induces ferroptosis and slows Tumor Growth in a mouse lymphoma model. Cell Chem Biol. 2019;26(5):623–33e9.
Hartert KT, Wenzl K, Krull JE, Manske M, Sarangi V, Asmann Y, et al. Targeting of inflammatory pathways with R2CHOP in high-risk DLBCL. Leukemia. 2021;35(2):522–33.
Melchardt T, Egle A, Greil R. How I treat diffuse large B-cell lymphoma. ESMO Open. 2023;8(1):100750.
Miao Y, Medeiros LJ, Li Y, Li J, Young KH. Genetic alterations and their clinical implications in DLBCL. Nat Reviews Clin Oncol. 2019;16(10):634–52.
Yang J, Li Y, Zhang Y, Fang X, Chen N, Zhou X, et al. Sirt6 promotes tumorigenesis and drug resistance of diffuse large B-cell lymphoma by mediating PI3K/Akt signaling. J Experimental Clin cancer Research: CR. 2020;39(1):142.
Liu Y, Barta SK. Diffuse large B-cell lymphoma: 2019 update on diagnosis, risk stratification, and treatment. Am J Hematol. 2019;94(5):604–16.
Chapuy B, Stewart C, Dunford AJ, Kim J, Kamburov A, Redd RA, et al. Molecular subtypes of diffuse large B cell lymphoma are associated with distinct pathogenic mechanisms and outcomes. Nat Med. 2018;24(5):679–90.
Sehn LH, Herrera AF, Flowers CR, Kamdar MK, McMillan A, Hertzberg M, et al. Polatuzumab Vedotin in relapsed or refractory diffuse large B-Cell lymphoma. J Clin Oncology: Official J Am Soc Clin Oncol. 2020;38(2):155–65.
Schuster SJ, Bishop MR, Tam CS, Waller EK, Borchmann P, McGuirk JP, et al. Tisagenlecleucel in adult relapsed or refractory diffuse large B-Cell lymphoma. N Engl J Med. 2019;380(1):45–56.
Zhou X, Chen N, Xu H, Zhou X, Wang J, Fang X, et al. Regulation of Hippo-YAP signaling by insulin-like growth factor-1 receptor in the tumorigenesis of diffuse large B-cell lymphoma. J Hematol Oncol. 2020;13(1):77.
Simithy J, Sidoli S, Yuan ZF, Coradin M, Bhanu NV, Marchione DM, et al. Characterization of histone acylations links chromatin modifications with metabolism. Nat Commun. 2017;8(1):1141.
Audia JE, Campbell RM. Histone modifications and Cancer. Cold Spring Harb Perspect Biol. 2016;8(4):a019521.
Sanchez GJ, Richmond PA, Bunker EN, Karman SS, Azofeifa J, Garnett AT, et al. Genome-wide dose-dependent inhibition of histone deacetylases studies reveal their roles in enhancer remodeling and suppression of oncogenic super-enhancers. Nucleic Acids Res. 2018;46(4):1756–76.
Razi S, Haghparast A, Chodari Khameneh S, Ebrahimi Sadrabadi A, Aziziyan F, Bakhtiyari M, et al. The role of tumor microenvironment on cancer stem cell fate in solid tumors. Cell Commun Signal. 2023;21(1):143.
de Visser KE, Joyce JA. The evolving tumor microenvironment: from cancer initiation to metastatic outgrowth. Cancer Cell. 2023;41(3):374–403.
West AC, Johnstone RW. New and emerging HDAC inhibitors for cancer treatment. J Clin Investig. 2014;124(1):30–9.
Zhang P, Brinton LT, Williams K, Sher S, Orwick SA-O, Tzung-Huei L et al. Targeting DNA damage repair functions of two histone deacetylases, HDAC8 and SIRT6, sensitizes Acute myeloid leukemia to NAMPT Inhibition. (1557–3265 (Electronic)).
Verza FA, Das U, Fachin AL, Dimmock JR, Marins M. Roles of histone deacetylases and inhibitors in Anticancer Therapy. 2020;12(6):1664.
Mondello P, Tadros S, Teater M, Fontan L, Chang AY, Jain N, et al. Selective inhibition of HDAC3 targets synthetic vulnerabilities and activates Immune Surveillance in Lymphoma. Cancer Discov. 2020;10(3):440–59.
Lasko LM, Jakob CG, Edalji RP, Qiu W, Montgomery D, Digiammarino EL, et al. Discovery of a selective catalytic p300/CBP inhibitor that targets lineage-specific tumours. Nature. 2017;550(7674):128–32.
Welti J, Sharp A, Brooks N, Yuan W, McNair C, Chand SN, et al. Targeting the p300/CBP Axis in Lethal prostate Cancer. Cancer Discov. 2021;11(5):1118–37.
Cai LY, Chen SJ, Xiao SH, Sun QJ, Ding CH, Zheng BN, et al. Targeting p300/CBP attenuates Hepatocellular Carcinoma Progression through Epigenetic Regulation of Metabolism. Cancer Res. 2021;81(4):860–72.
Veazey KJ, Cheng D, Lin K, Villarreal OD, Gao G, Perez-Oquendo M, et al. CARM1 inhibition reduces histone acetyltransferase activity causing synthetic lethality in CREBBP/EP300-mutated lymphomas. Leukemia. 2020;34(12):3269–85.
Li X, Su X, Liu R, Pan Y, Fang J, Cao L, et al. HDAC inhibition potentiates anti-tumor activity of macrophages and enhances anti-PD-L1-mediated tumor suppression. Oncogene. 2021;40(10):1836–50.
Buglio D, Khaskhely NM, Voo KS, Martinez-Valdez H, Liu Y-J, Younes A. HDAC11 plays an essential role in regulating OX40 ligand expression in Hodgkin lymphoma. Blood. 2011;117(10):2910–7.
Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.
Chen C, Grennan K, Badner J, Zhang D, Gershon E, Jin L, et al. Removing batch effects in analysis of expression microarray data: an evaluation of six batch adjustment methods. PLoS ONE. 2011;6(2):e17238.
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. 2022;13.
Guo J, Wang Z, Wu J, Liu M, Li M, Sun Y, et al. Endothelial SIRT6 is vital to prevent hypertension and Associated Cardiorenal Injury through Targeting Nkx3.2-GATA5 signaling. Circul Res. 2019;124(10):1448–61.
Liu X, Tan Y, Zhang C, Zhang Y, Zhang L, Ren P, et al. NAT10 regulates p53 activation through acetylating p53 at K120 and ubiquitinating Mdm2. EMBO Rep. 2016;17(3):349–66.
Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010;26(12):1572–3.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47–e.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.
Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–d92.
Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype Relationships and Predictors of response to checkpoint blockade. Cell Rep. 2017;18(1):248–62.
Sotiriou C, Wirapati P, Loi S, Harris A, Fox S, Smeds J, et al. Gene expression profiling in breast Cancer: understanding the molecular basis of histologic Grade to improve prognosis. JNCI: J Natl Cancer Inst. 2006;98(4):262–72.
Mariathasan S, Turley SJ, Nickles D, Castiglioni A, Yuen K, Wang Y, et al. TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature. 2018;554(7693):544–8.
Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang T-H, et al. The Immune Landscape of Cancer. Immunity. 2018;48(4):812–30e14.
Jie M, Wu Y, Gao M, Li X, Liu C, Ouyang Q, et al. CircMRPS35 suppresses gastric cancer progression via recruiting KAT7 to govern histone modification. Mol Cancer. 2020;19(1):56.
Zhou Y, Jin X, Ma J, Ding D, Huang Z, Sheng H, et al. HDAC5 loss impairs RB repression of pro-oncogenic genes and confers CDK4/6 inhibitor resistance in Cancer. Cancer Res. 2021;81(6):1486–99.
Yu Z, Zhou X, Wang X. Metabolic reprogramming in hematologic malignancies: advances and clinical perspectives. Cancer Res. 2022;82(17):2955–63.
Pradella D, Naro C, Sette C, Ghigna C. EMT and stemness: flexible processes tuned by alternative splicing in development and cancer progression. Mol Cancer. 2017;16(1):8.
Koushyar S, Powell AG, Vincan E, Phesse TJ. Targeting wnt signaling for the treatment of gastric Cancer. Int J Mol Sci. 2020;21(11).
Frontzek F, Staiger AM, Wullenkord R, Grau M, Zapukhlyak M, Kurz KS, et al. Molecular profiling of EBV associated diffuse large B-cell lymphoma. Leukemia. 2023;37(3):670–9.
Liu Y, Zhou X, Wang X. Targeting the tumor microenvironment in B-cell lymphoma: challenges and opportunities. J Hematol Oncol. 2021;14(1):125.
Ennishi D, Hsi ED, Steidl C, Scott DW. Toward a new molecular taxonomy of diffuse large B-cell lymphoma. Cancer Discov. 2020;10(9):1267–81.
Leivonen SK, Friman T, Autio M, Vaittinen S, Jensen AW, D’Amore F et al. Characterization and clinical impact of the tumor microenvironment in post-transplant aggressive B-cell lymphomas. Haematologica. 2023.
Tian X, Shen H, Li Z, Wang T, Wang S. Tumor-derived exosomes, myeloid-derived suppressor cells, and tumor microenvironment. J Hematol Oncol. 2019;12(1):84.
Kim HJ, Ji YR, Lee YM. Crosstalk between angiogenesis and immune regulation in the tumor microenvironment. Arch Pharm Res. 2022;45(6):401–16.
Tay RE, Richardson EK, Toh HC. Revisiting the role of CD4(+) T cells in cancer immunotherapy-new insights into old paradigms. Cancer Gene Ther. 2021;28(1–2):5–17.
Sungur CM, Murphy WJ. Positive and negative regulation by NK cells in cancer. Crit Rev Oncog. 2014;19(1–2):57–66.
Wang Y, Chen W, Lian J, Zhang H, Yu B, Zhang M, et al. The lncRNA PVT1 regulates nasopharyngeal carcinoma cell proliferation via activating the KAT2A acetyltransferase and stabilizing HIF-1α. Cell Death Differ. 2020;27(2):695–710.
Lu D, Song Y, Yu Y, Wang D, Liu B, Chen L, et al. KAT2A-mediated AR translocation into nucleus promotes abiraterone-resistance in castration-resistant prostate cancer. Cell Death Dis. 2021;12(8):787.
Haque ME, Jakaria M, Akther M, Cho DY, Kim IS, Choi DK. The GCN5: its biological functions and therapeutic potentials. Clin Sci (Lond). 2021;135(1):231–57.
Tzelepis K, Koike-Yusa H, De Braekeleer E, Li Y, Metzakopian E, Dovey OM, et al. A CRISPR dropout screen identifies genetic vulnerabilities and therapeutic targets in Acute myeloid leukemia. Cell Rep. 2016;17(4):1193–205.
We acknowledge TCGA, GTEx, and GEO database for providing their platforms and contributors to upload their significant datasets.
This study was supported by National Natural Science Foundation (No.82270200, No.82170189, No.82070203, No.81800194, No.81770210); Key Research and Development Program of Shandong Province (No.2018CXGC1213); China Postdoctoral Science Foundation (No. 2021T1404223); Translational Research Grant of NCRCH (No.2021WWB02, No.2020ZKMB01); Shandong Provincial Natural Science Foundation (ZR2021YQ51); Taishan Scholars Program of Shandong Province; Shandong Provincial Engineering Research Center of Lymphoma; Academic Promotion Programme of Shandong First Medical University (No. 2019QL018).
Ethics approval and consent to participate
This study was approved by the Medical Ethical Committee of Shandong Provincial Hospital, and written informed consent in accordance with the Declaration of Helsinki was obtained from each patient and healthy donors. Informed consent was obtained from all patients and healthy donors in this study. All methods were carried out in accordance with declaration of Helsinki.
Consent for publication
The authors have no relevant conflicts.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yu, Z., Ding, M., Cai, Y. et al. Histone regulator KAT2A acts as a potential biomarker related to tumor microenvironment and prognosis of diffuse large B cell lymphoma. BMC Cancer 23, 934 (2023). https://doi.org/10.1186/s12885-023-11401-4