Multimodel preclinical platform predicts clinical response of melanoma to immunotherapy

Stochastic modeling of single-cell gene expression adaptation reveals non-genomic contribution to evolution of tumor subclones

Cancer progression is an evolutionary process driven by the selection of cells adapted to gain growth advantage. We present a formal study on the adaptation of gene expression in subclonal evolution. We model evolutionary changes in gene expression as stochastic Ornstein-Uhlenbeck processes, jointly leveraging the evolutionary history of subclones and single-cell expression data. Applying our model to sublines derived from single cells of a mouse melanoma revealed that sublines with distinct phenotypes are underlined by different patterns of gene expression adaptation, indicating non-genetic mechanisms of cancer evolution. Sublines previously observed to be resistant to anti-CTLA4 treatment showed adaptive expression of genes related to invasion and non-canonical Wnt signaling, whereas sublines that responded to treatment showed adaptive expression of genes related to proliferation and canonical Wnt signaling. Our results suggest that clonal phenotypes emerge as the result of specific adaptivity patterns of gene expression. A record of this paper's transparent peer review process is included in the supplemental information.

Advancing cancer research through organoid technology

The complexity of tumors and the challenges associated with treatment often stem from the limitations of existing models in accurately replicating authentic tumors. Recently, organoid technology has emerged as an innovative platform for tumor research. This bioengineering approach enables researchers to simulate, in vitro, the interactions between tumors and their microenvironment, thereby enhancing the intricate interplay between tumor cells and their surroundings. Organoids also integrate multidimensional data, providing a novel paradigm for understanding tumor development and progression while facilitating precision therapy. Furthermore, advancements in imaging and genetic editing techniques have significantly augmented the potential of organoids in tumor research. This review explores the application of organoid technology for more precise tumor simulations and its specific contributions to cancer research advancements. Additionally, we discuss the challenges and evolving trends in developing comprehensive tumor models utilizing organoid technology.

In vivo dendritic cell reprogramming for cancer immunotherapy

Immunotherapy can lead to long-term survival for some cancer patients, yet generalized success has been hampered by insufficient antigen presentation and exclusion of immunogenic cells from the tumor microenvironment. Here, we developed an approach to reprogram tumor cells in vivo by adenoviral delivery of the transcription factors PU.1, IRF8, and BATF3, which enabled them to present antigens as type 1 conventional dendritic cells. Reprogrammed tumor cells remodeled their tumor microenvironment, recruited, and expanded polyclonal cytotoxic T cells; induced tumor regressions; and established long-term systemic immunity in multiple mouse melanoma models. In human tumor spheroids and xenografts, reprogramming to immunogenic dendritic-like cells progressed independently of immunosuppression, which usually limits immunotherapy. Our study paves the way for human clinical trials of in vivo immune cell reprogramming for cancer immunotherapy.

A machine learning model reveals expansive downregulation of ligand-receptor interactions that enhance lymphocyte infiltration in melanoma with developed resistance to immune checkpoint blockade

Immune checkpoint blockade (ICB) is a promising cancer therapy; however, resistance frequently develops. To explore ICB resistance mechanisms, we develop Immunotherapy Resistance cell-cell Interaction Scanner (IRIS), a machine learning model aimed at identifying cell-type-specific tumor microenvironment ligand-receptor interactions relevant to ICB resistance. Applying IRIS to deconvolved transcriptomics data of the five largest melanoma ICB cohorts, we identify specific downregulated interactions, termed resistance downregulated interactions (RDI), as tumors develop resistance. These RDIs often involve chemokine signaling and offer a stronger predictive signal for ICB response compared to upregulated interactions or the state-of-the-art published transcriptomics biomarkers. Validation across multiple independent melanoma patient cohorts and modalities confirms that RDI activity is associated with CD8 + T cell infiltration and highly manifested in hot/brisk tumors. This study presents a strongly predictive ICB response biomarker, highlighting the key role of downregulating chemotaxis-associated ligand-receptor interactions in inhibiting lymphocyte infiltration in resistant tumors.

IMMUNE AND MOLECULAR CORRELATES OF RESPONSE TO IMMUNOTHERAPY REVEALED BY BRAIN-METASTATIC MELANOMA MODELS

Despite the promising results of immune checkpoint blockade (ICB) therapy, outcomes for patients with brain metastasis (BrM) remain poor. Identifying resistance mechanisms has been hindered by limited access to patient samples and relevant preclinical models. Here, we developed two mouse melanoma BrM models that recapitulate the disparate responses to ICB seen in patients. We demonstrate that these models capture the cellular and molecular complexity of human disease and reveal key factors shaping the tumor microenvironment and influencing ICB response. BR1-responsive tumor cells express inflammatory programs that polarize microglia into reactive states, eliciting robust T cell recruitment. In contrast, BR3-resistant melanoma cells are enriched in neurological programs and exploit tolerance mechanisms to maintain microglia homeostasis and limit T cell infiltration. In humans, BR1 and BR3 expression signatures correlate positively or negatively with T cell infiltration and BrM patient outcomes, respectively. Our study provides clinically relevant models and uncovers mechanistic insights into BrM ICB responses, offering potential biomarkers and therapeutic targets to improve therapy efficacy.

Melanocyte lineage dynamics in development, growth and disease

Melanocytes evolved to produce the melanin that gives colour to our hair, eyes and skin. The melanocyte lineage also gives rise to melanoma, the most lethal form of skin cancer. The melanocyte lineage differentiates from neural crest cells during development, and most melanocytes reside in the skin and hair, where they are replenished by melanocyte stem cells. Because the molecular mechanisms necessary for melanocyte specification, migration, proliferation and differentiation are co-opted during melanoma initiation and progression, studying melanocyte development is directly relevant to human disease. Here, through the lens of advances in cellular omic and genomic technologies, we review the latest findings in melanocyte development and differentiation, and how these developmental pathways become dysregulated in disease.

Thin Amelanotic and Hypomelanotic Melanoma: Clinicopathological and Dermoscopic Features

Background and Objectives: Amelanotic/hypomelanotic melanomas (AHMs) account for 2-8% of all cutaneous melanomas. Due to their clinical appearance and the lack of specific dermoscopic indicators, AHMs are challenging to diagnose, particularly in thinner cutaneous lesions. The aim of our study was to evaluate the clinicopathological and dermoscopic features of thin AHMs. Identifying the baseline clinical-pathological features and dermoscopic aspects of thin AHMs is crucial to better understand this entity. Materials and Methods: We divided the AHM cohort into two groups based on Breslow thickness: thin (≤1.00 mm) and thick (>1.00 mm). This stratification helped identify any significant clinicopathological differences between the groups. For dermoscopic analysis, we employed the "pattern analysis" approach, which involves a simultaneous and subjective assessment of different criteria. Results: Out of the 2.800 melanomas analyzed for Breslow thickness, 153 were identified as AHMs. Among these, 65 patients presented with thin AHMs and 88 with thick AHMs. Red hair color and phototype II were more prevalent in patients with thin AHMs. The trunk was the most common anatomic site for thin AHMs. Patients with thin AHMs showed a higher number of multiple melanomas. Dermoscopic analysis revealed no significant difference between thin AHMs and thick AHMs, except for a more frequent occurrence of residual reticulum in thin AHMs. Conclusions: Thin AHMs typically affect individuals with lower phototypes and red hair color. These aspects can be related to the higher presence of pheomelanin, which provides limited protection against sun damage. This also correlates with the fact that the trunk, a site commonly exposed to intermittent sun exposure, is the primary anatomical location for thin AHMs. Multiple primary melanomas are more common in patients with thin AHMs, likely due to an intrinsic predisposition as well as greater periodic dermatologic follow-ups in this class of patients. Apart from the presence of residual reticulum, no other significant dermoscopic differences were observed, complicating the differential diagnosis between thin and thick AHMs based on dermoscopy alone.

Comparative Analysis of the Effect of the BRAF Inhibitor Dabrafenib in 2D and 3D Cell Culture Models of Human Metastatic Melanoma Cells

BACKGROUND/AIM

Melanoma, a variant of skin cancer, presents the highest mortality rates among all skin cancers. Despite advancements in targeted therapies, immunotherapies, and tissue culture techniques, the absence of an effective early treatment model remains a challenge. This study investigated the impact of dabrafenib on both 2D and 3D cell culture models with distinct molecular profiles.

MATERIALS AND METHODS

We developed a high-throughput workflow enabling drug screening on spheroids. Our approach involved cultivating 2D and 3D cultures derived from normal melanocytes and metastatic melanoma cells, treating them with dabrafenib and conducting viability, aggregation, migration, cell cycle, and apoptosis assays.

RESULTS

Dabrafenib exerted multifaceted influences, particularly on migration at concentrations of 10 and 25 μM. It induced a decrease in cell viability, impeded cellular adhesion to the matrix, inhibited cellular aggregation and spheroid formation, arrested the cell cycle in the G1 phase, and induced apoptosis.

CONCLUSION

These results confirm the therapeutic potential of dabrafenib in treating melanoma with the BRAF V600E mutation and that 3D models are validated models to study the potential of new molecules for therapeutic purposes. Furthermore, our study underscores the relevance of 3D models in simulating physiological in vivo microenvironments, providing insights into varied treatment responses between normal and tumor cells.

Malignant Melanoma: An Overview, New Perspectives, and Vitamin D Signaling

Melanoma, originating through malignant transformation of melanin-producing melanocytes, is a formidable malignancy, characterized by local invasiveness, recurrence, early metastasis, resistance to therapy, and a high mortality rate. This review discusses etiologic and risk factors for melanoma, diagnostic and prognostic tools, including recent advances in molecular biology, omics, and bioinformatics, and provides an overview of its therapy. Since the incidence of melanoma is rising and mortality remains unacceptably high, we discuss its inherent properties, including melanogenesis, that make this disease resilient to treatment and propose to use AI to solve the above complex and multidimensional problems. We provide an overview on vitamin D and its anticancerogenic properties, and report recent advances in this field that can provide solutions for the prevention and/or therapy of melanoma. Experimental papers and clinicopathological studies on the role of vitamin D status and signaling pathways initiated by its active metabolites in melanoma prognosis and therapy are reviewed. We conclude that vitamin D signaling, defined by specific nuclear receptors and selective activation by specific vitamin D hydroxyderivatives, can provide a benefit for new or existing therapeutic approaches. We propose to target vitamin D signaling with the use of computational biology and AI tools to provide a solution to the melanoma problem.

Transactivation of Met signaling by oncogenic Gnaq drives the evolution of melanoma in Hgf-Cdk4 mice

Recent pan-cancer genomic analyses have identified numerous oncogenic driver mutations that occur in a cell-type and tissue-specific distribution. For example, oncogenic mutations in Braf and Nras genes arise predominantly in melanocytic neoplasms of the epidermis, while oncogenic mutations in Gnaq/11 genes arise mostly in melanocytic lesions of the dermis or the uvea. The mechanisms promoting cell-type and tissue-specific oncogenic events currently remain poorly understood. Here, we report that Gnaq/11 hotspot mutations occur as early oncogenic drivers during the evolution of primary melanomas in Hgf-Cdk4 mice. Additional single base substitutions in the Trp53 gene and structural chromosomal aberrations favoring amplifications of the chromosomal region containing the Met receptor gene accumulate during serial tumor transplantation and in cell lines established in vitro. Mechanistically, we found that the GnaqQ209L mutation transactivates the Met receptor. Overexpression of oncogenic GnaqQ209L in the immortalized melanocyte cell line promoted in vivo growth that was enhanced by transgenic Hgf expression in the tumor microenvironment. This cross-signaling mechanism explains the selection of oncogenic Gnaq/11 in primary Hgf-Cdk4 melanomas and provides an example of how oncogenic driver mutations, intracellular signaling cascades, and microenvironmental cues cooperate to drive cancer development in a tissue-specific fashion.

Stochastic modelling of single-cell gene expression adaptation reveals non-genomic contribution to evolution of tumor subclones

Cancer progression is an evolutionary process driven by the selection of cells adapted to gain growth advantage. We present the first formal study on the adaptation of gene expression in subclonal evolution. We model evolutionary changes in gene expression as stochastic Ornstein-Uhlenbeck processes, jointly leveraging the evolutionary history of subclones and single-cell expression data. Applying our model to sublines derived from single cells of a mouse melanoma revealed that sublines with distinct phenotypes are underlined by different patterns of gene expression adaptation, indicating non-genetic mechanisms of cancer evolution. Interestingly, sublines previously observed to be resistant to anti-CTLA-4 treatment showed adaptive expression of genes related to invasion and non-canonical Wnt signaling, whereas sublines that responded to treatment showed adaptive expression of genes related to proliferation and canonical Wnt signaling. Our results suggest that clonal phenotypes emerge as the result of specific adaptivity patterns of gene expression.

Hypoxia is correlated with the tumor immune microenvironment: Potential application of immunotherapy in bladder cancer

OBJECTIVE

Hypoxia, which can considerably affect the tumor microenvironment, hinders the use of immunotherapy in bladder cancer (BLCA). Therefore, we aimed to identify reliable hypoxia-related biomarkers to guide clinical immunotherapy in BLCA.

METHODS

Using data downloaded from TCGA-BLCA cohort, we determined BLCA subtypes which divide 408 samples into different subtypes. Tumor immune infiltration levels of two clusters were quantified using ssGSEA, MCPcounter, EPIC, ESTIMATE, and TIMER algorithms. Next, we constructed a hypoxia score based on the expression of hypoxia-related genes. The IMvigor210 cohort and SubMap analysis were used to predict immunotherapeutic responses in patients with different hypoxia scores. Hub genes were screened using cytoscape, immunohistochemistry (IHC), and multispectral immunofluorescence were used to detect the spatial distribution of immune markers.

RESULTS

Patients with BLCA were categorized into cluster1 (n = 227) and Cluster2 (n = 181). Immune infiltration and expression of immune markers were higher in Cluster1. Immune infiltration was also more obvious in the high-hypoxia score group which related to a better predicted response to immunotherapy. IHC, and multispectral immunofluorescence confirmed the importance of TLR8 in immune infiltration and immune phenotype.

CONCLUSIONS

BLCA subtype can evaluate the infiltration of immune cells in the tumor microenvironment of different patients. Hypoxia score in this study could effectively predict immunotherapeutic responses in patients with BLCA. TLR8 may be a potential target for clinical immunotherapy.

High-throughput chemogenetic drug screening reveals PKC-RhoA/PKN as a targetable signaling vulnerability in GNAQ-driven uveal melanoma

Uveal melanoma (UM) is the most prevalent cancer of the eye in adults, driven by activating mutation of GNAQ/GNA11; however, there are limited therapies against UM and metastatic UM (mUM). Here, we perform a high-throughput chemogenetic drug screen in GNAQ-mutant UM contrasted with BRAF-mutant cutaneous melanoma, defining the druggable landscape of these distinct melanoma subtypes. Across all compounds, darovasertib demonstrates the highest preferential activity against UM. Our investigation reveals that darovasertib potently inhibits PKC as well as PKN/PRK, an AGC kinase family that is part of the "dark kinome." We find that downstream of the Gαq-RhoA signaling axis, PKN converges with ROCK to control FAK, a mediator of non-canonical Gαq-driven signaling. Strikingly, darovasertib synergizes with FAK inhibitors to halt UM growth and promote cytotoxic cell death in vitro and in preclinical metastatic mouse models, thus exposing a signaling vulnerability that can be exploited as a multimodal precision therapy against mUM.

Machine learning endometrial cancer risk prediction model: integrating guidelines of European Society for Medical Oncology with the tumor immune framework

OBJECTIVE

Current prognostic factors for endometrial cancer are not sufficient to predict recurrence in early stages. Treatment choices are based on the prognostic factors included in the risk classes defined by the ESMO-ESGO-ESTRO (European Society for Medical Oncology-European Society of Gynaecological Oncology-European Society for Radiotherapy and Oncology) consensus conference with the new biomolecular classification based on POLE, TP53, and microsatellite instability status. However, a minority of early stage cases relapse regardless of their low risk profiles. Integration of the immune context status to existing molecular based models has not been fully evaluated. This study aims to investigate whether the integration of the immune landscape in the tumor microenvironment could improve clinical risk prediction models and allow better profiling of early stages.

METHODS

Leveraging the potential of in silico deconvolution tools, we estimated the relative abundances of immune populations in public data and then applied feature selection methods to generate a machine learning based model for disease free survival probability prediction.

RESULTS

We included information on International Federation of Gynecology and Obstetrics (FIGO) stage, tumor mutational burden, microsatellite instability, POLEmut status, interferon γ signature, and relative abundances of monocytes, natural killer cells, and CD4+T cells to build a relapse prediction model and obtained a balanced accuracy of 69%. We further identified two novel early stage profiles that undergo different pathways of recurrence.

CONCLUSION

This study presents an extension of current prognostic factors for endometrial cancer by exploiting machine learning models and deconvolution techniques on available public biomolecular data. Prospective clinical trials are advisable to validate the early stage stratification.

Microsomal glutathione transferase 1 controls metastasis and therapeutic response in melanoma

While recent targeted and immunotherapies in malignant melanoma are encouraging, most patients acquire resistance, implicating a need to identify additional drug targets to improve outcomes. Recently, attention has been given to pathways that regulate redox homeostasis, especially the lipid peroxidase pathway that protects cells against ferroptosis. Here we identify microsomal glutathione S-transferase 1 (MGST1), a non-selenium-dependent glutathione peroxidase, as highly expressed in malignant and drug resistant melanomas and as a specific determinant of metastatic spread and therapeutic sensitivity. Loss of MGST1 in mouse and human melanoma enhanced cellular oxidative stress, and diminished glycolysis, oxidative phosphorylation, and pentose phosphate pathway. Gp100 activated pmel-1 T cells killed more Mgst1 KD than control melanoma cells and KD cells were more sensitive to cytotoxic anticancer drugs and ferroptotic cell death. When compared to control, mice bearing Mgst1 KD B16 tumors had more CD8+ T cell infiltration with reduced expression of inhibitory receptors and increased cytokine response, large reduction of lung metastases and enhanced survival. Targeting MGST1 alters the redox balance and limits metastases in melanoma, enhancing the therapeutic index for chemo- and immunotherapies.

Deactivation of ligand-receptor interactions enhancing lymphocyte infiltration drives melanoma resistance to Immune Checkpoint Blockade

Immune checkpoint blockade (ICB) is a promising cancer therapy; however, resistance often develops. To learn more about ICB resistance mechanisms, we developed IRIS (Immunotherapy Resistance cell-cell Interaction Scanner), a machine learning model aimed at identifying candidate ligand-receptor interactions (LRI) that are likely to mediate ICB resistance in the tumor microenvironment (TME). We developed and applied IRIS to identify resistance-mediating cell-type-specific ligand-receptor interactions by analyzing deconvolved transcriptomics data of the five largest melanoma ICB therapy cohorts. This analysis identifies a set of specific ligand-receptor pairs that are deactivated as tumors develop resistance, which we refer to as resistance deactivated interactions (RDI). Quite strikingly, the activity of these RDIs in pre-treatment samples offers a markedly stronger predictive signal for ICB therapy response compared to those that are activated as tumors develop resistance. Their predictive accuracy surpasses the state-of-the-art published transcriptomics biomarker signatures across an array of melanoma ICB datasets. Many of these RDIs are involved in chemokine signaling. Indeed, we further validate on an independent large melanoma patient cohort that their activity is associated with CD8+ T cell infiltration and enriched in hot/brisk tumors. Taken together, this study presents a new strongly predictive ICB response biomarker signature, showing that following ICB treatment resistant tumors turn inhibit lymphocyte infiltration by deactivating specific key ligand-receptor interactions.

Guadecitabine plus ipilimumab in unresectable melanoma: five-year follow-up and integrated multi-omic analysis in the phase 1b NIBIT-M4 trial

Association with hypomethylating agents is a promising strategy to improve the efficacy of immune checkpoint inhibitors-based therapy. The NIBIT-M4 was a phase Ib, dose-escalation trial in patients with advanced melanoma of the hypomethylating agent guadecitabine combined with the anti-CTLA-4 antibody ipilimumab that followed a traditional 3 + 3 design (NCT02608437). Patients received guadecitabine 30, 45 or 60 mg/m2/day subcutaneously on days 1 to 5 every 3 weeks starting on week 0 for a total of four cycles, and ipilimumab 3 mg/kg intravenously starting on day 1 of week 1 every 3 weeks for a total of four cycles. Primary outcomes of safety, tolerability, and maximum tolerated dose of treatment were previously reported. Here we report the 5-year clinical outcome for the secondary endpoints of overall survival, progression free survival, and duration of response, and an exploratory integrated multi-omics analysis on pre- and on-treatment tumor biopsies. With a minimum follow-up of 45 months, the 5-year overall survival rate was 28.9% and the median duration of response was 20.6 months. Re-expression of immuno-modulatory endogenous retroviruses and of other repetitive elements, and a mechanistic signature of guadecitabine are associated with response. Integration of a genetic immunoediting index with an adaptive immunity signature stratifies patients/lesions into four distinct subsets and discriminates 5-year overall survival and progression free survival. These results suggest that coupling genetic immunoediting with activation of adaptive immunity is a relevant requisite for achieving long term clinical benefit by epigenetic immunomodulation in advanced melanoma patients.

Tumors produce glucocorticoids by metabolite recycling, not synthesis, and activate Tregs to promote growth

Glucocorticoids are steroid hormones with potent immunosuppressive properties. Their primary source is the adrenals, where they are generated via de novo synthesis from cholesterol. In addition, many tissues have a recycling pathway in which glucocorticoids are regenerated from inactive metabolites by the enzyme 11β-hydroxysteroid dehydrogenase type 1 (11β-HSD1, encoded by Hsd11b1). Here, we find that multiple tumor types express Hsd11b1 and produce active glucocorticoids. Genetic ablation of Hsd11b1 in such cells had no effect on in vitro growth, but reduced in vivo tumor progression, which corresponded with increased frequencies of CD8+ tumor-infiltrating lymphocytes (TILs) expressing activation markers and producing effector cytokines. Tumor-derived glucocorticoids were found to promote signatures of Treg activation and suppress signatures of conventional T cell activation in tumor-infiltrating Tregs. Indeed, CD8+ T cell activation was restored and tumor growth reduced in mice with Treg-specific glucocorticoid receptor deficiency. Importantly, pharmacologic inhibition of 11β-HSD1 reduced tumor growth to the same degree as gene knockout and rendered immunotherapy-resistant tumors susceptible to PD-1 blockade. Given that HSD11B1 expression is upregulated in many human tumors and that inhibition of 11β-HSD1 is well tolerated in clinical studies, these data suggest that targeting 11β-HSD1 may be a beneficial adjunct in cancer therapy.

Activation of melanocortin-1 receptor signaling in melanoma cells impairs T cell infiltration to dampen antitumor immunity

Inhibition of T cell infiltration dampens antitumor immunity and causes resistance to immune checkpoint blockade (ICB) therapy. By in vivo CRISPR screening in B16F10 melanoma in female mice, here we report that loss of melanocortin-1 receptor (MC1R) in melanoma cells activates antitumor T cell response and overcomes resistance to ICB. Depletion of MC1R from another melanocytic melanoma model HCmel1274 also enhances ICB efficacy. By activating the GNAS-PKA axis, MC1R inhibits interferon-gamma induced CXCL9/10/11 transcription, thus impairing T cell infiltration into the tumor microenvironment. In human melanomas, high MC1R expression correlates with reduced CXCL9/10/11 expression, impaired T cell infiltration, and poor patient prognosis. Whereas MC1R activation is restricted to melanoma, GNAS activation by hotspot mutations is observed across diverse cancer types and is associated with reduced CXCL9/10/11 expression. Our study implicates MC1R as a melanoma immunotherapy target and suggests GNAS-PKA signaling as a pan-cancer oncogenic pathway inhibiting antitumor T cell response.

MacroH2A restricts inflammatory gene expression in melanoma cancer-associated fibroblasts by coordinating chromatin looping

MacroH2A has established tumour suppressive functions in melanoma and other cancers, but an unappreciated role in the tumour microenvironment. Using an autochthonous, immunocompetent mouse model of melanoma, we demonstrate that mice devoid of macroH2A variants exhibit increased tumour burden compared with wild-type counterparts. MacroH2A-deficient tumours accumulate immunosuppressive monocytes and are depleted of functional cytotoxic T cells, characteristics consistent with a compromised anti-tumour response. Single cell and spatial transcriptomics identify increased dedifferentiation along the neural crest lineage of the tumour compartment and increased frequency and activation of cancer-associated fibroblasts following macroH2A loss. Mechanistically, macroH2A-deficient cancer-associated fibroblasts display increased myeloid chemoattractant activity as a consequence of hyperinducible expression of inflammatory genes, which is enforced by increased chromatin looping of their promoters to enhancers that gain H3K27ac. In summary, we reveal a tumour suppressive role for macroH2A variants through the regulation of chromatin architecture in the tumour stroma with potential implications for human melanoma.

A Comprehensive Benchmark of Transcriptomic Biomarkers for Immune Checkpoint Blockades

Immune checkpoint blockades (ICBs) have revolutionized cancer therapy by inducing durable clinical responses, but only a small percentage of patients can benefit from ICB treatments. Many studies have established various biomarkers to predict ICB responses. However, different biomarkers were found with diverse performances in practice, and a timely and unbiased assessment has yet to be conducted due to the complexity of ICB-related studies and trials. In this study, we manually curated 29 published datasets with matched transcriptome and clinical data from more than 1400 patients, and uniformly preprocessed these datasets for further analyses. In addition, we collected 39 sets of transcriptomic biomarkers, and based on the nature of the corresponding computational methods, we categorized them into the gene-set-like group (with the self-contained design and the competitive design, respectively) and the deconvolution-like group. Next, we investigated the correlations and patterns of these biomarkers and utilized a standardized workflow to systematically evaluate their performance in predicting ICB responses and survival statuses across different datasets, cancer types, antibodies, biopsy times, and combinatory treatments. In our benchmark, most biomarkers showed poor performance in terms of stability and robustness across different datasets. Two scores (TIDE and CYT) had a competitive performance for ICB response prediction, and two others (PASS-ON and EIGS_ssGSEA) showed the best association with clinical outcome. Finally, we developed ICB-Portal to host the datasets, biomarkers, and benchmark results and to implement the computational methods for researchers to test their custom biomarkers. Our work provided valuable resources and a one-stop solution to facilitate ICB-related research.

Autoimmune alleles at the major histocompatibility locus modify melanoma susceptibility

Autoimmunity and cancer represent two different aspects of immune dysfunction. Autoimmunity is characterized by breakdowns in immune self-tolerance, while impaired immune surveillance can allow for tumorigenesis. The class I major histocompatibility complex (MHC-I), which displays derivatives of the cellular peptidome for immune surveillance by CD8+ T cells, serves as a common genetic link between these conditions. As melanoma-specific CD8+ T cells have been shown to target melanocyte-specific peptide antigens more often than melanoma-specific antigens, we investigated whether vitiligo- and psoriasis-predisposing MHC-I alleles conferred a melanoma-protective effect. In individuals with cutaneous melanoma from both The Cancer Genome Atlas (n = 451) and an independent validation set (n = 586), MHC-I autoimmune-allele carrier status was significantly associated with a later age of melanoma diagnosis. Furthermore, MHC-I autoimmune-allele carriers were significantly associated with decreased risk of developing melanoma in the Million Veteran Program (OR = 0.962, p = 0.024). Existing melanoma polygenic risk scores (PRSs) did not predict autoimmune-allele carrier status, suggesting these alleles provide orthogonal risk-relevant information. Mechanisms of autoimmune protection were neither associated with improved melanoma-driver mutation association nor improved gene-level conserved antigen presentation relative to common alleles. However, autoimmune alleles showed higher affinity relative to common alleles for particular windows of melanocyte-conserved antigens and loss of heterozygosity of autoimmune alleles caused the greatest reduction in presentation for several conserved antigens across individuals with loss of HLA alleles. Overall, this study presents evidence that MHC-I autoimmune-risk alleles modulate melanoma risk unaccounted for by current PRSs.

RNF8 enhances the sensitivity of PD-L1 inhibitor against melanoma through ubiquitination of galectin-3 in stroma

The failure of melanoma immunotherapy can be mediated by immunosuppression in the tumor microenvironment (TME), and insufficient activation of effector T cells against the tumor. Here, we show that inhibition of galectin-3 (gal-3) enhances the infiltration of T cells in TME and improves the sensitivity of anti-PD-L1 therapy. We identify that RNF8 downregulated the expression of gal-3 by K48-polyubiquitination and promoted gal-3 degradation via the ubiquitin-proteasome system. RNF8 deficiency in the host but sufficiency in implanted melanoma results in immune exclusion and tumor progression due to the upregulation of gal-3. Upregulation of gal-3 decreased the immune cell infiltration by restricting IL-12 and IFN-γ. Inhibition of gal-3 reverses immunosuppression and induces immune cell infiltration in the tumor microenvironment. Moreover, gal-3 inhibitor treatment can increase the sensitivity of PD-L1 inhibitors via increasing immune cell infiltration and enhancing immune response in tumors. This study reveals a previously unrecognized immunoregulation function of RNF8 and provides a promising strategy for the therapy of "cold" tumors. Tremendous effects of melanoma treatment can be achieved by facilitating immune cell infiltration combined with anti-PD-L1 treatment.

Late-Stage Metastatic Melanoma Emerges through a Diversity of Evolutionary Pathways

Understanding the evolutionary pathways to metastasis and resistance to immune-checkpoint inhibitors (ICI) in melanoma is critical for improving outcomes. Here, we present the most comprehensive intrapatient metastatic melanoma dataset assembled to date as part of the Posthumous Evaluation of Advanced Cancer Environment (PEACE) research autopsy program, including 222 exome sequencing, 493 panel-sequenced, 161 RNA sequencing, and 22 single-cell whole-genome sequencing samples from 14 ICI-treated patients. We observed frequent whole-genome doubling and widespread loss of heterozygosity, often involving antigen-presentation machinery. We found KIT extrachromosomal DNA may have contributed to the lack of response to KIT inhibitors of a KIT-driven melanoma. At the lesion-level, MYC amplifications were enriched in ICI nonresponders. Single-cell sequencing revealed polyclonal seeding of metastases originating from clones with different ploidy in one patient. Finally, we observed that brain metastases that diverged early in molecular evolution emerge late in disease. Overall, our study illustrates the diverse evolutionary landscape of advanced melanoma.

SIGNIFICANCE

Despite treatment advances, melanoma remains a deadly disease at stage IV. Through research autopsy and dense sampling of metastases combined with extensive multiomic profiling, our study elucidates the many mechanisms that melanomas use to evade treatment and the immune system, whether through mutations, widespread copy-number alterations, or extrachromosomal DNA. See related commentary by Shain, p. 1294. This article is highlighted in the In This Issue feature, p. 1275.

Identification of CircRNA signature associated with tumor immune infiltration to predict therapeutic efficacy of immunotherapy

Circular RNAs (circRNAs) play important roles in the regulation of cancer. However, the clinical implications and regulatory networks of circRNAs in cancer patients receiving immune checkpoint blockades (ICB) have not been fully elucidated. Here, we characterize circRNA expression profiles in two independent cohorts of 157 ICB-treated advanced melanoma patients and reveal overall overexpression of circRNAs in ICB non-responders in both pre-treatment and early during therapy. Then, we construct circRNA-miRNA-mRNA regulatory networks to reveal circRNA-related signaling pathways in the context of ICB treatment. Further, we construct an ICB-related circRNA signature (ICBcircSig) score model based on progression-free survival-related circRNAs to predict immunotherapy efficacy. Mechanistically, the overexpression of ICBcircSig circTMTC3 and circFAM117B could increase PD-L1 expression via the miR-142-5p/PD-L1 axis, thus reducing T cell activity and leading to immune escape. Overall, our study characterizes circRNA profiles and regulatory networks in ICB-treated patients, and highlights the clinical utility of circRNAs as predictive biomarkers of immunotherapy.

A Therapeutic Vaccine in Combination with Cyclic GMP-AMP Cures More Differentiated Melanomas in Mice

We have identified a combinational immunotherapy termed TheraVac vaccine (TheraVac) that can cure multiple large established mouse tumors, but it failed to cure melanoma in mice. TheraVac consists of an immunostimulating arm containing an agonist (HMGN1 [N1]) for TLR4 and an agonist (R848) for TLR7/8 that synergize to activate tumor-infiltrating dendritic cells (DCs) and promote Th1 immune responses. The second arm uses an immune checkpoint blockade, anti-PDL-1, to diminish tumor-associated immunosuppression. In this study, we investigated supplementation of TheraVac by a stimulator of IFN genes (STING) agonist, cyclic GMP-AMP (cGAMP), because together they synergize in activating DCs and produced more immunostimulating IL-12p70 and TNF-α cytokines. The synergistic activation and maturation of DCs is dependent on the activation of tank binding kinase-1 (TBK1). Treatment of three different melanin-producing mouse melanomas (B16F1, M3, and M4) with intratumoral delivery of cGAMP and TheraVac eradicated 60-80% of these melanomas. Immunoprofiling of M3 tumor treated with TheraVac plus cGAMP showed an increase in CD8+ CTLs and macrophages in the tumor. There was also a marked increase of CD4, CD8 effector and memory T cells and generation of functional tumor-specific CTLs in tumor-draining lymph nodes. The resultant tumor-free mice were selectively resistant to subsequent challenge with the same tumors, indicating long-term tumor-specific protective immunity. Overall, our findings have important implications for clinical trials with a combination of these immunotherapeutics to cure melanin-producing human melanomas, without the need for exogenous tumor Ags and no clear toxic effects in mice.

Preclinical models for drug discovery for metastatic disease

Despite many advances, metastatic disease remains essentially uncurable. Thus, there is an urgent need to better understand mechanisms that promote metastasis, drive tumor evolution, and underlie innate and acquired drug resistance. Sophisticated preclinical models that recapitulate the complex tumor ecosystem are key to this process. We begin with syngeneic and patient-derived mouse models that are the backbone of most preclinical studies. Second, we present some unique advantages of fish and fly models. Third, we consider the strengths of 3D culture models for resolving remaining knowledge gaps. Finally, we provide vignettes on multiplexed technologies to advance our understanding of metastatic disease.

Melanoma clonal subline analysis uncovers heterogeneity-driven immunotherapy resistance mechanisms

Intratumoral heterogeneity (ITH) can promote cancer progression and treatment failure, but the complexity of the regulatory programs and contextual factors involved complicates its study. To understand the specific contribution of ITH to immune checkpoint blockade (ICB) response, we generated single cell-derived clonal sublines from an ICB-sensitive and genetically and phenotypically heterogeneous mouse melanoma model, M4. Genomic and single cell transcriptomic analyses uncovered the diversity of the sublines and evidenced their plasticity. Moreover, a wide range of tumor growth kinetics were observed in vivo , in part associated with mutational profiles and dependent on T cell-response. Further inquiry into melanoma differentiation states and tumor microenvironment (TME) subtypes of untreated tumors from the clonal sublines demonstrated correlations between highly inflamed and differentiated phenotypes with the response to anti-CTLA-4 treatment. Our results demonstrate that M4 sublines generate intratumoral heterogeneity at both levels of intrinsic differentiation status and extrinsic TME profiles, thereby impacting tumor evolution during therapeutic treatment. These clonal sublines proved to be a valuable resource to study the complex determinants of response to ICB, and specifically the role of melanoma plasticity in immune evasion mechanisms.

The cancer epigenome: Non-cell autonomous player in tumor immunity

Dysregulation of the tumor-intrinsic epigenetic circuit is a key driver event for the development of cancer. Accumulating evidence suggests that epigenetic and/or genetic drivers stimulate intrinsic oncogenic pathways as well as extrinsic factors that modulate the immune system. These modulations indeed shape the tumor microenvironment (TME), allowing pro-oncogenic factors to become oncogenic, thereby contributing to cancer development and progression. Here we review the epigenetic dysregulation arising in cancer cells that disseminates throughout the TME and beyond. Recent CRISPR screening has elucidated key epigenetic drivers that play important roles in the proliferation of cancer cells (intrinsic) and inhibition of antitumor immunity (extrinsic), which lead to the development and progression of cancer. These epigenetic players can serve as promising targets for cancer therapy as a dual (two-in-one)-targeted approach. Considering the interplay between cancer and the immune system as a key determinant of immunotherapy, we discuss a novel lineage-tracing technology that enables longitudinal monitoring of cancer and immune phenotypic heterogeneity and fate paths during cancer development, progression, and therapeutic interventions.

Discovery of Targets for Immune-Metabolic Antitumor Drugs Identifies Estrogen-Related Receptor Alpha

Drugs that kill tumors through multiple mechanisms have the potential for broad clinical benefits. Here, we first developed an in silico multiomics approach (BipotentR) to find cancer cell-specific regulators that simultaneously modulate tumor immunity and another oncogenic pathway and then used it to identify 38 candidate immune-metabolic regulators. We show the tumor activities of these regulators stratify patients with melanoma by their response to anti-PD-1 using machine learning and deep neural approaches, which improve the predictive power of current biomarkers. The topmost identified regulator, ESRRA, is activated in immunotherapy-resistant tumors. Its inhibition killed tumors by suppressing energy metabolism and activating two immune mechanisms: (i) cytokine induction, causing proinflammatory macrophage polarization, and (ii) antigen-presentation stimulation, recruiting CD8+ T cells into tumors. We also demonstrate a wide utility of BipotentR by applying it to angiogenesis and growth suppressor evasion pathways. BipotentR (http://bipotentr.dfci.harvard.edu/) provides a resource for evaluating patient response and discovering drug targets that act simultaneously through multiple mechanisms.

SIGNIFICANCE

BipotentR presents resources for evaluating patient response and identifying targets for drugs that can kill tumors through multiple mechanisms concurrently. Inhibition of the topmost candidate target killed tumors by suppressing energy metabolism and effects on two immune mechanisms. This article is highlighted in the In This Issue feature, p. 517.

Research landscape and trends of melanoma immunotherapy: A bibliometric analysis

Background

Immunotherapy for lung cancer has been a hot research area for years. This bibliometric analysis was intended to present research trends on melanoma immunotherapy.

Method

On April 1, 2022, the authors identified 2,109 papers on melanoma immunotherapy using the Web of Science and extracted their general information and the total number of citations. The authors then conducted a bibliometric analysis to present the research landscape, clarify the research trends, and determine the most cited papers (top-papers) as well as major journals on melanoma immunotherapy. Subsequently, recent research hotspots were identified by analyzing the latest articles in major journals.

Results

The total and median number of citations of these 2,109 papers on melanoma immunotherapy was 137,686 and 11, respectively. "Improved survival with ipilimumab in patients with metastatic melanoma" by Hodi et al. was the most cited paper (9,824 citations). Among the journals, the top-paper number (16), average citations per paper (2,510.7), and top-papers rate (100%) of New England Journal of Medicine were the highest. Corresponding authors represented the USA took part in most articles (784). Since 2016, the hottest research area has changed from CTLA-4 to PD-1.

Conclusions

This bibliometric analysis comprehensively and quantitatively presents the research trends and hotspots based on 2,109 relevant publications, and further suggests future research directions. The researchers can benefit in selecting journals and in finding potential collaborators. This study can help researchers gain a comprehensive impression of the research landscape, historical development, and current hotspots in melanoma immunotherapy and can provide inspiration for future research.

From simplicity to complexity in current melanoma models

Despite the recent impressive clinical success of immunotherapy against melanoma, development of primary and adaptive resistance against immune checkpoint inhibitors remains a major issue in a large number of treated patients. This highlights the need for melanoma models that replicate the tumor's intricate dynamics in the tumor microenvironment (TME) and associated immune suppression to study possible resistance mechanisms in order to improve current and test novel therapeutics. While two-dimensional melanoma cell cultures have been widely used to perform functional genomics screens in a high-throughput fashion, they are not suitable to answer more complex scientific questions. Melanoma models have also been established in a variety of experimental (humanized) animals. However, due to differences in physiology, such models do not fully represent human melanoma development. Therefore, fully human three-dimensional in vitro models mimicking melanoma cell interactions with the TME are being developed to address this need for more physiologically relevant models. Such models include melanoma organoids, spheroids, and reconstructed human melanoma-in-skin cultures. Still, while major advances have been made to complement and replace animals, these in vitro systems have yet to fully recapitulate human tumor complexity. Lastly, technical advancements have been made in the organ-on-chip field to replicate functions and microstructures of in vivo human tissues and organs. This review summarizes advancements made in understanding and treating melanoma and specifically aims to discuss the progress made towards developing melanoma models, their applications, limitations, and the advances still needed to further facilitate the development of therapeutics.

Resistance to immune checkpoint blockade: Mechanisms, counter-acting approaches, and future directions

Immunotherapies seek to unleash the immune system against cancer cells. While a variety of immunotherapies exist, one of the most commonly used is immune checkpoint blockade, which refers to the use of antibodies to interfere with immunosuppressive signaling through immune checkpoint molecules. Therapies against various checkpoints have had success in the clinic across cancer types. However, the efficacy of checkpoint inhibitors has varied across different cancer types and non-responsive patient populations have emerged. Non-responders to these therapies have highlighted the importance of understanding underlying mechanisms of resistance in order to predict which patients will respond and to tailor individual treatment paradigms. In this review we discuss the literature surrounding tumor mediated mechanisms of immune checkpoint resistance. We also describe efforts to overcome resistance and combine checkpoint inhibitors with additional immunotherapies. Finally, we provide insight into the future of immune checkpoint blockade, including the need for improved preclinical modeling and predictive biomarkers to facilitate personalized cancer treatments for patients.

Melanoma central nervous system metastases: An update to approaches, challenges, and opportunities

Brain metastases are the most common brain malignancy. This review discusses the studies presented at the third annual meeting of the Melanoma Research Foundation in the context of other recent reports on the biology and treatment of melanoma brain metastases (MBM). Although symptomatic MBM patients were historically excluded from immunotherapy trials, efforts from clinicians and patient advocates have resulted in more inclusive and even dedicated clinical trials for MBM patients. The results of checkpoint inhibitor trials were discussed in conversation with current standards of care for MBM patients, including steroids, radiotherapy, and targeted therapy. Advances in the basic scientific understanding of MBM, including the role of astrocytes and metabolic adaptations to the brain microenvironment, are exposing new vulnerabilities which could be exploited for therapeutic purposes. Technical advances including single-cell omics and multiplex imaging are expanding our understanding of the MBM ecosystem and its response to therapy. This unprecedented level of spatial and temporal resolution is expected to dramatically advance the field in the coming years and render novel treatment approaches that might improve MBM patient outcomes.

SRC-RAC1 signaling drives drug resistance to BRAF inhibition in de-differentiated cutaneous melanomas

Rare gain-of-function mutations in RAC1 drive drug resistance to targeted BRAF inhibition in cutaneous melanoma. Here, we show that wildtype RAC1 is a critical driver of growth and drug resistance, but only in a subset of melanomas with elevated markers of de-differentiation. Similarly, SRC inhibition also selectively sensitized de-differentiated melanomas to BRAF inhibition. One possible mechanism may be the suppression of the de-differentiated state, as SRC and RAC1 maintained markers of de-differentiation in human melanoma cells. The functional differences between melanoma subtypes suggest that the clinical management of cutaneous melanoma can be enhanced by the knowledge of differentiation status. To simplify the task of classification, we developed a binary classification strategy based on a small set of ten genes. Using this gene set, we reliably determined the differentiation status previously defined by hundreds of genes. Overall, our study informs strategies that enhance the precision of BRAFi by discovering unique vulnerabilities of the de-differentiated cutaneous melanoma subtype and creating a practical method to resolve differentiation status.

A novel subtype to predict prognosis and treatment response with DNA driver methylation-transcription in ovarian cancer

Aims: To identify a novel subtype with DNA driver methylation-transcriptomic multiomics and predict prognosis and therapy response in serous ovarian cancer (SOC). Methods: SOC cohorts with both mRNA and methylation were collected, and DNA driver methylation (DNAme) was identified with the MithSig method. A novel prognostic subtype was developed by integrating the information on DNAme and prognosis-regulated DNAme-associated mRNA by similarity network fusion. Results: 43 overlapped DNAme were identified in three independent cohorts. SOC patients were categorized into three distinct subtypes by integrated multiomics. There were differences in prognosis, tumor microenvironment and response to therapy among the subtypes. Conclusion: This study identified 43 DNAmes and proposes a novel subtype toward personalized chemotherapy and immunotherapy for SOC patients based on multiomics.

Mutated processes predict immune checkpoint inhibitor therapy benefit in metastatic melanoma

Immune Checkpoint Inhibitor (ICI) therapy has revolutionized treatment for advanced melanoma; however, only a subset of patients benefit from this treatment. Despite considerable efforts, the Tumor Mutation Burden (TMB) is the only FDA-approved biomarker in melanoma. However, the mechanisms underlying TMB association with prolonged ICI survival are not entirely understood and may depend on numerous confounding factors. To identify more interpretable ICI response biomarkers based on tumor mutations, we train classifiers using mutations within distinct biological processes. We evaluate a variety of feature selection and classification methods and identify key mutated biological processes that provide improved predictive capability compared to the TMB. The top mutated processes we identify are leukocyte and T-cell proliferation regulation, which demonstrate stable predictive performance across different data cohorts of melanoma patients treated with ICI. This study provides biologically interpretable genomic predictors of ICI response with substantially improved predictive performance over the TMB.

Mechanisms of tumor resistance to immune checkpoint blockade and combination strategies to overcome resistance

Immune checkpoint blockade (ICB) has rapidly transformed the treatment paradigm for various cancer types. Multiple single or combinations of ICB treatments have been approved by the US Food and Drug Administration, providing more options for patients with advanced cancer. However, most patients could not benefit from these immunotherapies due to primary and acquired drug resistance. Thus, a better understanding of the mechanisms of ICB resistance is urgently needed to improve clinical outcomes. Here, we focused on the changes in the biological functions of CD8⁺ T cells to elucidate the underlying resistance mechanisms of ICB therapies and summarized the advanced coping strategies to increase ICB efficacy. Combinational ICB approaches and individualized immunotherapies require further in-depth investigation to facilitate longer-lasting efficacy and a more excellent safety of ICB in a broader range of patients.

Alpha radionuclide-chelated radioimmunotherapy promoters enable local radiotherapy/chemodynamic therapy to discourage cancer progression

BACKGROUND

Astatine-211 is an α-emitter with high-energy α-ray and high cytotoxicity for cancer cells. However, the targeted alpha therapy (TAT) also suffers from insufficient systematic immune activation, resulting in tumor metastasis and relapse. Combined immune checkpoint blockade (ICB) with chemodynamic therapy (CDT) could boost antitumor immunity, which may magnify the immune responses of TAT. This study aims to discourage tumor metastasis and relapse by tri-model TAT-CDT-ICB strategy.

METHODS

We successfully designed Mn-based radioimmunotherapy promoters (²¹¹At-ATE-MnO₂-BSA), which are consisting of ²¹¹At, MnO₂ and bovine serum albumin (BSA). The efficacy of ²¹¹At-ATE-MnO₂-BSA was studied as monotherapy or in combination with anti-PD-L1 in both metastatic and relapse models. The immune effects of radioimmunotherapy promoters on cytotoxic T lymphocytes and dendritic cells (DCs) were analyzed by flow cytometry. Enzyme-linked immunosorbent assay and immunofluorescence were used to explore the underlying mechanism.

RESULTS

Such radioimmunotherapy promoters could not only enhance the therapeutic outcomes of TAT and CDT, but also induce robust anti-cancer immune activity by activating dendritic cells. More intriguingly, ²¹¹At-ATE-MnO₂-BSA could effectively suppress the growths of primary tumors and distant tumors when combined with immune checkpoint inhibitors.

CONCLUSIONS

The tri-model TAT-CDT-ICB strategy provides a long-term immunological memory, which can protect against tumor rechallenge after eliminating original tumors. Therefore, this work presents a novel approach for TAT-CDT-ICB tri-modal cancer therapy with repressed metastasis and relapse in clinics.

The oxidoreductase CLIC4 is required to maintain mitochondrial function and resistance to exogenous oxidants in breast cancer cells

The chloride intracellular channel-4 (CLIC4) is one of the six highly conserved proteins in the CLIC family that share high structural homology with GST-omega in the GST superfamily. While CLIC4 is a multifunctional protein that resides in multiple cellular compartments, the discovery of its enzymatic glutaredoxin-like activity in vitro suggested that it could function as an antioxidant. Here, we found that deleting CLIC4 from murine 6DT1 breast tumor cells using CRISPR enhanced the accumulation of reactive oxygen species (ROS) and sensitized cells to apoptosis in response to H2O2 as a ROS-inducing agent. In intact cells, H2O2 increased the expression of both CLIC4 mRNA and protein. In addition, increased superoxide production in 6DT1 cells lacking CLIC4 was associated with mitochondrial hyperactivity including increased mitochondrial membrane potential and mitochondrial organelle enlargement. In the absence of CLIC4, however, H2O2-induced apoptosis was associated with low expression and degradation of the antiapoptotic mitochondrial protein Bcl2 and the negative regulator of mitochondrial ROS, UCP2. Furthermore, transcriptomic profiling of H2O2-treated control and CLIC4-null cells revealed upregulation of genes associated with ROS-induced apoptosis and downregulation of genes that sustain mitochondrial functions. Accordingly, tumors that formed from transplantation of CLIC4-deficient 6DT1 cells were highly necrotic. These results highlight a critical role for CLIC4 in maintaining redox-homeostasis and mitochondrial functions in 6DT1 cells. Our findings also raise the possibility of targeting CLIC4 to increase cancer cell sensitivity to chemotherapeutic drugs that are based on elevating ROS in cancer cells.

MCP-1/CCR2 axis inhibition sensitizes the brain microenvironment against melanoma brain metastasis progression

Development of resistance to chemo- and immunotherapies often occurs following treatment of melanoma brain metastasis (MBM). The brain microenvironment (BME), particularly astrocytes, cooperate toward MBM progression by upregulating secreted factors, among which we found that monocyte chemoattractant protein-1 (MCP-1) and its receptors, CCR2 and CCR4, were overexpressed in MBM compared with primary lesions. Among other sources of MCP-1 in the brain, we show that melanoma cells altered astrocyte secretome and evoked MCP-1 expression and secretion, which in turn induced CCR2 expression in melanoma cells, enhancing in vitro tumorigenic properties, such as proliferation, migration, and invasion of melanoma cells. In vivo pharmacological blockade of MCP-1 or molecular knockout of CCR2/CCR4 increased the infiltration of cytotoxic CD8⁺ T cells and attenuated the immunosuppressive phenotype of the BME as shown by decreased infiltration of Tregs and tumor-associated macrophages/microglia in several models of intracranially injected MBM. These in vivo strategies led to decreased MBM outgrowth and prolonged the overall survival of the mice. Our findings highlight the therapeutic potential of inhibiting interactions between BME and melanoma cells for the treatment of this disease.

Cross-species identification of cancer resistance-associated genes that may mediate human cancer risk

Cancer is a predominant disease across animals. We applied a comparative genomics approach to systematically characterize genes whose conservation levels correlate positively (PC) or negatively (NC) with cancer resistance estimates across 193 vertebrates. Pathway analysis reveals that NC genes are enriched for metabolic functions and PC genes in cell cycle regulation, DNA repair, and immune response, pointing to their corresponding roles in mediating cancer risk. We find that PC genes are less tolerant to loss-of-function (LoF) mutations, are enriched in cancer driver genes, and are associated with germline mutations that increase human cancer risk. Their relevance to cancer risk is further supported via the analysis of mouse functional genomics and cancer mortality of zoo mammals' data. In sum, our study describes a cross-species genomic analysis pointing to candidate genes that may mediate human cancer risk.

Multi-modal molecular programs regulate melanoma cell state

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.

Multimodal predictors for precision immunotherapy

Immune checkpoint blockade (ICB) unleashes immune cells to attack tumors, thereby inducing durable clinical responses in many cancer types. The number of patients responding to ICB is modest, however, and combination treatments are likely needed to overcome the multifaceted suppressive pathways active in the tumor microenvironment (TME). The development of precision immuno-oncology (IO) strategies allowing to identify the optimal treatment of each patient upfront is therefore a pivotal question in the field of cancer immunotherapy. Although single-parameter biomarkers can enrich for response to ICB, their predictive capacity is far from perfect and their clinical utility is complicated by their continuous nature and the difficulty to determine cut-offs that reliably distinguish responding patients from those without clinical benefit. The antitumor immune response that is induced or reinvigorated by immunotherapy is a complex cascade of events requiring the interplay of multiple cell types. To move towards precision IO, it is therefore essential to understand for each individual patient at which level(s) the antitumor immune response failed and how it can be therapeutically restored. Holistic approaches to profile human tumor microenvironments and treatment-induced responses may help to identify critical rate-limiting factors of antitumor immunity. These factors need to be translated into clinically applicable multimodal predictors that allow for the selection of the best IO treatment. This review discusses strategies to (i) create such holistic views of antitumor immunity, (ii) identify measurable parameters capturing the complexity of a patient's immune status, and (iii) facilitate the incorporation of precision IO research in the clinic.

Targeting WEE1/AKT restores p53-dependent NK cell activation to induce immune checkpoint blockade responses in 'cold' melanoma

Immunotherapy has revolutionized cancer treatment. Unfortunately, most tumor types do not respond to immunotherapy due to a lack of immune infiltration or 'cold' tumor microenvironment (TME), a contributing factor in treatment failure. Activation of the p53 pathway can increase apoptosis of cancer cells, leading to enhanced antigen presentation, and can stimulate natural killer (NK) cells through expression of stress ligands. Therefore, modulation of the p53 pathway in cancer cells with wildtype TP53 has the potential to enhance tumor immunogenicity to NK cells, produce an inflammatory TME, and ultimately lead to tumor regression. In this study, we report simultaneous targeting of the AKT/WEE1 pathways is a novel and tolerable approach to synergistically induce p53 activation to inhibit tumor development. This approach reduced the growth of melanoma cells and induced plasma membrane surface localization of the ER-resident protein calreticulin, an indicator of immunogenic cell death (ICD). Increase in ICD led to enhanced expression of stress ligands recognized by the activating NK cell receptor NKG2D, promoting tumor lysis. WEE1/AKT inhibition resulted in recruitment and activation of immune cells, including NK cells, in the TME, triggering an inflammatory cascade that transformed the 'cold' TME of B16F10 melanoma into a 'hot' TME that responded to anti-PD-1, resulting in complete regression of established tumors. These results suggest that AKT/WEE1 pathway inhibition is a potential approach to broaden the utility of class-leading anti-PD-1 therapies by enhancing p53-mediated, NK cell-dependent tumor inflammation and supports the translation of this novel approach to further improve response rates for metastatic melanoma.

Deconvolving Clinically Relevant Cellular Immune Cross-talk from Bulk Gene Expression Using CODEFACS and LIRICS Stratifies Patients with Melanoma to Anti-PD-1 Therapy

The tumor microenvironment (TME) is a complex mixture of cell types whose interactions affect tumor growth and clinical outcome. To discover such interactions, we developed CODEFACS (COnfident DEconvolution For All Cell Subsets), a tool deconvolving cell type-specific gene expression in each sample from bulk expression, and LIRICS (Ligand-Receptor Interactions between Cell Subsets), a statistical framework prioritizing clinically relevant ligand-receptor interactions between cell types from the deconvolved data. We first demonstrate the superiority of CODEFACS versus the state-of-the-art deconvolution method CIBERSORTx. Second, analyzing The Cancer Genome Atlas, we uncover cell type-specific ligand-receptor interactions uniquely associated with mismatch-repair deficiency across different cancer types, providing additional insights into their enhanced sensitivity to anti-programmed cell death protein 1 (PD-1) therapy compared with other tumors with high neoantigen burden. Finally, we identify a subset of cell type-specific ligand-receptor interactions in the melanoma TME that stratify survival of patients receiving anti-PD-1 therapy better than some recently published bulk transcriptomics-based methods.

SIGNIFICANCE

This work presents two new computational methods that can deconvolve a large collection of bulk tumor gene expression profiles into their respective cell type-specific gene expression profiles and identify cell type-specific ligand-receptor interactions predictive of response to immune-checkpoint blockade therapy. This article is highlighted in the In This Issue feature, p. 873.

Cancer Cell Phenotype Plasticity as a Driver of Immune Escape in Melanoma

Immunotherapies blocking negative immune checkpoints are now approved for the treatment of a growing number of cancers. However, even in metastatic melanoma, where sustained responses are observed, a significant number of patients still do not respond or display resistance. Increasing evidence indicates that non-genetic cancer cell-intrinsic alterations play a key role in resistance to therapies and immune evasion. Cancer cell plasticity, mainly associated with the epithelial-to-mesenchymal transition in carcinoma, relies on transcriptional, epigenetic or translational reprogramming. In melanoma, an EMT-like dedifferentiation process is characterized by the acquisition of invasive or neural crest stem cell-like features. Herein, we discuss recent findings on the specific roles of phenotypic reprogramming of melanoma cells in driving immune evasion and resistance to immunotherapies. The mechanisms by which dedifferentiated melanoma cells escape T cell lysis, mediate T cell exclusion or remodel the immune microenvironment will be detailed. The expanded knowledge on tumor cell plasticity in melanoma should contribute to the development of novel therapeutic combination strategies to further improve outcomes in this deadly metastatic cancer.

Targeting SOX10-deficient cells to reduce the dormant-invasive phenotype state in melanoma

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.

A Series of BRAF- and NRAS-Driven Murine Melanoma Cell Lines with Inducible Gene Modulation Capabilities

Murine cancer cell lines are powerful research tools to complement studies in genetically engineered mouse models. We have established 21 melanoma cell lines from embryonic stem cell-genetically engineered mouse models driven by alleles that model the most frequent genetic alterations in human melanoma. In addition, these cell lines harbor regulatory alleles for the genomic integration of transgenes and the regulation of expression of such transgenes. In this study, we report a comprehensive characterization of these cell lines. Specifically, we validated melanocytic origin, driver allele recombination and expression, and activation of the oncogenic MAPK and protein kinase B pathways. We further tested tumor formation in syngeneic immunocompetent recipients as well as the functionality of the integrated Tet-ON system and recombination-mediated cassette exchange homing cassette. Finally, by deleting the transcription factor MAFG with an inducible CRISPR/Cas9 approach, we show the utility of the regulatory alleles for candidate gene modulation. These cell lines will be a valuable resource for studying melanoma biology and therapy.

Establishing a Prognostic Model Based on Ulceration and Immune Related Genes in Melanoma Patients and Identification of EIF3B as a Therapeutic Target

Ulceration and immune status are independent prognostic factors for survival in melanoma patients. Herein univariate Cox regression analysis revealed 53 ulcer-immunity-related DEGs. We performed consensus clustering to divide The Cancer Genome Atlas (TCGA) cohort (n = 467) into three subtypes with different prognosis and biological functions, followed by validation in three merged Gene Expression Omnibus (GEO) cohorts (n = 399). Multiomics approach was used to assess differences among the subtypes. Cluster 3 showed relatively lesser amplification and expression of immune checkpoint genes. Moreover, Cluster 3 lacked immune-related pathways and immune cell infiltration, and had higher proportion of non-responders to immunotherapy. We also constructed a prognostic model based on ulceration and immune related genes in melanoma. EIF3B was a hub gene in the intersection between genes specific to Cluster 3 and those pivotal for melanoma growth (DepMap, https://depmap.org/portal/download/). High EIF3B expression in TCGA and GEO datasets was related to worst prognosis. In vitro models revealed that EIF3B knockdown inhibited melanoma cell migration and invasion, and decreased TGF-β1 level in supernatant compared with si-NC cells. EIF3B expression was negatively correlated with immune-related signaling pathways, immune cell gene signatures, and immune checkpoint gene expression. Moreover, its low expression could predict partial response to anti-PD-1 immunotherapy. To summarize, we established a prognostic model for melanoma and identified the role of EIF3B in melanoma progression and immunotherapy resistance development.