Abstract
Background:
The malignant progression of pancreatic cancer (PC) is frequently accompanied by intractable pain mediated through perineural invasion (PNI), yet the underlying epigenetic regulatory mechanisms remain elusive.This study aims to elucidate the role of DNA methylation in the pathogenesis of PC pain, including its interactive effects with the nervous and immune systems.
Methods:
Integrating multi-omics data from TCGA-PAAD (Pancreatic adenocarcinoma), we identified methylation driver genes (MDGs) using the MethylMix algorithm. By intersecting MDGs with pain-related gene sets and conducting multi-step regression modeling, we established a five-gene prognostic signature (PSMB8/COL17A1/BICC1/CTRC/TRIP13). Next, in order to elucidate the underlying mechanisms, we conducted differential expression analysis, protein-protein interaction network analysis, functional enrichment analysis, and single-cell sequencing. Additionally, we quantified immune infiltration using CIBERSORT and TIMER.
Results:
Pain-related MDGs are enriched in immune regulation, extracellular matrix reorganization, and cation channel activity, constituting the “immune-neural axis” of epigenetic regulation. The prognostic five-gene signature significantly stratifies patient survival risk (HR = 3.83, p = 1.4e−8), with its methylation levels positively correlated with CD4+ T cell infiltration and negatively correlated with dendritic cells. Model-derived differentially expressed genes exhibited dual immune-neural tropism at single-cell resolution, prominently enriched in presynaptic signaling and synaptic vesicle cycling. Mechanistically, MDGs orchestrate pain progression through PNI-associated neural remodeling and K+ channel-mediated neuronal hypersensitization.
Conclusion:
This study establishes a visceral pain model centered on pancreatic parenchymal nociception rather than secondary neural effects, and for the first time proposes an interconnected regulatory network linking epigenetic modifications, immune reprogramming, and neural plasticity, revealing dual pain pathogenesis mechanisms: (1) immune microenvironment reshaping that potentiates neuroinflammation, and (2) direct ion channel regulation enhancing neuronal excitability. These findings provide a mechanistic foundation for developing methylation-based prognostic biomarkers and multimodal analgesic therapeutic strategies targeting the immuno-neural nexus.
1 Introduction
The global incidence of PC continues to rise, characterized by its high malignancy and early-stage concealment. Over 80% of patients are diagnosed at an advanced stage, losing the opportunity for curative surgery. The high postoperative metastasis rate and a 5-year survival rate of less than 10% place PC as the fourth leading cause of cancer-related deaths worldwide (Kalra and Meltzer, 2025). Pain, as a primary symptom, affects approximately 90% of patients, significantly reducing quality of life and potentially triggering anxiety, depression, and even suicidal tendencies (Tarasiuk et al., 2023; Shrestha et al., 2024). In 2018, the WHO for the first time classified “chronic cancer pain” as a distinct disease category (), with effective pain management shown to prolong survival (Ozcan et al., 2019). However, the mechanisms underlying pain remain unclear, and clinical management relies heavily on empirical approaches, highlighting the urgent need for molecular mechanism research (; ).
The study found that the main cause of PC pain is neuropathic pain triggered by tumor cells invading nerves, which is considered a result of PNI (Zhu et al., 2024). PNI defined as cancer cells invading along nerves or the interstitial spaces of the neural sheath, perineurium, and endoneurium, is a special way for cancer to spread to distant sites (Zhu et al., 2024). The incidence of PNI in PC is nearly 100%, which is negatively related to the survival rate, and is a risk factor for R1 resection and recurrence (Selvaggi et al., 2022; Ozaki et al., 1999). The sensory nerve endings distributed within internal organs, such as the pancreas, and blood vessel walls are capable of perceiving stimuli including osmotic pressure, temperature, and pathological injuries (). Pain signals are transmitted via nerve fibers to the cell bodies of cranial/spinal ganglia (), and ultimately ascend to the central nervous system through afferent nerves. This conduction pathway confirms that visceral pain signals originate from pathological changes within the parenchyma of the organs (Ten Barge et al., 2025). However, previous studies on PNI-related pain models in PC have primarily focused on the sciatic nerve and dorsal root ganglia affected by tumor cells, concentrating on the pain signal transmission process while ignoring the initial site of pain generation in the viscera (Miura et al., 2018).
Recent studies have found that Schwann cells (key components of peripheral nerves) are significantly enriched in pancreatic tumor tissues (Xue et al., 2023). By interacting with tumor/immune/stromal cells, Schwann cells activate tumor-neural system dialogue, promoting metastasis and driving the process of carcinogenesis (Xu W. et al., 2024). On the one hand, tumor cells obtain more growth signals, remodel metabolism, and evade immune surveillance to promote survival by inducing their own innervation (). When PC cells invade nerves, they disrupt the normal structure of nerve fibers, leading to abnormal nerve signal transmission (Yang et al., 2024). Abnormal nerve impulses are transmitted along damaged nerve fibers, triggering pain sensations. During PNI, PC secretes nerve growth factors such as NGF, which may increase the sensitivity of nerve endings, known as hyperalgesia (Xu J. et al., 2024; Liu et al., 2024). PNI in PC may also lead to changes in neural plasticity, including nerve hyperplasia and hypertrophy, which further exacerbate abnormal nerve signal transmission and hyperalgesia (Xu et al., 2023). On the other hand, the sympathetic nerves, parasympathetic nerves, and sensory nerves within the pancreas voluntarily adapt to the progressive malignant process and establish bidirectional communication with the tumor to support its growth (). Besides direct effects, nerves regulate the progression of PC by actively regulating the functions of stroma and immune components, and the involvement of tumor microenvironment (TME) complicates the bidirectional interaction (). Taking the immune system as an example, immune cell activation and mediator release induce/maintain cancer pain (), and pain-related immune disorders (mainly inflammation) and immunosuppression delay pain resolution (Zhao et al., 2023), ultimately impairing patients’ immune function and prognosis. Combined with previous studies, it was found that the immune system and nervous system participate in regulating the occurrence and development of pain. However, little is known about how the two systems are related to the initial signal of pain generation.
In the study of pain-related mechanisms, epigenetic modifications involving DNA methylation have also garnered significant attention. In neuropathic pain models, DNMT-mediated hyper-methylation of the Oprm1/Kcna2 gene promoters leads to gene silencing, which triggers increased neuronal excitability and hyperalgesia (Sun et al., 2017; Sun et al., 2019). In rat dorsal root ganglion (DRG), DNMT1-driven downregulation of Cnr1 gene methylation can weaken its inhibitory effect on TRPV1, exacerbating visceral pain (). In CFA-induced chronic inflammatory pain, DNMT3b inhibition-induced de-methylation of the NGF promoter upregulates NGF expression and maintains the pain state by promoting C/EBPα binding (Yuan et al., 2020). Although oral cancer studies suggest that antitumor gene de-methylation and neurotrophin hyper-methylation are involved in the PNI process (), the intrinsic association between DNA methylation and PNI in PC pain remains unknown.
This study integrated TCGA PC methylation data with the MSigDB pain gene set, and identified 26 pain-related MDGs using the MethylMix algorithm to construct a pain risk scoring model. After stratifying patients based on this model and combining it with single-cell analysis, it was revealed that tumor-derived pain signals drive neuroplastic changes and immune microenvironment reprogramming through PNI, leading to a molecular axis of poor prognosis (Figure 1). This study innovatively focuses on the primary visceral pain site (distinct from secondary neural effects) and systematically analyzes the complete mechanism of pain signaling from its initial generation to transmission. It provides a theoretical foundation for the early diagnosis of PC pain and targeted multi-modal analgesic strategies for the “immune-neural hub”.
FIGURE 1
2 Materials and methods
2.1 Data collection
We downloaded RNA sequencing data, metabolomics data, somatic mutation data, and The Cancer Genome Atlas (TCGA; https://portal.gdc.cancer.gov/repository) (Zhang and Wang, 2015) from 178 PC tissues and four normal pancreatic tissues. The pain-related gene set HP_PAIN (M38128) was obtained by searching the keyword “pain” from the Molecular Signatures Database (MSigDB), which contains 835 genes related to pain (MSIGDB_URL:https://www.gsea-msigdb.org/gsea/msigdb/human/geneset/HP_PAIN). The mRNA expression profiles and clinical information of the GSE183795 dataset were downloaded from the GEO database as a validation cohort. This dataset includes microarray gene expression profiles of 139 pancreatic tumors, 102 adjacent non-tumor tissues, and three normal pancreases from donors with pancreatic ductal adenocarcinoma patients.
2.2 Identify pain - related MDGs in PC and conduct functional enrichment analysis
By using the beta hybrid model of the MethylMix package, we can identify sample subpopulations of PC with different DNA methylation compared with normal tissues, so as to identify differential and functional DNA methylation. Functional DNA methylation refers to 471 MDGs of PC based on significant negative correlation of matched gene expression data. By intersecting with pain-related genes in the MsigDB database, 26 pain-related MDGs were ultimately obtained.
Conduct Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis on 26 genes, and analyze gene characteristics from aspects of molecular function, biological process and cellular component as well as the interactions and regulations of their pathways in biological systems.
2.3 Construction and validation of risk scoring model
In this study, using the TCGA cohort, we employed the R package glmnet to integrate survival time, survival status, and gene expression data, and utilized the lasso-cox method for regression analysis. Additionally, we set up 10-fold cross-validation to obtain the optimal model. We set the λ value to be 0.0675911321389309, and finally obtained a model formula constructed by 10 genes as follows: RiskScore = 2.94730303909491e-05*CTRC+0.0038277720604039*TRIP13+0.00767681073907178*PSMB8-0.00156609272886264*IRF4+5.34857478531571e-05*GNE+0.00651372422293902*BICC1+0.00119961881628371*COL17A1-0.00169964562145739*FERMT3-0.00549988980019879*KCNJ2-0.000188824126427788*MAFB.
After multivariate survival analysis, five pain-related MDGs were found to be significantly associated with overall survival. These five genes are CTRC, TRIP13, PSMB8, BICC1, and COL17A1. We used the R package maxstat (Maximally selected rank statistics with several p-value approximations version: 0.7-25) to calculate the optimal cutoff value of RiskScore. We set the minimum number of samples in a group to be greater than 25% and the maximum number of samples in a group to be less than 75%. Finally, the optimal cutoff value was obtained as 1.81837560981286. Based on this, patients were divided into high-risk and low-risk groups. Then, we further used the survfit function in the R package survival to analyze the prognostic differences between the two groups, and the logrank test method was used to evaluate the significance of prognostic differences between different groups of samples. Eventually, we observed significant prognostic differences (p = 1.4e−8). Analyzed the relationship between risk scoring and pathological grading and staging, and observed the survival curves of the five main pain-related MGDs that make up the risk scoring model We selected the mRNA expression profiles and clinical information of the GSE183795 dataset downloaded from the GEO database as the validation cohort to observe the survival curves and prognostic analysis of high-risk and low-risk subgroups, thus verifying the risk score model.
2.4 Analysis of key pain-related MDGs
MethylMix is an algorithm for identifying highly methylated and hypo-methylated genes associated with diseases. MethylMix identifies methylation states based on the β-mixture model and compares them with normal DNA methylation states. MethylMix uses a new statistical quantity, namely, the difference methylation value or DM value, which is defined as the difference between the methylation state and the normal methylation state. Finally, the matched gene expression data are used to identify functional methylation states other than differences by focusing on the methylation changes that affect gene expression. This study, based on the TCGA PC gene expression matrix and DNA methylation, analyzed 471 PC MDGs by using the MethylMix package. The methylation states of five key pain-related MDGs and their correlations with mRNA expression levels were analyzed, and the gene expression differences in tumors and normal tissues were also analyzed.
2.5 Differential gene selection and weighted gene Co-expression network analysis (WGCNA)
Based on the risk score model, TCGA PC patients were divided into high-risk and low-risk subtypes. Limma is a differential expression screening method based on the generalized linear model. Here, we used the R package limma (version 3.40.6) for differential analysis to obtain differentially expressed genes between the high-risk and low-risk subtypes. A total of 5,962 differentially expressed genes were identified with a fold change of 1.5 and p < 0.05, including 3,886 upregulated genes and 2,076 downregulated genes. After analyzing their chromosomal locations, all differentially expressed genes underwent WGCNA to identify differentially expressed genes in the most relevant modules. Combined with survival time and survival status, the most relevant module cyan was determined, which encompassed 89 genes. Enrichment analysis of GO and KEGG was performed on the genes in this module.
2.6 Construction of protein-protein interaction (PPI) networks and identification of hub genes
To study the differentially expressed genes associated with pain risk subtypes in PC methylated driver genes, we input the selected cyan module genes into the String database to construct a PPI network, with a cutoff value set at 0.400, and visualized it using Cytoscape software. Subsequently, we utilized the Molecular Complex Detection (MCODE) tool in Cytoscape to analyze the gene interaction information. By applying criteria of degree cutoff = 2, node score cutoff = 0.2, and K-Core = 2, we identified hub modules within the PPI network, comprising eight genes, and conducted a visualization analysis of their gene functions.
2.7 Analysis of tumor mutation burden
Tumor Mutational Burden (TMB) as an emerging biomarker, can assist in predicting patients’ response to immunotherapy. In this study, we analyzed the mutational landscape of the top 20 most frequently mutated genes in the pain-related risk subtypes of PC patients, demonstrating the specific mutation conditions of all samples. With the help of “maftools” package in R software, we calculated the tumor mutation conditions of PC patients. The “survival” package was used to determine the optimal cut-off value, and the samples were divided into high TMB group and low TMB group based on TMB differential analysis. Subsequent survival analysis was conducted regarding the TMB.
2.8 Immune cell infiltration and gene set enrichment analysis
IOBR (Zeng et al., 2021) is a computing tool for immune tumor biology research. Here, based on the expression profile of TCGA PC patients, the R software package IOBR selected xCell method to calculate 64 immune cell infiltration scores and immune, matrix and microenvironment scores for each sample. The CIBERSORT method calculated 22 immune cell infiltration scores for each sample and analyzed the correlation between immune cell infiltration scores and 10 hub genes. For Gene Set Variation Analysis (GSVA), we used the R software package to calculate the enrichment score of each sample in the gene set from GSVA (version 1.40.1). We predefined the gene rank. Specifically, we first used gene expression to reach the spectrum, using the Hanzelmann et al. method, and downloaded the hallmark subset from the Molecular Signatures Database (http://www.gsea-msigdb.org/gsea/downloads.jsp) to evaluate relevant pathways and molecular mechanisms, Set the minimum gene set to 5 and the maximum gene set to 5,000, and calculate the enrichment scores of each sample in each gene set. Finally, the enrichment score matrix was obtained.
2.9 Analysis of the correlation between gene methylation, copy number variation and immune infiltration
GSCALite is a web-based platform for gene set cancer analysis. GSCALite integrates cancer genomics data of 33 cancer types from TCGA, drug response data from GDSC and CTRP, and normal tissue data from GTEx for an integrated data analysis workflow in gene set analysis. In this study, the key genes under investigation are input into the GSCALite platform to analyze their interrelationships with immune infiltration. TIMER is a comprehensive database whose primary function is to systematically analyze six types of tumor-infiltrating immune cells (B cells, CD4+ T cells, CD8+T cells, neutrophils, macrophages, and dendritic cells) in different cancer types through the TIMER algorithm. At the same time, it analyzes the relationship between gene expression and tumor purity. Genes that are highly expressed in the tumor microenvironment are negatively correlated with tumor purity, while genes that are highly expressed in tumor cells are positively correlated with tumor purity.
2.10 Single-cell sequencing analysis
CDCP is a comprehensive platform for single-cell data integration, sharing and analysis. Users can obtain detailed information about samples in the datasets included in the CDCP single-cell data platform online, and are allowed to use tSNE (t-Distributed Stochastic Neighbor Embedding) cell dimensionality reduction maps and clustering analysis maps of different cell types to visualize each single-cell dataset. In this study, the CDCP online platform was utilized to analyze the gene expression of key genes in single cells in the pancreas or nerves and visualize them for comparison.
3 Results
3.1 Construction and validation of risk scoring model for pain-related MDGs in PC
Emerging evidence has established epigenetic regulation, particularly DNA methylation, as a central modulator of chronic pain pathophysiology (). To systematically investigate methylation-driven mechanisms underlying PC-associated pain, we performed multi-dimensional bioinformatics analysis of PC samples from TCGA (n = 178). High-resolution methylation arrays and RNA-seq data were integrated through the MethylMix algorithm, revealing 471 aberrantly methylated driver genes (FDR < 0.05) with 333 hypermethylated and 138 hypomethylated candidates (Supplementary Table 1). Cross-referencing these epigenetic drivers with a curated database of 835 nociception-associated genes (PainGenesDB v4.0) identified 26 high-confidence pain-related MDGs (15 hypermethylated, 11 hypomethylated; Figure 2A; Supplementary Tables 2, 3). GO enrichment analysis revealed that these 26 genes are enriched in extracellular matrix structural constituents, DNA-binding transcription activator activity, ligand-gated cation channel activity, and other molecular functions. In biological processes, they are enriched in the regulation of immune system processes, immune function modulation, and cell differentiation regulation. In terms of cellular components, they are enriched in the extracellular matrix and integral components of the plasma membrane (Figure 2B). KEGG enrichment analysis showed that these genes are associated with signaling pathways such as T cell receptor signaling pathway and Th17 cell differentiation (Figure 2C). Thus, it can be seen that PC pain-related MDGs may play a role in immune response and cell differentiation and proliferation, suggesting an intrinsic link between methylation and pain in PC.
FIGURE 2
The lasso-cox regression analysis was conducted on the screened 26 pain-related MDGs, followed by a 10-fold cross-validation to obtain the optimal risk scoring model (Figures 2D,E). The λ value was determined to be 0.0675911321389309. The formula for the model constructed using the final 10 selected genes is as follows:
RiskScore = 2.94730303909491e-05*CTRC+0.0038277720604039*TRIP13+0.00767681073907178*PSMB8-0.00156609272886264*IRF4+5.34857478531571e-05*GNE+0.00651372422293902*BICC1+0.00119961881628371*COL17A1-0.00169964562145739*FERMT3-0.00549988980019879*KCNJ2-0.000188824126427788*MAFB.
After multivariate survival analysis, five pain-related MDGs were identified as significantly correlated with overall survival (Figure 2F). These five genes are: PSMB8, COL17A1, BICC1, CTRC, and TRIP13. The hazard ratio (HR) is 3.83. Using the R package “maxstat”, the optimal cutoff value for RiskScore was calculated to be 1.81837560981286. Based on this value, patients were divided into high and low groups, showing significant prognostic differences (p = 1.4e−8) (Figure 2G). The ROC curve revealed that the AUCs for patient survival at 1, 3, and 5 years were 0.76, 0.84, and 0.84, respectively (Figure 2H). The prognostic heatmap analyzed the relationship between different risk scores and patient follow-up time, events, and changes in gene expression. It was observed that as the risk score increased, the survival rate of patients decreased (Figure 2I). The risk scoring model was validated using the mRNA expression profiles and clinical information from the GSE183795 dataset as a validation cohort (Figure 2J).The HR of GSE183795 is 1.95 (p = 2.3e−3), and the ROC curve reveals that the AUCs for patient survival at 1, 3, and 5 years were 0.55, 0.64 and 0.67, respectively, which is significant (Supplementary Figure 1).
3.2 Risk model is significantly related to pathological staging of PC
Based on the tumor staging and TNM staging information of PC patients from TCGA, combined with the risk score derived from pain-related MDGs, the correlation between the risk score and the clinical information of PC was evaluated. In the tumor staging of PC, stage Ⅰ tumors generally do not have regional lymph node metastasis or distant metastasis, while stage Ⅱ and Ⅲ tumors have different degrees of regional lymph node metastasis, and stage Ⅳ tumors are generally accompanied by distant metastasis. As can be seen from the results in Figure 3A, the high-risk score is most correlated with PC stage Ⅱ. Stage Ⅱ tumors have 1–3 regional lymph node metastases. It is speculated that pain-related MDGs play an important role in the early metastasis process of tumors. Observing the results of the correlation between the TNM staging and the risk score of PC patients, it can be seen that pain is correlated with the size of the primary tumor (T staging) and the situation of regional lymph node metastasis to different degrees. No significant association with distant metastasis (Figures 2B–D). It can be seen that the risk score model constructed in this study is more suitable for the stage where PC begins to spread but has not yet had distant metastasis. A high-risk score indicates that the tumor has a tendency to metastasize. The Cox multivariate analysis obtained five key pain-related MDGs. The Kaplan-Meier curves of these five genes in TCGA PC patients were drawn respectively, and it can be seen that PSMB8, COL17A1, BICC1 and TRIP13 are significantly correlated with the prognosis (Figures 2E–I). The metastatic propensity signature (high COL17A1/PSMB8, low CTRC) suggests: PNI facilitation via collagen remodeling (COL17A1), immune-evasion priming through proteasomal antigen processing (PSMB8),protective trypsin depletion in tumor microenvironment (CTRC).
FIGURE 3
3.3 Characterization of key pain-related DNA MDGs
The MethyMix algorithm can identify genes’ low/high methylation states through the beta mixture model and find DNA MDGs in diseases by analyzing the correlation of the expression levels of corresponding genes. The expression patterns (left) and negative correlation maps (right) of key pain-related MDGs in different samples are shown in Figure 4 (PSMB8: Figure 4A, COL17A1: Figure 4C, BICC1: Figure 4E, CTRC: Figure 4G, TRIP13: Figure 4I). Combined with the differential methylation values of each driver gene (it’s the average difference in methylation between tumor samples and normal samples, Supplementary Table 4), it can be seen that PSMB8 and CTRC are not methylated in normal samples of PC in the TCGA database. Observing the mRNA expression levels of key pain-related MDGs in pancreatic tumor and normal samples, significant differences can be seen in the mRNA expression levels of PSMB8 and TRIP13 (PSMB8: Figure 4B, COL17A1: Figure 4D, BICC1: Figure 4F, CTRC: Figure 4H, TRIP13: Figure 4J).
FIGURE 4
3.4 Differential gene identification and WGCNA of risk scoring model
In the TCGA cohort, limma analysis identified a total of 3,886 upregulated and 2,076 downregulated differentially expressed genes (DEGs) (|logFC| > 1.5 and p-value < 0.05) between the two pain-related risk subtypes (Figure 5A; Supplementary Table 5). Chromosomal mapping revealed chromosome 19 enrichment (q = 0.003) of upregulated genes, while downregulated genes clustered on 7q21.3 (PAX/POU domain loci; Figure 5B). The heatmap visualization of the first 50 differentially expressed genes is shown in Figure 5C. The goodSamplesGenes method in the R package WGCNA was used to remove outlier genes and samples, and then the scale-free co-expression network was further constructed using WGCNA. Pearson correlation matrix and average linkage method were applied to all pairs of genes, and then a weighted adjacency matrix was constructed. The adjacency relationship was transformed into a topological overlap matrix (TOM), which can measure the network connectivity of a gene, defined as the sum of its adjacency relationships with all other genes, for network gene ratio, and the corresponding dissimilarity (1-TOM) was calculated (Figures 5D,E). Genes with similar expression profiles were classified into gene modules. According to the dissimilarity measure based on TOM, average linkage hierarchical clustering was carried out, and the minimum size (gene group) of the gene dendrogram was set to be 30 (Figures 5F,G). The sensitivity was set as 3. In order to further analyze the modules, the dissimilarity of characteristic genes of the modules was calculated, a cutting line was selected for the module dendrogram, and some modules were merged. In addition, modules with a distance less than 0.25 were also merged. Finally, 17 co-expression modules were obtained, among which the Gray module was considered as unable to be reasonably combined with any other model. Visualizing the correlations of the 17 co-expression gene modules resulted in Figure 5H, from which it can be seen that the darkgreen module has the highest correlation with the black module. By comprehensively analyzing the correlations of the 17 co-expression gene modules with tumor staging, it can be seen that the cyan module is most correlated with tumor stage and TNM staging. Therefore, it was selected as the key module for subsequent analysis, which contains a total of 89 genes (Figure 5I; Supplementary Table 6).
FIGURE 5
3.5 Enrichment analysis of key module genes and construction of protein-protein interaction networks
Subsequently, GO and KEGG enrichment analyses were performed on the crucial cyan module. The analysis results of biological processes indicated that the genes in the cyan module were concentrated on the establishment of extracellular protein localization and the secretion and transport of hormones, including insulin, neurotransmitters, and peptide hormones (Figure 6A). The analysis results of cellular components showed that the genes in the cyan module were focused on the process of signal transduction, including neuronal projection terminals, axon terminals, external vesicle membranes, synaptic vesicle membranes, transport vesicle membranes, potassium channel complexes, etc., all of which are closely related to signal transduction in the nervous system (Figure 6B). The molecular function analysis results of the cyan module genes revealed a close correlation with ATPase-coupled ion transmembrane transporter activity, potassium channel activity, voltage-gated cation channel activity, and synaptotagmin binding (Figure 6C). The KEGG analysis results demonstrated that the genes in the cyan module were concentrated on oxidative phosphorylation, mannose-type O-glycan biosynthesis, and diabetes or synaptic vesicle recycling-related aspects (Figure 6D). Among them, the pathogenic mechanisms related to sugar metabolism and diabetes are closely associated with pancreatic tissues, while synaptic vesicle recycling-related functions account for the largest proportion in the KEGG analysis of the cyan module, which is consistent with the results of GO enrichment analysis. By visualizing the synaptic vesicle recycling process and marking the genes involved in the cyan module in green, it can be seen that the genes in the cyan module participate in multiple key steps of synaptic vesicle recycling in signal transduction and are crucial for signal transduction (Figure 6E). The PPI network diagram of the cyan module genes was obtained by importing them into the Cytoscape software (Figure 6F). Using the MCODE module in the Cytoscape software to contract the genes, eight key differentially expressed genes were obtained, namely, PTPRN2, CPE, SNAP25, ABCC8, PTPRN, CHGB, SCG5, and PCSK2, these genes are all downregulated in the risk scoring model (Figure 6G).
FIGURE 6
3.6 Gene diversity of TCGA PC and single cell analysis of key MDGs
The diversity of biological DNA sequences plays a crucial role in the occurrence and development of diseases by affecting gene expression products or gene regulatory processes, resulting in different traits in individuals. Single nucleotide polymorphism (SNP) is a major form of genetic diversity and can occur in both coding and non-coding regions of the genome. We utilized the CMplot package to visualize the distribution of SNPs in the genomes of PC patients (Figure 7A). Subsequently, we identified the top 2 genes with the highest mutation rates among TCGA PC patients (Figure 7B). It is noteworthy that mutations in KRAS, TP53, SMAD4, and CDKN2A were the most common in both subgroups, among which KRAS and TP53 were the most dominant, with more than 5% in both groups. Missense mutations were the most common, followed by nonsense mutations. By summarizing the somatic mutation situation of cyan module genes, it was found that mutations were mainly SNPs, mainly concentrated in C>T and C>A (Figure 7C). The characteristic risk score was positively correlated with TMB (R = 0.2, p = 0.0096, Figure 7D left), while survival analysis showed that a higher level of TMB was significantly associated with a longer OS (p = 0.02, Figure 7D right).
FIGURE 7
Subsequently, we compared the copy number variations (CNVs) or mutations of five pain-related MDGs with various immune cells to observe the impact of gene changes in pain-related MDGs on immune infiltration in PC. Here, B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells (DCs) were selected to analyze the correlation between CNVs or mutations of genes and these immune cells. Among them, the CNVs of PSMB8, COL17A1, and BICC1 showed no obvious correlation with these six types of immune cells. The CNV of CTRC was significantly positively correlated with macrophages, and the TRIP13 was significantly positively correlated with neutrophils (Supplementary Figure 2A), indicating that the CNV changes of these five key pain-related MDGs had a relatively small impact on immune infiltration in PC. Combining the mutation changes of five pain-related MDGs, it can be seen that PSMB8 was significantly positively correlated with B cells, CD8+ T cells, and neutrophils; COL17A1 was significantly positively correlated with B cells, CD8+ T cells, and CD4+ T cells, and significantly negatively correlated with neutrophils and DCs; no significant correlation was found between the mutations of BICC1 and CTRC and immune cells; the mutation of TRIP13 was significantly positively correlated with CD8+ T cells and CD4+ T cells, and significantly negatively correlated with neutrophils and DCs (Supplementary Figure 2B). By comprehensively analyzing five MDGs and various immune cells, it was found that low CNVs of these five key MDGs were closely associated with CD4+ T cells, B cells, natural killer (NK) cells, etc., while high CNVs were significantly correlated with CD8-naive, natural regulatory T cells (nTreg), neutrophils, etc (Supplementary Figure 2C). Further analysis of the correlation between Single Nucleotide Variant (SNV) of five MDGs and immune infiltration showed that there was no significant correlation between them (Supplementary Figure 2D). At the same time, we also investigated the correlation between the mRNA expression levels of these genes and immune cells (Supplementary Figure 3). Combined with the above research findings, it can be concluded that the mutation is the most significant aspect in their DNA variations in terms of the impact on immune infiltration among these five key MDGs.
In order to observe the distribution of five pain-related MDGs in PC tissues or the immune microenvironment, we utilized the CDCP database and selected five public single-cell analysis data sets, namely, GSM6567161, GSM6567159, GSM663326, GSM728974 and GSM663329, to observe the gene distribution (Figure 7E). Among them, GSM6567161 and GSM6567159 mainly focus on the immune landscape of PC; GSM663326 and GSM663329 are concerned with the comparison between PC and corresponding normal tissues; GSM728974 focuses on tumor microenvironmental changes. It can be seen that the expressions of these genes mainly exist in PP cells, epithelial cells, mesenchymal cells, and acinar cells, which are several types of pancreatic tissue cells.
3.7 Analysis of immune function and drug sensitivity of differentially expressed genes
Based on the aforementioned research, it can be discovered that the pain risk subtypes constructed by the MDGs exhibit significant correlations with cellular immune functions. Through the GSEA-hallmark enrichment analysis of 5,962 differentially expressed genes in the risk scoring model, it was shown that the high-risk subgroup was associated with interferon-α response, interferon-β response, immune response, TGF-BETA signaling and IL6-JAK-STAT3 signaling, while the low-risk subgroup was associated with pancreas beta cell, spermatogenesis, myogenesis, oxidative phosphorylation and hedgehog signaling (Figure 8A). The analysis results of 22 types of immune cells for differentially expressed genes indicated that there were significant differences in dendritic cells resting, Macrophages M and Macrophages M1 cells (Figures 8B,C). After performing immune cell infiltration and GSVA enrichment scoring on the previously screened MCODE differentially expressed genes, it was found that they were closely related to multiple immune cells (Figure 8D). By analyzing the correlations between PC ESTIMATEScore, ImmuneScore, StromalScore and TumorPurity and the risk score, it can be seen that the high-risk score was significantly positively correlated with ESTIMATEScore and StromalScore, but negatively correlated with TumorPurity (Figure 8E). The analysis of MCODE gene immune cells showed a positive correlation between gene expression and immune cell function (Supplementary Figure 4). Due to the downregulation of gene expression in high-risk subtypes, it suggests a decrease in cellular immune function and an increase in the proportion of immune infiltration. This indicates that in the high-risk subgroup of PC, the proportions of immune cells and stromal components are higher, suggesting a higher likelihood of immune infiltration. Combined with the drug sensitivity analysis results from CTRP and GDSC websites, the drug sensitivity predictions of different tumor chemotherapy drugs can be seen (Figure 8F).
FIGURE 8
3.8 Single-cell analysis of key MCODE genes in the pancreas or nerves
In the development mechanism of PC pain, besides inflammatory pain caused by changes in immune function, neuropathic pain caused by tumor tissue compressing nerves is also crucial. Combined with previous studies, it can be seen that the risk scoring model constructed by pain-related MDGs may have certain correlations with early metastasis of pancreatic tumors. The selected differentially expressed genes play important roles in synapses related to nerve signal transmission and voltage-gated ion channels. Therefore, we further analyzed their expression distributions in the pancreas and nerves (Figure 9). GSM6567159 is a single-cell analysis data set of the immune landscape of PC, and GSM523652 is a single-cell transcription atlas of the human spinal cord. It can be seen that MCODE is expressed in different degrees not only in pancreatic tissues but also more in nerve cells, providing further evidence for the influence of perineural infiltration of PC on pain.
FIGURE 9
4 Discussion
In recent years, the global incidence of PC has shown a sharp increase. As one of the leading causes of cancer deaths, the diagnosis and treatment of PC are particularly important (Li et al., 2025). Pain is a significant factor affecting the quality of life for PC patients and is crucial for patient outcomes (Tarasiuk et al., 2023). The main mechanisms of PC pain include direct compression and infiltration of surrounding nerves by the tumor, as well as the release of inflammatory mediators (Zhu et al., 2024). Previous studies have indicated that cancer pain primarily involves genetic and epigenetic regulation. Epigenetic regulation, including DNA methylation, histone modifications, and non-coding RNAs, affects neuroinflammation, neuronal sensitization, and pain transmission by regulating the expression of pain-related genes or the activity of signaling pathways (Ni et al., 2024). Based on differences in the degree of PNI and inflammatory mediators/chemokines between PC and adjacent normal tissues, we extracted relevant data from TCGA for PC and adjacent normal tissues. Using the MethylMix software package, 471 MDGs were identified, and five pain-related MDGs (PSMB8/COL17A1/BICC1/CTRC/TRIP13) significantly associated with overall survival (OS) were finally selected for prognostic stratification of PC patients (HR = 3.83, p < 0.001). High-risk patients exhibit pronounced immunosuppression, with enriched functions in synaptic vesicle cycling and signal transmission. Combined with single-cell analysis, it was observed that differentially expressed genes are widely distributed in both the pancreas and nerves, suggesting that pain-related MDGs not only affect the immune microenvironment but also regulate signal transmission between tumors and nerves. This is a key process in initiating visceral nociception. Our study finds that the complex pathogenesis of PC pain extends beyond mechanical compression and inflammatory mediator release, involving intricate immuno-neural genetic crosstalk (Figure 10).
FIGURE 10
Epigenetic modifications can affect gene expression at multiple levels, including transcription, splicing, stability, and translation, thereby influencing the physiological and pathological functions of cells. These modifications primarily consist of DNA methylation, histone modifications, non-coding RNAs, RNA modifications, and chromatin remodeling (Zheng et al., 2025). Among these, DNA methylation is the earliest and most extensively studied epigenetic modification. Research has found that DNA methylation can inhibit or activate genes related to pain perception, conduction, and regulation, thereby altering nerve cells’ response to pain (Sun et al., 2019). In a model of spinal nerve ligation (SNL) in male Sprague Dawley rats, a persistently low level of hypomethylation was observed at CpG sites in the DRG. Reducing DNA methylation can cause pain hypersensitivity, while increasing DNA methylation can alleviate neuropathic pain (). DNA methylation can also regulates neural plasticity, affecting connections and communication between nerve cells (Pratt et al., 2022). In oral cancer, abnormal demethylation of various tumor suppressor genes and hypermethylation of nerve growth factors have been found to be involved in the development of PNI (). In this study, MDGs were identified using the MethylMix package combined with gene expression profiling and methylation data analysis. These genes showed a negative correlation between methylation and mRNA expression, and differences in DNA methylation levels were observed between cancer and normal samples, resulting in 471 MDGs. Further screening yielded 26 MDGs related to PC pain, followed by GO and KEGG enrichment analyses. GO enrichment analysis revealed three main functions of these epigenetic regulators: (1) structural regulation through extracellular matrix organization, (2) transcriptional reprogramming via DNA-binding activator activity, and (3) neuronal excitability modulation through ligand-gated cation channels. KEGG pathway dissection uncovered the convergence of the immune-pain axis, with significant enrichment in T cell receptor signaling and Th17 differentiation pathways. Notably, our findings maybe suggest a novel mechanistic link: (1) epigenetic silencing of ion channels may disrupt pain signal transduction, (2) methylation-mediated immune modulation (Th17 pathway genes) indicates tumor microenvironment crosstalk, and (3) Methylation of structural genes (COL family) may alter the PNI ability of tumors. The multi-scale integration of methylation drivers with pain biology demonstrates epigenetic coordination between cancer progression and nociceptive signaling for the first time. Further analysis of the 26 genes led to the construction of a risk scoring model, identifying five key genes: PSMB8, COL17A1, BICC1, CTRC, and TRIP13. Among them, PSMB8, an immunoproteasome component, plays a crucial role in inflammation regulation (Yang et al., 2009; Kitamura et al., 2011). Abnormal methylation of PSMB8 has been reported in ovarian, breast, and colorectal cancers (Liew et al., 2018; Siebenkas et al., 2017; Tian et al., 2022). COL17A1, a member of the collagen family, exhibits abnormal promoter methylation leading to overexpression in cervical and epithelial cancers, enhancing tumor invasiveness (Thangavelu et al., 2016). It has also been identified as a potential biomarker for PC prognosis (). BICC1, an RNA-binding protein, influences tumor progression in colon cancer through pathways such as extracellular matrix (ECM) receptor interaction and focal adhesion (Lv et al., 2020). CTRC (chymotrypsin C) encodes a protein crucial for regulating the activation and degradation of trypsinogen and procarboxypeptidase, protecting the pancreas from pancreatitis and influencing pain progression in chronic pancreatitis (; ). TRIP13 encodes a thyroid hormone receptor-interacting protein, also known as a hormone-dependent transcription factor, promoting the progression of pancreatic ductal adenocarcinoma by facilitating tumor tissue growth and metastasis (). Based on the speculation of regulatory mechanisms combining five genetic features and above, it can be observed that a crucial genetic driver coordinates the regulation among three parties. PSMB8 affects neuroinflammation by regulating the activity of immunoproteasome, while COL17A1 and BICC1 promote the progression of PNI in tumor cells by regulating matrix homeostasis and ECM. The disruption of CTRC destabilizes pancreatic protease, exacerbating substantive nociception, and TRIP13 may influence related pain progression through hormonal regulation. The above research results indirectly suggest that epigenetic dysregulation in PC may be involved in the regulation of cancer pain by affecting immune reprogramming and the nervous system. However, further exploration is needed to determine whether there is a direct connection between epigenetics, immunology, and neurology. The specific internal relationships and molecular mechanisms also require more clinical and functional validation.
In the context of cancer pain, previous studies have predominantly focused on the regulatory mechanisms of pain within the nerve conduction pathway, overlooking the role of the tumor itself and the complex tumor-nerve crosstalk (Xu et al., 2024a). PC tumors consist of malignant tumor cells, stromal cells, immune cells, and other extracellular matrix components. Past research suggests that dynamic remodeling of the tumor microenvironment is critical to the progression of cancer pain. During malignant proliferation, mechanical compression of surrounding tissues by the tumor mass causes injury, and the release of various inflammatory factors (such as prostaglandins, interleukins, and tumor necrosis factors) can damage nerve endings, triggering inflammatory pain (). Metabolites of cancer cells, such as lactic acid, induce local acidosis and activate ASICs/TRPV1 ion channels (Qian et al., 2021). Additionally, tumor cells can directly invade nerves, causing structural damage (). More importantly, studies have confirmed that tumor-associated macrophages (TAMs) can activate Schwann cells through the bFGF/PI3K/Akt/c-myc/GFP pathway. Schwann cells secrete IL-33, recruiting macrophages into the perineural environment and promoting their M2 tumor-promoting polarization (Zhang et al., 2024). The positive feedback loop between the two has a significant impact on the PNI process in PC, indicating crosstalk between the immune system and nervous system in the development of cancer pain. Differential expression analysis of 22 immune cell types in this study showed that pain risk subtypes are associated with dendritic cell dysfunction and M0/M1 macrophage imbalance. Immune cell infiltration and GSVA enrichment scores indicated that MCODE differentially expressed genes are negatively correlated with various immunosuppressive cells such as nTregs, neutrophils, and DC cells, and positively correlated with immune effector cells such as CD8 T, Tfh, NK, CD4 T, and MAIT. Previous research and analysis results suggest that the tumor parenchymal microenvironment, rather than secondary neural effects, regulates PC pain progression through neuro-immune interactions. However, the epigenetic regulatory mechanisms involved still require further exploration.
By conducting a WGCNA on the differential genes of different pain risk subgroups, the cyan module containing 89 differential genes most relevant to the T staging was obtained. The results from KEGG suggest that these differential genes are involved in synapse vesicle recycling, insulin secretion, diabetes, and other pancreas-related functions. GO analysis revealed that the biological processes of these differential genes are primarily enriched in signal release, hormone secretion and transport, and synaptic plasticity regulation. Synaptic plasticity refers to the dynamic changes in synaptic transmission efficiency, including long-term potentiation (LTP) and long-term depression (LTD) (). In pain pathways (such as the spinal dorsal horn, thalamus, and cortex), LTP can lead to hyperalgesia and allodynia by enhancing synaptic transmission efficiency (Rygh et al., 2002). In chronic pain, glutamic acid triggers calcium influx through NMDA receptors, activating downstream signals (such as CaMKII, PKC), inducing synaptic LTP, and amplifying pain signals (). The cellular components and molecular functions of the cyan module are mainly enriched in vesicle transport, potassium ion channels, and related transporter activity regulation. Vesicle transport is responsible for the synthesis, storage, and release of neurotransmitters (such as glutamate, substance P, and CGRP) (Tao et al., 2022). In nerve damage or inflammation, dysregulation of vesicle transport-related proteins (such as SNARE complexes, synaptotagmin) leads to excessive release of glutamate and substance P, enhancing postsynaptic neuronal excitability (Pan and Rutecki, 2014). CGRP and substance P participate in neurogenic inflammation through vesicle release, further activating peripheral and central pain pathways (Meseguer et al., 2014). Botulinum toxin (BoNT) inhibits vesicle release by cleaving SNARE proteins and has been used to treat migraine and neuropathic pain (Vacca et al., 2020). Potassium ion channels affect the generation and transmission of pain signals by regulating membrane potential repolarization and neuronal excitability (Loucif et al., 2018). Thus, it can be seen that gene expression differences among different pain subtypes are mainly concentrated in various steps of pain-related signal transmission.It can be seen that the differentially expressed genes of pain risk subtypes affect various nodes of pain signal transmission through dynamic regulation of neurotransmitters and synaptic efficacy reprogramming. Since differentially expressed genes are also involved in insulin secretion and diabetes-related pathways, it suggests to some extent the organ-specific coupling of neuroendocrine-nociceptive transmission. The above studies indicate to some extent the potential connection between epigenetic regulation and tumor-neural crosstalk, and the specific mechanism is also the focus of further exploration.
Using cytoscope to screen MCODE genes from the cyan module, eight most significant pain-related differential genes were identified, namely, PTPRN2, CPE, SNAP25, ABCC8, PTPRN, CHGB, SCG5, and PCSK2. TIMER analysis revealed that the expression levels of these genes in PC are mostly negatively correlated with tumor purity, indicating high expression in the tumor microenvironment. However, PTPRN2 showed a positive correlation with tumor purity, suggesting relatively higher expression in the tumor. Therefore, the GSM6567159 PC immune landscape dataset and the GSM5236520 human spinal cord dataset were selected for single-cell analysis of MCODE genes. Among them, PTPRN2 encodes a protein tyrosine phosphatase receptor expressed in endocrine and neuronal cells (Sorokin et al., 2015), playing a role in exocytosis and affecting synaptic plasticity in neurons (). The single-cell analysis in this study found that PTPRN2 is distributed in pancreatic delta cells and mesenchymal cells, while in the spinal cord, it is evenly distributed in astrocytes, neurons, and stem cells. Immune cell function analysis suggested a positive correlation between PTPRN2 expression levels and various cells such as B cells, CD8+ T cells, and CD4+ T cells in PC. In high pain risk subtypes, PTPRN2 expression levels decrease, leading to reduced cellular immune function. CPE encodes an enzyme widely present in neuroendocrine cells, playing a key role in the processing of peptides and hormones. Studies suggest that CPE can promote the entry of eosinophil cationic protein into neuroendocrine cells (Wu et al., 2004). In the pancreas, CPE may be involved in the post-processing of hormones such as insulin and glucagon, closely related to pancreatic function (). Single-cell analysis of the pancreas showed that CPE is abundantly distributed in pancreatic Alpha and Beta cells, as well as in astrocytes, neurons, and stem cells. In high pain risk subtypes, CPE expression levels decrease, positively correlating with the function of various immune cells such as B cells, CD8+ T cells, and CD4+ T cells, leading to reduced cellular immune function. Studies have found that SNAP25 increases the release of excitatory neurotransmitters (such as glutamic acid) in the dorsal horn of the spinal cord through presynaptic mechanisms, potentially leading to an imbalance of excitatory/inhibitory neurotransmitters, thus mediating the development of neuropathic pain (Tafoya et al., 2006). In this study, SNAP25 is less distributed in PC, mainly in neuronal cells, positively correlating with immune cell function. The ABCC8 gene encodes the SUR1 subunit of the ATP-sensitive K+ channel (KATP) in pancreatic beta cells. Mutations in this channel can lead to congenital hyperinsulinism (CHI), a disease associated with excessive or unregulated insulin secretion. Although it does not directly cause pain, it may have indirect links to the physiological functions of the pancreas and related pain states (). Single-cell analysis showed less distribution of ABCC8 in pancreatic tissue and neurons. In high pain risk subtypes, ABCC8 expression levels decrease, positively correlating with the function of various immune cells such as B cells, CD8+ T cells, and CD4+ T cells, leading to reduced cellular immune function. PTPRN is a receptor-type tyrosine phosphatase primarily involved in regulating secretory granule control in neuroendocrine cells and islet beta cells, affecting the release of hormones (such as insulin) and neurotransmitters () Chromogranin B (CHGB) is a marker protein of secretory granules, involved in regulating the storage, processing, and release of neuropeptides and catecholamines () SCG5 is a key component of secretory granules, involved in regulating the processing and secretion of neuropeptides and hormonesPC () SK2 is responsible for cleaving inactive precursor proteins (such as neuropeptide precursors) into active forms (such as endorphin, proinsulin)ese (Podvin et al., 2018) four gene-encoded proteins all play a role in the production and transmission of neurotransmitters. Single-cell analysis revealed that PTPRN and SCG5 are distributed in pancreatic beta cells, while the expression of CHGB and PCSK2 is relatively low. In spinal cord tissue, PTPRN expression is low, while CHGB, SCG5, and PCSK2 are all enriched in neuronal cells. TIMER immune infiltration analysis showed that these genes regulate cellular immune function through the tumor microenvironment. Thus, it can be inferred that in the risk subgroups of pain-related MDGs in PC, differential gene expression primarily regulates neurotransmitter processing and transmission steps in neural signals through the tumor microenvironment, further confirming the mutual crosstalk between tumor and nerves.Immune function and single-cell resolution reveal that PC pain-related MCODE genes constitute neurotransmitter processing hubs that mediate immune-neural interactions through the tumor microenvironment. Among them, PTPRN2/CPE/SCG5/PCSK2 regulate immune function based on the neuroendocrine system, SNAP25/CHGB/PTPRN affect neural signaling based on synaptic signal transmission, and ABCC8 indirectly participates in pain regulation based on the physiological characteristics of ion channels. This study proposes for the first time at the single-cell level that tumor epigenetics hijack both neuroendocrine and immune systems by affecting the microenvironment, leading to the development and progression of cancer pain. It focuses on pancreatic parenchymal damage rather than nerve transmission as the key to PC pain, providing a new perspective for the mechanistic study of cancer pain.
The current clinical research progress has indirectly validated the “DNA methylation-immunity-neuro” interaction model speculated in this study to some extent. This is specifically manifested in the combined therapy of DNA methyltransferase inhibitors (DNMTi) and immune checkpoint inhibitors (ICI), as well as analgesic research on drugs targeting MCODE genes. Studies have found that DNMTi can reverse tumor immunosuppression by reactivating antitumor immune signals or reshaping the expression of immune checkpoints. In ovarian cancer, Peng et al. discovered that DNMTi enhances the therapeutic effect of anti-PD-L1 by reactivating Th1-type chemokines and increasing T-cell infiltration (Peng et al., 2015). La ure Ricard et al. treated nine patients with relapsed/refractory Angioimmunoblastic T-cell lymphoma (AITL) using a combination of 5-azacytidine (one of the DNMTi) and nivolumab (anti-PD-1 checkpoint blockers). The overall response rate reached 78% with reduced chemotherapy toxicity (Ricard et al., 2024). In breast cancer, combined therapy with PD-1 antibodies and guadecitabine (one of the DNMTi) can enhance MHC-I expression and increase CD8 T-cell infiltration in TME, thereby enhancing the therapeutic effect of PD-1 antibodies (Luo et al., 2018). In mouse models of melanoma, colorectal, breast, and ovarian tumors, the combined treatment of DNMTis and anti-CTLA-4 ICI can extend survival to some extent (). Cl inical drug research on MCODE genes has found that potassium channel openers such as retigabin and flupirtine are effective in treating neuropathic pain caused by nerve ligation models, inflammatory pain caused by formaldehyde, and visceral pain caused by capsaicin (). Perineural injection of botulinum toxin type A (pBONT-A) targeting SNAP25 under ultrasound guidance can alleviate pain associated with PC nerve infiltration (Meyer-Frießem et al., 2019).
Through bioinformatics analysis, this study proposes an innovative regulatory framework: There is a significant correlation between DNA methylation changes in PC pain, the progression of PNI, and the remodeling of the immune microenvironment. At the mechanistic level, it has been found that methylated driver genes are associated with pathological processes such as neural plasticity and immunosuppression, suggesting that epigenetic reprogramming may participate in the pathogenesis of pain through PNI-related nerve damage and immune remodeling. However, the causal relationship between the two still needs experimental verification. The constructed “epigenetic-immune-neural” three-dimensional network provides a testable hypothesis for subsequent mechanistic studies. However, there are several limitations to this study: Firstly, studies based on public datasets have insufficient sample diversity and lack multi-center validation. Secondly, the lack of in vivo models prevents the analysis of the spatiotemporal dynamics of pain signals. Finally, the discovery of correlation has not been verified through functional experiments, and the specific molecular mechanism cannot be clarified. Future research should focus on establishing an in situ PC model for neural epigenome editing, tracking the dynamic evolution of methylation during the PNI process, and determining the proportional weight of immune and neural contributions to the pain phenotype.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. The data presented in the study are deposited in the TCGA repository, MSigDB and GEO database, the accession number are respectively: https://portal.gdc.cancer.gov/projects/TCGA-PAAD, HP_PAIN (M38128), GSE183795.
Ethics statement
The studies involving humans were approved by Since all the data used were published previously in the public databases, no additional ethical approval and patient consent were needed. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
CZ: Conceptualization, Methodology, Formal Analysis, Data curation, Writing – original draft, Visualization. WX: Methodology, Conceptualization, Data curation, Validation, Formal Analysis, Writing – original draft, Visualization. XuC: Methodology, Validation, Supervision, Writing – review and editing. XW: Data curation, Methodology, Formal Analysis, Writing – review and editing. XiC: Writing – review and editing, Investigation, Funding acquisition, Supervision. WZ: Funding acquisition, Project administration, Conceptualization, Writing – review and editing, Supervision.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. The work was funded by grants (82102784) from the National Natural Science Foundation of China, the Project of Science and Technology of Nantong (MSZ2023056) and funded by Health Committee of Nantong (QN2024030, MS2023056, MS2023055, MS2024016).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Correction note
This article has been corrected with minor changes. These changes do not impact the scientific content of the article.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2025.1600883/full#supplementary-material
References
1
AfaqF.AgarwalS.BajpaiP.DiffalhaS. A.KimH. G.PeterS.et al (2024). Targeting of oncogenic AAA-ATPase TRIP13 reduces progression of pancreatic ductal adenocarcinoma. Neoplasia47, 100951. 10.1016/j.neo.2023.100951
2
AvganN.SutherlandH. G.LeaR. A.HauptL. M.ShumD. H. K.GriffithsL. R. (2023). Association study of a comprehensive panel of neuropeptide-related polymorphisms suggest potential roles in verbal learning and memory. Genes (Basel)15 (1), 30. 10.3390/genes15010030
3
BennettM. I.KaasaS.BarkeA.KorwisiB.RiefW.TreedeR. D.et al (2019). The IASP classification of chronic pain for ICD-11: chronic cancer-related pain. Pain160 (1), 38–44. 10.1097/j.pain.0000000000001363
4
BetheaJ. R.FischerR. (2021). Role of peripheral immune cells for development and recovery of chronic pain. Front. Immunol.12, 641588. 10.3389/fimmu.2021.641588
5
BlissT. V.CollingridgeG. L.KaangB. K.ZhuoM. (2016). Synaptic plasticity in the anterior cingulate cortex in acute and chronic pain. Nat. Rev. Neurosci.17 (8), 485–496. 10.1038/nrn.2016.68
6
BusserollesJ.TsantoulasC.EschalierA.López GarcíaJ. A. (2016). Potassium channels in neuropathic pain: advances, challenges, and emerging ideas. Pain157 (Suppl. 1), S7–S14. 10.1097/j.pain.0000000000000368
7
CaiW.LiuF.WenY.HanC.PrasadM.XiaQ.et al (2021). Pain management in acute pancreatitis: a systematic review and meta-analysis of randomised controlled trials. Front. Med. (Lausanne)8, 782151. 10.3389/fmed.2021.782151
8
CapodannoY.HirthM. (2023). Targeting the cancer-neuronal crosstalk in the pancreatic cancer microenvironment. Int. J. Mol. Sci.24 (19), 14989. 10.3390/ijms241914989
9
ChenO.DonnellyC. R.JiR. R. (2020). Regulation of pain by neuro-immune interactions between macrophages and nociceptor sensory neurons. Curr. Opin. Neurobiol.62 (17-25), 17–25. 10.1016/j.conb.2019.11.006
10
ChenY. C.TaylorA. J.FulcherJ. M.SwensenA. C.DaiX. Q.KombaM.et al (2023). Deletion of carboxypeptidase E in beta-cells disrupts proinsulin processing but does not lead to spontaneous development of diabetes in mice. Diabetes72 (9), 1277–1288. 10.2337/db22-0945
11
CostenlaA. R.LopesL. V.de MendoncaA.RibeiroJ. A. (2001). A functional role for adenosine A3 receptors: modulation of synaptic plasticity in the rat hippocampus. Neurosci. Lett.302 (1), 53–57. 10.1016/s0304-3940(01)01633-0
12
DammM.WenigerM.KolschA. K.LampertC.CeyhanG. O.BeerS.et al (2020). The quality of pain management in pancreatic cancer: a prospective multi-center study. Pancreatology20 (7), 1511–1518. 10.1016/j.pan.2020.08.017
13
DemcsakA.ShariatzadehS.Sahin-TothM. (2024). Secretagogue-induced pancreatitis in mice devoid of chymotrypsin. Am. J. Physiol. Gastrointest. Liver Physiol.327 (3), G333–G344. 10.1152/ajpgi.00310.2023
14
DrewesA. M.OlesenA. E.FarmerA. D.SzigethyE.ReboursV.OlesenS. S. (2020). Gastrointestinal pain. Nat. Rev. Dis. Prim.6(1), 1. 10.1038/s41572-019-0135-7
15
DunbarE. K.GreerP. J.SalomanJ. L.AlbersK. M.YadavD.WhitcombD. C. (2025). Genetics of constant and severe pain in the NAPS2 cohort of recurrent acute and chronic pancreatitis patients. J. Pain27, 104754. 10.1016/j.jpain.2024.104754
16
FengQ.SongD.WangX. (2021). Pan-cancer analysis reveals that neurotrophin signaling correlates positively with anti-tumor immunity, clinical outcomes, and response to targeted therapies and immunotherapies in cancer. Life Sci.282, 119848. 10.1016/j.lfs.2021.119848
17
FungM. M.NguyenC.MehtaniP.SalemR. M.PerezB.ThomasB.et al (2008). Genetic variation within adrenergic pathways determines in vivo effects of presynaptic stimulation in humans. Circulation117 (4), 517–525. 10.1161/CIRCULATIONAHA.107.706317
18
GarrigaJ.LaumetG.ChenS. R.ZhangY.MadzoJ.IssaJ. J.et al (2018). Nerve injury-induced chronic pain is associated with persistent DNA methylation reprogramming in dorsal root ganglion. J. Neurosci.38 (27), 6090–6101. 10.1523/JNEUROSCI.2616-17.2018
19
GilZ.CavelO.KellyK.BraderP.ReinA.GaoS. P.et al (2010). Paracrine regulation of pancreatic cancer cell invasion by peripheral nerves. J. Natl. Cancer Inst.102 (2), 107–118. 10.1093/jnci/djp456
20
GloynA. L.WeedonM. N.OwenK. R.TurnerM. J.KnightB. A.HitmanG.et al (2003). Large-scale association studies of variants in genes encoding the pancreatic beta-cell KATP channel subunits Kir6.2 (KCNJ11) and SUR1 (ABCC8) confirm that the KCNJ11 E23K variant is associated with type 2 diabetes. Diabetes52 (2), 568–572. 10.2337/diabetes.52.2.568
21
GuoR.LiJ.HuJ.FuQ.YanY.XuS.et al (2023). Combination of epidrugs with immune checkpoint inhibitors in cancer immunotherapy: from theory to therapy. Int. Immunopharmacol.120, 110417. 10.1016/j.intimp.2023.110417
22
HarashimaS.HoriuchiT.WangY.NotkinsA. L.SeinoY.InagakiN. (2012). Sorting nexin 19 regulates the number of dense core vesicles in pancreatic beta-cells. J. Diabetes Investig.3 (1), 52–61. 10.1111/j.2040-1124.2011.00138.x
23
HirthM.GandlaJ.HoperC.GaidaM. M.AgarwalN.SimonettiM.et al (2020). CXCL10 and CCL21 promote migration of pancreatic cancer cells toward sensory neurons and neural remodeling in tumors in mice, associated with pain in patients. Gastroenterology159 (2), 665–681. 10.1053/j.gastro.2020.04.037
24
HongS.ZhengG.WileyJ. W. (2015). Epigenetic regulation of genes that modulate chronic stress-induced visceral pain in the peripheral nervous system. Gastroenterology148 (1), 148–157. 10.1053/j.gastro.2014.09.032
25
HuangW. L.WuS. F.HuangX.ZhouS. (2022). Integrated analysis of ECT2 and COL17A1 as potential biomarkers for pancreatic cancer. Dis. Markers2022, 9453549. 10.1155/2022/9453549
26
HurnikP.ChyraZ.SevcikovaT.StembirekJ.TrtkovaK. S.GaykalovaD. A.et al (2022). Epigenetic regulations of perineural invasion in head and neck squamous cell carcinoma. Front. Genet.13, 848557. 10.3389/fgene.2022.848557
27
JiangS.TullochA. G.KimT. A.FuY.RogersR.GaskellA.et al (1998). Characterization and chromosomal localization of PTP-NP-2, a new isoform of protein tyrosine phosphatase-like receptor, expressed on synaptic boutons. Gene215 (2), 345–359. 10.1016/s0378-1119(98)00282-0
28
JiangW.TanX. Y.LiJ. M.YuP.DongM. (2022). DNA methylation: a target in neuropathic pain. Front. Med. (Lausanne)9, 879902. 10.3389/fmed.2022.879902
29
JiangL.CaiS.WengZ.ZhangS.JiangS. H. (2025). Peripheral, central, and chemotherapy-induced neuropathic changes in pancreatic cancer. Trends Neurosci.48 (2), 124–139. 10.1016/j.tins.2024.11.008
30
KalraA.MeltzerS. J. (2025). The role of DNA methylation in gastrointestinal disease: an expanded review of malignant and nonmalignant gastrointestinal diseases. Gastroenterology168 (2), 245–266. 10.1053/j.gastro.2024.07.001
31
KitamuraA.MaekawaY.UeharaH.IzumiK.KawachiI.NishizawaM.et al (2011). A mutation in the immunoproteasome subunit PSMB8 causes autoinflammation and lipodystrophy in humans. J. Clin. Invest121 (10), 4150–4160. 10.1172/JCI58414
32
LiT.LinC.WangW. (2025). Global, regional, and national burden of pancreatic cancer from 1990 to 2021, its attributable risk factors, and projections to 2050: a systematic analysis of the global burden of disease study 2021. BMC Cancer25(1), 189. 10.1186/s12885-025-13597-z
33
LiewP. L.HuangR. L.WengY. C.FangC. L.Hui-Ming HuangT.LaiH. C. (2018). Distinct methylation profile of mucinous ovarian carcinoma reveals susceptibility to proteasome inhibitors. Int. J. Cancer143 (2), 355–367. 10.1002/ijc.31324
34
LiuZ.LiM.ZhangL.ShiX.LiaoT.JieL.et al (2024). NGF signaling exacerbates KOA peripheral hyperalgesia via the increased TRPV1-labeled synovial sensory innervation in KOA rats. Pain Res. Manag.2024, 1552594. 10.1155/2024/1552594
35
LoucifA. J. C.SaintotP. P.LiuJ.AntonioB. M.ZellmerS. G.YogerK.et al (2018). GI-530159, a novel, selective, mechanosensitive two-pore-domain potassium (K(2P)) channel opener, reduces rat dorsal root ganglion neuron excitability. Br. J. Pharmacol.175 (12), 2272–2283. 10.1111/bph.14098
36
LuoN.NixonM. J.Gonzalez-EricssonP. I.SanchezV.OpalenikS. R.LiH.et al (2018). DNA methyltransferase inhibition upregulates MHC-I to potentiate cytotoxic T lymphocyte responses in breast cancer. Nat. Commun.9 (1), 248. 10.1038/s41467-017-02630-w
37
LvJ.WangJ.ShangX.LiuF.GuoS. (2020). Survival prediction in patients with colon adenocarcinoma via multi-omics data integration using a deep learning algorithm. Biosci. Rep.40 (12), BSR20201482. 10.1042/BSR20201482
38
MeseguerV.AlpizarY. A.LuisE.TajadaS.DenlingerB.FajardoO.et al (2014). TRPA1 channels mediate acute neurogenic inflammation and pain produced by bacterial endotoxins. Nat. Commun.5, 3125. 10.1038/ncomms4125
39
Meyer-FrießemC. H.EitnerL. B.KaislerM.MaierC.VollertJ.WestermannA.et al (2019). Perineural injection of botulinum toxin-A in painful peripheral nerve injury - a case series: pain relief, safety, sensory profile and sample size recommendation. Curr. Med. Res. Opin.35 (10), 1793–1803. 10.1080/03007995.2019.1626228
40
MiuraT.MitsunagaS.IkedaM.OhnoI.TakahashiH.KuwataT.et al (2018). Neural invasion spreads macrophage-related allodynia via neural root in pancreatic cancer. Anesth. Analg.126 (5), 1729–1738. 10.1213/ANE.0000000000002803
41
NiC.ChenL.HuaB.HanZ.XuL.ZhouQ.et al (2024). Epigenetic mechanisms of bone cancer pain. Neuropharmacology261, 110164. 10.1016/j.neuropharm.2024.110164
42
OzakiH.HiraokaT.MizumotoR.MatsunoS.MatsumotoY.NakayamaT.et al (1999). The prognostic significance of lymph node metastasis and intrapancreatic perineural invasion in pancreatic cancer after curative resection. Surg. Today29(1), 16–22. 10.1007/BF02482964
43
OzcanS.AkinS.Yildiz AltunA.BolatE.OnalS. A. (2019). Effective pain therapy can improve the survival of patients with pancreatic cancer. Agri31 (3), 145–149. 10.14744/agri.2019.61224
44
PanY. Z.RuteckiP. A. (2014). Enhanced excitatory synaptic network activity following transient group I metabotropic glutamate activation. Neuroscience275, 22–32. 10.1016/j.neuroscience.2014.05.062
45
PengD.KryczekI.NagarshethN.ZhaoL.WeiS.WangW.et al (2015). Epigenetic silencing of TH1-type chemokines shapes tumour immunity and immunotherapy. Nature527 (7577), 249–253. 10.1038/nature15520
46
PodvinS.WojniczA.HookV. (2018). Human brain gene expression profiles of the cathepsin V and cathepsin L cysteine proteases, with the PC1/3 and PC2 serine proteases, involved in neuropeptide production. Heliyon4 (7), e00673. 10.1016/j.heliyon.2018.e00673
47
PrattK. J. B.SheaJ. M.Remesal-GomezL.BieriG.SmithL. K.CouthouisJ.et al (2022). Loss of neuronal Tet2 enhances hippocampal-dependent cognitive function. Cell Rep.41 (6), 111612. 10.1016/j.celrep.2022.111612
48
QianH. Y.ZhouF.WuR.CaoX. J.ZhuT.YuanH. D.et al (2021). Metformin attenuates bone cancer pain by reducing TRPV1 and ASIC3 expression. Front. Pharmacol.12, 713944. 10.3389/fphar.2021.713944
49
RicardL.CerveraP.StockerN.CorreE.Van de WyngaertZ.BanetA.et al (2024). A combination of 5-azacytidine and nivolumab is a potentially effective rescue therapy in relapsed/refractory AITL. Front. Immunol.15, 1410638. 10.3389/fimmu.2024.1410638
50
RyghL. J.TjolsenA.HoleK.SvendsenF. (2002). Cellular memory in spinal nociceptive circuitry. Scand. J. Psychol.43 (2), 153–159. 10.1111/1467-9450.00281
51
SelvaggiF.MelchiorreE.CasariI.CinalliS.CinalliM.AcetoG. M.et al (2022). Perineural invasion in pancreatic ductal adenocarcinoma: from molecules towards drugs of clinical relevance. Cancers (Basel)14 (23), 5793. 10.3390/cancers14235793
52
ShresthaS.SapkotaS.TeohS. L.KcB.PaudyalV.LeeS. W. H.et al (2024). Comprehensive assessment of pain characteristics, quality of life, and pain management in cancer patients: a multi-center cross-sectional study. Qual. Life Res.33 (10), 2755–2771. 10.1007/s11136-024-03725-w
53
SiebenkasC.ChiappinelliK. B.GuzzettaA. A.SharmaA.JeschkeJ.VatapalliR.et al (2017). Inhibiting DNA methylation activates cancer testis antigens and expression of the antigen processing and presentation machinery in colon and ovarian cancer cells. PLoS One12 (6), e0179501. 10.1371/journal.pone.0179501
54
SorokinA. V.NairB. C.WeiY.AzizK. E.EvdokimovaV.HungM. C.et al (2015). Aberrant expression of proPTPRN2 in cancer cells confers resistance to apoptosis. Cancer Res.75 (9), 1846–1858. 10.1158/0008-5472.CAN-14-2718
55
SunL.ZhaoJ. Y.GuX.LiangL.WuS.MoK.et al (2017). Nerve injury-induced epigenetic silencing of opioid receptors controlled by DNMT3a in primary afferent neurons. Pain158 (6), 1153–1165. 10.1097/j.pain.0000000000000894
56
SunL.GuX.PanZ.GuoX.LiuJ.AtianjohF. E.et al (2019). Contribution of DNMT1 to neuropathic pain genesis partially through epigenetically repressing Kcna2 in primary afferent neurons. J. Neurosci.39 (33), 6595–6607. 10.1523/JNEUROSCI.0695-19.2019
57
TafoyaL. C.MameliM.MiyashitaT.GuzowskiJ. F.ValenzuelaC. F.WilsonM. C. (2006). Expression and function of SNAP-25 as a universal SNARE component in GABAergic neurons. J. Neurosci.26 (30), 7826–7838. 10.1523/JNEUROSCI.1866-06.2006
58
TaoT.LiuY.ZhangJ.LaiW.LongH. (2022). NGF-induced upregulation of CGRP in orofacial pain induced by tooth movement is dependent on Atp6v0a1 and vesicle release. Int. J. Mol. Sci.23 (19), 11440. 10.3390/ijms231911440
59
TarasiukA.MirochaG.FichnaJ. (2023). Current status of complementary and alternative medicine interventions in the management of pancreatic cancer - an overview. Curr. Treat. Options Oncol.24 (12), 1852–1869. 10.1007/s11864-023-01146-4
60
Ten BargeJ. A.van den BoschG. E.SlaterR.van den HoogenN. J.ReissI. K. M.SimonsS. H. P. (2025). Visceral pain in preterm infants with necrotizing enterocolitis: underlying mechanisms and implications for treatment. Paediatr. Drugs27 (2), 201–220. 10.1007/s40272-024-00676-0
61
ThangaveluP. U.KrenacsT.DrayE.DuijfP. H. (2016). In epithelial cancers, aberrant COL17A1 promoter methylation predicts its misexpression and increased invasion. Clin. Epigenetics8, 120. 10.1186/s13148-016-0290-6
62
TianT.FuJ.LiD.LiuY.SunH.WangX.et al (2022). Methylation of immune-related genes in peripheral blood leukocytes and breast cancer. Front. Oncol.12, 817565. 10.3389/fonc.2022.817565
63
VaccaV.MadaroL.De AngelisF.ProiettiD.CobianchiS.OrsiniT.et al (2020). Revealing the therapeutic potential of botulinum neurotoxin type A in counteracting paralysis and neuropathic pain in spinally injured mice. Toxins (Basel)12 (8), 491. 10.3390/toxins12080491
64
WuC. M.ChangH. T.ChangM. D. (2004). Membrane-bound carboxypeptidase E facilitates the entry of eosinophil cationic protein into neuroendocrine cells. Biochem. J.382 (Pt 3), 841–848. 10.1042/BJ20040894
65
XuL.YangL.WuY.WanX.TangX.XuY.et al (2023). Rac1/PAK1 signaling contributes to bone cancer pain by regulation dendritic spine remodeling in rats. Mol. Pain19, 17448069231161031. 10.1177/17448069231161031
66
XuW.LiuJ.ZhangJ.LuJ.GuoJ. (2024a). Tumor microenvironment crosstalk between tumors and the nervous system in pancreatic cancer: molecular mechanisms and clinical perspectives. Biochim. Biophys. Acta Rev. Cancer1879 (1), 189032. 10.1016/j.bbcan.2023.189032
67
XuJ.ShangY.WangT.SongJ.ZhuW.ZengY.et al (2024b). Nerve growth factor from pancreatic cancer cells promotes the cancer progression by inducing nerve cell-secreted interleukin-6. J. Interferon Cytokine Res.44 (12), 541–549. 10.1089/jir.2024.0154
68
XueM.ZhuY.JiangY.HanL.ShiM.SuR.et al (2023). Schwann cells regulate tumor cells and cancer-associated fibroblasts in the pancreatic ductal adenocarcinoma microenvironment. Nat. Commun.14 (1), 4600. 10.1038/s41467-023-40314-w
69
YangZ.GagarinD.St LaurentG.3rdHammellN.TomaI.HuC. A.et al (2009). Cardiovascular inflammation and lesion cell apoptosis: a novel connection via the interferon-inducible immunoproteasome. Arterioscler. Thromb. Vasc. Biol.29 (8), 1213–1219. 10.1161/ATVBAHA.109.189407
70
YangY.YangW.ZhangR.WangY. (2024). Peripheral mechanism of cancer-induced bone pain. Neurosci. Bull.40 (6), 815–830. 10.1007/s12264-023-01126-6
71
YuanH.DuS.ChenL.XuX.WangY.JiF. (2020). Hypomethylation of nerve growth factor (NGF) promotes binding of C/EBPα and contributes to inflammatory hyperalgesia in rats. J. Neuroinflammation17 (1), 34. 10.1186/s12974-020-1711-1
72
ZengD.YeZ.ShenR.YuG.WuJ.XiongY.et al (2021). IOBR: multi-omics immuno-oncology biological research to decode tumor microenvironment and signatures. Front. Immunol.12, 687975. 10.3389/fimmu.2021.687975
73
ZhangK.WangH. (2015). Cancer genome atlas pan-cancer analysis Project. Zhongguo Fei Ai Za Zhi18 (4), 219–223. 10.3779/j.issn.1009-3419.2015.04.02
74
ZhangB.GuoX.HuangL.ZhangY.LiZ.SuD.et al (2024). Tumour-associated macrophages and Schwann cells promote perineural invasion via paracrine loop in pancreatic ductal adenocarcinoma. Br. J. Cancer130(4), 542–554. 10.1038/s41416-023-02539-w
75
ZhaoJ.HuhY.BortsovA.DiatchenkoL.JiR. R. (2023). Immunotherapies in chronic pain through modulation of neuroimmune interactions. Pharmacol. Ther.248, 108476. 10.1016/j.pharmthera.2023.108476
76
ZhengH.HuaY.YangS.LiuV.HuangN.LiJ.et al (2025). Epigenetic modification and tumor immunity: unraveling the interplay with the tumor microenvironment and its therapeutic vulnerability and implications. Cancer Lett.616, 217587. 10.1016/j.canlet.2025.217587
77
ZhuM.LuoF.XuB.XuJ. (2024). Research progress of neural invasion in pancreatic cancer. Curr. Cancer Drug Targets24 (4), 397–410. 10.2174/1568009623666230817105221
Summary
Keywords
pancreatic carcinoma, DNA methylation, perineural invasion (PNI), neuropathic pain, immune microenvironment, prognostic signature
Citation
Zhang C, Xu W, Chen X, Wang X, Cui X and Zhao W (2025) Pain - related methylation driver genes affect the prognosis of pancreatic cancer patients by altering immune function and perineural infiltration. Front. Genet. 16:1600883. doi: 10.3389/fgene.2025.1600883
Received
03 April 2025
Accepted
15 August 2025
Published
08 October 2025
Corrected
29 October 2025
Volume
16 - 2025
Edited by
Ragini Singh, National Institutes of Health (NIH), United States
Reviewed by
Boran Zhu, Nanjing University of Chinese Medicine, China
Manjit Kumar Srivastav, National Cancer Institute (NIH), United States
Updates
Copyright
© 2025 Zhang, Xu, Chen, Wang, Cui and Zhao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wenjing Zhao, zhaowenjing@ntu.edu.cn
† These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.