- Open Access
A novel inflammatory response-related signature predicts the prognosis of cutaneous melanoma and the effect of antitumor drugs
World Journal of Surgical Oncology volume 20, Article number: 263 (2022)
Cutaneous melanoma (CM) is a skin cancer that is highly metastatic and aggressive, with a dismal prognosis. This is the first study to use inflammatory response-related genes to build a model and evaluate their predictive significance in CM. This study used public databases to download CM patients’ mRNA expression profiles and clinical data to create multigene prognostic markers in the UCSC cohort. We compared overall survival (OS) between high- and low-risk groups using the Kaplan-Meier curve and determined independent predictors using Cox analysis. We also used enrichment analysis to assess immune cell infiltration fraction and immune pathway-related activity using KEGG enrichment analysis. Furthermore, we detected prognostic genes’ mRNA and protein expression in CM and normal skin tissues using qRT-PCR and immunohistochemistry. Finally, we developed a 5-gene predictive model that showed that patients in the high-risk group had a considerably shorter OS than those in the low-risk group. The analysis of the receiver operating characteristic (ROC) curve proved the model’s predictive ability. We also conducted a drug sensitivity analysis and discovered that the expression levels of prognostic genes were substantially linked with cancer cell sensitivity to antitumor medicines. The findings show that the model we developed, which consists of five inflammatory response-related genes, can be used to forecast the prognosis and immunological state of CM, giving personalized and precision medicine a new goal and direction.
Melanoma is a highly aggressive and invasive tumor that arises from melanocytes, often in the skin and mucous membranes . According to the American Joint Committee on Cancer (AJCC), between 2018 and 2021, there will be 91,270 and 106,110 new melanoma cases, respectively [2, 3]. According to statistics, the global incidence of melanoma was around 3.9/100,000 in 2017, representing a 41.2% increase since 1990 . Despite accounting for a small percentage of skin cancers, cutaneous melanoma (CM) has a significant recurrence rate, mortality, and medication resistance compared to other skin malignancies [5,6,7]. After complete resection, early-stage CM has a high cure rate, with a mortality rate of 70% and a 5-year survival rate of fewer than 16% in stages 3/4 patients . As a result, novel CM biomarkers must be investigated to guide clinical treatment and enhance the prognosis of CM patients. Gene-based biomarkers have gotten much attention since human gene sequencing technologies developed .
Researchers have been studying the role of inflammation in the genesis and progression of cancer [10,11,12]. Inflammation can both promote and prevent cancer [13, 14]. The link between cancer and indicators of the inflammatory response can be investigated by evaluating basic blood measurements. Many inflammatory response-related characteristics in the peripheral blood of patients with CM, such as thrombocytosis, leukocytosis, hypoproteinemia, and increased plasma fibrin, have been established in investigations [15, 16]. In the overall survival analysis of cancer, clinical systemic inflammatory markers such as neutrophil ratio, platelet-lymphocyte ratio, and lymphocyte-monocyte ratio revealed substantial predictive potential independent of previously identified prognostic variables for CM . The link between inflammatory response-related genes and CM prognosis, on the other hand, is uncertain .
In this study, we utilized the UCSC database to download the mRNA expression profiles of CM patients as well as the corresponding clinical data, used differential analysis to find the inflammatory response-related genes differentially expressed in CM, built a CM prognostic marker model, and validated the model’s stability and reliability using training and validation sets. We used functional enrichment analysis to investigate the model’s potential mechanisms of action and looked at the association between predictive gene expression and the type of immune infiltration. In addition, we looked at the relationship between prognostic gene expression and tumor chemoresistance, and our findings were then confirmed using qRT-PCR and IHC (data from the Human Protein Atlas), resulting in novel concepts for predicting CM prognosis.
Material and methods
We used the UCSC database (https://xena.ucsc.edu/) to find the GDC TCGA melanoma dataset and downloaded the HTSeq-fpkm dataset (n = 472) as the tumor group. As the control group (n = 323), we looked for the GTEX dataset on UCSC, retrieved the TOIL RSEM fpkm dataset, and extracted the skin data from it. We next joined the tumor and control groups to create our gene expression matrix. In addition, we obtained data on inflammation-related genes (n = 200) from the GSEA database (http://www.gsea-msigdb.org/). For the following study analysis, we shall merge the above datasets.
Screening for differential genes (DEGs) related to inflammatory response
We used the limma package in R language to identify differential genes (DEGs) associated with inflammatory responses, and we set the value of corFilter ≥ 0.4 and pvaluefilter ≤ 0.01. In the search for differential genes (DEGs), we set the significance threshold to |log2FC| ≥ 2 and the false-positive rate (FDR) ≤ 0.05.
Modeling of inflammatory response-related prognosis in the training and validation groups
Overall survival (OS) was used as a clinical endpoint in our study, and to improve the model’s accuracy, we screened patients with prognostic information in the tumor group patient group and then divided them into a training group (n = 228) and a validation group (n = 226). We then used one-way Cox analysis to screen for overall survival (OS) associated with inflammatory differential genes (n = 15) and then adjusted P-values using the Benjamini and Hochberg (BH) method. Hazard ratios (HR) were utilized to assess if an inflammatory factor was a risk factor (HR > 1) or a protective factor (HR < 1).
We used the least absolute shrinkage and selection operator (LASSO) Cox regression analysis to develop prognostic models for fitting overall survival (OS) in CM patients based on these inflammatory differential genes. We used the R language’s glmnet package to select and shrink variables such that some of the regression coefficients were strictly equal to 0 to obtain the best prognostic model . We used the normalized expression matrix of potential prognostic factors identified by DEGs as the independent variable and the overall survival (OS) and status of patients as the dependent variable in the regression analysis. We used tenfold cross-validation to determine the prediction model’s penalty parameter (λ) and followed the minimum criterion (the value of λ corresponds to the lowest likelihood deviation). We produced patient risk scores and risk profile models based on the expression levels of important inflammatory response genes and their accompanying regression coefficients. The risk score is calculated as follows: β (hub gene1) × expression (hub gene1) + β (hub gene2) × expression (hub gene2) + β (hub genen) × expression (hub genen)
We classified patients in the training and validation groups into low- and high-risk groups using the risk score formula and the median score value of the risk score in the training group as the cutoff point. We then used the R package survminer to plot Kaplan-Meier survival curves for overall survival (OS) in low- and high-risk groups, as well as ROC (receiver operating characteristic) curves to assess the model’s prediction usefulness and accuracy. We used the rtsne and ggplot2 packages in R to perform principal component analysis (PCA)  and t-SNE  analysis on the gene expression levels in the built models to minimize the dimensionality and analyze the distribution of the low- and high-risk groups.
Independent prognostic analysis
We employed univariate and multivariate Cox regression to determine if our prognostic model was independent of other commonly used clinical variables (e.g., age, gender) associated with OS in CM patients. Independent prognosis analyses were performed on all sample groups, the training group, and the validation group (P < 0.05). To predict the chance of survival in CM patients, we utilized a multifactor logistic model that included a risk score model and traditional clinical data (including age and gender). We used the model and examined its accuracy in the training and validation groups.
Functional enrichment: GSEA analysis
Gene set enrichment analysis (GSEA) ranks genes based on their differential expression in two types of samples and then determines whether the preset gene set is enriched at the top or bottom of this ranking table. We analyzed inflammatory response-related prognostic genes across high- and low-risk groups using GSEA analysis (version 4.2.0) to identify differential KEGG signaling networks and putative biological processes influencing tumor growth . We set the number of permutations to 1000 and the type of permutation to phenotypic, with a significance threshold of P < 0.05. The immune infiltration status of patients with CM and the immunological differences between high- and low-risk groups were determined using single-sample gene set enrichment analysis (ssGSEA). We calculated the infiltration fraction of 16 immune cells and the activation of 13 immune-related pathways in the high-risk and low-risk groups using the R language package GSVA.
Tumor microenvironment (TME) and immune scoring
Malignant tissues include not just tumor cells but also normal epithelial and stromal cells, immune cells, and vascular cells associated with tumors. Infiltrating stromal and immune cells, on the other hand, constitute the majority of normal cells in tumor tissue and play a critical role in tumor biology. We typically utilize immune, stem cell, and stromal scores to quantify the degree of immune and mesenchymal stromal cell infiltration in various tumor tissues. Spearman correlation analysis determined the link between risk and various scores. To examine the link between risk score and immune infiltration subtypes, we utilized a two-way analysis of variance (2-way ANOVA). Tumor stem cell characteristics are frequently employed to quantify the stem cell-like properties of tumors, and the link between tumor stemness and risk score was analyzed using Spearman correlation. The ESTIMATE algorithm in the R language ESTIMATE package was used to estimate the rate of the immune stromal component in TME for each sample, which is reported as three scores: stemnessScore, immuneScore, and stromalScore. Additionally, we employed computational approaches such as ESTIMATE , TIMER , MCP-counter , CIBERSORTx , and ssGSEA in the training group to synthesize immunological differences in the risk model using the R language’s pheatmap package to visualize the graphs.
Drug sensitivity analysis
To study the association between inflammation-related prognostic gene expression and drug sensitivity, we used the CellMiner database (https://discover.nci.nih.gov/cellminer/home.do) for drug sensitivity prediction and Pearson correlation analysis. After selecting the processed data set and downloading the RNA expression data and drug data (compound activity: DTP NCI-60), we read the drug-related data and screened the drug criteria, and after completing the preparation of the gene expression data, we completed the drug sensitivity analysis using the R packages impute and limma. Finally, we visualized and completed the scatter plot using the R packages ggplot2 and ggpubr.
Validation of key genes in CM and paraneoplastic tissues
Following permission from the Chinese PLA General Hospital’s Human Research Ethics Committee, we collected eight pairs of CM and para-cancerous normal tissue specimens. All patients were informed and signed an informed consent form. The relative expression of five key genes was determined using qRT-PCR. We isolated total RNA from cancer and paraneoplastic tissues using TRIzol and determined the quantity of RNA using a NanoDrop spectrophotometer. We then reverse-transcribed the RNA into cDNA using a 5 × RT Master Mix (BioRad) and analyzed mRNA expression levels using a 2 × SYBR Green PCR Kit. The expression of each gene was standardized to GAPDH. Using the 2−∆∆Ct technique, we quantified the real-time PCR analysis and determined the relative expression of genes linked with the validation reaction. The following primer sequences were used: C3AR1: forward, AAG CCA ATC TGG TGT CAG AAT C; reverse, CAG GAA TGC ACA TCA CAA AAG C; CXCL10: forward, GTG GCA TTC AAG GAG TAC CTC; reverse, TGA TGG CCT TCG ATT CTG GAT T; EIF2AK2: forward, GCC GCT AAA CTT GCA TAT CTT CA; reverse, TCA CAC GTA GTA GCA AAA GAA CC; EMP3: forward, CCT GAA TCT CTG GTA CGA CTG C; reverse, GCC ATT CTC GCT GAC ATT ACT G; ICAM1: forward, ATG CCC AGA CAT CTG TGT CC; and reverse, GGG GTC TCT ATG CCC AAC AA. The Huada Corporation synthesized all our primers (Beijing, China). We next used the online database HPA (https://www.proteinatlas.org/) to confirm the differential expression of immunohistochemical proteins associated with important genes in normal tissues and CM.
The Wilcoxon test was used to compare DEGs in CM and para-cancerous normal tissue, and the chi-square test was used to compare various proportions. The Mann-Whitney test was used to compare the ssGSEA scores of immune cells or immunological pathways across high- and low-risk groups, and the P-values were adjusted using the BH technique. We analyzed variations in OS between groups using Kaplan-Meier analysis and screened for independent determinants of OS using univariate and multivariate Cox analysis. The correlation study of prognostic model risk scores or prognostic gene expression levels with stemnessScore, immuneScore, stromalScore, and drug sensitivity was performed using Spearman and Pearson correlation analysis. All statistical analyses were conducted using the R programming language (version 4.1.1). All statistical tests were two sided, and a significance level of P < 0.05 was considered statistically significant.
Acquisition of genes related to inflammatory response
The study’s flow chart is depicted in Fig. 1. The study population included 472 patients with CM from the UCSC cohort and 323 patients with normal skin tissue expression data. The clinical characteristics of these patients are summarized in Table 1 (including age, gender, stage, and TNM stage). The screening of prognostic genes associated with inflammation was conducted with data on 200 prognostic genes associated with inflammation received from the GSEA database (Additional file 1).
Construction and validation of a predictive model related to inflammatory response
We removed samples with identical ID values and missing clinical data, leaving 454 samples with survival days and survival status. These samples were randomly divided into a training group (n = 228) (Additional file 2) and a validation group (n = 226) (Additional file 3). We initially screened using one-way Cox analysis and identified 15 inflammatory genes linked with OS (P < 0.01) (Fig. 2A) (Additional file 4) and demonstrated their correlation (Fig. 2B). We then performed LASSO regression analysis on these genes and screened inflammatory genes with more than 900 replicates in 1000 replacement samples. We then constructed a model containing five genes to predict the prognosis of CM patients based on their λ values (Fig. 3). The model consists of five key genes, including C3AR1, CXCL10, EIF2AK2, EMP3, and ICAM1. The hazard ratio (HR), 95% confidence interval (95% CI), and P-value of five hub gene could be seen at Table 2. The risk score is calculated as follows: (−0.00195702921106954) × expression (C3AR1) + (−0.0717393323519231) × expression (CXCL10) + (−0.178228168474953) × expression (EIF2AK2) + (0.0821050910327823) × expression (EMP3) + (−0.000851776035203931) × expression (ICAM1).
We split all patients in the training and validation groups into high- and low-risk groups based on the training group’s median risk score. In the training group, we discovered that patients’ chance of death increased, and their survival time reduced (Fig. 4A). Scatter plot analysis revealed that patients classified as high risk had a greater likelihood of dying sooner than patients classified as low risk (Fig. 4B). Kaplan-Meier curves revealed that OS was substantially worse in high-risk individuals than in low-risk patients (P < 0.01) (Fig. 4C). The risk heat map depicted the expression of many inflammatory genes associated with prognosis in high- and low-risk groups (Fig. 4D). We next revalidated our results in the validation group, where we noticed substantial variations in OS between the high-risk and low-risk groups (Fig. 4E–G) and inflammatory prognostic gene expression in the heat map (Fig. 4H) proving the model’s accuracy once again. We then used principal component analysis (PCA) (Fig. 5A–B) and t-SNE (Fig. 5C–D) analysis to demonstrate the discrete distribution of patients into different risk groups. Because patients in the high-risk group may die earlier and have shorter survival times than those in the low-risk group, we can intuitively assume that the model constructed using inflammatory prognosis-related genes can better differentiate the prognosis of patients with CM.
Independent prognostic analysis of OS
We utilized univariate and multivariate Cox regression analyses to determine if clinical parameters (age, gender), and risk score was significant independent predictors of OS. We showed that risk score, gender, and age were independent predictive predictors of OS in univariate and multivariate Cox regression analyses in the training group (Fig. 6A–B), validation group (Fig. 6C–D), and all samples (Fig. 6E–F). Additionally, the ROC areas at 1, 2, and 3 years were 0.754, 0.660, 0.654 (Fig. 7A), and 0.670, 0.709, 0.682 (Fig. 7B), respectively, for the training and validation groups. Our findings indicated that the risk score model was considerably better in predicting OS in patients with CM than other clinical factors, such as age and gender (Fig. 7C–D).
Immune status and tumor microenvironment analysis
To further synthesize the immune cell infiltration in CM patients and to explore the differences in immune status between high- and low-risk groups, we applied TIMER, CIBERSORT, CIBERSORT-ABS, QUANTISEQ, MCPCOUNTER, XCELL, and EPIC in the training group, and the results were presented in a heat map (Fig. 8). We then used ssGSEA to quantify the enrichment scores of different immune cell subpopulations, related functions, and pathways to explore the differences in immune cell enrichment scores between the high-risk and low-risk groups. In the training set, we found that the antigen presentation process contained significantly higher levels of aDCs, pDCs, APC-co inhibition, APC-co stimulation, HLA, and MHC class 1 in the low-risk group compared to the high-risk group. In addition, we found that CD8+T cells, T helper cells, Tfh cells, Th1 cells, Th2 cells, TIL cells, Treg cells, T-cell co-inhibition, and T-cell co-stimulation associated with T-cell regulation were all significantly higher in the low-risk group than in the high-risk group. The T helper cells, Tfh cells, Th1 cells, Th2 cells, TIL cells, Treg cells, T-cell co-inhibition, and T-cell co-stimulation were all significantly higher than those in the high-risk group, suggesting that there were differences in T-cell regulation between the low-risk and high-risk groups. In addition, we found that the score of CCR and checkpoint, cytolytic activity, B cells, NK cells, neutrophils, inflammation promoting, para-inflammation, type 1 IFN response, and other biological processes and cellular contents were significantly higher in the low-risk group. Response and other biological processes and cellular content were significantly higher than those in the high-risk group (Fig. 9A and C). We found similar results in the validation set as in the training set (Fig. 9B and D).
To examine the association between risk scores and immunological components in further detail, we conducted a correlation analysis between risk scores and immune infiltration. Six distinct immune infiltration types were identified in human tumors, including C1 (wound healing), C2 (IFN-γ dominant), C3 (inflammatory), C4 (lymphocyte depleted), C5 (immunologically quiet), and C6 (TGF-β dominant), although the C5 and C6 immune subtypes were excluded due to their absence in CM patients. We evaluated immune infiltration of CM in the training set and connected it with the risk score, finding that a high-risk score was strongly associated with C3, whereas a low-risk score was significantly associated with C2 (Fig. 9E).
Tumor stemness can be determined using RNA stemness based on mRNA expression (RNAss) or DNA stemness based on DNA methylation pattern (DNAss). Additionally, we evaluated tumor immunological microenvironment using the immuneScore and stromalScore. Correlation analysis examined the relationship between risk score and tumor stemness and immune microenvironment. The results indicated that the risk score was not significantly correlated with DNAss or RNAss but was significantly negatively correlated with immuneScore and stromalScore (P < 0.05) (Fig. 9F).
Biological function and pathway analysis
We used GSEA to conduct enrichment analysis on KEGG pathways in the high- and low-risk group (Fig. 10A). The results analysis revealed that the five pathways with the highest content in the high-risk group were Alzheimer’s disease, aminoacyl tRNA biosynthesis, chemokine signaling pathway, cytokine receptor interaction, and Huntington’s disease (P < 0.05, FDR < 25%) (Fig. 10B). The JAK-STAT signaling route, oxidative phosphorylation, RNA polymerase, toll-like receptor signaling pathway, and viral myocarditis were the five pathways with the highest levels in the low-risk group (P < 0.05, FDR < 25%) (Fig. 10C).
Prognostic gene expression and sensitivity of melanoma cells to chemotherapy
To make this study more therapeutically relevant, we examined the expression of prognostic genes in NCI-60 cell lines and the association between prognostic gene expression levels and drug sensitivity. The results indicated that all prognostic genes were significantly associated with drug susceptibility to certain chemotherapy agents (P < 0.01). The figure depicts the findings of the drug sensitivity analysis with the highest correlation. CXCL10 expression was associated with increased cancer cell resistance to LDK-378, brigatinib, alectinib, and PF-06463922, among others. C3AR1 expression was associated with higher resistance of cancer cells to denileukin diftitox ontak, isotretinoin, carmustine, estramustine, fluphenazine, nelfinavir, megestrol acetate, alectinib, cyclophosphamide, lomustine, and dromostanolone propiona, among others (Fig. 11).
Validation of key genes in melanoma and paraneoplastic tissues
To confirm the changes in expression of 5 important genes (C3AR1, CXCL10, EIF2AK2, EMP3, ICAM1) between CM and para-cancerous normal tissue, we first investigated their mRNA expression using quantitative real-time PCR (qRT-PCR). qRT-PCR analysis revealed that prognostic genes were expressed at a higher level in melanoma than in normal tissues adjacent to malignancy (P < 0.01) (Fig. 12A–E). We then validated the expression of key genes in CM and normal skin using the Human Protein Atlas (HPA) database; however, because the HPA database lacked immunohistochemical images of normal and tumor tissues for CXCL10 and EMP3, we validated only the expression of the remaining three key genes, which were all significantly different (Fig. 13A–C).
Cutaneous melanoma (CM), a highly malignant tumor that develops on the skin and mucosa, has increased incidence and increased clinical attention in recent years. Although CM treatment and prognosis have improved in recent years, clinical efficacy has not yet reached the predicted level due to the disease’s complicated genetic and molecular pathways. Apart from the conventional gold standard for pathological diagnosis, we frequently struggle with the early diagnosis and prognosis of CM due to the scarcity of melanoma biomarkers. S100B , VEGF (vascular endothelial growth factor) , IGF-1R (the insulin-like growth factor 1 receptor) , Wnt-5a , LDH (lactate dehydrogenase) , and MIA (melanoma inhibitory activity)  have all been shown to correlate positively with prognosis in CM. In contrast, indicators such as RUNX3 (human runt-related transcription factor 3)  have a high degree of accuracy in predicting the prognosis of CM, and traditional inflammation-related serum biomarkers such as neutrophil ratio and lymphocyte/monocyte ratio have also performed well in predicting the prognosis of CM, but no report of inflammatory response-related gene models as prognostic markers of CM has been made. Thus, in this investigation, we attempted to leverage genes associated with the inflammatory response to develop a predictive model for CM patients to enhance overall survival.
This study used multivariate Cox regression analysis to assess five prognostic genes related to inflammation and successfully developed a clinical prediction model for CM. We randomly assigned CM patients to train and validation groups based on median cutoff values and evaluated risk score as an independent risk factor for this prognostic model using univariate and multivariate data analysis of patients in high- and low-risk groups. In univariate and multivariate analyses, we observed that the risk score was more accurate than clinicopathological factors (such as age and gender) in predicting OS in CM patients. The AUC curves demonstrated our model’s excellent predictive performance in training and validation groups. The predictive model created in this work comprises five inflammatory response-related genes: C3AR1, CXCL10, EIF2AK2, EMP3, and ICAM1, all of which are all overexpressed in melanoma tissue and are associated with a poor prognosis. C3AR1 (complement component 3a receptor 1) is a protein-coding gene encodes a component of the complement system’s g protein-coupled transmembrane that spans the C3a receptor . C3AR1 was discovered to promote voluntary exercise behavior in melanoma mice by altering inflammatory and immunological responses in a study conducted by Zhi et al. . Although the effect of C3AR1 expression on melanoma is unknown, a study by Nabizadeh discovered that in the absence of complement C3aR (receptor for complement C3a), the development and growth of B16-F0 melanoma were inhibited in mice. In contrast, C3aR antagonists inhibited the growth of established melanoma . CXCL10 (chemokine ligand 10) is also known as interferon-γ-inducible protein 10 (IP 10). It is a member of the CXC chemokine superfamily and plays a role in immune response regulation, angiogenesis, apoptosis, the cell cycle, and cell proliferation [37, 38]. CXCL10 has been demonstrated to induce apoptosis in HeLa cells by suppressing the production of the HPV oncogenic proteins E6 and E7, consequently promoting the sustained expression of P53 in tumor cells . Additionally, CXCL10 has been shown to decrease estrogen-induced pro-tumor development by inhibiting VEGF production . CXCL10 was discovered to repress tumors and stimulate their development via a variety of mechanisms. Wennerberg et al. discovered that NK cells have a greater capacity to metastasis to melanomas that exhibit positive CXCL10 expression than melanomas that express negative CXCL10 expression, resulting in decreased tumor load and increased survival time . CXCL10 is a powerful angiogenesis inhibitor that binds to CXCR3 receptors and suppresses melanoma’s angiogenesis by decreasing intra-tumor vascular density and boosting apoptosis and necrosis of malignant tissue . EIF2AK2 (eukaryotic translation initiation factors 2 AK2), commonly known as PKR (dsRNA-dependent protein kinase), was initially discovered as a pathogen recognition molecule, with the N-terminal of this protein binding to double-stranded RNA to induce activation of the C-terminal catalytic region [43, 44]. Paola’s study demonstrated that when cancer cells are stimulated with pro-ICD (immunogenic cell death) drugs, PKR, a key mediator of eIF2 phosphorylation, can promote CTR translocation to the surface of melanoma cells by promoting cancer cell death and the release of damage-associated molecular patterns (DAMP) from dead cells . The EMP3 (epithelial membrane protein 3) gene is a member of the TMP22 (peripheral myelin protein) gene family . The EMP3 protein’s primary roles are assumed to be related to cell proliferation, differentiation, activation of the caspase apoptotic pathway, and intercellular connections . Exogenous expression of the EMP3 gene has been found to limit tumor cell proliferation and act as a tumor suppressor gene . It has been demonstrated that mRNAs associated with epithelial cell lineage interactions (EMP3 and EMP1) are frequently overexpressed in uveal melanoma (UM) with a high risk of metastasis, which is thought to promote plasticity in UM cells, thereby increasing their resistance to conventional chemotherapeutic agents [49, 50]. Although ICAM1 (intercellular adhesion molecule 1) is ubiquitous on the cell surface, its overexpression in several cancers promotes tumor cell invasion and metastasis and is inversely connected with patient prognosis . It has been demonstrated that human primary melanoma cells (T1) produce high levels of ICAM1, and that increased ICAM1 expression associated with PI3K/AKT pathway activation can be exploited by metastatic melanoma cells to resist CTL-mediated lysis . ICAM1 has been proposed as a novel target for a range of solid tumors in recent years, including melanoma and mesenchymal thyroid carcinoma (ATC), and has significant research potential .
We employed the GSEA technique to further investigate the biological pathways in high- and low-risk groups. Alzheimer’s disease, aminoacyl tRNA production, chemokine signaling, cytokine receptor interaction, and Huntington’s disease were all significantly enriched in the high-risk group. The aminoacyl-tRNA synthetase (ARS) is a critical enzyme that accurately translates the genetic information contained in messenger RNA into the amino acid sequence of proteins. It has a strong effect on cancer through its effects on apoptosis, RNA splicing, and angiogenesis . T-3861174, a prolyl-tRNA synthetase (PRS) inhibitor, was reported to promote apoptosis in melanoma SK-MEL-2 cells by activating the GCN2-ATF4 pathway, and the apoptosis-inducing activity of T-3861174 was eliminated when GCN2 was knocked down . This indicates that various ARS linked with cancers may be valuable therapeutic targets.
Numerous studies have demonstrated that the immune microenvironment plays a critical role in carcinogenesis, and that invading immune cells can operate as tumor growth promoters and inhibitors. To better understand the relationship between immunotherapy, the immune microenvironment, and neoplastic processes, we used the ssGSEA technique to determine the immune cell infiltration status of CM patients. We compared immune cell infiltration status and tumor microenvironment between high- and low-risk groups. We used correlation analysis to examine the relationship between risk scores and tumor stem cells and the immune microenvironment to understand better the prediction model’s potential mechanisms and predictive performance. CD8+ T cells and T helper cells were abundant in the low-risk group of patients. By secreting perforin, serine esterase, and lymphotoxin, CD8+T (CTL) lymphocytes contribute significantly to tumor clearance. CTL cells recognize target cells via the membrane receptor Ti-CD3 (TCR), which results in guanine and lymphotoxin binding to the CTL membrane. Phospholipase C activation in conjunction with nucleotide-binding proteins completes the activation of second messengers and information transduction, finally resulting in Ca2+ release and protein kinase activation, driving tumor cell death and breakdown . T helper cells are activated when they respond with peptide antigens presented by MHC II (major histocompatibility complex), which regulates or assists the immune response by secreting cytokines against tumor cells . The decreased amounts of immune cells such as CTL and T helper cells in the high-risk group show that immune control is disrupted, and antitumor immune function is weakened in this group, which may contribute significantly to their poor prognosis. As a result, we believe that CTL and T helper cells may be critical therapeutic targets for CM patients.
We explored the function of risk scores in the kind of immune infiltration better to understand the relationship between risk scores and immunological components. We discovered that high-risk scores were strongly linked with C3. In contrast, low-risk scores were significantly connected with C2, implying that C3 as a risk factor may promote carcinogenesis and progression, whilst C2 as a protective factor may retard tumorigenesis and progression. This conclusion is consistent with prior findings that C3 and C4, as immune infiltration types with moderate toxicity, are related to infiltration of suppressive immune cell populations, and patients have a poor prognosis . Many researchers now assume that the tumor microenvironment is formed of a network of stromal cells (fibroblasts, vascular cells, and inflammatory immune cells) . Tumor stem cells accelerate the clinical course of malignancies due to their high capacity for self-renewal and invasion, resulting in treatment resistance . As a result of our research, the high infiltration of tumor immune histiocytes in the high-risk group of patients was negatively connected with the immune and stromal scores, which was consistent with the risk score results.
To make this analysis more practically applicable, we examined the connection between important gene expression levels and chemotherapeutic drug sensitivity using data from the cellMiner database for the NCI-60 cell line. Numerous medications have been found in clinical trials to inhibit melanoma growth. Denileukin diftitox ontak, for example, is a diphtheria toxin-based fusion protein that has been licensed for the treatment of persistent skin malignancies such as cutaneous T-cell lymphoma and melanoma via CD25 depletion . Carmustine functions as a chloroethylating nitrosoureas (CNU) by alkylating DNA bases to generate interstrand cross-links (ICL), which inhibit DNA replication and transcription by covalently connecting DNA strands, resulting in severe cytotoxic DNA damage and hence anticancer activity. Because of its potent antitumor properties, Carmustine has been utilized clinically to treat malignant gliomas and melanomas . As a result of our findings, we hypothesize that some prognostic genes could be employed as therapeutic targets to overcome treatment resistance or adjuvant drug sensitivity. We also validated the model’s accuracy by demonstrating the expression levels of important genes in CM and normal skin tissues using qRT-PCR and IHC.
Since Virchow speculated in 1863 that tumors originate from chronic inflammation, there is growing evidence that inflammation plays an essential role in tumorigenesis, development, and evolution. A variety of cytokines produced during inflammation regulate the activation and migration of endothelial cells and their proliferation, survival, and apoptosis, thus playing an essential role in angiogenesis . From our pathway enrichment and immunoassay results, it is clear that patients in the high-risk group have less inflammation-associated cell infiltration, are in an overall low immune infiltration state, have less ability to target tumor clearance, which has more rapid tumor progression, and are more likely to metastasize to distant sites. Therefore, we hypothesize that in CM, as inflammatory infiltration deepens, tumor risk grade becomes higher, overall immune infiltration content becomes lower, and tumors are more likely to spread, metastasize, and recur. We suggest that the inflammatory environment and TME together constitute the tumor ecosphere. Inflammation in TME can accelerate alterations in epigenetic and TME components, jointly promoting tumor development . Hashimoto et al.’s study reported a rare case of undifferentiated pleomorphic sarcoma (UPS) in a high C-reactive protein (CRP) senior patient . Considering that both UPS and CM are in a state of hyper-immune infiltration compared to normal tissue and can be classified as “hot tumors,” we speculate that wide margin resection as soon as possible after exclusion of bacteremia-like infections may be a new option for the treatment of CM in a severe inflammatory state. Schuckmann et al.’s study demonstrated that people who regularly used anti-inflammatory drugs and statins had a much lower rate of CM than those who did not take such drugs . Dipak et al.’s study showed that ketorolac could indirectly stimulate some T cells of the immune system by preferentially inhibiting the COX-1 enzyme, enhancing the immune checkpoint inhibitor’s effect. Ketorolac can eradicate cancer metastasis and prolong survival in mouse models by administering the drug preoperatively . Therefore, we suggest that anti-inflammatory drugs can alter the prognosis of skin tumors. If used as a preventive effect, these drugs need to be taken regularly in daily life. If one already has CM, preoperative administration has been shown in animal models to eradicate cancer metastases and prolong survival, but further confirmation in clinical studies is needed.
This work creates a predictive model for CM patients based on inflammatory response genes. This study, however, has several drawbacks. First, we employed only the UCSC dataset to assess our model’s predictive ability, and we used limited validation approaches due to a lack of sufficient clinical samples. Our model may encourage others to conduct additional research on inflammatory prognostic variables in CM despite its shortcomings.
In summary, we established a new predictive model based on five inflammatory response genes associated with prognosis (C3AR1, CXCL10, EIF2AK2, EMP3, ICAM1). This prognostic model was demonstrated to be independently linked with overall survival (OS) in the UCSC cohort. It gives new directions and research value in the tumor microenvironment, functional analysis, immune response, and treatment sensitivity. Our research on inflammatory response-related prognostic genes not only elucidates their involvement in cancer but also lays the groundwork for the future development of personalized and precision medicine.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary files.
Wei CY, Zhu MX, Lu NH, Peng R, Yang X, Zhang PF, et al. Bioinformatics-based analysis reveals elevated MFSD12 as a key promoter of cell proliferation and a potential therapeutic target in melanoma. Oncogene. 2019;38:1876–91.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018;68:7–30.
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer Statistics, 2021. CA Cancer J Clin. 2021;71:7–33.
Wu Y, Wang Y, Wang L, Yin P, Lin Y, Zhou M. Burden of melanoma in China, 1990–2017: findings from the 2017 global burden of disease study. Int J Cancer. 2020;147:692–701.
Matias M, Pinho JO, Penetra MJ, Campos G, Reis CP, Gaspar MM. The Challenging Melanoma Landscape: From Early Drug Discovery to Clinical Approval. Cells. 2021;10(11):3088. https://doi.org/10.3390/cells10113088.
Yang L, Xu Y, Yan Y, Luo P, Chen S, Zheng B, Yan W, Chen Y, Wang C. Common Nevus and Skin Cutaneous Melanoma: Prognostic Genes Identified by Gene Co-Expression Network Analysis. Genes (Basel). 2019;10(10):747. https://doi.org/10.3390/genes10100747.
Karimkhani C, Green AC, Nijsten T, Weinstock MA, Dellavalle RP, Naghavi M, et al. The global burden of melanoma: results from the Global Burden of Disease Study 2015. Br J Dermatol. 2017;177:134–40.
Haydu LE, Lo SN, McQuade JL, Amaria RN, Wargo J, Ross MI, et al. Cumulative incidence and predictors of CNS metastasis for patients with American Joint Committee on Cancer 8th Edition Stage III Melanoma. J Clin Oncol. 2020;38:1429–41.
Liu Y, Jing R, Xu J, Liu K, Xue J, Wen Z, et al. Comparative analysis of oncogenes identified by microarray and RNA-sequencing as biomarkers for clinical prognosis. Biomark Med. 2015;9:1067–78.
Balkwill F, Mantovani A. Inflammation and cancer: back to Virchow? Lancet. 2001;357:539–45.
Bauer AT, Gorzelanny C, Gebhardt C, Pantel K, Schneider SW. Interplay between coagulation and inflammation in cancer: limitations and therapeutic opportunities. Cancer Treat Rev. 2021;102:102322.
Grivennikov SI, Greten FR, Karin M. Immunity, inflammation, and cancer. Cell. 2010;140:883–99.
Greten FR, Grivennikov SI. Inflammation and cancer: triggers, mechanisms, and consequences. Immunity. 2019;51:27–41.
Shalapour S, Karin M. Pas de Deux: Control of anti-tumor immunity by cancer-associated inflammation. Immunity. 2019;51:15–26.
Woo YR, Cho SH, Lee JD, Kim HS. The Human Microbiota and Skin Cancer. Int J Mol Sci. 2022;23(3):1813. https://doi.org/10.3390/ijms23031813.
Shao K, Hooper J, Feng H. Racial/Ethnic Health Disparities in Dermatology in the United States Part 2: Disease-specific Epidemiology, Characteristics, Management, and Outcomes. J Am Acad Dermatol. 2022;S0190-9622(22):00195-5. https://doi.org/10.1016/j.jaad.2021.12.062.
Hotz MJ, O'Halloran EA, Hill MV, Hayden K, Zaladonis AG, Deng M, et al. Tumor mutational burden and somatic mutation status to predict disease recurrence in advanced melanoma. Melanoma Res. 2022;32:112–9.
McMillan DC. The systemic inflammation-based Glasgow prognostic score: a decade of experience in patients with cancer. Cancer Treat Rev. 2013;39:534–40.
Gupta S, Lee REC, Faeder JR. Parallel tempering with Lasso for model reduction in systems biology. PLoS Comput Biol. 2020;16:e1007669.
Vidaurre D. A new model for simultaneous dimensionality reduction and time-varying functional connectivity estimation. PLoS Comput Biol. 2021;17:e1008580.
Habowski AN, Habowski TJ, Waterman ML. GECO: gene expression clustering optimization app for non-linear data visualization of patterns. BMC Bioinformatics. 2021;22:29.
Jin M, Li D, Liu W, Wang P, Xiang Z, Liu K. Pinin acts as a poor prognostic indicator for renal cell carcinoma by reducing apoptosis and promoting cell migration and invasion. J Cell Mol Med. 2021;25:4340–8.
Wu S, Liu S, Chen N, Zhang C, Zhang H, Guo X. Genome-wide identification of immune-related alternative splicing and splicing regulators involved in abdominal aortic aneurysm. Front Genet. 2022;13:816035.
Liu Y, Tu M, Wang L. Pan-cancer analysis predicts FOXS1 as a key target in prognosis and tumor immunotherapy. Int J General Med. 2022;15:2171–85.
Chen H, Zhang Z, Zhou L, Cai T, Liu B, Wang L, et al. Identification of CCL19 as a novel immune-related biomarker in diabetic nephropathy. Front Genet. 2022;13:830437.
Wu X, Zhao X, Xiong Y, Zheng M, Zhong C, Zhou Y. Deciphering cell-type-specific gene expression signatures of cardiac diseases through reconstruction of bulk transcriptomes. Front Cell Dev Biol. 2022;10:792774.
Ordóñez NG. Value of melanocytic-associated immunohistochemical markers in the diagnosis of malignant melanoma: a review and update. Hum Pathol. 2014;45:191–205.
Simonetti O, Lucarini G, Rubini C, Lazzarini R, DIPR, Offidani A. Clinical and prognostic significance of survivin, AKT and VEGF in primary mucosal oral melanoma. Anticancer Res. 2015;35:2113–20.
Park SL, Setiawan VW, Kanetsky PA, Zhang ZF, Wilkens LR, Kolonel LN, et al. Serum insulin-like growth factor-I and insulin-like growth factor binding protein-3 levels with risk of malignant melanoma. Cancer Causes Control. 2011;22:1267–75.
Prasad CP, Chaurasiya SK, Axelsson L, Andersson T. WNT-5A triggers Cdc42 activation leading to an ERK1/2 dependent decrease in MMP9 activity and invasive migration of breast cancer cells. Mol Oncol. 2013;7:870–83.
Partl R, Fastner G, Kaiser J, Kronhuber E, Cetin-Strohmer K, Steffal C, et al. KPS/LDH index: a simple tool for identifying patients with metastatic melanoma who are unlikely to benefit from palliative whole brain radiotherapy, Supportive care in cancer : official journal of the Multinational Association of. Support Care Cancer. 2016;24:523–8.
Henry L, Fabre C, Guiraud I, Bastide S, Fabbro-Peray P, Martinez J, et al. Clinical use of p-proteasome in discriminating metastatic melanoma patients: comparative study with LDH, MIA and S100B protein. Int J Cancer. 2013;133:142–8.
Kitago M, Martinez SR, Nakamura T, Sim MS, Hoon DS. Regulation of RUNX3 tumor suppressor gene expression in cutaneous melanoma, Clinical cancer research : an official journal of the American Association for. Cancer Res. 2009;15:2988–94.
Wang M, Busuttil RA, Pattison S, Neeson PJ, Boussioutas A. Immunological battlefield in gastric cancer and role of immunotherapies. World J Gastroenterol. 2016;22:6373–84.
Xia Z, Shang H, Cholewa J, Wang Q, Ding X, Su Q, et al. The effect of exercise on gene expression and signaling in mouse melanoma tumors. Med Sci Sports Exerc. 2020;52:1485–94.
Nabizadeh JA, Manthey HD, Steyn FJ, Chen W, Widiapradja A, Md Akhir FN, et al. The complement C3a receptor contributes to melanoma tumorigenesis by inhibiting neutrophil and CD4+ T cell responses. J Immunol. 1950;196(2016):4783–92.
Giuliani N, Bonomini S, Romagnani P, Lazzaretti M, Morandi F, Colla S, et al. CXCR3 and its binding chemokines in myeloma cells: expression of isoforms and potential relationships with myeloma cell proliferation and survival. Haematologica. 2006;91:1489–97.
Li G, Tian L, Hou JM, Ding ZY, He QM, Feng P, et al. Improved therapeutic effectiveness by combining recombinant CXC chemokine ligand 10 with cisplatin in solid tumors, Clinical cancer research : an official journal of the American Association for. Cancer Res. 2005;11:4217–24.
Zhang HM, Yuan J, Cheung P, Chau D, Wong BW, McManus BM, et al. Gamma interferon-inducible protein 10 induces HeLa cell apoptosis through a p53-dependent pathway initiated by suppression of human papillomavirus type 18 E6 and E7 expression. Mol Cell Biol. 2005;25:6247–58.
Aronica SM, Raiber L, Hanzly M, Kisela C. Antitumor/antiestrogenic effect of the chemokine interferon inducible protein 10 (IP-10) involves suppression of VEGF expression in mammary tissue. J Interf Cytokine Res. 2009;29:83–92.
Wennerberg E, Kremer V, Childs R, Lundqvist A. CXCL10-induced migration of adoptively transferred human natural killer cells toward solid tumors causes regression of tumor growth in vivo. Cancer Immunol Immunother. 2015;64:225–35.
Yang J, Richmond A. The angiostatic activity of interferon-inducible protein-10/CXCL10 in human melanoma depends on binding to CXCR3 but not to glycosaminoglycan. Mol Ther. 2004;9:846–55.
Bullido MJ, Martínez-García A, Tenorio R, Sastre I, Muñoz DG, Frank A, et al. Double stranded RNA activated EIF2 alpha kinase (EIF2AK2; PKR) is associated with Alzheimer's disease. Neurobiol Aging. 2008;29:1160–6.
Lu B, Nakamura T, Inouye K, Li J, Tang Y, Lundbäck P, et al. Novel role of PKR in inflammasome activation and HMGB1 release. Nature. 2012;488:670–4.
Giglio P, Gagliardi M, Tumino N, Antunes F, Smaili S, Cotella D, et al. PKR and GCN2 stress kinases promote an ER stress-independent eIF2α phosphorylation responsible for calreticulin exposure in melanoma cells. Oncoimmunology. 2018;7:e1466765.
Ben-Porath I, Benvenisty N. Characterization of a tumor-associated gene, a member of a novel family of genes encoding membrane glycoproteins. Gene. 1996;183:69–75.
Taylor V, Suter U. Epithelial membrane protein-2 and epithelial membrane protein-3: two novel members of the peripheral myelin protein 22 gene family. Gene. 1996;175:115–20.
Alaminos M, Dávalos V, Ropero S, Setién F, Paz MF, Herranz M, et al. EMP3, a myelin-related gene located in the critical 19q13.3 region, is epigenetically silenced and exhibits features of a candidate tumor suppressor in glioma and neuroblastoma. Cancer Res. 2005;65:2565–71.
Kaochar S, Dong J, Torres M, Rajapakshe K, Nikolos F, Davis CM, et al. ICG-001 exerts potent anticancer activity against uveal melanoma cells. Invest Ophthalmol Vis Sci. 2018;59:132–43.
Onken MD, Ehlers JP, Worley LA, Makita J, Yokota Y, Harbour JW. functional gene expression analysis uncovers phenotypic switch in aggressive uveal melanomas. Cancer Res. 2006;66:4602–9.
Zhang KE, Ge SJ, Lin XY, Lv BB, Cao ZX, Li JM, et al. Intercellular adhesion molecule 1 is a sensitive and diagnostically useful immunohistochemical marker of papillary thyroid cancer (PTC) and of PTC-like nuclear alterations in Hashimoto's thyroiditis. Oncol Lett. 2016;11:1722–30.
Hamaï A, Meslin F, Benlalam H, Jalil A, Mehrpour M, Faure F, et al. ICAM-1 has a critical role in the regulation of metastatic melanoma tumor susceptibility to CTL lysis by interfering with PI3K/AKT pathway. Cancer Res. 2008;68:9854–64.
Wei W, Jiang D, Lee HJ, Li M, Kutyreff CJ, Engle JW, et al. Development and characterization of CD54-targeted immunoPET imaging in solid tumors. Eur J Nucl Med Mol Imaging. 2020;47:2765–75.
Kim S, You S, Hwang D. Aminoacyl-tRNA synthetases and tumorigenesis: more than housekeeping. Nat Rev Cancer. 2011;11:708–18.
Arita T, Morimoto M, Yamamoto Y, Miyashita H, Kitazawa S, Hirayama T, et al. Prolyl-tRNA synthetase inhibition promotes cell death in SK-MEL-2 cells through GCN2-ATF4 pathway activation. Biochem Biophys Res Commun. 2017;488:648–54.
Basu R, Whitlock BM, Husson J, Le Floc'h A, Jin W, Oyler-Yaniv A, et al. Cytotoxic T cells use mechanical force to potentiate target cell killing. Cell. 2016;165:100–10.
Zhu J. T Helper Cell Differentiation, Heterogeneity, and Plasticity. Cold Spring Harb Perspect Biol. 2018;10(10):a030338. https://doi.org/10.1101/cshperspect.a030338.
Tamborero D, Rubio-Perez C, Muiños F, Sabarinathan R, Piulats JM, Muntasell A, et al. A pan-cancer landscape of interactions between solid tumors and infiltrating immune cell populations, Clinical cancer research : an official journal of the American Association for. Cancer Res. 2018;24:3717–28.
Bao S, Wu Q, McLendon RE, Hao Y, Shi Q, Hjelmeland AB, et al. Glioma stem cells promote radioresistance by preferential activation of the DNA damage response. Nature. 2006;444:756–60.
Cheung LS, Fu J, Kumar P, Kumar A, Urbanowski ME, Ihms EA, et al. Second-generation IL-2 receptor-targeted diphtheria fusion toxin exhibits antitumor activity and synergy with anti-PD-1 in melanoma. Proc Natl Acad Sci U S A. 2019;116:3100–5.
Nikolova T, Roos WP, Krämer OH, Strik HM, Kaina B. Chloroethylating nitrosoureas in cancer therapy: DNA damage, repair and cell death signaling, Biochimica et biophysica acta. Rev Cancer. 1868;2017:29–39.
Coussens LM, Werb Z. Inflammation and cancer. Nature. 2002;420:860–7.
Quail DF, Joyce JA. Microenvironmental regulation of tumor progression and metastasis. Nat Med. 2013;19:1423–37.
Hashimoto K, Nishimura S, Ito T, Oka N, Akagi M. Inflammatory Undifferentiated Pleomorphic Sarcoma Mimicking Bacteremia in an Elderly Patient: A Case Report. Medicina (Kaunas). 2021;57(2):175. https://doi.org/10.3390/medicina57020175.
von Schuckmann LA, Smith D, Hughes MCB, Malt M, van der Pols JC, Khosrotehrani K, et al. Associations of statins and diabetes with diagnosis of ulcerated cutaneous melanoma. J Invest Dermatol. 2017;137:2599–605.
Panigrahy D, Gartung A, Yang J, Yang H, Gilligan MM, Sulciner ML, et al. Preoperative stimulation of resolution and inflammation blockade eradicates micrometastases. J Clin Invest. 2019;129:2964–79.
Ethics approval and consent to participate
Following permission from the Chinese PLA General Hospital’s Human Research Ethics Committee, we collected eight pairs of CM and para-cancerous normal tissue specimens. Written informed consent was obtained from individual or guardian participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Xing, J., Li, Y., Chen, Y. et al. A novel inflammatory response-related signature predicts the prognosis of cutaneous melanoma and the effect of antitumor drugs. World J Surg Onc 20, 263 (2022). https://doi.org/10.1186/s12957-022-02726-8
- Cutaneous melanoma
- Inflammatory genes
- Prognostic model
- Immune states
- Tumor environment
- Drug sensitivity