Yuzhi Wang
Yunfei Xu2†
Yi Zhang- 1Department of Laboratory Medicine, Deyang People’s Hospital, Deyang, Sichuan, China
- 2Department of Laboratory Medicine, Chengdu Women’s and Children’s Central Hospital, Chengdu, Sichuan, China
- 3Department of Ultrasound, The First People’s Hospital of Yibin, Yibin, Sichuan, China
- 4Department of Blood Transfusion, Deyang People’s Hospital, Deyang, Sichuan, China
Background: Disulfidptosis is a newly identified variant of cell death characterized by disulfide accumulation, which is independent of ATP depletion. Accordingly, the latent influence of disulfidptosis on the prognosis of lung adenocarcinoma (LUAD) patients and the progression of tumors remains poorly understood.
Methods: We conducted a multifaceted analysis of the transcriptional and genetic modifications in disulfidptosis regulators (DRs) specific to LUAD, followed by an evaluation of their expression configurations to define DR clusters. Harnessing the differentially expressed genes (DEGs) identified from these clusters, we formulated an optimal predictive model by amalgamating 10 distinct machine learning algorithms across 101 unique combinations to compute the disulfidptosis score (DS). Patients were subsequently stratified into high and low DS cohorts based on median DS values. We then performed an exhaustive comparison between these cohorts, focusing on somatic mutations, clinical attributes, tumor microenvironment, and treatment responsiveness. Finally, we empirically validated the biological implications of a critical gene, KYNU, through assays in LUAD cell lines.
Results: We identified two DR clusters and there were great differences in overall survival (OS) and tumor microenvironment. We selected the "Least Absolute Shrinkage and Selection Operator (LASSO) + Random Survival Forest (RFS)" algorithm to develop a DS based on the average C-index across different cohorts. Our model effectively stratified LUAD patients into high- and low-DS subgroups, with this latter demonstrating superior OS, a reduced mutational landscape, enhanced immune status, and increased sensitivity to immunotherapy. Notably, the predictive accuracy of DS outperformed the published LUAD signature and clinical features. Finally, we validated the DS expression using clinical samples and found that inhibiting KYNU suppressed LUAD cells proliferation, invasiveness, and migration in vitro.
Conclusions: The DR-based scoring system that we developed enabled accurate prognostic stratification of LUAD patients and provides important insights into the molecular mechanisms and treatment strategies for LUAD.
Introduction
Lung cancer is a global health concern and one of the foremost sources of cancer patient morbidity and mortality, posing a grave threat to public health (1). In 2020, an estimated 2,206,771 new cases were diagnosed, and 1,796,144 fatalities occurred (2). Non-small cell lung cancer (NSCLC) stands as the prevalent pathological category among lung cancers, comprising about 85% of all cases, whereas lung adenocarcinoma (LUAD) is the most frequently occurring subtype of NSCLC (3). LUAD is plagued by a high incidence of invasive behavior and metastatic spread. However, significant advances in the treatment of LUAD have been achieved in the past few decades. Accordingly, targeted therapies and immunotherapies have been demonstrated to enhance the treatment efficacy and outcomes for patients with LUAD (4). However, the majority of patients develop drug resistance and relapse following initial treatment, resulting in no significant improvement in 5-year survival rates (5). Thus, to benefit more LUAD patients, there is an urgent need to identify new therapeutic targets and prognostic indicators for predicting survival and guiding clinical treatment in LUAD patients.
Programmed cell death (PCD) is a cellular death process governed by a molecular program that is controlled by particular genes. This process is essential for the normal development of organisms and the preservation of homeostasis (6). Exploration and characterization of these cell death mechanisms not only deepens our fundamental comprehension of cellular equilibrium, but also offers valuable perspectives for the therapeutic approach to a variety of diseases, including cancer. For instance, recent studies have progressively highlighted the tumor-inhibiting effects of ferroptosis, which is achieved through the deprivation of cysteine and the generation of reactive oxygen species (ROS) by p53 (7, 8). Moreover, the enhancement of lipid peroxidation by activated CD8+ T cells can induce ferroptosis, which contributes to the antitumor efficacy of immunotherapy (9). The role of autophagy in cancer can vary, depending on the specific tumor model and tumor stage. During the initial stages of cancer, autophagy functions as a protective mechanism, shielding normal cells from tumorigenesis by preventing DNA damage and mutations (10). However, in the context of fully-formed solid tumors, autophagy shifts its role and promotes tumor progression by promoting tumor growth, enhancing cell survival, enhancing resistance to platinum-based chemotherapy, and facilitating the formation of metastases (11). Autophagy inhibitors, therefore, comprise one of the treatment options for patients with advanced tumors. Recent research (Liu et al., 2023) demonstrates that excessive accumulation of disulfide induces a unique form of controlled cell death known as “disulfidptosis” that is distinct from apoptosis, necrosis, autophagy, and ferroptosis (12). SLC7A11 is an essential transport protein whose primary function is to facilitate the cellular uptake of cysteine (13). Cysteine is a necessary building block for the synthesis of glutathione and a crucial component for inhibiting oxidative stress in cells and regulating iron death pathways (14). However, it has also been shown to possess certain cytotoxic properties (15). Gan Bo et al., discovered that under conditions of glucose deprivation, high expression of SLC7A11 leads to a significant consumption of NADPH within cells, abnormal aggregation of disulfides like cysteine, inducing disulfide stress and rapid cell death (16). This form of cell death induced by glucose deprivation and high SLC7A11 expression in cancer cells cannot be prevented by inhibitors of cell death that act on other cells, nor is it caused by depletion of intracellular ATP. However, thiol-oxidizing agents, such as Diamide, can enhance this effect. Moreover, under glucose-deficient conditions, the number of disulfide bonds in the actin cytoskeleton increases significantly. Therefore, this study suggests that the induction of disulfide-dependent cell death by GLUT inhibitors may be an effective cancer treatment strategy. Consequently, focusing on disulfidptosis regulators (DRs) as potential targets provides new perspectives for understanding the complexities of the occurrence and development mechanism in LUAD. This approach is, therefore, of significant importance in enhancing the efficacy of treatment in patients with LUAD patients. However, the full scope of the impact of DRs on outcomes and treatments for LUAD patients has yet to be comprehensively explored.
In the current study, we divided 1569 LUAD samples into two disulfidptosis-associated subtypes according to 18 DRs and compared survival and immune infiltration between the subtypes. We also developed a disulfidptosis score (DS) to predict overall survival (OS) and to delineate the immunological landscape of LUAD. As indicated by the findings, a higher DS was associated with unfavorable prognostic outcomes and worse immunotherapy responses in LUAD, suggesting the potential clinical utility of DS as a tool for assessing prognosis and immunotherapy efficacy. Thus, the current study introduces an innovative approach for assessing the efficacy of immunotherapy and predicting the prognosis of LUAD patients based on DS.
Materials and methods
Data collection
The LUAD data were obtained from The Cancer Genome Atlas (TCGA)-LUAD (https://portal.gdc.cancer.gov/), GSE31210, GSE68465, and GSE72094 in Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/), which included RNA sequencing data, somatic mutation data, copy number variation (CNV) data, and corresponding clinical information. The TCGA-LUAD cohort, which consisted of 539 LUAD tissues and 59 cancer-adjacent tissue samples, served as a training cohort. Meanwhile, GSE31210, GSE68465, and GSE72094 were utilized for the validation cohort, which consisted of 1,150 LUAD patients. Additionally, the “sva” package was used to correct the batch effect between the different datasets by adopting the “combat” algorithm (17). Moreover, patients who lacked OS time were filtered out. Finally, 1569 eligible patients were encompassed in the study. The detailed clinical characteristics of all LUAD patients is listed in Supplementary Table 1. Additionally, 18 DRs were collected from the previous study (Supplementary Table 2) (12).
Consensus cluster analysis of DRs
A consensus clustering algorithm was employed to discern optimal subtypes founded on the expression of 18 DRs using the R package “ConsensusClusterPlus”. The number of clusters (K) and their stability (with 1,000 repeats for mast k = 9) were determined by the consensus clustering algorithm (18). The clustering was based on dividing centromeres with “Euclidean” distances (the most common and familiar distance measurement methods and correlation of K-Means clustering). Additionally, a T-distributed Stochastic Neighbor Embedding (tSNE) analysis was conducted to decrease the dimensions and differentiate the subtypes of information (tSNE can preserve local similarities between data points and is one of the most used unsupervised clustering visualization methods).
Differentially expressed genes and functional annotation
DEGs among DRs subtypes were determined using the “limma” R package with the filtering criteria of log (fold change) >1 and False-discovery rate (FDR) <0.05 (19). To further investigate the potential functions and enriched pathways of DEGs, functional enrichment analyses were conducted on DEGs employing the “clusterprofiler” R package (20).
Generation of DS
The DS was calculated to quantify the disulfidptosis patterns of the LUAD. Accordingly, a univariate cox regression analysis was conducted to identify the DEGs are related to prognosis. Subsequently, the patients were segregated into distinct gene cluster groups according to the expression of prognostic DEGs using the unsupervised clustering method. Based on prognostic DEGs, the 10 machine learning algorithms, including ‘Least Absolute Shrinkage and Selection Operator (LASSO, “glmnet” package)’, ‘Ridge (“glmnet” package)’, ‘Elastic network (“glmnet” package)’, ‘StepCox (“survival” package)’, ‘Survival support vector machine (survival-SVM, “survivalsvm” package)’, ‘CoxBoost (“CoxBoost” package)’ (21), ‘Supervised principal components, “superpc” package’ (22), ‘partial least squares regression for COX, “plsRcox” package’, ‘random survival forest (RSF, “randomForestSRC” package)’, ‘generalized boosted regression modeling, “gbm” package’ were used to constructed the models. Briefly, 101 combinations of 10 machine learning algorithms were used to build the models based on a leave-one-out cross-validation (LOOCV) framework. Models with <3 genes were excluded. Simultaneously, each patient’s linear score and concordance index (C-index) are calculated based on various models, and the optimal model is selected based on the average highest C-index in the training and testing cohorts. Utilizing the optimal model, the DS for each patient was determined. See the supplementary methods table for details. The patients were then separated into high- and low-DS groups using the median DS value. Subsequently, Kaplan–Meier survival analysis was utilized to compare the OS rates of patients in various DS groups.
Somatic mutation and CNV analysis
The “maftools” package was utilized to evaluate and visualize the mutation type and frequency of the genes (23). Correspondingly, the tumor mutation burden (TMB) of each LUAD sample was calculated based on the total count of somatic mutations per megabase (MB) in the exon coding region of the human genome. Different mutation types were classified as either synonymous or nonsynonymous mutations. The nonsynonymous variants included Frame_Shift_Del, Frame_Shift_Ins, In_Frame_Del, In_Frame_Ins, Missense, Nonsense, Nonstop, Splice_Site, and Translation_Start_Site. The maftools analysis focused on identifying significantly mutated genes (P <0.05) between the two groups and assessing the interaction effect of gene mutations. Only genes with mutations occurring 30 times or more in at least one group were considered for both analyses. GISTIC 2.0 was used to identify significant regions within CNV data (24). To quantify and compare CNVs, we calculated the fraction of altered genome (FGA), fraction of genome gained (FGG), and fraction of genome lost (FGL) for each sample. FGA reveals the proportion of genomic segments that have been altered. FGG/FGL considers specifically the genomic segments that have undergone gain or loss, respectively.
Tumor microenvironment landscape and hallmark pathway analyses
TME at LUAD were evaluated under four aspects. First, the immune score, stromal score, ESTIMATE score, and tumor purity were calculated using the ESTIMATE algorithm (25). Secondly, three algorithms, Single Sample Gene Set Enrichment Analysis (ssGSEA), Tumor Immune Estimation Resource (TIMER, https://cistrome.shinyapps.io/timer/), and “MCPcounter” were used to quantify the relative infiltration of immune cells in the entire cohort (26–28). Thirdly, the seven steps of cancer immunity cycle were analyzed using the Tracking Tumor Immunophenotype (TIP) website (http://biocc.hrbmu.edu.cn/TIP/) (29). In the fourth step, 35 inhibitory immune checkpoints were extracted from a prior study. Subsequently, Gene Set Enrichment Analysis (GSEA) was utilized to identify underlying mechanisms in hallmark gene sets with the recommended criteria (FDR <0.25 and NES >1), in order to determine the underlying hallmark pathways associated with DS (30). Meanwhile, “Gene Set Variation Analysis” (GSVA) package was applied to the two DS groups with an adjusted p value <0.01 (31). The “h.all.v7.4.symbols.gmt” hallmark gene sets from the MSigDB database were employed for GSVA implementation.
Assessment of immunotherapy and chemotherapy
To explore the predictive value of DS in LUAD patients after immunotherapy, we compared the immune dysfunction and exclusion (TIDE, http://tide.dfci.harvard.edu/) score; additionally, the submap algorithm was applied to compare the efficacy of immunotherapy among various DS subtypes (32, 33). In addition, four immunotherapy-treated cohorts, IMvigor210, GSE35640, GSE79671, and GSE173839, were collected to investigate the immunotherapy response ability of DS. The sensitivity of tumor cell lines to potential drugs was obtained from the Cancer Therapeutics Response Portal (CTRP, https://portals.broadinstitute.org/ctrp) and Profiling Relative Inhibition Simultaneously in Mixtures (PRISM, https://depmap.org/portal/prism/). The more sensitive a cell line is to a potential drug, the lower its area under the curve (AUC).
Single-cell RNA-sequencing analysis
ScRNA-seq data was downloaded from GSE131907. Subsequently, the “Seurat” R package was utilized to measure the gene expression levels by processing the raw data from each sample (34). Cells with fewer than 200 detected genes were removed. Accordingly, the top 2000 highly variable genes were selected for subsequent clustering analysis. Following this, single cells were classified into distinct subgroups via the application of the FindNeighbors and FindClusters functions (dim = 15 and resolution =0.2). In addition, the tSNE was constructed utilizing the top 15 primary components. Subsequently, immune cells and tumor cells were identified using the “single R” and “copyKAT” packages (35, 36). Pseudotime trajectory analysis is a method used to study the temporal order of cells during development, differentiation, or other biological processes. Its goal is to reduce the dimensionality of single-cell data from high-dimensional to one-dimensional, thus representing the temporal changes of cells in a pseudotime manner. This aids in understanding the timeline of processes like differentiation, development, and transformation in tumor cells. Hence, the cell trajectory cancer cell populations were ordered in pseudotime using the “Monocle” package (37).
Cell culture and transfection
Two human LUAD cells (PC-9 and H838) as well as a normal bronchial epithelial cell (BEAS-2B) were obtained from the American Type Culture Collection (Manassas, VA, USA). The PC-9 and H838 cells were grown in RPMI 1640 supplemented with 10% FBS and 1% penicillin-streptomycin, while BEAS-2B cells were grown in DMEM supplemented with 10% FBS and 1% penicillin-streptomycin in humidified air at 37°C and 5% CO2 (This condition simulates the physiological environment in the human body, aiding in maintaining normal cellular growth, metabolism, and function). To construct KYNU knockdown and negative control (NC) cell lines, H838 and PC-9 cells were seeded in 6-well plates at a density of 5 × 104 cells/well and transfected with 50 nM siRNAs-KYNU and siRNA negative control (siRNA-NC) using lipofectamine 3000 (This reagent possesses high efficiency, broad spectrum, and low toxicity) following the manufacturer’s guidelines (Hanheng, Shanghai, China). After 48 h of transfection, subsequent experiments were conducted (Typically, after 48 hours of transfecting siRNA, the target gene’s expression is effectively disrupted in the vast majority of cases, while cells continue to maintain a relatively healthy and appropriate growth state). The siRNA sequences are provided in Supplementary Table 3.
Tissue microarray and immunohistochemistry
For this study, a total of 15 TM tissue samples (HLugA030PG04-1) were utilized, including 15 LUAD and 15 adjacent non-tumor tissues. All the tissues were procured from Shanghai Outdo Biotech Co., Ltd. (Shanghai, China). TM was stained using IHC with KYNU. In brief, the antibodies were diluted to the suitable concentration and incubated overnight with the sections at 4°C. The avidin–biotin and streptavidin complex were then incubated with the biotinylated goat anti-rabbit IgG secondary antibody for 0.5 h. The cell nuclei were counterstained blue with Hematoxylin. Each specimen was graded based on the intensity (0: absent, 1: mild, 2: moderate, and 3: pronounced) and the proportion of positively stained cells (0: 0%, 1: 1%–25%, 2: 26%–50%, 3: 51%–75%, and 4: 76%–100%). The final IHC scores were then calculated as the product of the intensity score and the percentage score (38).
RNA-isolation and real-time reverse transcription polymerase chain reaction
RNA extraction was performed employing RNA extraction kit (AG, Changsha, China) in accordance with the manufacturer’s instructions. The complementary DNA (cDNA) was synthesized from total RNA for each sample by using a reverse transcription kit (Promega, Madison, Wisconsin). Subsequently, the mRNA expression was quantitatively analyzed via RT-PCR (Roche Light Cycler® 480II System). The relative gene expression was normalized to β-actin, acting as the control gene. The primer sequences for the genes are outlined in Supplementary Table 4. The fold change of the target gene was computed utilizing the 2−ΔΔCt method.
Proliferation assay
Cells in logarithmic phase were seeded into 96-well plates at density of 2,000 cells per 100 µl. They were then cultured for 0, 24, 48, 72, and 96 h. At the end of the incubation cycle, CCK-8 reagent (Saiku, Shanghai, China) was added to each well for 2 h at 37°C. Subsequently, the absorbance values of each well were measured at 450 nm. For colony formation assay, 500 cells were seeded into each well of a 12-well plate and incubated. After 10 days, the cells were stained with crystal violet solution, and fixed with 4% paraformaldehyde for 0.5 h. Finally, colonies consisting of more than 50 cells were counted.
Migration and invasion assay
The transwell experiment was executed using a transwell chamber (Transwell, Corning Costar, USA). For migration assays, a total of 4 × 104 cells were resuspended in 200 μl serum-free medium and placed to the upper chamber. The lower chamber was filled with 700 μl of medium containing 10% FBS. After 18 h, cells that migrated through the membrane were fixed and stained with hematoxylin. Following a 24 h-incubation period, cells that traversed the membrane were fixed with 4% paraformaldehyde for 30 minutes and stained with a crystal violet solution for 15 min. Subsequently, they were observed with a microscope, and five random fields of view were selected to count the cells. The protocol for the invasion assay was similar to the migration assay, with the exception that the upper chambers were coated with 70 μl of diluted Matrigel. In addition, the wound healing assay was utilized to examine cell migration. Serum-starved LUAD cells (2 × 105) were seeded in 6 well plates and then transfected with siRNA-NC or siRNA-KYNU. Upon reaching >95% cell confluence, a sterile 20 μl pipette tip was used to create a scratch in the monolayer of cells. After a PBS wash, the cells were incubated for 0 and 24 h in the corresponding basic culture medium. The monolayer cells were examined microscopically, and the gap distance was measured quantitatively to ascertain LUAD cell migration.
Statistical analysis
All statistical analyses and representations were conducted using R (version 4.2.1) and GraphPad Prism (version 9.00). The Chi-squared test was used to compare the proportion of individuals within two groups. Additionally, continuous variables in two or more groups were compared using Wilcoxon rank-sum test or Kruskal–Walli’s test. For correlation analysis, the Pearson correlation test was employed. A p-value less than 0.05 was deemed statistically significant.
Results
Genetic and transcriptional alterations of DRs in LUAD
Gene expression analysis using bulk RNA-seq demonstrated that the majority of DRs displayed relative higher expression levels in LUAD tissues compared to para-carcinoma tissues (Figure 1A). As depicted in Figure 1B, 104 of 616 (16.88%) LUAD samples possessed genetic mutations. A total of 13 of 18 DRs were found to be mutated, with CNOT1 exhibiting the highest rate of mutation. Among them, missense mutations were found to be the most frequent (Figure 1B). In order to unmask the genetic modifications in DRs, we presented an overview of the frequency of somatic and copy number mutations with malignancies. Moreover, analysis of these 18 DRs revealed that CNV alterations were common. NDUFS2, NUBPL, PPM1F, EPAS1, and LRPPRC displayed widespread CNV amplification, whereas CCNC, NDUFA11, OXSM, and GYS1 showed widespread CNV deletions (Figure 1C). Figure 1D depicts the locations of CNV alterations in LUAD DRs. The correlation network composed of 18 DRs is illustrated in Figure 1E (correlation coefficient >0.4; the positive correlation is represented by the red line, negative correlation is represented by the blue line).
Figure 1 The genomic features and expression of DRs in LUAD. (A) The differential expression of DRs between tumor and normal samples. (B) Mutation landscape of DRs in TCGA-LUAD. (C) The CNV mutation frequency of DRs. (D) Chromosome position and alteration of DRs. (E) Molecular interaction network map of DRs Negative correlations are illustrated in green, while positive correlations are denoted in pink. ns, not significant, *P < 0.05, **P < 0.01, ***P < 0.001.
Identification of DRs subtypes in the LUAD cohort
In order to delve further into understanding the DRs’ expression pattern in the oncogenesis of LUAD, 1569 patients from 4 independent LUAD cohorts were included. Accordingly, we performed unsupervised clustering and classification on the combined LUAD cohort based on DRs. Our results showed that k = 2 was the optimal choice (Figures S1). In addition, the tSNE results revealed significant differences between the two clusters in terms of DRs expression (Figure 2A). Moreover, the Kaplan–Meier survival analysis demonstrated that the DRcluster A had a greater survival advantage than DRcluster B (Figure 2B). In addition, the clinicopathological characteristics of the various DR clusters also revealed significant differences (Figure 2C). To further explore the biological behavioral difference between these two clusters, we also conducted a GSVA enrichment analysis (Figure 2D). The results demonstrated that DRcluster A was primarily enriched for carcinogenic pathways like focal adhesion, EMC receptor interaction, and others. Moreover, Figure 2E data also revealed significant differences in the relative expression of immune infiltration cells across two DRclusters.
Figure 2 The construction of DRclusters. (A) tSNE plot of two DRclusters. (B) Kaplan-Meier survival analysis between two DRclusters. (C) Heatmaps of the distribution of DRs in the two DRclusters. (D) GSEA analysis indicating significant enrichment of pathways in the two DRclusters. (E) The proportion of 24 kinds of immune cells in two DRclusters. ns, not significant, *P < 0.05, **P < 0.01, ***P < 0.001.
Identification of genes subtypes and establishment of DS
A total of 51 DEGs were identified from two DRclusters using the “limma” package. These DRs subtypes-related genes were significantly enriched in cellular metabolism (Figures 3A, B). Subsequently, we conducted a univariate Cox regression analysis for DEGs and identified 39 DEGs with prognostic importance for LUAD. To further investigate the heterogeneity of DRs subtypes, we performed unsupervised clustering on 39 DEGs. The 1569 patients with LUAD were also separated into two gene clusters, designated geneCluster A and B (Figure S2). Similar to DRcluster, these genes accurately differentiate LUAD patients, with distinct clusters of genes exhibiting variations in DRs, survival rates, and immune cell infiltration (Figures 3C–F). Following this, the LOOCV framework was used to fit 101 prediction models to both training and testing sets. The C-index for each model was calculated, and based on the model with the highest average C-index (0.713), “Lasso+RSF” was deemed the optimal model (Figures 4A–C, Supplementary Table 5). A DS was calculated for each patient based on the expression of 7 genes (KRT6A, NEIL3, KYNU, ABCC2, SFTPC, CPS1, and INSL4) weighted by their regression coefficients in the model, and patients were divided to high- or low‐DS groups based on the median cutoff point of DS (Supplementary Table 6). As evident from the K–M survival analysis, OS rates were significantly diminished in the high-DS group compared to the low-DS group (Figures 4D–G). Moreover, the relationship between different types of patients and their prognoses was analyzed (Figure 5), with results suggesting that a low DS was related to a better prognosis in all patient categories.
Figure 3 Identification of DRs gene clusters based on DEGs in DRclusters. (A, B) GO and KEGG enrichment analyses of DEGs among two DRclusters(C) Differences in the expression of DRs among the two genecluster (D) tSNE plot of two geneclusters. (E) Kaplan-Meier survival analysis between two geneclusters. (F) The proportion of 24 kinds of immune cells in two geneclusters. *p < 0.05; ***p < 0.001; ns, not significant.
Figure 4 Construction of prognostic signature based on DEGs. (A) A total of 101 kinds of prediction models via a leave-one-out cross-validation framework and further calculated the C-index of each model. (B) Cvfit and lambda curves of LASSO regression applied with minimum criteria. (C) The number of trees determined by minimum error and importance of the four most valuable genes based on the RSF algorithm. (D-G) Kaplan-Meier survival curves of OS for high- and low-DS groups of patients in the TCGA, GSE31210, GSE68465, GSE72094 cohorts, respectively.
Figure 5 Correlation analysis of DS. (A) Alluvial diagram of clusters distributions in groups with different DS and survival outcomes (B) Expression of DRs between high- and low-DS groups. (C) Differences in DS between DRclusters. (D) Differences in DS between gene clusters. (E) The circular pie chart for the proportion difference of clinical indices. *P < 0.05, ***P < 0.001.
Evaluation of the DS
A time-dependent receiver operating characteristic (ROC) curve was employed to assess the validity of DS, as the AUC value for TCGA (0.93–0.96), GSE31210 (0.68–0.79), GSE68465 (0.63–0.66), GSE72094 (0.63–0.72), and meta-cohort (0.71–0.78) (Figures 6A–E). In addition, the C-index of clinical factors in patients with LUAD was determined (Figures 6F–J). Notably, the DS had a higher predictive efficacy than the vast majority of clinical indicators. Subsequently, both DS and clinical indicators were subjected to univariate and multivariate Cox analyses. In all cohorts, the DS was determined to be an independent indicator of OS prognosis (Tables 1–4). In order to determine the prognostic efficacy of DS, we combined 56 previously published LUAD prognostic models and conducted a comparative analysis of each model’s C-index. These models were developed using a variety of biologically relevant features, including autophagy, EMT, ferroptosis, hypoxia, necroptosis, glycolysis, and m6A methylation. Accordingly, DS was found to exhibit superior performance relative to the vast majority of models across all cohorts (Figure 7, Supplementary Table 7). Cumulatively, these results, therefore, demonstrate that the DS would be a valuable LUAD prognostic model.
Figure 6 Evaluation of the DS. (A-E) Time-dependent ROC curves presented with the 1-5 year AUC in TCGA, GSE31210, GSE68465, GSE72094 and meta-cohort. (F) The C-index of the CDS for the TCGA, GSE31210, GSE68465, GSE72094 cohorts. (G-J) The C-index of the DS and other clinical factors in the TCGA, GSE31210, GSE68465, GSE72094 cohorts. *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; ns, not significant.
Table 1 Univariate and multivariate Cox analysis of the clinicopathological features and FA score with OS for TCGA cohort.
Table 2 Univariate and multivariate Cox analysis of the clinicopathological features and FA score with OS for GSE68465 cohort.
Table 3 Univariate and multivariate Cox analysis of the clinicopathological features and FA score with OS for GSE31210 cohort.
Table 4 Univariate and multivariate Cox analysis of the clinicopathological features and FA score with OS for GSE72094 cohort.
Figure 7 C-index analysis between the DS and 56 published signatures in TCGA, GSE31210, GSE68465, GSE72094 and meta-cohort. *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; ns, not significant.
Comparison of the mutations and CNV between DS groups
Using the “maftool” package, a comparison was made between the distribution differences of somatic mutations observed in high- and low-DS groups. (Figures 8A, B). Comparing the frequency of mutants between the high- and low-DS groups, more somatic mutations, both synonymous and nonsynonymous, were observed in the high-DS group (Figures 8C–E). In addition, maftools analysis results showed that 17 genes mutated more frequently in LUAD patients in the high-DS group, including KEAP1, STK11, SMARCA4, TTN, NCKAP5, COL5A2, ANKRD30A, TEX15, PCDH15, GRIN2B, AHNAK, FAT4, FMN2, FAT1, ZNF804B, DOCK2 and COL22A1 (Figure 8F), and there was extensive co-mutation between these genes (Figure 8G). This is consistent with the mutation analysis described above; accordingly, TMB was found to be higher in the high-DS group compared to the low-DS group (Figure 8H). Subsequently, LUAD patients were classified into two mutation groups based on their TMB score. When combining DS and TMB, we discovered that patients with low TMB from the high-DS group had the worst prognosis (Figure 8I). Subsequently, we used the GISTIC 2.0 software to decipher the amplification and deletion of CNA on chromosome. Compared to the high-DS group, the low-DS group had a greater burden of amplification and deletion at both the arm and focal levels (Figure 8J). However, no significant disparities were observed between the high- and low-DS groups in terms of FGA, FGG, and FGL (Figure 8K). These results indicate a certain degree of correlation between DS and mutations.
Figure 8 Integrated comparisons of somatic mutation and CNVs between high- and low-DS groups in the TCGA cohort. (A, B) Visual summary showing common genetic alterations in the high and low-DS groups. (C-E) Association between all mutation counts, synonymous mutation counts, nonsynonymous mutation counts, and DS and their distribution in the DS groups. (F) Forest plot of gene mutations in the patients. (G) Interaction effect of genes mutating differentially in patients. (H) Tumor mutation burden between high- and low-DS groups. (I) Comprehensive survival analysis based on DS and TMB. (J) Gene fragments profiles with amplification (red) and deletion (green) among the DS groups. (K) Comparison of the fraction of the genome altered, lost, and gained between the DS groups. *P < 0.05, **P < 0.01, ***P < 0.001.
TME analysis
In order to evaluate the discriminative potential of the DS subgroup for the TME and its applicability in immunotherapy, we simultaneously evaluated the abundance of immune cell infiltration across multiple samples using four distinct algorithms. Unsurprisingly, as the DS increased, the number of immune cells declined (Figure 9A). Apparently, the activation of key steps in the cancer immunity cycle, such as step 3 (priming_and_activation) and step 4 (CD4 T cell recruiting, Dendritic cell recruiting, Macrophage recruiting, Monocyte recruiting, and T cell recruiting) appeared to be significantly higher in the low-DS group than in the high-DS group (Figure 9B). Subsequently, the expression profile of immune checkpoints in the two DS groups were further evaluated. Accordingly, the analysis revealed that the low-DS group demonstrated elevated expression levels of immune checkpoints, including HHLA2 and CD48 (Figure 9C). Given the observed upregulation of immune-related characteristics in the group with low DS, its underlying biological mechanisms were investigated further. As evident from the findings, DS exhibited correlations with multiple metabolic pathways (Figure 9D). In addition, DS demonstrated a strong correlation with numerous immunotherapeutic strategies (Figure 9D). To investigate the cancer signaling pathways associated with DS, GSVA analysis was performed on high- and low-DS groups. Using a predetermined threshold, we discovered that 16 signature pathways were significantly upregulated in the high-DS group compared to the low-DS group (Figure 9E). GSEA validated that 12 of these pathways were upregulated in high-DS patients, the majority of which are known to be carcinogenic (Figures 9F–I).
Figure 9 Analysis of the TME in different DS groups. (A) Differences in immune infiltration status between two DS groups were evaluated by four algorithms. (B) The differences of cancer immunity cycle were showed in boxplot between two DS groups. (C) The differences of immune checkpoint related genes were showed in boxplot between two DS groups. (D) The correlations between the TIIClnc signature score and metabolic immune-related pathways, immune-related pathways based on GSVA of GO and KEGG terms were displayed in butterfly plot. (E) The difference in the hallmark gene sets between different DS groups based on GSVA. (F-I) The GSEA results for the 12 overlapping upregulated hallmark pathways in terms of the high-DS group. *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; ns, not significant.
Assessment of immunotherapy and chemotherapy
The potential of an immunotherapy response was subsequently predicted for each immune cluster utilizing the TIDE algorithm and submap analysis. Accordingly, lower TIDE scores were observed in the low-DS group, implying a higher sensitivity to immunotherapy in these patients (Figure 10A). Moreover, the submap results indicated that the group with a low DS level was more sensitive to CTLA4 inhibitors (Figure 10B). Although we evaluated an individual’s immunotherapy efficacy using two algorithms, it remains critical to directly compare the curative efficacy of immunotherapy cohorts across various DS groups. As a result, four immunotherapy cohorts were included for further analysis. In the IMvigor210 cohort, patients with DR/PR had significantly longer OS compared to patients with SD/PD, whereas the influence of DS on patient prognosis was minimal (Figure 10C). However, patients who responded better to immunotherapy had lower DS levels across all cohorts (Figures 10D–G). Although there were no significant differences in patient survival and DS between the two groups for some cohorts, the propensity for these results was consistent for the other cohorts. In addition, these results demonstrate that the DS was able to predict the efficacy of ICBs and can provide direction for the deployment of immunotherapy. Thus, to pinpoint candidate drugs that may exhibit heightened sensitivity in LUAD patients, we conducted drug response predictions using CTRP- and PRISM-derived data. Finally, the cross-correlation of the two pharmacogenomics databases allowed us to predict four drugs or compounds (including SB−743921, ispinesib, cabazitaxel, and gemcitabine) with therapeutic potential in patients (Figures 10H, I).
Figure 10 Prediction of immunotherapy and chemotherapy response. (A) A violin diagram illustrates the variance in TIDE scores among patients with diverse DS. (B) A comprehensive submap analysis of the meta-cohort and melanoma patients, inclusive of intricate immunotherapeutic data. (C) A Kaplan-Meier plot delineates the survival rates for patients categorized into high- and low-DS groups within the IMvigor cohort. (D-G) A box diagram depicts the disparity in DS among patients exhibiting immunotherapy responses in the IMvigor210, GSE35640, GSE79671, and GSE173839 cohorts. (H, I) The findings from the correlation study and differential drug response analysis of CTRP-derived pharmaceuticals and PRISM-derived pharmaceuticals are presented. ***P < 0.001.
Single-cell sequencing analysis
To analyze the expression of DS in TME, we used the LUAD single-cell dataset GSE131907 from the GEO database. All the cells were partitioned into 16 clusters using the k-nearest neighbor (KNN) clustering algorithm (Figure 11A). Subsequently, using the “single R” and “copycat” packages to annotate all cells, we were able to identify 9 distinct cellular subtypes, including B cells, endothelial cells, epithelial cells, cancer cells, macrophage cells, monocyte cells, smooth muscle cells, NK cells, and T cells (Figure 11B). Most of these cells are important components of the TME mentioned in the above results. Subsequently, we investigated the single-cell transcriptome localization of 7 genes in DS (Figure 11C). Concurrently, the DS for each cell was calculated, which showed that cells with a high DS predominately resided in the region of cancer cells (Figure 11D). Additionally, the temporal sequence of cancer cellular differentiation was revealed by the analysis of the pseudotime trajectory. Accordingly, low DS cancer cells appear at an earlier pseudotime than high DS cancer cells, which are primarily found in the earliest stages of differentiation (Figures 11E–H).
Figure 11 Exploration of DS in LUAD scRNA-seq data. (A) t-SNE plot colored by 16 cell subpopulations. (B) t-SNE plot of the distribution of 9 cell types. (C, D) Evaluation of DS gene expression and DS in scRNA-seq data in scRNA-seq data. (E-H) Pseudotime trajectory analysis in LUAD cells (Cells are colored based on states, pseudotime, cluster and DS groups.
KYNU evaluation in the LUAD cells
The qRT-PCR experiments were performed on LUAD cell lines to confirm the expression levels of DS genes in LUAD. Our finding discovered that ABCC2, NEIL3, KYNU, and CPS1 exhibited elevated expression in LUAD cell lines, while KRT6A and SFTPC were found to be underexpressed (Figure 12A). In addition, KYNU exhibited the most significant correlation between high expression and unfavorable poor patient prognosis. No reports have documented the role of KYNU in LUAD to date. As a result of these considerations, KYNU was selected as the focus of further experiments. Correspondingly, IHC staining analysis showed that the protein level of KYNU expression was elevated in LUAD tumor tissues relative to paracancerous tissues (Figure 12B). Subsequently, we investigated the function of KYNU in PC-9 and H838 LUAD cell lines through a series of cell-based experiments. Initially, the effect of the siRNA was confirmed via RT-qPCR (Figure 13A). As demonstrated by CCK-8 and clone formation assays, KYNU knockdown inhibits LUAD cell growth and their clone formation capacity (Figures 13B, C). In addition, wound healing and transwell assays confirmed that KYNU knockdown inhibited the cell migration and invasion capabilities of LUAD cells (Figures 13D, E). The precision of these findings corroborated that the expression of KYNU mirrored the variations anticipated through bioinformatic prediction.
Figure 12 Validation of expression levels of DS genes. (A) DS genes expression in LUAD and normal cell lines. (B) Protein expression levels of KYNU were assessed by IHC. *p < 0.05; **p < 0.01; ***p < 0.001; ns, not significant.
Figure 13 (A) The effect of siRNA to knockdown KYNU in LUAD cell lines was measured by RT–qPCR. (B, C) The CCK-8 and clone formation assays showed that knockdown of KYNU inhibited the proliferation of LUAD cells. (D, E) The wound healing and transwell assays showed that knockdown of KYNU inhibited the migration and invasion of LUAD cells. *P < 0.05; **P < 0.01; ***P < 0.001.
Discussion
LUAD persistently remains the principal contributor to cancer-related deaths among all cancer types and poses a substantial threat to global health (39). Prior research has revealed that the onset and progression of LUAD involve complex biological mechanisms, such as a multitude of genetic and epigenetic modifications (40, 41). Numerous staging systems have been proposed and utilized for clinical determinations to predict patient prognoses; however, these systems predominantly rely on clinicopathological features, ignoring the critical influence of complex molecular pathogenic processes in the oncogenesis and progression of LUAD (42, 43). As a result, there have been negligible improvements in patient outcomes. Consequently, the identification of superior predictive biomarkers for treatment responsiveness and patient outcomes could be advantageous in optimizing individualized therapeutic strategies and prognostic management for those afflicted with LUAD. Atypical accumulation of intracellular disulfides induces disulfide stress, which results in cellular toxicity and ultimately induces cell death (16, 44). Disulfide bonds are the most important redox-reactive covalent bonds between two cysteine residues within proteins. These bonds are regarded as cellular redox regulators and are intimately linked to the formation of disulfides. Recent studies have revealed that neoplastic cells may also experience disturbances in disulfide metabolism as a result of oxidative stress, a process that can potentially inhibit the proliferation of tumor cells and induce their apoptosis (45, 46). In addition, the disulfides inherent to neoplastic cells may act as conduits, modulating the responsiveness to chemotherapeutic agents and immunotherapy, and possibly serving as prognostic markers (47–79). This suggests that the application of disulfidptosis-focused translational medicine holds considerable promise as a candidate for clinical implementation across an array of human malignancies. Several DRs have been implicated in pathological and physiological processes of a variety of tumor. To counterbalance the oxidative stress induced by their heightened metabolic rate, tumor cells can upregulate the expression of the catalytic subunit SLC7A11 of the Xc−system, thereby maintaining high levels of glutathione (50). In addition, the overexpression of SLC7A11 in glioma cells improves their resistance to oxidative stress and decreases their sensitivity to temozolomide (51). In this regard, SLC3A2 is significantly upregulated in several types of malignant tumor cells, including those of the lung, breast, and prostate (52, 53). Furthermore, SLC3A2 is also an independent prognostic indicator for thymic epithelial tumors and NSCLC (54). Consequently, the overexpression of SLC3A2 contributes to radiotherapy resistance in tumors, indicating that SLC3A2 could surface as a promising clinical prospect in cancer treatment (55). Moreover, models composed of DRs have been established in certain tumor types, demonstrating their potential predictive value for patient prognosis and treatment efficacy (56, 57). These findings indicate that DRs have significant potential for elucidating the molecular mechanisms underlying LUAD and identifying novel biomarkers. Nonetheless, there is a dearth of pertinent research on how DRs influence prognosis, immune infiltration, and clinical response in LUAD.
In this study, we first analyzed the characteristics of DRs in LUAD, including extensive genetic and transcriptional level alterations. The majority of these genes are upregulated in LUAD patients and are associated with a poorer prognosis, suggesting a plausible role for DRs in the pathogenesis of LUAD. Using unsupervised clustering techniques on DRs transcriptomic expression data, we then divided LUAD patients from four distinct cohorts into two subgroups, designated DRcluster A and DRcluster B. In DRcluster B, the majority of DRs were significantly upregulated, indicating relatively active disulfidptosis. Compared to patients in DRcluster A, DRcluster B was associated with an increase in the number of immune cells that infiltrated the affected tissue. Significant infiltration levels of effector immune cells are critical for a successful immunotherapeutic response. Typically, a higher CD8+ T cell infiltration rate is indicative of a better prognosis for survival. As evidence, a higher concentration of cytotoxic CD8+ T cells permeating the tumor has been linked to superior outcomes for patients with NSCLC (58). This is consistent with our prognosis and analysis results of immune infiltration. As a result, we proceeded to identify 86 DEGs that distinguished the two DRclusters, and based on these DEGs, we formed a pair of gene clusters. Intriguingly, we discovered statistically significant differences in OS, DRs, and TME between gene clusters, revealing a strong correlation between DRclusters and gene clusters. In light of the lack of DRclusters for clinical application and the paucity of biomarkers for prognosis tracking, we developed a robust and effective model by transforming 10 machine learning algorithms into more than 101 combinations and selecting the best performing algorithm determined by the mean C-index across four LUAD cohorts (59). This facilitated the creation of a robust and efficient DR-based prognostic model, suitable for appraising the prognosis of LUAD patients. Ultimately, the combination of LASSO and RSF was deemed the superior model for constructing the DS. Survival analysis utilizing the median value of the DS revealed its association with LUAD prognosis, and concordant results were obtained from three independent cohorts. The AUC at various time points and the C-index suggest that the DS has exceptional clinical efficacy, surpassing the performance of a substantial number of other clinical attributes. Significantly, when compared to the 56 previously reported molecular signatures for LUAD, the predictive performance of the DS is consistently superior in nearly all cohorts examined.
The composition of our DS comprises 7 genes, including KRT6A, NEIL3, KYNU, ABCC2, SFTPC, CPS1, and INSL4. In addition, qRT-PCR analysis revealed that LUAD and human bronchial epithelial cell expression of the majority of genes differed significantly. Numerous identified genes exhibit a strong correlation with the onset and advancement of LC. For instance, overexpression of KRT6A in NSCLC is associated with poor prognosis (60). KRT6A, acting downstream of LSD1, upregulates G6PD and the pentose phosphate pathway flux via the MYC signaling cascade, thereby promoting NSCLC growth and invasion (61). The upregulation of NEIL3 expression in NSCLC tissues and cell lines correlates with clinical progression and a poor prognosis. By partially activating the PI3K/AKT/mTOR signaling pathway, NEIL3 contributes to the progression of NSCLC (62). A recent study demonstrated that SFTPC expression is suppressed in human LUAD tissues and cell lines, and its overexpression inhibits LCcell proliferation in vitro and in vivo (63). INSL4, via autocrine or paracrine effects, promotes the proliferation and invasion of NSCLC by enhancing the MAPK and AKT signaling pathways. Moreover, INSL4 serves as a detrimental prognostic indicator for patients suffering from NSCLC. Among the 7 genes, IHC confirmed the elevated expression of KYNU in LUAD tissue samples, and the association between high KYNU expressions was notably associated with a shorter OS in the cohort. Moreover, cellular experiments suggest that the knockdown of KYNU can curb the proliferation, invasion, and migration capabilities of LUAD cells, suggesting its oncogenic role in LUAD.
TMB has recently emerged as a promising prognostic biomarker for numerous tumor types. A higher TMB is frequently associated with improved survival outcomes (9). For instance, a study of NSCLC patients revealed that those with elevated TMB levels experienced prolonged OS when subjected to PD-1/PD-L1 antibody therapy (64). In this study, LUAD patients with an escalated DS had elevated TMB. This could be the result of patient heterogeneity or a small sample size. Recently, the proliferation and efficacy of targeted immunotherapies have begun to transform the landscape of cancer treatment (65, 66). Given the complex interaction between the tumor immune microenvironment and host immune responses, there is an urgent need for predictive biomarkers that facilitate individualized therapy. Increased concentrations of CD8+ T cells in the tumor microenvironment correlate with an improved prognosis and increased survival rates among patients with NSCLC (67). Furthermore, the existence of dysfunctional CD8+ T cells within lung tumors and malignant pleural effusions has been documented, thereby diminishing their capacity to mount an effective antitumor response (68). In preclinical studies, NK cell-based therapies have demonstrated the ability to prevent the development of pulmonary metastases (69). Several studies demonstrate that extracorporeal stimulation of autologous Natural Killer (NK) cells with Interleukin-2 (IL-2) in conjunction with adoptive transfer and subcutaneous IL-2 infusions increased overall survival (OS) in a subset of patients with advanced cancers (70). Within the pulmonary environment, neoplasm-associated B lymphocytes can differentiate into plasma cells, thereby producing tumor-specific antibodies capable of recognizing and reacting to tumor-associated antigens (71). Accordingly, the presence of both follicular B cells and tumor-infiltrating plasma cells has been positively correlated with increased longevity in patients with NSCLC, highlighting the protective contribution of plasma cells and antibodies in combating tumor proliferation (72). Moreover, single-cell sequencing analysis results indicate that high DS cells are primarily located in the tumor cell and B cell regions, suggesting a possible interaction between them. Consistent with previous findings, LUAD patients with a low DS demonstrated a high level of immune cell infiltration, including CD8 cells, NK cells, B cells, dendritic cells, mast cells, and central memory T cells, all of which play crucial roles in either bolster or counter tumor immunity during immunotherapy. Tumor cells characterized by reduced differentiation levels frequently demonstrate accelerated growth rates and heightened invasiveness, often correlating with unfavorable prognoses (73). Pseudotime analysis outcomes revealed a spatial disposition wherein tumor cells exhibiting diminished CDS levels occupied the initial phase of the differentiation trajectory. Conversely, those with elevated CDS levels were situated at the concluding stage of differentiation. Consequently, it becomes evident that CDS levels could potentially be linked to the extent of differentiation as well as invasiveness in tumor cells. In accordance with the tumor immunoediting hypothesis, the high-DS group exhibited greater immunosuppression but decreased immunoreactivity compared to the low-DS group (74). Disparities in the immune infiltrating microenvironment could potentially contribute to cancer progression and result in a poorer prognosis. Elevated IDO1 expression contributes to the development of an immunosuppressive TME by promoting T cell and NK cell inactivation and activating and expanding Tregs and DCs (75, 76). The role of IDO1-mediated tryptophan (TRP) metabolism in resistance to therapies targeting CTLA-4 or PD-1 demonstrates its potential as a promising target to augment existing immunotherapy approaches (77). CD40 is an essential co-stimulatory protein involved in the pro-inflammatory immune activation of antigen-presenting cells like dendritic cells and immunosuppressive macrophages within the cancer landscape (78, 79). Earlier investigations demonstrate that CD40 stimulation, in addition to activating tumor-associated immunosuppressive macrophages and T cells and inhibiting tumor progression (80, 81), also remodels the TME and heightens the tumor’s responsiveness to checkpoint blockade therapies in various types of cancer (82, 83). The findings of the current study show that the expression of IDO1 was downregulated in the low-DS group relative to the high-DS group, while CD40 expression was upregulated. In addition, differences in the expression of other immune checkpoints between the two groups suggested immunotherapy would have divergent effects (84). To determine the degree of immunotherapy response in LUAD patients, we evaluated the TIDE and submap algorithms to discovered that low-DS individuals identified by the model may be suitable candidates for immune checkpoint blockade therapies targeting CTLA-4 and PD-1. The scores based on immunotherapy algorithm scores are merely a reflection of theoretical hypotheses and cannot represent the actual efficacy within actual cohorts. For a more comprehensive analysis of the predictive efficacy of DS in immunotherapy, we thus incorporated multiple immunotherapy cohorts. Thus, in predicting therapeutic responses within immunotherapy cohorts, DS exhibits a trend that is consistent with immunotherapy algorithms, as indicated by our findings. These results suggest that a lower DS score may be a potent indicator of immunotherapy response in patients with LUAD.
The efficacy of pharmaceutical interventions is closely correlated with drug sensitivity and individual variation, indicating that personalized therapies based on specific subtypes could reduce the prevalence of ineffective therapies among LUAD patients. An analysis of drug sensitivity differences among LUAD patients with varying DS revealed 4 pharmaceuticals with significantly divergent sensitivities. Subsequent sensitivity projections demonstrated that SB-743921, cabazitaxel, gemcitabine, and ispinesib could potentially serve as superior therapeutic options for individuals with a high DS. Through signaling pathways such as G2M_CHECKPOINT, DNA_REPAIR, and PI3K_AKT_MTOR_SIGNALING, some of these pharmaceuticals exert their antitumor effects, as revealed by our enrichment analyses. Cabazitaxel functions by closely interacting with microtubule proteins, resulting in the inhibition of their depolymerization and consequently the inhibition of cellular mitosis. This mechanism induces cell cycle cessation, resulting in the programmed death or apoptosis of neoplastic cells (85). Recent research revealed that cabazitaxel can induce G2/M phase block and autophagy in LUAD cells by inhibiting the PI3K-AKT-mTOR pathway, indicating its potential as a chemotherapy drug for LUAD patients (86). Gemcitabine is a pyrimidine nucleoside analogue antimetabolite that can inhibit the synthesis and repair of DNA, thereby inducing cellular autophagy and apoptosis (87). When used as a stand-alone treatment, gemcitabine has consistently demonstrated response rates greater than 20% while maintaining a favorable tolerability profile. Moreover, its therapeutic efficacy can be augmented through combination regimens with platinum-based compounds such as cisplatin, thereby synergistically enhancing its overall efficacy (88, 89).
Despite the use of potent open-source data to elucidate two distinct characteristics of LUAD disulfidptosis subtypes and to develop a robust DS evaluation model, this investigation is limited by certain factors. First, this study relies on patient data obtained from publicly available retrospective cohorts and lacks the prospective real-world data required to validate the clinical applicability of the proposed scoring system. In addition, due to financial and resource limitations, we conducted preliminary in vitro experiments to investigate the functionality of KYNU in LUAD. The further experimentation and investigation are required for a more complete comprehension of the molecular mechanisms. This study focuses on bioinformatics analysis and preliminary functional investigations to identify possible biomarkers. These limitations will be addressed in future research.
Conclusions
In conclusion, this study identified disulfidptosis-related subgroups and developed a DS for evaluating the prognosis, immune infiltration, mutations, and treatment sensitivity of LUAD patients. Studying disparities in disulfidptosis patterns has deepened our comprehension of both tumor heterogeneity and the intricate complexities within the TME. We have also constructed and validated a DS that accurately predicts patient prognosis and treatment efficacy assessment, offering a potentially powerful new tool for clinical decision-making, patient outcomes, and individualized treatment strategies. Besides, this study has contributed to advancing the understanding of molecular complexity in LUAD and provides directions and potential avenues for future LUAD research.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by Ethics Committee and Institutional Review Board of the Outdo Biotech. Co., Ltd. (SHYJS-CP-2210035). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
CY and YZ conceived and designed the study. CL and YW performed the collection and assembly of data. YX and CL analyzed the data. YZ and YW performed experiments and wrote the manuscript. All authors read and approved the final manuscript.
Acknowledgments
We thank Bullet Edits Limited for the linguistic editing and proofreading of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1233260/full#supplementary-material
References
1. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin (2015) 65(2):87–108. doi: 10.3322/caac.21262
2. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2022. CA Cancer J Clin (2022) 72(1):7–33. doi: 10.3322/caac.21708
3. Song C, Wu Z, Wang Q, Wang Y, Guo Z, Li S, et al. A combined two-mRNA signature associated with PD-L1 and tumor mutational burden for prognosis of lung adenocarcinoma. Front Cell Dev Biol (2021) 9:634697. doi: 10.3389/fcell.2021.634697
4. Reck M, Heigener DF, Mok T, Soria JC, Rabe KF. Management of non-small-cell lung cancer: recent developments. Lancet (2013) 382(9893):709–19. doi: 10.1016/s0140-6736(13)61502-0
5. Miller KD, Nogueira L, Devasia T, Mariotto AB, Yabroff KR, Jemal A, et al. Cancer treatment and survivorship statistics, 2022. CA Cancer J Clin (2022) 72(5):409–36. doi: 10.3322/caac.21731
6. Tower J. Programmed cell death in aging. Ageing Res Rev (2015) 23(Pt A):90–100. doi: 10.1016/j.arr.2015.04.002
7. Liang C, Zhang X, Yang M, Dong X. Recent progress in ferroptosis inducers for cancer therapy. Adv Mater (2019) 31(51):e1904197. doi: 10.1002/adma.201904197
8. Xie Y, Zhu S, Song X, Sun X, Fan Y, Liu J, et al. The tumor suppressor p53 limits ferroptosis by blocking DPP4 activity. Cell Rep (2017) 20(7):1692–704. doi: 10.1016/j.celrep.2017.07.055
9. Wang W, Green M, Choi JE, Gijón M, Kennedy PD, Johnson JK, et al. CD8(+) T cells regulate tumour ferroptosis during cancer immunotherapy. Nature (2019) 569(7755):270–4. doi: 10.1038/s41586-019-1170-y
10. Rangel M, Kong J, Bhatt V, Khayati K, Guo JY. Autophagy and tumorigenesis. FEBS J (2022) 289(22):7177–98. doi: 10.1111/febs.16125
11. Mulcahy Levy JM, Thorburn A. Autophagy in cancer: moving from understanding mechanism to improving therapy responses in patients. Cell Death Differ (2020) 27(3):843–57. doi: 10.1038/s41418-019-0474-7
12. Liu X, Nie L, Zhang Y, Yan Y, Wang C, Colic M, et al. Actin cytoskeleton vulnerability to disulfide stress mediates disulfidptosis. Nat Cell Biol (2023) 25(3):404–414. doi: 10.1038/s41556-023-01091-2
13. Shin CS, Mishra P, Watrous JD, Carelli V, D’Aurelio M, Jain M, et al. The glutamate/cystine xCT antiporter antagonizes glutamine metabolism and reduces nutrient flexibility. Nat Commun (2017) 8:15074. doi: 10.1038/ncomms15074
14. Koppula P, Zhang Y, Shi J, Li W, Gan B. The glutamate/cystine antiporter SLC7A11/xCT enhances cancer cell dependency on glucose by exporting glutamate. J Biol Chem (2017) 292(34):14240–9. doi: 10.1074/jbc.M117.798405
15. Goji T, Takahara K, Negishi M, Katoh H. Cystine uptake through the cystine/glutamate antiporter xCT triggers glioblastoma cell death under glucose deprivation. J Biol Chem (2017) 292(48):19721–32. doi: 10.1074/jbc.M117.814392
16. Liu X, Olszewski K, Zhang Y, Lim EW, Shi J, Zhang X, et al. Cystine transporter regulation of pentose phosphate pathway dependency and disulfide stress exposes a targetable metabolic vulnerability in cancer. Nat Cell Biol (2020) 22(4):476–86. doi: 10.1038/s41556-020-0496-x
17. Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics (2012) 28(6):882–3. doi: 10.1093/bioinformatics/bts034
18. Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics (2010) 26(12):1572–3. doi: 10.1093/bioinformatics/btq170
19. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43(7):e47. doi: 10.1093/nar/gkv007
20. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118
21. Binder H, Allignol A, Schumacher M, Beyersmann J. Boosting for high-dimensional time-to-event data with competing risks. Bioinformatics (2009) 25(7):890–6. doi: 10.1093/bioinformatics/btp088
22. Bair E, Tibshirani R. Semi-supervised methods to predict patient survival from gene expression data. PloS Biol (2004) 2(4):E108. doi: 10.1371/journal.pbio.0020108
23. Mayakonda A, Lin DC, Assenov Y, Plass C, Koeffler HP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res (2018) 28(11):1747–56. doi: 10.1101/gr.239244.118
24. Mermel CH, Schumacher SE, Hill B, Meyerson ML, Beroukhim R, Getz G. GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol (2011) 12(4):R41. doi: 10.1186/gb-2011-12-4-r41
25. Chakraborty H, Hossain A. R package to estimate intracluster correlation coefficient with confidence interval for binary data. Comput Methods Programs BioMed (2018) 155:85–92. doi: 10.1016/j.cmpb.2017.10.023
26. Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep (2017) 18(1):248–62. doi: 10.1016/j.celrep.2016.12.019
27. 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. doi: 10.1186/s13059-016-1070-5
28. Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity (2013) 39(4):782–95. doi: 10.1016/j.immuni.2013.10.003
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. doi: 10.1158/0008-5472.Can-18-0689
30. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U.S.A. (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102
31. Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinf (2013) 14:7. doi: 10.1186/1471-2105-14-7
32. Jiang P, Gu S, Pan D, Fu J, Sahu A, Hu X, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med (2018) 24(10):1550–8. doi: 10.1038/s41591-018-0136-1
33. Roh W, Chen PL, 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. doi: 10.1126/scitranslmed.aah3560
34. Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM3, et al. Comprehensive integration of single-cell data. Cell (2019) 177(7):1888–1902.e21. doi: 10.1016/j.cell.2019.05.031
35. Aran D, Looney AP, Liu L, Wu E, Fong V, Hsu A, et al. Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage. Nat Immunol (2019) 20(2):163–72. doi: 10.1038/s41590-018-0276-y
36. Gao R, Bai S, Henderson YC, Lin Y, Schalck A, Yan Y, et al. Delineating copy number and clonal substructure in human tumors from single-cell transcriptomes. Nat Biotechnol (2021) 39(5):599–608. doi: 10.1038/s41587-020-00795-2
37. Qiu X, Mao Q, Tang Y, Wang L, Chawla R, Pliner HA, et al. Reversed graph embedding resolves complex single-cell trajectories. Nat Methods (2017) 14(10):979–82. doi: 10.1038/nmeth.4402
38. Allred DC, Harvey JM, Berardo M, Clark GM. Prognostic and predictive factors in breast cancer by immunohistochemical analysis. Mod Pathol (1998) 11(2):155–68.
39. Yang D, Liu Y, Bai C, Wang X, Powell CA. Epidemiology of lung cancer and lung cancer screening programs in China and the United States. Cancer Lett (2020) 468:82–7. doi: 10.1016/j.canlet.2019.10.009
40. Kabbout M, Garcia MM, Fujimoto J, Liu DD, Woods D, Chow CW, et al. ETS2 mediated tumor suppressive function and MET oncogene inhibition in human non-small cell lung cancer. Clin Cancer Res (2013) 19(13):3383–95. doi: 10.1158/1078-0432.Ccr-13-0341
41. Saji H, Tsuboi M, Shimada Y, Kato Y, Hamanaka W, Kudo Y, et al. Gene expression profiling and molecular pathway analysis for the identification of early-stage lung adenocarcinoma patients at risk for early recurrence. Oncol Rep (2013) 29(5):1902–6. doi: 10.3892/or.2013.2332
42. Westaway DD, Toon CW, Farzin M, Sioson L, Watson N, Brady PW, et al. The International Association for the Study of Lung Cancer/American Thoracic Society/European Respiratory Society grading system has limited prognostic significance in advanced resected pulmonary adenocarcinoma. Pathology (2013) 45(6):553–8. doi: 10.1097/PAT.0b013e32836532ae
43. Rosell R, Carcereny E, Gervais R, Vergnenegre A, Massuti B, Felip E, et al. Erlotinib versus standard chemotherapy as first-line treatment for European patients with advanced EGFR mutation-positive non-small-cell lung cancer (EURTAC): a multicentre, open-label, randomised phase 3 trial. Lancet Oncol (2012) 13(3):239–46. doi: 10.1016/s1470-2045(11)70393-x
44. Joly JH, Delfarah A, Phung PS, Parrish S, Graham NA. A synthetic lethal drug combination mimics glucose deprivation-induced cancer cell death in the presence of glucose. J Biol Chem (2020) 295(5):1350–65. doi: 10.1074/jbc.RA119.011471
45. Hogg PJ. Biological regulation through protein disulfide bond cleavage. Redox Rep (2002) 7(2):71–7. doi: 10.1179/135100002125000299
46. Daly EB, Wind T, Jiang XM, Sun L, Hogg PJ. Secretion of phosphoglycerate kinase from tumour cells is controlled by oxygen-sensing hydroxylases. Biochim Biophys Acta (2004) 1691(1):17–22. doi: 10.1016/j.bbamcr.2003.11.004
47. Zhang S, Guo N, Wan G, Zhang T, Li C, Wang Y, et al. pH and redox dual-responsive nanoparticles based on disulfide-containing poly(β-amino ester) for combining chemotherapy and COX-2 inhibitor to overcome drug resistance in breast cancer. J Nanobiotechnol (2019) 17(1):109. doi: 10.1186/s12951-019-0540-9
48. Zhang W, Quan Y, Ma X, Zeng L, Li J, Chen S, et al. Synergistic effect of glutathione and IgG4 in immune evasion and the implication for cancer immunotherapy. Redox Biol (2023) 60:102608. doi: 10.1016/j.redox.2023.102608
49. Karatas F, Acat M, Sahin S, Inci F, Karatas G, Neselioglu S, et al. The prognostic and predictive significance of serum thiols and disulfide levels in advanced non-small cell lung cancer. Aging Male (2020) 23(5):619–28. doi: 10.1080/13685538.2018.1559805
50. Koppula P, Zhuang L, Gan B. Cystine transporter SLC7A11/xCT in cancer: ferroptosis, nutrient dependency, and cancer therapy. Protein Cell (2021) 12(8):599–620. doi: 10.1007/s13238-020-00789-5
51. Polewski MD, Reveron-Thornton RF, Cherryholmes GA, Marinov GK, Aboody KS. SLC7A11 overexpression in glioblastoma is associated with increased cancer stem cell-like properties. Stem Cells Dev (2017) 26(17):1236–46. doi: 10.1089/scd.2017.0123
52. Digomann D, Kurth I, Tyutyunnykova A, Chen O, Löck S, Gorodetska I, et al. The CD98 heavy chain is a marker and regulator of head and neck squamous cell carcinoma radiosensitivity. Clin Cancer Res (2019) 25(10):3152–63. doi: 10.1158/1078-0432.Ccr-18-2951
53. Hasegawa K, Ikeda S, Yaga M, Watanabe K, Urakawa R, Iehara A, et al. Selective targeting of multiple myeloma cells with a monoclonal antibody recognizing the ubiquitous protein CD98 heavy chain. Sci Transl Med (2022) 14(632):eaax7706. doi: 10.1126/scitranslmed.aax7706
54. Kanai Y. Amino acid transporter LAT1 (SLC7A5) as a molecular target for cancer diagnosis and therapeutics. Pharmacol Ther (2022) 230:107964. doi: 10.1016/j.pharmthera.2021.107964
55. Bajaj J, Konuma T, Lytle NK, Kwon HY, Ablack JN, Cantor JM, et al. CD98-mediated adhesive signaling enables the establishment and propagation of acute myelogenous leukemia. Cancer Cell (2016) 30(5):792–805. doi: 10.1016/j.ccell.2016.10.003
56. Zhao S, Wang L, Ding W, Ye B, Cheng C, Shao J, et al. Crosstalk of disulfidptosis-related subtypes, establishment of a prognostic signature and immune infiltration characteristics in bladder cancer based on a machine learning survival framework. Front Endocrinol (Lausanne) (2023) 14:1180404. doi: 10.3389/fendo.2023.1180404
57. Wang T, Guo K, Zhang D, Wang H, Yin J, Cui H, et al. Disulfidptosis classification of hepatocellular carcinoma reveals correlation with clinical prognosis and immune profile. Int Immunopharmacol (2023) 120:110368. doi: 10.1016/j.intimp.2023.110368
58. Barnes TA, Amir E. HYPE or HOPE: the prognostic value of infiltrating immune cells in cancer. Br J Cancer (2017) 117(4):451–60. doi: 10.1038/bjc.2017.220
59. Liu Z, Liu L, Weng S, Guo C, Dang Q, Xu H, et al. Machine learning-based integration develops an immune-derived lncRNA signature for improving outcomes in colorectal cancer. Nat Commun (2022) 13(1):816. doi: 10.1038/s41467-022-28421-6
60. Xiao J, Lu X, Chen X, Zou Y, Liu A, Li W, et al. Eight potential biomarkers for distinguishing between lung adenocarcinoma and squamous cell carcinoma. Oncotarget (2017) 8(42):71759–71. doi: 10.18632/oncotarget.17606
61. Che D, Wang M, Sun J, Li B, Xu T, Lu Y, et al. KRT6A promotes lung cancer cell growth and invasion through MYC-regulated pentose phosphate pathway. Front Cell Dev Biol (2021) 9:694071. doi: 10.3389/fcell.2021.694071
62. Huang H, Hua Q. NEIL3 mediates lung cancer progression and modulates PI3K/AKT/mTOR signaling: A potential therapeutic target. Int J Genomics (2022) 2022:8348499. doi: 10.1155/2022/8348499
63. Li B, Meng YQ, Li Z, Yin C, Lin JP, Zhu DJ, et al. MiR-629-3p-induced downregulation of SFTPC promotes cell proliferation and predicts poor survival in lung adenocarcinoma. Artif Cells Nanomed Biotechnol (2019) 47(1):3286–96. doi: 10.1080/21691401.2019.1648283
64. Rizvi NA, Hellmann MD, Snyder A, Kvistborg P, Makarov V, Havel JJ, et al. Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science (2015) 348(6230):124–8. doi: 10.1126/science.aaa1348
65. Zhou F, Qiao M, Zhou C. The cutting-edge progress of immune-checkpoint blockade in lung cancer. Cell Mol Immunol (2021) 18(2):279–93. doi: 10.1038/s41423-020-00577-5
66. Bejarano L, Jordāo MJC, Joyce JA. Therapeutic targeting of the tumor microenvironment. Cancer Discovery (2021) 11(4):933–59. doi: 10.1158/2159-8290.Cd-20-1808
67. Fumet JD, Richard C, Ledys F, Klopfenstein Q, Joubert P, Routy B, et al. Prognostic and predictive role of CD8 and PD-L1 determination in lung tumor tissue of patients under anti-PD-1 therapy. Br J Cancer (2018) 119(8):950–60. doi: 10.1038/s41416-018-0220-9
68. Prado-Garcia H, Romero-Garcia S, Aguilar-Cazares D, Meneses-Flores M, Lopez-Gonzalez JS. Tumor-induced CD8+ T-cell dysfunction in lung cancer patients. Clin Dev Immunol (2012) 2012:741741. doi: 10.1155/2012/741741
69. Kim S, Iizuka K, Aguila HL, Weissman IL, Yokoyama WM. In vivo natural killer cell activities revealed by natural killer cell-deficient mice. Proc Natl Acad Sci U.S.A. (2000) 97(6):2731–6. doi: 10.1073/pnas.050588297
70. Rosenberg SA, Lotze MT, Yang JC, Topalian SL, Chang AE, Schwartzentruber DJ, et al. Prospective randomized trial of high-dose interleukin-2 alone or in conjunction with lymphokine-activated killer cells for the treatment of patients with advanced cancer. J Natl Cancer Inst (1993) 85(8):622–32. doi: 10.1093/jnci/85.8.622
71. Germain C, Gnjatic S, Tamzalit F, Knockaert S, Remark R, Goc J, et al. Presence of B cells in tertiary lymphoid structures is associated with a protective immunity in patients with lung cancer. Am J Respir Crit Care Med (2014) 189(7):832–44. doi: 10.1164/rccm.201309-1611OC
72. Lohr M, Edlund K, Botling J, Hammad S, Hellwig B, Othman A, et al. The prognostic relevance of tumour-infiltrating plasma cells and immunoglobulin kappa C indicates an important role of the humoral immune response in non-small cell lung cancer. Cancer Lett (2013) 333(2):222–8. doi: 10.1016/j.canlet.2013.01.036
73. Ana Kátia S, Teixeira JL, Vasconcelos A. Histopathological profile of patients diagnosed with Malignant tumors assisted in a hospital of reference of Agreste Pernambucano. Jornal Brasileiro Patologia e Medicina Laboratorial (2019) 55(1). doi: 10.5935/1676-2444.20190002
74. Hellmann MD, Nathanson T, Rizvi H, Creelan BC, Sanchez-Vega F, Ahuja A, et al. Genomic features of response to combination immunotherapy in patients with advanced non-small-cell lung cancer. Cancer Cell (2018) 33(5):843–852.e4. doi: 10.1016/j.ccell.2018.03.018
75. Frumento G, Rotondo R, Tonetti M, Damonte G, Benatti U, Ferrara GB. Tryptophan-derived catabolites are responsible for inhibition of T and natural killer cell proliferation induced by indoleamine 2,3-dioxygenase. J Exp Med (2002) 196(4):459–68. doi: 10.1084/jem.20020121
76. Mellor AL, Keskin DB, Johnson T, Chandler P, Munn DH. Cells expressing indoleamine 2,3-dioxygenase inhibit T cell responses. J Immunol (2002) 168(8):3771–6. doi: 10.4049/jimmunol.168.8.3771
77. Holmgaard RB, Zamarin D, Munn DH, Wolchok JD, Allison JP. Indoleamine 2,3-dioxygenase is a critical resistance mechanism in antitumor T cell immunotherapy targeting CTLA-4. J Exp Med (2013) 210(7):1389–402. doi: 10.1084/jem.20130066
78. Alexandroff AB, Jackson AM, Paterson T, Haley JL, Ross JA, Longo DL, et al. Role for CD40-CD40 ligand interactions in the immune response to solid tumours. Mol Immunol (2000) 37(9):515–26. doi: 10.1016/s0161-5890(00)00079-1
79. Vonderheide RH. CD40 agonist antibodies in cancer immunotherapy. Annu Rev Med (2020) 71:47–58. doi: 10.1146/annurev-med-062518-045435
80. von Leoprechting A, van der Bruggen P, Pahl HL, Aruffo A, Simon JC. Stimulation of CD40 on immunogenic human Malignant melanomas augments their cytotoxic T lymphocyte-mediated lysis and induces apoptosis. Cancer Res (1999) 59(6):1287–94.
81. French RR, Chan HT, Tutt AL, Glennie MJ. CD40 antibody evokes a cytotoxic T-cell response that eradicates lymphoma and bypasses T-cell help. Nat Med (1999) 5(5):548–53. doi: 10.1038/8426
82. Ma HS, Poudel B, Torres ER, Sidhom JW, Robinson TM, Christmas B, et al. A CD40 agonist and PD-1 antagonist antibody reprogram the microenvironment of nonimmunogenic tumors to allow T-cell-mediated anticancer activity. Cancer Immunol Res (2019) 7(3):428–42. doi: 10.1158/2326-6066.Cir-18-0061
83. Zhang J, Li Y, Yang S, Zhang L, Wang W. Anti-CD40 mAb enhanced efficacy of anti-PD1 against osteosarcoma. J Bone Oncol (2019) 17:100245. doi: 10.1016/j.jbo.2019.100245
84. Hiam-Galvez KJ, Allen BM, Spitzer MH. Systemic immunity in cancer. Nat Rev Cancer (2021) 21(6):345–59. doi: 10.1038/s41568-021-00347-z
85. Mellado B, Jimenez N, Marin-Aguilera M, Reig O. Diving into cabazitaxel’s mode of action: more than a taxane for the treatment of castration-resistant prostate cancer patients. Clin Genitourin Cancer (2016) 14(4):265–70. doi: 10.1016/j.clgc.2015.12.030
86. Huo R, Wang L, Liu P, Zhao Y, Zhang C, Bai B, et al. Cabazitaxel-induced autophagy via the PI3K/Akt/mTOR pathway contributes to A549 cell death. Mol Med Rep (2016) 14(4):3013–20. doi: 10.3892/mmr.2016.5648
87. Noble S, Goa KL. Gemcitabine. A review of its pharmacology and clinical potential in non-small cell lung cancer and pancreatic cancer. Drugs (1997) 54(3):447–72. doi: 10.2165/00003495-199754030-00009
88. Crinò L, Cappuzzo F. Gemcitabine in non-small cell lung cancer. Expert Opin Pharmacother (2002) 3(6):745–53. doi: 10.1517/14656566.3.6.745
Keywords: disulfidptosis, lung adenocarcinoma, prognosis, tumor microenvironment, immunotherapy
Citation: Wang Y, Xu Y, Liu C, Yuan C and Zhang Y (2023) Identification of disulfidptosis-related subgroups and prognostic signatures in lung adenocarcinoma using machine learning and experimental validation. Front. Immunol. 14:1233260. doi: 10.3389/fimmu.2023.1233260
Received: 01 June 2023; Accepted: 04 September 2023;
Published: 20 September 2023.
Edited by:
Ines Zidi, Tunis El Manar University, TunisiaReviewed by:
Fereshteh Yazdanpanah, University of Pennsylvania, United StatesQiangzhe Zhang, Nankai University, China
Marwa Hasni, Higher National Engineering School of Tunis, Tunisia
Copyright © 2023 Wang, Xu, Liu, Yuan and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yi Zhang, NDcyMTg5OTI2QHFxLmNvbQ==; Chengliang Yuan, MTM1NjgyMzU2MjhAMTYzLmNvbQ==
†These authors have contributed equally to this work and share first authorship