A TCF4-dependent gene regulatory network confers resistance to immunotherapy in melanoma.

To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.

Toward Minimal Residual Disease-Directed Therapy in Melanoma.

Many patients with advanced cancers achieve dramatic responses to a panoply of therapeutics yet retain minimal residual disease (MRD), which ultimately results in relapse. To gain insights into the biology of MRD, we applied single-cell RNA sequencing to malignant cells isolated from BRAF mutant patient-derived xenograft melanoma cohorts exposed to concurrent RAF/MEK-inhibition. We identified distinct drug-tolerant transcriptional states, varying combinations of which co-occurred within MRDs from PDXs and biopsies of patients on treatment. One of these exhibited a neural crest stem cell (NCSC) transcriptional program largely driven by the nuclear receptor RXRG. An RXR antagonist mitigated accumulation of NCSCs in MRD and delayed the development of resistance. These data identify NCSCs as key drivers of resistance and illustrate the therapeutic potential of MRD-directed therapy. They also highlight how gene regulatory network architecture reprogramming may be therapeutically exploited to limit cellular heterogeneity, a key driver of disease progression and therapy resistance.

Decoding the interplay between genetic and non-genetic drivers of metastasis.

Metastasis is a multistep process by which cancer cells break away from their original location and spread to distant organs, and is responsible for the vast majority of cancer-related deaths. Preventing early metastatic dissemination would revolutionize the ability to fight cancer. Unfortunately, the relatively poor understanding of the molecular underpinnings of metastasis has hampered the development of effective anti-metastatic drugs. Although it is now accepted that disseminating tumour cells need to acquire multiple competencies to face the many obstacles they encounter before reaching their metastatic site(s), whether these competencies are acquired through an accumulation of metastasis-specific genetic alterations and/or non-genetic events is often debated. Here we review a growing body of literature highlighting the importance of both genetic and non-genetic reprogramming events during the metastatic cascade, and discuss how genetic and non-genetic processes act in concert to confer metastatic competencies. We also describe how recent technological advances, and in particular the advent of single-cell multi-omics and barcoding approaches, will help to better elucidate the cross-talk between genetic and non-genetic mechanisms of metastasis and ultimately inform innovative paths for the early detection and interception of this lethal process.

CD4+ T cell-induced inflammatory cell death controls immune-evasive tumours.

Most clinically applied cancer immunotherapies rely on the ability of CD8+ cytolytic T cells to directly recognize and kill tumour cells1-3. These strategies are limited by the emergence of major histocompatibility complex (MHC)-deficient tumour cells and the formation of an immunosuppressive tumour microenvironment4-6. The ability of CD4+ effector cells to contribute to antitumour immunity independently of CD8+ T cells is increasingly recognized, but strategies to unleash their full potential remain to be identified7-10. Here, we describe a mechanism whereby a small number of CD4+ T cells is sufficient to eradicate MHC-deficient tumours that escape direct CD8+ T cell targeting. The CD4+ effector T cells preferentially cluster at tumour invasive margins where they interact with MHC-II+CD11c+ antigen-presenting cells. We show that T helper type 1 cell-directed CD4+ T cells and innate immune stimulation reprogramme the tumour-associated myeloid cell network towards interferon-activated antigen-presenting and iNOS-expressing tumouricidal effector phenotypes. Together, CD4+ T cells and tumouricidal myeloid cells orchestrate the induction of remote inflammatory cell death that indirectly eradicates interferon-unresponsive and MHC-deficient tumours. These results warrant the clinical exploitation of this ability of CD4+ T cells and innate immune stimulators in a strategy to complement the direct cytolytic activity of CD8+ T cells and natural killer cells and advance cancer immunotherapies.

A cellular hierarchy in melanoma uncouples growth and metastasis.

Although melanoma is notorious for its high degree of heterogeneity and plasticity1,2, the origin and magnitude of cell-state diversity remains poorly understood. Equally, it is unclear whether growth and metastatic dissemination are supported by overlapping or distinct melanoma subpopulations. Here, by combining mouse genetics, single-cell and spatial transcriptomics, lineage tracing and quantitative modelling, we provide evidence of a hierarchical model of tumour growth that mirrors the cellular and molecular logic underlying the cell-fate specification and differentiation of the embryonic neural crest. We show that tumorigenic competence is associated with a spatially localized perivascular niche, a phenotype acquired through an intercellular communication pathway established by endothelial cells. Consistent with a model in which only a fraction of cells are fated to fuel growth, temporal single-cell tracing of a population of melanoma cells with a mesenchymal-like state revealed that these cells do not contribute to primary tumour growth but, instead, constitute a pool of metastatic initiating cells that switch cell identity while disseminating to secondary organs. Our data provide a spatially and temporally resolved map of the diversity and trajectories of melanoma cell states and suggest that the ability to support growth and metastasis are limited to distinct pools of cells. The observation that these phenotypic competencies can be dynamically acquired after exposure to specific niche signals warrant the development of therapeutic strategies that interfere with the cancer cell reprogramming activity of such microenvironmental cues.

PHGDH heterogeneity potentiates cancer cell dissemination and metastasis.

Cancer metastasis requires the transient activation of cellular programs enabling dissemination and seeding in distant organs1. Genetic, transcriptional and translational heterogeneity contributes to this dynamic process2,3. Metabolic heterogeneity has also been observed4, yet its role in cancer progression is less explored. Here we find that the loss of phosphoglycerate dehydrogenase (PHGDH) potentiates metastatic dissemination. Specifically, we find that heterogeneous or low PHGDH expression in primary tumours of patients with breast cancer is associated with decreased metastasis-free survival time. In mice, circulating tumour cells and early metastatic lesions are enriched with Phgdhlow cancer cells, and silencing Phgdh in primary tumours increases metastasis formation. Mechanistically, Phgdh interacts with the glycolytic enzyme phosphofructokinase, and the loss of this interaction activates the hexosamine-sialic acid pathway, which provides precursors for protein glycosylation. As a consequence, aberrant protein glycosylation occurs, including increased sialylation of integrin αvß3, which potentiates cell migration and invasion. Inhibition of sialylation counteracts the metastatic ability of Phgdhlow cancer cells. In conclusion, although the catalytic activity of PHGDH supports cancer cell proliferation, low PHGDH protein expression non-catalytically potentiates cancer dissemination and metastasis formation. Thus, the presence of PHDGH heterogeneity in primary tumours could be considered a sign of tumour aggressiveness.

Targeted chemotherapy overcomes drug resistance in melanoma.

The emergence of drug resistance is a major obstacle for the success of targeted therapy in melanoma. Additionally, conventional chemotherapy has not been effective as drug-resistant cells escape lethal DNA damage effects by inducing growth arrest commonly referred to as cellular dormancy. We present a therapeutic strategy termed "targeted chemotherapy" by depleting protein phosphatase 2A (PP2A) or its inhibition using a small molecule inhibitor (1,10-phenanthroline-5,6-dione [phendione]) in drug-resistant melanoma. Targeted chemotherapy induces the DNA damage response without causing DNA breaks or allowing cellular dormancy. Phendione treatment reduces tumor growth of BRAFV600E-driven melanoma patient-derived xenografts (PDX) and diminishes growth of NRASQ61R-driven melanoma, a cancer with no effective therapy. Remarkably, phendione treatment inhibits the acquisition of resistance to BRAF inhibition in BRAFV600E PDX highlighting its effectiveness in combating the advent of drug resistance.

Melanoma plasticity and phenotypic diversity: therapeutic barriers and opportunities.

An incomplete view of the mechanisms that drive metastasis, the primary cause of cancer-related death, has been a major barrier to development of effective therapeutics and prognostic diagnostics. Increasing evidence indicates that the interplay between microenvironment, genetic lesions, and cellular plasticity drives the metastatic cascade and resistance to therapies. Here, using melanoma as a model, we outline the diversity and trajectories of cell states during metastatic dissemination and therapy exposure, and highlight how understanding the magnitude and dynamics of nongenetic reprogramming in space and time at single-cell resolution can be exploited to develop therapeutic strategies that capitalize on nongenetic tumor evolution.

LifeTime and improving European healthcare through cell-based interceptive medicine.

Here we describe the LifeTime Initiative, which aims to track, understand and target human cells during the onset and progression of complex diseases, and to analyse their response to therapy at single-cell resolution. This mission will be implemented through the development, integration and application of single-cell multi-omics and imaging, artificial intelligence and patient-derived experimental disease models during the progression from health to disease. The analysis of large molecular and clinical datasets will identify molecular mechanisms, create predictive computational models of disease progression, and reveal new drug targets and therapies. The timely detection and interception of disease embedded in an ethical and patient-centred vision will be achieved through interactions across academia, hospitals, patient associations, health data management systems and industry. The application of this strategy to key medical challenges in cancer, neurological and neuropsychiatric disorders, and infectious, chronic inflammatory and cardiovascular diseases at the single-cell level will usher in cell-based interceptive medicine in Europe over the next decade.

Mitochondrial MDM2 Regulates Respiratory Complex I Activity Independently of p53.

Accumulating evidence indicates that the MDM2 oncoprotein promotes tumorigenesis beyond its canonical negative effects on the p53 tumor suppressor, but these p53-independent functions remain poorly understood. Here, we show that a fraction of endogenous MDM2 is actively imported in mitochondria to control respiration and mitochondrial dynamics independently of p53. Mitochondrial MDM2 represses the transcription of NADH-dehydrogenase 6 (MT-ND6) in vitro and in vivo, impinging on respiratory complex I activity and enhancing mitochondrial ROS production. Recruitment of MDM2 to mitochondria increases during oxidative stress and hypoxia. Accordingly, mice lacking MDM2 in skeletal muscles exhibit higher MT-ND6 levels, enhanced complex I activity, and increased muscular endurance in mild hypoxic conditions. Furthermore, increased mitochondrial MDM2 levels enhance the migratory and invasive properties of cancer cells. Collectively, these data uncover a previously unsuspected function of the MDM2 oncoprotein in mitochondria that play critical roles in skeletal muscle physiology and may contribute to tumor progression.

Identification of the tumour transition states occurring during EMT.

In cancer, the epithelial-to-mesenchymal transition (EMT) is associated with tumour stemness, metastasis and resistance to therapy. It has recently been proposed that, rather than being a binary process, EMT occurs through distinct intermediate states. However, there is no direct in vivo evidence for this idea. Here we screen a large panel of cell surface markers in skin and mammary primary tumours, and identify the existence of multiple tumour subpopulations associated with different EMT stages: from epithelial to completely mesenchymal states, passing through intermediate hybrid states. Although all EMT subpopulations presented similar tumour-propagating cell capacity, they displayed differences in cellular plasticity, invasiveness and metastatic potential. Their transcriptional and epigenetic landscapes identify the underlying gene regulatory networks, transcription factors and signalling pathways that control these different EMT transition states. Finally, these tumour subpopulations are localized in different niches that differentially regulate EMT transition states.

Codon-specific translation reprogramming promotes resistance to targeted therapy.

Reprogramming of mRNA translation has a key role in cancer development and drug resistance 1 . However, the molecular mechanisms that are involved in this process remain poorly understood. Wobble tRNA modifications are required for specific codon decoding during translation2,3. Here we show, in humans, that the enzymes that catalyse modifications of wobble uridine 34 (U34) tRNA (U34 enzymes) are key players of the protein synthesis rewiring that is induced by the transformation driven by the BRAF V600E oncogene and by resistance to targeted therapy in melanoma. We show that BRAF V600E -expressing melanoma cells are dependent on U34 enzymes for survival, and that concurrent inhibition of MAPK signalling and ELP3 or CTU1 and/or CTU2 synergizes to kill melanoma cells. Activation of the PI3K signalling pathway, one of the most common mechanisms of acquired resistance to MAPK therapeutic agents, markedly increases the expression of U34 enzymes. Mechanistically, U34 enzymes promote glycolysis in melanoma cells through the direct, codon-dependent, regulation of the translation of HIF1A mRNA and the maintenance of high levels of HIF1α protein. Therefore, the acquired resistance to anti-BRAF therapy is associated with high levels of U34 enzymes and HIF1α. Together, these results demonstrate that U34 enzymes promote the survival and resistance to therapy of melanoma cells by regulating specific mRNA translation.

Chromatin-Bound MDM2 Regulates Serine Metabolism and Redox Homeostasis Independently of p53.

The mouse double minute 2 (MDM2) oncoprotein is recognized as a major negative regulator of the p53 tumor suppressor, but growing evidence indicates that its oncogenic activities extend beyond p53. Here, we show that MDM2 is recruited to chromatin independently of p53 to regulate a transcriptional program implicated in amino acid metabolism and redox homeostasis. Identification of MDM2 target genes at the whole-genome level highlights an important role for ATF3/4 transcription factors in tethering MDM2 to chromatin. MDM2 recruitment to chromatin is a tightly regulated process that occurs during oxidative stress and serine/glycine deprivation and is modulated by the pyruvate kinase M2 (PKM2) metabolic enzyme. Depletion of endogenous MDM2 in p53-deficient cells impairs serine/glycine metabolism, the NAD(+)/NADH ratio, and glutathione (GSH) recycling, impacting their redox state and tumorigenic potential. Collectively, our data illustrate a previously unsuspected function of chromatin-bound MDM2 in cancer cell metabolism.

Cross-species analysis of enhancer logic using deep learning.

Deciphering the genomic regulatory code of enhancers is a key challenge in biology because this code underlies cellular identity. A better understanding of how enhancers work will improve the interpretation of noncoding genome variation and empower the generation of cell type-specific drivers for gene therapy. Here, we explore the combination of deep learning and cross-species chromatin accessibility profiling to build explainable enhancer models. We apply this strategy to decipher the enhancer code in melanoma, a relevant case study owing to the presence of distinct melanoma cell states. We trained and validated a deep learning model, called DeepMEL, using chromatin accessibility data of 26 melanoma samples across six different species. We show the accuracy of DeepMEL predictions on the CAGI5 challenge, where it significantly outperforms existing models on the melanoma enhancer of IRF4 Next, we exploit DeepMEL to analyze enhancer architectures and identify accurate transcription factor binding sites for the core regulatory complexes in the two different melanoma states, with distinct roles for each transcription factor, in terms of nucleosome displacement or enhancer activation. Finally, DeepMEL identifies orthologous enhancers across distantly related species, where sequence alignment fails, and the model highlights specific nucleotide substitutions that underlie enhancer turnover. DeepMEL can be used from the Kipoi database to predict and optimize candidate enhancers and to prioritize enhancer mutations. In addition, our computational strategy can be applied to other cancer or normal cell types.

Long non-coding RNA TINCR suppresses metastatic melanoma dissemination by preventing ATF4 translation.

Transition from proliferative-to-invasive phenotypes promotes metastasis and therapy resistance in melanoma. Reversion of the invasive phenotype, however, is challenged by the poor understanding of mechanisms underlying its maintenance. Here, we report that the lncRNA TINCR is down-regulated in metastatic melanoma and its silencing increases the expression levels of invasive markers, in vitro migration, in vivo tumor growth, and resistance to BRAF and MEK inhibitors. The critical mediator is ATF4, a central player of the integrated stress response (ISR), which is activated in TINCR-depleted cells in the absence of starvation and eIF2α phosphorylation. TINCR depletion increases global protein synthesis and induces translational reprogramming, leading to increased translation of mRNAs encoding ATF4 and other ISR proteins. Strikingly, re-expression of TINCR in metastatic melanoma suppresses the invasive phenotype, reduces numbers of tumor-initiating cells and metastasis formation, and increases drug sensitivity. Mechanistically, TINCR interacts with mRNAs associated with the invasive phenotype, including ATF4, preventing their binding to ribosomes. Thus, TINCR is a suppressor of the melanoma invasive phenotype, which functions in nutrient-rich conditions by repressing translation of selected ISR RNAs.

Dynamic reversal of random X-Chromosome inactivation during iPSC reprogramming.

Induction and reversal of chromatin silencing is critical for successful development, tissue homeostasis, and the derivation of induced pluripotent stem cells (iPSCs). X-Chromosome inactivation (XCI) and reactivation (XCR) in female cells represent chromosome-wide transitions between active and inactive chromatin states. Although XCI has long been studied, providing important insights into gene regulation, the dynamics and mechanisms underlying the reversal of stable chromatin silencing of X-linked genes are much less understood. Here, we use allele-specific transcriptomics to study XCR during mouse iPSC reprogramming in order to elucidate the timing and mechanisms of chromosome-wide reversal of gene silencing. We show that XCR is hierarchical, with subsets of genes reactivating early, late, and very late during reprogramming. Early genes are activated before the onset of late pluripotency genes activation. Early genes are located genomically closer to genes that escape XCI, unlike genes reactivating late. Early genes also show increased pluripotency transcription factor (TF) binding. We also reveal that histone deacetylases (HDACs) restrict XCR in reprogramming intermediates and that the severe hypoacetylation state of the inactive X Chromosome (Xi) persists until late reprogramming stages. Altogether, these results reveal the timing of transcriptional activation of monoallelically repressed genes during iPSC reprogramming, and suggest that allelic activation involves the combined action of chromatin topology, pluripotency TFs, and chromatin regulators. These findings are important for our understanding of gene silencing, maintenance of cell identity, reprogramming, and disease.