SOX10 promotes melanoma cell invasion by regulating melanoma inhibitory activity

Specific SOX10 enhancer elements modulate phenotype plasticity and drug resistance in melanoma

Recent studies indicate that the development of drug resistance and increased invasiveness in melanoma is largely driven by transcriptional plasticity rather than canonical coding mutations. Understanding the mechanisms behind cell identity shifts in oncogenic transformation and cancer progression is crucial for advancing our understanding of melanoma and other aggressive cancers. While distinct melanoma phenotypic states have been well characterized, the processes and transcriptional controls that enable cells to shift between these states remain largely unknown. In this study, we initially leverage the well-established zebrafish melanoma model as a high-throughput system to dissect and analyze transcriptional control elements that are hijacked by melanoma. We identify key characteristics of these elements, making them translatable to human enhancer identification despite the lack of direct sequence conservation. Building on our identification of a zebrafish sox10 enhancer necessary for melanoma initiation, we extend these findings to human melanoma, identifying two human upstream enhancer elements that are critical for full SOX10 expression. Stable biallelic deletion of these enhancers using CRISPR-Cas9 induces a distinct phenotype shift across multiple human melanoma cell lines from a melanocytic phenotype towards an undifferentiated phenotype and is also characterized by an increase in drug resistance that mirrors clinical data including an upregulation of NTRK1, a tyrosine kinase, and potential therapeutic target. These results provide new insights into the transcriptional regulation of SOX10 in human melanoma and underscore the role of individual enhancer elements and potentially NTRK1 in driving melanoma phenotype plasticity and drug resistance. Our work lays the groundwork for future gene-based and combination kinase-inhibitor therapies targeting SOX10 regulation and NTRK1 as a potential avenue for enhancing the efficacy of current melanoma treatments.

DrugMap: A quantitative pan-cancer analysis of cysteine ligandability

Cysteine-focused chemical proteomic platforms have accelerated the clinical development of covalent inhibitors for a wide range of targets in cancer. However, how different oncogenic contexts influence cysteine targeting remains unknown. To address this question, we have developed "DrugMap," an atlas of cysteine ligandability compiled across 416 cancer cell lines. We unexpectedly find that cysteine ligandability varies across cancer cell lines, and we attribute this to differences in cellular redox states, protein conformational changes, and genetic mutations. Leveraging these findings, we identify actionable cysteines in NF-κB1 and SOX10 and develop corresponding covalent ligands that block the activity of these transcription factors. We demonstrate that the NF-κB1 probe blocks DNA binding, whereas the SOX10 ligand increases SOX10-SOX10 interactions and disrupts melanoma transcriptional signaling. Our findings reveal heterogeneity in cysteine ligandability across cancers, pinpoint cell-intrinsic features driving cysteine targeting, and illustrate the use of covalent probes to disrupt oncogenic transcription-factor activity.

Dormancy of cutaneous melanoma

Many cancer-related deaths including melanoma result from metastases that develop months or years after the initial cancer therapy. Even the most effective drugs and immune therapies rarely eradicate all tumor cells. Instead, they strongly reduce cancer burden, permitting dormant cancer cells to persist in niches, where they establish a cellular homeostasis with their host without causing clinical symptoms. Dormant cancers respond poorly to most drugs and therapies since they do not proliferate and hide in niches. It therefore remains a major challenge to develop novel therapies for dormant cancers. In this review we focus on the mechanisms regulating the initiation of cutaneous melanoma dormancy as well as those which are involved in reawakening of dormant cutaneous melanoma cells. In recent years the role of neutrophils and niche components in reawakening of melanoma cells came into focus and indicate possible future therapeutic applications. Sophisticated in vitro and in vivo melanoma dormancy models are needed to make progress in this field and are discussed.

RNA sequencing analysis reveals that missense mutation in SOX10 is associated with iris color phenotype in quail

To investigate the molecular mechanisms underlying the differences in iris color in quail, the transcriptome of iris tissue from black quail and Korean quail at day 10 of hatching was RNA sequenced in this study. The differentially expressed genes (DEGs) were screened, functionally annotated and enriched after the quality control and mapping of the raw data. RT-qPCR validation was performed using EIF2S3 as an internal reference gene. The screened SNPs were studied by bioinformatics analysis and iris color correlation analysis. The results showed that there were 425 upregulated genes and 364 downregulated genes in 789 DEGs. Gene Ontology (GO) enrichment analysis revealed that 139 DEGs were significantly enriched in 154 GO terms. The Kyoto Encyclopedia of Genes and Genomes enrichment results showed that the Notch signaling pathway, melanogenesis and tyrosine metabolism were associated with pigment synthesis (p < 0.05). The expression levels of the ASIP, MLPH, PMEL, TYR and SOX10 genes were significantly different in black quail iris and Korean quail iris, as verified by RT-qPCR. The SOX10 gene c.324G>C mutation, which caused the replacement of p.Glu108Asp, had a highly significant correlation with iris color in black quail and Korean quail, which may be one of the reasons for different in iris color between these two quail species.

SOX10 deficiency-mediated LAMB3 upregulation determines the invasiveness of MAPKi-resistant melanoma

Melanoma that develops adaptive resistance to MAPK inhibitors (MAPKi) through transcriptional reprograming-mediated phenotype switching is associated with enhanced metastatic potential, yet the underlying mechanism of this improved invasiveness has not been fully elucidated. In this study, we show that MAPKi-resistant melanoma cells are more motile and invasive than the parental cells. We further show that LAMB3, a β subunit of the extracellular matrix protein laminin-332 is upregulated in MAPKi-resistant melanoma cells and that the LAMB3-Integrin α3/α6 signaling mediates the motile and invasive phenotype of resistant cells. In addition, we demonstrate that SOX10 deficiency in MAPKi-resistant melanoma cells drives LAMB3 upregulation through TGF-β signaling. Transcriptome profiling and functional studies further reveal a FAK/MMPs axis mediates the pro-invasiveness effect of LAMB3. Using a mouse lung metastasis model, we demonstrate LAMB3 depletion inhibits the metastatic potential of MAPKi-resistant cells in vivo. In summary, this study identifies a SOX10low/TGF-β/LAMB3/FAK/MMPs signaling pathway that determines the migration and invasion properties of MAPKi-resistant melanoma cells and provide rationales for co-targeting LAMB3 to curb the metastasis of melanoma cells in targeted therapy.

Sox10-Deficient Drug-Resistant Melanoma Cells Are Refractory to Oncolytic RNA Viruses

Targeted therapy resistance frequently develops in melanoma due to intratumor heterogeneity and epigenetic reprogramming. This also typically induces cross-resistance to immunotherapies. Whether this includes additional modes of therapy has not been fully assessed. We show that co-treatments of MAPKi with VSV-based oncolytics do not function in a synergistic fashion; rather, the MAPKis block infection. Melanoma resistance to vemurafenib further perturbs the cells' ability to be infected by oncolytic viruses. Resistance to vemurafenib can be induced by the loss of SOX10, a common proliferative marker in melanoma. The loss of SOX10 promotes a cross-resistant state by further inhibiting viral infection and replication. Analysis of RNA-seq datasets revealed an upregulation of interferon-stimulated genes (ISGs) in SOX10 knockout populations and targeted therapy-resistant cells. Interestingly, the induction of ISGs appears to be independent of type I IFN production. Overall, our data suggest that the pathway mediating oncolytic resistance is due to the loss of SOX10 during acquired drug resistance in melanoma.

Synthetic Lethality Screening with Recursive Feature Machines

Synthetic lethality refers to a genetic interaction where the simultaneous perturbation of gene pairs leads to cell death. Synthetically lethal gene pairs (SL pairs) provide a potential avenue for selectively targeting cancer cells based on genetic vulnerabilities. The rise of large-scale gene perturbation screens such as the Cancer Dependency Map (DepMap) offers the opportunity to identify SL pairs automatically using machine learning. We build on a recently developed class of feature learning kernel machines known as Recursive Feature Machines (RFMs) to develop a pipeline for identifying SL pairs based on CRISPR viability data from DepMap. In particular, we first train RFMs to predict viability scores for a given CRISPR gene knockout from cell line embeddings consisting of gene expression and mutation features. After training, RFMs use a statistical operator known as average gradient outer product to provide weights for each feature indicating the importance of each feature in predicting cellular viability. We subsequently apply correlation-based filters to re-weight RFM feature importances and identify those features that are most indicative of low cellular viability. Our resulting pipeline is computationally efficient, taking under 3 minutes for analyzing all 17, 453 knockouts from DepMap for candidate SL pairs. We show that our pipeline more accurately recovers experimentally verified SL pairs than prior approaches. Moreover, our pipeline finds new candidate SL pairs, thereby opening novel avenues for identifying genetic vulnerabilities in cancer.

DrugMap: A quantitative pan-cancer analysis of cysteine ligandability

Cysteine-focused chemical proteomic platforms have accelerated the clinical development of covalent inhibitors of a wide-range of targets in cancer. However, how different oncogenic contexts influence cysteine targeting remains unknown. To address this question, we have developed DrugMap , an atlas of cysteine ligandability compiled across 416 cancer cell lines. We unexpectedly find that cysteine ligandability varies across cancer cell lines, and we attribute this to differences in cellular redox states, protein conformational changes, and genetic mutations. Leveraging these findings, we identify actionable cysteines in NFκB1 and SOX10 and develop corresponding covalent ligands that block the activity of these transcription factors. We demonstrate that the NFκB1 probe blocks DNA binding, whereas the SOX10 ligand increases SOX10-SOX10 interactions and disrupts melanoma transcriptional signaling. Our findings reveal heterogeneity in cysteine ligandability across cancers, pinpoint cell-intrinsic features driving cysteine targeting, and illustrate the use of covalent probes to disrupt oncogenic transcription factor activity.

Dynamical modeling of proliferative-invasive plasticity and IFNγ signaling in melanoma reveals mechanisms of PD-L1 expression heterogeneity

BACKGROUND

Phenotypic heterogeneity of melanoma cells contributes to drug tolerance, increased metastasis, and immune evasion in patients with progressive disease. Diverse mechanisms have been individually reported to shape extensive intra-tumor and inter-tumor phenotypic heterogeneity, such as IFNγ signaling and proliferative to invasive transition, but how their crosstalk impacts tumor progression remains largely elusive.

METHODS

Here, we integrate dynamical systems modeling with transcriptomic data analysis at bulk and single-cell levels to investigate underlying mechanisms behind phenotypic heterogeneity in melanoma and its impact on adaptation to targeted therapy and immune checkpoint inhibitors. We construct a minimal core regulatory network involving transcription factors implicated in this process and identify the multiple 'attractors' in the phenotypic landscape enabled by this network. Our model predictions about synergistic control of PD-L1 by IFNγ signaling and proliferative to invasive transition were validated experimentally in three melanoma cell lines-MALME3, SK-MEL-5 and A375.

RESULTS

We demonstrate that the emergent dynamics of our regulatory network comprising MITF, SOX10, SOX9, JUN and ZEB1 can recapitulate experimental observations about the co-existence of diverse phenotypes (proliferative, neural crest-like, invasive) and reversible cell-state transitions among them, including in response to targeted therapy and immune checkpoint inhibitors. These phenotypes have varied levels of PD-L1, driving heterogeneity in immunosuppression. This heterogeneity in PD-L1 can be aggravated by combinatorial dynamics of these regulators with IFNγ signaling. Our model predictions about changes in proliferative to invasive transition and PD-L1 levels as melanoma cells evade targeted therapy and immune checkpoint inhibitors were validated in multiple RNA-seq data sets from in vitro and in vivo experiments.

CONCLUSION

Our calibrated dynamical model offers a platform to test combinatorial therapies and provide rational avenues for the treatment of metastatic melanoma. This improved understanding of crosstalk among PD-L1 expression, proliferative to invasive transition and IFNγ signaling can be leveraged to improve the clinical management of therapy-resistant and metastatic melanoma.

Sox10 escalates vascular inflammation by mediating vascular smooth muscle cell transdifferentiation and pyroptosis in neointimal hyperplasia

Vascular smooth muscle cells (VSMCs) can transdifferentiate into macrophage-like cells in the context of sustained inflammatory injury, which drives vascular hyperplasia and atherosclerotic complications. Using single-cell RNA sequencing, we identify that macrophage-like VSMCs are the key cell population in mouse neointimal hyperplasia. Sex-determining region Y (SRY)-related HMG-box gene 10 (Sox10) upregulation is associated with macrophage-like VSMC accumulation and pyroptosis in vitro and in the neointimal hyperplasia of mice. Tumor necrosis factor α (TNF-α)-induced Sox10 lactylation in a phosphorylation-dependent manner by PI3K/AKT signaling drives transcriptional programs of VSMC transdifferentiation, contributing to pyroptosis. The regulator of G protein signaling 5 (RGS5) interacts with AKT and blocks PI3K/AKT signaling and Sox10 phosphorylation at S24. Sox10 silencing mitigates vascular inflammation and forestalls neointimal hyperplasia in RGS5 knockout mice. Collectively, this study shows that Sox10 is a regulator of vascular inflammation and a potential control point in inflammation-related vascular disease.

Dynamical modelling of proliferative-invasive plasticity and IFNγ signaling in melanoma reveals mechanisms of PD-L1 expression heterogeneity

Phenotypic heterogeneity of melanoma cells contributes to drug tolerance, increased metastasis, and immune evasion in patients with progressive disease. Diverse mechanisms have been individually reported to shape extensive intra- and inter-tumoral phenotypic heterogeneity, such as IFNγ signaling and proliferative to invasive transition, but how their crosstalk impacts tumor progression remains largely elusive. Here, we integrate dynamical systems modeling with transcriptomic data analysis at bulk and single-cell levels to investigate underlying mechanisms behind phenotypic heterogeneity in melanoma and its impact on adaptation to targeted therapy and immune checkpoint inhibitors. We construct a minimal core regulatory network involving transcription factors implicated in this process and identify the multiple "attractors" in the phenotypic landscape enabled by this network. Our model predictions about synergistic control of PD-L1 by IFNγ signaling and proliferative to invasive transition were validated experimentally in three melanoma cell lines - MALME3, SK-MEL-5 and A375. We demonstrate that the emergent dynamics of our regulatory network comprising MITF, SOX10, SOX9, JUN and ZEB1 can recapitulate experimental observations about the co-existence of diverse phenotypes (proliferative, neural crest-like, invasive) and reversible cell-state transitions among them, including in response to targeted therapy and immune checkpoint inhibitors. These phenotypes have varied levels of PD-L1, driving heterogeneity in immune-suppression. This heterogeneity in PD-L1 can be aggravated by combinatorial dynamics of these regulators with IFNγ signaling. Our model predictions about changes in proliferative to invasive transition and PD-L1 levels as melanoma cells evade targeted therapy and immune checkpoint inhibitors were validated in multiple data sets from in vitro and in vivo experiments. Our calibrated dynamical model offers a platform to test combinatorial therapies and provide rational avenues for the treatment of metastatic melanoma. This improved understanding of crosstalk among PD-L1 expression, proliferative to invasive transition and IFNγ signaling can be leveraged to improve the clinical management of therapy-resistant and metastatic melanoma.

SM22α Deletion Contributes to Neurocognitive Impairment in Mice through Modulating Vascular Smooth Muscle Cell Phenotypes

Considerable evidence now indicates that cognitive impairment is primarily a vascular disorder. The depletion of smooth muscle 22 alpha (SM22α) contributes to vascular smooth muscle cells (VSMCs) switching from contractile to synthetic and proinflammatory phenotypes in the context of inflammation. However, the role of VSMCs in the pathogenesis of cognitive impairment remains undetermined. Herein, we showed a possible link between VSMC phenotypic switching and neurodegenerative diseases via the integration of multi-omics data. SM22α knockout (Sm22α^-/-) mice exhibited obvious cognitive impairment and cerebral pathological changes, which were visibly ameliorated by the administration of AAV-SM22α. Finally, we confirmed that SM22α disruption promotes the expression of SRY-related HMG-box gene 10 (Sox10) in VSMCs, thereby aggravating the systemic vascular inflammatory response and ultimately leading to cognitive impairment in the brain. Therefore, this study supports the idea of VSMCs and SM22α as promising therapeutic targets in cognitive impairment to improve memory and cognitive decline.

Cooperative induction of receptor tyrosine kinases contributes to adaptive MAPK drug resistance in melanoma through the PI3K pathway

Vemurafenib-induced drug resistance in melanoma has been linked to receptor tyrosine kinase (RTK) upregulation. The MITF and SOX10 genes play roles as master regulators of melanocyte and melanoma development. Here, we aimed to explore the complex mechanisms behind the MITF/SOX10-controlled RTK-induced drug resistance in melanoma. To achieve this, we used a number of molecular techniques, including melanoma patient data from TCGA, vemurafenib-resistant melanoma cell lines, and knock-down studies. The melanoma cell lines were classified as proliferative or invasive based upon their MITF/AXL expression activity. We measured the change of expression activity for MITF/SOX10 and their receptor (AXL/ERBB3) and ligand (NRG1/GAS6) targets known to be involved in RTK-induced drug resistance after vemurafenib treatment. We find that melanoma cell lines characterized as proliferative (high MITF low AXL), transform into an invasive (low MITF, high AXL) cell state after vemurafenib resistance, indicating novel feedback loops and advanced compensatory regulation mechanisms between the master regulators, receptors, and ligands involved in vemurafenib-induced resistance. Together, our data disclose fine-tuned mechanisms involved in RTK-facilitated vemurafenib resistance that will be challenging to overcome by using single drug targeting strategies against melanoma.

A cis-regulatory lexicon of DNA motif combinations mediating cell-type-specific gene regulation

Gene expression is controlled by transcription factors (TFs) that bind cognate DNA motif sequences in cis-regulatory elements (CREs). The combinations of DNA motifs acting within homeostasis and disease, however, are unclear. Gene expression, chromatin accessibility, TF footprinting, and H3K27ac-dependent DNA looping data were generated and a random-forest-based model was applied to identify 7,531 cell-type-specific cis-regulatory modules (CRMs) across 15 diploid human cell types. A co-enrichment framework within CRMs nominated 838 cell-type-specific, recurrent heterotypic DNA motif combinations (DMCs), which were functionally validated using massively parallel reporter assays. Cancer cells engaged DMCs linked to neoplasia-enabling processes operative in normal cells while also activating new DMCs only seen in the neoplastic state. This integrative approach identifies cell-type-specific cis-regulatory combinatorial DNA motifs in diverse normal and diseased human cells and represents a general framework for deciphering cis-regulatory sequence logic in gene regulation.

Improved NSGA-II algorithms for multi-objective biomarker discovery

MOTIVATION

In modern translational research, the development of biomarkers heavily relies on use of omics technologies, but implementations with basic data mining algorithms frequently lead to false positives. Non-dominated Sorting Genetic Algorithm II (NSGA2) is an extremely effective algorithm for biomarker discovery but has been rarely evaluated against large-scale datasets. The exploration of the feature search space is the key to NSGA2 success but in specific cases NSGA2 expresses a shallow exploration of the space of possible feature combinations, possibly leading to models with low predictive performances.

RESULTS

We propose two improved NSGA2 algorithms for finding subsets of biomarkers exhibiting different trade-offs between accuracy and feature number. The performances are investigated on gene expression data of breast cancer patients. The results are compared with NSGA2 and LASSO. The benchmarking dataset includes internal and external validation sets. The results show that the proposed algorithms generate a better approximation of the optimal trade-offs between accuracy and set size. Moreover, validation and test accuracies are better than those provided by NSGA2 and LASSO. Remarkably, the GA-based methods provide biomarkers that achieve a very high prediction accuracy (>80%) with a small number of features (<10), representing a valid alternative to known biomarker models, such as Pam50 and MammaPrint.

AVAILABILITY AND IMPLEMENTATION

The software is publicly available on GitHub at github.com/UEFBiomedicalInformaticsLab/BIODAI/tree/main/MOO.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Recent advances in the regulation mechanism of SOX10

Neural crest (NC) is the primitive neural structure in embryonic stage, which develops from ectodermal neural plate cells and epithelial cells. When the neural fold forms into neural tube, neural crest also forms a cord like structure above the neural tube and below the ectoderm. Neural crest cells (NCC) have strong migration and proliferation abilities. A number of tissue cells differentiate from neural crest cells, such as melanocytes, central and peripheral neurons, glial cells, craniofacial cells, osteoblasts, chondrocytes and smooth muscle cells. The migration and differentiation of neural crest cells are regulated by a gene network where a variety of genes, transcriptional factors, signal pathways and growth factors are involved.

DEPDC1B Promotes Melanoma Angiogenesis and Metastasis through Sequestration of Ubiquitin Ligase CDC16 to Stabilize Secreted SCUBE3

The ability of melanoma to acquire metastasis through the induction of angiogenesis is one of the major causes of skin cancer death. Here, it is found that high transcript levels of DEP domain containing 1B (DEPDC1B) in cutaneous melanomas are significantly associated with a poor prognosis. Tissue microarray analysis indicates that DEPDC1B expression is positively correlated with SOX10 in the different stages of melanoma. Consistently, DEPDC1B is both required and sufficient for melanoma growth, metastasis, angiogenesis, and functions as a direct downstream target of SOX10 to partly mediate its oncogenic activity. In contrast to other tumor types, the DEPDC1B-mediated enhancement of melanoma metastatic potential is not dependent on the activities of RHO GTPase signaling and canonical Wnt signaling, but is acquired through secretion of signal peptide, CUB domain and EGF like domain containing 3 (SCUBE3), which is crucial for promoting angiogenesis in vitro and in vivo. Mechanistically, DEPDC1B regulates SCUBE3 protein stability through the competitive association with ubiquitin ligase cell division cycle 16 (CDC16) to prevent SCUBE3 from undergoing degradation via the ubiquitin-proteasome pathway. Importantly, expression of SOX10, DEPDC1B, and SCUBE3 are positively correlated with microvessel density in the advanced stage of melanomas. In conclusion, it is revealed that a SOX10-DEPDC1B-SCUBE3 regulatory axis promotes melanoma angiogenesis and metastasis, which suggests that targeting secreted SCUBE3 can be a therapeutic strategy against metastatic melanoma.

Targeting SOX10-deficient cells to reduce the dormant-invasive phenotype state in melanoma

Cellular plasticity contributes to intra-tumoral heterogeneity and phenotype switching, which enable adaptation to metastatic microenvironments and resistance to therapies. Mechanisms underlying tumor cell plasticity remain poorly understood. SOX10, a neural crest lineage transcription factor, is heterogeneously expressed in melanomas. Loss of SOX10 reduces proliferation, leads to invasive properties, including the expression of mesenchymal genes and extracellular matrix, and promotes tolerance to BRAF and/or MEK inhibitors. We identify the class of cellular inhibitor of apoptosis protein-1/2 (cIAP1/2) inhibitors as inducing cell death selectively in SOX10-deficient cells. Targeted therapy selects for SOX10 knockout cells underscoring their drug tolerant properties. Combining cIAP1/2 inhibitor with BRAF/MEK inhibitors delays the onset of acquired resistance in melanomas in vivo. These data suggest that SOX10 mediates phenotypic switching in cutaneous melanoma to produce a targeted inhibitor tolerant state that is likely a prelude to the acquisition of resistance. Furthermore, we provide a therapeutic strategy to selectively eliminate SOX10-deficient cells.

Transcriptome analysis of genes potentially associated with white and black plumage formation in Chinese indigenous ducks (Anas platyrhynchos)

1. Plumage colour is an important recognisable characteristic of duck (Anas platyrhynchos), but the coloration mechanisms remain largely unknown. To elucidate the molecular mechanisms underlying the formation of black and white plumage, the following study applied RNA sequencing (RNA-Seq) to catalogue the global gene expression profiles in the duck feather bulbs of black and white colours.2. Black feather bulbs were collected from Putian Black ducks (B-PTB) and black Longsheng Jade-green ducks (B-LS), while white feather bulbs were collected from Putian White ducks (W-PTW), Putian Black ducks (W-PTB) and Longsheng Jade-green ducks (W-LS). Sixteen cDNA libraries were constructed and sequenced for transcriptional analysis. Three comparison groups were employed to analyse differentially expressed genes (DEGs), including W-PTB versus B-PTB, W-PTW versus B-PTB and W-LS versus B-LS.3. The results showed 180 DEGs between W-PTB and B-PTB, 303 DEGs between W-PTW and B-PTB, and 108 DEGs between W-LS and B-LS. Further analysis showed that 18 DEGs were directly involved in the pigmentation process and melanogenesis signalling pathway. Additionally, the distribution of DEGs varied amongst groups whereby ASIP appeared only in the W-LS versus B-LS group, GNAI1 and ZEB2 appeared only in the W-PTW versus B-PTB group, and KITLG, EDN3 and FZD4 appeared only in W-PTB versus B-PTB.4. The findings suggested that the mechanism of feather albinism may differ between duck breeds. This study provided new information for discovering genes that are important for feather pigmentation and helps elucidate molecular mechanisms involved in black and white plumage in ducks.

In Vivo Melanoma Cell Morphology Reflects Molecular Signature and Tumor Aggressiveness

Melanoma is the deadliest type of skin cancer, characterized by high cellular heterogeneity which contributes to therapy resistance and unpredictable disease outcome. Recently, by correlating Reflectance-Confocal-Microscopy (RCM) morphology with histopathological type, we identified four distinct melanoma-subtypes: dendritic-cell (DC), round-cell (RC), dermal-nest (DN), and combined-type (CT) melanomas. In the present study, each RCM-melanoma subtype expressed a specific biomolecular profile and biological behavior in vitro. Markers of tumor aggressiveness, including Ki67, MERTK, nestin and stemness markers, were highest in the most invasive CT and DN melanomas, as compared to DC and RC. This was also confirmed in multicellular tumor spheroids. Transcriptomic analysis showed a modulation of cancer progression-associated genes from DC to CT melanomas. The switch from E- to N-cadherin expression proved the epithelial-to-mesenchymal transition from DC to CT subtypes. The DN melanoma was predominantly located in the dermis, as also shown in skin reconstructs. It displayed a unique behavior and a molecular profile associated with a high degree of aggressiveness. Altogether, our results demonstrate that each RCM-melanoma subtype has a distinct biological and gene expression profile, related to tumor aggressiveness, confirming that RCM can be a dependable tool for in vivo detecting different types of melanoma and for early diagnostic screening.

Gene Expression and Mutational Profile in BAP-1 Inactivated Melanocytic Lesions of Progressive Malignancy from a Patient with Multiple Lesions

BAP-1 (BRCA1-associated protein 1) inactivated melanocytic lesions are a group of familial or sporadic lesions with unique histology and molecular features. They are of great clinical interest, at least in part due to the potential for malignant transformation and association with a familial cancer predisposition syndrome. Here, we describe a patient with multiple spatially and temporally distinct melanocytic lesions with loss of BAP1 expression by immunohistochemistry. RNA sequencing was performed on three independent lesions spanning the morphologic spectrum: a benign nevus, an atypical tumor, and a melanoma arising from a pre-existing BAP1-inactivated nevus. The three lesions demonstrated largely distinct gene expression and mutational profiles. Gene expression analysis revealed that genes involved in receptor protein kinase pathways were progressively upregulated from nevus to melanoma. Moreover, a clear enrichment of genes regulated in response to UV radiation was found in the melanoma from this patient, as well as upregulation of MAPK pathway-related genes and several transcription factors related to melanomagenesis.

Signal pathways of melanoma and targeted therapy

Melanoma is the most lethal skin cancer that originates from the malignant transformation of melanocytes. Although melanoma has long been regarded as a cancerous malignancy with few therapeutic options, increased biological understanding and unprecedented innovations in therapies targeting mutated driver genes and immune checkpoints have substantially improved the prognosis of patients. However, the low response rate and inevitable occurrence of resistance to currently available targeted therapies have posed the obstacle in the path of melanoma management to obtain further amelioration. Therefore, it is necessary to understand the mechanisms underlying melanoma pathogenesis more comprehensively, which might lead to more substantial progress in therapeutic approaches and expand clinical options for melanoma therapy. In this review, we firstly make a brief introduction to melanoma epidemiology, clinical subtypes, risk factors, and current therapies. Then, the signal pathways orchestrating melanoma pathogenesis, including genetic mutations, key transcriptional regulators, epigenetic dysregulations, metabolic reprogramming, crucial metastasis-related signals, tumor-promoting inflammatory pathways, and pro-angiogenic factors, have been systemically reviewed and discussed. Subsequently, we outline current progresses in therapies targeting mutated driver genes and immune checkpoints, as well as the mechanisms underlying the treatment resistance. Finally, the prospects and challenges in the development of melanoma therapy, especially immunotherapy and related ongoing clinical trials, are summarized and discussed.

SOX10 regulates melanoma immunogenicity through an IRF4-IRF1 axis

Loss-of-function mutations of JAK1/2 impair cancer cell responsiveness to IFN-γ and immunogenicity. Therefore, an understanding of compensatory pathways to activate IFN-γ signaling in cancer cells is clinically important for the success of immunotherapy. Here we demonstrate that the transcription factor SOX10 hinders immunogenicity of melanoma cells through the IRF4-IRF1 axis. Genetic and pharmacological approaches revealed that SOX10 repressed IRF1 transcription via direct induction of a negative regulator, IRF4. The SOX10-IRF4-IRF1 axis regulated PD-L1 expression independently of JAK-STAT pathway activity, and suppression of SOX10 increased the efficacy of combination therapy with an anti-PD-1 antibody and HDAC inhibitor against a clinically relevant melanoma model. Thus, the SOX10-IRF4-IRF1 axis serves as a potential target that can bypass JAK-STAT signaling to immunologically warm up melanoma with a "cold" tumor immune microenvironment.

Therapeutic combination silencing VEGF and SOX10 increases the antiangiogenic effect in the mouse melanoma model B16-F10 - in vitro and in vivo studies

INTRODUCTION

Gene therapy is an innovative form of treatment of genetic diseases, in which psiRNA molecules silencing specific genes are applied.

AIM

The study evaluated the anti-tumour effect of psiRNA silencing preparations of the vascular endothelial growth factor (VEGF) and Sry-related HMG-Box gene 10 (SOX10) on melanoma (B16-F10) by inhibiting angiogenesis.

MATERIAL AND METHODS

The preparations based on plasmid vectors psiRNA silencing the gene SOX10 and VEGF that form complexes with cationic lipid (psiRNA/carrier) have been developed. psiRNA preparations were tested on the mouse melanoma cell line B16-F10, both in vitro and in vivo. The silencing activity of transfected melanoma cells with the obtained psiRNA preparations was examined using the qPCR and Western blot methods. The anti-tumour activity of psiRNA preparations on melanoma tumour cells was then evaluated in a mouse in vivo model.

RESULTS

In vitro studies have shown that the B16-F10 cells efficiently transfect non-viral preparations - psiRNA: Lyovec (74-89%). Worth mentioning is the fact that silencing SOX10 in B16-F10 melanoma cells increases the expression of the COL18A1 gene (compared to the preparation inhibiting only VEGF), which codes the endostatin to stop angiogenesis. In vivo results show that the level of haemoglobin in tumours of mice treated with psiRNA formulations was over 6 times lower than controls and tumour mass was 60-80% lower.

CONCLUSIONS

The novel study proves that simultaneous inhibition of SOX10 and VEGF enhances the antiangiogenic action and thus contributes to a significant halt of disease development. In addition, these data expand knowledge about SOX10 regulation and functions.

SOX10 Knockdown Inhibits Melanoma Cell Proliferation via Notch Signaling Pathway

Purpose

Melanoma is a serious and malignant disease worldwide. Seeking diagnostic markers and potential therapeutic targets is urgent for melanoma treatment. SOX10, a member of the SoxE family of genes, is a transcription factor which can regulate the transcription of a wide variety of genes in multiple cellular processes.

Methods

The mRNA level and protein expression of SOX10 is confirmed by bioinformatic analysis and IHC staining. MTT, clone formation and EdU analysis showed that SOX10 knockdown (KD) could significantly inhibit melanoma cell proliferation. FACS analysis showed that SOX10 KD could markedly enhance the level of cell apoptosis. The downstream target signaling pathway is predicted by RNA-seq based on the public GEO database. The activation of Notch signaling mediated by SOX10 is tested by qPCR and Western blot.

Results

Ectopic upregulation of SOX10 was found in melanoma patient tissues compared to normal nevus tissues in mRNA and protein levels. Furthermore, both mRNA and protein level of SOX10 were negatively correlated with melanoma patient's prognosis. SOX10 knockdown could obviously suppress the proliferation ability of melanoma cells by inactivating Notch signaling pathway.

Conclusion

Our study confirmed that SOX10 is an oncogene and activate Notch signaling pathway, which suggests the potential treatment for melanoma patients by target SOX10/Notch axis.

PITX1 inhibits the growth and proliferation of melanoma cells through regulation of SOX family genes

Melanoma is one of the most aggressive types of cancer wherein resistance to treatment prevails. Therefore, it is important to discover novel molecular targets of melanoma progression as potential treatments. Here we show that paired-like homeodomain transcription factor 1 (PITX1) plays a crucial role in the inhibition of melanoma progression through regulation of SRY-box transcription factors (SOX) gene family mRNA transcription. Overexpression of PITX1 in melanoma cell lines resulted in a reduction in cell proliferation and an increase in apoptosis. Additionally, analysis of protein levels revealed an antagonistic cross-regulation between SOX9 and SOX10. Interestingly, PITX1 binds to the SOX9 promoter region as a positive regulatory transcription factor; PITX1 mRNA expression levels were positively correlated with SOX9 expression, and negatively correlated with SOX10 expression in melanoma tissues. Furthermore, transcription of the long noncoding RNA (lncRNA), survival-associated mitochondrial melanoma-specific oncogenic noncoding RNA (SAMMSON), was decreased in PITX1-overexpressing cells. Taken together, the findings in this study indicate that PITX1 may act as a negative regulatory factor in the development and progression of melanoma via direct targeting of the SOX signaling.

SFPQ promotes an oncogenic transcriptomic state in melanoma

The multifunctional protein, splicing factor, proline- and glutamine-rich (SFPQ) has been implicated in numerous cancers often due to interaction with coding and non-coding RNAs, however, its role in melanoma remains unclear. We report that knockdown of SFPQ expression in melanoma cells decelerates several cancer-associated cell phenotypes, including cell growth, migration, epithelial to mesenchymal transition, apoptosis, and glycolysis. RIP-seq analysis revealed that the SFPQ-RNA interactome is reprogrammed in melanoma cells and specifically enriched with key melanoma-associated coding and long non-coding transcripts, including SOX10, AMIGO2 and LINC00511 and in most cases SFPQ is required for the efficient expression of these genes. Functional analysis of two SFPQ-enriched lncRNA, LINC00511 and LINC01234, demonstrated that these genes independently contribute to the melanoma phenotype and a more detailed analysis of LINC00511 indicated that this occurs in part via modulation of the miR-625-5p/PKM2 axis. Importantly, analysis of a large clinical cohort revealed that elevated expression of SFPQ in primary melanoma tumours may have utility as a prognostic biomarker. Together, these data suggest that SFPQ is an important driver of melanoma, likely due to SFPQ-RNA interactions promoting the expression of numerous oncogenic transcripts.

How Neural Crest Transcription Factors Contribute to Melanoma Heterogeneity, Cellular Plasticity, and Treatment Resistance

Cutaneous melanoma represents one of the deadliest types of skin cancer. The prognosis strongly depends on the disease stage, thus early detection is crucial. New therapies, including BRAF and MEK inhibitors and immunotherapies, have significantly improved the survival of patients in the last decade. However, intrinsic and acquired resistance is still a challenge. In this review, we discuss two major aspects that contribute to the aggressiveness of melanoma, namely, the embryonic origin of melanocytes and melanoma cells and cellular plasticity. First, we summarize the physiological function of epidermal melanocytes and their development from precursor cells that originate from the neural crest (NC). Next, we discuss the concepts of intratumoral heterogeneity, cellular plasticity, and phenotype switching that enable melanoma to adapt to changes in the tumor microenvironment and promote disease progression and drug resistance. Finally, we further dissect the connection of these two aspects by focusing on the transcriptional regulators MSX1, MITF, SOX10, PAX3, and FOXD3. These factors play a key role in NC initiation, NC cell migration, and melanocyte formation, and we discuss how they contribute to cellular plasticity and drug resistance in melanoma.

Epithelial-to-mesenchymal-like transition events in melanoma

Epithelial-to-mesenchymal transition (EMT), a process through which epithelial tumor cells acquire mesenchymal phenotypic properties, contributes to both metastatic dissemination and therapy resistance in cancer. Accumulating evidence indicates that nonepithelial tumors, including melanoma, can also gain mesenchymal-like properties that increase their metastatic propensity and decrease their sensitivity to therapy. In this review, we discuss recent findings, illustrating the striking similarities-but also knowledge gaps-between the biology of mesenchymal-like state(s) in melanoma and mesenchymal state(s) from epithelial cancers. Based on this comparative analysis, we suggest hypothesis-driven experimental approaches to further deepen our understanding of the EMT-like process in melanoma and how such investigations may pave the way towards the identification of clinically relevant biomarkers for prognosis and new therapeutic strategies.

LACTB suppresses melanoma progression by attenuating PP1A and YAP interaction

Very limited progress has been made in the management of advanced melanoma, especially melanoma of uveal origin. Lactamase β (LACTB) is a novel tumor suppressor; however, its biological function in melanoma remains unknown. Herein we demonstrated markedly lower LACTB expression levels in melanoma tissues and cell lines. Overexpression of LACTB suppressed the proliferation, migration and invasion of melanoma cells in vitro. Mechanistically, LACTB inhibited the activity of yes-associated protein (YAP). We showed that the level of phospho-YAP (Serine 127) was increased upon LACTB overexpression, which prevented the translocation of YAP to the nucleus. Further, LACTB could directly bind to PP1A and attenuate the interaction between PP1A and YAP, resulting in decreased YAP dephosphorylation and inactivation in a LATS1-independent manner. Additionally, transfection of phosphorylation-defective YAP mutants reversed LACTB-induced tumor suppression. Upstream, we demonstrated that SOX10 binds to the LACTB promoter and negatively regulates its transcription. Overexpression of LACTB also suppressed the tumorigenicity and lung metastasis of MUM2B uveal melanoma cells in vivo. Taken together, our findings indicate a novel SOX10/LACTB/PP1A signaling cascade that renders YAP inactive and modulates melanoma progression, offering a new therapeutic target for melanoma treatment.

Sex-Determining Region Y Chromosome-Related High-Mobility-Group Box 10 in Cancer: A Potential Therapeutic Target

Sex-determining region Y-related high mobility group-box 10 (SOX10), a member of the SOX family, has recently been highlighted as an essential transcriptional factor involved in developmental biology. Recently, the functionality of SOX 10 has been increasingly revealed by researchers worldwide. It has been reported that SOX10 significantly regulates the proliferation, migration, and apoptosis of tumors and is closely associated with the progression of cancer. In this review, we first introduce the basic background of the SOX family and SOX10 and then discuss the pathophysiological roles of SOX10 in cancer. Besides, we enumerate the application of SOX10 in the pathological diagnosis and therapeutic potential of cancer. Eventually, we summarize the potential directions and perspectives of SOX10 in neoplastic theranostics. The information compiled herein may assist in additional studies and increase the potential of SOX10 as a therapeutic target for cancer.

Cross-species analysis of enhancer logic using deep learning

Deciphering the genomic regulatory code of enhancers is a key challenge in biology because this code underlies cellular identity. A better understanding of how enhancers work will improve the interpretation of noncoding genome variation and empower the generation of cell type-specific drivers for gene therapy. Here, we explore the combination of deep learning and cross-species chromatin accessibility profiling to build explainable enhancer models. We apply this strategy to decipher the enhancer code in melanoma, a relevant case study owing to the presence of distinct melanoma cell states. We trained and validated a deep learning model, called DeepMEL, using chromatin accessibility data of 26 melanoma samples across six different species. We show the accuracy of DeepMEL predictions on the CAGI5 challenge, where it significantly outperforms existing models on the melanoma enhancer of IRF4 Next, we exploit DeepMEL to analyze enhancer architectures and identify accurate transcription factor binding sites for the core regulatory complexes in the two different melanoma states, with distinct roles for each transcription factor, in terms of nucleosome displacement or enhancer activation. Finally, DeepMEL identifies orthologous enhancers across distantly related species, where sequence alignment fails, and the model highlights specific nucleotide substitutions that underlie enhancer turnover. DeepMEL can be used from the Kipoi database to predict and optimize candidate enhancers and to prioritize enhancer mutations. In addition, our computational strategy can be applied to other cancer or normal cell types.

Transcriptional characterization of conjunctival melanoma identifies the cellular tumor microenvironment and prognostic gene signatures

This study characterizes the transcriptome and the cellular tumor microenvironment (TME) of conjunctival melanoma (CM) and identifies prognostically relevant biomarkers. 12 formalin-fixed and paraffin-embedded CM were analyzed by MACE RNA sequencing, including six cases each with good or poor clinical outcome, the latter being defined by local recurrence and/or systemic metastases. Eight healthy conjunctival specimens served as controls. The TME of CM, as determined by bioinformatic cell type enrichment analysis, was characterized by the enrichment of melanocytes, pericytes and especially various immune cell types, such as plasmacytoid dendritic cells, natural killer T cells, B cells and mast cells. Differentially expressed genes between CM and control were mainly involved in inhibition of apoptosis, proteolysis and response to growth factors. POU3F3, BIRC5 and 7 were among the top expressed genes associated with inhibition of apoptosis. 20 genes, among them CENPK, INHA, USP33, CASP3, SNORA73B, AAR2, SNRNP48 and GPN1, were identified as prognostically relevant factors reaching high classification accuracy (area under the curve: 1.0). The present study provides new insights into the TME and the transcriptional profile of CM and additionally identifies new prognostic biomarkers. These results add new diagnostic tools and may lead to new options of targeted therapy for CM.

Effect of Konjac Glucomannan (KGM) on the Reconstitution of the Dermal Environment against UVB-Induced Condition

Skin layers serve as a barrier against unexpected critical changes in the body due to environmental factors. Excessive ultraviolet (UV) B exposure increases the levels of age-related factors, leading to senescent cells and damaged skin tissues. Widely used as a dietary supplement, konjac (Amorphophallus konjac) glucomannan (KGM) has shown skin regeneration potential in patch or sheet form with anti-inflammatory or immunosuppressive effects. However, the ability of KGM to reconstitute senescent/damaged skin following UV radiation has not been explored. Here, we demonstrate that KGM alleviates skin damage by increasing the proportion of young cell populations in UVB-exposed senescent human epidermal primary melanocytes. Young cell numbers increased depending on KGM dosage, but the senescent cells were not removed. Real-time quantitative polymerase chain reaction (RT-qPCR) and Western blot analysis showed that mRNA and protein levels of age- and pigmentation-related factors decreased in a manner dependent on the rate at which new cells were generated. Moreover, an analysis of mRNA and protein levels indicated that KGM facilitated youth by increasing cell proliferation in UVB-damaged human fibroblasts. Thus, KGM is a highly effective natural agent for maintaining skin homeostasis by promoting the reconstitution of the dermal environment against UVB-induced acute senescence or skin damage.

HOXA9 mediates and marks premalignant compartment size expansion in colonic adenomas

The transformation of normal colonic epithelium to colorectal cancer (CRC) involves a relatively ordered progression, and understanding the molecular alterations involved may aid rational design of strategies aimed at preventing or counteracting disease. Homeobox A9 (HOXA9) is an oncogene in leukemia and has been implicated in CRC pathology, although its role in disease etiology remains obscure at best. We observe that HOXA9 expression is increased in colonic adenomas compared with location-matched healthy colon epithelium. Its forced expression results in dramatic genetic and signaling changes, with increased expression of growth factors IGF1 and FLT3, super-activity of the AKT survival pathway and a concomitant increase in compartment size. Furthermore, a reduced mRNA expression of the epithelial to mesenchymal transition marker N-cadherin as well as reduced activity of the actin cytoskeletal mediator PAK was seen, which is in apparent agreement with an observed reduced migratory response in HOXA9-overexpressing cells. Thus, HOXA9 appears closely linked with adenoma growth while impairing migration and metastasis and hence is both a marker and driver of premalignant polyp growth. Colonic polyps grow but remain premalignant for up to decades. Here, we show that HOXA9 drives growth in premalignant polyps, but simultaneously prevents further transformation.

Nuclear DLC1 exerts oncogenic function through association with FOXK1 for cooperative activation of MMP9 expression in melanoma

A Rho GTPase-activating protein (RhoGAP), deleted in liver cancer 1 (DLC1), is known to function as a tumor suppressor in various cancer types; however, whether DLC1 is a tumor-suppressor gene or an oncogene in melanoma remains to be clarified. Here we revealed that high DLC1 expression was detected in most of the melanoma tissues where it was localized in both the nuclei and the cytoplasm. Functional studies unveiled that DLC1 was both required and sufficient for melanoma growth and metastasis. These tumorigenic events were mediated by nuclear-localized DLC1 in a RhoGAP-independent manner. Mechanistically, mass spectrometry analysis identified a DLC1-associated protein, FOXK1 transcription factor, which mediated oncogenic events in melanoma by translocating and retaining DLC1 into the nucleus. RNA-sequencing profiling studies further revealed MMP9 as a direct target of FOXK1 through DLC1-regulated promoter occupancy for cooperative activation of MMP9 expression to promote melanoma invasion and metastasis. Concerted action of DLC1-FOXK1 in MMP9 gene regulation was further supported by their highly correlated expression in melanoma patients' samples and cell lines. Together, our results not only unravel a mechanism by which nuclear DLC1 functions as an oncogene in melanoma but also suggest an unexpected role of RhoGAP protein in transcriptional regulation.

A Review of the Molecular Pathways Involved in Resistance to BRAF Inhibitors in Patients with Advanced-Stage Melanoma

Melanoma is an aggressive malignancy of melanocytes and most commonly arises in the skin. In 2002, BRAF gene mutations were identified in melanoma, and this finding resulted in the development of several small-molecule molecular inhibitors that specifically targeted the BRAF V600E mutation. The development of targeted therapies for advanced-stage melanoma, including tyrosine kinase inhibitors (TKIs) of the BRAF (V600E) kinase, vemurafenib and dabrafenib, have been approved for the treatment of advanced melanoma leading to improved clinical outcomes. However, the development of BRAF inhibitor (BRAFi) resistance has significantly reduced the therapeutic efficacy after prolonged treatment. Recent studies have identified the molecular mechanisms for BRAFi resistance. This review aims to describe the impact of BRAFi resistance on the pathogenesis of melanoma, the current status of molecular pathways involved in BRAFi resistance, including intrinsic resistance, adaptive resistance, and acquired resistance. This review will discuss how an understanding of the mechanisms associated with BRAFi resistance may aid the identification of useful strategies for overcoming the resistance to BRAF-targeted therapy in patients with advanced-stage melanoma.

Cardiac myxoma and neural crests: a tense relationship

In cardiac myxomas, the malignant transformation process, selecting incidental gene mutations and leading to loss of proliferation control, has not a so drastic effects in terms of growth rate of tumor mass, but frequently the particular location of lesion engrosses the high risk for health. For accurate cancer cell profiling, it is important to establish the embryologic origin of malignant cells and their initial commitments, above all, in the sight of therapeutic strategies and solutions. Here, we advance, for cardiac myxoma, the hypothesis of an origin from cardiac neural crest cells and we attempt to support it by an integrated discussion of current knowledge about embryological characteristics of neural crest cells and most recent studies focusing cardiac myxomas. We discuss the relationship between the basic plasticity of cardiac neural crest cells and some typical mutations arising in neoplastic lesions as well as the expression of typical cell markers of neural crests derivatives. Dysfunctions in proliferative and migratory programs, focused in other studies, are evaluated in the context of the topological and histopathological characteristics of cardiac myxomas.

Complex Formation with Monomeric α-Tubulin and Importin 13 Fosters c-Jun Protein Stability and Is Required for c-Jun's Nuclear Translocation and Activity

Microtubules are highly dynamic structures, which consist of α- and β-tubulin heterodimers. They are essential for a number of cellular processes, including intracellular trafficking and mitosis. Tubulin-binding chemotherapeutics are used to treat different types of tumors, including malignant melanoma. The transcription factor c-Jun is a central driver of melanoma development and progression. Here, we identify the microtubule network as a main regulator of c-Jun activity. Monomeric α-tubulin fosters c-Jun protein stability by protein-protein interaction. In addition, this complex formation is necessary for c-Jun's nuclear localization sequence binding to importin 13, and consequent nuclear import and activity of c-Jun. A reduction in monomeric α-tubulin levels by treatment with the chemotherapeutic paclitaxel resulted in a decline in the nuclear accumulation of c-Jun in melanoma cells in an experimental murine model and in patients' tissues. These findings add important knowledge to the mechanism of the action of microtubule-targeting drugs and indicate the newly discovered regulation of c-Jun by the microtubule cytoskeleton as a novel therapeutic target for melanoma and potentially also other types of cancer.

Crosstalk of intracellular post-translational modifications in cancer

Post-translational modifications (PTMs) have been reported to play pivotal roles in numerous cellular biochemical and physiological processes. Multiple PTMs can influence the actions of each other positively or negatively, termed as PTM crosstalk or PTM code. During recent years, development of identification strategies for PTMs co-occurrence has revealed abundant information of interplay between PTMs. Increasing evidence demonstrates that deregulation of PTMs crosstalk is involved in the genesis and development of various diseases. Insight into the complexity of PTMs crosstalk will help us better understand etiology and provide novel targets for drug therapy. In the present review, we will discuss the important functional roles of PTMs crosstalk in proteins associated with cancer diseases.

SOX9 is a dose-dependent metastatic fate determinant in melanoma

BACKGROUND

In this research, we aimed to resolve contradictory results whether SOX9 plays a positive or negative role in melanoma progression and determine whether SOX9 and its closely related member SOX10 share the same or distinct targets in mediating their functions in melanoma.

METHODS

Immunofluorescence, TCGA database and qPCR were used to analyze the correlation between the expression patterns and levels of SOX9, SOX10 and NEDD9 in melanoma patient samples. AlamarBlue, transwell invasion and colony formation assays in melanoma cell lines were conducted to investigate the epistatic relationship between SOX10 and NEDD9, as well as the effects of graded SOX9 expression levels. Lung metastasis was determined by tail vein injection assay. Live cell imaging was conducted to monitor dynamics of melanoma migratory behavior. RHOA and RAC1 activation assays measured the activity of Rho GTPases.

RESULTS

High SOX9 expression was predominantly detected in patients with distant melanoma metastases whereas SOX10 was present in the different stages of melanoma. Both SOX9 and SOX10 exhibited distinct but overlapping expression patterns with metastatic marker NEDD9. Accordingly, SOX10 was required for NEDD9 expression, which partly mediated its oncogenic functions in melanoma cells. Compensatory upregulation of SOX9 expression in SOX10-inhibited melanoma cells reduced growth and migratory capacity, partly due to elevated expression of cyclin-dependent kinase inhibitor p21 and lack of NEDD9 induction. Conversely, opposite phenomenon was observed when SOX9 expression was further elevated to a range of high SOX9 expression levels in metastatic melanoma specimens, and that high levels of SOX9 can restore melanoma progression in the absence of SOX10 both in vitro and in vivo. In addition, overexpression of SOX9 can also promote invasiveness of the parental melanoma cells by modulating the expression of various matrix metalloproteinases. SOX10 or high SOX9 expression regulates melanoma mesenchymal migration through the NEDD9-mediated focal adhesion dynamics and Rho GTPase signaling.

CONCLUSIONS

These results unravel NEDD9 as a common target for SOX10 or high SOX9 to partly mediate their oncogenic events, and most importantly, reconcile previous discrepancies that suboptimal level of SOX9 expression is anti-metastatic whereas high level of SOX9 is metastatic in a heterogeneous population of melanoma.

The myelin protein PMP2 is regulated by SOX10 and drives melanoma cell invasion

The transcription factor sex determining region Y-box 10 (SOX10) plays a key role in the development of melanocytes and glial cells from neural crest precursors. SOX10 is involved in melanoma initiation, proliferation, invasion, and survival. However, specific mediators which impart its oncogenic properties remain widely unknown. To identify target genes of SOX10, we performed RNA sequencing after ectopic expression of SOX10 in human melanoma cells. Among nine differentially regulated genes, peripheral myelin protein 2 (PMP2) was consistently upregulated in several cell lines. Direct regulation of PMP2 by SOX10 was shown by chromatin immunoprecipitation, electrophoretic mobility shift, and luciferase reporter assays. Moreover, a coregulation of PMP2 by SOX10 and early growth response 2 in melanoma cells was found. Phenotypical investigation demonstrated that PMP2 expression can increase melanoma cell invasion. As PMP2 protein was detected only in a subset of melanoma cell lines, it might contribute to melanoma heterogeneity.

UBE2N Promotes Melanoma Growth via MEK/FRA1/SOX10 Signaling

UBE2N is a K63-specific ubiquitin conjugase linked to various immune disorders and cancer. Here, we demonstrate that UBE2N and its partners UBE2V1 and UBE2V2 are highly expressed in malignant melanoma. Silencing of UBE2N and its partners significantly decreased melanoma cell proliferation and subcutaneous tumor growth. This was accompanied by increased expression of E-cadherin, p16, and MC1R and decreased expression of melanoma malignancy markers including SOX10, Nestin, and ABCB5. Mass spectrometry-based phosphoproteomic analysis revealed that UBE2N loss resulted in distinct alterations to the signaling landscape: MEK/ERK signaling was impaired, FRA1 and SOX10 gene regulators were downregulated, and p53 and p16 tumor suppressors were upregulated. Similar to inhibition of UBE2N and MEK, silencing FRA1 decreased SOX10 expression and cell proliferation. Conversely, exogenous expression of active FRA1 increased pMEK and SOX10 expression, and restored anchorage-independent cell growth of cells with UBE2N loss. Systemic delivery of NSC697923, a small-molecule inhibitor of UBE2N, significantly decreased melanoma xenograft growth. These data indicate that UBE2N is a novel regulator of the MEK/FRA1/SOX10 signaling cascade and is indispensable for malignant melanoma growth. Our findings establish the basis for targeting UBE2N as a potential treatment strategy for melanoma.Significance: These findings identify ubiquitin conjugase UBE2N and its variant partners as novel regulators of MAPK signaling and potential therapeutic targets in melanoma. Cancer Res; 78(22); 6462-72. ©2018 AACR.

High MITF Expression Is Associated with Super-Enhancers and Suppressed by CDK7 Inhibition in Melanoma

Cutaneous melanoma is an aggressive tumor that accounts for most skin cancer deaths. Among the physiological barriers against therapeutic success is a strong survival program driven by genes such as MITF that specify melanocyte identity, a phenomenon known in melanoma biology as lineage dependency. MITF overexpression is occasionally explained by gene amplification, but here we show that super-enhancers are also important determinants of MITF overexpression in some melanoma cell lines and tumors. Although compounds that directly inhibit MITF are unavailable, a covalent CDK7 inhibitor, THZ1, has recently been shown to potently suppress the growth of various cancers through the depletion of master transcription-regulating oncogenes and the disruption of their attendant super-enhancers. We also show that melanoma cells are highly sensitive to CDK7 inhibition both in vitro and in vivo and that THZ1 can dismantle the super-enhancer apparatus at MITF and SOX10 in some cell lines, thereby extinguishing their intracellular levels. Our results show a dimension to MITF regulation in melanoma cells and point to CDK7 inhibition as a potential strategy to deprive oncogenic transcription and suppress tumor growth in melanoma.

BRN2, a POUerful driver of melanoma phenotype switching and metastasis

The POU domain family of transcription factors play a central role in embryogenesis and are highly expressed in neural crest cells and the developing brain. BRN2 is a class III POU domain protein that is a key mediator of neuroendocrine and melanocytic development and differentiation. While BRN2 is a central regulator in numerous developmental programs, it has also emerged as a major player in the biology of tumourigenesis. In melanoma, BRN2 has been implicated as one of the master regulators of the acquisition of invasive behaviour within the phenotype switching model of progression. As a mediator of melanoma cell phenotype switching, it coordinates the transition to a dedifferentiated, slow cycling and highly motile cell type. Its inverse expression relationship with MITF is believed to mediate tumour progression and metastasis within this model. Recent evidence has now outlined a potential epigenetic switching mechanism in melanoma cells driven by BRN2 expression that induces melanoma cell invasion. We summarize the role of BRN2 in tumour cell dissemination and metastasis in melanoma, while also examining it as a potential metastatic regulator in other tumour models.

A unique gene expression signature is significantly differentially expressed in tumor-positive or tumor-negative sentinel lymph nodes in patients with melanoma

The purpose of this study was to learn whether molecular characterization through gene expression profiling of node-positive and node-negative sentinel lymph nodes (SLNs) in patients with clinical stage I and II melanoma may improve the understanding of mechanisms of metastasis and identify gene signatures for SLNs/SLNs that correlate with diagnosis or clinical outcome. Gene expression profiling was performed on SLN biopsies of 48 (24 SLN and 24 SLN) patients (T3a/b-T4a/b) who underwent staging of SLNs using transcriptome profiling analysis on 5 μm sections of fresh SLNs. U133A 2.0 Affymetrix gene chips were used. Significance analysis of microarrays was used to test the association between gene expression level and SLN status. Genes with fold change more than 1.5 and q value less than 0.05 were considered differentially expressed. Pathway analysis was performed using Ingenuity Pathway Analysis. The Benjamini and Hochberg method was used to adjust for multiple testing in pathway analysis. We identified 89 probe sets that were significantly differentially expressed (1.5-27-fold; q<0.05). Upon performing the pathway analysis, it was found that 25 genes were common among the most significant and biologically relevant canonical pathways. The molecules and pathways that achieved differential expression of highest statistical significance were notably related to melanoma and its microenvironment and to signaling pathways implicated in immunosuppression and development of cancer. A 25-gene signature is significantly differentially expressed between SLN and SLN and is related to melanoma oncogenesis and immunosuppression. The identified expression profile provides a signature of melanoma nodal involvement. These findings warrant further investigation into the mechanisms of metastasis, melanoma metastasis diagnosis, and prediction of outcome.