Abstract:OBJECTIVES:Malignant melanoma is a highly malignant and heterogeneous skin cancer. Although immunotherapy has improved survival rates, the inhibitory effect of tumor microenvironment has weakened its efficacy. To improve survival and treatment strategies, we need to develop immune-related prognostic models. Based on the analysis of the Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and Sequence Read Archive (SRA) database, this study aims to establish an immune-related prognosis prediction model, and to evaluate the tumor immune microenvironment by risk score to guide immunotherapy.METHODS:Skin cutaneous melanoma (SKCM) transcriptome sequencing data and corresponding clinical information were obtained from the TCGA database, differentially expressed genes were analyzed, and prognostic models were developed using univariate Cox regression, the LASSO method, and stepwise regression. Differentially expressed genes in prognostic models confirmed by real-time reverse transcription PCR (real-time RT-PCR) and Western blotting. Survival analysis was performed by using the Kaplan-Meier method, and the effect of the model was evaluated by time-dependent receiver operating characteristic curve as well as multivariate Cox regression, and the prognostic model was validated by 2 GEO melanoma datasets. Furthermore, correlations between risk score and immune cell infiltration, Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) score, immune checkpoint mRNA expression levels, tumor immune cycle, or tumor immune micro-environmental pathways were analyzed. Finally, we performed association analysis for risk score and the efficacy of immunotherapy.RESULTS:We identified 4 genes that were differentially expressed in TCGA-SKCM datasets, which were mainly associated with the tumor immune microenvironment. A prognostic model was also established based on 4 genes. Among 4 genes, the mRNA and protein levels of killer cell lectin like receptor D1 (KLRD1), leukemia inhibitory factor (LIF), and cellular retinoic acid binding protein 2 (CRABP2) genes in melanoma tissues differed significantly from those in normal skin (all P<0.01). The prognostic model was a good predictor of prognosis for patients with SKCM. The patients with high-risk scores had significantly shorter overall survival than those with low-risk scores, and consistent results were achieved in the training cohort and multiple validation cohorts (P<0.001). The risk score was strongly associated with immune cell infiltration, ESTIMATE score, immune checkpoint mRNA expression levels, tumor immune cycle, and tumor immune microenvironmental pathways (P<0.001). The correlation analysis showed that patients with the high-risk scores were in an inhibitory immune microenvironment based on the prognostic model (P<0.01).CONCLUSIONS:The immune-related SKCM prognostic model constructed in this study can effectively predict the prognosis of SKCM patients. Considering its close correlation to the tumor immune microenvironment, the model has some reference value for clinical immunotherapy of SKCM.

Identification of Immune-Related Biomarkers Associated with Tumorigenesis and Prognosis in Skin Cutaneous Melanoma.

Characterization of M6a-Related Genes Landscape in Skin Cutaneous Melanoma to Aid Immunotherapy and Assess Prognosis

New Prognostic Biomarkers and Drug Targets for Skin Cutaneous Melanoma Via Comprehensive Bioinformatic Analysis and Validation

Identification of Genes Related to Immune Infiltration in the Tumor Microenvironment of Cutaneous Melanoma

Identification and Validation of Three Risk Models in Skin Melanoma Based on Bioinformatics Analysis

Identification of Potential Prognostic Biomarkers Associated With Cancerometastasis in Skin Cutaneous Melanoma

Identification of Therapeutic Targets and Prognostic Biomarkers Among Integrin Subunits in the Skin Cutaneous Melanoma Microenvironment

Analysis of immune subtypes based on immunogenomic profiling identifies prognostic signature for cutaneous melanoma

Analyzing Prognostic Hub Genes in the Microenvironment of Cutaneous Melanoma by Computer Integrated Bioinformatics

Analysis of prognostic immune and stromal- related genes in melanoma microenvironment

A Potential Diagnostic and Prognostic Biomarker TMEM176B and Its Relationship With Immune Infiltration in Skin Cutaneous Melanoma

The Identification of the Metabolism Subtypes of Skin Cutaneous Melanoma Associated with the Tumor Microenvironment and the Immunotherapy

Prognostic Role of Tumor Mutation Burden Combined With Immune Infiltrates in Skin Cutaneous Melanoma Based on Multi-Omics Analysis

Analysis on Tumor Immune Microenvironment and Construction of a Prognosis Model for Immune-Related Skin Cutaneous Melanoma.

Comprehensive analysis of prognostic immune‐related genes in the tumor microenvironment of cutaneous melanoma

Development and Validation of an Immune-Related Gene Pair Signature in Skin Cutaneous Melanoma

Characterization of coagulation-related gene signature to predict prognosis and tumor immune microenvironment in skin cutaneous melanoma

Multi-omics analysis of tumor mutation burden combined with immune infiltrates in melanoma

C1QA, C1QB, and GZMB are novel prognostic biomarkers of skin cutaneous melanoma relating tumor microenvironment

CXC Chemokines as Therapeutic Targets and Prognostic Biomarkers in Skin Cutaneous Melanoma Microenvironment

Based on Scrna-Seq and Bulk RNA-seq to Establish Tumor Immune Microenvironment-Associated Signature of Skin Melanoma and Predict Immunotherapy Response