This is a PDF file of a peer-reviewed paper that has been accepted for publication. Although unedited, the content has been subjected to preliminary formatting. Nature is providing this early version of the typeset paper as a service to our authors and readers. The text and figures will undergo copyediting and a proof review before the paper is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers apply.

Adoptive cell transfer therapies (ACTs) with cytotoxic T cells that target melanocytic antigens can achieve remissions in patients with metastatic melanomas, but tumours frequently relapse. Hypotheses explaining the acquired resistance to ACTs include the selection of antigen-deficient tumour cell variants and the induction of T-cell tolerance. However, the lack of appropriate experimental melanoma models has so far impeded clear insights into the underlying mechanisms. Here we establish an effective ACT protocol in a genetically engineered mouse melanoma model that recapitulates tumour regression, remission and relapse as seen in patients. We report the unexpected observation that melanomas acquire ACT resistance through an inflammation-induced reversible loss of melanocytic antigens. In serial transplantation experiments, melanoma cells switch between a differentiated and a dedifferentiated phenotype in response to T-cell-driven inflammatory stimuli. We identified the proinflammatory cytokine tumour necrosis factor (TNF)-α as a crucial factor that directly caused reversible dedifferentiation of mouse and human melanoma cells. Tumour cells exposed to TNF-α were poorly recognized by T cells specific for melanocytic antigens, whereas recognition by T cells specific for non-melanocytic antigens was unaffected or even increased. Our results demonstrate that the phenotypic plasticity of melanoma cells in an inflammatory microenvironment contributes to tumour relapse after initially successful T-cell immunotherapy. On the basis of our work, we propose that future ACT protocols should simultaneously target melanocytic and non-melanocytic antigens to ensure broad recognition of both differentiated and dedifferentiated melanoma cells, and include strategies to sustain T-cell effector functions by blocking immune-inhibitory mechanisms in the tumour microenvironment.

Intermittent intense ultraviolet (UV) exposure represents an important aetiological factor in the development of malignant melanoma. The ability of UV radiation to cause tumour-initiating DNA mutations in melanocytes is now firmly established, but how the microenvironmental effects of UV radiation influence melanoma pathogenesis is not fully understood. Here we report that repetitive UV exposure of primary cutaneous melanomas in a genetically engineered mouse model promotes metastatic progression, independent of its tumour-initiating effects. UV irradiation enhanced the expansion of tumour cells along abluminal blood vessel surfaces and increased the number of lung metastases. This effect depended on the recruitment and activation of neutrophils, initiated by the release of high mobility group box 1 (HMGB1) from UV-damaged epidermal keratinocytes and driven by Toll-like receptor 4 (TLR4). The UV-induced neutrophilic inflammatory response stimulated angiogenesis and promoted the ability of melanoma cells to migrate towards endothelial cells and use selective motility cues on their surfaces. Our results not only reveal how UV irradiation of epidermal keratinocytes is sensed by the innate immune system, but also show that the resulting inflammatory response catalyses reciprocal melanoma-endothelial cell interactions leading to perivascular invasion, a phenomenon originally described as angiotropism in human melanomas by histopathologists. Angiotropism represents a hitherto underappreciated mechanism of metastasis that also increases the likelihood of intravasation and haematogenous dissemination. Consistent with our findings, ulcerated primary human melanomas with abundant neutrophils and reactive angiogenesis frequently show angiotropism and a high risk for metastases. Our work indicates that targeting the inflammation-induced phenotypic plasticity of melanoma cells and their association with endothelial cells represent rational strategies to specifically interfere with metastatic progression.

Malignant abdominal fluid (ascites) frequently develops in women with advanced high-grade serous ovarian cancer (HGSOC) and is associated with drug resistance and a poor prognosis 1 . To comprehensively characterize the HGSOC ascites ecosystem, we used single-cell RNA-seq (scRNA-seq) to profile ~11,000 cells from 22 ascites specimens from 11 HGSOC patients. We found significant inter-patient variability in the composition and functional programs of ascites cells, including immunomodulatory fibroblast sub-populations and dichotomous macrophage populations. We find that the previously described "immunoreactive" and "mesenchymal" subtypes of HGSOC, which have prognostic implications, reflect the abundance of immune infiltrates and fibroblasts rather than distinct subsets of malignant cells 2 . Malignant cell variability was partly explained by heterogeneous copy number alterations (CNA) patterns or expression of a stemness program. Malignant cells shared expression of inflammatory programs that were largely recapitulated in scRNA-seq of ~35,000 cells from additionally collected samples, including three ascites, two primary HGSOC tumors and three patient-ascites-derived xenograft models. Inhibition of the JAK/STAT-pathway, which was expressed in both malignant cells and CAFs, had potent anti-tumor activity in primary short-term cultures and PDX models. Our work contributes to resolving the HSGOC landscape 3-5 and provides a resource for the development of novel therapeutic approaches.

Inhibitors of the receptor tyrosine kinase c-MET are currently used in the clinic to target oncogenic signaling in tumor cells. We found that concomitant c-MET inhibition promoted adoptive T cell transfer and checkpoint immunotherapies in murine cancer models by increasing effector T cell infiltration in tumors. This therapeutic effect was independent of tumor cell-intrinsic c-MET dependence. Mechanistically, c-MET inhibition impaired the reactive mobilization and recruitment of neutrophils into tumors and draining lymph nodes in response to cytotoxic immunotherapies. In the absence of c-MET inhibition, neutrophils recruited to T cell-inflamed microenvironments rapidly acquired immunosuppressive properties, restraining T cell expansion and effector functions. In cancer patients, high serum levels of the c-MET ligand HGF correlated with increasing neutrophil counts and poor responses to checkpoint blockade therapies. Our findings reveal a role for the HGF/c-MET pathway in neutrophil recruitment and function and suggest that c-MET inhibitor co-treatment may improve responses to cancer immunotherapy in settings beyond c-MET-dependent tumors.

Resistance to immune checkpoint inhibitors (ICI) that activate T cell mediated anti-tumor immunity is a key challenge in cancer therapy, yet the underlying mechanisms remain poorly understood. To further elucidate those, we developed a new approach, Perturb-CITE-seq, for pooled CRISPR perturbation screens with multi-modal RNA and protein single-cell profiling readout and applied it to screen patient-derived autologous melanoma and tumor infiltrating lymphocyte (TIL) co-cultures. We profiled RNA and 20 surface proteins in over 218,000 cells under ~750 perturbations, chosen by their membership in an immune evasion program that is associated with immunotherapy resistance in patients. Our screen recovered clinically-relevant resistance mechanisms concordantly reflected in RNA, protein and perturbation effects on susceptibility to T cell mediated killing. These were organized in eight co-functional modules whose perturbation distinctly affect four co-regulated programs associated with immune evasion. Among these were defects in the IFNγ-JAK/STAT pathway and in antigen presentation, and several novel mechanisms, including loss or downregulation of CD58, a surface protein without known mouse homolog. Leveraging the rich profiles in our screen, we found that loss of CD58 did not compromise MHC protein expression and that CD58 was not transcriptionally induced by the IFNγ pathway, allowing us to distinguish it as a novel mechanism of immune resistance.We further show that loss of CD58 on cancer cells conferred immune evasion across multiple T cell and Natural Killer cell patient co-culture models. Notably, CD58 is downregulated in tumors with resistance to immunotherapy in melanoma patients. Our work identifies novel mechanisms at the nexus of immune evasion and drug resistance and provides a general framework for deciphering complex mechanisms by large-scale perturbation screens with multi-modal singlecell profiles, including in systems consisting of multiple cell types.

We propose an extension of a standard stochastic individual-based model in population dynamics which broadens the range of biological applications. Our primary motivation is modelling of immunotherapy of malignant tumours. In this context the different actors, T-cells, cytokines or cancer cells, are modelled as single particles (individuals) in the stochastic system. The main expansions of the model are distinguishing cancer cells by phenotype and genotype, including environment-dependent phenotypic plasticity that does not affect the genotype, taking into account the effects of therapy and introducing a competition term which lowers the reproduction rate of an individual in addition to the usual term that increases its death rate. We illustrate the new setup by using it to model various phenomena arising in immunotherapy. Our aim is twofold: on the one hand, we show that the interplay of genetic mutations and phenotypic switches on different timescales as well as the occurrence of metastability phenomena raise new mathematical challenges. On the other hand, we argue why understanding purely stochastic events (which cannot be obtained with deterministic models) may help to understand the resistance of tumours to therapeutic approaches and may have non-trivial consequences on tumour treatment protocols. This is supported through numerical simulations.

No abstract