Integrative genome comparison of primary and metastatic melanomas

CDKN2A, a key gene in copper-induced cell death model, influencing melanoma invasion and apoptosis

Skin cutaneous melanoma (SKCM) is one of the most lethal cancers translating into 75% of skin cancer-related deaths. Despite the advances in SKCM management and treatment strategies, the overall survival of patients remains unsatisfactory due to the metastatic properties of SKCM as well as the absence of effective prognostic biomarkers. Recent studies have shown that overload copper renders accumulation of mitochondrial proteins and fuels a form of cell death at odds with known death mechanisms and is hinged on mitochondrial respiration, the so-called cuproptosis. However, the exact role of cuproptosis in SKCM development and progression is unknown, and painting a clear picture of its functions in SKCM is fraught with challenges. A more systematic investigation is justified. In this study, we were posed to dissect the clout and latent regulatory mechanisms of cuproptosis-related genes (CRGs) in reining in SKCM progression. Also, we identified three CRGs that stood out were used to construct a prognostic model, which could be employed to predict the prognosis of patients with SKCM. Finally, through pan-cancer analysis, we found that the four cuproptosis key genes play a role in multiple tumors, suggesting that cuproptosis may impact tumor progression at the pan-cancer level. Taken together, these findings may not only contribute to the development of treatment strategies but also provide clues for treatment decision-making.

Exploring new mechanisms in cancer molecular pathways and pathogenic cell transformation: PIP4K2A as a prognostic marker and therapeutic target in cutaneous malignant melanoma

BACKGROUND

Cutaneous malignant melanoma is a very aggressive and metastatic form of skin cancer, typically linked with poor outcomes. Advances in genomic analysis have underscored the crucial role of T cells in tumor immunity. Immune checkpoint inhibitors have notably transformed melanoma treatment by boosting T cell activity. Studies of gene expression have found that the phosphatidylinositol-4-phosphate kinase 2A (PIP4K2A) gene is abnormally expressed in various tumors, indicating its potential role in tumor progression. Utilizing single-cell sequencing and machine learning, researchers can now explore the complex interactions between T cells and melanoma cells at a genomic level. This study aimed to investigate the role of the PIP4K2A gene in cutaneous malignant melanoma, with a focus on its influence on T cell-mediated immune responses.

METHODS

Samples from cutaneous melanoma patients were analysed by single-cell transcriptome for differentially expressed genes and signalling pathways associated with cutaneous melanoma. Then, genes were identified and predictive models were built based on the transcriptomic data using machine learning models to assess whether the expression level of PIP4K2A could effectively predict the malignancy and prognosis of cutaneous melanoma. In addition, we also performed drug therapy predictive analysis and immunotherapy analysis.Finally, the critical role of PIP4K2A in cutaneous melanoma was further confirmed by immunohistochemistry.

RESULTS

The PIP4K2A gene exhibited a significantly elevated expression level in cutaneous malignant melanoma, showing a strong correlation with the clinical stage and patient prognosis. At the therapeutic level, high PIP4K2A expression is less responsive to immunotherapy, and this gene is a risk factor for drug therapy in cutaneous malignant melanoma. Additionally, our experimental outcomes validated this observation.

CONCLUSIONS

The PIP4K2A gene could be a crucial prognostic marker for cutaneous malignant melanoma, as it significantly affects T cell activity within the tumor microenvironment. This study offers essential insights into melanoma pathogenesis and assists in pinpointing new early diagnostic markers and therapeutic targets. Utilizing advanced genomic tools and computational techniques, the research enhances our understanding of T cell dynamics in melanoma, facilitating the development of personalized medicine and more effective immunotherapy strategies.

Protein kinase A and local signaling in cancer

Protein kinase A (PKA) is a basophilic kinase implicated in the modulation of many cell-signaling and physiological processes. PKA also contributes to cancer-relevant events such as growth factor action, cell cycle control, cell migration and tumor metabolism. Germline and somatic mutations in PKA, gene amplifications, and chromosome rearrangements that encode kinase fusions, are linked to a growing number of malignant neoplasms. Mislocalization of PKA by exclusion from A-Kinase Anchoring Protein (AKAP) signaling islands further underlies cancer progression. This article highlights the influence of AKAP signaling and local kinase action in selected hallmarks of cancer. We also feature the utility of kinase inhibitor drugs as frontline and future anti-cancer therapies.

The classification of melanocytic gene signatures

Gene expression profiling technologies have revolutionized cell biology, enabling researchers to identify gene signatures linked to various biological attributes of melanomas, such as pigmentation status, differentiation state, proliferative versus invasive capacity, and disease progression. Although the discovery of gene signatures has significantly enhanced our understanding of melanocytic phenotypes, reconciling the numerous signatures reported across independent studies and different profiling platforms remains a challenge. Current methods for classifying melanocytic gene signatures depend on exact gene overlap and comparison with unstandardized baseline transcriptomes. In this study, we aimed to categorize published gene signatures into clusters based on their similar patterns of expression across clinical cutaneous melanoma specimens. We analyzed nearly 800 melanoma samples from six gene expression repositories and developed a classification framework for gene signatures that is resilient against biases in gene identification across profiling platforms and inconsistencies in baseline standards. Using 39 frequently cited published gene signatures, our analysis revealed seven principal classes of gene signatures that correlate with previously identified phenotypes: Differentiated, Mitotic/MYC, AXL, Amelanotic, Neuro, Hypometabolic, and Invasive. Each class is consistent with the phenotypes that the constituent gene signatures represent, and our classification method does not rely on overlapping genes between signatures. To facilitate broader application, we created WIMMS (what is my melanocytic signature, available at https://wimms.tanlab.org/), a user-friendly web application. WIMMS allows users to categorize any gene signature, determining its relationship to predominantly cited signatures and its representation within the seven principal classes.

A Predictive Network-Based Immune Checkpoint Blockade Immunotherapeutic Signature Optimizing Patient Selection and Treatment Strategies

Immune checkpoint blockade (ICB) therapy has brought significant advancements to the field of oncology. However, the diverse responses among patients highlight the need for more accurate predictive tools. In this study, insights are drawn from tumor-immunology pathways, and a novel network-based ICB immunotherapeutic signature, termed ICBnetIS, is constructed. The signature is derived from advanced biological network-based computational strategies involving co-expression networks and molecular interactions networks. The efficacy of ICBnetIS is established through its association with enhanced patient survival and a robust immune response characterized by diverse immune cell infiltration and active anti-tumor immune pathways. The validation process positions ICBnetIS as an effective tool in predicting responses to ICB therapy, analyzing ICB data from a broad collection of over 700 samples from multiple cancer types of more than 15 datasets. It achieves an aggregated prediction AUC of 0.784, which outperforms the other nine renowned immunotherapeutic signatures, indicating the superior predictive capability of ICBnetIS. To sum up, the findings suggest ICBnetIS as a potent tool in predicting ICB therapy responses, offering significant implications for patient selection and treatment optimization in oncology. The study highlights the role of ICBnetIS in advancing personalized treatment strategies, potentially transforming the clinical landscape of ICB therapy.

A comprehensive identification of potential molecular targets and small drugs candidate for melanoma cancer using bioinformatics and network-based screening approach

Melanoma is the third most common malignant skin tumor and has increased in morbidity and mortality over the previous decade due to its rapid spread into the bloodstream or lymphatic system. This study used integrated bioinformatics and network-based methodologies to reliably identify molecular targets and small molecular medicines that may be more successful for Melanoma diagnosis, prognosis and treatment. The statistical LIMMA approach utilized for bioinformatics analysis in this study found 246 common differentially expressed genes (cDEGs) between case and control samples from two microarray gene-expression datasets (GSE130244 and GSE15605). Protein-protein interaction network study revealed 15 cDEGs (PTK2, STAT1, PNO1, CXCR4, WASL, FN1, RUNX2, SOCS3, ITGA4, GNG2, CDK6, BRAF, AGO2, GTF2H1 and AR) to be critical in the development of melanoma (KGs). According to regulatory network analysis, the most important transcriptional and post-transcriptional regulators of DEGs and hub-DEGs are ten transcription factors and three miRNAs. We discovered the pathogenetic mechanisms of MC by studying DEGs' biological processes, molecular function, cellular components and KEGG pathways. We used molecular docking and dynamics modeling to select the four most expressed genes responsible for melanoma malignancy to identify therapeutic candidates. Then, utilizing the Connectivity Map (CMap) database, we analyzed the top 4-hub-DEGs-guided repurposable drugs. We validated four melanoma cancer drugs (Fisetin, Epicatechin Gallate, 1237586-97-8 and PF 431396) using molecular dynamics simulation with their target proteins. As a result, the results of this study may provide resources to researchers and medical professionals for the wet-lab validation of MC diagnosis, prognosis and treatments.Communicated by Ramaswamy H. Sarma.

Integrative and comprehensive pan-cancer analysis of ubiquitin specific peptidase 11 (USP11) as a prognostic and immunological biomarker

The significance of USP11 as a critical regulator in cancer has garnered substantial attention, primarily due to its catalytic activity as a deubiquitinating enzyme. Nonetheless, a thorough evaluation of USP11 across various cancer types in pan-cancer studies remains absent. Our analysis integrates data from a variety of sources, including five immunotherapy cohorts, thirty-three cohorts from The Cancer Genome Atlas (TCGA), and sixteen cohorts from the Gene Expression Omnibus (GEO), two of which involve single-cell transcriptomic data. Our findings indicate that aberrant USP11 expression is predictive of survival outcomes across various cancer types. The highest frequency of genomic alterations was observed in uterine corpus endometrial carcinoma (UCEC), with single-cell transcriptome analysis revealing significantly higher USP11 expression in plasmacytoid dendritic cells and mast cells. Notably, USP11 expression was associated with the infiltration levels of CD8+ T cells and natural killer (NK) activated cells. Additionally, in the skin cutaneous melanoma (SKCM) phs000452 cohort, patients with higher USP11 mRNA levels during immunotherapy experienced a significantly shorter median progression-free survival. USP11 emerges as a promising molecular biomarker with significant potential for predicting patient prognosis and immunoreactivity across various cancer types.

Proteogenomic insights into the biology and treatment of pan-melanoma

Melanoma is one of the most prevalent skin cancers, with high metastatic rates and poor prognosis. Understanding its molecular pathogenesis is crucial for improving its diagnosis and treatment. Integrated analysis of multi-omics data from 207 treatment-naïve melanomas (primary-cutaneous-melanomas (CM, n = 28), primary-acral-melanomas (AM, n = 81), primary-mucosal-melanomas (MM, n = 28), metastatic-melanomas (n = 27), and nevi (n = 43)) provides insights into melanoma biology. Multivariate analysis reveals that PRKDC amplification is a prognostic molecule for melanomas. Further proteogenomic analysis combined with functional experiments reveals that the cis-effect of PRKDC amplification may lead to tumor proliferation through the activation of DNA repair and folate metabolism pathways. Proteome-based stratification of primary melanomas defines three prognosis-related subtypes, namely, the ECM subtype, angiogenesis subtype (with a high metastasis rate), and cell proliferation subtype, which provides an essential framework for the utilization of specific targeted therapies for particular melanoma subtypes. The immune classification identifies three immune subtypes. Further analysis combined with an independent anti-PD-1 treatment cohort reveals that upregulation of the MAPK7-NFKB signaling pathway may facilitate T-cell recruitment and increase the sensitivity of patients to immunotherapy. In contrast, PRKDC may reduce the sensitivity of melanoma patients to immunotherapy by promoting DNA repair in melanoma cells. These results emphasize the clinical value of multi-omics data and have the potential to improve the understanding of melanoma treatment.

IGSF3 is a homophilic cell adhesion molecule that drives lung metastasis of melanoma by promoting adhesion to vascular endothelium

The immunoglobulin superfamily (IgSF) is one of the largest families of cell-surface molecules involved in various cell-cell interactions, including cancer-stromal interactions. In this study, we undertook a comprehensive RT-PCR-based screening for IgSF molecules that promote experimental lung metastasis in mice. By comparing the expression of 325 genes encoding cell-surface IgSF molecules between mouse melanoma B16 cells and its highly metastatic subline, B16F10 cells, we found that expression of the immunoglobulin superfamily member 3 gene (Igsf3) was significantly enhanced in B16F10 cells than in B16 cells. Knockdown of Igsf3 in B16F10 cells significantly reduced lung metastasis following intravenous injection into C57BL/6 mice. IGSF3 promoted adhesion of B16F10 cells to vascular endothelial cells and functioned as a homophilic cell adhesion molecule between B16F10 cells and vascular endothelial cells. Notably, the knockdown of IGSF3 in either B16F10 cells or vascular endothelial cells suppressed the transendothelial migration of B16F10 cells. Moreover, IGSF3 knockdown suppressed the extravasation of B16F10 cells into the lungs after intravenous injection. These results suggest that IGSF3 promotes the metastatic potential of B16F10 cells in the lungs by facilitating their adhesion to vascular endothelial cells.

Acute myeloid leukemia with a novel AKAP9::PDGFRA fusion transformed from essential thrombocythemia: A case report and mini review

Acute myeloid leukemia (AML) is a heterogeneous hematological malignancy associated with various combinations of gene mutations, epigenetic abnormalities, and chromosome rearrangement-related gene fusions. Despite the significant degree of heterogeneity in its pathogenesis, many gene fusions and point mutations are recurrent in AML and have been employed in risk stratification over the last several decades. Gene fusions have long been recognized for understanding tumorigenesis and their proven roles in clinical diagnosis and targeted therapies. Advances in DNA sequencing technologies and computational biology have contributed significantly to the detection of known fusion genes as well as for the discovery of novel ones. Several recurring gene fusions in AML have been linked to prognosis, treatment response, and disease progression. In this report, we present a case with a long history of essential thrombocythemia and hallmark CALR mutation transforming to AML characterized by a previously unreported AKAP9::PDGFRA fusion gene. We propose mechanisms by which this fusion may contribute to the pathogenesis of AML and its potential as a molecular target for tyrosine kinase inhibitors.

Integrated in vivo functional screens and multi-omics analyses identify α-2,3-sialylation as essential for melanoma maintenance

Glycosylation is a hallmark of cancer biology, and altered glycosylation influences multiple facets of melanoma growth and progression. To identify glycosyltransferases, glycans, and glycoproteins essential for melanoma maintenance, we conducted an in vivo growth screen with a pooled shRNA library of glycosyltransferases, lectin microarray profiling of benign nevi and melanoma patient samples, and mass spectrometry-based glycoproteomics. We found that α-2,3 sialyltransferases ST3GAL1 and ST3GAL2 and corresponding α-2,3-linked sialosides are upregulated in melanoma compared to nevi and are essential for melanoma growth in vivo and in vitro. Glycoproteomics revealed that glycoprotein targets of ST3GAL1 and ST3GAL2 are enriched in transmembrane proteins involved in growth signaling, including the amino acid transporter Solute Carrier Family 3 Member 2 (SLC3A2/CD98hc). CD98hc suppression mimicked the effect of ST3GAL1 and ST3GAL2 silencing, inhibiting melanoma cell proliferation. We found that both CD98hc protein stability and its pro-survival effect in melanoma are dependent upon α-2,3 sialylation mediated by ST3GAL1 and ST3GAL2. In summary, our studies reveal that α-2,3-sialosides functionally contribute to melanoma maintenance, supporting ST3GAL1 and ST3GAL2 as novel therapeutic targets in these tumors.

A novel method to assess copy number variations in melanocytic neoplasms: Droplet digital PCR for precise quantitation of MYC and MYB genes

INTRODUCTION

While most melanocytic neoplasms can be classified as either benign or malignant by histopathology alone, ancillary molecular diagnostic tests can be necessary to establish the correct diagnosis in challenging cases. Currently, the detection of copy number variations (CNVs) by fluorescence in situ hybridization and chromosomal microarray (CMA) are the most popular methods, but remain expensive and inaccessible. We aim to develop a relatively inexpensive, fast, and accessible molecular assay to detect CNVs relevant to melanoma using droplet digital polymerase chain reaction (ddPCR) technology.

METHODS

In this proof-of-concept study, we evaluated CNVs in MYC and MYB genes from 73 cases of benign nevi, borderline melanocytic lesions, and primary and metastatic melanoma at our institution from 2015 to 2022. A multiplexed ddPCR assay and CMA were performed on each sample, and the results were compared.

RESULTS

Concordance analysis of ddPCR with CMA for quantification of MYC and MYB CNVs revealed a sensitivity and specificity of 89% and 86% for MYC and 83% and 74% for MYB, respectively.

CONCLUSION

We demonstrate the first use of a multiplexed ddPCR assay to identify CNVs in melanocytic neoplasms. With further improvement and validation, ddPCR may represent a low-cost and rapid tool to aid in the diagnosis of histopathologically ambiguous melanocytic tumors.

Super-enhancer-driven expression of BAHCC1 promotes melanoma cell proliferation and genome stability

Super-enhancers (SEs) are stretches of enhancers ensuring a high level of expression of key genes associated with cell function. The identification of cancer-specific SE-driven genes is a powerful means for the development of innovative therapeutic strategies. Here, we identify a MITF/SOX10/TFIIH-dependent SE promoting the expression of BAHCC1 in a broad panel of melanoma cells. BAHCC1 is highly expressed in metastatic melanoma and is required for tumor engraftment, growth, and dissemination. Integrative genomics analyses reveal that BAHCC1 is a transcriptional regulator controlling expression of E2F/KLF-dependent cell-cycle and DNA-repair genes. BAHCC1 associates with BRG1-containing remodeling complexes at the promoters of these genes. BAHCC1 silencing leads to decreased cell proliferation and delayed DNA repair. Consequently, BAHCC1 deficiency cooperates with PARP inhibition to induce melanoma cell death. Our study identifies BAHCC1 as an SE-driven gene expressed in melanoma and demonstrates how its inhibition can be exploited as a therapeutic target.

Logic-based modeling and drug repurposing for the prediction of novel therapeutic targets and combination regimens against E2F1-driven melanoma progression

Melanoma presents increasing prevalence and poor outcomes. Progression to aggressive stages is characterized by overexpression of the transcription factor E2F1 and activation of downstream prometastatic gene regulatory networks (GRNs). Appropriate therapeutic manipulation of the E2F1-governed GRNs holds the potential to prevent metastasis however, these networks entail complex feedback and feedforward regulatory motifs among various regulatory layers, which make it difficult to identify druggable components. To this end, computational approaches such as mathematical modeling and virtual screening are important tools to unveil the dynamics of these signaling networks and identify critical components that could be further explored as therapeutic targets. Herein, we integrated a well-established E2F1-mediated epithelial-mesenchymal transition (EMT) map with transcriptomics data from E2F1-expressing melanoma cells to reconstruct a core regulatory network underlying aggressive melanoma. Using logic-based in silico perturbation experiments of a core regulatory network, we identified that simultaneous perturbation of Protein kinase B (AKT1) and oncoprotein murine double minute 2 (MDM2) drastically reduces EMT in melanoma. Using the structures of the two protein signatures, virtual screening strategies were performed with the FDA-approved drug library. Furthermore, by combining drug repurposing and computer-aided drug design techniques, followed by molecular dynamics simulation analysis, we identified two potent drugs (Tadalafil and Finasteride) that can efficiently inhibit AKT1 and MDM2 proteins. We propose that these two drugs could be considered for the development of therapeutic strategies for the management of aggressive melanoma.

Exploring and validating the metastasis mechanism of pathenolide interfering with cutaneous melanoma through ER stress-dependent apoptosis based on the network pharmacology

BACKGROUND

This study aimed to explore the mechanism of parthenolide in inhibiting melanoma metastasis through network pharmacology and cell experiment.

MATERIALS AND METHODS

This research obtained the targets of the drug from the HERB database and PubChem database, the differential expression gene of metastatic cutaneous melanoma was obtained by differentially expression gene analysis of four Gene Expression Omnibus (GEO) datasets. The intersection of drug targets and differentially expression genes were considered to be related to drugs that inhibit metastasis of cutaneous melanoma. The STRING database was used to construct the protein-protein interaction (PPI) network, and cytohubba package in Cytoscape software was used to rank the PPI network targets. The enrichment analysis was used to screen out the relevance Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and gene ontology to explain the underlying mechanism of drug inhibiting the cutaneous melanoma metastatic; cell viability, apoptosis, cell migration and protein levels were assessed using cell counting kit-8 (CCK-8) assay, Annexin V-FITC/PI assay, wound healing assay, respectively. Finally, combining pathway maps and literature, we detected ATF4 and proteins upstream and downstream of ATF4 through Western blot.

RESULTS

A total of 87 targets were screened out from the drug databases, and a total of 1635 differentially expression genes was obtained from the differentially expression genes analysis of GEO datasets, a total of nine targets (VEGFA, ANXA5, ICAM1, SELE, NFKBIA, ATF4, CTNNB1, SELP and HPGDS) were considered to be related to drugs that inhibit metastasis of cutaneous melanoma. The result of enrichment analysis showed that the drug inhibits the metastatic of cutaneous melanoma through multiple pathways such as TNF signalling pathway, lipid and atherosclerosis and fluid shear stress and atherosclerosis, relevance multiple biological processes, cellular components and molecular function; cell experiments showed that parthenolide could inhibit tumour cell migration and induce a decrease of cell viability. Flow cytometry results showed that parthenolide induced tumour cell apoptosis. Western blot results suggested that parthenolide exerted therapeutic effects by regulating ATF4 protein and its upstream and downstream proteins, namely endoplasmic reticulum (ER) stress signalling pathway.

CONCLUSION

Parthenolide induces ER stress-dependent apoptosis in melanoma cells.

Combining single-cell and transcriptomic analysis revealed the immunomodulatory effect of GOT2 on a glutamine-dependent manner in cutaneous melanoma

Background: Reprogramming in glutamine metabolism is a hallmark of cancers, while its role in cutaneous melanoma has not been studied at great length. Methods: Here, we constructed a glutamine metabolism-related prognostic signature in cutaneous melanoma with a variety of bioinformatics methods according to the glutamine metabolism regulatory molecules. Moreover, experimental verification was carried out for the key gene. Results: We have identified two subgroups of cutaneous melanoma patients, each with different prognoses, immune characteristics, and genetic mutations. GOT2 was the most concerned key gene among the model genes. We verified its role in promoting tumor cell proliferation by CCK-8 and clone formation assays. Conclusion: Our study cast new light on the prognosis of cutaneous melanoma, and the internal mechanism regulating glutamine metabolism of GOT2 may provide a new avenue for treating the cutaneous melanoma disease precisely.

A Novel Method to Detect Copy Number Variation in Melanoma: Droplet Digital PCR for Quantitation of the CDKN2A Gene, a Proof-of-Concept Study

ABSTRACT

A definitive diagnosis of nevus or melanoma is not always possible for histologically ambiguous melanocytic neoplasms. In such cases, ancillary molecular testing can support a diagnosis of melanoma if certain chromosomal aberrations are detected. Current technologies for copy number variation (CNV) detection include chromosomal microarray analysis (CMA) and fluorescence in situ hybridization. Although CMA and fluorescence in situ hybridization are effective, their utilization can be limited by cost, turnaround time, and inaccessibility outside of large reference laboratories. Droplet digital polymerase chain reaction (ddPCR) is a rapid, automated, and relatively inexpensive technology for CNV detection. We investigated the ability of ddPCR to quantify CNV in cyclin-dependent kinase inhibitor 2A ( CDKN2A ), the most commonly deleted tumor suppressor gene in melanoma. CMA data were used as the gold standard. We analyzed 57 skin samples from 52 patients diagnosed with benign nevi, borderline lesions, primary melanomas, and metastatic melanomas. In a training cohort comprising 29 randomly selected samples, receiver operator characteristic curve analysis revealed an optimal ddPCR cutoff value of 1.73 for calling CDKN2A loss. In a validation cohort comprising the remaining 28 samples, ddPCR detected CDKN2A loss with a sensitivity and specificity of 94% and 90%, respectively. Significantly, ddPCR could also identify whether CDKN2A losses were monoallelic or biallelic. These pilot data suggest that ddPCR can detect CDKN2A deletions in melanocytic tumors with accuracy comparable with CMA. With further validation, ddPCR could provide an additional CNV assay to aid in the diagnosis of challenging melanocytic neoplasms.

Transcriptome Analysis Identifies Oncogenic Tissue Remodeling during Progression from Common Nevi to Early Melanoma

Early detection and treatment of melanoma, the most aggressive skin cancer, improves the median 5-year survival rate of patients from 25% to 99%. Melanoma development involves a stepwise process during which genetic changes drive histologic alterations within nevi and surrounding tissue. Herein, a comprehensive analysis of publicly available gene expression data sets of melanoma, common or congenital nevi (CN), and dysplastic nevi (DN), assessed molecular and genetic pathways leading to early melanoma. The results demonstrate several pathways reflective of ongoing local structural tissue remodeling activity likely involved during the transition from benign to early-stage melanoma. These processes include the gene expression of cancer-associated fibroblasts, collagens, extracellular matrix, and integrins, which assist early melanoma development and the immune surveillance that plays a substantial role at this early stage. Furthermore, genes up-regulated in DN were also overexpressed in melanoma tissue, supporting the notion that DN may serve as a transitional phase toward oncogenesis. CN collected from healthy individuals exhibited different gene signatures compared with histologically benign nevi tissue located adjacent to melanoma (adjacent nevi). Finally, the expression profile of microdissected adjacent nevi tissue was more similar to melanoma compared with CN, revealing the melanoma influence on this annexed tissue.

Proteomic profiling of cutaneous melanoma explains the aggressiveness of distant organ metastasis

Despite recent developments in managing metastatic melanomas, patients' overall survival remains low. Therefore, the current study aims to understand better the proteome-wide changes associated with melanoma metastasis that will assist with identifying targeted therapies. The latest development in mass spectrometry-based proteomics, together with extensive bioinformatics analysis, was used to investigate the molecular changes in 60 formalin-fixed and paraffin-embedded samples of primary and lymph nodes (LN) and distant organ metastatic melanomas. A total of 4631 proteins were identified, of which 72 and 453 were significantly changed between the LN and distant organ metastatic melanomas compared to the primary lesions (adj. p-value <0.05). An increase in proteins such as SLC9A3R1, CD20 and GRB2 and a decrease in CST6, SERPINB5 and ARG1 were associated with regional LN metastasis. By contrast, increased metastatic activities in distant organ metastatic melanomas were related to higher levels of CEACAM1, MC1R, AKT1 and MMP3-9 and decreased levels of CDKN2A, SDC1 and SDC4 proteins. Furthermore, machine learning analysis classified the lesions with up to 92% accuracy based on their metastatic status. The findings from this study provide up to date proteome-level information about the progression of melanomas to regional LN and distant organs, leading to the identification of protein signatures with potential for clinical translation.

Age-related next-generation sequencing mutational analysis in 1196 melanomas

BACKGROUND AND OBJECTIVES

Melanoma mutational burden is high and approximately 50% have oncogenic mutations in BRAF. We sought to evaluate age-related mutational differences in melanoma.

METHODS

We analyzed melanoma samples in the Genomics Evidence Neoplasia Information Exchange database. Targetable mutations were identified using the Precision Oncology Knowledge Base (OncoKB).

RESULTS

We found 1194 patients with a common set of 30 genes. The top mutated genes in patients <40 years old (y/o) (n = 98) were BRAF (59%), TP53 (31%), NRAS (17%), and PTEN (14%); in 40-59 y/o (n = 354) were BRAF (51%), NRAS (30%), TP53 (26%), and APC (13%); and in ≥60 y/o (n = 742) were BRAF (38%), NRAS (33%), TP53 (26%), and KDR (19%). BRAF mutations were almost mutually exclusive from NRAS mutations in <40 y/o (58/59). Mutational burden increased with age, with means of 2.39, 2.92, and 3.67 mutations per sample in patients <40, 40-59, and ≥60 y/o, respectively (p < 0.0001). There were 10 targetable mutations meeting OncoKB criteria for melanoma: BRAF (level 1), RET (level 1), KIT (level 2), NRAS (level 3A), TP53 (level 3A), and FGFR2, MET, PTEN, PIK3CA, and KRAS (level 4).

CONCLUSIONS

Mutations in melanoma have age-related differences and demonstrates potential targetable mutations for personalized therapies.

SMG7-AS1 as a prognostic biomarker and predictor of immunotherapy responses for skin cutaneous melanoma

Skin cutaneous melanoma (SKCM) is the most life-threatening skin cancer and lacks early detection and effective treatment strategies. Many long noncoding RNAs are associated with the development of tumors and may serve as potential immunotherapeutic targets. In this study, microarray analysis was performed to screen for differentially expressed lncRNAs between SKCM and normal tissues, and SMG7-AS1 was identified as an upregulated lncRNA in SKCM. Subsequently, bioinformatic analysis revealed that dysregulation of SMG7-AS1 influences metastasis and immune infiltration. qRT-PCR of clinical samples demonstrated that the expression of SMG7-AS1 was higher in melanoma tissues. Flow cytometry showed that SMG7-AS1 plays a vital role in the cell cycle. Additionally, SMG7-AS1 was found to be associated with immunotherapy responses. To the best of our knowledge, this study is the first to report that SMG7-AS1 is associated with SKCM and may serve as a prognostic biomarker and predictor of immunotherapy responses in SKCM.

Collectin-11 promotes cancer cell proliferation and tumor growth

Collectin-11 (CL-11) is a recently described soluble C-type lectin that has distinct roles in embryonic development, host defence, autoimmunity, and fibrosis. Here we report that CL-11 also plays an important role in cancer cell proliferation and tumor growth. Melanoma growth was found to be suppressed in Colec11-/- mice in a s.c. B16 melanoma model. Cellular and molecular analyses revealed that CL-11 is essential for melanoma cell proliferation, angiogenesis, establishment of more immunosuppressive tumor microenvironment, and the reprogramming of macrophages to M2 phenotype within melanomas. In vitro analysis revealed that CL-11 can activate tyrosine kinase receptors (EGFR, HER3) and ERK, JNK, and AKT signaling pathways and has a direct stimulatory effect on murine melanoma cell proliferation. Furthermore, blockade of CL-11 (treatment with L-fucose) inhibited melanoma growth in mice. Analysis of open data sets revealed that COLEC11 gene expression is upregulated in human melanomas and that high COLEC11 expression has a trend toward poor survival. CL-11 also had direct stimulatory effects on human tumor cell proliferation in melanoma and several other types of cancer cells in vitro. Overall, our findings provide the first evidence to our knowledge that CL-11 is a key tumor growth-promoting protein and a promising therapeutic target in tumor growth.

Ethanol Metabolism and Melanoma

Malignant melanoma is the deadliest form of skin cancer. Despite significant efforts in sun protection education, melanoma incidence is still rising globally, drawing attention to other socioenvironmental risk factors for melanoma. Ethanol and acetaldehyde (AcAH) are ubiquitous in our diets, medicines, alcoholic beverages, and the environment. In the liver, ethanol is primarily oxidized to AcAH, a toxic intermediate capable of inducing tumors by forming adducts with proteins and DNA. Once in the blood, ethanol and AcAH can reach the skin. Although, like the liver, the skin has metabolic mechanisms to detoxify ethanol and AcAH, the risk of ethanol/AcAH-associated skin diseases increases when the metabolic enzymes become dysfunctional in the skin. This review highlights the evidence linking cutaneous ethanol metabolism and melanoma. We summarize various sources of skin ethanol and AcAH and describe how the reduced activity of each alcohol metabolizing enzyme affects the sensitivity threshold to ethanol/AcAH toxicity. Data from the Gene Expression Omnibus database also show that three ethanol metabolizing enzymes (alcohol dehydrogenase 1B, P450 2E1, and catalase) and an AcAH metabolizing enzyme (aldehyde dehydrogenase 2) are significantly reduced in melanoma tissues.

The mechanosensitive TRPV2 calcium channel promotes human melanoma invasiveness and metastatic potential

Melanoma is a highly aggressive cancer endowed with a unique capacity of rapidly metastasizing, which is fundamentally driven by aberrant cell motility behaviors. Discovering "migrastatics" targets, specifically controlling invasion and dissemination of melanoma cells during metastasis, is therefore of primary importance. Here, we uncover the prominent expression of the plasma membrane TRPV2 calcium channel as a distinctive feature of melanoma tumors, directly related to melanoma metastatic dissemination. In vitro as well as in vivo, TRPV2 activity is sufficient to confer both migratory and invasive potentials, while conversely TRPV2 silencing in highly metastatic melanoma cells prevents aggressive behavior. In invasive melanoma cells, TRPV2 channel localizes at the leading edge, in dynamic nascent adhesions, and regulates calcium-mediated activation of calpain and the ensuing cleavage of the adhesive protein talin, along with F-actin organization. In human melanoma tissues, TRPV2 overexpression correlates with advanced malignancy and poor prognosis, evoking a biomarker potential. Hence, by regulating adhesion and motility, the mechanosensitive TRPV2 channel controls melanoma cell invasiveness, highlighting a new therapeutic option for migrastatics in the treatment of metastatic melanoma.

Impact of Green Gold Nanoparticle Coating on Internalization, Trafficking, and Efficiency for Photothermal Therapy of Skin Cancer

Skin cancer is a global health issue and mainly composed of melanoma and nonmelanoma cancers. For the first clinical proof of concept on humans, we decided to study good prognosis skin cancers, i.e., carcinoma basal cell. In UE, the first-line treatment remains surgical resection, healing most of the tumors, but presents aesthetic disadvantages with a high reoccurrence rate on exposed areas. Moreover, the therapeutic indications could extend to melanoma and metastasis, which is a different medical strategy that could combine this treatment. Indeed, patients with late-stage melanoma are in a therapeutic impasse, despite recent targeted and immunological therapies. Photothermal therapy using gold nanoparticles is the subject of many investigations due to their strong potential to treat cancers by physical, thermal destruction. We developed gold nanoparticles synthesized by green chemistry (gGNPs), using endemic plant extract from Reunion Island, which have previously showed their efficiency at a preclinical stage. Here, we demonstrate that these gGNPs are less cytotoxic than gold nanoparticles synthesized by Turkevich's method. Furthermore, our work describes the optimization of gGNP coating and stabilization, also taking into consideration the gGNP path in cells (endocytosis, intracellular trafficking, and exocytosis), their specificity toward cancerous cells, their cytotoxicity, and their in vivo efficiency. Finally, based on the metabolic switch of cancerous cells overexpressing Glut transporters in skin cancers, we demonstrated that glucose-stabilized gGNP (gGNP@G) enables a quick internalization, fourfold higher in cancerous cells in contrast to healthy cells with no side cytotoxicity, which is particularly relevant to target and treat cancer.

LAP1 supports nuclear adaptability during constrained melanoma cell migration and invasion

Metastasis involves dissemination of cancer cells away from a primary tumour and colonization at distal sites. During this process, the mechanical properties of the nucleus must be tuned since they pose a challenge to the negotiation of physical constraints imposed by the microenvironment and tissue structure. We discovered increased expression of the inner nuclear membrane protein LAP1 in metastatic melanoma cells, at the invasive front of human primary melanoma tumours and in metastases. Human cells express two LAP1 isoforms (LAP1B and LAP1C), which differ in their amino terminus. Here, using in vitro and in vivo models that recapitulate human melanoma progression, we found that expression of the shorter isoform, LAP1C, supports nuclear envelope blebbing, constrained migration and invasion by allowing a weaker coupling between the nuclear envelope and the nuclear lamina. We propose that LAP1 renders the nucleus highly adaptable and contributes to melanoma aggressiveness.

Biomarker Discovery for Meta-Classification of Melanoma Metastatic Progression Using Transfer Learning

Melanoma is considered to be the most serious and aggressive type of skin cancer, and metastasis appears to be the most important factor in its prognosis. Herein, we developed a transfer learning-based biomarker discovery model that could aid in the diagnosis and prognosis of this disease. After applying it to the ensemble machine learning model, results revealed that the genes found were consistent with those found using other methodologies previously applied to the same TCGA (The Cancer Genome Atlas) data set. Further novel biomarkers were also found. Our ensemble model achieved an AUC of 0.9861, an accuracy of 91.05, and an F1 score of 90.60 using an independent validation data set. This study was able to identify potential genes for diagnostic classification (C7 and GRIK5) and diagnostic and prognostic biomarkers (S100A7, S100A7, KRT14, KRT17, KRT6B, KRTDAP, SERPINB4, TSHR, PVRL4, WFDC5, IL20RB) in melanoma. The results show the utility of a transfer learning approach for biomarker discovery in melanoma.

Melanoma: Molecular genetics, metastasis, targeted therapies, immunotherapies, and therapeutic resistance

Cutaneous melanoma is a common cancer and cases have steadily increased since the mid 70s. For some patients, early diagnosis and surgical removal of melanomas is lifesaving, while other patients typically turn to molecular targeted therapies and immunotherapies as treatment options. Easy sampling of melanomas allows the scientific community to identify the most prevalent mutations that initiate melanoma such as the BRAF, NRAS, and TERT genes, some of which can be therapeutically targeted. Though initially effective, many tumors acquire resistance to the targeted therapies demonstrating the need to investigate compensatory pathways. Immunotherapies represent an alternative to molecular targeted therapies. However, inter-tumoral immune cell populations dictate initial therapeutic response and even tumors that responded to treatment develop resistance in the long term. As the protocol for combination therapies develop, so will our scientific understanding of the many pathways at play in the progression of melanoma. The future direction of the field may be to find a molecule that connects all of the pathways. Meanwhile, noncoding RNAs have been shown to play important roles in melanoma development and progression. Studying noncoding RNAs may help us to understand how resistance - both primary and acquired - develops; ultimately allow us to harness the true potential of current therapies. This review will cover the basic structure of the skin, the mutations and pathways responsible for transforming melanocytes into melanomas, the process by which melanomas metastasize, targeted therapeutics, and the potential that noncoding RNAs have as a prognostic and treatment tool.

Genomic landscape of the immunogenicity regulation in skin melanomas with diverse tumor mutation burden

Skin melanoma cells are tightly interconnected with their tumor microenvironment (TME), which influences their initiation, progression, and sensitivity/resistance to therapeutic interventions. An immune-active TME favors patient response to immune checkpoint inhibition (ICI), but not all patients respond to therapy. Here, we assessed differential gene expression in primary and metastatic tumors from the TCGA-SKCM dataset, compared to normal skin samples from the GTEx project and validated key findings across 4 independent GEO datasets, as well as using immunohistochemistry in independent patient cohorts. We focused our attention on examining the expression of various immune receptors, immune-cell fractions, immune-related signatures and mutational signatures across cutaneous melanomas with diverse tumor mutation burdens (TMB). Globally, the expression of most immunoreceptors correlated with patient survival, but did not differ between TMBhigh and TMBlow tumors. Melanomas were enriched in "naive T-cell", "effector memory T-cell", "exhausted T-cell", "resting Treg T-cell" and "Th1-like" signatures, irrespective of their BRAF, NF1 or RAS mutational status. Somatic mutations in IDO1 and HLA-DRA were frequent and could be involved in hindering patient response to ICI therapies. We finally analyzed transcriptome profiles of ICI-treated patients and associated their response with high levels of IFNγ, Merck18, CD274, CD8, and low levels of myeloid-derived suppressor cells (MDSCs), cancer-associated fibroblasts (CAFs) and M2 macrophages, irrespective of their TMB status. Overall, our findings highlight the importance of pre-existing T-cell immunity in ICI therapeutic outcomes in skin melanoma and suggest that TMBlow patients could also benefit from such therapies.

CEBPB is associated with active tumor immune environment and favorable prognosis of metastatic skin cutaneous melanoma

Metastatic skin cutaneous melanoma (SKCM) is a common malignancy that accounts for low morbidity but high mortality of skin cancer. SKCM is characterized by high lymphocytic infiltration, whereas the states of infiltrated cells are variable in patients leading to a heterogeneous prognosis and hindering appropriate clinical decisions. It is therefore urgent to identify markers associated with lymphocytic infiltration, cellular conditions, and the prognosis of SKCM. In this study, we report that CEBPB, a transcriptional factor, is mainly expressed in macrophages in metastatic SKCM and associated with an active tumor immune environment and a favorable prognosis through integrated analysis of single-cell and bulk RNA-seq datasets. High CEBPB expression is significantly associated with active inflammation and immune response pathways in both macrophages and bulk SKCM tumor tissues. A signature based on CEBPB-associated genes that are specifically expressed in macrophages could robustly and prognostically separate different metastatic SKCM patients. In addition, the associations between the metastatic SKCM tumor signature and microenvironment with respect to T-cell recruitment and state, inflammation response, angiogenesis, and so on were also determined. In conclusion, we present here the first report on CEBPB in tumor immune environment and prognosis regulation in metastatic SKCM and construct a reliable signature, which should provide a useful biomarker for stratification of the patient's prognosis and therapeutic selection.

Exploration and validation of metastasis-associated genes for skin cutaneous melanoma

Skin cutaneous melanoma is a malignant and highly metastatic skin tumor, and its morbidity and mortality are still rising worldwide. However, the molecular mechanisms that promote melanoma metastasis are unclear. Two datasets (GSE15605 and GSE46517) were retrieved to identify the differentially expressed genes (DEGs), including 23 normal skin tissues (N), 77 primary melanoma tissues (T) and 85 metastatic melanoma tissues (M). Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were performed to explore the functions of the DEGs. We constructed protein–protein interaction network using the STRING database and Cytoscape software. Using the cytoHubba plugin of Cytoscape, we identified the most significant hub genes by five analytical methods (Degree, Bottleneck, MCC, MNC, and EPC). Hub gene expression was validated using the UALCAN website. Clinical relevance was investigated using The Cancer Genome Atlas resources. Finally, we explored the association between metastasis-associated genes and immune infiltrates through the Tumor Immune Estimation Resource (TIMER) database and performed drug–gene interaction analysis using the Drug-Gene Interaction database. A total of 294 specific genes were related to melanoma metastasis and were mainly involved in the positive regulation of locomotion, mitotic cell cycle process, and epithelial cell differentiation. Four hub genes (CDK1, FOXM1, KIF11, and RFC4) were identified from the cytoHubba plugin of Cytoscape. CDK1 was significantly upregulated in metastatic melanoma compared with primary melanoma, and high CDK1 expression was positively correlated with worse overall survival. Immune infiltration analysis revealed that CDK1 expression negatively correlated with macrophage infiltration (Rho = − 0.164, P = 2.02e−03) and positively correlated with neutrophil cells (Rho = 0.269, P = 2.72e−07) in SKCM metastasis. In addition, we identified that CDK1 had a close interaction with 10 antitumor drugs. CDK1 was identified as a hub gene involved in the progression of melanoma metastasis and may be regarded as a therapeutic target for melanoma patients to improve prognosis and prevent metastasis in the future.

Identification of a Seven-Differentially Expressed Gene-Based Recurrence-Free Survival Model for Melanoma Patients

Melanoma is a malignant tumor that originates in melanocytes of the skin or mucous membrane, which has a high mortality rate and worse prognosis. Therefore, perspective prognosis evaluation seems more important for patients' treatment. Gene expression profiles of melanoma were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, respectively. 130 consistent differentially expressed genes (DEGs) were identified between melanoma and nevus tissues from two GEO cohorts. Prognostic genes were identified by univariate analysis, and 20 of them were regarded to be associated with the recurrence-free survival (RFS) of melanoma patients. Then, the LASSO Cox regression analysis chose seven of them to establish a seven-DEG-based RFS predicting signature. We demonstrated that this model was more powerful to predict RFS risk than other individual clinical features and was able to independently predict the RFS outcomes in different subsets of patients. We attempted to search for the underlying mechanisms by analyzing the coexpression genes of the seven candidates, and the pathway enrichment analyses indicated that immune response-related pathways might play a critical role in melanoma progression. Finally, we establish a robust seven-DEG-based RFS predicting signature, which will facilitate the personalized treatment of melanoma patients.

NRF3 Decreases during Melanoma Carcinogenesis and Is an Independent Prognostic Marker in Melanoma

The prognostic significance of the major redox regulator, nuclear factor erythroid-2-related factor 2 (NRF2), is recognized in many cancers, but the role of NRF3 is not studied. Analysis from the Gene Expression Omnibus datasets showed that NRF3 mRNA levels increased from benign to dysplastic naevi (p = 0.04). We characterized the immunohistochemical expression of NRF3 in 81 naevi, 67 primary skin melanomas, and 51 lymph node metastases. The immunohistochemical expression of cytoplasmic NRF3 decreased from benign to dysplastic naevi (p < 0.001) and further to primary melanomas (p < 0.001). High cytoplasmic NRF3 protein expression in pigment cells of the primary melanomas associated with worse melanoma-specific survival in multivariate analysis, specifically in the subgroup of patients with the lymph node metastases at the time of diagnosis (hazard ratio 3.179; 95% confidence interval 1.065-9.493; p = 0.038). Intriguingly, we did not observe associations between NRF3 and the traditional prognostic factors such as Breslow thickness, ulceration, or stage. Together, this data represents the primary description about the role of NRF3 in pigment tumours that is worthy of further explorations.

The histone demethylase PHF8 regulates TGFβ signaling and promotes melanoma metastasis

The contribution of epigenetic dysregulation to metastasis remains understudied. Through a meta-analysis of gene expression datasets followed by a mini-screen, we identified Plant Homeodomain Finger protein 8 (PHF8), a histone demethylase of the Jumonji C protein family, as a previously unidentified prometastatic gene in melanoma. Loss- and gain-of-function approaches demonstrate that PHF8 promotes cell invasion without affecting proliferation in vitro and increases dissemination but not subcutaneous tumor growth in vivo, thus supporting its specific contribution to the acquisition of metastatic potential. PHF8 requires its histone demethylase activity to enhance melanoma cell invasion. Transcriptomic and epigenomic analyses revealed that PHF8 orchestrates a molecular program that directly controls the TGFβ signaling pathway and, as a consequence, melanoma invasion and metastasis. Our findings bring a mechanistic understanding of epigenetic regulation of metastatic fitness in cancer, which may pave the way for improved therapeutic interventions.

Identification of Keratinocyte Differentiation-Involved Genes for Metastatic Melanoma by Gene Expression Profiles

Background

Melanoma is the deadliest type of skin cancer. Until now, its pathological mechanisms, particularly the mechanism of metastasis, remain largely unknown. Our study on the identification of genes in association with metastasis for melanoma provides a novel understanding of melanoma.

Methods

From the Gene Expression Omnibus (GEO) database, the gene expression microarray datasets GSE46517, GSE7553, and GSE8401 were downloaded. We made use of R aiming at analyzing the differentially expressed genes (DEGs) between metastatic and nonmetastatic melanoma. R was also used in differentially expressed miRNA (DEM) data mining from GSE18509, GSE19387, GSE24996, GSE34460, GSE35579, GSE36236, and GSE54492 datasets referring to Li's study. Based on the DEG and DEM data, we performed functional enrichment analysis through the application of the DAVID database. Furthermore, we constructed the protein-protein interaction (PPI) network and established functional modules by making use of the STRING database. Through making use of Cytoscape, the PPI results were visualized. We predicted the targets of the DEMs through applying TargetScan, miRanda, and PITA databases and identified the overlapping genes between DEGs and predicted targets, followed by the construction of DEM-DEG pair network. The expressions of these keratinocyte differentiation-involved genes in Module 1 were identified based on the data from TCGA.

Results

239 DEGs were screened out in all 3 datasets, which were inclusive of 21 positively regulated genes and 218 negatively regulated genes. Based on these 239 DEGs, we finished constructing the PPI network which was formed from 225 nodes and 846 edges. We finished establishing 3 functional modules. And we analyzed 92 overlapping genes and 26 miRNA, including 11 upregulated genes targeted by 11 negatively regulated DEMs and 81 downregulated genes targeted by 15 positively regulated DEMs. As proof of the differential expression of metastasis-associated genes, eleven keratinocyte differentiation-involved genes, including LOR, EVPL, SPRR1A, FLG, SPRR1B, SPRR2B, TGM1, DSP, CSTA, CDSN, and IVL in Module 1, were obviously downregulated in metastatic melanoma tissue in comparison with primary melanoma tissue based on the data from TCGA.

Conclusion

239 melanoma metastasis-associated genes and 26 differentially expressed miRNA were identified in our study. The keratinocyte differentiation-involved genes may take part in melanoma metastasis, providing a latent molecular mechanism for this disease.

AKAP9 Upregulation Predicts Unfavorable Prognosis in Pediatric Acute Myeloid Leukemia and Promotes Stemness Properties via the Wnt/β-Catenin Pathway

Background

PRKA kinase anchor protein 9 (AKAP9) is a scaffold protein involved in various cellular processes, including cell adhesion, proliferation, differentiation, and apoptosis. Although the oncogenic role of AKAP9 in solid tumors is well elucidated, the functions and mechanisms of AKAP9 in acute myeloid leukemia (AML) are still not understood.

Methods

We used the gene expression omnibus (GEO) database (GSE2191) to determine the mRNA expression of AKAP9 in the bone marrow of pediatric AML and healthy patients. We further used the therapeutically available research to generate effective treatments (TARGET) database to elucidate the relationship between AKAP9 expression and clinical outcomes in pediatric patients with AML. In addition, cell proliferation, cell cycle, apoptosis, RT-PCR, and Western blotting assays were applied to reveal the functions of AKAP9 and the underlying mechanisms of AKAP9 silencing in THP1 and HL60 cell lines.

Results

AKAP9 is overexpressed in the bone marrow of pediatric AML patients as compared with that of healthy patients. High expression of AKAP9 was found to be a predictor of poor overall survival (OS) and event-free survival (EFS). Using univariate and multivariate survival analyses, we found that high AKAP9 expression is an independent predictor of a worse OS and EFS. Functionally, AKAP9 silencing significantly inhibited AML cell proliferation, and cell cycle progression and promoted apoptosis. Moreover, AKAP9 silencing significantly downregulated the expression of stemness markers and β-catenin.

Conclusion

AKAP9 upregulation is a predictor of unfavorable prognosis, promotes stemness, and activates the Wnt/β-catenin pathway in AML patients. AKAP9 may act as a prognostic biomarker of AML in pediatric patients and a future therapeutic target.

Genomic and Transcriptomic Underpinnings of Melanoma Genesis, Progression, and Metastasis

Melanoma is a deadly skin cancer with rapidly increasing incidence worldwide. The discovery of the genetic drivers of melanomagenesis in the last decade has led the World Health Organization to reclassify melanoma subtypes by their molecular pathways rather than traditional clinical and histopathologic features. Despite this significant advance, the genomic and transcriptomic drivers of metastatic progression are less well characterized. This review describes the known molecular pathways of cutaneous and uveal melanoma progression, highlights recently identified pathways and mediators of metastasis, and touches on the influence of the tumor microenvironment on metastatic progression and treatment resistance. While targeted therapies and immune checkpoint blockade have significantly aided in the treatment of advanced disease, acquired drug resistance remains an unfortunately common problem, and there is still a great need to identify potential prognostic markers and novel therapeutic targets to aid in such cases.

Identification of Genes Potentially Associated with Melanoma Tumorigenesis Through Co-Expression Network Analysis

Background

Melanoma is one of the most malignant and aggressive skin tumors, and its incidence is increasing worldwide. However, few studies have investigated the process of tumorigenesis from normal skin to melanoma.

Methods

Several bioinformatics analyses, including GEO databases, Oncomine database, TCGA database, STRING, MCODE and cytoHubba plug-ins, GEPIA, TIMER and TRRUST and DGIdb, were performed to disclose the hub genes and immunology implicated in primary melanoma tumorigenesis. Finally, immunohistochemistry (IHC) and quantitative real-time PCR (qPCR) were used to validate the results of bioinformatics analysis in vitro.

Results

A total of 295 overlapping DEGs (ODEGs) (157 upregulated and 138 downregulated) and 9 hub genes were identified between primary melanoma and normal skin tissues. Functional analysis of these 9 hub genes indicated that the genes were primarily enriched in cell chemotaxis, the chemokine-mediated signaling pathway, the extracellular region, the extracellular space, chemokine activity and CXCR3 chemokine receptor binding. KEGG pathway enrichment showed that these genes were primarily involved in the chemokine signaling pathway, cytokine–cytokine receptor interaction, the toll-like receptor signaling pathway, the cytosolic DNA-sensing pathway and the TNF signaling pathway. Upregulated CCL5, CCL4, CXCL9 and CXCL10 demonstrated good overall survival (OS), and most of them have a higher expression in stage 0 and 1 of melanoma. Moreover, immune infiltration analysis showed that the above hub genes showed a strong positive correlation between their expression and infiltration of the six immune cell subsets. Transcription factor regulation network suggested that RELA and NFKB1 are the transcription factors of CCL4, CCL5, CXCL10 and CXCL2, while IRF7, IRF3 and IRF1 are the transcription factors of CCL5 and CXCL10. Drug–gene interaction analysis identified 46 drug–gene interactions. In vitro data demonstrated that the level of CCL4, CCL5, CXCL9 and CXCL10 is higher in melanoma than that in normal skin tissues, either at tissue or cell lines level.

Conclusion

In summary, we identified 4 key chemokine members related to tumorigenesis and progression in primary melanoma, and these results may help to elucidate melanoma tumorigenesis and facilitate its treatment.

A preclinical pipeline to evaluate migrastatics as therapeutic agents in metastatic melanoma

BACKGROUND

Metastasis is a hallmark of cancer and responsible for most cancer deaths. Migrastatics were defined as drugs interfering with all modes of cancer cell invasion and thus cancers' ability to metastasise. First anti-metastatic treatments have recently been approved.

METHODS

We used bioinformatic analyses of publicly available melanoma databases. Experimentally, we performed in vitro target validation (including 2.5D cell morphology analysis and mass spectrometric analysis of RhoA binding partners), developed a new traceable spontaneously metastasising murine melanoma model for in vivo validation, and employed histology (haematoxylin/eosin and phospho-myosin II staining) to confirm drug action in harvested tumour tissues.

RESULTS

Unbiased and targeted bioinformatic analyses identified the Rho kinase (ROCK)-myosin II pathway and its various components as potentially relevant targets in melanoma. In vitro validation demonstrated redundancy of several RhoGEFs upstream of RhoA and confirmed ROCK as a druggable target downstream of RhoA. The anti-metastatic effects of two ROCK inhibitors were demonstrated through in vivo melanoma metastasis tracking and inhibitor effects also confirmed ex vivo by digital pathology.

CONCLUSIONS

We proposed a migrastatic drug development pipeline. As part of the pipeline, we provide a new traceable spontaneous melanoma metastasis model for in vivo quantification of metastasis and anti-metastatic effects by non-invasive imaging.

GJB5 association with BRAF mutation and survival in cutaneous malignant melanoma

BACKGROUND

Gap junctional intercellular communication is crucial for epidermal cellular homeostasis. Inability to establish melanocyte-keratinocytes contacts and loss of intercellular junction's integrity may contribute to melanoma development. Connexins, laminins and desmocollins have been implicated in the control of melanoma growth, where their reduced expression has been reported in metastatic lesions.

OBJECTIVES

The aim of this study was to investigate Connexin 31.1 (GJB5) expression and identify any association with BRAF mutational status, melanoma patient prognosis and MAPK inhibitors (MAPKi) treatment.

MATERIAL AND METHODS

GJB5 expression was measured at RNA and protein level in melanoma clinical samples and established cell lines treated or not with BRAF and MEK inhibitors, as well as in cell lines which developed MAPK inhibitors resistance. Findings were further validated and confirmed by analysis of independent datasets.

RESULTS

Our analysis reveals significant downregulation of GJB5 expression in metastatic melanoma lesions compared to primary ones and in BRAF mutated versus BRAF wild-type melanomas. Likewise, GJB5 expression is significantly lower in BRAF^V600E compared with BRAF^WT cell lines and increases upon MAPKi treatment. MAPKi-resistant melanoma cells display a similar expression pattern compared to BRAF^WT cells, with increased GJB5 expression associated with morphological changes. Enhancement of BRAF^V600E expression in BRAF^WT melanoma cells significantly upregulates miR-335-5p expression with consequent downregulation of GJB5, one of its targets. Furthermore, overexpression of miR-335-5p in two BRAF^WT cell lines confirms specific GJB5 protein downregulation. RT-qPCR analysis also revealed upregulation of miR-335 in BRAF^V600E melanoma cells, which is significantly downregulated in cells resistant to MEK inhibitors. Our data were further validated using the TCGA-SKCM dataset, where BRAF mutations associate with increased miR-335 expression and inversely correlate with GJB5 expression. In clinical samples, GJB5 underexpression is also associated with patient overall worse survival, especially at early stages.

CONCLUSION

We identified a significant association between metastases / BRAF mutation and low GJB5 expression in melanoma. Our results identify a novel mechanism of Gap-junctional protein regulation, suggesting a prognostic role for GJB5 in cutaneous melanoma.

Simultaneous silencing of lysophosphatidylcholine acyltransferases 1-4 by nucleic acid nanoparticles (NANPs) improves radiation response of melanoma cells

Radiation induces the generation of platelet-activating factor receptor (PAF-R) ligands, including PAF and oxidized phospholipids. Alternatively, PAF is also synthesized by the biosynthetic enzymes lysophosphatidylcholine acyltransferases (LPCATs) which are expressed by tumor cells including melanoma. The activation of PAF-R by PAF and oxidized lipids triggers a survival response protecting tumor cells from radiation-induced cell death, suggesting the involvement of the PAF/PAF-R axis in radioresistance. Here, we investigated the role of LPCATs in the melanoma cell radiotherapy response. LPCAT is a family of four enzymes, LPCAT1-4, and modular nucleic acid nanoparticles (NANPs) allowed for the simultaneous silencing of all four LPCATs. We found that the in vitro simultaneous silencing of all four LPCAT transcripts by NANPs enhanced the therapeutic effects of radiation in melanoma cells by increasing cell death, reducing long-term cell survival, and activating apoptosis. Thus, we propose that NANPs are an effective strategy for improving radiotherapy efficacy in melanomas.

The Research on the Treatment of Metastatic Skin Cutaneous Melanoma by Huanglian Jiedu Decoction Based on the Analysis of Immune Infiltration Analysis

Objective

To explore the potential mechanism of Huanglian Jiedu Decoction (HJD) treatment and prevention of metastatic Cutaneous Melanoma (CM) occurrence and metastasis based on network pharmacological methods and immune infiltration analysis.

Methods

The GEO database was used to obtain metastatic CM disease targets, the TCMSP database and the HERB database were used to obtain HJD action targets, core genes were screened by protein interaction network, and the potential mechanism of HJD in the treatment of metastatic CM was explored by enrichment analysis, prognostic analysis and immune infiltration analysis.

Results

HJD treatment of metastatic CM involved 60 targets, enrichment analysis showed that HJD treatment of metastatic CM involved Chemokine signaling pathway, NF-kappa B signaling pathway, and Fluid shear stress and atherosclerosis, etc. Prognostic analysis revealed that HJD had a certain ability to improve the prognosis of metastatic CM patients. Immune infiltration analysis showed that HJD could inhibit the immune cell infiltration of metastatic CM patients by acting on related targets.

Conclusions

Our study identified the potential mechanism of HJD in the treatment of metastatic CM through network pharmacology, and revealed the mechanism of HJD in the prevention of Skin Cutaneous Melanoma metastasis through immune infiltration analysis and prognostic analysis.

Prognostic score model-based signature genes for predicting the prognosis of metastatic skin cutaneous melanoma

PURPOSE

Cutaneous melanoma (SKCM) is the most invasive malignancy of skin cancer. Metastasis to distant lymph nodes or other system is an indicator of poor prognosis in melanoma patients. The aim of this study was to identify reliable prognostic biomarkers for SKCMs.

METHODS

Four RNA-sequencing datasets associated with SKCMs were downloaded from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) database as well as corresponding clinical information. Differentially expressed genes (DEGs) were screened between primary and metastatic samples by using MetaDE tool. Weighted gene co-expression network analysis (WGCNA) was conducted to screen functional modules. A prognostic score (PS)-based predictive model and nomogram model were constructed to identify signature genes and independent clinicopathologic factors.

RESULTS

Based on MetaDE analysis and WGCNA, a total of 456 overlapped genes were identified as hub genes related to SKCMs progression. Functional enrichment analysis revealed these genes were mainly involved in the hippo signaling pathway, signaling pathways regulating pluripotency of stem cells, pathways in cancer. In addition, eight optimal DEGs (RFPL1S, CTSV, EGLN3, etc.) were identified as signature genes by using PS model. Cox regression analysis revealed that pathologic stage T, N and recurrence were independent prognostic factors. Three clinical factors and PS status were incorporated to construct a nomogram predictive model for estimating the three years and five-year survival probability of individuals.

CONCLUSIONS

The prognosis prediction model of this study may provide a promising method for decision making in clinic and prognosis predicting of SKCM patients.

Zwint facilitates melanoma progression by promoting c-Myc expression

ZW10 interactor (Zwint) is upregulated in various types of tumors and exerts a carcinogenic effect. However, little is known about the expression profile, function and molecular mechanisms of action of Zwint in melanoma. Therefore, the aim of the present study was to investigate the expression levels of Zwint in melanoma cell lines and tissues. It was revealed that Zwint was highly expressed in melanoma samples. Functional experiments indicated that Zwint knockdown suppressed the proliferation and migration of A375 melanoma cells. Further mechanistic studies demonstrated that Zwint knockdown decreased the protein expression levels of c-Myc, MMP-2, Slug, mTOR, phosphorylated (p)-mTOR, p-p38 and fibronectin, while it increased the protein expression levels of E-cadherin and MMP-9. Among these genes, c-Myc, MMP-2 and Slug were overexpressed to investigate their effects on cell proliferation following Zwint knockdown. The results demonstrated that overexpression of c-Myc, but not MMP-2 or Slug, rescued the effects of Zwint knockdown on melanoma cell proliferation and migration. Taken together, the results of the present study suggested that Zwint may act as an oncogene in melanoma by regulating c-Myc expression.

A MyD88/IL1R Axis Regulates PD-1 Expression on Tumor-Associated Macrophages and Sustains Their Immunosuppressive Function in Melanoma

Macrophages are critical mediators of tissue homeostasis, cell proliferation, and tumor metastasis. Tumor-associated macrophages (TAM) are generally associated with tumor-promoting immunosuppressive functions in solid tumors. Here, we examined the transcriptional landscape of adaptor molecules downstream of Toll-like receptors in human cancers and found that higher expression of MYD88 correlated with tumor progression. In murine melanoma, MyD88, but not Trif, was essential for tumor progression, angiogenesis, and maintaining the immunosuppressive phenotype of TAMs. In addition, MyD88 expression in myeloid cells drove melanoma progression. The MyD88/IL1 receptor (IL1R) axis regulated programmed cell death (PD)-1 expression on TAMs by promoting recruitment of NF-κBp65 to the Pdcd1 promoter. Furthermore, a combinatorial immunotherapy approach combining the MyD88 inhibitor with anti-PD-1 blockade elicited strong antitumor effects. Thus, the MyD88/IL1R axis maintains the immunosuppressive function of TAMs and promotes tumor growth by regulating PD-1 expression. SIGNIFICANCE: These findings indicate that MyD88 regulates TAM-immunosuppressive activity, suggesting that macrophage-mediated immunotherapy combining MYD88 inhibitors with PD-1 blockade could result in better treatment outcomes in a wide variety of cancers. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/81/9/2358/F1.large.jpg.

L Antigen Family Member 3 Serves as a Prognostic Biomarker for the Clinical Outcome and Immune Infiltration in Skin Cutaneous Melanoma

L Antigen Family Member 3 (LAGE3) is an important RNA modification-related protein. Whereas few studies have interrogated the LAGE3 protein, there is limited data on its role in tumors. Here, we analyzed and profiled the LAGE3 protein in skin cutaneous melanoma (CM) using TCGA, GTEx, or GEO databases. Our data showed an upregulation of LAGE3 in melanoma cell lines compared to normal skin cell lines. Besides, the Kaplan–Meier curves and Cox proportional hazard model revealed that LAGE3 was an independent survival indicator for CM, especially in metastatic CM. Moreover, LAGE3 was negatively associated with multiple immune cell infiltration levels in CM, especially CD8⁺ T cells in metastatic CM. Taken together, our study suggests that LAGE3 could be a potential prognostic biomarker and might be a potential target for the development of novel CM treatment strategies.

Long noncoding RNA LINC00518 induces radioresistance by regulating glycolysis through an miR-33a-3p/HIF-1α negative feedback loop in melanoma

The long noncoding RNA, LINC00518, is highly expressed in various types of cancers and is involved in cancer progression. Although LINC00518 promotes the metastasis of cutaneous malignant melanoma (CMM), the mechanism underlaying its effects on CMM radiosensitivity remains unclear. In this study, LINC00518 expression was significantly upregulated in CMM samples, and LINC00518 levels were associated with poor prognosis of patients with CMM. Knockdown of LINC00518 in CMM cells significantly inhibited cell invasion, migration, proliferation, and clonogenicity. LINC00518-mediated invasion, migration, proliferation, and clonogenicity were negatively regulated by the microRNA, miR-33a-3p, in vitro, which increased sensitivity to radiotherapy via inhibition of the hypoxia-inducible factor 1α (HIF-1α)/lactate dehydrogenase A glycolysis axis. Additionally, HIF-1α recognized the miR-33a-3p promoter region and recruited histone deacetylase 2, which decreased the expression of miR-33a-3p and formed an LINC00518/miR-33a-3p/HIF-1α negative feedback loop. Furthermore, signaling with initially activated glycolysis and radioresistance in CMM cells was impaired by Santacruzamate A, a histone deacetylase inhibitor, and 2-deoxy-D-glucose, a glycolytic inhibitor. Lastly, knockdown of LINC00518 expression sensitized CMM cancer cells to radiotherapy in an in vivo subcutaneously implanted tumor model. In conclusion, LINC00518 was confirmed to be an oncogene in CMM, which induces radioresistance by regulating glycolysis through an miR-33a-3p/HIF-1α negative feedback loop. Our study, may provide a potential strategy to improve the treatment outcome of radiotherapy in CMM.

Novel oxidative stress-related prognostic biomarkers for melanoma associated with tumor metastasis

ABSTRACT

Skin cutaneous melanoma (SKCM) is a prevalent skin cancer whose metastatic form is dangerous due to its high morbidity and mortality. Previous studies have systematically established the vital role of oxidative stress (OS) in melanoma progression. This study aimed to identify prognostic OS genes closely associated with SKCM and illustrate their potential mechanisms. Transcriptome data and corresponding clinical traits of patients with SKCM were retrieved from The Cancer Genome Atlas and Gene Expression Omnibus databases. A weighted gene co-expression network analysis was conducted to identify relationships between clinical features and OS genes in specific modules. Subsequently, Cox regression analysis was performed on candidate OS genes; four hub prognosis-associated OS genes (AKAP9, VPS13C, ACSL4, and HMOX2) were identified to construct a prognostic model. After a series of bioinformatics analysis, our prognostic model was identified significantly associated with the overall survival of patients with SKCM and metastatic ability of the cancer. Furthermore, our risk model demonstrated improved diagnostic accuracy in the Cancer Genome Atlas and Gene Expression Omnibus cohorts. In addition, we established 2 nomograms based on either risk score or hub genes, which displayed favorable discriminating ability for SKCM. Our results provide novel insight into the potential applications of OS-associated genes in SKCM.

Do Molecular Profiles of Primary Versus Metastatic Radioiodine Refractory Differentiated Thyroid Cancer Differ?

Management of metastatic radioiodine refractory differentiated thyroid cancer (DTC) can be a therapeutic challenge. Generally, little is known about the paired molecular profile of the primary tumor and the metastases and whether they harbor the same genetic abnormalities. The present study compared the molecular profile of paired tumor specimens (primary tumor/metastatic sites) from patients with radioiodine refractory DTC in order to gain insight into a possible basis for resistance to radioiodine. Twelve patients with radioiodine refractory metastases were studied; median age at diagnosis of 61 years (range, 25-82). Nine patients had papillary TC (PTC), one had follicular TC (FTC), and two had Hürthle cell TC (HTC). Distant metastases were present in the lungs (n = 10), bones (n = 4), and liver (n = 1). The molecular profiling of paired tumors was performed with a panel of 592 genes for Next Generation Sequencing, RNA-sequencing, and immunohistochemistry. Digital microfluidic PCR was used to investigate TERT promoter mutations. The genetic landscape of all paired sites comprised BRAF, NRAS, HRAS, TP53, ATM, MUTYH, POLE, and NTRK genes, including BRAF and NTRK fusions. BRAF V600E was the most common point mutation in the paired specimens (5/12). TERT promoter mutation C228T was detected in one case. PD-L1 expression at metastatic sites was highly positive (95%) for one patient with HTC. All specimens were stable for microsatellite instability testing, and the tumor mutation burden was low to intermediate. Therefore, the molecular profile of DTC primary and metastatic lesions can show heterogeneity, which may help explain some altered responses to therapeutic intervention.

Transcript levels of keratin 1/5/6/14/15/16/17 as potential prognostic indicators in melanoma patients

Keratins (KRTs), the intermediate filament-forming proteins of epithelial cells, are extensively used as diagnostic biomarkers in cancers and associated with tumorigenesis and metastasis in multiple cancers. However, the diverse expression patterns and prognostic values of KRTs in melanoma have yet to be elucidated. In the current study, we examined the transcriptional and clinical data of KRTs in patients with melanoma from GEO, TCGA, ONCOMINE, GEPIA, cBioPortal, TIMER and TISIDB databases. We found that the mRNA levels of KRT1/2/5/6/8/10/14/15/16/17 were significantly differential expressed between primary melanoma and metastatic melanoma. The expression levels of KRT1/2/5/6/10/14/15/16/17 were correlated with advanced tumor stage. Survival analysis revealed that the high transcription levels of KRT1/5/6/14/15/16/17 were associated with low overall survival in melanoma patients. GSEA analysis indicated that the most involved hallmarks pathways were P53 pathway, KRAS signaling, estrogen response early and estrogen response late. Furthermore, we found some correlations among the expression of KRTs and the infiltration of immune cells. Our study may provide novel insights for the selection of prognostic biomarkers for melanoma.