Melanomas can have multiple coexisting cell states, including proliferative (PRO) versus invasive (INV) subpopulations that represent a "go or grow" trade-off; however, how these populations interact is poorly understood. Using a combination of zebrafish modeling and analysis of patient samples, we show that INV and PRO cells form spatially structured heterotypic clusters and cooperate in the seeding of metastasis, maintaining cell state heterogeneity. INV cells adhere tightly to each other and form clusters with a rim of PRO cells. Intravital imaging demonstrated cooperation in which INV cells facilitate dissemination of less metastatic PRO cells. We identified the TFAP2 neural crest transcription factor as a master regulator of clustering and PRO/INV states. Isolation of clusters from patients with metastatic melanoma revealed a subset with heterotypic PRO-INV clusters. Our data suggest a framework for the co-existence of these two divergent cell populations, in which heterotypic clusters promote metastasis via cell-cell cooperation.

Much of our understanding of the function of histone post-translational modifications in metazoans is inferred from their genomic localization and / or extrapolated from yeast studies. For example, acetylation of histone H3 lysine 56 (H3 K56Ac) is assumed to be important for transcriptional regulation in metazoan cells based on its occurrence at promoters and its function in yeast. Here we directly assess the function of H3 K56Ac during chromatin disassembly from gene regulatory regions during transcriptional induction in human cells by using mutations that either mimic or prevent H3 K56Ac. Although there is rapid histone H3 disassembly during induction of some estrogen receptor responsive genes, depletion of the histone chaperone ASF1A/B, which is required for H3 K56 acetylation, has no effect on chromatin disassembly at these regions. During the course of this work, we found that all the commercially available antibodies to H3 K56Ac are non-specific in human cells and in Drosophila. We used H3-YFP fusions to show that the H3 K56Q mutation can promote chromatin disassembly from regulatory regions of some estrogen responsive genes in the context of transcriptional induction. However, neither the H3 K56R nor K56Q mutation significantly altered chromatin disassembly dynamics by FRAP analysis. These results indicate that unlike the situation in yeast, human cells do not use H3 K56Ac to promote chromatin disassembly from regulatory regions or from the genome in general. Furthermore, our work highlights the need for rigorous characterization of the specificity of antibodies to histone post-translational modifications in vivo.

Melanoma exhibits numerous transcriptional cell states including neural crest-like cells as well as pigmented melanocytic cells. How these different cell states relate to distinct tumorigenic phenotypes remains unclear. Here, we use a zebrafish melanoma model to identify a transcriptional program linking the melanocytic cell state to a dependence on lipid droplets, the specialized organelle responsible for lipid storage. Single-cell RNA-sequencing of these tumors show a concordance between genes regulating pigmentation and those involved in lipid and oxidative metabolism. This state is conserved across human melanoma cell lines and patient tumors. This melanocytic state demonstrates increased fatty acid uptake, an increased number of lipid droplets, and dependence upon fatty acid oxidative metabolism. Genetic and pharmacologic suppression of lipid droplet production is sufficient to disrupt cell cycle progression and slow melanoma growth in vivo. Because the melanocytic cell state is linked to poor outcomes in patients, these data indicate a metabolic vulnerability in melanoma that depends on the lipid droplet organelle.

In melanoma, predicting which tumors will ultimately metastasize guides treatment decisions. Transcriptional signatures of primary tumors have been utilized to predict metastasis, but which among these are driver or passenger events remains unclear. We used data from the adjuvant AVAST-M trial to identify a predictive gene signature in localized tumors that ultimately metastasized. Using a zebrafish model of primary melanoma, we interrogated the top genes from the AVAST-M signature in vivo. This identified GRAMD1B, a cholesterol transfer protein, as a bona fide metastasis suppressor, with a majority of knockout animals rapidly developing metastasis. Mechanistically, excess free cholesterol or its metabolite 27-hydroxycholesterol promotes invasiveness via activation of an AP-1 program, which is associated with increased metastasis in humans. Our data demonstrate that the transcriptional seeds of metastasis are embedded within localized tumors, suggesting that early targeting of these programs can be used to prevent metastatic relapse.

Oncogenes only transform cells under certain cellular contexts, a phenomenon called oncogenic competence. Using a combination of a human pluripotent stem cell–derived cancer model along with zebrafish transgenesis, we demonstrate that the transforming ability of BRAFV600E along with additional mutations depends on the intrinsic transcriptional program present in the cell of origin. In both systems, melanocytes are less responsive to mutations, whereas both neural crest and melanoblast populations are readily transformed. Profiling reveals that progenitors have higher expression of chromatin-modifying enzymes such as ATAD2, a melanoma competence factor that forms a complex with SOX10 and allows for expression of downstream oncogenic and neural crest programs. These data suggest that oncogenic competence is mediated by regulation of developmental chromatin factors, which then allow for proper response to those oncogenes.

The access-repair-restore model for the role of chromatin in DNA repair infers that chromatin is a mere obstacle to DNA repair. However, here we show that blocking chromatin assembly, via knockdown of the histone chaperones ASF1 or CAF-1 or a mutation that prevents ASF1A binding to histones, hinders Rad51 loading onto ssDNA during homologous recombination. This is a consequence of reduced recruitment of the Rad51 loader MMS22L-TONSL to ssDNA, resulting in persistent RPA foci, extensive DNA end resection, persistent activation of the ATR-Chk1 pathway, and cell cycle arrest. In agreement, histones occupy ssDNA during DNA repair in yeast. We also uncovered DNA-PKcs-dependent DNA damage-induced ASF1A phosphorylation, which enhances chromatin assembly, promoting MMS22L-TONSL recruitment and, hence, Rad51 loading. We propose that transient assembly of newly synthesized histones onto ssDNA serves to recruit MMS22L-TONSL to efficiently form the Rad51 nucleofilament for strand invasion, suggesting an active role of chromatin assembly in homologous recombination.