Direct reprogramming provides a fundamentally new approach for the generation of patient-specific cells. Here, by screening a pool of candidate transcription factors, we identify that a combination of the three factors, MITF, SOX10 and PAX3, directly converts mouse and human fibroblasts to functional melanocytes. Induced melanocytes (iMels) activate melanocyte-specific networks, express components of pigment production and delivery system and produce melanosomes. Human iMels properly integrate into the dermal-epidermal junction and produce and deliver melanin pigment to surrounding keratinocytes in a 3D organotypic skin reconstruct. Human iMels generate pigmented epidermis and hair follicles in skin reconstitution assays in vivo. The generation of iMels has important implications for studies of melanocyte lineage commitment, pigmentation disorders and cell replacement therapies.

Tumor associated inflammation predicts response to immune checkpoint blockade in human melanoma. Current theories on regulation of inflammation center on anti-tumor T cell responses. Here we show that tumor associated B cells are vital to melanoma associated inflammation. Human B cells express pro- and anti-inflammatory factors and differentiate into plasmablast-like cells when exposed to autologous melanoma secretomes in vitro. This plasmablast-like phenotype can be reconciled in human melanomas where plasmablast-like cells also express T cell-recruiting chemokines CCL3, CCL4, CCL5. Depletion of B cells in melanoma patients by anti-CD20 immunotherapy decreases tumor associated inflammation and CD8+ T cell numbers. Plasmablast-like cells also increase PD-1+ T cell activation through anti-PD-1 blockade in vitro and their frequency in pretherapy melanomas predicts response and survival to immune checkpoint blockade. Tumor associated B cells therefore orchestrate and sustain melanoma inflammation and may represent a predictor for survival and response to immune checkpoint blockade therapy.

The lack of tumor infiltration by CD8+ T cells is associated with poor patient response to anti-PD-1 therapy. Understanding how tumor infiltration is regulated is key to improving treatment efficacy. Here, we report that phosphorylation of HRS, a pivotal component of the ESCRT complex involved in exosome biogenesis, restricts tumor infiltration of cytolytic CD8+ T cells. Following ERK-mediated phosphorylation, HRS interacts with and mediates the selective loading of PD-L1 to exosomes, which inhibits the migration of CD8+ T cells into tumors. In tissue samples from patients with melanoma, CD8+ T cells are excluded from the regions where tumor cells contain high levels of phosphorylated HRS. In murine tumor models, overexpression of phosphorylated HRS increases resistance to anti-PD-1 treatment, whereas inhibition of HRS phosphorylation enhances treatment efficacy. Our study reveals a mechanism by which phosphorylation of HRS in tumor cells regulates anti-tumor immunity by inducing PD-L1+ immunosuppressive exosomes, and suggests HRS phosphorylation blockade as a potential strategy to improve the efficacy of cancer immunotherapy.

Anti-PD-1 therapy is used as a front-line treatment for many cancers, but mechanistic insight into this therapy resistance is still lacking. Here we generate a humanized (Hu)-mouse melanoma model by injecting fetal liver-derived CD34+ cells and implanting autologous thymus in immune-deficient NOD-scid IL2Rγnull (NSG) mice. Reconstituted Hu-mice are challenged with HLA-matched melanomas and treated with anti-PD-1, which results in restricted tumor growth but not complete regression. Tumor RNA-seq, multiplexed imaging and immunohistology staining show high expression of chemokines, as well as recruitment of FOXP3+ Treg and mast cells, in selective tumor regions. Reduced HLA-class I expression and CD8+/Granz B+ T cells homeostasis are observed in tumor regions where FOXP3+ Treg and mast cells co-localize, with such features associated with resistance to anti-PD-1 treatment. Combining anti-PD-1 with sunitinib or imatinib results in the depletion of mast cells and complete regression of tumors. Our results thus implicate mast cell depletion for improving the efficacy of anti-PD-1 therapy.

Melanoma cells display distinct intrinsic phenotypic states. Here, we seek to characterize the molecular regulation of these states using multi-omic analyses of whole exome, transcriptome, microRNA, long non-coding RNA and DNA methylation data together with reverse-phase protein array data on a panel of 68 highly annotated early passage melanoma cell lines. We demonstrate that clearly defined cancer cell intrinsic transcriptomic programs are maintained in melanoma cells ex vivo and remain highly conserved within melanoma tumors, are associated with distinct immune features within tumors, and differentially correlate with checkpoint inhibitor and adoptive T cell therapy efficacy. Through integrative analyses we demonstrate highly complex multi-omic regulation of melanoma cell intrinsic programs that provide key insights into the molecular maintenance of phenotypic states. These findings have implications for cancer biology and the identification of new therapeutic strategies. Further, these deeply characterized cell lines will serve as an invaluable resource for future research in the field.

Melanoma is a highly plastic tumor characterized by dynamic interconversion of different cell identities depending on the biological context. Melanoma cells with high expression of the H3K4 demethylase KDM5B (JARID1B) rest in a slow-cycling, yet reversible persister state. Over time, KDM5Bhigh cells can promote rapid tumor repopulation with equilibrated KDM5B expression heterogeneity. The cellular identity of KDM5Bhigh persister cells has not been studied so far, missing an important cell state-directed treatment opportunity in melanoma. Here, we have established a doxycycline-titratable system for genetic induction of permanent intratumor expression of KDM5B and screened for chemical agents that phenocopy this effect. Transcriptional profiling and cell functional assays confirmed that the dihydropyridine 2-phenoxyethyl 4-(2-fluorophenyl)-2,7,7-trimethyl-5-oxo-1,4,5,6,7,8-hexa-hydro-quinoline-3-carboxylate (termed Cpd1) supports high KDM5B expression and directs melanoma cells towards differentiation along the melanocytic lineage and to cell cycle-arrest. The high KDM5B state additionally prevents cell proliferation through negative regulation of cytokinetic abscission. Moreover, treatment with Cpd1 promoted the expression of the melanocyte-specific tyrosinase gene specifically sensitizing melanoma cells for the tyrosinase-processed antifolate prodrug 3-O-(3,4,5-trimethoxybenzoyl)-(-)-epicatechin (TMECG). In summary, our study provides proof-of-concept for a dual hit strategy in melanoma, in which persister state-directed transitioning limits tumor plasticity and primes melanoma cells towards lineage-specific elimination.

Protein levels and function are poorly predicted by genomic and transcriptomic analysis of patient tumours. Therefore, direct study of the functional proteome has the potential to provide a wealth of information that complements and extends genomic, epigenomic and transcriptomic analysis in The Cancer Genome Atlas (TCGA) projects. Here we use reverse-phase protein arrays to analyse 3,467 patient samples from 11 TCGA 'Pan-Cancer' diseases, using 181 high-quality antibodies that target 128 total proteins and 53 post-translationally modified proteins. The resultant proteomic data are integrated with genomic and transcriptomic analyses of the same samples to identify commonalities, differences, emergent pathways and network biology within and across tumour lineages. In addition, tissue-specific signals are reduced computationally to enhance biomarker and target discovery spanning multiple tumour lineages. This integrative analysis, with an emphasis on pathways and potentially actionable proteins, provides a framework for determining the prognostic, predictive and therapeutic relevance of the functional proteome.

To realize the full potential of immunotherapy, it is critical to understand the drivers of tumor infiltration by immune cells. Previous studies have linked immune infiltration with tumor neoantigen levels, but the broad applicability of this concept remains unknown. Here, we find that while this observation is true across cancers characterized by recurrent mutations, it does not hold for cancers driven by recurrent copy number alterations, such as breast and pancreatic tumors. To understand immune invasion in these cancers, we developed an integrative multi-omics framework, identifying the DNA damage response protein ATM as a driver of cytokine production leading to increased immune infiltration. This prediction was validated in numerous orthogonal datasets, as well as experimentally in vitro and in vivo by cytokine release and immune cell migration. These findings demonstrate diverse drivers of immune cell infiltration across cancer lineages and may facilitate the clinical adaption of immunotherapies across diverse malignancies.

Gastric cancer is a heterogeneous cancer, making treatment responses difficult to predict. Here we show that we identify two distinct molecular subtypes, mesenchymal phenotype (MP) and epithelial phenotype (EP), by analyzing genomic and proteomic data. Molecularly, MP subtype tumors show high genomic integrity characterized by low mutation rates and microsatellite stability, whereas EP subtype tumors show low genomic integrity. Clinically, the MP subtype is associated with markedly poor survival and resistance to standard chemotherapy, whereas the EP subtype is associated with better survival rates and sensitivity to chemotherapy. Integrative analysis shows that signaling pathways driving epithelial-to-mesenchymal transition and insulin-like growth factor 1 (IGF1)/IGF1 receptor (IGF1R) pathway are highly activated in MP subtype tumors. Importantly, MP subtype cancer cells are more sensitive to inhibition of IGF1/IGF1R pathway than EP subtype. Detailed characterization of these two subtypes could identify novel therapeutic targets and useful biomarkers for prognosis and therapy response.

Development of candidate cancer treatments is a resource-intensive process, with the research community continuing to investigate options beyond static genomic characterization. Toward this goal, we have established the genomic landscapes of 536 patient-derived xenograft (PDX) models across 25 cancer types, together with mutation, copy number, fusion, transcriptomic profiles, and NCI-MATCH arms. Compared with human tumors, PDXs typically have higher purity and fit to investigate dynamic driver events and molecular properties via multiple time points from same case PDXs. Here, we report on dynamic genomic landscapes and pharmacogenomic associations, including associations between activating oncogenic events and drugs, correlations between whole-genome duplications and subclone events, and the potential PDX models for NCI-MATCH trials. Lastly, we provide a web portal having comprehensive pan-cancer PDX genomic profiles and source code to facilitate identification of more druggable events and further insights into PDXs' recapitulation of human tumors.

Both genomic and transcriptomic signatures have been developed to predict responses of metastatic melanoma to immune checkpoint blockade (ICB) therapies; however, most of these signatures are derived from pre-treatment biopsy samples. Here, we build pathway-based super signatures in pre-treatment (PASS-PRE) and on-treatment (PASS-ON) tumor specimens based on transcriptomic data and clinical information from a large dataset of metastatic melanoma treated with anti-PD1-based therapies as the training set. Both PASS-PRE and PASS-ON signatures are validated in three independent datasets of metastatic melanoma as the validation set, achieving area under the curve (AUC) values of 0.45-0.69 and 0.85-0.89, respectively. We also combine all test samples and obtain AUCs of 0.65 and 0.88 for PASS-PRE and PASS-ON signatures, respectively. When compared with existing signatures, the PASS-ON signature demonstrates more robust and superior predictive performance across all four datasets. Overall, we provide a framework for building pathway-based signatures that is highly and accurately predictive of response to anti-PD1 therapies based on on-treatment tumor specimens. This work would provide a rationale for applying pathway-based signatures derived from on-treatment tumor samples to predict patients' therapeutic response to ICB therapies.

Primary liver cancer is a heterogeneous disease in terms of its etiology, histology, and therapeutic response. Concurrent proteomic and genomic characterization of a large set of clinical liver cancer samples can help elucidate the molecular basis of heterogeneity and thus serve as a valuable resource for personalized liver cancer treatment. In this study, we perform proteomic profiling of ~300 proteins on 259 primary liver cancer tissues with reverse-phase protein arrays, mutational analysis using whole genome sequencing and transcriptional analysis with RNA-Seq. Patients are of Japanese ethnic background and mainly HBV or HCV positive, providing insight into this important liver cancer subtype. Unsupervised classification of tumors based on protein expression profiles reveal three proteomic subclasses R1, R2, and R3. The R1 subclass is immunologically hot and demonstrated a good prognosis. R2 contains advanced proliferative tumor with TP53 mutations, high expression of VEGF receptor 2 and the worst prognosis. R3 is enriched with CTNNB1 mutations and elevated mTOR signaling pathway activity. Twenty-two proteins, including CDK1 and CDKN2A, are identified as potential prognostic markers. The proteomic classification presented in this study can help guide therapeutic decision making for liver cancer treatment.

The circadian clock regulator Bmal1 modulates tumorigenesis, but its reported effects are inconsistent. Here, we show that Bmal1 has a context-dependent role in mouse melanoma tumor growth. Loss of Bmal1 in YUMM2.1 or B16-F10 melanoma cells eliminates clock function and diminishes hypoxic gene expression and tumorigenesis, which could be rescued by ectopic expression of HIF1α in YUMM2.1 cells. By contrast, over-expressed wild-type or a transcriptionally inactive mutant Bmal1 non-canonically sequester myosin heavy chain 9 (Myh9) to increase MRTF-SRF activity and AP-1 transcriptional signature, and shift YUMM2.1 cells from a Sox10high to a Sox9high immune resistant, mesenchymal cell state that is found in human melanomas. Our work describes a link between Bmal1, Myh9, mouse melanoma cell plasticity, and tumor immunity. This connection may underlie cancer therapeutic resistance and underpin the link between the circadian clock, MRTF-SRF and the cytoskeleton.

In melanoma, therapies with inhibitors to oncogenic BRAFV600E are highly effective but responses are often short-lived due to the emergence of drug-resistant tumor subpopulations. We describe here a mechanism of acquired drug resistance through the tumor microenvironment, which is mediated by human tumor-associated B cells. Human melanoma cells constitutively produce the growth factor FGF-2, which activates tumor-infiltrating B cells to produce the growth factor IGF-1. B-cell-derived IGF-1 is critical for resistance of melanomas to BRAF and MEK inhibitors due to emergence of heterogeneous subpopulations and activation of FGFR-3. Consistently, resistance of melanomas to BRAF and/or MEK inhibitors is associated with increased CD20 and IGF-1 transcript levels in tumors and IGF-1 expression in tumor-associated B cells. Furthermore, first clinical data from a pilot trial in therapy-resistant metastatic melanoma patients show anti-tumor activity through B-cell depletion by anti-CD20 antibody. Our findings establish a mechanism of acquired therapy resistance through tumor-associated B cells with important clinical implications.Resistance to BRAFV600E inhibitors often occurs in melanoma patients. Here, the authors describe a potential mechanism of acquired drug resistance mediated by tumor-associated B cells-derived IGF-1.

Copy number loss of PIK3R1 (p85α) most commonly occurs in ovarian cancer among all cancer types. Here we report that ovarian cancer cells manifest a spectrum of tumorigenic phenotypes upon knockdown of PIK3R1. PIK3R1 loss activates AKT and p110-independent JAK2/STAT3 signaling through inducing changes in the phosphorylation of the docking protein Gab2, thereby relieving the negative inhibition on AKT and promoting the assembly of JAK2/STAT3 signalosome, respectively. Additional mechanisms leading to AKT activation include enhanced p110α kinase activity and a decrease in PTEN level. PIK3R1 loss renders ovarian cancer cells vulnerable to inhibition of AKT or JAK2/STAT3. The combination of AKT and STAT3 inhibitors significantly increases the anti-tumor effect compared to single-agent treatments. Together, our findings provide a rationale for mechanism-based therapeutic approach that targets tumors with loss of PIK3R1.

Checkpoint recovery, the process that checkpoint-arrested cells with normal DNA repair capacity resume cell cycle progression, is essential for genome stability. However, the signaling network of the process has not been clearly defined. Here, we combine functional proteomics, mathematical modeling, and molecular biology to identify mTORC1, the nutrient signaling integrator, as the determinant for G2/M checkpoint recovery. Inhibition of the mTORC1 pathway delays mitotic entry after DNA damage through KDM4B-mediated regulation of CCNB1 and PLK1 transcription. Cells with hyper-mTORC1 activity caused by TSC2 depletion exhibit accelerated G2/M checkpoint recovery. Those Tsc2-null cells are sensitive to WEE1 inhibition in vitro and in vivo by driving unscheduled mitotic entry and inducing mitotic catastrophe. These results reveal that mTORC1 functions as a mediator between nutrition availability sensing and cell fate determination after DNA damage, suggesting that checkpoint inhibitors may be used to treat mTORC1-hyperactivated tumors such as those associated with tuberous sclerosis complex.

Prolyl hydroxylase domain protein 2 (PHD2) is a well-known master oxygen sensor. However, the role of PHD2 in tumor initiation remains controversial. We find that during the transition of human nevi to melanoma, the expression of PHD2 protein is significantly decreased and lower expression PHD2 in melanoma is associated with worse clinical outcome. Knockdown of PHD2 leads to elevated Akt phosphorylation in human melanocytes. Mice with conditional melanocyte-specific expression of Phd2lox/lox (Tyr::CreER;Phd2lox/lox) fail to develop pigmented lesions. However, deletion of Phd2 in combination with expression of BRafV600E in melanocytes (Tyr::CreER;Phd2lox/lox;BRafCA) leads to the development of melanoma with 100% penetrance and frequent lymph node metastasis. Analysis of tumor tissues using reverse phase protein arrays demonstrates that Phd2 deletion activates the AKT-mTOR-S6 signaling axis in the recovered tumors. These data indicate that PHD2 is capable of suppressing tumor initiation largely mediated through inhibiting of the Akt-mTOR signaling pathway in the melanocyte lineage.

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.

PIK3R2 encodes the p85β regulatory subunit of phosphatidylinositol 3-kinase and is frequently amplified in cancers. The signaling mechanism and therapeutic implication of p85β are poorly understood. Here we report that p85β upregulates the protein level of the receptor tyrosine kinase AXL to induce oncogenic signaling in ovarian cancer. p85β activates p110 activity and AKT-independent PDK1/SGK3 signaling to promote tumorigenic phenotypes, which are all abolished upon inhibition of AXL. At the molecular level, p85β alters the phosphorylation of TRIM2 (an E3 ligase) and optineurin (an autophagy receptor), which mediate the selective regulation of AXL by p85β, thereby disrupting the autophagic degradation of the AXL protein. Therapeutically, p85β expression renders ovarian cancer cells vulnerable to inhibitors of AXL, p110, or PDK1. Conversely, p85β-depleted cells are less sensitive to these inhibitors. Together, our findings provide a rationale for pharmacological blockade of the AXL signaling axis in PIK3R2-amplified ovarian cancer.