
==== Front
Int J Med Sci
Int J Med Sci
ijms
International Journal of Medical Sciences
1449-1907
Ivyspring International Publisher Sydney

10.7150/ijms.96134
ijmsv21p2158
Research Paper
Multiple Comprehensive Analyses Identify Lysine Demethylase KDM as a Potential Therapeutic Target for Pancreatic Cancer
Shen Wan-Jou 1#
Kao Hsuan-Min 2#
Wang Chih-Yang 345
Kousar Rubina 6
Lin Jing-Shan 6
Ko Ching-Chung 789
Lin Hung-Yun 35101112
Ta Hoang Dang Khoa 4
Anuraga Gangga 413
Xuan Do Thi Minh 14
Kumar Sachin 315
Dey Sanskriti 3
Ly Ngoc Phung 31617
Wang Wei-Jan 618✉
1 Graduate Institute of Biomedical Sciences, College of Medicine, China Medical University, Taichung 40402, Taiwan.
2 Department of Geriatric, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Chiayi City 600566, Taiwan.
3 Graduate Institute of Cancer Biology and Drug Discovery, College of Medical Science and Technology, Taipei Medical University, Taipei 11031, Taiwan
4 Ph.D. Program for Cancer Molecular Biology and Drug Discovery, College of Medical Science and Technology, Taipei Medical University and Academia Sinica, Taipei 11031, Taiwan
5 TMU Research Center of Cancer Translational Medicine, Taipei Medical University, Taipei 11031, Taiwan.
6 Department of Biological Science and Technology, China Medical University, Taichung 40676, Taiwan
7 Department of Medical Imaging, Chi-Mei Medical Center, Tainan, Taiwan
8 Department of Health and Nutrition, Chia Nan University of Pharmacy and Science, Tainan, Taiwan
9 School of Medicine, College of Medicine, National Sun Yat-Sen University, Kaohsiung, Taiwan.
10 Cancer Center, Wan Fang Hospital, Taipei Medical University, Taipei 11031, Taiwan.
11 Traditional Herbal Medicine Research Center of Taipei Medical University Hospital, Taipei Medical University, Taipei 11031, Taiwan.
12 Pharmaceutical Research Institute, Albany College of Pharmacy and Health Sciences, Rensselaer, NY 12144, USA.
13 Department of Statistics, Faculty of Science and Technology, Universitas PGRI Adi Buana, Surabaya, East Java 60234, Indonesia.
14 Faculty of Pharmacy, Van Lang University, 69/68 Dang Thuy Tram Street, Ward 13, Binh Thanh District, Ho Chi Minh City 70000, Vietnam.
15 Faculty of Biotechnology and Applied Sciences, Shoolini University of Biotechnology and Management Sciences, Himachal Pradesh, India.
16 Natural Product Research Center, Korea Institute of Science and Technology (KIST), Gangneung 25451, Republic of Korea.
17 Division of Bio-Medical Science and Technology, KIST School, University of Science and Technology (UST), Seoul 02792, Republic of Korea.
18 Research Center for Cancer Biology, China Medical University, Taichung 40402, Taiwan.
✉ Corresponding author: cvcsky@cmu.edu.tw (W.-J.W.).
# These authors equally contributed to this study.

Competing Interests: The authors have declared that no competing interest exists.

2024
12 8 2024
21 11 21582169
9 3 2024
22 7 2024
© The author(s)
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Pancreatic cancer (PC) is a challenging and heterogeneous disease with a high mortality rate. Despite advancements in treatment, the prognosis for PC patients remains poor, with a high chance of disease recurrence. Biomarkers are crucial for diagnosing cancer, predicting patient prognosis and selecting treatments. However, the current lack of effective biomarkers for PC could contribute to the insufficiency of existing treatments. These findings underscore the urgent need to develop novel strategies to fight this disease. This study utilized multiple comprehensive bioinformatic analyses to identify potential therapeutic target genes in PC, focusing on histone lysine demethylases (KDMs). We found that high expression levels of KDM family genes, particularly KDM1A, KDM5A and KDM5B, were associated with improved overall survival in the cohort. Furthermore, the infiltration of various immune cells, including B cells, neutrophils, CD8+ T cells, dendritic cells, and macrophages, was positively correlated with KDM1A, KDM5A, and KDM5B expression. Moreover, MetaCore pathway analysis revealed interesting connections between KDM1A and the cell cycle and proliferation, between KDM5A and DNA damage and double-strand break repair through homologous recombination, and between KDM5B and WNT/β-catenin signaling. These findings suggest that KDM1A, KDM5A and KDM5B may serve as promising biomarkers and therapeutic targets for PC, a disease of high importance due to its aggressive nature and urgent need for novel biomarkers to improve diagnosis and treatment.
==== Body
pmcIntroduction

Pancreatic cancer (PC), a type of cancer that affects the digestive system, is associated with a 5-year survival rate of less than 5% and a median survival duration of 6 months in patients in whom the cancer cannot be surgically removed 1. Moreover, the majority of PCs are of the exocrine type, with ductal adenocarcinoma being the most common, accounting for 80-90% of all cases. 2. PC ranks as the 13th most common type of cancer globally, with a total incidence of 458,000 cases, and is the seventh most lethal cancer type, accounting for approximately 50,000 deaths in the United States alone 3. The symptoms of PC may not be distinct, making early diagnosis challenging with current screening techniques 4. Therefore, understanding the molecular mechanisms that may induce treatment resistance, morbidity, and mortality may aid in devising relevant clinical strategies for treating PC 5.

Biomarkers are crucial in cancer diagnosis, prognosis prediction, and treatment selection. Samples of tissue, plasma, urine, and saliva can be analyzed to identify relevant biomarkers, providing a noninvasive, cost-effective, and efficient screening method 6. The development of novel biomarkers can improve early detection, prognosis prediction, and treatment selection, ultimately leading to better patient outcomes 7-11. Since PC is often diagnosed at an advanced stage when treatment options are limited, the discovery of new biomarkers could allow for earlier intervention and improved survival rates. Moreover, new biomarkers may also help monitor the effectiveness of treatments and develop more targeted therapies for this challenging disease 12-14.

Histone lysine demethylases (KDMs) have been demonstrated to remove methylation marks from lysine residues in histone tails, affecting gene transcription during development and cancer progression. Moreover, their activity can either promote or suppress tumor growth, depending on the cellular mechanism and specific isoforms of KDMs. The KDM family consists of 19 members: KDM1A, KDM1B, KDM2A, KDM2B, KDM3A, KDM3B, KDM4A, KDM4B, KDM4C, KDM4D, KDM4E, KDM5A, KDM5B, KDM5C, KDM5D, KDM6A, KDM6B, KDM7A, and KDM8 15. KDM1A is highly expressed in various human cancers, such as triple-negative breast cancer and cancer stem cells (CSCs). In prostate cancer, the overexpression of KDM1B is associated with poor prognosis. Inhibitors targeting KDM1B significantly reduce the development of prostate tumors 16. Additionally, in breast cancer, KDM2A is overexpressed and associated with poor prognosis 17. Similarly, in ovarian cancer, overexpression of KDM2B is associated with poor prognosis, whereas a decrease in KDM2B expression is correlated with reduced proliferation and migration of cancer cells 18. In lung cancer, the overexpression of KDM1A and KDM3A has been associated with lymphatic metastasis, pathological grade, and clinical stage. Another study reported that high KDM3B expression in breast cancer cells and tissues leads to rapid growth and invasion of cancer cells 19. Furthermore, KDM4A overexpression has been observed in malignant pleural mesothelioma, an aggressive malignancy with a poor prognosis 20. Moreover, a study indicated that KDM4C expression substantially increased in gastric CSCs and correlated favorably with relevant biomarkers 21. In addition, KDM4D promotes leukemia progression by activating the expression of myeloid leukemia proteins 22. Although evidence suggests that the KDM family plays a crucial role in cancers, an overview of KDMs in PC is still lacking.

In this study, we investigated the overexpression patterns of different members of the KDM family in PC and focused on all KDM family members through comprehensive bioinformatics analysis (Fig. 1.). To understand the importance of the KDM family in PC development and prognosis, we analyzed relevant data from the Genotype-Tissue Expression databases and The Cancer Genome Atlas (TCGA) 23. The differential expression levels of various KDM family genes in PC were observed via data obtained from the Gene Expression Profiling Interactive Analysis 2 (GEPIA2) and University of Alabama at Birmingham CANcer (UALCAN) databases. This observation indicated their potential oncogenic and tumor-suppressive activities in PC. We further evaluated the diagnostic and prognostic value of KDM family members and assessed their protein expression levels in PC. Overall, we found that KDM1A, KDM5A, and KDM5B were highly expressed and critical for prognostic assessment and immune infiltration in PC. These findings enhance our understanding of the role of KDMs in PC and provide novel biomarkers for PC treatment, enabling the prediction of prognosis on the basis of KDM functions.

Materials and Methods

GEPIA2

We collected data on the expression of hub genes in both normal and malignant tissues from the GEPIA2 database (http://gepia2.cancer-pku.cn/index.html). We then investigated the differences in hub gene expression to validate the data obtained from TCGA. Additionally, we analyzed the expression patterns of the target genes in various cancers and the correlations between target gene expression and tumor grade via data from the GEPIA2 database 23-26.

UALCAN

The UALCAN online portal (http://ualcan.path.uab.edu) provides data on patient survival, cancer transcriptomes, and proteomics. Using data from TCGA, UALCAN users can analyze the expression of protein-coding genes and their impact on patient survival across 33 types of cancers 27-29. In this study, we examined the expression profiles of all 19 members of the KDM family in patients with PC.

Survival analysis

The KM survival curve estimates the likelihood of patient survival over time by dividing time into short intervals. We utilized the KM plotter (http://kmplot.com/analysis/) tool to assess the predictive value of the mRNA expression levels of various KDM family genes in PC. Patients with PC were classified into high- and low-expression groups on the basis of the median values of mRNA expression, and their overall survival (OS) was evaluated. The KM database provides data on hazard ratios, 95% confidence intervals, median mRNA expression, and p values, making it a valuable resource for evidence-based survival data. Additionally, we determined the log-rank p values and hazard ratios (HRs) with 95% confidence intervals 30-32.

Tumor Immune Estimation Resource

The Tumor Immune Estimation Resource (TIMER: http://timer.comp-genomics.org) online server is a comprehensive tool for systematically assessing infiltrating immune cells in different cancers. We utilized this tool to estimate the abundance of six types of immune cells: CD4+ T cells, neutrophils, B cells, dendritic cells, CD8+ T cells, and macrophages. Using the TIMER database, we examined the correlations between KDM1A, KDM5A, and KDM5B expression and the infiltration of various immune cells. TCGA data were analyzed via TIMER to assess KDM1A, KDM5A, and KDM5B expression patterns in both PC and normal tissues 33-35. The gene symbols for KDM1A, KDM5A, and KDM5B are plotted on the x-axis, whereas the related marker genes' symbols are plotted on the y-axis. Gene expression levels are represented as log2 RSEM values.

Protein-protein interaction network and gene coexpression

We utilized the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING: https://string-db.org/cgi/input) database to investigate protein‒protein interactions (PPIs). KDM1A, KDM5A, and KDM5B sequences were queried in the STRING database to identify their PPIs. On the basis of the overexpression patterns of these proteins, we constructed a PPI network 36-38.

Functional component and potential pathway analyses

We employed cBioPortal for Cancer Genomics (https://www.cbioportal.org) software to analyze the coexpression of selected genes with other genes. The relevant data were subsequently uploaded to the Database for Annotation, Visualization, and Integrated Discovery (DAVID: https://david.ncifcrf.gov), a tool for systematically gathering biological data from multiple genes 39-41, and an online platform for data analysis and visualization (http://www.bioinformatics.com.cn/). Additionally, we obtained data from the Kyoto Encyclopedia of Genes and Genomes (KEGG), which assigns biological functions to genes and genomes, whereas DAVID provides information about gene products, processes, and functions. To gain a deeper understanding of the biological processes associated with the differentially expressed genes (DEGs), we employed a curated pathway enrichment analysis using MetaCore™ (Clarivate Analytics, UK) (https://clarivate.com/products/biopharma/research-development/early-research-intelligence-solutions/). For each time point, we uploaded the DEGs relative to control cells into MetaCore™ for pathway analysis. The Pathway Maps tool was then used to identify significantly enriched pathways associated with these DEGs. An adjusted p value threshold of less than 0.05 was applied for the analysis, and all expression changes were included. The results are presented as -log(p value) values, ranked by their statistical significance 42-44.

Immunohistochemical (IHC) analysis

PC tissue microarrays, sourced from TissueArray (PA2082a), were analyzed to explore the expression of KDM1A, KDM5A, and KDM5B via IHC staining. The microarrays included normal and tumor samples from PC patients 45-47. For the IHC staining, we used antibodies against KDM1A (Abcam, ab37165), KDM5A (Abcam, ab92533), and KDM5B (Abcam, ab244220). The staining procedures were performed according to the TnAlink Polymer Detection System (DAB) manufacturers' protocols (BioTnA, TAHC04D).

Statistical analysis

We conducted survival analysis via the KM plotter tool with its default settings. Recurrence-free survival was the primary focus, utilizing the automatically calculated best cutoff values and the J best probe set. Furthermore, we utilized the TCGA Pan-Cancer Atlas dataset from the cBioPortal database to acquire patient data and examine the impact of the expression of various KDM family genes on overall survival (OS). We determined all potential cutoff values between the lower and upper quartiles, selecting the best-presenting threshold as the final cutoff value. A log-rank p value less than 0.05 was considered statistically significant 48-50.

Results

Differential expression analysis of the KDM family genes in PC

KDMs are critical enzymes that regulate gene expression by removing methyl groups from lysine residues on histones. This process is essential for maintaining the dynamic balance of genomic expression, ensuring genome integrity, and supporting epigenetic inheritance. Aberrant KDM expression and activity in cancer are pivotal in tumor progression and drug resistance. However, no study has explored the overall role and expression of the entire KDM family in PC. This study utilized an extensive range of publicly available databases and web-based analysis tools. The bioinformatics method we employed is depicted in Figure 1. To identify the role of KDMs in PC, we investigated the expression patterns of various KDM family genes in both normal and cancerous tissues from PC patients via the GEPIA2 database. Our findings revealed significantly upregulated expression levels (p < 0.5) of KDM1A, KDM1B, KDM2A, KDM2B, KDM3B, KDM4A, KDM5A, KDM5B, and KDM6B in PC tissue, suggesting their potential oncogenic role (Fig. 2). However, KDM3A, KDM4B, KDM4C, KDM4D, KDM4E, KDM5C, KDM6A, KDM7A, and KDM8 were not significantly different in the tumor tissue. Furthermore, immunohistochemical data analysis of the expression of KDM5B in human tissue, which is available in the Human Protein Atlas database, revealed greater staining intensity in PC tissue than in normal tissue. These findings suggest that certain members of the KDM family could serve as promising biomarkers for PC.

The prognostic relevance of KDM family gene expression in PC and its correlation with the infiltration of immune cells

We then examined whether the prognosis of cancer patients was related to the overexpression of KDM1A, KDM1B, KDM2A, KDM2B, KDM3B, KDM4A, KDM5A, KDM5B, KDM5D, and KDM6B in PC tumor tissue. The results revealed that high expression levels of KDM family genes were associated with two opposite scenarios in PC patients: the high expression levels were linked to either extended or poorer overall survival (OS) in patients with PC. Notably, poor overall survival (OS) was observed for KDM1A, KDM5A, and KDM5B (p < 0.05). Conversely, KDM2B, KDM4B, KDM4C, KDM4D, KDM4E, KDM5D, KDM6B, and KDM8 were associated with prolonged OS (Fig. 3A). Furthermore, to assess the protein expression levels of KDM1A, KDM5A, and KDM5B in human tumor tissue, we conducted IHC on a tumor tissue microarray. In clinical PC samples, significantly greater expression of KDM1A, KDM5A, and KDM5B was observed than in normal pancreatic tissue (Fig. 3B). These results suggest that the KDM family, particularly KDM1A, KDM5A, and KDM5B, plays an oncogenic role and is associated with poor prognosis in PC patients. In addition, accumulating studies have demonstrated that immune infiltration in PC is crucial in the development of more effective immunotherapies, as it influences disease progression, treatment response, and patient outcomes 51. Therefore, we were interested in investigating whether KDM1A, KDM5A, and KDM5B expression was correlated with the enrichment of tumor-infiltrating lymphocytes in PC. Using the TIMER database, we found that the expression of the KDM1A gene is associated with the infiltration of CD8+ T cells, macrophages, neutrophils, and dendritic cells. Furthermore, KDM5A was related to the infiltration of B cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells. Additionally, KDM5B was associated with the infiltration of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells in patients with PC (Fig. 4). These results suggest that KDM1A, KDM5A, and KDM5B may serve as oncogenic biomarkers and correlate with antitumor immunity in PC.

Analysis of gene interactions and molecular functions involving KDM1A, KDM5A, and KDM5B

To gain further insight into the molecular regulation of KDMs, we used the STRING database to identify PPIs within the KDM family. The results indicated that KDM1A and KDM5A significantly interact with HDAC1/2, which plays a critical role in cancer by regulating gene expression through histone deacetylation, thereby influencing cell proliferation, apoptosis, and metastasis 52. Moreover, we found that KDM5B interacts with HDAC1/4 to mediate miRNA expression and enhance docetaxel resistance in lung adenocarcinoma (Supplementary Fig. 1A) 53. Furthermore, we utilized gene sets coexpressed with KDM1A, KDM5A, and KDM5B from the TCGA database. Our analysis revealed several pathways associated with KDM1A that are significantly associated with pancreatic cancer, including ephrin receptor signaling, cell cycle regulation, carbon metabolism, progesterone-mediated oocyte maturation, and glycolysis/glycogenesis, as identified by MetaCore analysis and KEGG database (Fig. 5 and Supplementary Fig. 1B). Our findings provide additional insights into the functional roles of KDM5A, revealing its association with pathways related to DNA damage, the IL-6-mediated immune response, the DNA damage response, and double-strand break repair through homologous recombination (Fig. 6 and Supplementary Fig. 1B). Moreover, our study demonstrated that KDM5B is associated with WNT/β-catenin signaling in the cytoplasm, HGF signaling in cancer, sphingolipid metabolism, and the Hippo pathway (Fig. 7 and Supplementary Fig. 1B). Furthermore, we identified the top 50 pathways related to KDM1A, KDM5A and KDM5B, and the results revealed that the KDM family may contribute to multiple oncogenic pathways in PC (Supplementary Tables 1-3). These correlations offer additional evidence of their potential participation in crucial cellular processes and signaling pathways relevant to cancer biology. Taken together, in this study, we conducted a comprehensive analysis via multiple credible databases, which implied that KDMs, particularly KDM1A, 5A, and 5B, may play critical roles in the molecular pathways underlying the development and progression of PC (Fig. 7).

Discussion

While PC has a lower incidence than other malignancies, such as breast, lung, colorectal, and prostate cancers, it was the third leading cause of cancer-related death in 2021 54. By 2030, it is projected to become the second leading cause of cancer-related death in the United States 55. Despite this, there has been limited success in improving the outcomes of PC over the past two decades, underscoring the need for effective treatment and diagnostic approaches.

Although key genetic mutations, such as those in KRAS, TP53, SMAD4, and CDKN2A, have been identified in PC patients, the overall outcomes for these patients have not been markedly improved. Therefore, identifying novel biomarkers for PC treatment is imperative as the key to improving patient outcomes and extending survival. By meticulously analyzing these biomarkers and their correlation with treatment response, researchers can pave the way for more effective therapeutic strategies. This holds particular promise for KDMs, enzymes that regulate gene expression through histone methylation, offering hope for improved patient care and increased survival rates. Although their involvement in cancer is well known, the molecular functions, expression patterns, and significance of KDM family genes in patient survival in PC are still unclear. In this study, we used multiple comprehensive bioinformatics approaches to analyze KDM family proteins in PC. Our findings suggest that KDM1A, KDM5A, and KDM5B could serve as biomarkers for PC and may be targeted to improve long-term patient survival.

We found that KDM1A, also known as LSD1, is overexpressed in PC and plays a significant role in disease development and progression, which is also correlated with the glycolysis pathway. Several KDM genes, including KDM1A, KDM5A, and KDM5B, were upregulated in PC tissues compared with normal controls, indicating their potential as oncogenes. Furthermore, high KDM1A, KDM5A and KDM5B expression correlated with poorer overall survival in PC patients. These results highlight these KDM family members as potential diagnostic and prognostic markers for PC. Studies have demonstrated that glycolysis is a critical regulator in cancer and plays a role in cancer stemness and chemoresistance 56. Chemotherapy is known to trigger cancer stemness, which is characterized by enhanced DNA repair capacity and contributes to temozolomide (TMZ) resistance in glioblastoma (GBM). High expression of KDM1A has been linked to therapy resistance and recurrence in GBM. Inhibiting KDM1A has been shown to sensitize glioma stem cells (GSCs) to TMZ, indicating its involvement in TMZ resistance 57. Furthermore, KDM1A acts as a coregulator for hormone receptors (AR, ER, GR, MR) and may control ERβ expression in ovarian cancer (OCa). KDM1A knockdown enhances ERβ expression, and combination therapy with a KDM1A inhibitor and ERβ agonist reduces OCa cell viability and invasion and promotes apoptosis in various OCa models 58.

In PC, inducing KDM2B protein expression has been shown to promote tumor progression through interaction with Polycomb group proteins to repress developmental genes. Furthermore, researchers have demonstrated that KDM2B collaborates with KDM5A and that Myc activates genes involved in ribosomal and mitochondrial functions, which are crucial for the metabolism and growth of PC​ 59​. Our study revealed that KDM2B gene expression is correlated with favorable prognosis in patients with PC. Therefore, gene expression and protein expression may be regulated by different mechanisms in PC, which requires further investigation. Intriguingly, our data revealed a link between KDM expression and immune cell infiltration in PC. Specifically, KDM1A, KDM5A, and KDM5B expression correlated with the infiltration of multiple immune cell types, including CD8+ T cells, B cells, macrophages and dendritic cells. These findings suggest that targeting these KDM genes may target cancer cells directly and modulate the antitumor immune response. Further investigations are needed to explore this possibility and its potential for immunotherapeutic strategies in PC. These enzymes play significant roles in a wide range of physiological and pathological processes, indicating their complex and context-dependent functions. Moreover, our results revealed that KDM5A and KDM5B are correlated with DNA damage repair and the WNT/β-catenin signaling pathway. Consistently, recent studies have shown that KDM5A can be upregulated in tumor tissues compared with normal prostate tissues and may promote malignancy through the regulation of chromosome instability 60, 61. Furthermore, research has demonstrated that KDM5B regulates the androgen receptor and interacts with the PI3K/AKT pathway, contributing to the progression of prostate cancer 62. In addition to KDM2B, we also found that KDM5D and KDM6B are correlated with good prognosis in PC patients. According to a previous study, KDM5D plays a critical role in prostate cancer by regulating the response to androgen deprivation therapy (ADT) and taxane treatments. Loss of KDM5D increases the expression of MYBL2, a transcription factor for which elevated levels are associated with poorer clinical outcomes, leading to resistance to ADT and taxanes 63. However, the role of KDM5D in PC remains unclear. In addition, KDM6B acts as a tumor suppressor in PC by demethylating histone H3 lysine 27, thereby controlling the expression of CEBPA, a key gene associated with cancer progression. Inhibition of KDM6B expression enhances PC cell aggressiveness, increases metastasis, and disrupts cellular senescence mediated by oncogenic KRAS, highlighting its critical role in maintaining cellular control in early cancer stages 64. Moreover, KDM6B has also been demonstrated to promote neuroblastoma cell differentiation by removing the repressive chromatin marker H3K27me3. KDM6B expression is induced by retinoic acid via HOXC9 and inhibits tumor cell growth by promoting cell differentiation 65. These findings highlight the complex roles of KDMs in cancer biology, particularly in PC, where they can influence both the suppression and activation of different pathways critical for cancer progression.

In addition, previous studies have demonstrated that KDM5A and KDM5C are overexpressed in gemcitabine-resistant pancreatic cancer cells and that their expression is regulated by CD44, a stem cell marker linked to chemoresistance. Moreover, studies have shown that KDM5 family members are involved in several critical oncogenic pathways, including epithelial‒mesenchymal transition (EMT) 66. These findings suggest that members of the KDM family could serve as new biomarkers or potential therapeutic targets for overcoming drug resistance in PC. In our study, MetaCore analysis suggested that KDM5A may contribute to the regulation of the IL-6 pathway. Consistent with this finding, KDM5A has been demonstrated to regulate osteosarcoma cell viability through the IL-6/JAK/STAT3 pathway, highlighting its potential role in modulating this signaling axis 67. Moreover, a previous study showed that KDM5B, by interacting with TIEG1, enhances the repression of SMAD7 thereby augmenting TGF-β signaling and potentially suppressing tumorigenesis 68. Interestingly, our bioinformatics analysis revealed that KDM5B is involved in multiple pathways, such as the WNT, HGF, and TGF-β pathways, in PC. Taken together, these findings showed that both KDM5A and KDM5B might act as protumorigenic factors in cancer, highlighting their potential as targets for therapeutic intervention. However, further research is needed to fully elucidate the molecular mechanisms underlying their roles, especially in PC progression, and to explore their potential as therapeutic targets in PC treatment.

Conclusions

We conducted comprehensive bioinformatics analyses using data from the TCGA and the Genotype-Tissue Expression databases. We observed high expression levels of KDM family genes in PC tissue. Our prognostic analysis suggested that KDM1A, KDM5A, and KDM5B could be promising biomarkers and therapeutic targets for PC. In addition, we discovered that the elevated expression of KDM1A, KDM5A, and KDM5B in PC may be correlated with poor prognosis and significantly decreased immune cell infiltration. These data suggest that targeting KDMs could increase the effectiveness of immunotherapy for PC.

Supplementary Material

Supplementary figure and tables.

Bioinformatical analyses and data mining were performed by the staff members of the Bioinformatics Core Facility at Taipei Medical University. This manuscript was edited by Wallace Academic Editing. The authors gratefully acknowledge the support provided by Daniel P. Chamberlin from the Office of Research and Development at Taipei Medical University in enhancing the linguistic quality of the manuscript.

Funding

This research was funded by grants from the Taiwanese National Science and Technology Council (grant numbers: NSTC 112-2320-B-039 -012 -MY3 to W-J. W.), China Medical University (grant number: CMU109-MF-55 to W-J. W.), Taipei Medical University (TMU; grant number: TMU-108-AE1-B16 to C-Y.W.), Ditmanson Medical Foundation Chia-Yi Christian Hospital (grant numbers: R111-056). This work was additionally supported by the TMU Research Center of Cancer Translational Medicine (Featured Areas Research Center Program) within the framework of the Higher Education Sprout Project of the Taiwanese Ministry of Education.

Ethics approval and consent to participate

This study is approved by the ethics committee of China Medical University.

Data availability statement

cBioPortal for Cancer Genomics: https://cbioportal.org, Human Protein Atlas: https://www.proteinatlas.org, Kaplan-Meier plotter: https://kmplot.com, and MetaCore: https://portal.genego.com. The datasets used and analyzed in the present study are available from the corresponding author on reasonable request.

Figure 1 The bioinformatics analysis flowchart depicts the processing of publicly available data from TCGA and GTEx databases.

Figure 2 The expression levels of KDM family proteins in pancreatic cancer (PC). KDM family mRNA expression analysis in PCs from the GEPIA2 database. Box plot (GEPIA2): normal tissue, gray; tumor tissue, red. *, p<0.01.

Figure 3 The prognostic value of the mRNA levels of the KDM family in patients with pancreatic cancer (PC) was determined via Kaplan‒Meier plotter analysis. (A) The hazard ratio (HR) indicates the prognostic value for PC patients. The log [rank p] test was used to determine the level of prognostic significance, with a value of p<0.05 considered significant. High expression of KDM1A/5A/5B was significantly associated with poor prognosis, indicating poorer outcomes. Conversely, the HRs of KDM2B/5D/6B were significantly lower, suggesting better prognostic outcomes in patients with pancreatic cancer. (B) Immunohistochemical patterns of KDM1A, KDM5A, and KDM5B expression in normal and tumor tissues from patients with PC.

Figure 4 Correlation analysis was performed to assess the relationships between differentially expressed KDM family genes and immune cell infiltration via the TIMER database. The figure illustrates the associations of KDM1A, KDM5A, and KDM5B gene expression with tumor purity and markers of tumor-infiltrating immune cells, including B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells. Spearman correlations were used to assess the relationships between these KDM genes and the mentioned immune cells, with statistical significance set at p<0.05.

Figure 5 MetaCore pathway analysis of the coexpression gene network involving KDM1A in pancreatic cancer. We extracted the top 1000 genes coexpressed with KDM1A from the TCGA database. We conducted a pathway analysis and generated a pathway list ordered by the -log p value. The "Ephrin receptors pathway" ranked highest in the biological process category. The figure illustrates interactions between genes and proteins, with symbols representing proteins and arrows depicting protein interactions (green for activation, red for inhibition).

Figure 6 Expression of the KDM5A signaling pathway in pancreatic cancer (MetaCore). We used the MetaCore platform to analyze genes coexpressed with KDM5A from the associated TCGA database. Our analysis revealed that "DNA damage" was the most common biological process.

Figure 7 MetaCore enrichment pathway analysis of genes coexpressed with KDM5B in pancreatic cancer. We used the MetaCore platform to analyze genes coexpressed with KDM5B from the associated TCGA database. Our analysis revealed that "WNT/β-catenin signaling" ranked highest in the biological process analysis.
==== Refs
1 Balaban EP Mangu PB Yee NS Locally Advanced Unresectable Pancreatic Cancer: American Society of Clinical Oncology Clinical Practice Guideline Summary J Oncol Pract 2017 13 265 9 28399382
2 Turanli B Yildirim E Gulfidan G Arga KY Sinha RJJoPM Current state of “omics” biomarkers in pancreatic cancer J Pers Med 2021 11 127 33672926
3 Moffat GT Epstein AS O'Reilly EM Pancreatic cancer-A disease in need: Optimizing and integrating supportive care Cancer 2019 125 3927 35 31381149
4 Wolrab D Jirásko R Cífková E Höring M Mei D Chocholoušková M Lipidomic profiling of human serum enables detection of pancreatic cancer Nat Commun 2022 13 1 16 34983933
5 Murugan NJ Voutsadakis IA Proteasome regulators in pancreatic cancer World J Gastrointest Oncol 2022 14 38 35116102
6 Armstrong EA Beal EW Chakedis J Paredes AZ Moris D Pawlik TM Exosomes in pancreatic cancer: From early detection to treatment J Gastrointest Surg 2018 22 737 50 29423813
7 Lin JC Liu TP Yang PM CDKN2A-Inactivated Pancreatic Ductal Adenocarcinoma Exhibits Therapeutic Sensitivity to Paclitaxel: A Bioinformatics Study J Clin Med 2020 9 4019 33322698
8 Lin TY Wang PW Huang CH Yang PM Pan TL Characterizing the Relapse Potential in Different Luminal Subtypes of Breast Cancers with Functional Proteomics Int J Mol Sci 2020 21 6077 32846884
9 Liu LW Hsieh YY Yang PM Bioinformatics Data Mining Repurposes the JAK2 (Janus Kinase 2) Inhibitor Fedratinib for Treating Pancreatic Ductal Adenocarcinoma by Reversing the KRAS (Kirsten Rat Sarcoma 2 Viral Oncogene Homolog)-Driven Gene Signature J Pers Med 2020 10 130 32947833
10 Yang PM Hsieh YY Du JL Yen SC Hung CF Sequential Interferon β-Cisplatin Treatment Enhances the Surface Exposure of Calreticulin in Cancer Cells via an Interferon Regulatory Factor 1-Dependent Manner Biomolecules 2020 10 643 32326356
11 Yang PM Lin LS Liu TP Sorafenib Inhibits Ribonucleotide Reductase Regulatory Subunit M2 (RRM2) in Hepatocellular Carcinoma Cells Biomolecules 2020 10 117 31936661
12 Wang S Zheng Y Yang F Zhu L Zhu XQ Wang ZF The molecular biology of pancreatic adenocarcinoma: translational challenges and clinical perspectives Signal Transduct Target Ther 2021 6 249 34219130
13 Ta HDK Minh Xuan DT Tang WC Anuraga G Ni YC Pan SR Novel Insights into the Prognosis and Immunological Value of the SLC35A (Solute Carrier 35A) Family Genes in Human Breast Cancer Biomedicines 2021 9 1804 34944621
14 Ta HDK Wang WJ Phan NN An Ton NT Anuraga G Ku SC Potential Therapeutic and Prognostic Values of LSM Family Genes in Breast Cancer Cancers (Basel) 2021 13 4902 34638387
15 Sterling J Menezes SV Abbassi RH Munoz L Histone lysine demethylases and their functions in cancer Int J Cancer 2021 148 2375 88 33128779
16 Tang D He J Dai Y Geng X Leng Q Jiang H Targeting KDM1B-dependent miR-215-AR-AGR2-axis promotes sensitivity to enzalutamide-resistant prostate cancer Cancer Gene Ther 2022 29 543 557 33854217
17 Chen J-Y Li C-F Lai Y-S Hung W-C Lysine demethylase 2A expression in cancer-associated fibroblasts promotes breast tumour growth Br J Cancer 2021 124 484 93 33024266
18 Kuang Y Xu H Lu F Meng J Yi Y Yang H Inhibition of microRNA let-7b expression by KDM2B promotes cancer progression by targeting EZH2 in ovarian cancer Cancer Sci 2021 112 231 42 33091189
19 Hu A Hong F Li D Xie Q Chen K Zhu L KDM3B-ETF1 fusion gene downregulates LMO2 via the WNT/β-catenin signaling pathway, promoting metastasis of invasive ductal carcinoma Cancer Gene Ther 2022 29 215 24 33828234
20 Lapidot M Case AE Weisberg EL Meng C Walker SR Garg S Essential role of the histone lysine demethylase KDM4A in the biology of malignant pleural mesothelioma (MPM) Br J Cancer 2021 125 582 92 34088988
21 Lang T Xu J Zhou L Zhang Z Ma X Gu J Disruption of KDM4C-ALDH1A3 feed-forward loop inhibits stemness, tumorigenesis and chemoresistance of gastric cancer stem cells Signal Transduct Target Ther 2021 6 1 3 33384407
22 Wu W Cao X Mo L Overexpression of KDM4D promotes acute myeloid leukemia cell development by activating MCL-1 Am J Transl Res 2021 13 2308 19 34017391
23 Wang CC Lin WH Ku SC Shen WJ Ta HDK Anuraga G Circadian rhythm-related factors of PER and CRY family genes function as novel therapeutic targets and prognostic biomarkers in lung adenocarcinoma Aging (Albany NY) 2022 14 9056 89 36385012
24 Chou CW Hsieh YH Ku SC Shen WJ Anuraga G Khoa Ta HD Potential Prognostic Biomarkers of OSBPL Family Genes in Patients with Pancreatic Ductal Adenocarcinoma Biomedicines 2021 9 1601 34829830
25 Wang CC Shen WJ Anuraga G Hsieh YH Khoa Ta HD Xuan DTM Penetrating Exploration of Prognostic Correlations of the FKBP Gene Family with Lung Adenocarcinoma J Pers Med 2022 13 49 36675710
26 Kao TJ Wang CY Hsu TI Wu YH Chuang JY Huang CC TDP-43 is essential for Eph receptor-class-specific spinal motor axon trajectory into the limb Neurosci Res 2022 176 31 9 34687812
27 Chandrashekar DS Karthikeyan SK Korla PK Patel H Shovon AR Athar M UALCAN: An update to the integrated cancer data analysis platform Neoplasia 2022 25 18 27 35078134
28 Ku SC Liu HL Su CY Yeh IJ Yen MC Anuraga G Comprehensive analysis of prognostic significance of cadherin (CDH) gene family in breast cancer Aging (Albany NY) 2022 14 8498 567 36315446
29 Hagerling C Gonzalez H Salari K Wang CY Lin C Robles I Immune effector monocyte-neutrophil cooperation induced by the primary tumor prevents metastatic progression of breast cancer Proc Natl Acad Sci U S A 2019 116 21704 14 31591235
30 Choy TK Wang CY Phan NN Khoa Ta HD Anuraga G Liu YH Identification of Dipeptidyl Peptidase (DPP) Family Genes in Clinical Breast Cancer Patients via an Integrated Bioinformatics Approach Diagnostics (Basel) 2021 11 1204 34359286
31 Chen WC Hsu HP Li CY Yang YJ Hung YH Cho CY Cancer stem cell marker CD90 inhibits ovarian cancer formation via β3 integrin Int J Oncol 2016 49 1881 9 27633757
32 Chen YL Lee KT Wang CY Shen CH Chen SC Chung WP Low expression of cytosolic NOTCH1 predicts poor prognosis of breast cancer patients Am J Cancer Res 2022 12 2084 101 35693094
33 Li T Fan J Wang B Traugh N Chen Q Liu JS TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells Cancer Res 2017 77 e108 e10 29092952
34 Wang CY Xuan DTM Ye PH Li CY Anuraga G Ta HDK Synergistic suppressive effects on triple-negative breast cancer by the combination of JTC-801 and sodium oxamate Am J Cancer Res 2023 13 4661 77 37970352
35 Xuan DTM Wu CC Wang WJ Hsu HP Ta HDK Anuraga G Glutamine synthetase regulates the immune microenvironment and cancer development through the inflammatory pathway Int J Med Sci 2023 20 35 49 36619229
36 Szklarczyk D Gable AL Lyon D Junge A Wyder S Huerta-Cepas J STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets Nucleic Acids Research 2018 47 D607 D13
37 Anuraga G Lang J Xuan DTM Ta HDK Jiang JZ Sun Z Integrated bioinformatics approaches to investigate alterations in transcriptomic profiles of monkeypox infected human cell line model J Infect Public Health 2024 17 60 9 37992435
38 Xuan DTM Yeh IJ Wu CC Su CY Liu HL Chiao CC Comparison of Transcriptomic Signatures between Monkeypox-Infected Monkey and Human Cell Lines J Immunol Res 2022 2022 3883822 36093436
39 Long T Liu Z Zhou X Yu S Tian H Bao Y Identification of differentially expressed genes and enriched pathways in lung cancer using bioinformatics analysis Mol Med Rep 2019 19 2029 40 30664219
40 Wu YH Yeh IJ Phan NN Yen MC Liu HL Wang CY Severe acute respiratory syndrome coronavirus (SARS-CoV)-2 infection induces dysregulation of immunity: in silico gene expression analysis Int J Med Sci 2021 18 1143 52 33526974
41 Cheng LC Kao TJ Phan NN Chiao CC Yen MC Chen CF Novel signaling pathways regulate SARS-CoV and SARS-CoV-2 infectious disease Medicine (Baltimore) 2021 100 e24321 33607766
42 Wang CY Chiao CC Phan NN Li CY Sun ZD Jiang JZ Gene signatures and potential therapeutic targets of amino acid metabolism in estrogen receptor-positive breast cancer Am J Cancer Res 2020 10 95 113 32064155
43 Wu YH Yeh IJ Phan NN Yen MC Hung JH Chiao CC Gene signatures and potential therapeutic targets of Middle East respiratory syndrome coronavirus (MERS-CoV)-infected human lung adenocarcinoma epithelial cells J Microbiol Immunol Infect 2021 54 845 57 34176764
44 Kao TJ Wu CC Phan NN Liu YH Ta HDK Anuraga G Prognoses and genomic analyses of proteasome 26S subunit, ATPase (PSMC) family genes in clinical breast cancer Aging (Albany NY) 2021 13 17970 34329194
45 Lottini T Duranti C Iorio J Martinelli M Colasurdo R D'Alessandro FN Combination therapy with a bispecific antibody targeting the hERG1/β1 integrin complex and gemcitabine in pancreatic ductal adenocarcinoma Cancers 2023 15 2013 37046674
46 Tram VTN Khoa Ta HD Anuraga G Dung PVT Xuan DTM Dey S Dysbindin Domain-Containing 1 in Prostate Cancer: New Insights into Bioinformatic Validation of Molecular and Immunological Features Int J Mol Sci 2023 24 11930 37569304
47 Weng TY Huang SS Yen MC Lin CC Chen YL Lin CM A novel cancer therapeutic using thrombospondin 1 in dendritic cells Mol Ther 2014 22 292 302 24127010
48 Wang CC Shen WJ Anuraga G Khoa Ta HD Xuan DTM Chen ST Novel Potential Therapeutic Targets of PTPN Families for Lung Cancer J Pers Med 2022 12 1947 36556168
49 Weng TY Wu HF Li CY Hung YH Chang YW Chen YL Homoharringtonine induced immune alteration for an Efficient Anti-tumor Response in Mouse Models of Non-small Cell Lung Adenocarcinoma Expressing Kras Mutation Sci Rep 2018 8 8216 29844447
50 Shahi P Slorach EM Wang CY Chou J Lu A Ruderisch A The Transcriptional Repressor ZNF503/Zeppo2 Promotes Mammary Epithelial Cell Proliferation and Enhances Cell Invasion J Biol Chem 2015 290 3803 13 25538248
51 Looi CK Chung FF Leong CO Wong SF Rosli R Mai CW Therapeutic challenges and current immunomodulatory strategies in targeting the immunosuppressive pancreatic tumor microenvironment J Exp Clin Cancer Res 2019 38 162 30987642
52 Kuang J Li P Zhai Z Fan Y Xu H Zhao C Exclusion of HDAC1/2 complexes by oncogenic nuclear condensates Mol Cancer 2024 23 85 38678233
53 Jiang M Qi F Zhang K Zhang X Ma J Xia S MARCKSL1-2 reverses docetaxel-resistance of lung adenocarcinoma cells by recruiting SUZ12 to suppress HDAC1 and elevate miR-200b Mol Cancer 2022 21 150 35864549
54 Di Federico A Mosca M Pagani R Carloni R Frega G De Giglio A Immunotherapy in Pancreatic Cancer: Why Do We Keep Failing? A Focus on Tumor Immune Microenvironment, Predictive Biomarkers and Treatment Outcomes 2022 14 2429
55 Padrón LJ Maurer DM O'Hara MH O'Reilly EM Wolff RA Wainberg ZA Sotigalimab and/or nivolumab with chemotherapy in first-line metastatic pancreatic cancer: clinical and immunologic analyses from the randomized phase 2 PRINCE trial Nat Med 2022 28 1167 1177 35662283
56 Jiang YX Siu MKY Wang JJ Leung THY Chan DW Cheung ANY PFKFB3 Regulates Chemoresistance, Metastasis and Stemness via IAP Proteins and the NF-kappaB Signaling Pathway in Ovarian Cancer Front Oncol 2022 12 748403 35155224
57 Alejo S Palacios BE Venkata PP He Y Li W Johnson JD Lysine-specific histone demethylase 1A (KDM1A/LSD1) inhibition attenuates DNA double-strand break repair and augments the efficacy of temozolomide in glioblastoma Neuro Oncol 2023 25 1249 61 36652263
58 Venkata PP Jayamohan S He Y Alejo S Johnson JD Palacios BE Pharmacological inhibition of KDM1A/LSD1 enhances estrogen receptor beta-mediated tumor suppression in ovarian cancer Cancer Lett 2023 575 216383 37714256
59 Tzatsos A Paskaleva P Ferrari F Deshpande V Stoykova S Contino G KDM2B promotes pancreatic cancer via Polycomb-dependent and -independent transcriptional programs J Clin Invest 2013 123 727 39 23321669
60 Du C Lv C Feng Y Yu S Activation of the KDM5A/miRNA-495/YTHDF2/m6A-MOB3B axis facilitates prostate cancer progression J Exp Clin Cancer Res 2020 39 223 33087165
61 Domingo-Reines J Montes R Garcia-Moreno A Gallardo A Sanchez-Manas JM Ellson I The pediatric leukemia oncoprotein NUP98-KDM5A induces genomic instability that may facilitate malignant transformation Cell Death Dis 2023 14 357 37301844
62 Li G Kanagasabai T Lu W Zou MR Zhang SM Celada SI KDM5B Is Essential for the Hyperactivation of PI3K/AKT Signaling in Prostate Tumorigenesis Cancer Res 2020 80 4633 43 32868382
63 Yoshikawa Y Stopsack KH Wang XV Chen YH Mazzu YZ Burton F Increased MYBL2 expression in aggressive hormone-sensitive prostate cancer Mol Oncol 2022 16 3994 4010 36087093
64 Yamamoto K Tateishi K Kudo Y Sato T Yamamoto S Miyabayashi K Loss of histone demethylase KDM6B enhances aggressiveness of pancreatic cancer through downregulation of C/EBPalpha Carcinogenesis 2014 35 2404 14 24947179
65 Yang L Zha Y Ding J Ye B Liu M Yan C Histone demethylase KDM6B has an anti-tumorigenic function in neuroblastoma by promoting differentiation Oncogenesis 2019 8 3 30631055
66 Wang D Zhang Y Liao Z Ge H Gungor C Li Y KDM5 family of demethylases promotes CD44-mediated chemoresistance in pancreatic adenocarcinomas Sci Rep 2023 13 18250 37880235
67 Peng D Lin B Xie M Zhang P Guo Q Li Q Histone demethylase KDM5A promotes tumorigenesis of osteosarcoma tumor Cell Death Discov 2021 7 9 33436536
68 Kim J Shin S Subramaniam M Bruinsma E Kim TD Hawse JR Histone demethylase JARID1B/KDM5B is a corepressor of TIEG1/KLF10 Biochem Biophys Res Commun 2010 401 412 6 20863814
