Dataset Information

Identification and validation of an epithelial mesenchymal transition-related gene pairs signature for prediction of overall survival in patients with skin cutaneous melanoma

ABSTRACT:

Background

We aimed to construct a novel epithelial-mesenchymal transition (EMT)-related gene pairs (ERGPs) signature to predict overall survival (OS) in skin cutaneous melanoma (CM) patients.

Methods

Expression data of the relevant genes, corresponding clinicopathological parameters, and follow-up data were obtained from The Cancer Genome Atlas database. Univariate Cox regression analysis was utilized to identify ERGPs significantly associated with OS, and LASSO analysis was used to identify the genes used for the construction of the ERGPs signature. The optimal cutoff value determined by the receiver operating characteristic curve was used to classify patients into high-risk and low-risk groups. Survival curves were generated using the Kaplan–Meier method, and differences between the two groups were estimated using the log-rank test. The independent external datasets GSE65904 and GSE19234 were used to verify the performance of the ERGPs signature using the area under the curve (AUC) values. In addition, we also integrated clinicopathological parameters and risk scores to develop a nomogram that can individually predict the prognosis of patients with CM.

Results

A total of 104 ERGPs related to OS were obtained, of which 21 ERGPs were selected for the construction of the signature. All CM patients were stratified into high-and low-risk groups based on an optimal risk score cutoff value of 0.281. According to the Kaplan–Meier analysis, the mortality rate in the low-risk group was lower than that in the high-risk group in the TCGA cohort (P < 0.001), GSE65904 cohort (P = 0.006), and GSE19234 cohort (P = 0.002). Multivariate Cox regression analysis indicated that our ERGP signature was an independent risk factor for OS in CM patients in the three cohorts (for TCGA: HR, 2.560; 95% CI [1.907–3.436]; P < 0.001; for GSE65904: HR = 2.235, 95% CI [1.492–3.347], P < 0.001; for GSE19234: HR = 2.458, 95% CI [1.065–5.669], P = 0.035). The AUC value for predicting the 5-year survival rate of patients with CM of our developed model was higher than that of two previously established prognostic signatures. Both the calibration curve and the C-index (0.752, 95% CI [0.678–0.826]) indicated that the developed nomogram was highly accurate. Most importantly, the decision curve analysis results showed that the nomogram had a higher net benefit than that of the American Joint Committee on Cancer stage system.

Conclusion

Our study established an ERGPs signature that could be potentially used in a clinical setting as a genetic biomarker for risk stratification of CM patients. In addition, the ERGPs signature could also predict which CM patients will benefit from PD-1 and PD-L1 inhibitors.

SUBMITTER: Shi Y

PROVIDER: S-EPMC8785661 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:The prognosis of skin cutaneous melanoma (SKCM) remains poor, and patients with SKCM show a poor response to immunotherapy. Thus, we aimed to identify necroptosis-related biomarkers, which can help predict the prognosis of SKCM and improve the effectiveness of precision medicine. Data of SKCM were obtained from The Cancer Genome Atlas (TCGA) and GEO databases. TCGA samples were classified into two clusters by consensus clustering of necroptosis-related genes. Univariate Cox and least absolute shrinkage and selection operator regression analyses led to the identification of 11 genes, which were used to construct a prognostic model. GSE65904 was used as the test set. Principal component, t-distributed stochastic neighbor embedding, and Kaplan–Meier survival analyses indicated that samples in the train and test sets could be divided into two groups, with the high-risk group showing a worse prognosis. Univariate and multivariate Cox regression analyses were performed, and a nomogram, calibration curve, and time-dependent receiver operating characteristic curve were constructed to verify the efficacy of our model. The 1-, 3-, and 5-year areas under the receiver operating characteristic curves for the train set were 0.702, 0.663, and 0.701 and for the test set were 0.613, 0.627, and 0.637, respectively. Moreover, we performed Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses between the high- and low-risk groups. Single sample gene set enrichment analysis, immune cell infiltration analysis, tumor microenvironment scores, immune checkpoint analysis, and half-maximal inhibitory concentration prediction indicated that the high-risk group showed weaker antitumor immunity; further, the response to immune checkpoint inhibitors was worse, and the high-risk group was sensitive to fewer antitumor drugs. Tumor mutational burden analysis, Kaplan–Meier survival analysis, and correlation analysis between risk score and RNA stemness score revealed that the high-risk group with low tumor mutational burden and high RNA stemness score was potentially associated with poor prognosis. To conclude, our model, which was based on 11 necroptosis-related genes, could predict the prognosis of SKCM; in addition, it has guiding significance for the selection of clinical treatment and provides new research directions to enhance necroptosis against SKCM.