Ligand-activated BMP signaling inhibits cell differentiation and death to promote melanoma

GDF6 in gastric cancer upregulated by helicobacter pylori induces epithelial-mesenchymal translation via the TGF-β/SMAD3 signaling pathway

OBJECTIVE

To investigate the association between Helicobacter pylori infection and GDF6 expression in gastric cancer patients, and to determine its influence on prognosis and resistance to capecitabine.

METHODS

Tumor and adjacent non-tumor tissues were collected from 148 gastric cancer patients who underwent surgery in our department from October 2019 to June 2022. Of these patients, 78 tested positive for Helicobacter pylori and 70 tested negative. Hematoxylin-eosin (HE) and immunofluorescence staining were utilized to quantify GDF6 expression in cancerous and adjacent tissues. Patient prognosis was monitored via follow-up. Western blotting analyzed GDF6 expression in common gastric cancer cell lines. HGC27 cells exhibiting high GDF6 expression and BGC823 cells with low expression were used to create GDF6-silenced and overexpressed cell lines. The impact of GDF6 on the proliferation, migration, invasion, and cloning abilities of gastric cancer cells was evaluated using the CCK-8 assay, scratch test, Transwell assay, and plate colony formation assay. Fluorescent quantitative PCR and Western blotting assessed the effects of GDF6 levels on epithelial-mesenchymal transition (EMT) and tumor cell stemness.

RESULTS

GDF6 expression in gastric cancer tissues was significantly correlated with cancer grading and staging (P<0.05). Helicobacter pylori-positive tissues exhibited significantly higher GDF6 expression levels than negative samples (P<0.05). Kaplan-Meier survival analysis indicated that high GDF6 expression was associated with poor survival prognosis. Overexpressed GDF6 enhanced the proliferation, migration, and invasion abilities of gastric cancer cells, while silencing GDF6 yielded opposite results. Increased GDF6 expression upregulated TGF-β expression and the phosphorylation levels of SMAD3, leading to an elevation in mesenchymal cell markers N-cadherin, vimentin, and a reduction in epithelial cell markers cytokeratins, E-cadherin. Moreover, high GDF6 levels contributed to increased resistance to capecitabine and enhanced the expression of tumor stem cell markers Nanog, Sox-2, Oct-4, CD44, amplifying tumor cell stemness.

CONCLUSION

Helicobacter pylori infection is associated with increased GDF6 expression in gastric cancer tissue, correlating with poor survival prognosis. Elevated GDF6 expression promotes the proliferation, migration, and invasion abilities of gastric cancer cells, facilitates EMT via the TGF-β/SMAD3 pathway, and intensifies cell stemness and capecitabine resistance. Consequently, GDF6 presents itself as a potential new target for gastric cancer treatment.

DATA AVAILABILITY STATEMENT

The data that support the findings of this study are available from the corresponding author upon reasonable request.

The PIEZO1/miR-155-5p/GDF6/SMAD2/3 signaling axis is involved in inducing the occurrence and progression of osteoarthritis under excessive mechanical stress

OBJECTIVE

To elucidate the molecular mechanism of overloading-induced osteoarthritis (OA) and to find a novel therapeutic target.

METHODS

We utilized human cartilage specimens, mouse chondrocytes, a destabilization of the medial meniscus (DMM) mouse model, and a mouse hindlimb weight-bearing model to validate the role of overloading on chondrocyte senescence and OA development. Then, we observed the effect of PIEZO1-miR-155-5p-GDF6-SMAD2/3 signaling axis on the preservation of joint metabolic homeostasis under overloading in vivo, in vitro and ex vivo by qPCR, Western blot, enzyme-linked immunosorbent assay (ELISA), immunohistochemistry, immunofluorescence, SA-β-gal staining, CCK8 assay, et al. Finally, we verified the therapeutic effects of intra-articular injection of miR-155-5p inhibitor or recombinant GDF6 on the murine overloading-induced OA models.

RESULTS

Chondrocytes sensesed the mechanical overloading through PIEZO1 and up-regulated miR-155-5p expression. MiR-155-5p mimics could copy the effects of overloading-induced chondrocyte senescence and OA. Additionally, miR-155-5p could suppress the mRNA expression of Gdf6-Smad2/3 in various tissues within the joint. Overloading could disrupt joint metabolic homeostasis by downregulating the expression of anabolism indicators and upregulating the expression of catabolism indicators in the chondrocytes and synoviocytes, while miR-155-5p inhibition or GDF6 supplementation could exert an antagonistic effect by preserving the joint homeostasis. Finally, in the in vivo overloading models, intra-articular injection of miR-155-5p inhibitor or recombinant GDF6 could significantly mitigate the severity of impending OA and lessened the progression of existing OA.

CONCLUSION

GDF6 overexpression or miR-155-5p inhibition could attenuate overloading-induced chondrocyte senescence and OA through the PIEZO1-miR-155-5p-GDF6-SMAD2/3 signaling pathway. Our study provides a new therapeutic target for the treatment of overloading-induced OA.

Immune escape and metastasis mechanisms in melanoma: breaking down the dichotomy

Melanoma is one of the most lethal neoplasms of the skin. Despite the revolutionary introduction of immune checkpoint inhibitors, metastatic spread, and recurrence remain critical problems in resistant cases. Melanoma employs a multitude of mechanisms to subvert the immune system and successfully metastasize to distant organs. Concerningly, recent research also shows that tumor cells can disseminate early during melanoma progression and enter dormant states, eventually leading to metastases at a future time. Immune escape and metastasis have previously been viewed as separate phenomena; however, accumulating evidence is breaking down this dichotomy. Recent research into the progressive mechanisms of melanoma provides evidence that dedifferentiation similar to classical epithelial to mesenchymal transition (EMT), genes involved in neural crest stem cell maintenance, and hypoxia/acidosis, are important factors simultaneously involved in immune escape and metastasis. The likeness between EMT and early dissemination, and differences, also become apparent in these contexts. Detailed knowledge of the mechanisms behind "dual drivers" simultaneously promoting metastatically inclined and immunosuppressive environments can yield novel strategies effective in disabling multiple facets of melanoma progression. Furthermore, understanding progression through these drivers may provide insight towards novel treatments capable of preventing recurrence arising from dormant dissemination or improving immunotherapy outcomes.

A novel DNA methylation signature revealed GDF6 and RCC1 as potential prognostic biomarkers correlated with cell proliferation in clear cell renal cell carcinoma

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) accounts for the majority (80%-90%) of renal cell carcinoma (RCC) patients at the time of diagnosis, and approximately 15% of ccRCC patients will develop distant metastasis or recurrence during their lifetime. Increasing number of studies have revealed that the aberrant DNA methylations is closely correlated with the tumorigenesis in ccRCC.

RESULTS

In this study, we utilized a LASSO (least absolute shrinkage and selection operator) model to identify a combination of 13 probes-based DNA methylation signature that associated with the progression-free survival (PFS) of ccRCC patients. First, differentially methylated regions (CpGs) related to PFS and phenotypes were identified. Next, prognostic DNA methylation probes were selected from the differentially methylated probes (DMPs) and calculated risk scores to stratify patients with ccRCC. The performance of this signature was validated in an independent testing set using various analyses, including Kaplan-Meier analysis for PFS and receiver operating characteristic (ROC) curve analysis. Based on our 13-DNA methylation probes signature, ccRCC patients were successfully stratified into high- and low-risk groups. Combining DNA methylation signature with clinical variables such as T stage, M stage and tumor grade could further improve the accuracy of prediction. Moreover, we highlight two molecular biomarkers (RCC1 and GDF6) corresponding to our probes. Invitro experiments showed that knockdown of RCC1 or GDF6 in ccRCC cell lines reduced cell proliferation, which indicated that both biomarkers are associated with tumorigenesis.

CONCLUSIONS

The 13-probes-based DNA methylation signature has the potential to serve as an independent tool for survival outcome improvement and treatment strategy selection for ccRCC patients. In addition, our findings suggest that RCC1 and GDF6 may serve as promising markers for ccRCC.

BMP signaling in cancer stemness and differentiation

The BMP (Bone morphogenetic protein) signaling pathway plays a central role in metazoan biology, intricately shaping embryonic development, maintaining tissue homeostasis, and influencing disease progression. In the context of cancer, BMP signaling exhibits context-dependent dynamics, spanning from tumor suppression to promotion. Cancer stem cells (CSCs), a modest subset of neoplastic cells with stem-like attributes, exert substantial influence by steering tumor growth, orchestrating therapy resistance, and contributing to relapse. A comprehensive grasp of the intricate interplay between CSCs and their microenvironment is pivotal for effective therapeutic strategies. Among the web of signaling pathways orchestrating cellular dynamics within CSCs, BMP signaling emerges as a vital conductor, overseeing CSC self-renewal, differentiation dynamics, and the intricate symphony within the tumor microenvironment. Moreover, BMP signaling's influence in cancer extends beyond CSCs, intricately regulating cellular migration, invasion, and metastasis. This multifaceted role underscores the imperative of comprehending BMP signaling's contributions to cancer, serving as the foundation for crafting precise therapies to navigate multifaceted challenges posed not only by CSCs but also by various dimensions of cancer progression. This article succinctly encapsulates the diverse roles of the BMP signaling pathway across different cancers, spanning glioblastoma multiforme (GBM), diffuse intrinsic pontine glioma (DIPG), colorectal cancer, acute myeloid leukemia (AML), lung cancer, prostate cancer, and osteosarcoma. It underscores the necessity of unraveling underlying mechanisms and molecular interactions. By delving into the intricate tapestry of BMP signaling's engagement in cancers, researchers pave the way for meticulously tailored therapies, adroitly leveraging its dualistic aspects-whether as a suppressor or promoter-to effectively counter the relentless march of tumor progression.

Functional analysis of recurrent CDC20 promoter variants in human melanoma

Small nucleotide variants in non-coding regions of the genome can alter transcriptional regulation, leading to changes in gene expression which can activate oncogenic gene regulatory networks. Melanoma is heavily burdened by non-coding variants, representing over 99% of total genetic variation, including the well-characterized TERT promoter mutation. However, the compendium of regulatory non-coding variants is likely still functionally under-characterized. We developed a pipeline to identify hotspots, i.e. recurrently mutated regions, in melanoma containing putatively functional non-coding somatic variants that are located within predicted melanoma-specific regulatory regions. We identified hundreds of statistically significant hotspots, including the hotspot containing the TERT promoter variants, and focused on a hotspot in the promoter of CDC20. We found that variants in the promoter of CDC20, which putatively disrupt an ETS motif, lead to lower transcriptional activity in reporter assays. Using CRISPR/Cas9, we generated an indel in the CDC20 promoter in human A375 melanoma cell lines and observed decreased expression of CDC20, changes in migration capabilities, increased growth of xenografts, and an altered transcriptional state previously associated with a more proliferative and less migratory state. Overall, our analysis prioritized several recurrent functional non-coding variants that, through downregulation of CDC20, led to perturbation of key melanoma phenotypes.

GABA Regulates Electrical Activity and Tumor Initiation in Melanoma

Oncogenes can initiate tumors only in certain cellular contexts, which is referred to as oncogenic competence. In melanoma, whether cells in the microenvironment can endow such competence remains unclear. Using a combination of zebrafish transgenesis coupled with human tissues, we demonstrate that GABAergic signaling between keratinocytes and melanocytes promotes melanoma initiation by BRAFV600E. GABA is synthesized in melanoma cells, which then acts on GABA-A receptors in keratinocytes. Electron microscopy demonstrates specialized cell-cell junctions between keratinocytes and melanoma cells, and multielectrode array analysis shows that GABA acts to inhibit electrical activity in melanoma/keratinocyte cocultures. Genetic and pharmacologic perturbation of GABA synthesis abrogates melanoma initiation in vivo. These data suggest that GABAergic signaling across the skin microenvironment regulates the ability of oncogenes to initiate melanoma.

SIGNIFICANCE

This study shows evidence of GABA-mediated regulation of electrical activity between melanoma cells and keratinocytes, providing a new mechanism by which the microenvironment promotes tumor initiation. This provides insights into the role of the skin microenvironment in early melanomas while identifying GABA as a potential therapeutic target in melanoma. See related commentary by Ceol, p. 2128. This article is featured in Selected Articles from This Issue, p. 2109.

Bone Morphogenetic Protein 13 Has Protumorigenic Effects on Hepatocellular Carcinoma Cells In Vitro

Activated hepatic stellate cells (HSCs) play a key role in hepatic fibrosis and, thus, build the "soil" for hepatocarcinogenesis. Furthermore, HSCs are known to promote the progression of hepatocellular carcinoma (HCC), but the molecular mechanisms are only incompletely understood. Recently, we newly described the expression of bone morphogenetic protein 13 (BMP13) by HSCs in fibrotic liver tissue. In addition, BMP13 has mostly been studied in the context of cartilage and bone repair, but not in liver disease or cancer. Thus, we aimed to analyze the expression and function of BMP13 in HCC. Expression analyses revealed high BMP13-expression in activated human HSCs, but not in human HCC-cell-lines. Furthermore, analysis of human HCC tissues showed a significant correlation between BMP13 and α-smooth muscle actin (α-SMA), and immunofluorescence staining confirmed the co-localization of BMP13 and α-SMA, indicating activated HSCs as the cellular source of BMP13 in HCC. Stimulation of HCC cells with recombinant BMP13 increased the expression of the inhibitors of differentiation 1 (ID1) and 2 (ID2), which are known targets of BMP-signaling and cell-cycle promotors. In line with this, BMP13-stimulation caused an induced SMAD 1/5/9 and extracellular signal-regulated kinase (ERK) phosphorylation, as well as reduced expression of cyclin-dependent kinase inhibitors 1A (CDKN1A) and 2A (CDKN2A). Furthermore, stimulation with recombinant BMP13 led to increased proliferation and colony size formation of HCC cells in clonogenicity assays. The protumorigenic effects of BMP13 on HCC cells were almost completely abrogated by the small molecule dorsomorphin 1 (DMH1), which selectively blocks the intracellular kinase domain of ALK2 and ALK3, indicating that BMP13 acts via these BMP type I receptors on HCC cells. In summary, this study newly identifies stroma-derived BMP13 as a potential new tumor promotor in HCC and indicates this secreted growth-factor as a possible novel therapeutic target in HCC.

Comparative Tumor Microenvironment Analysis of Primary and Recurrent Ovarian Granulosa Cell Tumors

Adult-type granulosa cell tumors (aGCT) are rare ovarian sex cord tumors with few effective treatments for recurrent disease. The objective of this study was to characterize the tumor microenvironment (TME) of primary and recurrent aGCTs and to identify correlates of disease recurrence. Total RNA sequencing (RNA-seq) was performed on 24 pathologically confirmed, cryopreserved aGCT samples, including 8 primary and 16 recurrent tumors. After read alignment and quality-control filtering, DESeq2 was used to identify differentially expressed genes (DEG) between primary and recurrent tumors. Functional enrichment pathway analysis and gene set enrichment analysis was performed using "clusterProfiler" and "GSVA" R packages. TME composition was investigated through the analysis and integration of multiple published RNA-seq deconvolution algorithms. TME analysis results were externally validated using data from independent previously published RNA-seq datasets. A total of 31 DEGs were identified between primary and recurrent aGCTs. These included genes with known function in hormone signaling such as LHCGR and INSL3 (more abundant in primary tumors) and CYP19A1 (more abundant in recurrent tumors). Gene set enrichment analysis revealed that primarily immune-related and hormone-regulated gene sets expression was increased in recurrent tumors. Integrative TME analysis demonstrated statistically significant depletion of cancer-associated fibroblasts in recurrent tumors. This finding was confirmed in multiple independent datasets.

IMPLICATIONS

Recurrent aGCTs exhibit alterations in hormone pathway gene expression as well as decreased infiltration of cancer-associated fibroblasts, suggesting dual roles for hormonal signaling and TME remodeling underpinning disease relapse.

BORIS/CTCFL-mediated chromatin accessibility alterations promote a pro-invasive transcriptional signature in melanoma cells

Melanoma is the deadliest form of skin cancer, due to its tendency to metastasize early. Brother of regulator of imprinted sites (BORIS), also known as CCCTC binding factor-like (CTCFL), is a transcription regulator that becomes ectopically expressed in melanoma. We recently showed that BORIS contributes to melanoma phenotype switching by altering the gene expression program of melanoma cells from an intermediate melanocytic state toward a more mesenchymal-like state. However, the mechanism underlying this transcriptional switch remains unclear. Here, ATAC-seq was used to study BORIS-mediated chromatin accessibility alterations in melanoma cells harboring an intermediate melanocytic state. The gene set that gained promoter accessibility, following ectopic BORIS expression, showed enrichment for biological processes associated with melanoma invasion, while promoters of genes associated with proliferation showed reduced accessibility. Integration of ATAC-seq and RNA-seq data demonstrated that increased chromatin accessibility was associated with transcriptional upregulation of genes involved in tumor progression processes, and the aberrant activation of oncogenic transcription factors, while reduced chromatin accessibility and downregulated genes were associated with repressed activity of tumor suppressors and proliferation factors. Together, these findings indicate that BORIS mediates transcriptional reprogramming in melanoma cells by altering chromatin accessibility and gene expression, shifting the cellular transcription landscape of melanoma cells toward a mesenchymal-like genetic signature.

Mesenchymal/stromal stem cells: necessary factors in tumour progression

Mesenchymal/stromal stem cells (MSCs) are a crucial component of the tumour microenvironment (TME). They can be recruited from normal tissues into the TME and educated by tumour cells to transform into tumour-associated MSCs, which are oncogenic cells that promote tumour development and progression by impacting or transforming into various kinds of cells, such as immune cells and endothelial cells. Targeting MSCs in the TME is a novel strategy to prevent malignant processes. Exosomes, as communicators, carry various RNAs and proteins and thus link MSCs and the TME, which provides options for improving outcomes and developing targeted treatment.

CDK12 is hyperactivated and a synthetic-lethal target in BRAF-mutated melanoma

Melanoma is the deadliest form of skin cancer and considered intrinsically resistant to chemotherapy. Nearly all melanomas harbor mutations that activate the RAS/mitogen-activated protein kinase (MAPK) pathway, which contributes to drug resistance via poorly described mechanisms. Herein we show that the RAS/MAPK pathway regulates the activity of cyclin-dependent kinase 12 (CDK12), which is a transcriptional CDK required for genomic stability. We find that melanoma cells harbor constitutively high CDK12 activity, and that its inhibition decreases the expression of long genes containing multiple exons, including many genes involved in DNA repair. Conversely, our results show that CDK12 inhibition promotes the expression of short genes with few exons, including many growth-promoting genes regulated by the AP-1 and NF-κB transcription factors. Inhibition of these pathways strongly synergize with CDK12 inhibitors to suppress melanoma growth, suggesting promising drug combinations for more effective melanoma treatment.

Construction of mRNA prognosis signature associated with differentially expressed genes in early stage of stomach adenocarcinomas based on TCGA and GEO datasets

Background

Stomach adenocarcinomas (STAD) are the most common malignancy of the human digestive system and represent the fourth leading cause of cancer-related deaths. As early-stage STAD are generally mild or asymptomatic, patients with advanced STAD have short overall survival. Early diagnosis of STAD has a considerable influence on clinical outcomes.

Methods

The mRNA expression data and clinical indicators of STAD and normal tissues were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. The gene expression differences were analyzed by R packages, and gene function enrichment analysis was performed. Kaplan–Meier method and univariate Cox proportional risk regression analysis were used to screen differential expressed genes (DEGs) related to survival of STAD patients. Multivariate Cox proportional risk regression analysis was used to further screen and determine the prognostic DEGs in STAD patients, and to construct a multigene prognostic prediction signature. The accuracy of predictive signature was tested by receiver operating characteristic (ROC) curve software package, and the nomogram of patients with STAD was drawn. Cox regression was used to investigate the correlation between multigene prognostic signature and clinical factors. The predictive performance of this model was compared with two other models proposed in previous studies using KM survival analysis, ROC curve analysis, Harrell consistency index and decision curve analysis (DCA). qRT-PCR and Western blot were used to verify the expression levels of prognostic genes. The pathways and functions of possible involvement of features were predicted using the GSEA method.

Results

A total of 569 early-stage specific DEGs were retrieved from TCGA-STAD dataset, including 229 up-regulated genes and 340 down-regulated genes. Enrichment analysis showed that the early-stage specific DEGs were associated with cytokine–cytokine receptor interaction, neuroactive ligand–receptor interaction, and calcium signaling pathway. Multiple Cox regression algorithm was used to identify 10 early-stage specific DEGs associated with overall survival (P < 0.01) of STAD patients, and a multi-mRNA prognosis signature was established. The patients were divided into high-risk group and low-risk group according to the risk score. In the training set, the prognostic signature was positively correlated with tumor size and stage (P < 0.05), survival curve (P < 0.001) and time-dependent ROC (AUC = 0.625). In the training dataset and test dataset, the both signatures had good predictive efficiencies. Cox regression and DCA analysis revealed that the prognostic signature was an independent factor and had a better predict effect than the conventional TNM stage classification method and the earlier published biomarkers on the prognosis of STAD patients.

Conclusion

In this study, based on the early-stage specifically expressed genes, the prognostic signature constructed through TCGA and GEO datasets may become an indicator for clinical prognosis assessment of STAD and a new strategy for targeted therapy in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40001-022-00827-4.

Oncogenic BRAF induces whole-genome doubling through suppression of cytokinesis

Melanomas and other solid tumors commonly have increased ploidy, with near-tetraploid karyotypes being most frequently observed. Such karyotypes have been shown to arise through whole-genome doubling events that occur during early stages of tumor progression. The generation of tetraploid cells via whole-genome doubling is proposed to allow nascent tumor cells the ability to sample various pro-tumorigenic genomic configurations while avoiding the negative consequences that chromosomal gains or losses have in diploid cells. Whereas a high prevalence of whole-genome doubling events has been established, the means by which whole-genome doubling arises is unclear. Here, we find that BRAFV600E, the most common mutation in melanomas, can induce whole-genome doubling via cytokinesis failure in vitro and in a zebrafish melanoma model. Mechanistically, BRAFV600E causes decreased activation and localization of RhoA, a critical cytokinesis regulator. BRAFV600E activity during G1/S phases of the cell cycle is required to suppress cytokinesis. During G1/S, BRAFV600E activity causes inappropriate centriole amplification, which is linked in part to inhibition of RhoA and suppression of cytokinesis. Together these data suggest that common abnormalities of melanomas linked to tumorigenesis - amplified centrosomes and whole-genome doubling events - can be induced by oncogenic BRAF and other mutations that increase RAS/MAPK pathway activity.

Single-Cell RNA Sequencing Unravels Upregulation of Immune Cell Crosstalk in Relapsed Pediatric Ependymoma

Ependymoma (EPN) is a malignant glial tumor occurring throughout central nervous system, which commonly presents in children. Although recent studies have characterized EPN samples at both the bulk and single-cell level, intratumoral heterogeneity across subclones remains a confounding factor that impedes understanding of EPN biology. In this study, we generated a high-resolution single-cell dataset of pediatric ependymoma with a particular focus on the comparison of subclone differences within tumors and showed upregulation of cilium-associated genes in more highly differentiated subclone populations. As a proxy to traditional pseudotime analysis, we applied a novel trajectory scoring method to reveal cellular compositions associated with poor survival outcomes across primary and relapsed patients. Furthermore, we identified putative cell–cell communication features between relapsed and primary samples and showed upregulation of pathways associated with immune cell crosstalk. Our results revealed both inter- and intratumoral heterogeneity in EPN and provided a framework for studying transcriptomic signatures of individual subclones at single-cell resolution.

Oxidative stress from DGAT1 oncoprotein inhibition in melanoma suppresses tumor growth when ROS defenses are also breached

Dysregulated cellular metabolism is a cancer hallmark for which few druggable oncoprotein targets have been identified. Increased fatty acid (FA) acquisition allows cancer cells to meet their heightened membrane biogenesis, bioenergy, and signaling needs. Excess FAs are toxic to non-transformed cells but surprisingly not to cancer cells. Molecules underlying this cancer adaptation may provide alternative drug targets. Here, we demonstrate that diacylglycerol O-acyltransferase 1 (DGAT1), an enzyme integral to triacylglyceride synthesis and lipid droplet formation, is frequently up-regulated in melanoma, allowing melanoma cells to tolerate excess FA. DGAT1 over-expression alone transforms p53-mutant zebrafish melanocytes and co-operates with oncogenic BRAF or NRAS for more rapid melanoma formation. Antagonism of DGAT1 induces oxidative stress in melanoma cells, which adapt by up-regulating cellular reactive oxygen species defenses. We show that inhibiting both DGAT1 and superoxide dismutase 1 profoundly suppress tumor growth through eliciting intolerable oxidative stress.

Long-term non-invasive drug treatments in adult zebrafish that lead to melanoma drug resistance

Zebrafish embryos are widely used for drug discovery, however, administering drugs to adult zebrafish is limited by current protocols that can cause stress. Here, we developed a drug formulation and administration method for adult zebrafish by producing food-based drug pellets that are consumed voluntarily. We applied this to zebrafish with BRAF-mutant melanoma, a model that has significantly advanced our understanding of melanoma progression, but not of drug resistance due to the limitations of current treatment methods. Zebrafish with melanomas responded to short-term, precise and daily dosing with drug pellets made with the BRAF^V600E inhibitor, vemurafenib. On-target drug efficacy was determined by phospho-Erk staining. Continued drug treatment led to the emergence, for the first time in zebrafish, of acquired drug resistance and melanoma relapse, modelling the responses seen in melanoma patients. This method presents a controlled, non-invasive approach that permits long-term drug studies and can be widely applied to adult zebrafish models.

Summary: We developed a precise drug-pellet dosing method for adult zebrafish melanoma models that can be widely applied to drug discovery in adult zebrafish.

Construction and Validation of a Ferroptosis-Related Prognostic Signature for Melanoma Based on Single-Cell RNA Sequencing

Melanoma, the deadliest type of skin cancer, is on the rise globally. The generally poor prognosis makes melanoma still an enormous public health problem. Ferroptosis is a newly emerging form of iron-dependent regulated cell death, which has been implicated in the development and treatment of several tumors. However, whether there is a connection between ferroptosis-related genes and the prognosis of melanoma patients remains an enigma. In the present study, we identified a ferroptosis-related genes signature to predict the prognosis of melanoma patients by analyzing single-cell RNA-sequencing data from the Gene Expression Omnibus (GEO). Single-cell trajectory analysis was performed to explore malignant differentiation. CellChat was used to investigate intercellular communications in melanoma. Collectively, a novel four-gene signature (CP, MAP1LC3A, transferrin, and TP53) was constructed for prognosis prediction. COX proportional hazards regression analysis showed that the established ferroptosis-associated risk model was an independent prognostic predictor for melanoma patients (HR = 2.3293; 95%CI 1.1528–4.706) (p < 0.018). Patients with low-risk scores had significantly better overall survival (OS) than those with high-risk scores in The Cancer Genome Atlas, GSE59455, and GSE22153 dataset (p = 0.0015, p = 0.031, p = 0.077). Furthermore, the gene expression level of the four genes were verified in multistrain melanoma cell lines and normal human epidermal melanocytes (NHEM). The protein expression level of the four genes in clinical samples were further verified in the Human Protein Atlas (HPA) databases. Taken together, our study identified the prognostic significance of the ferroptosis-related genes in melanoma and developed a novel four-gene prognostic signature, which may shed light on the prognostic assessment and clinical decision making for melanoma patients.

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.

Shifting the focus of zebrafish toward a model of the tumor microenvironment

Cancer cells exist in a complex ecosystem with numerous other cell types in the tumor microenvironment (TME). The composition of this tumor/TME ecosystem will vary at each anatomic site and affects phenotypes such as initiation, metastasis, and drug resistance. A mechanistic understanding of the large number of cell-cell interactions between tumor and TME requires models that allow us to both characterize as well as genetically perturb this complexity. Zebrafish are a model system optimized for this problem, because of the large number of existing cell-type-specific drivers that can label nearly any cell in the TME. These include stromal cells, immune cells, and tissue resident normal cells. These cell-type-specific promoters/enhancers can be used to drive fluorophores to facilitate imaging and also CRISPR cassettes to facilitate perturbations. A major advantage of the zebrafish is the ease by which large numbers of TME cell types can be studied at once, within the same animal. While these features make the zebrafish well suited to investigate the TME, the model has important limitations, which we also discuss. In this review, we describe the existing toolset for studying the TME using zebrafish models of cancer and highlight unique biological insights that can be gained by leveraging this powerful resource.

TRANSCRIPTIONAL PROFILE AND CHROMATIN ACCESSIBILITY IN ZEBRAFISH MELANOCYTES AND MELANOMA TUMORS

Transcriptional and epigenetic characterization of melanocytes and melanoma cells isolated from their in vivo context promises to unveil key differences between these developmentally related normal and cancer cell populations. We therefore engineered an enhanced Danio rerio (zebrafish) melanoma model with fluorescently labeled melanocytes to allow for isolation of normal (wild type) and premalignant (BRAF^V600E-mutant) populations for comparison to fully transformed BRAF^V600E-mutant, p53 loss-of-function melanoma cells. Using fluorescence-activated cell sorting to isolate these populations, we performed high-quality RNA- and ATAC-seq on sorted zebrafish melanocytes vs. melanoma cells, which we provide as a resource here. Melanocytes had consistent transcriptional and accessibility profiles, as did melanoma cells. Comparing melanocytes and melanoma, we note 4128 differentially expressed genes and 56,936 differentially accessible regions with overall gene expression profiles analogous to human melanocytes and the pigmentation melanoma subtype. Combining the RNA- and ATAC-seq data surprisingly revealed that increased chromatin accessibility did not always correspond with increased gene expression, suggesting that though there is widespread dysregulation in chromatin accessibility in melanoma, there is a potentially more refined gene expression program driving cancerous melanoma. These data serve as a resource to identify candidate regulators of the normal vs. diseased states in a genetically controlled in vivo context.

Bone morphogenetic protein induces bone invasion of melanoma by epithelial-mesenchymal transition via the Smad1/5 signaling pathway

Oral malignant melanoma, which frequently invades the hard palate or maxillary bone, is extremely rare and has a poor prognosis. Bone morphogenetic protein (BMP) is abundantly expressed in bone matrix and is highly expressed in malignant melanoma, inducing an aggressive phenotype. We examined the role of BMP signaling in the acquisition of an aggressive phenotype in melanoma cells in vitro and in vivo. In five cases, immunohistochemistry indicated the phosphorylation of Smad1/5 (p-Smad1/5) in the nuclei of melanoma cells. In the B16 mouse and A2058 human melanoma cell lines, BMP2, BMP4, or BMP7 induces morphological changes accompanied by the downregulation of E-cadherin, and the upregulation of N-cadherin and Snail, markers of epithelial-mesenchymal transition (EMT). BMP2 also stimulates cell invasion by increasing matrix metalloproteinase activity in B16 cells. These effects were canceled by the addition of LDN193189, a specific inhibitor of Smad1/5 signaling. In vivo, the injection of B16 cells expressing constitutively activated ALK3 enhanced zygoma destruction in comparison to empty B16 cells by increasing osteoclast numbers. These results suggest that the activation of BMP signaling induces EMT, thus driving the acquisition of an aggressive phenotype in malignant melanoma.

Immunogenomic Analyses of the Prognostic Predictive Model for Patients With Renal Cancer

Background

Renal cell carcinoma (RCC) is associated with poor prognostic outcomes. The current stratifying system does not predict prognostic outcomes and therapeutic benefits precisely for RCC patients. Here, we aim to construct an immune prognostic predictive model to assist clinician to predict RCC prognosis.

Methods

Herein, an immune prognostic signature was developed, and its predictive ability was confirmed in the kidney renal clear cell carcinoma (KIRC) cohorts based on The Cancer Genome Atlas (TCGA) dataset. Several immunogenomic analyses were conducted to investigate the correlations between immune risk scores and immune cell infiltrations, immune checkpoints, cancer genotypes, tumor mutational burden, and responses to chemotherapy and immunotherapy.

Results

The immune prognostic signature contained 14 immune-associated genes and was found to be an independent prognostic factor for KIRC. Furthermore, the immune risk score was established as a novel marker for predicting the overall survival outcomes for RCC. The risk score was correlated with some significant immunophenotypic factors, including T cell infiltration, antitumor immunity, antitumor response, oncogenic pathways, and immunotherapeutic and chemotherapeutic response.

Conclusions

The immune prognostic, predictive model can be effectively and efficiently used in the prediction of survival outcomes and immunotherapeutic responses of RCC patients.

Isotope tracing in adult zebrafish reveals alanine cycling between melanoma and liver

The cell-intrinsic nature of tumor metabolism has become increasingly well characterized. The impact that tumors have on systemic metabolism, however, has received less attention. Here, we used adult zebrafish harboring BRAFV600E-driven melanoma to study the effect of cancer on distant tissues. By applying metabolomics and isotope tracing, we found that melanoma consume ~15 times more glucose than other tissues measured. Despite this burden, circulating glucose levels were maintained in disease animals by a tumor-liver alanine cycle. Excretion of glucose-derived alanine from tumors provided a source of carbon for hepatic gluconeogenesis and allowed tumors to remove excess nitrogen from branched-chain amino acid catabolism, which we found to be activated in zebrafish and human melanoma. Pharmacological inhibition of the tumor-liver alanine cycle in zebrafish reduced tumor burden. Our findings underscore the significance of metabolic crosstalk between tumors and distant tissues and establish the adult zebrafish as an attractive model to study such processes.

The impact of single-nucleotide polymorphisms on intravesical recurrence after bacillus Calmette-Guérin therapy for non-muscle invasive bladder cancer in a genome-wide association study

OBJECTIVE

Bacillus Calmette-Guérin (BCG) instillation therapy is widely used to reduce intravesical recurrence in non-muscle invasive bladder cancer (NMIBC). In this study, we aimed to reveal the genetic variations associated with intravesical recurrence after BCG therapy for NMIBC in a genome-wide association study (GWAS).

MATERIALS AND METHODS

This study included Japanese patients with NMIBC, in whom genomic DNA was obtained from whole blood samples. The association between genetic variation and treatment failure was analyzed by GWAS in 44 patients treated with BCG instillation as a discovery cohort. Candidate single-nucleotide polymorphisms (SNPs) were examined separately in 47 patients treated with BCG instillation and in 62 patients treated with chemotherapeutic agent instillation as validation studies.

RESULTS

Among the 44 patients in the discovery cohort, 14 cases experienced intravesical recurrent diseases. GWAS identified 12 candidate SNPs (rs9374832, rs35176001, rs363765, rs2127120, rs4277759, rs73664140, rs1607282, rs12141654, rs4541358, rs6986852, rs12373386, and rs17637903). In the validation study, a genetic risk stratification model using the number of risk alleles in rs363765 and rs6986852 discriminated the risk of intravesical recurrence after BCG therapy, but not after non-BCG therapy.

CONCLUSION

This study suggested that several SNPs were associated with intravesical recurrence after BCG therapy for NMIBC. A genetic risk model may be useful to predict intravesical recurrence after BCG therapy, warranting further research and development for clinical application.

Evidence of Cooperation between Hippo Pathway and RAS Mutation in Thyroid Carcinomas

Thyroid cancer incidences have been steadily increasing worldwide and are projected to become the fourth leading cancer diagnosis by 2030. Improved diagnosis and prognosis predictions for this type of cancer depend on understanding its genetic bases and disease biology. RAS mutations have been found in a wide range of thyroid tumors, from benign to aggressive thyroid carcinomas. Based on that and in vivo studies, it has been suggested that RAS cooperates with other driver mutations to induce tumorigenesis. This study aims to identify genetic alterations or pathways that cooperate with the RAS mutation in the pathogenesis of thyroid cancer. From a cohort of 120 thyroid carcinomas, 11 RAS-mutated samples were identified. The samples were subjected to RNA-Sequencing analyses. The mutation analysis in our eleven RAS-positive cases uncovered that four genes that belong to the Hippo pathway were mutated. The gene expression analysis revealed that this pathway was dysregulated in the RAS-positive samples. We additionally explored the mutational status and expression profiling of 60 RAS-positive papillary thyroid carcinomas (PTC) from The Cancer Genome Atlas (TCGA) cohort. Altogether, the mutational landscape and pathway enrichment analysis (gene set enrichment analysis (GSEA) and Kyoto Encyclopedia of Genes and Genome (KEGG)) detected the Hippo pathway as dysregulated in RAS-positive thyroid carcinomas. Finally, we suggest a crosstalk between the Hippo and other signaling pathways, such as Wnt and BMP.

Melanoma models for the next generation of therapies

There is a lack of appropriate melanoma models that can be used to evaluate the efficacy of novel therapeutic modalities. Here, we discuss the current state of the art of melanoma models including genetically engineered mouse, patient-derived xenograft, zebrafish, and ex vivo and in vitro models. We also identify five major challenges that can be addressed using such models, including metastasis and tumor dormancy, drug resistance, the melanoma immune response, and the impact of aging and environmental exposures on melanoma progression and drug resistance. Additionally, we discuss the opportunity for building models for rare subtypes of melanomas, which represent an unmet critical need. Finally, we identify key recommendations for melanoma models that may improve accuracy of preclinical testing and predict efficacy in clinical trials, to help usher in the next generation of melanoma therapies.

Bone Morphogenic Protein Signaling and Melanoma

OPINION STATEMENT

Malignant melanoma is a deadly form of skin cancer caused by neoplastic transformation of melanocytic cells. Despite recent progress in melanoma therapy, by inhibition of activated oncogenes or immunotherapy, survival rate for metastatic melanoma patients remains low. The remarkable phenotypic plasticity of melanoma cells allows for rapid development of invasive properties and metastatic tumors, the main cause of mortality in melanoma patients. Phenotypic and molecular analyses of developing tumors revealed that epithelial-mesenchymal transition (EMT), a cellular and molecular mechanism, controls transition from mature melanocyte to less differentiated melanocyte lineage progenitor cells forming melanoma tumors. This transition is facilitated by persistence of transcriptional regulatory circuit characteristic of embryonic stage in mature melanocytes. Switching of the developmental program of mature melanocyte to EMT is induced by accumulated mutations, especially targeting BRAF, N-RAS, or MEK1/2 signaling pathways, and further promoted by dynamic stimuli from local environment including hypoxia, interactions with extracellular matrix and growth factors or cytokines. Recent reports demonstrate that signaling mediated by transforming growth factor-β (TGF-β) and bone morphogenic proteins (BMPs) play critical roles in inducing EMT by controlling expression of critical transcription factors. BMPs are essential modulators of differentiation, proliferation, apoptosis, invasiveness, and metastases in developing melanoma tumors. They control transcription and epigenetic landscape of melanoma cells. Better understanding of the role of BMPs may lead to new strategies to control EMT processes in melanocyte cell lineage and to achieve clinical benefits for the patients.

Reeling in the Zebrafish Cancer Models

Zebrafish are rapidly becoming a leading model organism for cancer research. The genetic pathways driving cancer are highly conserved between zebrafish and humans, and the ability to easily manipulate the zebrafish genome to rapidly generate transgenic animals makes zebrafish an excellent model organism. Transgenic zebrafish containing complex, patient-relevant genotypes have been used to model many cancer types. Here we present a comprehensive review of transgenic zebrafish cancer models as a resource to the field and highlight important areas of cancer biology that have yet to be studied in the fish. The ability to image cancer cells and niche biology in an endogenous tumor makes zebrafish an indispensable model organism in which we can further understand the mechanisms that drive tumorigenesis and screen for potential new cancer therapies.

SATB2 induction of a neural crest mesenchyme-like program drives melanoma invasion and drug resistance

Recent genomic and scRNA-seq analyses of melanoma demonstrated a lack of recurrent genetic drivers of metastasis, while identifying common transcriptional states correlating with invasion or drug resistance. To test whether transcriptional adaptation can drive melanoma progression, we made use of a zebrafish mitfa:BRAFV600E;tp53-/- model, in which malignant progression is characterized by minimal genetic evolution. We undertook an overexpression-screen of 80 epigenetic/transcriptional regulators and found neural crest-mesenchyme developmental regulator SATB2 to accelerate aggressive melanoma development. Its overexpression induces invadopodia formation and invasion in zebrafish tumors and human melanoma cell lines. SATB2 binds and activates neural crest-regulators, including pdgfab and snai2. The transcriptional program induced by SATB2 overlaps with known MITF^lowAXL^high and AQP1⁺NGFR1^high drug-resistant states and functionally drives enhanced tumor propagation and resistance to Vemurafenib in vivo. In summary, we show that melanoma transcriptional rewiring by SATB2 to a neural crest mesenchyme-like program can drive invasion and drug resistance in autochthonous tumors.

XGSEA: CROSS-species gene set enrichment analysis via domain adaptation

MOTIVATION

Gene set enrichment analysis (GSEA) has been widely used to identify gene sets with statistically significant difference between cases and controls against a large gene set. GSEA needs both phenotype labels and expression of genes. However, gene expression are assessed more often for model organisms than minor species. Also, importantly gene expression are not measured well under specific conditions for human, due to high risk of direct experiments, such as non-approved treatment or gene knockout, and then often substituted by mouse. Thus, predicting enrichment significance (on a phenotype) of a given gene set of a species (target, say human), by using gene expression measured under the same phenotype of the other species (source, say mouse) is a vital and challenging problem, which we call CROSS-species gene set enrichment problem (XGSEP).

RESULTS

For XGSEP, we propose the CROSS-species gene set enrichment analysis (XGSEA), with three steps of: (1) running GSEA for a source species to obtain enrichment scores and $p$-values of source gene sets; (2) representing the relation between source and target gene sets by domain adaptation; and (3) using regression to predict $p$-values of target gene sets, based on the representation in (2). We extensively validated the XGSEA by using five regression and one classification measurements on four real data sets under various settings, proving that the XGSEA significantly outperformed three baseline methods in most cases. A case study of identifying important human pathways for T -cell dysfunction and reprogramming from mouse ATAC-Seq data further confirmed the reliability of the XGSEA.

AVAILABILITY

Source code of the XGSEA is available through https://github.com/LiminLi-xjtu/XGSEA.

Deciphering Melanoma Cell States and Plasticity with Zebrafish Models

Dynamic cellular heterogeneity underlies melanoma progression and therapy resistance. Advances in single-cell technologies have revealed an increasing number of tumor and microenvironment cell states in melanoma, but little is understood about their function in vivo. Zebrafish models are a powerful system for discovery, live imaging, and functional investigation of cell states throughout melanoma progression and treatment. By capturing dynamic melanoma states in living animals, zebrafish have the potential to resolve the complexity of melanoma heterogeneity from a single cell through disease processes within the context of the whole body, revealing novel cancer biology and therapeutic targets.

Making a melanoma: Molecular and cellular changes underlying melanoma initiation

Melanoma arises from the melanocyte lineage and is the most aggressive and lethal form of skin cancer. There are several genetic, genomic, and cellular changes associated with melanoma initiation. Here, we discuss these alterations and the melanoma cells of origin in which they are proposed to promote melanomagenesis.

BMP signaling: at the gate between activated melanocyte stem cells and differentiation

Through recurrent bouts synchronous with the hair cycle, quiescent melanocyte stem cells (McSCs) become activated to generate proliferative progeny that differentiate into pigment-producing melanocytes. The signaling factors orchestrating these events remain incompletely understood. Here, we use single-cell RNA sequencing with comparative gene expression analysis to elucidate the transcriptional dynamics of McSCs through quiescence, activation, and melanocyte maturation. Unearthing converging signs of increased WNT and BMP signaling along this progression, we endeavored to understand how these pathways are integrated. Employing conditional lineage-specific genetic ablation studies in mice, we found that loss of BMP signaling in the lineage leads to hair graying due to a block in melanocyte maturation. We show that interestingly, BMP signaling functions downstream from activated McSCs and maintains WNT effector, transcription factor LEF1. Employing pseudotime analysis, genetics, and chromatin landscaping, we show that following WNT-mediated activation of McSCs, BMP and WNT pathways collaborate to trigger the commitment of proliferative progeny by fueling LEF1- and MITF-dependent differentiation. Our findings shed light upon the signaling interplay and timing of cues that orchestrate melanocyte lineage progression in the hair follicle and underscore a key role for BMP signaling in driving complete differentiation.

Loss of prdm1a accelerates melanoma onset and progression

Melanoma is an aggressive, deadly skin cancer derived from melanocytes, a neural crest cell derivative. Melanoma cells mirror the developmental program of neural crest cells in that they exhibit the same gene expression patterns and utilize similar cellular mechanisms, including increased cell proliferation, epithelial-mesenchymal transition, and migration. Here we studied the role of neural crest regulator PRDM1 in melanoma onset and progression. In development, Prdm1a functions to promote neural crest progenitor fate, and in melanoma, we found that PRDM1 has reduced copy number and is recurrently deleted in both zebrafish and humans. When examining expression of neural crest and melanocyte development genes, we show that sox10 progenitor expression is high in prdm1a-/- mutants, while more differentiated melanocyte markers are reduced, suggesting that normally Prdm1a is required for differentiation. Data mining of human melanoma datasets indicates that high PRDM1 expression in human melanoma is correlated with better patient survival and decreased PRDM1 expression is common in metastatic tumors. When one copy of prdm1a is lost in the zebrafish melanoma model Tg(mitfa:BRAFV600E );p53-/- ;prdm1a+/- , melanoma onset occurs more quickly, and the tumors that form have a larger area with increased expression of sox10. These data demonstrate a novel role for PRDM1 as a tumor suppressor in melanoma.

From Tank to Treatment: Modeling Melanoma in Zebrafish

Melanoma is the deadliest form of skin cancer and one of few cancers with a growing incidence. A thorough understanding of its pathogenesis is fundamental to developing new strategies to combat mortality and morbidity. Zebrafish—due in large part to their tractable genetics, conserved pathways, and optical properties—have emerged as an excellent system to model melanoma. Zebrafish have been used to study melanoma from a single tumor initiating cell, through metastasis, remission, and finally into relapse. In this review, we examine seminal zebrafish studies that have advanced our understanding of melanoma.

Regulation of zebrafish melanocyte development by ligand-dependent BMP signaling

Preventing terminal differentiation is important in the development and progression of many cancers including melanoma. Recent identification of the BMP ligand GDF6 as a novel melanoma oncogene showed GDF6-activated BMP signaling suppresses differentiation of melanoma cells. Previous studies have identified roles for GDF6 orthologs during early embryonic and neural crest development, but have not identified direct regulation of melanocyte development by GDF6. Here, we investigate the BMP ligand gdf6a, a zebrafish ortholog of human GDF6, during the development of melanocytes from the neural crest. We establish that the loss of gdf6a or inhibition of BMP signaling during neural crest development disrupts normal pigment cell development, leading to an increase in the number of melanocytes and a corresponding decrease in iridophores, another neural crest-derived pigment cell type in zebrafish. This shift occurs as pigment cells arise from the neural crest and depends on mitfa, an ortholog of MITF, a key regulator of melanocyte development that is also targeted by oncogenic BMP signaling. Together, these results indicate that the oncogenic role ligand-dependent BMP signaling plays in suppressing differentiation in melanoma is a reiteration of its physiological roles during melanocyte development.

Embryonic bone morphogenetic protein and nodal induce invasion in melanocytes and melanoma cells

Despite recent progress in melanoma therapy via inhibition of activated oncogenes or immune stimulation, most stage IV melanoma patients still have limited survival times. Existing therapeutic approaches eventually fail to prevent further invasion and metastasis, which is driven by a morphological process termed epithelial-mesenchymal transition (EMT). We previously demonstrated that inhibition of EMT in melanoma cells via antagonizing the bone morphogenetic protein (BMP)-pathway abrogated EMT and neural crest migration of melanoma cells in chick embryos. Here, we show that BMP-2 is highly expressed in invasive melanoma cells and is elevated in the serum of stage IV melanoma patients compared to stage IB-IIC patients and healthy controls. Highly BMP-2-expressing melanoma cells display enhanced invasion in the rhombencephalon of the chick embryo. In addition to driving neural crest migration in the zebrafish embryo, the agonists BMP-2, BMP-7 and nodal induce EMT/invasion in radial growth phase melanoma cells and in human melanocytes in skin reconstructs. Blocking either BMP or nodal signaling by antagonists (noggin, lefty), or the Alk4/5/7-receptor inhibitor SB431542, decreases EMT and invasion of melanoma cells in human epidermal skin reconstructs. Together, our data suggest that inhibition of EMT-inducing pathways in melanoma might be a therapeutic approach to attenuate melanoma cell invasiveness.

Summary: We show that bone morphogenetic protein and nodal drive epithelial-mesenchymal transition (EMT) and invasiveness in melanoma cells, induce EMT and a melanoma-like invasive phenotype in melanocytes.