Construction and validation of a novel pyroptosis-related signature to predict prognosis in patients with cutaneous melanoma

Multiomics characterization of pyroptosis in the tumor microenvironment and therapeutic relevance in metastatic melanoma

BACKGROUND

Pyroptosis, mediated by gasdermins with the release of multiple inflammatory cytokines, has emerged as playing an important role in targeted therapy and immunotherapy due to its effectiveness at inhibiting tumor growth. Melanoma is one of the most commonly used models for immunotherapy development, though an inadequate immune response can occur. Moreover, the development of pyroptosis-related therapy and combinations with other therapeutic strategies is limited due to insufficient understanding of the role of pyroptosis in the context of different tumor immune microenvironments (TMEs).

METHODS

Here, we present a computational model (pyroptosis-related gene score, PScore) to assess the pyroptosis status. We applied PScore to 1388 melanoma samples in our in-house cohort and eight other publicly available independent cohorts and then calculated its prognostic power of and potential as a predictive marker of immunotherapy efficacy. Furthermore, we performed association analysis for PScore and the characteristics of the TME by using bulk, single-cell, and spatial transcriptomics and assessed the association of PScore with mutation status, which contributes to targeted therapy.

RESULTS

Pyroptosis-related genes (PRGs) showed distinct expression patterns and prognostic predictive ability in melanoma. Most PRGs were associated with better survival in metastatic melanoma. Our PScore model based on genes associated with prognosis exhibits robust performance in survival prediction in multiple metastatic melanoma cohorts. We also found PScore to be associated with BRAF mutation and correlate positively with multiple molecular signatures, such as KRAS signaling and the IFN gamma response pathway. Based on our data, melanoma with an immune-enriched TME had a higher PScore than melanoma with an immune-depleted or fibrotic TME. Additionally, monocytes had the highest PScore and malignant cells and fibroblasts the lowest PScore based on single-cell and spatial transcriptome analyses. Finally, a higher PScore was associated with better therapeutic efficacy of immune checkpoint blockade, suggesting the potential of pyroptosis to serve as a marker of immunotherapy response.

CONCLUSIONS

Collectively, our findings indicate that pyroptosis is a prognostic factor and is associated with the immune response in metastatic melanoma, as based on multiomics data. Our results provide a theoretical basis for drug combination and reveal potential immunotherapy response markers.

Emerging prognostic biomarkers in advanced cutaneous melanoma: a literature update

INTRODUCTION

Over the past two years, the scientific community has witnessed an exponential growth in research focused on identifying prognostic biomarkers for melanoma, both in pre-clinical and clinical settings. This surge in studies reflects the need of developing effective prognostic indicators in the field of melanoma.

AREAS COVERED

The aim of this work is to review the scientific literature on the most recent findings on the development or validation of prognostic biomarkers in melanoma, in the attempt of providing both clinicians and researchers with an updated broad synopsis of prognostic biomarkers in cutaneous melanoma.

EXPERT OPINION

While the field of prognostic biomarkers in melanoma appears promising, there are several complexities and limitations to address. The interdependence of clinical, histological, and molecular features requires accurate classification of different biomarker families. Correlation does not imply causation, and adjustments for confounding factors are often overlooked. In this scenario, large-scale studies based on high-quality clinical trial data can provide more reliable evidence. It is essential to avoid oversimplification by focusing on a single biomarker, as the interactions among multiple factors contribute to define the disease course and patient's outcome. Furthermore, implementing well-supported evidence in real-life settings can help advance prognostic biomarker research in melanoma.

Novel pyroptosis-related lncRNAs and ceRNAs predict osteosarcoma prognosis and indicate immune microenvironment signatures

Objective

To study pyroptosis-related biomarkers that are associated with the prognosis and immune microenvironment characteristics of osteosarcoma (OS). The goal is to establish a foundation for the prognosis and treatment of OS.

Methods

We retrieved transcriptome and clinical data from The Cancer Genome Atlas (TCGA) database for 88 OS patients. Using this data, we constructed a prognostic model to identify pyroptosis-related genes (PRGs) associated with OS prognosis. To further explore the biological function of these PRGs, we performed enrichment analysis. To identify pyroptosis-related long non-coding RNAs (PRLncs) associated with the prognosis of OS, we performed co-expression analysis. Subsequently, a risk prognostic model was constructed using these PRLncs to generate a risk score, termed as PRLncs-score, thereby obtaining PRLncs associated with the prognosis of OS. The accuracy of the prognostic model was verified through survival analysis, risk curve, independent prognostic analysis, receiver operating characteristic (ROC) curve, difference analysis between high- and low-risk groups, and clinical correlation analysis. And to determine whether PRLncs-score is independent prognostic factor for OS. In addition, we further conducted external and internal validation for the risk prognosis model. Further analyses of immune cell infiltration and tumor microenvironment were performed. A pyroptosis-related competitive endogenous RNA (PRceRNA) network was constructed to obtain PRceRNAs associated with the prognosis of OS and performed gene set enrichment analysis (GSEA) on PRceRNA genes.

Results

We obtained five PRGs (CHMP4C, BAK1, GSDMA, CASP1, and CASP6) that predicted OS prognosis and seven PRLncs (AC090559.1, AP003119.2, CARD8-AS1, AL390728.4, SATB2-AS1, AL133215.2, and AC009495.3) and one PRceRNA (CARD8-AS1-hsa-miR-21-5p-IL1B) that predicted OS prognosis and indicated characteristics of the OS immune microenvironment. The PRLncs-score, in combination with other clinical features, was established as an independent prognostic factor for OS patients. Subsequent scrutiny of the tumor microenvironment and immune infiltration indicated that patients with low-PRLncs-scores were associated with reduced metastatic risk, improved survival rates, heightened levels of immune cells and stroma, and increased immune activity compared to those with high-PRLncs-scores.

Conclusion

The study's findings offer insight into the prognosis of OS and its immune microenvironment, and hold promise for improving early diagnosis and immunotherapy.

Clinical-mediated discovery of pyroptosis in CD8+ T cell and NK cell reveals melanoma heterogeneity by single-cell and bulk sequence

Histologically, melanoma tissues had fewer positive cells percentage of pyroptosis-related genes (PRGs), GZMA, GSDMB, NLRP1, IL18, and CHMP4A in epidermal than in normal skin. Pyroptosis, a new frontier in cancer, affects the tumor microenvironment and tumor immunotherapy. Nevertheless, the role of pyroptosis remains controversial, which reason is partly due to the heterogeneity of the cellular composition in melanoma. In this study, we present a comprehensive analysis of the single-cell transcriptome landscape of pyroptosis in melanoma specimens. Our findings reveal dysregulation in the expression of PRGs, particularly in immune cells, such as CD8+ cells (representing CD8+ T cells) and CD57+ cells (representing NK cells). Additionally, the immunohistochemical and multiplex immunofluorescence staining experiments results further confirmed GZMA+ cells and GSDMB+ cells were predominantly expressed in immune cells, especially in CD8 + T cells and NK cells. Melanoma specimens secreted a minimal presence of GZMA+ merged CD8+ T cells (0.11%) and GSDMB+ merged CD57+ cells (0.08%), compared to the control groups exhibiting proportions of 4.02% and 0.62%, respectively. The aforementioned findings indicate that a reduced presence of immune cells within tumors may play a role in diminishing the ability of pyroptosis, consequently posing a potential risk to the anti-melanoma properties. To quantify clinical relevance, we constructed a prognostic risk model and an individualized nomogram (C-index=0.58, P = 0.002), suggesting a potential role of PRGs in malignant melanoma prevention. In conclusion, our integrated single-cell and bulk RNA-seq analysis identified immune cell clusters and immune gene modules with experiment validation, contributing to our better understanding of pyroptosis in melanoma.

A comprehensive signature based on endoplasmic reticulum stress-related genes in predicting prognosis and immunotherapy response in melanoma

Melanoma is considered as one of the most invasion types of skin cancer with high mortality rates. Although combination of immune checkpoint therapy with local surgical excision provide a novel promising therapeutic strategies, the overall prognosis of melanoma patients remains unsatisfactory. Endoplasmic reticulum (ER) stress, a process of protein misfolding and undue accumulation, has been proven to play an indispensable regulatory role in tumor progression and tumor immunity. However, whether the signature based ER genes has predictive value for the prognosis and immunotherapy of melanoma has not been systematically manifested. In this study, the LASSO regression and multivariate Cox regression were applied to construct a novel signature for predicting melanoma prognosis both in the training and testing set. Intriguingly, we found that patients endowed with high- and low-risk scores displayed differences in clinicopathologic classification, immune cell infiltration level, tumor microenvironment, and immune checkpoint treatment response. Subsequently, based on molecular biology experiments, we validated that silencing the expression of RAC1, an ERG composed of the risk signature, could restrain the proliferation and migration, promote apoptosis, as well as increase the expression of PD-1/PD-L1 and CTLA4 in melanoma cells. Taken together, the risk signature was regarded as promising predictors for melanoma prognosis and might provide prospective strategies to ameliorate patients' response to immunotherapy.

The Therapeutic Potential of Pyroptosis in Melanoma

Pyroptosis is a programmed cell death characterized by the rupture of the plasma membranes and release of cellular content leading to inflammatory reaction. Four cellular mechanisms inducing pyroptosis have been reported thus far, including the (i) caspase 1-mediated canonical, (ii) caspase 4/5/11-mediated non-canonical, (iii) caspase 3/8-mediated and (iv) caspase-independent pathways. Although discovered as a defense mechanism protecting cells from infections of intracellular pathogens, pyroptosis plays roles in tumor initiation, progression and metastasis of tumors, as well as in treatment response to antitumor drugs and, consequently, patient outcome. Pyroptosis induction following antitumor therapies has been reported in several tumor types, including lung, colorectal and gastric cancer, hepatocellular carcinoma and melanoma. This review provides an overview of the cellular pathways of pyroptosis and discusses the therapeutic potential of pyroptosis induction in cancer, particularly in melanoma.

A combined aging and immune prognostic signature predict prognosis and responsiveness to immunotherapy in melanoma

Background: Melanoma is the most lethal, and one of the most aggressive forms of cutaneous malignancies, which poor response to treatment has always puzzled clinicians. As is known to all, aging and immune microenvironment are two crucial factors impacting melanoma biological progress through the tumor microenvironment (TME). However, reliable biomarkers for predicting melanoma prognosis based on aging and immune microenvironment and therapeutic efficacy of immune checkpoints remain to be determined.

Methods: The aging-related genes (ARGs) were obtained from the Human Ageing Genomic Resources and immune-related genes (IRGs) were downloaded from the Immunology database as well as Analysis Portal (ImmPort) database. Next, we initially performed LASSO regression and multivariate Cox regression to identify prognostic ARGs and IRGs in the TCGA and GSE65904 datasets, and firstly constructed a novel comprehensive index of aging and immune (CIAI) signature. Finally, in vitro molecular biology experiments were performed to assess the regulatory role of CNTFR in melanoma cell lines proliferation and migration, macrophage recruitment, and M2 polarization.

Results: This novel CIAI signature consisted of 7 genes, including FOXM1, TP63, ARNTL, KIR2DL4, CCL8, SEMA6A, and CNTFR, in which melanoma patients in the high-CIAI group had shorter OS, DSS, and PFI, indicating CIAI model served as an independent prognostic index. Moreover, we found the CIAI score was potentially correlated with immune scores, estimate score, immune cell infiltration level, tumor microenvironment, immunotherapy effect, and drug sensitivity. Finally, CNTFR might function as oncogenes in melanoma cell lines and the silencing of CNTFR reduced macrophage recruitment and M2 polarization.

Conclusion: In this study, we have first presented a novel prognostic CIAI model applied to assess immune checkpoint therapy and the efficacy of conventional chemotherapy agents in melanoma patients. Thus providing a new insight for combating melanoma.

Melanoma Molecular Subtypes and Development of Prognostic and Immunotherapy-Related Genetic Characteristics by Ferroptosis Gene Analysis

The dissimilarity is a major problem in clinical therapy of skin cutaneous melanoma (SKCM). Objective and reproducible classification systems may help decode SKCM heterogeneity. ConsensusClusterPlus was used to establish a stable immune molecular classification based on ferroptosis-related genes that had been acquired from FerrDb. Moreover, the prognosis, somatic mutations, immune microenvironment characteristics, functional enrichment, and clinical responsiveness to the immune checkpoint blockade of different subtypes in two independent melanin datasets were compared. Kaplan-Meier curves, univariate, multivariate, least absolute contraction, and selection operator (LASSO) Cox regression analysis were used to develop a molecular model for predicting survival, which was verified by a nomogram on the basis of independent prognostic indicators. Two molecular subtypes (C1 and C2) for SKCM were first identified according to ferroptosis-related genes; C1 showed a poor prognosis, with lower infiltration degree of immune cells and TIED score and higher homologous recombination defects, fraction altered, the number of segments, and copy number amplification and deletion. These characteristics of C2 were the opposite of C1. A ferroptosis-related prognosis risk score (FPRS) model was constructed using 6 of 463 genes with differential expression between C1 and C2. This model splits patients into low- and high-risk cohorts. There were significant differences in the infiltration and proportion of immune cells, immune checkpoint gene expression, responsiveness to immune checkpoint therapy, and sensitivity to chemotherapeutic medications between low- and high-risk cohorts. This model was an independent prognostic marker for SKCM and has a high AUC. In summary, we have identified two subtypes of SKCM with different molecular and immune characteristics on the basis of ferroptosis-related genes and further developed and verified an FPRS model, which might independently serve as a prognostic marker for SKCM.

Pan-cancer analyses of pyroptosis with functional implications for prognosis and immunotherapy in cancer

Background

Programmed cell death is an active and orderly form of cell death regulated by intracellular genes that plays an important role in the normal occurrence and development of the immune system, and pyroptosis has been found to be involved in tumorigenesis and development. However, compressive analysis and biological regulation of pyroptosis genes are lacking in cancers.

Methods

Using data from The Cancer Genome Atlas, we established a score level model to quantify the pyroptosis level in cancer. Multiomics bioinformatic analyses were performed to assess pyroptosis-related molecular features and the effect of pyroptosis on immunotherapy in cancer.

Results

In the present study, we performed a comprehensive analysis of pyroptosis and its regulator genes in cancers. Most pyroptosis genes were aberrantly expressed in different types of cancer, attributed to the CAN frequency and differences in DNA methylation levels. We established a pyroptosis level model and found that pyroptosis had dual roles across cancers, while the pyroptosis levels were different among multiple cancers and were significantly associated with clinical prognosis. The dual role of pyroptosis was also shown to affect immunotherapeutic efficacy in several cancers. Multiple pyroptosis genes showed close associations with drug sensitivity across cancers and may be considered therapeutic targets in cancer.

Conclusions

Our comprehensive analyses provide new insight into the functions of pyroptosis in the initiation, development, progression and treatment of cancers, suggesting corresponding prognostic and therapeutic utility.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03313-x.