Skip to main content

Comprehensive analysis identifies cuproptosis-related gene DLAT as a potential prognostic and immunological biomarker in pancreatic adenocarcinoma

Abstract

Background

Cuproptosis is a regulated cell death form associated with tumor progression, clinical outcomes, and immune response. However, the role of cuproptosis in pancreatic adenocarcinoma (PAAD) remains unclear. This study aims to investigate the implications of cuproptosis-related genes (CRGs) in PAAD by integrated bioinformatic methods and clinical validation.

Methods

Gene expression data and clinical information were downloaded from UCSC Xena platform. We analyzed the expression, mutation, methylation, and correlations of CRGs in PAAD. Then, based on the expression profiles of CRGs, patients were divided into 3 groups by consensus clustering algorithm. Dihydrolipoamide acetyltransferase (DLAT) was chosen for further exploration, including prognostic analysis, co-expression analysis, functional enrichment analysis, and immune landscape analysis. The DLAT-based risk model was established by Cox and LASSO regression analysis in the training cohort, and then verified in the validation cohort. Quantitative reverse transcriptase polymerase chain reaction (RT-qPCR) and immunohistochemistry (IHC) assays were performed to examine the expression levels of DLAT in vitro and in vivo, respectively.

Results

Most CRGs were highly expressed in PAAD. Among these genes, increased DLAT could serve as an independent risk factor for survival. Co-expression network and functional enrichment analysis indicated that DLAT was engaged in multiple tumor-related pathways. Moreover, DLAT expression was positively correlated with diverse immunological characteristics, such as immune cell infiltration, cancer-immunity cycle, immunotherapy-predicted pathways, and inhibitory immune checkpoints. Submap analysis demonstrated that DLAT-high patients were more responsive to immunotherapeutic agents. Notably, the DLAT-based risk score model possessed high accuracy in predicting prognosis. Finally, the upregulated expression of DLAT was verified by RT-qPCR and IHC assays.

Conclusions

We developed a DLAT-based model to predict patients’ clinical outcomes and demonstrated that DLAT was a promising prognostic and immunological biomarker in PAAD, thereby providing a new possibility for tumor therapy.

Peer Review reports

Introduction

Pancreatic cancer is a highly malignant tumor with extremely poor prognosis. In 2023, approximately 64,050 new cases will be detected, leading to estimated 50,550 deaths [1]. Pancreatic adenocarcinoma (PAAD) is the most common histological subtype of pancreatic cancer, accounting for approximately more than 85% of all cases [2]. The prevalence of PAAD is relevant to smoking, obesity, diabetes, alcohol use, chronic pancreatitis, and genetic factors [3]. Patients with PAAD are commonly unable to receive early diagnosis and treatments due to their non-specific symptoms and aggressive behaviors [4]. They generally result in poor outcomes with a 5-year survival rate as low as 9% [5]. At initial diagnosis, more than half of the patients suffer from metastatic diseases, exhibiting an average survival time of less than one year [6]. Multiple therapeutic approaches, such as surgery, chemotherapy, radiotherapy, targeted therapy, and immunotherapy have been developed to increase patients’ survival rates [7]. Nevertheless, PAAD continues to threaten human life and health seriously. Therefore, it is imperative to detect novel prognostic biomarkers to overcome the clinical dilemma in PAAD treatment.

Immunotherapy has provided more therapeutic opinions for cancer patients in recent years. Immune checkpoints inhibitors that target programmed cell death-1 (PD-1), programmed cell death-ligand 1 (PD-L1), and cytotoxic T-lymphocyte antigen-4 (CTLA-4) have been approved for the treatment of PAAD [8, 9]. Unfortunately, the immunosuppressive tumor microenvironment (TME) in PAAD that consists of dense fibrotic stromal, suppressed effector T cells, activated myeloid-derived suppressor cells, and insufficient secreted cytokines strongly limits the efficacy of immunotherapy [10]. Thus, it is crucial to identify useful immunological indicators to promote individualized immunotherapy in PAAD.

Cuproptosis, a newly defined cell death program, has attracted much attention due to its unique mechanism [11]. Unlike other cell death forms, such as apoptosis, autophagy, pyroptosis, and ferroptosis, cuproptosis is a copper-dependent pathway related to the tricarboxylic acid (TCA) cycle in mitochondria [12]. As an essential trace metal, copper is vital in regulating enzyme functions and body homeostasis. However, copper accumulation will induce cell death [13]. Tsvetkov et al. elucidated this lethal mechanism that excessive copper ions could directly bind to the lipoylated mitochondrial proteins, leading to lipoylated proteins aggregation and iron-sulfur cluster proteins reduction, eventually causing cell death [11]. Most researchers made explorations on the role of cuproptosis in tumors, and discovered that cuproptosis-related genes (CRGs) could function as prognostic biomarkers with the potential to guide immunotherapy in hepatocellular carcinoma [14], kidney cancer [15], breast cancer [16], and other tumors [17,18,19]. It was considered that cuproptosis might open the area of tumor treatment and provide new therapeutic targets in the future. Thus, we performed this analysis to explore the prognostic and therapeutic implications of cuproptosis in PAAD.

Herein, we conducted a comprehensive analysis to assess the function of CRGs in PAAD, including expression profiles, mutation and methylation status, correlation analysis, and consensus clustering analysis. We chose DLAT for further exploration and performed survival analysis to determine its prognostic value. To further investigate its biological function, we established the protein–protein interaction (PPI), gene–gene interaction, and co-expression network, then identified differentially expressed genes (DEGs) for enrichment analysis. We also uncovered the positive correlations between DLAT and immunological characteristics. In summary, our research indicated that DLAT could be a potential prognostic and immunological biomarker in PAAD.

Materials and methods

Data collection

Thirteen CRGs were obtained from the research of Tsvetkov et al., including 7 positive regulatory genes (FDX1, LIPT1, LIAS, DLD, DLAT, PDHA1, and PDHB), 3 negative regulatory genes (MTF1, GLS, and CDKN2A), and 3 copper transport-related genes (SLC31A1, ATP7A and ATP7B) [11]. The gene expression data of PAAD samples collected in the Cancer Genome Atlas (TCGA) database and normal samples saved in the genotype-tissue expression (GTEx) database were obtained from the UCSC Xena platform (http://xena.ucsc.edu/). The fragments per kilobase million (FPKM) sequencing data were transformed to transcripts per million (TPM) data and normalized by log2 (TPM + 1). Then, the TPM data of TCGA and GTEx databases were integrated using the "normalize between array" function of "limma" package (v3.52.4). To verify the differential expression of DLAT between tumor and normal tissues, GSE62452, GSE71729, GSE15471, and GSE16515 datasets were downloaded from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/).

Mutation and methylation analysis

Copy number variation (CNV) data were collected using the UCSC Xena platform. The location of CNV changes on chromosomes was visualized through "RCircos" package (v1.2.2). The promoter methylation levels of CRGs were analyzed by the UALCAN database (http://ualcan.path.uab.edu/) [20].

Consensus clustering analysis

Patients in TCGA cohort were clustered using "ConsensusClusterPlus" package (v1.60.0) according to their expression levels of CRGs. In the consensus clustering analysis, the maximum number of clusters was defined as 6, and 80% of the total samples were drawn 1000 times, clusterAlg = "pam", innerLinkage = "ward.D". The optimal number of clusters was determined by consensus matrix heat map, consensus cumulative distribution function, and delta area plot together. Following that, we compared the age, gender, histologic grade, clinical stage, and survival status among different clusters. Survival differences were analyzed based on Kaplan–Meier curves using "survminer" package (v0.4.9). The "DESeq2" package (v1.36.0) was employed to identify DEGs between cluster A and cluster C on the basis of count data, with the threshold set at fold change > 2 and adjusted p < 0.05.

Functional enrichment analysis

The "ClusterProfiler" package (v4.4.4) was utilized to explore the function of potential targets, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways [21], and Gene Set Enrichment Analysis (GSEA). The "HALLMARK" gene set was downloaded from the Molecular Signatures Database. The "GSVA" package (v1.44.5) was used to perform Gene Set Variation Analysis (GSVA) for analyzing pathways activities between different groups. The volcano map, box plot, and bubble plot were completed by "ggplot2" package (v3.3.6). The "pheatmap" package (v1.0.12) was used to draw the heatmap.

Survival analysis

To assess the overall survival (OS) differences, we conducted Kaplan–Meier analysis based on log-rank test through "survival" package (v3.4–0) and "survminer" package (v0.4.9). The comparison of objective response rates between DLAT-high and DLAT-low groups was carried out by chi-square test. We then evaluated hazard ratio (HR) with 95% confidence interval (CI) using univariate and multivariate Cox proportional hazards regression. Forest plot was performed by "forestplot" package (v.3.0.0) to visualize p value, HR, and 95% CI of each variable.

Construction of PPI, gene–gene interaction, and co-expression network

The STRING database (https://string-db.org/) was used to create the PPI network [22]. The GeneMANIA database (http://genemania.org/) was utilized to establish the gene–gene interaction network [23]. Pearson’s correlation analysis was applied to evaluate the correlations between DLAT and other genes for co-expression network construction.

Immunological characteristics evaluation

Immunological characteristics were composed of immune cell infiltration, TME score, immunotherapy-predicted pathways, cancer-immunity cycle, and inhibitory immune checkpoints. We first analyzed the correlation between DLAT CNV and immune cell infiltration in TIMER database (https://cistrome.shinyapps.io/timer/) [24]. Immune infiltrating cells included B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils, and dendritic cells. Following that, we evaluated the relationship between DLAT expression and the abundance of immune cells in the TIMER database. We also assessed the association between DLAT expression and the infiltration levels of 24 immune cells, which were generated by single sample Gene Set Enrichment Analysis (ssGSEA). The "IOBR" package (v0.99.9) was used to calculate the TME score [25].

The relationships between DLAT and immunotherapy-predicted pathways, as well as cancer-immunity cycle were also investigated. The 18 immunotherapy-predicted pathways included IFN-gamma signature, APM signal, base excision repair, cell cycle, DNA replication, Fanconi anemia pathway, and other pathways. The signature genes were acquired from previous studies [26, 27]. Moreover, the 7 steps of cancer-immunity cycle were described in 2013, including release of cancer antigens, cancer antigen presentation, priming and activation, trafficking of T cells to tumors, infiltration of T cells into tumors, recognition of cancer cells by T cells, and killing of cancer cells [28]. The signature gene set was obtained from the TIP website (http://biocc.hrbmu.edu.cn/TIP/) [29]. All the correlations were calculated through Pearson’s correlation analysis.

We integrated 5 frequently used immunological algorithms, including TIMER [30], CIBERSORT [31], MCPCOUNTER [32], QUANTISEQ [33], and EPIC [34] to obtain the immune cell infiltration matrix. The heatmap was carried out by "pheatmap" package (v1.0.12) to visualize the distribution of tumor-infiltrating immune cells.

Furthermore, we assessed the associations between DLAT and common inhibitory immune checkpoints, such as PD-L1, PD-L2, CTLA-4, TIM-3, TIGIT and so on. To investigate the role of DLAT in guiding immunotherapy, we analyzed the response to immune checkpoint blockade (ICB) by comparing our data with a meaningful melanoma dataset of 47 patients who received anti-PD-1 and anti-CTLA-4 treatments by subclass mapping algorithm [35,36,37]. Moreover, IMvigor210 [27], CheckMate [38], GSE78220 [39], and GSE91061 [40] cohorts were also used to predict the clinical prognosis of patients under immunotherapy.

Construction of risk score model

The "DESeq2" package (v1.36.0) was used to determine the DEGs (fold change = 1.5, adjusted p < 0.05) between DLAT-high and DLAT-low groups. Next, we conducted the univariable Cox regression analysis to obtain the DEGs with prognostic value. The risk score model was then developed on the basis of the prognostic DEGs. The "survival" package (v3.4–0) and "glmnet" package (v4.1–4) were employed to perform Cox and LASSO regression analysis, respectively. All patients in the TCGA-PAAD cohort were included in the training set. They were classified into low- or high-risk groups using the median risk score. Survival differences were evaluated by log-rank test. The accuracy of model was estimated by time-dependent receiver operating characteristic (ROC) curves with corresponding areas under the curve (AUC) values using "survivalROC" package (v1.0.3).

Cell culture

Human pancreatic cancer cell lines (BxPC-3 and PANC-1) and normal pancreatic cell line (HPDE6-C7) were obtained from the Chinese Academy of Sciences (Shanghai, China) and cultured with complete growth medium, as recommended by the manufacturer. All cells were incubated at 37℃ in a 5% of CO2 environment.

Quantitative reverse transcriptase polymerase chain reaction (RT-qPCR)

Total RNA was extracted of cells using TRIzol reagent. Next, the RNA was reverse transcribed into cDNA using a PrimeScript RT Reagent Kit with gDNA Eraser (Takara, Japan) according to the manufacturer’s instructions. RT-qPCR was conducted using TB Green Premix Ex Taq II (Takara, Japan). The reaction conditions were set as follows: 95℃ for 30 s, 40 cycles of 95 ℃ for 3 s and 60 ℃ for 30 s. Relative mRNA levels were determined by the 2−∆∆Ct method using GAPDH mRNA as an internal control. The primer sequences were shown as follows. DLAT-forward: ACTCCCCAGCCTTTAGCTC, DLAT-reverse: CAATCCCTTTCTCTACTGCCAAC, GAPDH forward: GCTGAGAACGGGAAGCTTGT, GAPDH reverse: GCCAGGGGTGCTAAGCAGTT.

Immunohistochemistry (IHC)

Human pancreatic cancer tissues were collected from 12 patients in Harbin Medical University Cancer Hospital. Patients have not received any anti-cancer treatment before tissue collection. Tissues were fixed in formalin, embedded in paraffin, and cut into 3 μm sections. The 2-step plus Poly-HRP Anti-Rabbit IgG detection system with DAB solution (Elabscience, E-IR-R217) was used to perform IHC staining. Sections were incubated with primary antibody against DLAT (1:500, Proteintech, 13,426-1AP) overnight at 4℃, followed by secondary antibody polyperoxidase-anti-rabbit IgG for 30 min at 37℃, and finally stained by DAB reagent.

Statistical analysis

The comparison between 2 groups was analyzed by Wilcoxon test, and the comparison among 3 groups was evaluated by Kruskal–Wallis test. Chi-square test was used for the comparison of objective response rate. Kaplan–Meier analysis was utilized to assess the survival differences through log-rank test. Pearson’s correlation analysis was employed to examine the correlation between 2 variables. All statistical analyses were performed using R (v4.2.1). p < 0.05 was considered statistically significant.

Results

Expression profiles and mutation status of CRGs in PAAD

Tsvetkov et al. have defined 13 CRGs which were expected to serve as novel targets in cancer treatment. We analyzed their expression levels in PAAD tissues and normal tissues, and found that all genes were differently expressed (Fig. 1A-B). In detail, FDX1, LIAS, LIPT1, DLD, DLAT, PDHB, MTF1, GLS, CDKN2A, SLC31A1, ATP7A, and ATP7B were highly expressed in tumor tissues, while PDHA1 was upregulated in normal tissues. Moreover, we examined the CNV alterations and revealed that GLS had the most significant copy number gain, whereas CDKN2A showed the most significant copy number loss (Fig. 1C). The location of CNV alterations on chromosomes was visualized in Fig. 1D. Promoter methylation levels of FDX1, LIAS, GLS, CDKN2A, and ATP7B were upregulated in tumor tissues, whereas the level of DLD was upregulated in healthy tissues (Fig. 1E). In addition, most CRGs were positively correlated with each other (Fig. 1F). The log-rank analysis depicted 5 prognostic genes, of which FDX1, DLAT, and ATP7A were associated with poor outcomes, while LIAS and LIPT1 were related to favorable outcomes (Fig. 1F).

Fig. 1
figure 1

Characteristics of CRGs in PAAD. A Heat map and B box plot of 13 CRGs in tumor tissues and normal tissues. C CNV frequency of CRGs in PAAD. D Location of CNV alterations on chromosomes. E Promoter methylation levels of CRGs in tumor tissues and normal tissues. F Correlation analysis and prognostic values of CRGs in PAAD. Positive correlation, red; negative correlation, green (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001)

Identification and biological characteristics of CRG clusters in PAAD

Based on the expression of CRGs, patients were clustered using consensus clustering algorithm. The results showed that consensus matrix k = 3 can produce the best result (Fig. 2A-C). Patients were classified into 3 subgroups, termed cluster A (n = 37), cluster B (n = 80), and cluster C (n = 54). The associations between 3 clusters and clinical features were presented in Fig. 2D. Five genes (DLAT, PDHA1, GLS, CDKN2A, and ATP7B) were differentially expressed among clusters. PDHA1 and ATP7B were increased in cluster A, CDKN2A was upregulated in cluster B, while DLAT and GLS were highly expressed in cluster C (Fig. 2E). Survival analysis indicated a significant difference on OS (log-rank test, p = 0.042) (Fig. 2F). Patients in cluster A exhibited the best median OS of 44.4 months, whereas those in cluster C showed the worst OS of 17.2 months, implying that patients in cluster C had higher degree of malignancy than those in cluster A.

Fig. 2
figure 2

Consensus clustering analysis of CRGs in PAAD. A-C Consensus matrix heatmap defining 3 clusters and their correlation area. D Associations between 3 clusters and clinical features. E Box plot of 5 differentially expressed CRGs among 3 clusters. F Kaplan–Meier curves for patients in different subgroups. G Volcano plot and H heatmap of DEGs between cluster A and C. I GO, J KEGG, and K GSVA of upregulated DEGs between cluster A and C (*p < 0.05; **p < 0.01; ****p < 0.0001)

We identified 224 upregulated genes (including ITLN1, CLCA1, SPINK4, NPHS2, GAST, and MUC2) and 194 downregulated genes between cluster A and cluster C (Fig. 2G-H). To explore their function, we conducted the enrichment analysis based on upregulated DEGs. The results of GO, KEGG, and GSVA together uncovered that upregulated DEGs were closely associated with tumor metabolic regulation, such as olefinic compound metabolic process, retinol metabolism, and xenobiotic metabolism (Fig. 2I-K).

Expression profiles and prognostic values of DLAT in PAAD

In this research, DLAT possessed the lowest p value (p < 0.001) in prognostic analysis, indicating that DLAT was strongly associated with clinical outcomes in PAAD. (Fig. 1F). Moreover, we found that the survival time was increasingly shortened in cluster A, B, and C (Fig. 2F). Meanwhile, the expression levels of DLAT were getting higher (Fig. 2E). This result further suggested that high expression of DLAT was correlated with poor prognosis. Thus, we selected DLAT for further analysis.

We downloaded the external datasets GSE62452, GSE71729, GSE15471, and GSE16515 to confirm the differential expression of DLAT. In accordance with previous results, DLAT was significantly increased in tumor tissues than normal tissues (Fig. 3A). To further examine its prognostic signification, patients in TCGA cohort with the survival time of more than 30 days were separated into DLAT-high (n = 71) and DLAT-low (n = 100) groups. We found that DLAT-high patients exhibited a shorter OS (log-rank test, p < 0.001) (Fig. 3B) and a lower objective response rate (chi-square test, p = 0.043) (Fig. 3C) than DLAT-low patients. Moreover, increased DLAT expression was considered as an independent risk factor for survival (univariable Cox, HR = 2.100, p < 0.001; multivariable Cox, HR = 2.287, p < 0.001) (Fig. 3D-E).

Fig. 3
figure 3

Expression profiles and prognostic values of DLAT in PAAD. A Box plot of DLAT expression between tumor and normal tissues in GSE62452, GSE71729, GSE15471, and GSE16515 datasets. B Kaplan–Meier curves for patients in DLAT-high and DLAT-low groups. C Comparison of clinical response rates (CR, complete response; PR, partial response; SD, stable disease; PD, progressive disease) in DLAT-high and DLAT-low groups. D Univariate and E multivariate Cox analysis of DLAT expression and clinical factors (*p < 0.05; **p < 0.01; ***p < 0.001)

Construction of PPI, gene–gene interaction, and co-expression network

In order to analyze the biological function of DLAT, we first established the PPI network in the STRING database (Fig. 4A). Then, the DLAT-associated gene–gene interaction network was built using the GeneMANIA database (Fig. 4B). It was observed that the CRGs, including PDHB, PDHA1, and DLD were involved in the PPI and gene–gene interaction network. Additionally, we performed correlation analysis to construct DLAT co-expression network. A total of 5474 genes were found to be positively related to DLAT (Fig. 4C) (Supplementary Table 1). Among them, 6 genes exhibited the strongest correlation, namely RAB6A (r = 0.73, p < 0.001), ARCN1 (r = 0.70, p < 0.001), PRKDC (r = 0.67, p < 0.001), RBM7 (r = 0.67, p < 0.001), ARMC1 (r = 0.66, p < 0.001), and ZW10 (r = 0.66, p < 0.001) (Fig. 4D). GSEA of DLAT co-expression network indicated the top 6 enriched hallmark terms, including G2/M checkpoint (NES = 3.255, p < 0.001), E2F targets (NES = 3.153, p < 0.001), mTORC1 signaling (NES = 2.948, p < 0.001), MYC targets (NES = 2.943, p < 0.001), epithelial-mesenchymal transition (NES = 2.908, p < 0.001), and protein secretion (NES = 2.903, p < 0.001), which were associated with cell proliferation and tumor metastasis (Fig. 4E). The detailed description of GSEA was shown in Supplementary Table 2.

Fig. 4
figure 4

Integrated analysis of PPI, gene–gene interaction, and co-expression network. A Construction of PPI network. B Construction of gene–gene interaction network. C Volcano map of co-expression genes associated with DLAT. D Correlation analysis of DLAT expression with RAB6A, ARCN1, PRKDC, RBM7, ARMC1, and ZW10. E Hallmark terms enriched in G2/M checkpoint, E2F targets, mTORC1 signaling, MYC targets, epithelial-mesenchymal transition, and protein secretion based on GSEA

Identification of DEGs and enrichment analysis

We identified 144 upregulated genes (such as LCT, MUC21, PAX7, and C6orf15) and 680 downregulated genes (such as C6orf58, GAST, DBET, FGL1, and OXT) in DLAT-high group when compared with DLAT-low group (Fig. 5A) (Supplementary Table 3). The top 200 DEGs were shown in Fig. 5B. We then carried out enrichment analysis to clarify the function of DLAT. GO analysis consisted of biological process, cell component, and molecular function. The results showed that DEGs were involved in nuclear pore organization, condensed chromosome outer kinetochore, and annealing activity (Fig. 5C-E). KEGG analysis demonstrated that DEGs were primarily engaged in cell proliferation and DNA repair, such as mismatch repair, DNA replication, and cell cycle (Fig. 5F). The following GSVA revealed the enrichments in G2/M checkpoint, E2F targets, mTORC1 signaling, TGF-β signaling, and complement, which were related to cell proliferation and immune response (Fig. 5G). In conclusion, DLAT may promote the tumor progression through these tumor-related pathways.

Fig. 5
figure 5

Identification of DEGs between DLAT-high and DLAT-low groups and functional enrichment analysis. A Volcano plot of DEGs. B Heatmap plot of the top 200 DEGs. C-F GO and KEGG pathway analysis of DEGs. G Differences in pathway activities scored per patient by GSVA between DLAT-high and DLAT-low groups

Immunological analysis of DLAT in PAAD

We first investigated the impact of DLAT on immune cell infiltration in the TIMER database, and observed that DLAT CNV has a significant association with the infiltration levels of CD4 + T cell and neutrophil (Fig. 6A). Meanwhile, DLAT expression was positively correlated with the infiltration levels of B cell (r = 0.352, p < 0.001), CD8 + T cell (r = 0.589, p < 0.001), macrophage (r = 0.482, p < 0.001), neutrophil (r = 0.409, p < 0.001), and dendritic cell (r = 0.548, p < 0.001) (Fig. 6B). In addition, DLAT expression exhibited positive correlations with most immune cells analyzed by ssGSEA algorithm (Fig. 6C). We then calculated the TME scores based on "ESTIMATE" package, and revealed close associations between DLAT and stromal score (r = 0.26, p < 0.001), ESTIMATE score (r = 0.18, p = 0.02), as well as tumor purity (r = -0.17, p = 0.03) (Fig. 6D). Moreover, DLAT was positively related to steps 1–6 of cancer-immunity cycle and the enrichment scores of immunotherapy-predicted pathways gene signatures, such as cell cycle, Fanconi anemia pathway, homologous recombination, and others (Fig. 6E).

Fig. 6
figure 6

Correlation analysis between DLAT expression and immune cell infiltration. A Immune cell infiltration under various copy numbers of DLAT. B Correlations between DLAT and immune cell infiltration based on TIMER database. C Correlations between DLAT and immune cell infiltration calculated by ssGSEA algorithm. D Correlations between DLAT and stromal, immune, ESTIMATE scores, as well as tumor purity. E Correlations between DLAT and immunotherapy-predicted pathways (left) and cancer-immunity cycle (right) (***p < 0.001)

According to previous results, DLAT was likely to produce a favorable effect on anticancer immune response in PAAD. Immune cell infiltration in DLAT-high and DLAT-low groups was shown in Fig. 7A. We integrated 5 main stream immunoinformatic methods, including TIMER, CIBERSORT, MCPCOUNTER, QUANTISEQ, and EPIC, to elucidate tumor-infiltrating immune cell types in different groups. The infiltration levels of immune cells were different between DLAT-high and DLAT-low groups (Fig. 7B).

Fig. 7
figure 7

Immunological characteristics in DLAT-high and DLAT-low groups. A Proportions of 22 immune cell types in DLAT-high and DLAT-low groups. B Heatmap of immune cell infiltration in DLAT-high and DLAT-low groups. C Correlations between DLAT and inhibitory immune checkpoints. D Correlations between DLAT and PD-L1 expression. E Comparison of PD-L1 expression in DLAT-high and DLAT-low groups. F Submap algorithm evaluated the expression similarity between TCGA-PAAD patients and melanoma patients receiving anti-PD-1 and anti-CTLA-4 therapies (R represents responsive, while noR represents no responsive; ***p < 0.001)

Immune checkpoints played an immunosuppressive role in immune response by inhibiting immune cell activation. In this research, we found that DLAT expression was positively related to common inhibitory immune checkpoints, such as PD-L1, PD-L2, CTLA-4, TIM-3, and TIGIT (Fig. 7C). PD-L1 expression was an important biomarker to predict immunotherapy response. We, therefore, examined the relationship between DLAT and PD-L1, and demonstrated a positive correlation between them (r = 0.57, p < 0.001) (Fig. 7D). Compared to DLAT-low group, PD-L1 was statistically upregulated in DLAT-high group (Fig. 7E). Subsequently, we compared the expression profile of TCGA-PAAD patients with a melanoma dataset containing 47 patients who received anti-PD-1 and anti-CTLA-4 treatments [35,36,37]. It was observed that DLAT-high patients were more responsive to anti-CTLA-4 therapy than DLAT-low patients (Nominal p = 0.008) (Fig. 7F). In addition, patients in DLAT-high group were more promising to respond to anti-PD-L1 therapy (Supplementary Fig. 1A). However, when comparing the expression profile of DLAT-high group with PD-1-response group, the results showed insignificant differences, indicating that DLAT exhibited little impact on the anti-PD-1 therapeutic response (Supplementary Fig. 1B-D). The diverse results may be attributed to the stronger correlation between DLAT and PD-L1/CTLA-4 than it between DLAT and PD-1 (Supplementary Fig. 1E).

Development and validation of DLAT-based risk score model

To create the DLAT-based model, we defined 427 prognostic DEGs by univariate Cox regression analysis (Supplementary Table 4) and selected 13 optimal genes through LASSO regression analysis (Fig. 8A-B). According to the 13-gene model, patients were scored using the following formula and then split into high-risk and low-risk groups based on the median risk score. Risk score = TSPOAP1 * (-0.02367876) + ANLN * 0.011937625 + TLE2 * (-0.08933642) + GNB3 * (-0.036416417) + GLTPD2 * (-0.018670796) + BPIFB4 * (-0.003192766) + FAM111B * 0.102912035 + INSYN2B * 0.070370212 + RNU6-892P * (-0.017035372) + LINC01940 * 0.044119945 + ANKRD18B * 0.034310234 + CCDC188 * (-0.012306307) + HNRNPD-DT * (-0.037777293).

Fig. 8
figure 8

Construction and validation of risk score model. A, B Thirteen modeling genes determined by LASSO Cox regression. C Distribution of risk score and survival status in TCGA training cohort. D Kaplan–Meier curves for high-risk and low-risk patients in TCGA training cohort. E Time-dependent ROC curves for predicting the 1-, 3-, and 5-year survival status in TCGA training cohort. F Distribution of risk score, G Kaplan–Meier curves, and H Time-dependent ROC curves in TCGA validation cohort. I Sankey diagram of different CRG clusters, DLAT subgroups, risk scores, and survival outcomes. J Comparison of risk scores in cluster A, B, and C. K Comparison of risk scores between DLAT-high and DLAT-low groups (ns represents no significance, ****p < 0.0001)

The risk score increased along with the death rates in the training cohort (the whole dataset, n = 171) (Fig. 8C). Moreover, the OS was significantly prolonged in low-risk group (log-rank test, p < 0.0001) (Fig. 8D). The ROC analysis demonstrated that this model had a high accuracy in predicting 1-, 3-, 5- year survival status, with the AUC values of 0.746, 0.801, and 0.878, respectively (Fig. 8E). Additionally, this risk model performed well in the validation cohort (randomly selected, n = 85) (Fig. 8F-H).

The Sankey diagram displayed the proportions of patients with different CRG clusters, DLAT expression, risk scores, and survival outcomes (Fig. 8I). Patients in cluster C showed higher risk scores than those in other clusters, however, the difference was not significant (Fig. 8J). Furthermore, DLAT-high patients exhibited statistically higher risk scores than DLAT-low patients (Fig. 8K).

Validation of the expression of DLAT in PAAD

We performed RT-qPCR to further confirm the differential expression of DLAT. The results exhibited that DLAT mRNA was significantly increased in tumor cell lines compared to normal cell line (Supplementary Fig. 2). Moreover, we collected pancreatic tumor tissues and adjacent normal tissues from 12 patients. Their clinical characteristics were summarized in Supplementary Table 5. The IHC assays demonstrated that DLAT was highly expressed in pancreatic cancer tissues (Fig. 9).

Fig. 9
figure 9

Representative IHC images of DLAT staining in pancreatic tumor and adjacent normal tissues

Discussion

PAAD is a common malignant tumor worldwide with an augmented incidence over the past decades. Despite advances in the treatment of PAAD, there was very little improvement in OS, with a 5-year survival rate of 9% [41]. Thus, it is ultimate to discover novel therapeutic targets to improve the clinical outcomes of PAAD patients.

Tsvetkov et al. revealed a novel cell death pathway termed "cuproptosis" induced by intracellular copper. Unlike other forms, cuproptosis occurs when excessive copper ions bind to mitochondrial lipoylated components, causing cell death by proteotoxic stress. However, the role that cuproptosis plays in PAAD remains to be illuminated.

In this work, we carried out the bioinformatic analysis of CRGs using public data acquired from TCGA, GTEx, and GEO databases. We examined the expression of CRGs and found that most genes exhibited increased expression levels in PAAD tissues. Subsequently, patients were classified into 3 subgroups using consensus clustering analysis. Patients in cluster A showed the most favorable behaviors with the longest OS time. The enrichment analysis indicated that upregulated DEGs in cluster A were mainly enriched in cellular metabolism pathways (such as retinol metabolism, drug metabolism by cytochrome P450, and xenobiotics metabolism by cytochrome P450), and tumor-related signaling pathways (such as Wnt/β catenin pathways, DNA repair, p53 pathways, and unfolded protein response). Previous researches have fully uncovered that these pathways played a significant role in cancer development by regulating cellular metabolism [42], cell growth [43], apoptosis [44], and tumor immunity [45]. Accordingly, CRGs may participate in the development of PAAD via these classical pathways.

The "Warburg effect" is a unique metabolic feature in tumors, that is, the energy required for cancer cellular processes is primarily generated by aerobic glycolysis rather than mitochondrial oxidative phosphorylation [46]. A newly developed viewpoint showed that the changes in tumor metabolism may be partly connected to impaired mitochondrial function caused by the inhibition of pyruvate dehydrogenase complex (PDC) [47]. PDC is an enzyme complex responsible to regulate glucose oxidative metabolism. It is composed of pyruvate dehydrogenase (E1), dihydrolipoamide acetyltransferase (E2), and dihydrolipoamide dehydrogenase (E3) [48]. These components work orderly in the inner mitochondrial membrane to catalyze the conversion of oxidative decarboxylation of pyruvate to acetyl coenzyme A, thereby interconnecting glycolysis and TCA cycle. Nevertheless, the inactivation of PDC controlled by pyruvate dehydrogenase kinase limits the synthesis of acetyl coenzyme A, so that glucose oxidation is blocked and metabolic disorders are developed [49].

The dihydrolipoamide acetyltransferase (PDC-E2) is encoded by DLAT gene. Dihydrolipoamide acetyltransferase has been widely described in an autoimmune liver disease, primary biliary cholangitis, in which patients lost their immune tolerance to PDC-E2 [50, 51]. Autoreactive T lymphocytes were increased in primary biliary cholangitis and destroyed bile duct epithelial cells, leading to destructive lymphocytic cholangitis, ultimately resulting in cirrhosis and even liver failure. Moreover, pyruvate dehydrogenase deficiency, a disease characterized by primary lactic acidosis and neurological dysfunction, was related to DLAT gene mutation [52]. As for its role in tumor development, Shan et al. reported that DLAT emerged as an upstream acetyltransferase for K76 and activated oxidative pentose phosphate pathway, finally promoting cell growth in cancer cells [53]. In addition, DLAT was increased in gastric cancer, and siRNA-mediated knockdown of DLAT could increase the pyruvate levels [54]. DLAT-mediated glycolysis reprograming also contributed to the carcinogenesis of non-small cell lung cancer [55]. Furthermore, DLAT exhibited prognostic values in hepatocellular carcinoma patients based on public databases and clinical samples [56].

Although CRGs have been studied in PAAD, DLAT has hardly been considered as a special gene for systemic analysis. Researchers generally analyzed the CRGs as a whole rather than focusing on one gene [57,58,59,60]. Compared to other genes, DLAT possessed the lowest p value in prognostic analysis, indicating that it may play a more important role in the prognosis of PAAD. Therefore, a comprehensive analysis of DLAT was needed. In this work, we found that DLAT served as an independent risk factor for OS in PAAD. KEGG pathway analysis revealed that DLAT was primarily correlated to TCA cycle and pancreatic cancer pathways, which was in accordance with the fact that cuproptosis was related to lipoylated TCA cycle proteins, further demonstrating the rationality and feasibility of our work. Functional enrichment analysis also revealed that DLAT was engaged in cell division, indicating that DLAT might promote tumor progression by activating proliferative processes, such as mismatch repair, DNA replication, nucleotide excision repair, and cell cycle.

As an essential trace element, copper is required for maintaining enzyme functions and body homeostasis [13]. The levels of copper were usually increased in malignant tumors [61, 62]. Copper was involved in the regulation of cellular and humoral immunity, and both excessive or insufficient copper ions will impair cellular functions [63, 64]. We evaluated the relationship between DLAT and multiple immunological characteristics in PAAD. As a result, DLAT showed significantly positive correlations with immune cell infiltration, most cancer-immunity cycle steps, multiple immunotherapy-predicted pathways, and common inhibitory immune checkpoints. Voli et al. reported that copper could enhance PD-L1 expression and promote cancer immune evasion in vitro, while copper chelators could increase CD8 + T cell infiltration and prolong survival time in vivo [65]. Similarly, DLAT-high patients with increased PD-L1 expression had worse survival status than DLAT-low patients. Furthermore, DLAT-high patients were more sensitive to respond to anti-CTLA-4/PD-L1 treatment than DLAT-low patients, implying that DLAT may function as a biomarker to predict immunotherapy response.

In actual clinical practice, with the help of risk score model, doctors can better comprehend the status of patients to make precise treatment decisions. Since the introduction of cuproptosis, many researchers have developed the risk score models to predict the clinical outcomes of patients with PAAD. However, the reported models exhibited relatively lower AUC values (range from 0.68–0.84) [57,58,59]. Therefore, a more effective model was needed. Considering the potential effect of DLAT on clinical outcomes, we constructed a risk model based on DLAT to assess the individual survival risks. It was found that the model we built showed higher AUC values than others (Supplementary Table 6). Thus, the DLAT-based model possessed better effectiveness in predicting prognosis and could be more useful in clinical practice.

Overall, we comprehensively evaluated the prognostic values and biological function of DLAT in PAAD. In addition, we revealed the relationship between DLAT and multiple immunological characteristics, which could provide new ideas for tumor immunotherapy. We also constructed a DLAT-based model to predict patients’ prognosis. However, our study still has certain limitations. Firstly, the prognostic and immunotherapeutic implications of DLAT were analyzed using public databases, which may cause a certain bias. Secondly, the impact of DLAT on tumor progression needs to be further validated by cell function experiments. In conclusion, more experimental and clinical data are required in the future to verify the above conclusions.

Conclusion

We conducted a systematic analysis of DLAT in PAAD by integrated bioinformatic methods and clinical validation. It was identified that DLAT could emerge as an independent risk factor for survival, with the potential to guide personalized immunotherapy in PAAD. Moreover, the DLAT-based model exhibited high accuracy in prognosis prediction. Our findings could help illuminate the function of DLAT and contribute to accurate treatment of PAAD.

Availability of data and materials

RNA-sequencing data can be obtained from the UCSC Xena (https://xenabrowser.net/datapages/) (GDC TCGA PAAD and GTEx). The study’s GEO database (GSE62452, GSE71729, GSE15471, GSE16515, GSE78220, and GSE91061) can be downloaded from https://www.ncbi.nlm.nih.gov/gds/. The original data presented in the study were included in Supplementary Information. Further inquiries can be directed to the corresponding author.

Abbreviations

PAAD:

Pancreatic adenocarcinoma

CRGs:

Cuproptosis-related genes

DLAT:

Dihydrolipoamide acetyltransferase

RT-qPCR:

Quantitative reverse transcriptase polymerase chain reaction

IHC:

Immunohistochemistry

PD-1:

Programmed cell death-1

PD-L1:

Programmed cell death-ligand 1

CTLA-4:

Cytotoxic T-lymphocyte antigen-4

TME:

Tumor microenvironment

TCA:

Tricarboxylic acid

PPI:

Protein-protein interaction

DEGs:

Differentially expressed genes

TCGA:

The Cancer Genome Atlas

GTEx:

Genotype-tissue expression

GEO:

Gene Expression Omnibus

FPKM:

Fragments per kilobase million

TPM:

Transcripts per million

CNV:

Copy number variation

GO:

Gene ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

GSEA:

Gene set enrichment analysis

GSVA:

Gene set variation analysis

ssGSEA:

Single sample gene set enrichment analysis

OS:

Overall survival

HR:

Hazard ratio

CI:

Confidence interval

ICB:

Immune checkpoint blockade

ROC:

Receiver operating characteristic

AUC:

Areas under the curve

PDC:

Pyruvate dehydrogenase complex

References

  1. Siegel RL, Miller KD, Wagle NS, Jemal A. Cancer statistics, 2023. CA Cancer J Clin. 2023;73(1):17–48.

    Article  PubMed  Google Scholar 

  2. Li Y, Liu Z, Zhang Y. Expression and prognostic impact of FZDs in pancreatic adenocarcinoma. BMC Gastroenterol. 2021;21(1):79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Klein AP. Pancreatic cancer epidemiology: understanding the role of lifestyle and inherited risk factors. Nat Rev Gastroenterol Hepatol. 2021;18(7):493–502.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Park W, Chawla A, O’Reilly EM. Pancreatic Cancer: A Review. JAMA. 2021;326(9):851–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Rawla P, Sunkara T, Gaduputi V. Epidemiology of Pancreatic Cancer: Global Trends, Etiology and Risk Factors. World J Oncol. 2019;10(1):10–27.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Mizrahi JD, Surana R, Valle JW, Shroff RT. Pancreatic cancer. Lancet. 2020;395(10242):2008–20.

    Article  CAS  PubMed  Google Scholar 

  7. Neoptolemos JP, Kleeff J, Michl P, Costello E, Greenhalf W, Palmer DH. Therapeutic developments in pancreatic cancer: current and future perspectives. Nat Rev Gastroenterol Hepatol. 2018;15(6):333–48.

    Article  PubMed  Google Scholar 

  8. Galluzzi L, Chan TA, Kroemer G, Wolchok JD, López-Soto A. The hallmarks of successful anticancer immunotherapy. Sci Transl Med. 2018;10(459):eaat 7807.

    Article  Google Scholar 

  9. Schizas D, Charalampakis N, Kole C, Economopoulou P, Koustas E, Gkotsis E, et al. Immunotherapy for pancreatic cancer: A 2020 update. Cancer Treat Rev. 2020;86:102016.

    Article  CAS  PubMed  Google Scholar 

  10. Fan J, Wang M, Chen H, Shang D, Das JK, Song J. Current advances and outlooks in immunotherapy for pancreatic ductal adenocarcinoma. Mol Cancer. 2020;19(1):32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Tsvetkov P, Coy S, Petrova B, Dreishpoon M, Verma A, Abdusamad M, et al. Copper induces cell death by targeting lipoylated TCA cycle proteins. Science. 2022;375(6586):1254–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Tang D, Kang R, Berghe TV, Vandenabeele P, Kroemer G. The molecular machinery of regulated cell death. Cell Res. 2019;29(5):347–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Kim B, Nevitt T, Thiele DJ. Mechanisms for copper acquisition, distribution and regulation. Nat Chem Biol. 2008;4(3):176–85.

    Article  CAS  PubMed  Google Scholar 

  14. Wang Y, Zhang Y, Wang L, Zhang N, Xu W, Zhou J, et al. Development and experimental verification of a prognosis model for cuproptosis-related subtypes in HCC. Hepatol Int. 2022;16(6):1435–47.

    Article  PubMed  Google Scholar 

  15. Zhang G, Chen X, Fang J, Tai P, Chen A, Cao K. Cuproptosis status affects treatment options about immunotherapy and targeted therapy for patients with kidney renal clear cell carcinoma. Front Immunol. 2022;13:954440.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Sha S, Si L, Wu X, Chen Y, Xiong H, Xu Y, et al. Prognostic analysis of cuproptosis-related gene in triple-negative breast cancer. Front Immunol. 2022;13:922780.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Yang M, Zheng H, Xu K, Yuan Q, Aihaiti Y, Cai Y, et al. A novel signature to guide osteosarcoma prognosis and immune microenvironment: Cuproptosis-related lncRNA. Front Immunol. 2022;13:919231.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Jiang R, Huan Y, Li Y, Gao X, Sun Q, Zhang F, et al. Transcriptional and genetic alterations of cuproptosis-related genes correlated to malignancy and immune-infiltrate of esophageal carcinoma. Cell Death Discov. 2022;8(1):370.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Li G, Luo Q, Wang X, Zeng F, Feng G, Che G. Deep learning reveals cuproptosis features assist in predict prognosis and guide immunotherapy in lung adenocarcinoma. Front Endocrinol (Lausanne). 2022;13:970269.

    Article  PubMed  Google Scholar 

  20. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi BVSK, et al. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia. 2017;19(8):649–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43(Database issue):D447–52.

    Article  CAS  PubMed  Google Scholar 

  23. Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(Web Server issue):W214–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 2017;77(21):e108–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, TorresGarcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.

    Article  PubMed  Google Scholar 

  26. Hu J, Yu A, Othmane B, Qiu D, Li H, Li C, et al. Siglec15 shapes a non-inflamed tumor microenvironment and predicts the molecular subtype in bladder cancer. Theranostics. 2021;11(7):3089–108.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Mariathasan S, Turley SJ, Nickles D, Castiglioni A, Yuen K, Wang Y, et al. TGFβ attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature. 2018;554(7693):544–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Chen DS, Mellman I. Oncology meets immunology: the cancer-immunity cycle. Immunity. 2013;39(1):1–10.

    Article  PubMed  Google Scholar 

  29. Xu L, Deng C, Pang B, Zhang X, Liu W, Liao G, et al. TIP: A Web Server for Resolving Tumor Immunophenotype Profiling. Cancer Res. 2018;78(23):6575–80.

    Article  CAS  PubMed  Google Scholar 

  30. Li B, Severson E, Pignon J, Zhao H, Li T, Novak J, et al. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol. 2016;17(1):174.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016;17(1):218.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Finotello F, Mayer C, Plattner C, Laschober G, Rieder D, Hackl H, et al. Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data. Genome Med. 2019;11(1):34.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Racle J, de Jonge K, Baumgaertner P, Speiser DE, Gfeller D. Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression data. Elife. 2017;6:e26476.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Hoshida Y, Brunet J, Tamayo P, Golub TR, Mesirov JP. Subclass mapping: identifying common subtypes in independent disease data sets. PLoS One. 2007;2(11):e1195.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Lu X, Jiang L, Zhang L, Zhu Y, Hu W, Wang J, et al. Immune Signature-Based Subtypes of Cervical Squamous Cell Carcinoma Tightly Associated with Human Papillomavirus Type 16 Expression, Molecular Features, and Clinical Outcome. Neoplasia. 2019;21(6):591–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Roh W, Chen P, Reuben A, Spencer CN, Prieto PA, Miller JP, et al. Integrated molecular analysis of tumor biopsies on sequential CTLA-4 and PD-1 blockade reveals markers of response and resistance. Sci Transl Med. 2017;9(379):eaah3560.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Braun DA, Hou Y, Bakouny Z, Ficial M, Sant’ Angelo M, Forman J, et al. Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell renal cell carcinoma. Nat Med. 2020;26(6):909–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Hugo W, Zaretsky JM, Sun L, Song C, Moreno BH, Hu-Lieskovan S, et al. Genomic and Transcriptomic Features of Response to Anti-PD-1 Therapy in Metastatic Melanoma. Cell. 2016;165(1):35–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Riaz N, Havel JJ, Makarov V, Desrichard A, Urba WJ, Sims JS, et al. Tumor and Microenvironment Evolution during Immunotherapy with Nivolumab. Cell. 2017;171(4):934-949.e916.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Khalaf N, ElSerag HB, Abrams HR, Thrift AP. Burden of Pancreatic Cancer: From Epidemiology to Practice. Clin Gastroenterol Hepatol. 2021;19(5):876–84.

    Article  PubMed  Google Scholar 

  42. Chen X, CubillosRuiz JR. Endoplasmic reticulum stress signals in the tumour and its microenvironment. Nat Rev Cancer. 2021;21(2):71–88.

    Article  CAS  PubMed  Google Scholar 

  43. Jeggo PA, Pearl LH, Carr AM. DNA repair, genome stability and cancer: a historical perspective. Nat Rev Cancer. 2016;16(1):35–42.

    Article  CAS  PubMed  Google Scholar 

  44. Vazquez A, Bond EE, Levine AJ, Bond GL. The genetics of the p53 pathway, apoptosis and cancer therapy. Nat Rev Drug Discov. 2008;7(12):979–87.

    Article  CAS  PubMed  Google Scholar 

  45. Gattinoni L, Ji Y, Restifo NP. Wnt/beta-catenin signaling in T-cell immunity and cancer immunotherapy. Clin Cancer Res. 2010;16(19):4695–701.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009;324(5930):1029–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Saunier E, Benelli C, Bortoli S. The pyruvate dehydrogenase complex in cancer: An old metabolic gatekeeper regulated by new pathways and pharmacological agents. Int J Cancer. 2016;138(4):809–17.

    Article  CAS  PubMed  Google Scholar 

  48. Patel MS, Nemeria NS, Furey W, Jordan F. The pyruvate dehydrogenase complexes: structure-based function and regulation. J Biol Chem. 2014;289(24):16615–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Harris RA, BowkerKinley MM, Huang B, Wu P. Regulation of the activity of the pyruvate dehydrogenase complex. Adv Enzyme Regul. 2002;42:249–59.

    Article  CAS  PubMed  Google Scholar 

  50. Joplin R, Lindsay JG, Johnson GD, Strain A, Neuberger J. Membrane dihydrolipoamide acetyltransferase (E2) on human biliary epithelial cells in primary biliary cirrhosis. Lancet. 1992;339(8785):93–4.

    Article  CAS  PubMed  Google Scholar 

  51. Hirschfield GM, Gershwin ME. The immunobiology and pathophysiology of primary biliary cirrhosis. Annu Rev Pathol. 2013;8:303–30.

    Article  CAS  PubMed  Google Scholar 

  52. Head RA, Brown RM, Zolkipli Z, Shahdadpuri R, King MD, Clayton PT, et al. Clinical and genetic spectrum of pyruvate dehydrogenase deficiency: dihydrolipoamide acetyltransferase (E2) deficiency. Ann Neurol. 2005;58(2):234–41.

    Article  CAS  PubMed  Google Scholar 

  53. Shan C, Elf S, Ji Q, Kang H, Zhou L, Hitosugi T, et al. Lysine acetylation activates 6-phosphogluconate dehydrogenase to promote tumor growth. Mol Cell. 2014;55(4):552–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Goh WQJ, Ow GS, Kuznetsov VA, Chong S, Lim YP. DLAT subunit of the pyruvate dehydrogenase complex is upregulated in gastric cancer-implications in cancer therapy. Am J Transl Res. 2015;7(6):1140–51.

    PubMed  PubMed Central  Google Scholar 

  55. Chen Q, Wang Y, Yang L, Sun L, Wen Y, Huang Y, et al. PM2.5 promotes NSCLC carcinogenesis through translationally and transcriptionally activating DLAT-mediated glycolysis reprograming. J Exp Clin Cancer Res. 2022;41(1):229.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Bai W, Liu J, Li M, Yang X, Wang Y, Wang G, et al. A Novel Cuproptosis-Related Signature Identified DLAT as a Prognostic Biomarker for Hepatocellular Carcinoma Patients. World J Oncol. 2022;13(5):299–310.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Huang X, Zhou S, Tóth J, Hajdu A. Cuproptosis-related gene index: A predictor for pancreatic cancer prognosis, immunotherapy efficacy, and chemosensitivity. Front Immunol. 2022;13:978865.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Xu Y, Li H, Lan A, Wu Q, Tang Z, Shu D, et al. Cuprotosis-Related Genes: Predicting Prognosis and Immunotherapy Sensitivity in Pancreatic Cancer Patients. J Oncol. 2022;2022:2363043.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Li J, Yin J, Li W, Wang H, Ni B. Molecular subtypes based on cuproptosis-related genes and tumor microenvironment infiltration characteristics in pancreatic adenocarcinoma. Cancer Cell Int. 2023;23(1):7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Chen Y, Zou X, Ma M, Liu Y, Wang R, Dai Z, et al. Expression Profiles of Cuproptosis-Related Genes Determine Distinct Subtypes of Pancreatic Ductal Adenocarcinoma. Curr Oncol. 2023;30(2):1648–62.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Aubert L, Nandagopal N, Steinhart Z, Lavoie G, Nourreddine S, Berman J, et al. Copper bioavailability is a KRAS-specific vulnerability in colorectal cancer. Nat Commun. 2020;11(1):3701.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Kaiafa GD, Saouli Z, Diamantidis MD, Kontoninas Z, Voulgaridou V, Raptaki M, et al. Copper levels in patients with hematological malignancies. Eur J Intern Med. 2012;23(8):738–41.

    Article  CAS  PubMed  Google Scholar 

  63. Percival SS. Copper and immunity. Am J Clin Nutr. 1998;67(5 Suppl):1064S-1068S.

    Article  CAS  PubMed  Google Scholar 

  64. Arredondo M, Núñez MT. Iron and copper metabolism. Mol Aspects Med. 2005;26(4–5):313–27.

    Article  CAS  PubMed  Google Scholar 

  65. Voli F, Valli E, Lerra L, Kimpton K, Saletta F, Giorgi FM, et al. Intratumoral Copper Modulates PD-L1 Expression and Influences Tumor Immune Evasion. Cancer Res. 2020;80(19):4129–44.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (No. 82073301), Haiyan Foundation of Harbin Medical University Cancer Hospital (No. JJMS2021-03, JJMS2021-22, and JJMS2021-23), and Beijing Medical Award Foundation (No. YXJL-2021–0268-0332).

Author information

Authors and Affiliations

Authors

Contributions

RX and XZ designed the study. XZ, YZ, and JH performed the data analysis, and XZ wrote the manuscript. XY, HX, and ZB provided software support. HZ, YS, and RW downloaded the data. XD, RM, and XL contributed to the experiments of this study. RX and JZ revised the manuscript. RX provided the foundation and support. All authors contributed to the article and approved the final manuscript.

Corresponding authors

Correspondence to Jiuxin Zhu or Rui Xie.

Ethics declarations

Ethics approval and consent to participate

The study protocol was met with the declaration of Helsinki and was approved by the Ethics Committee of Harbin Medical University Cancer Hospital (No. KY2023-32). All methods were performed in accordance with relevant guidelines and regulations, and written informed consent was obtained from all the study participants.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1:

Supplementary Table 1. Correlation between DLAT and other genes in pancreatic cancer. Supplementary Table 2. GSEA of DLAT co-expression network based on "HALLMARK" gene sets. Supplementary Table 3. DEGs between DLAT-high and DLAT-low groups. Supplementary Table 4. Univariate Cox regression analysis identified 427 prognostic DEGs between DLAT-high and DLAT-low groups. Supplementary Table 5. Clinical characteristics of 12 pancreatic cancer patients. Supplementary Table 6. Comparison of AUC values between DLAT-based model and other models.

Additinal file 2:

Supplementary Figure 1. Prediction of immunotherapy response differences between DLAT-high and DLAT-low groups in (A) IMvigor210, (B) CheckMate, (C) GSE78220, and (D) GSE91061 cohorts (R represents responsive, while noR represents no responsive). (E) Correlation analysis between DLAT and PD-L1, CTLA-4, as well as PD-1. Supplementary Figure 2. The mRNA expression of DLAT in pancreatic cancer cell lines (BxPC-3 and PANC-1) and normal cell line (HPDE6-C7) (***p < 0.001).

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, X., Zhou, Y., Hu, J. et al. Comprehensive analysis identifies cuproptosis-related gene DLAT as a potential prognostic and immunological biomarker in pancreatic adenocarcinoma. BMC Cancer 23, 560 (2023). https://doi.org/10.1186/s12885-023-11042-7

Download citation

  • Received:

  • Accepted:

  • Published:

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

Keywords