The Identification and Validation of a Robust Immune-Associated Gene Signature in Cutaneous Melanoma

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.

Explore the impact of hypoxia-related genes (HRGs) in Cutaneous melanoma

BACKGROUND

Cutaneous melanoma (CM) has an overall poor prognosis due to a high rate of metastasis. This study aimed to explore the role of hypoxia-related genes (HRGs) in CM.

METHODS

We first used on-negative matrix factorization consensus clustering (NMF) to cluster CM samples and preliminarily analyzed the relationship of HRGs to CM prognosis and immune cell infiltration. Subsequently, we identified prognostic-related hub genes by univariate COX regression analysis and the least absolute shrinkage and selection operator (LASSO) and constructed a prognostic model. Finally, we calculated a risk score for patients with CM and investigated the relationship between the risk score and potential surrogate markers of response to immune checkpoint inhibitors (ICIs), such as TMB, IPS values, and TIDE scores.

RESULTS

Through NMF clustering, we identified high expression of HRGs as a risk factor for the prognosis of CM patients, and at the same time, increased expression of HRGs also indicated a poorer immune microenvironment. Subsequently, we identified eight gene signatures (FBP1, NDRG1, GPI, IER3, B4GALNT2, BGN, PKP1, and EDN2) by LASSO regression analysis and constructed a prognostic model.

CONCLUSION

Our study identifies the prognostic significance of hypoxia-related genes in melanoma and shows a novel eight-gene signature to predict the potential efficacy of ICIs.

Construction and Identification of an NLR-Associated Prognostic Signature Revealing the Heterogeneous Immune Response in Skin Cutaneous Melanoma

Background

Skin cutaneous melanoma (SKCM) is the deadliest dermatology tumor. Ongoing researches have confirmed that the NOD-like receptors (NLRs) family are crucial in driving carcinogenesis. However, the function of NLRs signaling pathway-related genes in SKCM remains unclear.

Objective

To establish and identify an NLRs-related prognostic signature and to explore its predictive power for heterogeneous immune response in SKCM patients.

Methods

Establishment of the predictive signature using the NLRs-related genes by least absolute shrinkage and selection operator-Cox regression analysis (LASSO-COX algorithm). Through univariate and multivariate COX analyses, NLRs signature's independent predictive effectiveness was proven. CIBERSORT examined the comparative infiltration ratios of 22 distinct types of immune cells. RT-qPCR and immunohistochemistry implemented expression validation for critical NLRs-related prognostic genes in clinical samples.

Results

The prognostic signature, including 7 genes, was obtained by the LASSO-Cox algorithm. In TCGA and validation cohorts, SKCM patients with higher risk scores had remarkably poorer overall survival. The independent predictive role of this signature was confirmed by multivariate Cox analysis. Additionally, a graphic nomogram demonstrated that the risk score of the NLRs signature has high predictive accuracy. SKCM patients in the low-risk group revealed a distinct immune microenvironment characterized by the significantly activated inflammatory response, interferon-α/γ response, and complement pathways. Indeed, several anti-tumor immune cell types were significantly accumulated in the low-risk group, including M1 macrophage, CD8 T cell, and activated NK cell. It is worth noting that our NLRs prognostic signature could serve as one of the promising biomarkers for predicting response rates to immune checkpoint blockade (ICB) therapy. Furthermore, the results of expression validation (RT-qPCR and IHC) were consistent with the previous analysis.

Conclusion

A promising NLRs signature with excellent predictive efficacy for SKCM was developed.

Machine learning-based construction of immunogenic cell death-related score for improving prognosis and response to immunotherapy in melanoma

BACKGROUND

Immunogenic cell death (ICD) is a form of regulated cell death (RCD) which could drive the activation of the innate and adaptive immune responses. In this work, we aimed to develop an ICD-related signature to facilitate the assessment of prognosis and immunotherapy response for melanoma patients.

METHODS

A set of machine learning methods, including consensus clustering, non-negative matrix factorization (NMF) method and least absolute shrinkage and selection operator (LASSO) logistic regression model, and bioinformatics analytic tools were integrated to construct an ICD-related risk score (ICDscore). CIBERSORT and ESTIMATE algorithm were used to evaluate the infiltration of immune cells. The 'pRRophetic' package in R and 6 cohorts of melanoma patients receiving immunotherapy were used for therapy sensitivity analyses. The predictive performance between ICDscore with other mRNA signatures were also compared.

RESULTS

The ICDscore could predict prognosis and immunotherapy response in multiple cohorts, and displayed superior performance than other forms of cell death-related signatures or 52 published signatures. The melanoma patients with low ICDscore were marked with high infiltration of immune cells, high expression of immune checkpoint inhibitor-related genes, and increased tumor mutation burden.

CONCLUSIONS

In conclusion, we constructed a stable and robust ICD-related signature for evaluating the prognosis and benefits of immunotherapy, and it could serve as a promising tool to guide decision-making and surveillance for individual melanoma patients.

Crosstalk of four kinds of cell deaths defines subtypes of cutaneous melanoma for precise immunotherapy and chemotherapy

Background

Cell death patterns can give therapeutic and biological clues that facilitate the development of individualized treatments for this lethal form of skin cancer.

Methods

We employed unsupervised clustering to establish robust classifications based on the four kinds of cell death-associated gene expression of 462 melanoma patients in the Cancer Genome Atlas (TCGA) and tested their reproducibility in two independent melanoma cohorts of 558 patients. We then used dimensionality reduction of graph learning to display the different characteristics of cell death patterns and immune microenvironments.

Results

We examined 570 cell death-associated gene expression data of melanoma patients for exploration, independent verification, and comprehensive classification of five reproducible melanoma subtypes (CS1 to CS5) with different genomic and clinical features. Patients in death-inactive subtypes (CS1, CS2, and CS5) had the least immune and stromal cell infiltration, and their prognosis was the poorest. A death-active subtype (CS4), on the other hand, had the highest infiltrated immune and stromal cells and elevated immune-checkpoints. As a result, these patients had the highest response to immunotherapy and the best prognosis. An additional subtype (CS3) had more diversified cell death and immune characteristics with moderate prognoses. Based on graph learning, we successfully divided the CS3 subtype into two subgroups (group A and group B) with distinct survival outcomes and immune features. Finally, we identified eight potential chemical drugs that were specifically targeted for the therapy of melanoma subtypes.

Conclusions

This research defines the intrinsic subtypes of melanoma based on the crosstalk of four kinds of cell deaths, which affords a blueprint for clinical strategies and guiding precise immunotherapy and chemotherapy for melanoma patients.

What genes are differentially expressed in individuals with schizophrenia? A systematic review

Schizophrenia is a severe, complex mental disorder characterized by a combination of positive symptoms, negative symptoms, and impaired cognitive function. Schizophrenia is highly heritable (~80%) with multifactorial etiology and complex polygenic genetic architecture. Despite the large number of genetic variants associated with schizophrenia, few causal variants have been established. Gaining insight into the mechanistic influences of these genetic variants may facilitate our ability to apply these findings to prevention and treatment. Though there have been more than 300 studies of gene expression in schizophrenia over the past 15 years, none of the studies have yielded consistent evidence for specific genes that contribute to schizophrenia risk. The aim of this work is to conduct a systematic review and synthesis of case-control studies of genome-wide gene expression in schizophrenia. Comprehensive literature searches were completed in PubMed, EmBase, and Web of Science, and after a systematic review of the studies, data were extracted from those that met the following inclusion criteria: human case-control studies comparing the genome-wide transcriptome of individuals diagnosed with schizophrenia to healthy controls published between January 1, 2000 and June 30, 2020 in the English language. Genes differentially expressed in cases were extracted from these studies, and overlapping genes were compared to previous research findings from the genome-wide association, structural variation, and tissue-expression studies. The transcriptome-wide analysis identified different genes than those previously reported in genome-wide association, exome sequencing, and structural variation studies of schizophrenia. Only one gene, GBP2, was replicated in five studies. Previous work has shown that this gene may play a role in immune function in the etiology of schizophrenia, which in turn could have implications for risk profiling, prevention, and treatment. This review highlights the methodological inconsistencies that impede valid meta-analyses and synthesis across studies. Standardization of the use of covariates, gene nomenclature, and methods for reporting results could enhance our understanding of the potential mechanisms through which genes exert their influence on the etiology of schizophrenia. Although these results are promising, collaborative efforts with harmonization of methodology will facilitate the identification of the role of genes underlying schizophrenia.

Prognostic modeling of patients with metastatic melanoma based on tumor immune microenvironment characteristics

Most of the malignant melanomas are already in the middle and advanced stages when they are diagnosed, which is often accompanied by the metastasis and spread of other organs. Besides, the prognosis of patients is bleak. The characteristics of the local immune microenvironment in metastatic melanoma have important implications for both tumor progression and tumor treatment. In this study, data on patients with metastatic melanoma from the TCGA and GEO datasets were selected for immune, stromal, and estimate scores, and overlapping differentially expressed genes were screened. A nine-IRGs prognostic model (ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22) was established by univariate COX regression, LASSO and multivariate COX regression. Receiver operating characteristic curves were used to test the predictive accuracy of the model. Immune infiltration was analyzed by using CIBERSORT and Xcell in high-risk and low-risk groups. The immune infiltration of the high-risk group was significantly lower than that of the low-risk group. Immune checkpoint analysis revealed that the expression of PDCD1, CTLA4, TIGIT, CD274, HAVR2 and LAG3 demonstrated the visible difference in groups with different levels of risk scores. WGCNA analysis found that the yellow-green module contained seven genes from the nine-IRG prognostic model, and the yellow-green module had the highest correlation with risk scores. The results of GO and KEGG suggested that the genes in the yellow-green module were mainly enriched in immune-related biological processes. Finally, the expression characteristics of ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22 were analyzed between metastatic melanoma and normal samples. Overall, a prognostic model for metastatic melanoma based on the tumor immune microenvironment characteristics was established, which left plenty of space for further studies. It could function well in helping people to understand characteristics of the immune microenvironment in metastatic melanoma.

Development and Validation of a CD8+ T Cell Infiltration-Related Signature for Melanoma Patients

Aim

Immunotherapy shows efficacy in only a subset of melanoma patients. Here, we intended to construct a risk score model to predict melanoma patients’ sensitivity to immunotherapy.

Methods

Integration analyses were performed on melanoma patients from high-dimensional public datasets. The CD8+ T cell infiltration related genes (TIRGs) were selected via TIMER and CIBERSORT algorithm. LASSO Cox regression was performed to screen for the crucial TIRGs. Single sample gene set enrichment analysis (ssGSEA) and ESTIMATE algorithm were used to evaluate the immune activity. The prognostic value of the risk score was determined by univariate and multivariate Cox regression analysis.

Results

184 candidate TIRGs were identified in melanoma patients. Based on the candidate TIRGs, melanoma patients were classified into three clusters which were characterized by different immune activity. Six signature genes were further screened out of 184 TIRGs and a representative risk score for patient survival was constructed based on these six signature genes. The risk score served as an indicator for the level of CD8+ T cell infiltration and acted as an independent prognostic factor for the survival of melanoma patients. By using the risk score, we achieved a good predicting result for the response of cancer patients to immunotherapy. Moreover, pan-cancer analysis revealed the risk score could be used in a wide range of non-hematologic tumors.

Conclusions

Our results showed the potential of using signature gene-based risk score as an indicator to predict melanoma patients’ sensitivity to immunotherapy.