BackgroundNecroptosis is crucial for organismal development and pathogenesis. To date, the role of necroptosis in skin cutaneous melanoma (SKCM) is yet unveiled. In addition, the part of melanin pigmentation was largely neglected in the bioinformatic analysis. In this study, we aimed to construct a novel prognostic model based on necroptosis-related genes and analysis the pigmentation phenotype of patients to provide clinically actionable information for SKCM patients.MethodsWe downloaded the SKCM data from the TCGA and GEO databases in this study and identified the differently expressed and prognostic necroptosis-related genes. Patients’ pigmentation phenotype was evaluated by the GSVA method. Then, using Lasso and Cox regression analysis, a novel prognostic model was constructed based on the intersected genes. The risk score was calculated and the patients were divided into two groups. The survival differences between the two groups were compared using Kaplan-Meier analysis. The ROC analysis was performed and the area under curves was calculated to evaluate the prediction performances of the model. Then, the GO, KEGG and GSEA analyses were performed to elucidate the underlying mechanisms. Differences in the tumor microenvironment, patients’ response to immune checkpoint inhibitors (ICIs) and pigmentation phenotype were analyzed. In order to validate the mRNA expression levels of the selected genes, quantitative real-time PCR (qRT-PCR) was performed.ResultsAltogether, a novel prognostic model based on four genes (BOK, CD14, CYLD and FASLG) was constructed, and patients were classified into high and low-risk groups based on the median risk score. Low-risk group patients showed better survival status. The model showed high accuracy in the training and the validation cohort. Pathway and functional enrichment analysis indicated that immune-related pathways were differently activated in the two groups. In addition, immune cells infiltration patterns and sensitivity of ICIs showed a significant difference between patients from two risk groups. The pigmentation score was positively related to the risk score in pigmentation phenotype analysis.ConclusionIn conclusion, this study established a novel prognostic model based on necroptosis-related genes and revealed the possible connections between necroptosis and melanin pigmentation. It is expected to provide a reference for clinical treatment.

PurposeFerroptosis-related lncRNAs are promising biomarkers for predicting the prognosis of many cancers. However, a ferroptosis-related signature to predict the prognosis of cutaneous melanoma (CM) has not been identified. The purpose of this study was to construct a ferroptosis-related lncRNA signature to predict prognosis and immunotherapy efficacy in CM.MethodsFerroptosis-related differentially expressed genes (FDEGs) and lncRNAs (FDELs) were identified using TCGA, GTEx, and FerrDb datasets. We performed Cox and LASSO regressions to identify key FDELs, and constructed a risk score to stratify patients into high- and low-risk groups. The lncRNA signature was evaluated using the areas under the receiver operating characteristic curves (AUCs) and Kaplan-Meier analyses in the training, testing, and entire cohorts. Multivariate Cox regression analyses including the lncRNA signature and common clinicopathological characteristics were performed to identify independent predictors of overall survival (OS). A nomogram was developed for clinical use. We performed gene set enrichment analyses (GSEA) to identify significantly enriched pathways. Differences in the tumor microenvironment (TME) between the 2 groups were assessed using 7 algorithms. To predict the efficacy of immune checkpoint inhibitors (ICI), we analyzed the association between PD1 and CTLA4 expression and the risk score. Finally, differences in Tumor Mutational Burden (TMB) and molecular drugs Sensitivity between the 2 groups were performed.ResultsWe identified 5 lncRNAs (AATBC, AC145423.2, LINC01871, AC125807.2, and AC245041.1) to construct the risk score. The AUC of the lncRNA signature was 0.743 in the training cohort and was validated in the testing and entire cohorts. Kaplan-Meier analyses revealed that the high-risk group had poorer prognosis. Multivariate Cox regression showed that the lncRNA signature was an independent predictor of OS with higher accuracy than traditional clinicopathological features. The 1-, 3-, and 5-year survival probabilities for CM patients were 92.7%, 57.2%, and 40.2% with an AUC of 0.804, indicating a good accuracy and reliability of the nomogram. GSEA showed that the high-risk group had lower ferroptosis and immune response. TME analyses confirmed that the high-risk group had lower immune cell infiltration (e.g., CD8+ T cells, CD4+ memory-activated T cells, and M1 macrophages) and lower immune functions (e.g., immune checkpoint activation). Low-risk patients whose disease expressed PD1 or CTLA4 were likely to respond better to ICIs. The analysis demonstrated that the TMB had significantly difference between low- and high- risk groups. Chemotherapy drugs, such as sorafenib, Imatinib, ABT.888 (Veliparib), Docetaxel, and Paclitaxel showed Significant differences in the estimated IC50 between the two risk groups.ConclusionOur novel ferroptosis-related lncRNA signature was able to accurately predict the prognosis and ICI outcomes of CM patients. These ferroptosis-related lncRNAs might be potential biomarkers and therapeutic targets for CM.

PurposeThe purpose of this study was to construct a gene signature comprising genes related to both inflammation and pyroptosis (GRIPs) to predict the prognosis of patients with cutaneous melanoma patients and the efficacy of immunotherapy, chemotherapy, and targeted therapy in these patients.MethodsGene expression profiles were collected from The Cancer Genome Atlas. Weighted gene co-expression network analysis was performed to identify GRIPs. Univariable Cox regression and Lasso regression further selected key prognostic genes. Multivariable Cox regression was used to construct a risk score, which stratified patients into high- and low-risk groups. Areas under the ROC curves (AUCs) were calculated, and Kaplan-Meier analyses were performed for the two groups, following validation in an external cohort from Gene Expression Omnibus (GEO). A nomogram including the GRIP signature and clinicopathological characteristics was developed for clinical use. Gene set enrichment analysis illustrated differentially enriched pathways. Differences in the tumor microenvironment (TME) between the two groups were assessed. The efficacies of immune checkpoint inhibitors (ICIs), chemotherapeutic agents, and targeted agents were predicted for both groups. Immunohistochemical analyses of the GRIPs between the normal and CM tissues were performed using the Human Protein Atlas data. The qRT-PCR experiments validated the expression of genes in CM cell lines, Hacat, and PIG1 cell lines.ResultsA total of 185 GRIPs were identified. A novel gene signature comprising eight GRIPs (TLR1, CCL8, EMP3, IFNGR2, CCL25, IL15, RTP4, and NLRP6) was constructed. The signature had AUCs of 0.714 and 0.659 for predicting 3-year overall survival (OS) in the TCGA entire and GEO validation cohorts, respectively. Kaplan-Meier analyses revealed that the high-risk group had a poorer prognosis. Multivariable Cox regression showed that the GRIP signature was an independent predictor of OS with higher accuracy than traditional clinicopathological features. The nomogram showed good accuracy and reliability in predicting 3-year OS (AUC = 0.810). GSEA and TME analyses showed that the high-risk group had lower levels of pyroptosis, inflammation, and immune response, such as lower levels of CD8+ T-cell infiltration, CD4+ memory-activated T-cell infiltration, and ICI. In addition, low-risk patients whose disease expressed PD-1 or CTLA-4 were likely to respond better to ICIs, and several chemotherapeutic and targeted agents. Immunohistochemical analysis confirmed the distinct expression of five out of the eight GRIPs between normal and CM tissues.ConclusionOur novel 8-GRIP signature can accurately predict the prognosis of patients with CM and the efficacies of multiple anticancer therapies. These GRIPs might be potential prognostic biomarkers and therapeutic targets for CM.

<abstract> <p>Skin cutaneous melanoma (SKCM) is one of the most malignant skin cancers and remains a health concern worldwide. Pyroptosis is a newly recognized form of programmed cell death and plays a vital role in cancer progression. We aim to construct a prognostic model for SKCM patients based on pyroptosis-related genes (PRGs). SKCM patients from The Cancer Genome Atlas (TCGA) were divided into training and validation cohorts. We used GSE65904 downloaded from GEO database as an external validation cohort. We performed Cox regression and the least absolute shrinkage and selection operator (LASSO) regression to identify prognostic genes and built a risk score. Patients were divided into high- and low-risk groups based on the risk score. Differently expressed genes (DEGs), immune cell infiltration and immune-related pathways activation were compared between the two groups. We established a model containing 4 PRGs, i.e., GSDMA, GSDMC, AIM2 and NOD2. The overall survival (OS) time was significantly different between the 2 groups. The risk score was an independent predictor for prognosis in both the uni- and multi-variable Cox regressions. Gene ontology (GO) and Kyoto Encylopedia of Genes and Genomes (KEGG) analyses showed that DEGs were enriched in immune-related pathways. Most types of immune cells were highly expressed in the low risk group. All immune pathways were significantly up-regulated in the low-risk group. In addition, low-risk patients had a better response to immune checkpoint inhibitors. Our novel pyroptosis-related gene signature could predict the prognosis of SKCM patients and their response to immune checkpoint inhibitors.</p> </abstract>

PurposeFatty acid metabolism (FAM) affects the immune phenotype in a metabolically dynamic tumor microenvironment (TME), but the use of FAM-related genes (FAMGs) to predict the prognosis and immunotherapy response of cutaneous melanoma (CM) patients has not been investigated. In this study, we aimed to construct FAM molecular subtypes and identify key prognostic biomarkers in CM.MethodsWe used a CM dataset in The Cancer Genome Atlas (TCGA) to construct FAM molecular subtypes. We performed Kaplan–Meier (K-M) analysis, gene set enrichment analysis (GSEA), and TME analysis to assess differences in the prognosis and immune phenotype between subtypes. We used weighted gene co-expression network analysis (WGCNA) to identify key biomarkers that regulate tumor metabolism and immunity between the subtypes. We compared overall survival (OS), progression-free survival (PFS), and disease-specific survival (DSS) between CM patients with high or low biomarker expression. We applied univariable and multivariable Cox analyses to verify the independent prognostic value of the FAM biomarkers. We used GSEA and TME analysis to investigate the immune-related regulation mechanism of the FAM subtype biomarker. We evaluated the immune checkpoint inhibition (ICI) response and chemotherapy sensitivity between CM patients with high or low biomarker expression. We performed real-time fluorescent quantitative PCR (qRT-PCR) and semi-quantitative analysis of the immunohistochemical (IHC) data from the Human Protein Atlas to evaluate the mRNA and protein expression levels of the FAM biomarkers in CM.ResultsWe identified 2 FAM molecular subtypes (cluster 1 and cluster 2). K-M analysis showed that cluster 2 had better OS and PFS than cluster 1 did. GSEA showed that, compared with cluster 1, cluster 2 had significantly upregulated immune response pathways. The TME analysis indicated that immune cell subpopulations and immune functions were highly enriched in cluster 2 as compared with cluster 1. WGCNA identified 6 hub genes (ACSL5, ALOX5AP, CD1D, CD74, IL4I1, and TBXAS1) as FAM biomarkers. CM patients with high expression levels of the six biomarkers had better OS, PFS, and DSS than those with low expression levels of the biomarkers. The Cox regression analyses verified that the 6 FAM biomarkers can be independent prognostic factors for CM patients. The single-gene GSEA showed that the high expression levels of the 6 genes were mainly enriched in T-cell antigen presentation, the PD-1 signaling pathway, and tumor escape. The TME analysis confirmed that the FAM subtype biomarkers were not only related to immune infiltration but also highly correlated with immune checkpoints such as PD-1, PD-L1, and CTLA-4. TIDE scores confirmed that patients with high expression levels of the 6 biomarkers had worse immunotherapy responses. The 6 genes conveyed significant sensitivity to some chemotherapy drugs. qRT-PCR and IHC analyses verified the expression levels of the 6 biomarkers in CM cells.ConclusionOur FAM subtypes verify that different FAM reprogramming affects the function and phenotype of infiltrating immune cells in the CM TME. The FAM molecular subtype biomarkers can be independent predictors of prognosis and immunotherapy response in CM patients.