Understanding tumor heterogeneity is an important challenge in current cancer research. Transcription and epigenetic profiling of cultured melanoma cells have defined at least two distinct cell phenotypes characterized by distinctive gene expression signatures associated with high or low/absent expression of microphthalmia-associated transcription factor (MITF). Nevertheless, heterogeneity of cell populations and gene expression in primary human tumors is much less well characterized. We performed single-cell gene expression analyses on 472 cells isolated from needle biopsies of 5 primary human melanomas, 4 superficial spreading, and one acral melanoma. The expression of -high and-low signature genes was assessed and compared to investigate intra- and intertumoral heterogeneity and correlated gene expression profiles. Single-cell gene expression analyses revealed varying degrees of intra- and intertumor heterogeneity conferred by the variable expression of distinct sets of genes in different tumors. Expression of partially correlated with that of its known target genes, while expression correlated best with and Nevertheless, cells simultaneously expressing -high and-low signature genes were observed both by single-cell analyses and RNAscope. Single-cell analyses can be performed on limiting numbers of cells from primary human melanomas revealing their heterogeneity. Although tumors comprised variable proportions of cells with the -high and-low gene expression signatures characteristic of melanoma cultures, primary tumors also comprised cells expressing markers of both signatures defining a novel cell state in tumors .

MIcrophthalmia-associated Transcription Factor (MITF) regulates melanocyte and melanoma physiology. We show that MITF associates the NURF chromatin-remodelling factor in melanoma cells. ShRNA-mediated silencing of the NURF subunit BPTF revealed its essential role in several melanoma cell lines and in untransformed melanocytes in vitro. Comparative RNA-seq shows that MITF and BPTF co-regulate overlapping gene expression programs in cell lines in vitro. Somatic and specific inactivation of Bptf in developing murine melanoblasts in vivo shows that Bptf regulates their proliferation, migration and morphology. Once born, Bptf-mutant mice display premature greying where the second post-natal coat is white. This second coat is normally pigmented by differentiated melanocytes derived from the adult melanocyte stem cell (MSC) population that is stimulated to proliferate and differentiate at anagen. An MSC population is established and maintained throughout the life of the Bptf-mutant mice, but these MSCs are abnormal and at anagen, give rise to reduced numbers of transient amplifying cells (TACs) that do not express melanocyte markers and fail to differentiate into mature melanin producing melanocytes. MSCs display a transcriptionally repressed chromatin state and Bptf is essential for reactivation of the melanocyte gene expression program at anagen, the subsequent normal proliferation of TACs and their differentiation into mature melanocytes.

Somatic oncogenic mutation of BRAF coupled with inactivation of PTEN constitute a frequent combination of genomic alterations driving the development of human melanoma. Mice genetically engineered to conditionally express oncogenic Braf V600E and inactivate Pten in melanocytes following tamoxifen treatment rapidly develop melanoma. While early-stage melanomas comprised melanin-pigmented Mitf and Dct-expressing cells, expression of these and other melanocyte identity genes was lost in later stage tumours that showed histological and molecular characteristics of de-differentiated neural crest type cells. Melanocyte identity genes displayed loss of active chromatin marks and RNA polymerase II and gain of heterochromatin marks, indicating epigenetic reprogramming during tumour progression. Nevertheless, late-stage tumour cells grown in culture re-expressed Mitf, and melanocyte markers and Mitf together with Sox10 coregulated a large number of genes essential for their growth. In this melanoma model, somatic inactivation that the catalytic Brg1 (Smarca4) subunit of the SWI/SNF complex and the scaffolding Bptf subunit of the NuRF complex delayed tumour formation and deregulated large and overlapping gene expression programs essential for normal tumour cell growth. Moreover, we show that Brg1 and Bptf coregulated many genes together with Mitf and Sox10. Together these transcription factors and chromatin remodelling complexes orchestrate essential gene expression programs in mouse melanoma cells.

Chromodomain helicase DNA binding protein 4 (CHD4) is an ATPase subunit of the Nucleosome Remodelling and Deacetylation (NuRD) complex that regulates gene expression. CHD4 is essential for growth of multiple patient derived melanoma xenografts and for breast cancer. Here we show that CHD4 regulates expression of PADI1 (Protein Arginine Deiminase 1) and PADI3 in multiple cancer cell types modulating citrullination of arginine residues of the allosterically-regulated glycolytic enzyme pyruvate kinase M2 (PKM2). Citrullination of PKM2 R106 reprogrammes cross-talk between PKM2 ligands lowering its sensitivity to the inhibitors Tryptophan, Alanine and Phenylalanine and promoting activation by Serine. Citrullination thus bypasses normal physiological regulation by low Serine levels to promote excessive glycolysis and reduced cell proliferation. We further show that PADI1 and PADI3 expression is up-regulated by hypoxia where PKM2 citrullination contributes to increased glycolysis. We provide insight as to how conversion of arginines to citrulline impacts key interactions within PKM2 that act in concert to reprogramme its activity as an additional mechanism regulating this important enzyme.

In melanoma, initiating oncogenic mutations in BRAF or NRAS are detected in premalignant lesions that accumulate additional mutations and genomic instability as the tumor evolves to the metastatic state. Here we investigate evolution of clonal composition and neoantigen landscape in an atypical melanoma displaying recurrent cutaneous lesions over a 6-year period without development of extracutaneous metastases. Whole exome sequencing of four cutaneous lesions taken during the 6-year period identified a collection of single nucleotide variants and small insertions and deletions shared among all tumors, along with progressive selection of subclones displaying fewer single nucleotide variants. Later tumors also displayed lower neoantigen burden compared to early tumors, suggesting that clonal evolution was driven, at least in part, by counter selection of subclones with high neoantigen burdens. Among the selected mutations are a missense mutation in MAP2K1 (F53Y) and an inversion on chromosome 7 generating a AKAP9-BRAF fusion. The mutant proteins cooperatively activate the MAPK signaling pathway confirming they are potential driver mutations of this tumor. We therefore describe the long-term genetic evolution of cutaneous metastatic melanoma characterized by an unexpected phenotypic stability and neoantigen-driven clonal selection.

CHD3 and CHD4 are mutually exclusive ATPase subunits of the Nucleosome Remodelling and Deacetylation (NuRD) complex that regulates gene expression. CHD4 is essential for growth of multiple patient derived melanoma xenografts and for breast cancer. Here we show that CHD4 regulates expression of PADI1 (Protein Arginine Deiminase 1) and PADI3 in multiple cancer cell types modulating citrullination of three arginines of the allostericallyregulated glycolytic enzyme pyruvate kinase M2 (PKM2). PKM2 citrullination lowers sensitivity to the inhibitors Tryptophan, Alanine and Phenylalanine shifting the equilibrium towards the activator Serine. Citrullination thus bypasses normal physiological regulation by low Serine levels to promote excessive glycolysis defining a novel pathway regulating proliferation of melanoma and other cancer cells. We provide unique insight as to how conversion of arginines to citrulline impacts key interactions within PKM2 reprogramming its activity as an additional layer of complexity to the mechanisms regulating this important enzyme.Significance. This study describes the mechanism by which citrullination of key arginines in pyruvate kinase M2 reprograms its allosteric regulation by the activating amino acid serine.PKM2 citrullination defines a novel pathway regulating glycolysis and cancer cell proliferation.3 Introduction.

Conversion of peptidyl-arginine to peptidyl citrulline, known as citrullination, is a posttranslational protein modification catalysed by the PADI (Protein Arginine Deiminase) family of enzymes. PADI1 and PADI3 catalyse citrullination of arginine 106 in the glycolytic enzyme pyruvate kinase M2 modulating its allosteric regulation, glycolysis and cancer cell proliferation.

Malignant melanoma is a highly aggressive cancer with a propensity for early metastasis. Melanocyte transformation results predominantly from oncogenic mutations in BRAF, NRAS or NF1 leading to constitutive activation of the MAP kinase pathway driving cell proliferation and second site mutations such as loss of CDKN1A, or PTEN or activating mutations in the beta-catenin pathway that allow escape from oncogene induced senescence. Nevertheless, irrespective of the nature of the driver mutations, melanoma cell physiology is strongly regulated by transcription factors and epigenetic mechanisms. MITF (Microphthalmia-associated Transcription Factor) and SOX10 are two major transcription factors that regulate both normal melanocyte and melanoma cell physiology. Using a combination of mouse genetics, biochemistry and high throughput genomics we have identified cofactors for MITF and addressed the mechanisms by which MITF, SOX10 and their cofactors regulate gene expression in melanocytes and melanoma.