Xu H, Du QC, Wang XY, Zhou L, Wang J, Ma YY, Liu MY, Yu H. Comprehensive analysis of the relationship between cuproptosis-related genes and esophageal cancer prognosis. World J Clin Cases 2022; 10(33): 12089-12103 [DOI: 10.12998/wjcc.v10.i33.12089]
Corresponding Author of This Article
Hua Yu, MM, Associate Chief Physician, General Surgery, Shanghai Fourth People’s Hospital, School of Medicine, Tongji University, No. 1279 Sanmen Road, Shanghai 200434, China. firstname.lastname@example.org
Checklist of Responsibilities for the Scientific Editor of This Article
This article is an open-access article which was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Author contributions: Xu H and Du QC analyzed the data and wrote the manuscript; Yu H designed the study; Zhou L, Wang J, Ma YY, and Liu MY collected the data and corrected the paper; all authors have read and approved the final manuscript.
Supported bythe Shanghai Fourth People’s Hospital Affiliated to Tongji University School of Medicine Talent Boosting Plan, No. SY-XKZT-2020-3007.
Institutional review board statement: The data for the study came from public databases and did not involve blood or tissue samples from humans or animals. Therefore, there were no ethical issues involved in this study.
Conflict-of-interest statement: The authors declare no conflicts of interest.
Data sharing statement: No additional data are available.
STROBE statement: The authors have read the STROBE Statement—checklist of items, and the manuscript was prepared and revised according to the STROBE Statement—checklist of items.
Open-Access: This article is an open-access article that was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution NonCommercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: https://creativecommons.org/Licenses/by-nc/4.0/
Corresponding author: Hua Yu, MM, Associate Chief Physician, General Surgery, Shanghai Fourth People’s Hospital, School of Medicine, Tongji University, No. 1279 Sanmen Road, Shanghai 200434, China. email@example.com
Received: June 18, 2022 Peer-review started: June 18, 2022 First decision: August 6, 2022 Revised: August 18, 2022 Accepted: October 9, 2022 Article in press: October 9, 2022 Published online: November 26, 2022
Esophageal cancer is one of the most common malignant tumors of the digestive system, with a 5-year survival rate of 15% to 50%. Cuproptosis, a unique kind of cell death driven by protein lipoylation, is strongly connected to mitochondrial metabolism. The clinical implications of cuproptosis-related genes in esophageal cancer, however, are mainly unknown.
To identify cuprotosis-related genes that are differentially expressed in esophageal cancer and investigate their prognostic significance.
With |log fold change| > 1 and false discovery rate < 0.05 as criteria, the Wilcox test was used to evaluate the differentially expressed genes between 151 tumor tissues and 151 normal esophageal tissues. Cuproptosis-related genes were selected to be linked with prognosis using univariate Cox regression analysis. Genes were separated into high- and low- expression groups based on their cutoff value of gene expression, and the correlation between the two groups and overall survival or progression-free survival was investigated using the log-rank test. The C-index, calibration curve, and receiver operator characteristic (ROC) curve were used to assess a nomogram containing clinicopathological characteristics and cuproptosis-related genes.
Pyruvate dehydrogenase A1 (PDHA1) was found to be highly correlated with prognosis in univariate Cox regression analysis (hazard ratio = 22.96, 95% confidence interval = 3.09-170.73; P = 0.002). According to Kaplan-Meier survival curves, low expression of PDHA1 was associated with a better prognosis (log-rank P = 0.0007). There was no significant correlation between PDHA1 expression and 22 different types of immune cells. Tumor necrosis factor superfamily member 15 (TNFSF15) (P = 3.2 × 10-6; r = 0.37), TNFRSF14 (P = 8.1 × 10-8; r = 0.42), H long terminal repeat-associating 2 (P = 6.0 × 10-8; r = 0.42) and galectin 9 (P = 3.1 × 10-6; r = 0.37) were all found to be considerably greater in the high PDHA1 expression group, according to an analysis of genes related to 47 immunological checkpoints. The low PDHA1 expression group had significantly lower levels of cluster of differentiation 44 (CD44) (P = 0.00028; R = -0.29), TNFRSF18 (P = 1.2 × 10-5; R = -0.35), programmed cell death 1 ligand 2 (P = 0.0032; R = -0.24), CD86 (P = 0.018; R = -0.19), and CD40 (P = 0.0047; R = -0.23), and the differences were statistically significant. We constructed a prognostic nomogram incorporating pathological type, tumor-node-metastasis stage, and PDHA1 expression, and the C-index, calibration curve, and ROC curve revealed that the nomogram’s predictive performance was good.
Cuproptosis-related genes can be used as a prognostic predictor for esophageal cancer patients, providing novel insights into cancer treatment.
Core Tip: Esophageal carcinoma has a poor prognosis and is one of the major causes of cancer-related deaths worldwide. Despite recent advancements in the surgical and pharmacological treatment of esophageal cancer, the prognosis remains poor. Copper toxicity has been linked to the incidence and progression of esophageal cancer in numerous studies. At the gene level, however, the probable biochemical mechanism is unknown. We included 19 cuproptosis-related genes and screened a gene that could successfully predict the prognosis of esophageal cancer by statistical analysis to further elucidate the role of cuproptosis-related genes in impacting the prognosis of esophageal cancer.
Citation: Xu H, Du QC, Wang XY, Zhou L, Wang J, Ma YY, Liu MY, Yu H. Comprehensive analysis of the relationship between cuproptosis-related genes and esophageal cancer prognosis. World J Clin Cases 2022; 10(33): 12089-12103
Esophageal cancer is a major global health issue, and its incidence is rapidly increasing. Esophageal cancer is classified into two types: esophageal squamous cell carcinoma (ESCC), which accounts for 90% of all occurrences, and esophageal adenocarcinoma (EAC). In recent years, epidemiological research has revealed that the incidence of EAC has grown 3-4 times, with the proportion increasing[2,3]. Despite significant advances in the diagnosis and treatment of esophageal cancer, the mortality rate varies from 15% to 20%, placing it fourth among all cancer-related deaths.
Esophageal cancer is the result of a complex process involving various causes and polygene alterations. Using high-throughput sequencing technologies, a comprehensive mutation catalog was evaluated, and substantial genetic alterations have been discovered in the malignancies. Gene alterations are often linked to aberrant expression and are becoming more essential for the early diagnosis and prognosis of esophageal cancer. Currently, several gene expression products are employed as indicators for esophageal cancer diagnosis and prognosis[5,6]. Somatic mutations in tumor protein p53 (TP53) have been found in more than 83% of ESCCs. Adenocarcinoma and squamous cell carcinoma both have TP53 point mutations[7,8]. In addition, numerous cell cycle-controlling genes are overexpressed in ESCC. For example, cyclin-dependent kinase 4/cyclin-dependent kinase 6 accounted for 23.6%, murine double minute 25.7%, and cyclin D1 46.4%, showing that the above components are implicated in the incidence and development of ESCC. As a result, there is a pressing need to uncover genetic anomalies in esophageal cancer and understand their molecular basis to enhance early diagnosis and minimize esophageal cancer mortality.
Previous studies have reported several types of precisely programmed cell death, including apoptosis, pyroptosis, necrosis, and ferroptosis. Similar to iron, copper is a trace metal in cells that plays an integral role in maintaining protein functions. Excessive copper can cause cytotoxicity, but the exact mechanism is unclear. Tsvetkov et al discovered that the copper carrier elesclomol, which was originally used to treat cancer, kills cells in excess of copper. Elesclomol does not trigger cell death on its own, suggesting that copper toxicity is to blame. Dar et al found that patients with esophageal cancer had considerably higher mean blood copper levels than controls, with a mean copper concentration of 169 g/dL in the cancer group and 149 g/dL in the control group. Therefore, we hypothesized that copper shortage or excess is linked to the occurrence and progression of esophageal cancer. Copper, however, has been linked to the development and progression of esophageal cancer in few studies.
There has not been any research on the role of cuproptosis-related genes in the genesis, progression, and metastasis of esophageal cancer to date. Cuproptosis-related genes and their processes need to be better understood in order to improve the prognosis of malignant tumors and uncover novel treatment targets. Bioinformatics analysis was used to evaluate the expression profile of cuproptosis-related genes and its predictive significance in esophageal cancer in this study.
MATERIALS AND METHODS
The study did not include any human participants, data, or tissue, nor did it include any animals. All of the information was gathered from a public database.
The Cancer Genome Atlas (https://portal.gdc.cancer.gov/) provided gene expression data and clinical information for 151 esophageal cancer samples, whereas the Genotype-Tissue Expression database (https://xenabrowser.net/) provided gene expression levels for 151 healthy tissue samples. The "limma" package in R was used to conduct matrix normalization. We selected 19 genes from the scientific literature that have been linked to cuproptosis in prior studies.
Screening of cuproptosis-related prognosis genes
In 151 patients with esophageal cancer and 151 healthy controls, a total of 55185 genes were acquired. The log2 (x+1) scale was used to standardize all of the expression data. Differentially expressed genes (DEGs) were screened between tumor tissue and normal esophageal tissue using the "limma" package and Wilcox-test, with the conditions of |log fold change [FC]| > 1 and false discovery rate (FDR) < 0.05. Using the "limma" package, we evaluated those genes that were differently expressed between esophageal cancer and healthy control tissues based on 19 cuproptosis-related genes. The FDR less than 0.05 was used as a criterion for further investigation. The researcher employed a univariate Cox regression analysis to find predictive markers linked to overall survival (OS) and progression-free survival (PFS). The optimum cutoff value for a prognostic gene was identified using a receiver operating characteristic (ROC) curve. Patients were classified into high-risk (less than the cutoff) and low-risk (greater than the cutoff) esophageal cancer groups based on gene expression cutoff values.
Function enrichment analysis
The R package "ClusterProfiler" was used to analyze gene ontology (GO) and pathway enrichment [Kyoto Encyclopedia of Genes and Genomes (KEGG)]. It was deemed substantially enriched when the P value and adjusted P value were both less than 0.05. To identify the functional role of the genes, GO analysis was done on the significantly expressed genes, and the expression levels of the genes were displayed in GO circle plots using the R package "GOplot". To determine tumor-infiltrating immune-cell fractions in esophageal cancer patients, we employed the CIBERSORT algorithm. On samples with a CIBERSORT result of P value less than 0.05, further analysis was performed. Then for reference, we collected 547 gene expression profiles from the CIBERSORT website (http://cibersort.stanford.edu/). The Pearson’s test and the "corrplot" program were used to correlate infiltrated immune cells.
Development of a nomogram for predicting OS
A cuproptosis-related gene and numerous clinicopathological variables were also subjected to univariate and multivariate Cox regression analyses. For predicting the OS, we utilized the R "rms" package to generate a nomogram including pathological categories, American Joint Committee on Cancer (AJCC)-tumor-node-metastasis (TNM) stages, and differential cuproptosis-related genes. The uniformity of the nomogram was assessed using calibration curves to anticipate various OS results. Harrell’s concordance index was used to generate the C-index. The area under the receiver operating characteristic (AUC) curve and ROC were both utilized to evaluate how predictive our nomogram was.
All statistical analyses were performed using R software, version 4.1.0 (http://www.rproject.org/). The Students t-test was used to compare variables such as age at diagnosis, sex, AJCC-TNM stage, mutational status of v-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog (KRAS), v-Raf murine sarcoma viral oncogene homolog B1 genes (BRAF), epidermal growth factor receptor (EGFR), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA, unpaired t-test). The Cox proportional hazards model was used to find genes and clinicopathological variables linked to survival. Based on the optimal cutoff of differential cuproptosis-related genes, the Youden index method determined that esophageal cancer patients should be divided into high- and low-risk groups. For assessing survival in the high- and low-risk groups, Kaplan-Meier survival curves were employed, and log-rank tests were performed to compare survival rates. For two-sided P values, statistical significance was considered as a value less than 0.05.
Demographics of the 151 esophagus cancer patients
Our study comprised 151 esophageal cancer patients and 151 healthy controls. Table 1 shows the characteristics of patients with EAC (n = 74) and ESCC (n = 77). All esophageal cancer patients had a median survival duration of 13.4 mo. EAC had a median survival time of 14.8 mo, whereas ESCC had a median survival time of 13 mo.
Table 1 Demographic characteristics and clinicopathological characteristics.
Relationship between clinical parameters and PFS/OS
Table 2 summarizes the relationships between clinical pathological features and these individuals’ OS or PFS. When compared to ESCC, EAC was related to a lower OS (P = 0.011). AJCC-TNM stages III and IV were linked to OS (P = 0.021) and PFS (P = 0.013). When compared to females, males had a lower PFS (P = 0.009). Age, EGFR, BRAF, KRAS, or PIK3CA status, among other clinicopathological variables, were not linked to OS or PFS.
Table 2 Relationship between the clinicopathological characteristics and survival outcome.
Identification of prognostic cuproptosis-related genes
With |log FC| > 1 and FDR < 0.05, a total of 7055 DEGs were identified in the healthy control group, comprising 3494 upregulated genes and 3561 downregulated genes. Then, using the FDR < 0.05 threshold, gene expression analysis was done to identify cuproptosis-related genes and 18 genes fulfilled our requirements (Supplementary Table 1). In the intersection of two gene sets, there were five genes: high-affinity copper uptake protein 1 (SLC31A1), ferredoxin, lipoyl(octanoyl) transferase 2, pyruvate dehydrogenase A1 (PDHA1), and programmed cell death 6 interacting protein (DEGs and cuproptosis-related genes). Furthermore, using univariate Cox regression analysis, we discovered that one gene (PDHA1) was strongly linked to prognosis [hazard ratio (HR) = 22.96, 95% confidence interval (CI) = 3.09-170.73, P = 0.002] (Table 3).
Table 3 The results of the univariate Cox regression analysis.
HR: Hazard ratio; HR.95L: Low 95% confidence interval of HR; HR.95H: High 95% confidence interval of HR.
Functional enrichment analysis
Using GO enrichment and KEGG pathway analyses, we explored the biological processes, cellular components, and molecular activities of DEGs. Energy metabolism and glycolysis/gluconeogenesis signaling pathways were highly enriched among DEGs with cuproptosis relevance. The phrases "tricarboxylic acid cycle," "acetylCoA metabolic process," and "acetylCoA biosynthetic process from pyruvate" are considerably enriched in Figure 1A. "Carbon metabolism," "tricarboxylic acid cycle," "pyruvate metabolism," "glycolysis/gluconeogenesis," "platinum drug resistance," "biosynthesis of cofactors," "mineral absorption," and "central carbon metabolism in cancer" were the most substantially enriched pathways (Figure 1B).
Figure 1 Esophageal cancer functional enrichment analysis of cuproptosis-related genes.
A: The expression of cuproptosis-related genes is represented by the outer circle in each enriched Gene Ontology (GO) terms; red dots on each GO term shows upregulated cuproptosis-related genes, and the inner circle indicates the importance of GO terms (log10 adjusted P value). The downregulated cuproptosis-related genes are shown by blue dots, and the chart on the right shows the distribution of cuproptosis-related genes in important GO terms; B: The top 16 Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis findings are shown. The GeneRatio denotes the number of differentially expressed genes (DEGs) detected in one GO pathway in proportion to the total number of DEGs.
Prognostic value of PDHA1
The 151 esophageal cancer patients were divided into two groups based on the cutoff value: Those with high PDHA1 expression and those with low PDHA1 expression. In the analysis of OS, high PDHA1 expression was linked to considerably lower OS rates than low expression (P = 0.007; Figure 2A). However, there was no link between PDHA1 expression and PFS (Figure 2B). Within a multivariate context, we investigated the link between PDHA1 expression and survival result (Table 4). When we controlled for clinical prognostic indicators that were significant (P < 0.05) in univariate Cox regression models, we discovered that overexpression of the PDHA1 gene might predict poor clinical outcomes. The expression of PDHA1 (HR: 1.67, 95 CI: 1.03-2.73; P = 0.0386) and AJCC-TNM stage (HR: 2.30, 95%CI: 1.58-3.35; P < 0.001) were independent risk factors for OS, but not for PFS, according to multivariate Cox regression analysis.
Figure 2 The Kaplan-Meier survival of patients with esophageal cancer who were classified according to pyruvate dehydrogenase A1 expression.
A and B: In esophageal cancer, patients with low pyruvate dehydrogenase A1 (PDHA1) expression levels had better overall survival (A) and progression-free survival (B) than patients with high PDHA1 expression levels. The significance of the difference between high and low expression was calculated using the log-rank test.
Table 4 Univariate and multivariate Cox regression analyses of pyruvate dehydrogenase A1 expression and the clinicopathological factors in the esophagus cancer patients.
Different infiltration levels of immune cells between low and high PDHA1 groups
Our study used the CIBERSORT algorithm to examine the ratio of tumor-infiltrating immune cells (TICs) in esophageal cancer to further verify the correlation between PDHA1 and immune cells (Figure 3A). P < 0.05 was used to classify samples as statistically different. The ratio of TICs in the low-risk group was represented by the first 108 of 151 esophageal cancer patients, whereas the ratio of TICs in the high-risk group was represented by the last 43 samples. The correlation between 22 immune cells was represented using a heatmap (Figure 3B). Neutrophils and activated mast cells, macrophages M1 and activated natural killer cells, plasma cells and naive B cells, activated mast cells and activated dendritic cells, and macrophages M1 and naive CD4+ T cells were the top five results with a positive correlation. CD8+ T cells and resting memory CD4+ T cells, on the other hand, were the immune cells that were most negatively associated. Furthermore, the bar graph revealed that the immune infiltration of 22 immune cell types in high and low PDHA1 esophageal cancer patients did not statistically significantly differ (Figure 3C).
Figure 3 Tumor-infiltrating immune cell profile.
A: A barplot shows the ratio of 22 different tumor-infiltrating immune cells in esophageal cancer patients. Column names (X-axis): sample ID; B: A heatmap shows the association between 22 different types of tumor-infiltrating immune cells. Each dot represents the P value of the correlation between two different cell types, and Pearson’s correlation coefficient was used to determine significance; C: A bar graph shows the difference between esophageal cancer patients with high or low expression of pyruvate dehydrogenase A1 in terms of 22 different types of tumor-infiltrating immune cells. The Wilcoxon rank-sum test was used for the significance test.
PDHA1 expression and immune checkpoint correlation analysis
We also investigated the relationships between PDHA1 and 47 genes associated with immunological checkpoints, which have been reported in the literature. The results showed that tumor necrosis factor superfamily member 15 (TNFSF15) (P = 3.2 × 10-6; r = 0.37), TNFRSF14 (P = 8.1 × 10-8; r = 0.42), H long terminal repeat-associating 2 (HHLA2) (P = 6.0 × 10-8; r = 0.42) and galectin 9 (LGALS9) (P = 3.1 × 10-6; r = 0.37) were significantly higher than that of the low expression group of PDHA1. Cluster of differentiation 44 (CD44) (P = 0.00028; R = -0.29), TNFRSF18 (P = 1.2 × 10-5; R = -0.35), programmed cell death 1 ligand 2 (PDCD1LG2) (P = 0.0032; R = -0.24), CD86 (P = 0.018; R = -0.19), and CD40 (P = 0.0047; R = -0.23) were significantly lower than that of the low expression group of PDHA1 (Figure 4A). The expression profiles of 47 immune checkpoint genes involved in cuproptosis were investigated further. In addition, we categorized the above genes into high and low expression groups based on their cutoff values. Kaplan-Meier analysis and log-rank testing were used to assess OS. The findings revealed that low HHLA2, TNFRSF18, and CD44 overexpression was linked to a significantly shorter OS and a worse prognosis (Figure 4B-D).
Figure 4 Correlation between pyruvate dehydrogenase A1 expression and immune checkpoint-related genes.
A: Barplot reveals expression of 47 immune checkpoint-related genes in high- and low-pyruvate dehydrogenase A1 (PDHA1) esophageal cancer patients; B-D: H long terminal repeat-associating 2 (B), tumor necrosis factor superfamily member 18 (C), and cluster of differentiation 44 (D) Kaplan-Meier survival curves of high and low expression of immune checkpoint-related genes. aP < 0.05; bP < 0.01; cP < 0.001.
Constructing and evaluating a predictive nomogram
The 1-, 2-, and 3-year OS probability was calculated using a nomogram that included the pathological type, AJCC-TNM stage, and PDHA1 expression (Figure 5A). The actual observed vs anticipated rates of the 1-, 2-, and 3-year OS demonstrated near to the ideal 45° oblique line on the calibration curve (Figure 5B). Furthermore, the AUC of ROC curves for 1-, 2-, and 3-year survival were 0.725, 0.776, 0.619, and 0.810, respectively (Figure 5C). With a C-index of 0.703 for OS, the nomogram showed promising discrimination.
Figure 5 Developing a predictive nomogram.
A: A nomogram that predicts patients with esophageal cancer’s 1-, 2-, and 3-year survival; B: A calibration curve of 1-, 2-, and 3-year survival in the nomogram and ideal model; C: Receiver operating characteristic curves and area under the curve for the nomogram’s 1, 2, and 3-year survival. AUC: Area under the curve; PDHA1: Pyruvate dehydrogenase A1.
In the present study, we investigated the expression signature of 19 cuproptosis-related genes in esophageal cancer tissues and explored their relationships with OS and PFS. A prognostic nomogram involving gene expression and clinicopathological parameters was constructed for the first time. Functional analysis showed that DEGs were enriched in energy metabolism, especially in pathways related to the tricarboxylic acid cycle. Cuproptosis-related genes were also confirmed to be associated with immune checkpoint genes.
In our study, 151 patients with esophageal cancer and 151 healthy controls were included. In univariate Cox regression analysis, PDHA1 was selected from 19 cuproptosis-related genes to be associated with the prognosis of esophageal cancer. At the same time, PDHA1 expression was also different in tumor tissues and healthy tissues. The results of the present study suggested that the PDHA1 expression was relatively low in cancer tissues, which was consistent with previous studies[16,17]. Li et al also reported that high expression of PDHA1 in ovarian cancer cells was significantly correlated with better OS and PFS. Our results suggested that the high expression of PDHA1 expression was associated with poor OS, but not with PFS.
Li et al showed that PDHA1 knockout inhibited glucose entering the tricarboxylic acid cycle, resulting in the reconnection of glutamine metabolism by increasing glutaminase 1 and glutamate dehydrogenase 1 expression, thus increasing the survival rate of glutamine-dependent cells. Consistent with our results of functional enrichment, these DEGs were significantly enriched in tricarboxylic acid cyclic-related pathways. We hypothesized that PDHA1 gene knockdown or low expression causes mitochondrial malfunction, resulting in aberrant generation of intracellular reactive oxygen species (ROS) and adenosine triphosphate, which is compatible with the cuproptosis mechanism.
TNM staging is linked to OS and PFS in esophageal cancer, whereas pathological types are associated with OS in esophageal cancer, according to our findings. Through univariate and multivariate Cox analyses, Zhang et al also claimed that pathological stage is an independent risk factor for OS and PFS in patients with operable esophageal cancer. Therefore, pathological types and TNM stages were included in the prognostic nomogram. Men and women had different rates of esophageal cancer, with men having nearly twice as many cases as women. We did not include sex in our nomogram because, while our analysis indicated that sex impacted PFS, the OS of esophageal cancer was not connected to sex.
CIBERSORT and ssGSEA techniques were used to analyze the composition of tumor-infiltrating immune cells (TICs) in each sample. Regulatory T cells and CD8+ T cells have been shown to play an important role in anti-tumor immunity in previous studies[21-23]. Studies have demonstrated that PDHA1 mediates metabolic reprogramming in macrophages[24,25]. However, there was no significant difference between the 22 types of immune cells in the high and low PDHA1 expression groups, according to our findings.
Immune checkpoint inhibitors have recently been studied in a variety of cancers, and they provide a novel therapy option. However, no evidence of a link between cuproptosis and immune checkpoint genes has been found. TNFSF15, TNFRSF14, HHLA2, LGALS9, CD44, TNFRSF18, PDCD1LG2, CD86, and CD40 all have a strong relationship with PDHA1 expression in our recent research. The link between immune checkpoint-related genes expression and prognosis in certain esophageal cancer patients is still debated[27,28]. The immune checkpoint-related genes associated with PDHA1 expression were analyzed using the log-rank test, and the results revealed that HHLA2, TNFRSF18, and CD44 were substantially correlated with prognosis.
Our study had a number of limitations. First, the sample size must be increased to analyze EAC and ESCC individually. Second, if the sample size is high enough, the treatment procedures and stages of the esophagus must be unified. Finally, given that prognostic characteristics were generated and analyzed using data from public databases, further biological evidence, in addition to the statistical evidence we present, was required.
This study analyzed the association between cuproprosis-related genes and the prognosis of esophageal cancer in a systematic method. Cuproprosis-related genes, especially PDHA1, can be used as prognostic predictors in esophageal cancer patients, providing additional information on how to treat the disease.
Despite many breakthroughs in treatment, the general prognosis for esophageal cancer, one of the least responsive malignancies to cancer therapy, remains poor. As a result, identifying biomarkers and understanding the molecular mechanisms of esophageal cancer were critical for improving patient outcomes.
A nomogram for predicting the prognosis of esophageal cancer would be developed by evaluating cuproprosis-related genes features and their correlation with prognosis in order to predict the prognosis of esophageal cancer.
Considering cuproprosis-related genes expression was linked to patient prognosis, we intended to develop a nomogram to predict prognosis based on cuproprosis-related genes characteristics and evaluate its prediction performance.
Cuproprosis-related genes were found to be linked with esophageal cancer prognosis using univariate COX regression analysis on 151 esophageal cancer samples. The C-index, calibration curve, and receiver operator characteristic (ROC) curve were used to evaluate the prediction ability of a prognostic nomogram created by combining clinicopathological variables and cuproprosis-related genes.
Univariate COX regression analysis of 19 Cuproprosis-related genes revealed that the expression of pyruvate dehydrogenase A1 (PDHA1) was associated with the prognosis of esophageal cancer. The low PDHA1 expression group had a better prognosis of esophageal cancer, according to the log-rank test. There was no statistical correlation between PDHA1 expression and 22 immune cells; however there was a correlation between PDHA1 expression and several immune checkpoint genes. The C-index, calibration curve, and ROC curve were used to confirm the predictive ability of the esophageal cancer prognostic nomogram, which was developed by combining pathological type, tumor-node-metastasis stage, and PDHA1 expression.
Cuproprosis-related genes are correlated to esophageal cancer prognosis, and a deep understanding of its molecular mechanism might contribute to novel cancer treatments in the future.
To enhance the overall survival of esophageal cancer patients, researchers must investigate cuproprosis biomarkers and anticipate possible therapy targets.
Tokunaga M, Okimoto K, Akizue N, Ishikawa K, Hirotsu Y, Amemiya K, Ota M, Matsusaka K, Nishimura M, Matsushita K, Ishikawa T, Nagashima A, Shiratori W, Kaneko T, Oura H, Kanayama K, Ohta Y, Taida T, Saito K, Matsumura T, Chiba T, Mochizuki H, Arai M, Kato J, Ikeda JI, Omata M, Kato N. Genetic profiles of Barrett's esophagus and esophageal adenocarcinoma in Japanese patients.Sci Rep. 2021;11:17671.
[PubMed] [DOI][Cited in This Article: ][Reference Citation Analysis (0)]
Li Y, Huang R, Li X, Yu D, Zhang M, Wen J, Goscinski MA, Trope CG, Nesland JM, Suo Z. Decreased expression of pyruvate dehydrogenase A1 predicts an unfavorable prognosis in ovarian carcinoma.Am J Cancer Res. 2016;6:2076-2087.
[PubMed] [DOI][Cited in This Article: ]
Saloura V, Izumchenko E, Zuo Z, Bao R, Korzinkin M, Ozerov I, Zhavoronkov A, Sidransky D, Bedi A, Hoque MO, Koeppen H, Keck MK, Khattri A, London N, Kotlov N, Fatima A, Vougiouklakis T, Nakamura Y, Lingen M, Agrawal N, Savage PA, Kron S, Kline J, Kowanetz M, Seiwert TY. Immune profiles in primary squamous cell carcinoma of the head and neck.Oral Oncol. 2019;96:77-88.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 36][Cited by in F6Publishing: 39][Article Influence: 12.0][Reference Citation Analysis (0)]
Semba H, Takeda N, Isagawa T, Sugiura Y, Honda K, Wake M, Miyazawa H, Yamaguchi Y, Miura M, Jenkins DM, Choi H, Kim JW, Asagiri M, Cowburn AS, Abe H, Soma K, Koyama K, Katoh M, Sayama K, Goda N, Johnson RS, Manabe I, Nagai R, Komuro I. HIF-1α-PDK1 axis-induced active glycolysis plays an essential role in macrophage migratory capacity.Nat Commun. 2016;7:11635.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 171][Cited by in F6Publishing: 173][Article Influence: 28.5][Reference Citation Analysis (0)]
Samstein RM, Lee CH, Shoushtari AN, Hellmann MD, Shen R, Janjigian YY, Barron DA, Zehir A, Jordan EJ, Omuro A, Kaley TJ, Kendall SM, Motzer RJ, Hakimi AA, Voss MH, Russo P, Rosenberg J, Iyer G, Bochner BH, Bajorin DF, Al-Ahmadie HA, Chaft JE, Rudin CM, Riely GJ, Baxi S, Ho AL, Wong RJ, Pfister DG, Wolchok JD, Barker CA, Gutin PH, Brennan CW, Tabar V, Mellinghoff IK, DeAngelis LM, Ariyan CE, Lee N, Tap WD, Gounder MM, D'Angelo SP, Saltz L, Stadler ZK, Scher HI, Baselga J, Razavi P, Klebanoff CA, Yaeger R, Segal NH, Ku GY, DeMatteo RP, Ladanyi M, Rizvi NA, Berger MF, Riaz N, Solit DB, Chan TA, Morris LGT. Tumor mutational load predicts survival after immunotherapy across multiple cancer types.Nat Genet. 2019;51:202-206.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 1654][Cited by in F6Publishing: 1764][Article Influence: 551.3][Reference Citation Analysis (0)]