Determination of prognosis in metastatic melanoma through integration of clinico-pathologic, mutation, mRNA, microRNA, and protein information

The MITF/mir-579-3p regulatory axis dictates BRAF-mutated melanoma cell fate in response to MAPK inhibitors

Therapy of melanoma has improved dramatically over the last years thanks to the development of targeted therapies (MAPKi) and immunotherapies. However, drug resistance continues to limit the efficacy of these therapies. Our research group has provided robust evidence as to the involvement of a set of microRNAs in the development of resistance to target therapy in BRAF-mutated melanomas. Among them, a pivotal role is played by the oncosuppressor miR-579-3p. Here we show that miR-579-3p and the microphthalmia-associated transcription factor (MITF) influence reciprocally their expression through positive feedback regulatory loops. In particular we show that miR-579-3p is specifically deregulated in BRAF-mutant melanomas and that its expression levels mirror those of MITF. Luciferase and ChIP studies show that MITF is a positive regulator of miR-579-3p, which is located in the intron 11 of the human gene ZFR (Zink-finger recombinase) and is co-transcribed with its host gene. Moreover, miR-579-3p, by targeting BRAF, is able to stabilize MITF protein thus inducing its own transcription. From biological points of view, early exposure to MAPKi or, alternatively miR-579-3p transfection, induce block of proliferation and trigger senescence programs in BRAF-mutant melanoma cells. Finally, the long-term development of resistance to MAPKi is able to select cells characterized by the loss of both miR-579-3p and MITF and the same down-regulation is also present in patients relapsing after treatments. Altogether these findings suggest that miR-579-3p/MITF interplay potentially governs the balance between proliferation, senescence and resistance to therapies in BRAF-mutant melanomas.

Construction and optimization of multi-platform precision pathways for precision medicine

In the enduring challenge against disease, advancements in medical technology have empowered clinicians with novel diagnostic platforms. Whilst in some cases, a single test may provide a confident diagnosis, often additional tests are required. However, to strike a balance between diagnostic accuracy and cost-effectiveness, one must rigorously construct the clinical pathways. Here, we developed a framework to build multi-platform precision pathways in an automated, unbiased way, recommending the key steps a clinician would take to reach a diagnosis. We achieve this by developing a confidence score, used to simulate a clinical scenario, where at each stage, either a confident diagnosis is made, or another test is performed. Our framework provides a range of tools to interpret, visualize and compare the pathways, improving communication and enabling their evaluation on accuracy and cost, specific to different contexts. This framework will guide the development of novel diagnostic pathways for different diseases, accelerating the implementation of precision medicine into clinical practice.

Comprehensive analysis of necroptotic patterns and associated immune landscapes in individualized treatment of skin cutaneous melanoma

Skin cutaneous melanoma (SKCM) constitutes a malignant cutaneous neoplasm characterized by an exceedingly unfavorable prognosis. Over the past years, necroptosis, a manifestation of inflammatory programmed cell demise, has gained substantial traction in its application. However, a conclusive correlation between the expression of necroptosis-related genes (NRGs) and SKCM patient's prognosis remains elusive. In this endeavor, we have undertaken an integrative analysis of genomic data, aiming to provide an exhaustive evaluation of the intricate interplay between melanoma necroptosis and immune-infiltration nuances within the tumor microenvironment. Through meticulous scrutiny, we have endeavored to discern the prognostic potency harbored by individual necroptosis-associated genes. Our efforts culminated in the establishment of a risk stratification framework, allowing for the appraisal of necroptosis irregularities within each afflicted cutaneous melanoma patient. Notably, those SKCM patients classified within the low-risk cohort exhibited a markedly elevated survival quotient, in stark contrast to their high-risk counterparts (p < 0.001). Remarkably, the low-risk cohort not only displayed a more favorable survival rate but also exhibited an enhanced responsiveness to immunotherapeutic interventions, relative to their high-risk counterparts. The outcomes of this investigation proffer insights into a conceivable mechanistic underpinning linking necroptosis-related attributes to the intricacies of the tumor microenvironment. This prompts a conjecture regarding the plausible association between necroptosis characteristics and the broader tumor microenvironmental milieu. However, it is imperative to emphasize that the pursuit of discerning whether the expression profiles of NRG genes can indeed be regarded as viable therapeutic targets necessitates further comprehensive exploration and scrutiny. In conclusion, our study sheds light on the intricate interrelationship between necroptosis-related factors and the tumor microenvironment, potentially opening avenues for therapeutic interventions. However, the prospect of translating these findings into clinical applications mandates rigorous investigation.

Machine learning-derived identification of tumor-infiltrating immune cell-related signature for improving prognosis and immunotherapy responses in patients with skin cutaneous melanoma

BACKGROUND

Immunoblockade therapy based on the PD-1 checkpoint has greatly improved the survival rate of patients with skin cutaneous melanoma (SKCM). However, existing anti-PD-1 therapeutic efficacy prediction markers often exhibit a poor situation of poor reliability in identifying potential beneficiary patients in clinical applications, and an ideal biomarker for precision medicine is urgently needed.

METHODS

10 multicenter cohorts including 4 SKCM cohorts and 6 immunotherapy cohorts were selected. Through the analysis of WGCNA, survival analysis, consensus clustering, we screened 36 prognostic genes. Then, ten machine learning algorithms were used to construct a machine learning-derived immune signature (MLDIS). Finally, the independent data sets (GSE22153, GSE54467, GSE59455, and in-house cohort) were used as the verification set, and the ROC index standard was used to evaluate the model.

RESULTS

Based on computing framework, we found that patients with high MLDIS had poor overall survival and has good prediction performance in all cohorts and in-house cohort. It is worth noting that MLDIS performs better in each data set than almost all models which from 51 prognostic signatures for SKCM. Meanwhile, high MLDIS have a positive prognostic impact on patients treated with anti-PD-1 immunotherapy by driving changes in the level of infiltration of immune cells in the tumor microenvironment. Additionally, patients suffering from SKCM with high MLDIS were more sensitive to immunotherapy.

CONCLUSIONS

Our study identified that MLDIS could provide new insights into the prognosis of SKCM and predict the immunotherapy response in patients with SKCM.

TimiGP: Inferring cell-cell interactions and prognostic associations in the tumor immune microenvironment through gene pairs

Determining the prognostic association of different immune cell types in the tumor microenvironment is critical for understanding cancer biology and developing new therapeutic strategies. However, this is challenging in certain cancer types, where the abundance of different immune subsets is highly correlated. In this study, we develop a computational method named TimiGP to overcome this challenge. Based on bulk gene expression and survival data, TimiGP infers cell-cell interactions that reveal the association between immune cell relative abundance and prognosis. As demonstrated in metastatic melanoma, TimiGP prioritizes immune cells critical in prognosis based on the identified cell-cell interactions. Highly consistent results are obtained by TimiGP when applied to seven independent melanoma datasets and when different cell-type marker sets are used as inputs. Additionally, TimiGP can leverage single-cell RNA sequencing data to delineate the tumor immune microenvironment at high resolutions across a wide range of cancer types.

Multi-omics profiling approach in food allergy

The prevalence of food allergy (FA) among children is increasing, affecting nearly 8% of children, and FA is the most common cause of anaphylaxis and anaphylaxis-related emergency department visits in children. Importantly, FA is a complex, multi-system, multifactorial disease mediated by food-specific immunoglobulin E (IgE) and type 2 immune responses and involving environmental and genetic factors and gene-environment interactions. Early exposure to external and internal environmental factors largely influences the development of immune responses to allergens. Genetic factors and gene-environment interactions have established roles in the FA pathophysiology. To improve diagnosis and identification of FA therapeutic targets, high-throughput omics approaches have emerged and been applied over the past decades to screen for potential FA biomarkers, such as genes, transcripts, proteins, and metabolites. In this article, we provide an overview of the current status of FA omics studies, namely genomic, transcriptomic, epigenomic, proteomic, exposomic, and metabolomic. The current development of multi-omics integration of FA studies is also briefly discussed. As individual omics technologies only provide limited information on the multi-system biological processes of FA, integration of population-based multi-omics data and clinical data may lead to robust biomarker discovery that could translate into advances in disease management and clinical care and ultimately lead to precision medicine approaches.

Identification of Pyroptosis-Relevant Signature in Tumor Immune Microenvironment and Prognosis in Skin Cutaneous Melanoma Using Network Analysis

Background

Pyroptosis is closely related to the programmed death of cancer cells as well as the tumor immune microenvironment (TIME) via the host-tumor crosstalk. However, the role of pyroptosis-related genes as prognosis and TIME-related biomarkers in skin cutaneous melanoma (SKCM) patients remains unknown.

Methods

We evaluated the expression profiles, copy number variations, and somatic mutations (CNVs) of 27 genes obtained from MSigDB database regulating pyroptosis among TCGA-SKCM patients. Thereafter, we conducted single-sample gene set enrichment analysis (ssGSEA) for evaluating pyroptosis-associated expression patterns among cases and for exploring the associations with clinicopathological factors and prognostic outcome. In addition, a prognostic pyroptosis-related signature (PPRS) model was constructed by performing Cox regression, weighted gene coexpression network analysis (WGCNA), and least absolute shrinkage and selection operator (LASSO) analysis to score SKCM patients. On the other hand, we plotted the ROC and survival curves for model evaluation and verified the robustness of the model through external test sets (GSE22153, GSE54467, and GSE65904). Meanwhile, we examined the relations of clinical characteristics, oncogene mutations, biological processes (BPs), tumor stemness, immune infiltration degrees, immune checkpoints (ICs), and treatment response with PPRS via multiple methods, including immunophenoscore (IPS) analysis, gene set variation analysis (GSVA), ESTIMATE, and CIBERSORT. Finally, we constructed a nomogram incorporating PPRS and clinical characteristics to improve risk evaluation of SKCM.

Results

Many pyroptosis-regulated genes showed abnormal expression within SKCM. TP53, TP63, IL1B, IL18, IRF2, CASP5, CHMP4C, CHMP7, CASP1, and GSDME were detected with somatic mutations, among which, a majority displayed CNVs at high frequencies. Pyroptosis-associated profiles established based on pyroptosis-regulated genes showed markedly negative relation to low stage and superior prognostic outcome. Blue module was found to be highly positively correlated with pyroptosis. Later, this study established PPRS based on the expression of 8 PAGs (namely, GBP2, HPDL, FCGR2A, IFITM1, HAPLN3, CCL8, TRIM34, and GRIPAP1), which was highly associated with OS, oncogene mutations, tumor stemness, immune infiltration degrees, IC levels, treatment responses, and multiple biological processes (including cell cycle and immunoinflammatory response) in training and test set samples.

Conclusions

Based on our observations, analyzing modification patterns associated with pyroptosis among diverse cancer samples via PPRS is important, which can provide more insights into TIME infiltration features and facilitate immunotherapeutic development as well as prognosis prediction.

Immune-related risk score: An immune-cell-pair-based prognostic model for cutaneous melanoma

Background

Melanoma is among the most malignant immunologic tumor types and is associated with high mortality. However, a considerable number of melanoma patients cannot benefit from immunotherapy owing to individual differences. This study attempts to build a novel prediction model of melanoma that fully considers individual differences in the tumor microenvironment.

Methods

An immune-related risk score (IRRS) was constructed based on cutaneous melanoma data from The Cancer Genome Atlas (TCGA). Single-sample gene set enrichment analysis (ssGSEA) was used to calculate immune enrichment scores of 28 immune cell signatures. We performed pairwise comparisons to obtain scores for cell pairs based on the difference in the abundance of immune cells within each sample. The resulting cell pair scores, in the form of a matrix of relative values of immune cells, formed the core of the IRRS.

Results

The area under the curve (AUC) for the IRRS was over 0.700, and when the IRRS was combined with clinical information, the AUC reached 0.785, 0.817, and 0.801 for the 1-, 3-, and 5-year survival, respectively. Differentially expressed genes between the two groups were enriched in staphylococcal infection and estrogen metabolism pathway. The low IRRS group showed a better immunotherapeutic response and exhibited more neoantigens, richer T-cell receptor and B-cell receptor diversity, and higher tumor mutation burden.

Conclusion

The IRRS enables a good prediction of prognosis and immunotherapy effect, based on the difference in the relative abundance of different types of infiltrating immune cells, and could provide support for further research in melanoma.

Omics technologies in allergy and asthma research: An EAACI position paper

Allergic diseases and asthma are heterogenous chronic inflammatory conditions with several distinct complex endotypes. Both environmental and genetic factors can influence the development and progression of allergy. Complex pathogenetic pathways observed in allergic disorders present a challenge in patient management and successful targeted treatment strategies. The increasing availability of high-throughput omics technologies, such as genomics, epigenomics, transcriptomics, proteomics, and metabolomics allows studying biochemical systems and pathophysiological processes underlying allergic responses. Additionally, omics techniques present clinical applicability by functional identification and validation of biomarkers. Therefore, finding molecules or patterns characteristic for distinct immune-inflammatory endotypes, can subsequently influence its development, progression, and treatment. There is a great potential to further increase the effectiveness of single omics approaches by integrating them with other omics, and nonomics data. Systems biology aims to simultaneously and longitudinally understand multiple layers of a complex and multifactorial disease, such as allergy, or asthma by integrating several, separated data sets and generating a complete molecular profile of the condition. With the use of sophisticated biostatistics and machine learning techniques, these approaches provide in-depth insight into individual biological systems and will allow efficient and customized healthcare approaches, called precision medicine. In this EAACI Position Paper, the Task Force "Omics technologies in allergic research" broadly reviewed current advances and applicability of omics techniques in allergic diseases and asthma research, with a focus on methodology and data analysis, aiming to provide researchers (basic and clinical) with a desk reference in the field. The potential of omics strategies in understanding disease pathophysiology and key tools to reach unmet needs in allergy precision medicine, such as successful patients' stratification, accurate disease prognosis, and prediction of treatment efficacy and successful prevention measures are highlighted.

Photothermal Nano-Vaccine Promoting Antigen Presentation and Dendritic Cells Infiltration for Enhanced Immunotherapy of Melanoma via Transdermal Microneedles Delivery

Immunotherapy has demonstrated the potential to cure melanoma, while the current response rate is still unsatisfactory in clinics. Extensive evidence indicates the correlation between the efficacy and pre-existing T-cell in tumors, whereas the baseline T-cell infiltration is lacking in low-response melanoma patients. Herein, we demonstrated the critical contribution of dendritic cells (DCs) on melanoma survival and baseline T-cell level, as well as the efficacy of immunotherapy. Capitalized on this fact, we developed a photothermal nano-vaccine to simultaneously promote tumor antigens presentation and DCs infiltration for enhanced immunotherapy. The nano-vaccine was composed of polyserotonin (PST) core and tannic acid (TA)/Mn2+ coordination-based metal-organic-framework (MOF) shell for β-catenin silencing DNAzyme loading, which was further integrated into dissolving microneedles to allow noninvasive and transdermal administration at melanoma skin. The nano-vaccine could rapidly penetrate skin upon microneedles insertion and exert a synergistically amplified photothermal effect to induce immunogenic cell death (ICD). The MOF shell then dissociated and released Mn2+ as a cofactor to self-activate DNAzyme for β-catenin suppression, which in turn caused a persistent CCL4 excretion to promote the infiltration of DCs into the tumor. Meanwhile, the liberated PST core could effectively capture and facilitate tumor antigens presentation to DCs. As a result, potent antitumor efficacies were achieved for both primary and distal tumors without any extra treatment, indicating the great promise of such a nano-vaccine for on-demand personalized immunotherapy of melanoma.

Core immune cell infiltration signatures identify molecular subtypes and promote precise checkpoint immunotherapy in cutaneous melanoma

Yutao Wang, China Medical University, ChinaThe tumor microenvironment (TME) has been shown to impact the prognosis of tumors in patients including cutaneous melanoma (CM); however, not all components of TME are important. Given the aforementioned situation, the functional immune cell contents correlated with CM patient prognosis are needed to optimize present predictive models and reflect the overall situation of TME. We developed a novel risk score named core tumor-infiltrating immune cell score (cTICscore), which showed certain advantages over existing biomarkers or TME-related signatures in predicting the prognosis of CM patients. Furthermore, we explored a new gene signature named cTILscore−related module gene score (cTMGs), based on four identified TME-associated genes (GCH1, GZMA, PSMB8, and PLAAT4) showing a close correlation with the cTICscore, which was generated by weighted gene co-expression network analysis and least absolute shrinkage and selection operator analysis to facilitate clinical application. Patients with low cTMGs had significantly better overall survival (OS, P = 0.002,< 0.001, = 0.002, and = 0.03, respectively) in the training and validating CM datasets. In addition, the area under the curve values used to predict the immune response in four CM cohorts were 0.723, 0.723, 0.754, and 0.792, respectively, and that in one gastric cohort was 0.764. Therefore, the four-gene signature, based on cTICscore, might improve prognostic information, serving as a predictive tool for CM patients receiving immunotherapy.cutaneous melanoma, tumor microenvironment, prognosis, immunotherapy, cTICscore

Cross-Platform Omics Prediction procedure: a statistical machine learning framework for wider implementation of precision medicine

In this modern era of precision medicine, molecular signatures identified from advanced omics technologies hold great promise to better guide clinical decisions. However, current approaches are often location-specific due to the inherent differences between platforms and across multiple centres, thus limiting the transferability of molecular signatures. We present Cross-Platform Omics Prediction (CPOP), a penalised regression model that can use omics data to predict patient outcomes in a platform-independent manner and across time and experiments. CPOP improves on the traditional prediction framework of using gene-based features by selecting ratio-based features with similar estimated effect sizes. These components gave CPOP the ability to have a stable performance across datasets of similar biology, minimising the effect of technical noise often generated by omics platforms. We present a comprehensive evaluation using melanoma transcriptomics data to demonstrate its potential to be used as a critical part of a clinical screening framework for precision medicine. Additional assessment of generalisation was demonstrated with ovarian cancer and inflammatory bowel disease studies.

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.

Identification of Lactate-Related Gene Signature for Prediction of Progression and Immunotherapeutic Response in Skin Cutaneous Melanoma

Skin cutaneous melanoma (SKCM) is a skin cancer type characterized by a high degree of immune cell infiltration. The potential function of lactate, a main metabolic product in the tumor microenvironment (TME) of SKCM, remains unclear. In this study, we systemically analyzed the predictive value of lactate-related genes (LRGs) for prognosis and response to immune checkpoint inhibitors (ICIs) in SKCM patients included from The Cancer Genome Atlas (TCGA) database. Cluster 3, by consensus clustering for 61 LRGs, manifested a worse clinical outcome, attributed to the overexpression of malignancy marks. In addition, we created a prognostic prediction model for high- and low-risk patients and verified its performance in a validation cohort, GSE65904. Between TME and the risk model, we found a negative relation of the immunocyte infiltration levels with patients’ risk scores. The low-risk cases had higher ICI expression and could benefit better from ICIs relative to the high-risk cases. Thus, the lactate-related prognosis risk signature may comprehensively provide a basis for future investigations on immunotherapeutic treatment for SKCM.

Pan-Cancer Analysis Reveals the Signature of TMC Family of Genes as a Promising Biomarker for Prognosis and Immunotherapeutic Response

Transmembrane Channel-like (TMC) genes are critical in the carcinogenesis, proliferation, and cell cycle of human cancers. However, the multi-omics features of TMCs and their role in the prognosis and immunotherapeutic response of human cancer have not been explored. We discovered that TMCs 4-8 were commonly deregulated and correlated with patient survival in a variety of cancers. For example, TMC5 and TMC8 were correlated with the relapse and overall survival rates of breast cancer and skin melanoma, respectively. These results were validated by multiple independent cohorts. TMCs were regulated by DNA methylation and somatic alterations, such as TMC5 amplification in breast cancer (523/1062, 49.2%). Six algorithms concordantly uncovered the critical role of TMCs in the tumor microenvironment, potentially regulating immune cell toxicity and lymphocytes infiltration. Moreover, TMCs 4-8 were correlated with tumor mutation burden and expression of PD-1/PD-L1/CTLA4 in 33 cancers. Thus, we established an immunotherapy response prediction (IRP) score based on the signature of TMCs 4-8. Patients with higher IRP scores showed higher immunotherapeutic responses in five cohorts of skin melanoma (area under curve [AUC] = 0.90 in the training cohort, AUCs range from 0.70 to 0.83 in the validation cohorts). Together, our study highlights the great potential of TMCs as biomarkers for prognosis and immunotherapeutic response, which can pave the way for further investigation of the tumor-infiltrating mechanisms and therapeutic potentials of TMCs in cancer.

Identification of an m6A Regulators-Mediated Prognosis Signature For Survival Prediction and Its Relevance to Immune Infiltration in Melanoma

Despite robust evidence for the role of m6A in cancer development and progression, its association with immune infiltration and survival outcomes in melanoma remains obscure. Here, we aimed to develop an m6A-related risk signature to improve prognostic and immunotherapy responder prediction performance in the context of melanoma. We comprehensively analyzed the m6A cluster and immune infiltration phenotypes of public datasets. The TCGA (n = 457) and eleven independent melanoma cohorts (n = 758) were used as the training and validation datasets, respectively. We identified two m6A clusters (m6A-clusterA and m6A-clusterB) based on the expression pattern of m6A regulators via unsupervised consensus clustering. IGF2BP1 (7.49%), KIAA1429 (7.06%), and YTHDC1 (4.28%) were the three most frequently mutated genes. There was a correlation between driver genes mutation statuses and the expression of m6A regulators. A significant difference in tumor-associated immune infiltration between two m6A clusters was detected. Compared with m6A-clusterA, the m6A-clusterB was characterized by a lower immune score and immune cell infiltration but higher mRNA expression-based stemness index (mRNAsi). An m6A-related risk signature consisting of 12 genes was determined via Cox regression analysis and divided the patients into low- and high-risk groups (IL6ST, MBNL1, NXT2, EIF2A, CSGALNACT1, C11orf58, CD14, SPI1, NCCRP1, BOK, CD74, PAEP). A nomogram was developed for the prediction of the survival rate. Compared with the high-risk group, the low-risk group was characterized by high expression of immune checkpoints and immunophenoscore (IPS), activation of immune-related pathways, and more enriched in immune cell infiltrations. The low-risk group had a favorable prognosis and contained the potential beneficiaries of the immune checkpoint blockade therapy and verified by the IMvigor210 cohort (n = 298). The m6A-related signature we have determined in melanoma highlights the relationships between m6A regulators and immune cell infiltration. The established risk signature was identified as a promising clinical biomarker of melanoma.

Development and validation of an immune gene set-based prognostic signature in cutaneous melanoma

We aimed to fully understand the landscape of the skin cutaneous melanoma (SKCM) microenvironment and develop an immune prognostic signature that can predict the prognosis for SKCM patients. RNA sequencing data and clinical information were downloaded from the Cancer Genome Atlas and Gene Expression Omnibus databases. The immune-prognostic signature was constructed by LASSO Cox regression analysis. We calculated the relative abundance of 29 immune-related gene sets based on the mRNA expression profiles of 314 SKCM patients in the Cancer Genome Atlas training set. Hierarchical clustering was performed to classify SKCM patients into three clusters: immunity-high, -medium and -low. The values of our prognostic model in predicting disease progression, metastasis and immunotherapeutic responses were also validated. In conclusion, the prognostic model demonstrated a powerful ability to distinguish and predict SKCM patients' prognosis.

Pyroptosis-Related Gene Signatures Can Robustly Diagnose Skin Cutaneous Melanoma and Predict the Prognosis

Skin cutaneous melanoma (SKCM) is a chronically malignant tumor with a high mortality rate. Pyroptosis, a kind of pro-inflammatory programmed cell death, has been linked to cancer in recent studies. However, the value of pyroptosis in the diagnosis and prognosis of SKCM is not clear. In this study, it was discovered that 20 pyroptosis-related genes (PRGs) differed in expression between SKCM and normal tissues, which were related to diagnosis and prognosis. Firstly, based on these genes, nine machine-learning algorithms were shown to perform well in constructing diagnostic classifiers, including K-Nearest Neighbor (KNN), logistic regression, Support Vector Machine (SVM), Artificial Neural Network (ANN), decision tree, random forest, XGBoost, LightGBM, and CatBoost. Secondly, the least absolute shrinkage and selection operator (LASSO) Cox regression analysis was applied and the prognostic model was constructed based on 9 PRGs. Subgroups in low and high risks determined by the prognostic model were shown to have different survival. Thirdly, functional enrichment analyses were performed by applying the gene set enrichment analysis (GSEA), and results suggested that the risk was related to immune response. In conclusion, the expression signatures of pyroptosis-related genes are effective and robust in the diagnosis and prognosis of SKCM, which is related to immunity.

Clinical and Molecular Correlates of NLRC5 Expression in Patients With Melanoma

NLRC5 is an important regulator in antigen presentation and inflammation, and its dysregulation promotes tumor progression. In melanoma, the impact of NLRC5 expression on molecular phenotype, clinical characteristics, and tumor features is largely unknown. In the present study, public datasets from the Cancer Cell Line Encyclopedia (CCLE), Gene Expression Omnibus (GEO), The Cancer Genome Atlas (TCGA), and cBioPortal were used to address these issues. We identify that NLRC5 is expressed in both immune cells and melanoma cells in melanoma samples and its expression is regulated by SPI1 and DNA methylation. NLRC5 expression is closely associated with Breslow thickness, Clark level, recurrence, pathologic T stage, and ulceration status in melanoma. Truncating/splice mutations rather than missense mutations in NLRC5 could compromise the expression of downstream genes. Low expression of NLRC5 is associated with poor prognosis, low activity of immune-related signatures, low infiltrating level of immune cells, and low cytotoxic score in melanoma. Additionally, NLRC5 expression correlates with immunotherapy efficacy in melanoma. In summary, these findings suggest that NLRC5 acts as a tumor suppressor in melanoma via modulating the tumor immune microenvironment. Targeting the NLRC5 related pathway might improve efficacy of immunotherapy for melanoma patients.

Conserved pan-cancer microenvironment subtypes predict response to immunotherapy

The clinical use of molecular targeted therapy is rapidly evolving but has primarily focused on genomic alterations. Transcriptomic analysis offers an opportunity to dissect the complexity of tumors, including the tumor microenvironment (TME), a crucial mediator of cancer progression and therapeutic outcome. TME classification by transcriptomic analysis of >10,000 cancer patients identifies four distinct TME subtypes conserved across 20 different cancers. The TME subtypes correlate with patient response to immunotherapy in multiple cancers, with patients possessing immune-favorable TME subtypes benefiting the most from immunotherapy. Thus, the TME subtypes act as a generalized immunotherapy biomarker across many cancer types due to the inclusion of malignant and microenvironment components. A visual tool integrating transcriptomic and genomic data provides a global tumor portrait, describing the tumor framework, mutational load, immune composition, anti-tumor immunity, and immunosuppressive escape mechanisms. Integrative analyses plus visualization may aid in biomarker discovery and the personalization of therapeutic regimens.

Potential prognostic value of PD-L1 and NKG2A expression in Indonesian patients with skin nodular melanoma

OBJECTIVE

Biomarker mRNA levels have been suggested to be predictors of patient survival and therapy response in melanoma cases. This study aimed to investigate the correlations between the mRNA expression levels of PD-L1 and NKG2A in melanoma tissue with clinicopathologic characteristics and survival in Indonesian primary nodular melanoma patients.

RESULTS

Thirty-one tissue samples were obtained; two were excluded from survival analysis due to Breslow depth of less than 4 mm. The median survival of upregulated and normoregulated PD-L1-patients were 15.800 ± 2.345 and 28.945 ± 4.126 months, respectively. However, this difference was not significant statistically (p = 0.086). Upregulated and normoregulated NKG2A patients differed very little in median survival time (25.943 ± 7.415 vs 26.470 ± 3.854 months; p = 0.981). Expression of PD-L1 and NKG2A were strongly correlated (r_s: 0.787, p < 0.001). No clinicopathologic associations with PD-L1 and NKG2A mRNA levels were observed. These results suggest that PD-L1 may have potential as a prognostic factor. Although an unlikely prognostic factor, NKG2A may become an adjunct target for therapy. The strong correlation between PD-L1 and NKG2A suggests that anti-PD-1 and anti-NKG2A agents could be effective in patients with PD-L1 upregulation. The mRNA levels of these two genes may help direct choice of immunotherapy and predict patient outcomes.

AutoEncoder-Based Computational Framework for Tumor Microenvironment Decomposition and Biomarker Identification in Metastatic Melanoma

Melanoma is one of the most aggressive cancer types whose prognosis is determined by both the tumor cell-intrinsic and -extrinsic features as well as their interactions. In this study, we performed systematic and unbiased analysis using The Cancer Genome Atlas (TCGA) melanoma RNA-seq data and identified two gene signatures that captured the intrinsic and extrinsic features, respectively. Specifically, we selected genes that best reflected the expression signals from tumor cells and immune infiltrate cells. Then, we applied an AutoEncoder-based method to decompose the expression of these genes into a small number of representative nodes. Many of these nodes were found to be significantly associated with patient prognosis. From them, we selected two most prognostic nodes and defined a tumor-intrinsic (TI) signature and a tumor-extrinsic (TE) signature. Pathway analysis confirmed that the TE signature recapitulated cytotoxic immune cell related pathways while the TI signature reflected MYC pathway activity. We leveraged these two signatures to investigate six independent melanoma microarray datasets and found that they were able to predict the prognosis of patients under standard care. Furthermore, we showed that the TE signature was also positively associated with patients' response to immunotherapies, including tumor vaccine therapy and checkpoint blockade immunotherapy. This study developed a novel computational framework to capture the tumor-intrinsic and -extrinsic features and identified robust prognostic and predictive biomarkers in melanoma.

A novel prognostic biomarker LCP2 correlates with metastatic melanoma-infiltrating CD8+ T cells

Lymphocyte cytosolic protein 2 (LCP2) is one of the SLP-76 family of adapters, which are critical intermediates in signal cascades downstream of several receptors. LCP2 regulates immunoreceptor signaling (such as T-cell receptors) and is also required for integrin signaling in neutrophils and platelets. However, the role of LCP2 in the tumor microenvironment is still unknown. In this study, we found a significant increase of mRNA and protein expression of LCP2 in metastatic skin cutaneous melanoma compared to normal skin. The upregulation of LCP2 was associated with good overall survival of patients with metastatic skin cutaneous melanoma, who received pharmacotherapy and radiation. GSEA signaling pathways analysis showed that LCP2 was involved in multiple pathways of immune response and correlation analysis revealed LCP2 was positively correlated with molecules in TCR signaling and 11 immune checkpoints, while LCP2 negatively correlated with 2 immune checkpoints in the metastatic skin cutaneous melanoma. According to the different expressions of LCP2, high LCP2 expression was positively correlated with more tumor-infiltrating CD8+ T cells. Furthermore, Kaplan-Meier plot indicated that LCP2 acted as a prognostic biomarker for progression-free survival of patients with metastatic skin cutaneous melanoma receiving anti-PD1 immunotherapy. In conclusion, our results integrated both the expression and function of LCP2 in melanoma using multiple tools, shedding light on the potential role of LCP2 in melanoma, and suggesting LCP2 serves as a prognostic biomarker and therapeutic target in anti-tumor immunity.

Systematic characterization of mutations altering protein degradation in human cancers

The ubiquitin-proteasome system (UPS) is the primary route for selective protein degradation in human cells. The UPS is an attractive target for novel cancer therapies, but the precise UPS genes and substrates important for cancer growth are incompletely understood. Leveraging multi-omics data across more than 9,000 human tumors and 33 cancer types, we found that over 19% of all cancer driver genes affect UPS function. We implicate transcription factors as important substrates and show that c-Myc stability is modulated by CUL3. Moreover, we developed a deep learning model (deepDegron) to identify mutations that result in degron loss and experimentally validated the prediction that gain-of-function truncating mutations in GATA3 and PPM1D result in increased protein stability. Last, we identified UPS driver genes associated with prognosis and the tumor microenvironment. This study demonstrates the important role of UPS dysregulation in human cancer and underscores the potential therapeutic utility of targeting the UPS.

Integration of Proteomics and Other Omics Data

In recent biomedical studies, multidimensional profiling, which collects proteomics as well as other types of omics data on the same subjects, is getting increasingly popular. Proteomics, transcriptomics, genomics, epigenomics, and other types of data contain overlapping as well as independent information, which suggests the possibility of integrating multiple types of data to generate more reliable findings/models with better classification/prediction performance. In this chapter, a selective review is conducted on recent data integration techniques for both unsupervised and supervised analysis. The main objective is to provide the "big picture" of data integration that involves proteomics data and discuss the "intuition" beneath the recently developed approaches without invoking too many mathematical details. Potential pitfalls and possible directions for future developments are also discussed.

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

INTRODUCTION

Skin cutaneous melanoma (SKCM) is a common skin malignancy worldwide, and its metastasis and mortality rates are high. The molecular characteristics exhibited by tumor-immune interactions have drawn the attention from researchers. Therefore, increased knowledge and new strategies to identify effective immune-related biomarkers may improve the clinical management of SKCM by providing more accurate prognostic information.

PATIENTS AND METHODS

In this study, we established a prognostic immune-related gene pair (IRGP) signature for predicting the survival of SKCM patients. The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases provided gene expression profiles together with clinical information, and the samples were randomly divided into three groups including the training, testing, and validation datasets. The regression model of least absolute shrinkage and selection operator (LASSO) helped to identify a 13-IRGP signature with a significant relation to the survival of SKCM patients.

RESULTS

The training, TCGA, and independent sets have an average value of area under the curve of 0.79, 0.76, and 0.82, respectively. In addition, this 13-IRGP signature can noticeably divide SKCM patients into high-risk group and low-risk group with significantly different prognoses. Many biological activities such as gene family were enriched among the genes in our IRGP signature. While analyzing the risk signature and clinical characteristics, there was a large difference in the risk score between T stage and tumor stage grouping. Finally, we constructed a nomogram and forest plots of the risk score and clinical features.

CONCLUSION

In summary, we developed a robust 13-IRGP prognostic signature in SKCM, which can identify and provide new insights into immunological biomarkers.

Association of HERV-K and LINE-1 hypomethylation with reduced disease-free survival in melanoma patients

Aim: To evaluate CpG methylation of long interspersed nuclear elements 1 (LINE-1) and human endogenous retrovirus K (HERV-K) retroelements as potential prognostic biomarkers in cutaneous melanoma. Materials & methods: Methylation of HERV-K and LINE-1 retroelements was assessed in resected melanoma tissues from 82 patients ranging in age from 14 to 88 years. In addition, nevi from eight patients were included for comparison with nonmalignant melanocytic lesions. Results: Methylation levels were lower in melanomas than in nevi. HERV-K and LINE-1 methylation were decreased in melanoma patients with clinical parameters associated with adverse prognosis, while they were independent of age and gender. Hypomethylation of HERV-K (but not LINE-1) was an independent predictor of reduced disease-free survival. Conclusion: HERV-K hypomethylation can be a potential independent biomarker of melanoma recurrence.

BRN2 and MITF together impact AXL expression in melanoma

The inverse relationship between transcription factor MITF and receptor tyrosine kinase AXL has received much attention recently. It is thought that melanoma tumors showing AXL^high /MITF^low levels are resistant to therapy. We show here that a population of cells within melanoma tumors with extremely high expression of AXL are negative/low for both MITF and the transcription factor BRN2. Depletion of both transcription factors from cultured melanoma cell lines produced an increase in AXL expression greater than depletion of MITF alone. Further, re-expression of BRN2 led to decreased AXL expression, indicating a role for BRN2 in regulation of AXL levels unrelated to effects on MITF level. As AXL has been recognized as a marker of therapy resistance, these cells may represent a population of cells responsible for disease relapse and as potential targets for therapeutic treatment.

A novel immune checkpoint-related seven-gene signature for predicting prognosis and immunotherapy response in melanoma

BACKGROUND

New emergence of immunotherapy has significantly improved clinical outcome of melanoma patients with advanced and metastatic diseases. We aimed to develop a gene signature based on the expression of PD-1/PD-L1 signaling pathway genes to predict prognosis and immunotherapy response in melanoma patients.

METHODS

Melanoma samples from The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO) were used as the training set and external validation sets respectively. Prognostic genes for overall survival (OS) were identified by univariate Cox regression analysis. Then a multi-gene risk signature was established with the Least Absolute Shrinkage and Selector Operation (LASSO) regression and multivariate Cox regression. The predictive and prognostic value of gene signature was evaluated by Kaplan Meier curve, Time-dependent receiver operating characteristic (ROC) curve, and area under curve (AUC). Gene set enrichment analysis (GSEA) was performed to investigate the discrepantly enriched biological processes between low-risk and high-risk group of melanoma patients.

RESULTS

A seven-gene risk signature (BATF2, CTLA4, EGFR, HLA-DQB1, IKBKG, PIK3R2, PPP3CA) was constructed. The signature was an independent risk factor for OS (hazard ratio = 1.544, p < 0.001) and it could robustly predict OS in both training and validation sets. Besides, high risk scores indicated advanced clinical stage and no response to immunotherapy for melanoma patients. GSEA demonstrated that high risk score was intimately associated with immune response and immune regulation. In conclusion, the novel seven-gene signature could serve as a robust biomarker for prognosis and a potential indicator of immunotherapy response in melanoma.

Combining texture features of whole slide images improves prognostic prediction of recurrence-free survival for cutaneous melanoma patients

BACKGROUND

Accurate prediction of recurrence-free survival (RFS) is important for the prognosis of cutaneous melanoma patients. The image-based pathological examination remains as the gold standard for diagnosis. It is of clinical interest to account for computer-aided processing of pathology image when performing prognostic analysis.

METHODS

We enrolled in this study a total of 152 patients from TCGA-SKCM (The Cancer Genome Atlas Skin Cutaneous Melanoma project) with complete information in recurrence-related survival time, baseline variables (clinicopathologic variables, mutation status of BRAF and NRAS genes), gene expression data, and whole slide image (WSI) features. We preprocessed WSI to segment global or nucleus areas, and extracted 3 types of texture features from each region. We performed cross validation and used multiple evaluation metrics including C-index and time-dependent AUC to determine the best model of predicting recurrence events. We further performed differential gene expression analysis between the higher and lower-risk groups within AJCC pathologic tumor stage III patients to explore the underlying molecular mechanisms driving risk stratification.

RESULTS

The model combining baseline variables and WSI features had the best performance among models with any other types of data integration. The prognostic risk score generated by this model could provide a higher-resolution risk stratification within pathologically defined subgroups. We found the selected image features captured important immune-related variations, such as the aberration of expression in T cell activation and proliferation gene sets, and therefore contributed to the improved prediction.

CONCLUSIONS

Our study provided a prognostic model based on the combination of baseline variables and computer-processed WSI features. This model provided more accurate prediction than models based on other types of data combination in recurrence-free survival analysis.

TRIAL REGISTRATION

This study was based on public open data from TCGA and hence the study objects were retrospectively registered.

miRNAs in the Diagnosis and Prognosis of Skin Cancer

Skin cancer is, at present, the most common type of malignancy in the Caucasian population. Its incidence has increased rapidly in the last decade for both melanoma and non-melanoma skin cancer. Differential expression profiles of microRNAs (miRNAs) have been reported for a variety of different cancers, including skin cancers. Since miRNAs’ discovery as regulators of gene expression, their importance grew in the field of oncology. miRNAs can post-transcriptionally regulate gene expression, tumor initiation, development progression, and aggressiveness. Nowadays, these short regulatory RNAs are perceived as one of the epigenetic markers for the identification of new diagnostic and/or prognostic molecular markers. Moreover, as miRNAs can drive tumorigenesis, they might eventually represent new therapy targets. Some miRNAs are pleiotropic, such as miR-214, which was found deregulated in several other tumors besides skin cancers. Some others are specific for one or more skin cancer types, like miR-21 and miR-221 for cutaneous melanoma and cutaneous squamous carcinoma or miR-155 for melanoma and cutaneous lymphoma. The goal of this review was to summarize some of the main miRNA detection technologies that are used to evaluate miRNAs in tissues and body fluids. Furthermore, their quantification limits, conformity, and robustness are discussed. Aberrant miRNA expression is analyzed for cutaneous melanoma, cutaneous squamous cell carcinoma (CSCC), skin lymphomas, cutaneous lymphoma, and Merkel cell carcinoma (MCC). In this type of disease, miRNAs are described as potential biomarkers to diagnose early lesion and/or early metastatic disease. In the future, whether in tissue or circulating in body fluids, miRNAs will gain their place in skin cancer diagnosis, prognosis, and future therapeutic targets.

MicroRNA heterogeneity in melanoma progression

The altered expression of miRNAs has been linked with neocarcinogenesis or the development of human malignancies including melanoma. Of significance, multiple clinical studies have documented that distinct sets of microRNAs (miRNAs) could be utilized as prognostic biomarkers for cancer development or predict the outcomes of treatment responses. To that end, an in-depth validation of such differentially expressed miRNAs is necessary in diverse settings of cancer patients in order to devise novel approaches to control tumor growth and/or enhance the efficacy of clinically-relevant therapeutic options. Moreover, considering the heterogeneity and sophisticated regulation of miRNAs, the precise delineation of their cellular targets could also be explored to design personalized medicine. Given the significance of miRNAs in regulating several key cellular processes of tumor cells including cell cycle progression and apoptosis, we review the findings of such miRNAs implicated in melanoma tumorigenesis. Understanding the novel mechanistic insights of such miRNAs will be useful for developing diagnostic or prognostic biomarkers or devising future therapeutic intervention for malignant melanoma.

Challenges in the Integration of Omics and Non-Omics Data

Omics data integration is already a reality. However, few omics-based algorithms show enough predictive ability to be implemented into clinics or public health domains. Clinical/epidemiological data tend to explain most of the variation of health-related traits, and its joint modeling with omics data is crucial to increase the algorithm’s predictive ability. Only a small number of published studies performed a “real” integration of omics and non-omics (OnO) data, mainly to predict cancer outcomes. Challenges in OnO data integration regard the nature and heterogeneity of non-omics data, the possibility of integrating large-scale non-omics data with high-throughput omics data, the relationship between OnO data (i.e., ascertainment bias), the presence of interactions, the fairness of the models, and the presence of subphenotypes. These challenges demand the development and application of new analysis strategies to integrate OnO data. In this contribution we discuss different attempts of OnO data integration in clinical and epidemiological studies. Most of the reviewed papers considered only one type of omics data set, mainly RNA expression data. All selected papers incorporated non-omics data in a low-dimensionality fashion. The integrative strategies used in the identified papers adopted three modeling methods: Independent, conditional, and joint modeling. This review presents, discusses, and proposes integrative analytical strategies towards OnO data integration.

Characterisation and validation of Mel38; A multi-tissue microRNA signature of cutaneous melanoma

Background

Histopathologic examination of melanocytic neoplasms can be challenging and subjective, with no specific circulating or tissue-based biomarkers currently available. Recently, a circulating 38-microRNA profile of melanoma (Mel38) was described. In this study, Mel38 expression and its impact on downstream mRNA regulation in solid tissue is examined.

Methods

Mel38 was applied to archival, clinically-annotated, solid-tissue genomic datasets representing benign naevi, primary and metastatic melanoma. Statistical analysis of the signature in relation to disease status, patient outcome and molecular pathways was performed.

Results

Mel38 is able to stratify genomic data from solid tissue biopsies on the basis of disease status and differences in melanoma-specific survival. Experimentally-verified messenger-RNA targets of Mel38 also exhibit prognostic expression patterns and represent key molecular pathways and events in melanoma development and progression.

Conclusion

The Mel38 microRNA profile may have diagnostic and prognostic utility in solid tissue as well as being a robust circulating biomarker of melanoma.

The Prognostic Significance of Low-Frequency Somatic Mutations in Metastatic Cutaneous Melanoma

Background: Little is known about the prognostic significance of somatically mutated genes in metastatic melanoma (MM). We have employed a combined clinical and bioinformatics approach on tumor samples from cutaneous melanoma (SKCM) as part of The Cancer Genome Atlas project (TCGA) to identify mutated genes with potential clinical relevance. Methods: After limiting our DNA sequencing analysis to MM samples (n = 356) and to the CANCER CENSUS gene list, we filtered out mutations with low functional significance (snpEFF). We performed Cox analysis on 53 genes that were mutated in ≥3% of samples, and had ≥50% difference in incidence of mutations in deceased subjects versus alive subjects. Results: Four genes were potentially prognostic [RAC1, FGFR1, CARD11, CIITA; false discovery rate (FDR) < 0.2]. We identified 18 additional genes (e.g., SPEN, PDGFRB, GNAS, MAP2K1, EGFR, TSC2) that were less likely to have prognostic value (FDR < 0.4). Most somatic mutations in these 22 genes were infrequent (< 10%), associated with high somatic mutation burden, and were evenly distributed across all exons, except for RAC1 and MAP2K1. Mutations in only 9 of these 22 genes were also identified by RNA sequencing in >75% of the samples that exhibited corresponding DNA mutations. The low frequency, UV signature type and RNA expression of the 22 genes in MM samples were confirmed in a separate multi-institution validation cohort (n = 413). An underpowered analysis within a subset of this validation cohort with available patient follow-up (n = 224) showed that somatic mutations in SPEN and RAC1 reached borderline prognostic significance [log-rank favorable (p = 0.09) and adverse (p = 0.07), respectively]. Somatic mutations in SPEN, and to a lesser extent RAC1, were not associated with definite gene copy number or RNA expression alterations. High (>2+) nuclear plus cytoplasmic expression intensity for SPEN was associated with longer melanoma-specific overall survival (OS) compared to lower (≤ 2+) nuclear intensity (p = 0.048). We conclude that expressed somatic mutations in infrequently mutated genes beyond the well-characterized ones (e.g., BRAF, RAS, CDKN2A, PTEN, TP53), such as RAC1 and SPEN, may have prognostic significance in MM.

A Leukocyte Infiltration Score Defined by a Gene Signature Predicts Melanoma Patient Prognosis

Melanoma is the most aggressive type of skin cancer in the United States with an increasing incidence. Melanoma lesions often exhibit high immunogenicity, with infiltrating immune cells playing important roles in regression of tumors occurring spontaneously or caused by therapeutic treatment. Computational and experimental methods have been used to estimate the abundance of immune cells in tumors, but their applications are limited by the requirement of large gene sets or multiple antibodies. Although the prognostic role of immune cells has been appreciated, a systematic investigation of their association with clinical factors, genomic features, prognosis and treatment response in melanoma is still lacking. This study, identifies a 25-gene signature based on RNA-seq data from The Cancer Genome Atlas (TCGA)-Skin Cutaneous Melanoma (TCGA-SKCM) dataset. This signature was used to calculate sample-specific Leukocyte Infiltration Scores (LIS) in six independent melanoma microarray datasets and scores were found to vary substantially between different melanoma lesion sites and molecular subtypes. For metastatic melanoma, LIS was prognostic in all datasets with high LIS being associated with good survival. The current approach provided additional prognostic information over established clinical factors, including age, tumor stage, and gender. In addition, LIS was predictive of patient survival in stage III melanoma, and treatment efficacy of tumor-specific antigen vaccine. IMPLICATIONS: This study identifies a 25-gene signature that effectively estimates the level of immune cell infiltration in melanoma, which provides a robust biomarker for predicting patient prognosis.

Using machine learning tools for protein database biocuration assistance

Biocuration in the omics sciences has become paramount, as research in these fields rapidly evolves towards increasingly data-dependent models. As a result, the management of web-accessible publicly-available databases becomes a central task in biological knowledge dissemination. One relevant challenge for biocurators is the unambiguous identification of biological entities. In this study, we illustrate the adequacy of machine learning methods as biocuration assistance tools using a publicly available protein database as an example. This database contains information on G Protein-Coupled Receptors (GPCRs), which are part of eukaryotic cell membranes and relevant in cell communication as well as major drug targets in pharmacology. These receptors are characterized according to subtype labels. Previous analysis of this database provided evidence that some of the receptor sequences could be affected by a case of label noise, as they appeared to be too consistently misclassified by machine learning methods. Here, we extend our analysis to recent and quite substantially modified new versions of the database and reveal their now extremely accurate labeling using several machine learning models and different transformations of the unaligned sequences. These findings support the adequacy of our proposed method to identify problematic labeling cases as a tool for database biocuration.

Altered integrin expression patterns shown by microarray in human cutaneous melanoma

A large variety of molecular pathways in melanoma progression suggests that no individual molecular alteration is crucial in itself. Our aim was to define the molecular alterations underlying metastasis formation. Gene expression profiling was performed using microarray and qRT-PCR to define alterations between matched primary and metastatic melanoma cell lines. These data were integrated with publicly available unmatched tissue data. The invasiveness of cell lines was determined by Matrigel invasion assays and invasive clones from primary melanoma-derived cell lines were also selected. Two metastatic cell line models were created: the regional lymph node WM983A-WM983A-WM983B and the distant lung WM793B-WM793B-1205Lu metastatic models. The majority of metastasis genes were downregulated and enriched in adhesion and ITGA6-B4 pathways. Upregulation of immune pathways was characteristic of distant metastases, whereas increased Rap1 signaling was specific for regional (sub)cutaneous metastases. qRT-PCR analysis of selected integrins (A2, A3, A4, A9, B5, B8, A6, B1, and B3) highlighted the possible importance of ITGA3/4 and B8 in the metastatic process, distinguishing regional and distant metastases. We identified functionally relevant gene clusters that influenced metastasis formation. Our data provide further evidence that integrin expression patterns may be important in distant metastasis formation.

A multi-step classifier addressing cohort heterogeneity improves performance of prognostic biomarkers in three cancer types

Cancer research continues to highlight the extensive genetic diversity that exists both between and within tumors. This intrinsic heterogeneity poses one of the central challenges to predicting patient clinical outcome and the personalization of treatments. Despite progress in some individual tumor types, it is not yet possible to prospectively, accurately classify patients by expected survival. One hypothesis proposed to explain this is that the prognostic classifiers developed to date are insufficiently sensitive and specific; however it is also possible that patients are not equally easy to classify by any given biomarker. We demonstrate in a cohort of 45 AJCC stage III melanoma patients that clinico-pathologic biomarkers can identify those patients that are most likely to be misclassified by a molecular biomarker. The process of modelling the classifiability of patients was then replicated in a cohort of 49 stage II breast cancer patients and 53 stage III colon cancer patients. A multi-step procedure incorporating this information not only improved classification accuracy but also indicated the specific clinical attributes that had made classification problematic in each cohort. These findings show that, even when cohorts are of moderate size, including features that explain the patient-specific performance of a prognostic biomarker in a classification framework can improve the modelling and estimation of survival.

microRNA-7-5p inhibits melanoma cell proliferation and metastasis by suppressing RelA/NF-κB

microRNA-7-5p (miR-7-5p) is a tumor suppressor in multiple cancer types and inhibits growth and invasion by suppressing expression and activity of the epidermal growth factor receptor (EGFR) signaling pathway. While melanoma is not typically EGFR-driven, expression of miR-7-5p is reduced in metastatic tumors compared to primary melanoma. Here, we investigated the biological and clinical significance of miR-7-5p in melanoma. We found that augmenting miR-7-5p expression in vitro markedly reduced tumor cell viability, colony formation and induced cell cycle arrest. Furthermore, ectopic expression of miR-7-5p reduced migration and invasion of melanoma cells in vitro and reduced metastasis in vivo. We used cDNA microarray analysis to identify a subset of putative miR-7-5p target genes associated with melanoma and metastasis. Of these, we confirmed nuclear factor kappa B (NF-κB) subunit RelA, as a novel direct target of miR-7-5p in melanoma cells, such that miR-7-5p suppresses NF-κB activity to decrease expression of canonical NF-κB target genes, including IL-1β, IL-6 and IL-8. Importantly, the effects of miR-7-5p on melanoma cell growth, cell cycle, migration and invasion were recapitulated by RelA knockdown. Finally, analysis of gene array datasets from multiple melanoma patient cohorts revealed an association between elevated RelA expression and poor survival, further emphasizing the clinical significance of RelA and its downstream signaling effectors. Taken together, our data show that miR-7-5p is a potent inhibitor of melanoma growth and metastasis, in part through its inactivation of RelA/NF-κB signaling. Furthermore, miR-7-5p replacement therapy could have a role in the treatment of this disease.

sPAGM: inferring subpathway activity by integrating gene and miRNA expression-robust functional signature identification for melanoma prognoses

MicroRNAs (miRNAs) regulate biological pathways by inhibiting gene expression. However, most current analytical methods fail to consider miRNAs, when inferring functional or pathway activities. In this study, we developed a model called sPAGM to infer subpathway activities by integrating gene and miRNA expressions. In this model, we reconstructed subpathway graphs by embedding miRNA components, and characterized subpathway activity (sPA) scores by simultaneously considering the expression levels of miRNAs and genes. The results showed that the sPA scores could distinguish different samples across tumor types, as well as samples between tumor and normal conditions. Moreover, the sPAGM model displayed more specificities than the entire pathway-based analyses. This model was applied to melanoma tumors to perform a prognosis analysis, which identified a robust 55-subpathway signature. By using The Cancer Genome Atlas and independently verified data sets, the subpathway-based signature significantly predicted the patients’ prognoses, which were independent of clinical variables. In the prognostic performance comparison, the sPAGM model was superior to the gene-only and miRNA-only methods. Finally, we dissected the functional roles and interactions of components within the subpathway signature. Taken together, the sPAGM model provided a framework for inferring subpathway activities and identifying functional signatures for clinical applications.

Innovative array-based assay for omics pattern in melanoma

Cutaneous melanoma remains a major health issue and still an important challenge for research. Thus, omics complex evaluation can provide a more specific molecular classification for this heterogeneous disease. Complex omics analysis based on genomic and proteomic microarrays can identify disease markers that prognosticate disease evolution or can monitor therapies efficacy. Among the technologies that gained momentum in the last years, array-based comparative genomic hybridization offered the possibility to analyze chromosomal numerical aberrations within cutaneous melanomas providing important support for molecular classification of melanoma tumors. This technology can identify new chromosomal alterations and discover new deregulated melanoma genes that can be further used as therapy targets. Integrating genetic profiling with clinical and pathological parameters would lead to seminal improvements in diagnosis, prognosis, and therapy.

Integrated analysis of multidimensional omics data on cutaneous melanoma prognosis

Multiple types of genetic, epigenetic, and genomic changes have been implicated in cutaneous melanoma prognosis. Many of the existing studies are limited in analyzing a single type of omics measurement and cannot comprehensively describe the biological processes underlying prognosis. As a result, the obtained prognostic models may be less satisfactory, and the identified prognostic markers may be less informative. The recently collected TCGA (The Cancer Genome Atlas) data have a high quality and comprehensive omics measurements, making it possible to more comprehensively and more accurately model prognosis. In this study, we first describe the statistical approaches that can integrate multiple types of omics measurements with the assistance of variable selection and dimension reduction techniques. Data analysis suggests that, for cutaneous melanoma, integrating multiple types of measurements leads to prognostic models with an improved prediction performance. Informative individual markers and pathways are identified, which can provide valuable insights into melanoma prognosis.

Identification, Review, and Systematic Cross-Validation of microRNA Prognostic Signatures in Metastatic Melanoma

In metastatic melanoma, it is vital to identify and validate biomarkers of prognosis. Previous studies have systematically evaluated protein biomarkers or mRNA-based expression signatures. No such analyses have been applied to microRNA (miRNA)-based prognostic signatures. As a first step, we identified two prognostic miRNA signatures from publicly available data sets (Gene Expression Omnibus/The Cancer Genome Atlas) of global miRNA expression profiling information. A 12-miRNA signature predicted longer survival after surgery for resection of American Joint Committee on Cancer stage III disease (>4 years, no sign of relapse) and outperformed American Joint Committee on Cancer standard-of-care prognostic markers in leave-one-out cross-validation analysis (error rates 34% and 38%, respectively). A similar 15-miRNA biomarker derived from The Cancer Genome Atlas miRNA-seq data performed slightly worse (39%) than these current biomarkers. Both signatures were then assessed for replication in two independent data sets and subjected to systematic cross-validation together with the three other miRNA-based prognostic signatures proposed in the literature to date. Five miRNAs (miR-142-5p, miR-150-5p, miR-342-3p, miR-155-5p, and miR-146b-5p) were reproducibly associated with patient outcome and have the greatest potential for application in the clinic. Our extensive validation approach highlighted among multiple independent cohorts the translational potential and limitations of miRNA signatures, and pointed to future directions in the analysis of this emerging class of markers.

Differential distribution improves gene selection stability and has competitive classification performance for patient survival

A consistent difference in average expression level, often referred to as differential expression (DE), has long been used to identify genes useful for classification. However, recent cancer studies have shown that when transcription factors or epigenetic signals become deregulated, a change in expression variability (DV) of target genes is frequently observed. This suggests that assessing the importance of genes by either differential expression or variability alone potentially misses sets of important biomarkers that could lead to improved predictions and treatments. Here, we describe a new approach for assessing the importance of genes based on differential distribution (DD), which combines information from differential expression and differential variability into a unified metric. We show that feature ranking and selection stability based on DD can perform two to three times better than DE or DV alone, and that DD yields equivalent error rates to DE and DV. Finally, assessing genes via differential distribution produces a complementary set of selected genes to DE and DV, potentially opening up new categories of biomarkers.

Recurrent inactivating RASA2 mutations in melanoma

Analysis of 501 melanoma exomes identified RASA2, encoding a RasGAP, as a tumor-suppressor gene mutated in 5% of melanomas. Recurrent loss-of-function mutations in RASA2 were found to increase RAS activation, melanoma cell growth and migration. RASA2 expression was lost in ≥30% of human melanomas and was associated with reduced patient survival. These findings identify RASA2 inactivation as a melanoma driver and highlight the importance of RasGAPs in cancer.

Label noise in subtype discrimination of class C G protein-coupled receptors: A systematic approach to the analysis of classification errors

BACKGROUND

The characterization of proteins in families and subfamilies, at different levels, entails the definition and use of class labels. When the adscription of a protein to a family is uncertain, or even wrong, this becomes an instance of what has come to be known as a label noise problem. Label noise has a potentially negative effect on any quantitative analysis of proteins that depends on label information. This study investigates class C of G protein-coupled receptors, which are cell membrane proteins of relevance both to biology in general and pharmacology in particular. Their supervised classification into different known subtypes, based on primary sequence data, is hampered by label noise. The latter may stem from a combination of expert knowledge limitations and the lack of a clear correspondence between labels that mostly reflect GPCR functionality and the different representations of the protein primary sequences.

RESULTS

In this study, we describe a systematic approach, using Support Vector Machine classifiers, to the analysis of G protein-coupled receptor misclassifications. As a proof of concept, this approach is used to assist the discovery of labeling quality problems in a curated, publicly accessible database of this type of proteins. We also investigate the extent to which physico-chemical transformations of the protein sequences reflect G protein-coupled receptor subtype labeling. The candidate mislabeled cases detected with this approach are externally validated with phylogenetic trees and against further trusted sources such as the National Center for Biotechnology Information, Universal Protein Resource, European Bioinformatics Institute and Ensembl Genome Browser information repositories.

CONCLUSIONS

In quantitative classification problems, class labels are often by default assumed to be correct. Label noise, though, is bound to be a pervasive problem in bioinformatics, where labels may be obtained indirectly through complex, many-step similarity modelling processes. In the case of G protein-coupled receptors, methods capable of singling out and characterizing those sequences with consistent misclassification behaviour are required to minimize this problem. A systematic, Support Vector Machine-based method has been proposed in this study for such purpose. The proposed method enables a filtering approach to the label noise problem and might become a support tool for database curators in proteomics.