LORIS robustly predicts patient outcomes with immune checkpoint blockade therapy using common clinical, pathologic and genomic features.

Despite the revolutionary impact of immune checkpoint blockade (ICB) in cancer treatment, accurately predicting patient responses remains challenging. Here, we analyzed a large dataset of 2,881 ICB-treated and 841 non-ICB-treated patients across 18 solid tumor types, encompassing a wide range of clinical, pathologic and genomic features. We developed a clinical score called LORIS (logistic regression-based immunotherapy-response score) using a six-feature logistic regression model. LORIS outperforms previous signatures in predicting ICB response and identifying responsive patients even with low tumor mutational burden or programmed cell death 1 ligand 1 expression. LORIS consistently predicts patient objective response and short-term and long-term survival across most cancer types. Moreover, LORIS showcases a near-monotonic relationship with ICB response probability and patient survival, enabling precise patient stratification. As an accurate, interpretable method using a few readily measurable features, LORIS may help improve clinical decision-making in precision medicine to maximize patient benefit. LORIS is available as an online tool at https://loris.ccr.cancer.gov/ .

Decoupling the correlation between cytotoxic and exhausted T lymphocyte states enhances melanoma immunotherapy response prediction.

Cytotoxic T lymphocyte (CTL) and terminal exhausted T lymphocyte (ETL) activities crucially influence immune checkpoint inhibitor (ICI) response. Despite this, the efficacy of ETL and CTL transcriptomic signatures for response prediction remains limited. Investigating this across the TCGA and publicly available single-cell cohorts, we find a strong positive correlation between ETL and CTL expression signatures in most cancers. We hence posited that their limited predictability arises due to their mutually canceling effects on ICI response. Thus, we developed DETACH, a computational method to identify a gene set whose expression pinpoints to a subset of melanoma patients where the CTL and ETL correlation is low. DETACH enhances CTL's prediction accuracy, outperforming existing signatures. DETACH signature genes activity also demonstrates a positive correlation with lymphocyte infiltration and the prevalence of reactive T cells in the tumor microenvironment (TME), advancing our understanding of the CTL cell state within the TME.

ASS1 metabolically contributes to the nuclear and cytosolic p53-mediated DNA damage response.

Downregulation of the urea cycle enzyme argininosuccinate synthase (ASS1) in multiple tumors is associated with a poor prognosis partly because of the metabolic diversion of cytosolic aspartate for pyrimidine synthesis, supporting proliferation and mutagenesis owing to nucleotide imbalance. Here, we find that prolonged loss of ASS1 promotes DNA damage in colon cancer cells and fibroblasts from subjects with citrullinemia type I. Following acute induction of DNA damage with doxorubicin, ASS1 expression is elevated in the cytosol and the nucleus with at least a partial dependency on p53; ASS1 metabolically restrains cell cycle progression in the cytosol by restricting nucleotide synthesis. In the nucleus, ASS1 and ASL generate fumarate for the succination of SMARCC1, destabilizing the chromatin-remodeling complex SMARCC1-SNF5 to decrease gene transcription, specifically in a subset of the p53-regulated cell cycle genes. Thus, following DNA damage, ASS1 is part of the p53 network that pauses cell cycle progression, enabling genome maintenance and survival. Loss of ASS1 contributes to DNA damage and promotes cell cycle progression, likely contributing to cancer mutagenesis and, hence, adaptability potential.

Prediction of DNA methylation-based tumor types from histopathology in central nervous system tumors with deep learning.

Precision in the diagnosis of diverse central nervous system (CNS) tumor types is crucial for optimal treatment. DNA methylation profiles, which capture the methylation status of thousands of individual CpG sites, are state-of-the-art data-driven means to enhance diagnostic accuracy but are also time consuming and not widely available. Here, to address these limitations, we developed Deep lEarning from histoPathoLOgy and methYlation (DEPLOY), a deep learning model that classifies CNS tumors to ten major categories from histopathology. DEPLOY integrates three distinct components: the first classifies CNS tumors directly from slide images ('direct model'), the second initially generates predictions for DNA methylation beta values, which are subsequently used for tumor classification ('indirect model'), and the third classifies tumor types directly from routinely available patient demographics. First, we find that DEPLOY accurately predicts beta values from histopathology images. Second, using a ten-class model trained on an internal dataset of 1,796 patients, we predict the tumor categories in three independent external test datasets including 2,156 patients, achieving an overall accuracy of 95% and balanced accuracy of 91% on samples that are predicted with high confidence. These results showcase the potential future use of DEPLOY to assist pathologists in diagnosing CNS tumors within a clinically relevant short time frame.

PERCEPTION predicts patient response and resistance to treatment using single-cell transcriptomics of their tumors.

Tailoring optimal treatment for individual cancer patients remains a significant challenge. To address this issue, we developed PERCEPTION (PERsonalized Single-Cell Expression-Based Planning for Treatments In ONcology), a precision oncology computational pipeline. Our approach uses publicly available matched bulk and single-cell (sc) expression profiles from large-scale cell-line drug screens. These profiles help build treatment response models based on patients' sc-tumor transcriptomics. PERCEPTION demonstrates success in predicting responses to targeted therapies in cultured and patient-tumor-derived primary cells, as well as in two clinical trials for multiple myeloma and breast cancer. It also captures the resistance development in patients with lung cancer treated with tyrosine kinase inhibitors. PERCEPTION outperforms published state-of-the-art sc-based and bulk-based predictors in all clinical cohorts. PERCEPTION is accessible at https://github.com/ruppinlab/PERCEPTION . Our work, showcasing patient stratification using sc-expression profiles of their tumors, will encourage the adoption of sc-omics profiling in clinical settings, enhancing precision oncology tools based on sc-omics.

Outcome differences by sex in oncology clinical trials.

Identifying sex differences in outcomes and toxicity between males and females in oncology clinical trials is important and has also been mandated by National Institutes of Health policies. Here we analyze the Trialtrove database, finding that, strikingly, only 472/89,221 oncology clinical trials (0.5%) had curated post-treatment sex comparisons. Among 288 trials with comparisons of survival, outcome, or response, 16% report males having statistically significant better survival outcome or response, while 42% reported significantly better survival outcome or response for females. The strongest differences are in trials of EGFR inhibitors in lung cancer and rituximab in non-Hodgkin's lymphoma (both favoring females). Among 44 trials with side effect comparisons, more trials report significantly lesser side effects in males (N = 22) than in females (N = 13). Thus, while statistical comparisons between sexes in oncology trials are rarely reported, important differences in outcome and toxicity exist. These considerable outcome and toxicity differences highlight the need for reporting sex differences more thoroughly going forward.

Programming a Ferroptosis-to-Apoptosis Transition Landscape Revealed Ferroptosis Biomarkers and Repressors for Cancer Therapy.

Ferroptosis and apoptosis are key cell-death pathways implicated in several human diseases including cancer. Ferroptosis is driven by iron-dependent lipid peroxidation and currently has no characteristic biomarkers or gene signatures. Here a continuous phenotypic gradient between ferroptosis and apoptosis coupled to transcriptomic and metabolomic landscapes is established. The gradual ferroptosis-to-apoptosis transcriptomic landscape is used to generate a unique, unbiased transcriptomic predictor, the Gradient Gene Set (GGS), which classified ferroptosis and apoptosis with high accuracy. Further GGS optimization using multiple ferroptotic and apoptotic datasets revealed highly specific ferroptosis biomarkers, which are robustly validated in vitro and in vivo. A subset of the GGS is associated with poor prognosis in breast cancer patients and PDXs and contains different ferroptosis repressors. Depletion of one representative, PDGFA-assaociated protein 1(PDAP1), is found to suppress basal-like breast tumor growth in a mouse model. Omics and mechanistic studies revealed that ferroptosis is associated with enhanced lysosomal function, glutaminolysis, and the tricarboxylic acid (TCA) cycle, while its transition into apoptosis is attributed to enhanced endoplasmic reticulum(ER)-stress and phosphatidylethanolamine (PE)-to-phosphatidylcholine (PC) metabolic shift. Collectively, this study highlights molecular mechanisms underlying ferroptosis execution, identified a highly predictive ferroptosis gene signature with prognostic value, ferroptosis versus apoptosis biomarkers, and ferroptosis repressors for breast cancer therapy.

Transcriptomic Profiling of Plasma Extracellular Vesicles Enables Reliable Annotation of the Cancer-Specific Transcriptome and Molecular Subtype.

Longitudinal monitoring of patients with advanced cancers is crucial to evaluate both disease burden and treatment response. Current liquid biopsy approaches mostly rely on the detection of DNA-based biomarkers. However, plasma RNA analysis can unleash tremendous opportunities for tumor state interrogation and molecular subtyping. Through the application of deep learning algorithms to the deconvolved transcriptomes of RNA within plasma extracellular vesicles (evRNA), we successfully predicted consensus molecular subtypes in patients with metastatic colorectal cancer. Analysis of plasma evRNA also enabled monitoring of changes in transcriptomic subtype under treatment selection pressure and identification of molecular pathways associated with recurrence. This approach also revealed expressed gene fusions and neoepitopes from evRNA. These results demonstrate the feasibility of using transcriptomic-based liquid biopsy platforms for precision oncology approaches, spanning from the longitudinal monitoring of tumor subtype changes to the identification of expressed fusions and neoantigens as cancer-specific therapeutic targets, sans the need for tissue-based sampling. SIGNIFICANCE: The development of an approach to interrogate molecular subtypes, cancer-associated pathways, and differentially expressed genes through RNA sequencing of plasma extracellular vesicles lays the foundation for liquid biopsy-based longitudinal monitoring of patient tumor transcriptomes.

DNA and RNA base editors can correct the majority of pathogenic single nucleotide variants.

The majority of human genetic diseases are caused by single nucleotide variants (SNVs) in the genome sequence. Excitingly, new genomic techniques known as base editing have opened efficient pathways to correct erroneous nucleotides. Due to reliance on deaminases, which have the capability to convert A to I(G) and C to U, the direct applicability of base editing might seem constrained in terms of the range of mutations that can be reverted. In this evaluation, we assess the potential of DNA and RNA base editing methods for treating human genetic diseases. Our findings indicate that 62% of pathogenic SNVs found within genes can be amended by base editing; 30% are G>A and T>C SNVs that can be corrected by DNA base editing, and most of them by RNA base editing as well, and 29% are C>T and A>G SNVs that can be corrected by DNA base editing directed to the complementary strand. For each, we also present several factors that affect applicability such as bystander and off-target occurrences. For cases where editing the mismatched nucleotide is not feasible, we introduce an approach that calculates the optimal substitution of the deleterious amino acid with a new amino acid, further expanding the scope of applicability. As personalized therapy is rapidly advancing, our demonstration that most SNVs can be treated by base editing is of high importance. The data provided will serve as a comprehensive resource for those seeking to design therapeutic base editors and study their potential in curing genetic diseases.

Chromosome 7 to the rescue: overcoming chromosome 10 loss in gliomas.

The co-occurrence of chromosome 10 loss and chromosome 7 gain in gliomas is the most frequent loss-gain co-aneuploidy pair in human cancers, a phenomenon that has been investigated without resolution since the late 1980s. Expanding beyond previous gene-centric studies, we investigate the co-occurrence in a genome-wide manner taking an evolutionary perspective. First, by mining large tumor aneuploidy data, we predict that the more likely order is 10 loss followed by 7 gain. Second, by analyzing extensive genomic and transcriptomic data from both patients and cell lines, we find that this co-occurrence can be explained by functional rescue interactions that are highly enriched on 7, which can possibly compensate for any detrimental consequences arising from the loss of 10. Finally, by analyzing transcriptomic data from normal, non-cancerous, human brain tissues, we provide a plausible reason why this co-occurrence happens preferentially in cancers originating in certain regions of the brain.

A systematic analysis of the landscape of synthetic lethality-driven precision oncology.

BACKGROUND: Synthetic lethality (SL) denotes a genetic interaction between two genes whose co-inactivation is detrimental to cells. Because more than 25 years have passed since SL was proposed as a promising way to selectively target cancer vulnerabilities, it is timely to comprehensively assess its impact so far and discuss its future. METHODS: We systematically analyzed the literature and clinical trial data from the PubMed and Trialtrove databases to portray the preclinical and clinical landscape of SL oncology. FINDINGS: We identified 235 preclinically validated SL pairs and found 1,207 pertinent clinical trials, and the number keeps increasing over time. About one-third of these SL clinical trials go beyond the typically studied DNA damage response (DDR) pathway, testifying to the recently broadening scope of SL applications in clinical oncology. We find that SL oncology trials have a greater success rate than non-SL-based trials. However, about 75% of the preclinically validated SL interactions have not yet been tested in clinical trials. CONCLUSIONS: Dissecting the recent efforts harnessing SL to identify predictive biomarkers, novel therapeutic targets, and effective combination therapy, our systematic analysis reinforces the hope that SL may serve as a key driver of precision oncology going forward. FUNDING: Funded by the Samsung Research Funding & Incubation Center of Samsung Electronics, the Institute of Information & Communications Technology Planning & Evaluation (IITP) grant funded by the Republic of Korea government (MSIT), the Kwanjeong Educational Foundation, the Intramural Research Program of the National Institutes of Health (NIH), National Cancer Institute (NCI), and Center for Cancer Research (CCR).

Transcriptional control of leukemogenesis by the chromatin reader SGF29.

ABSTRACT: Aberrant expression of stem cell-associated genes is a common feature in acute myeloid leukemia (AML) and is linked to leukemic self-renewal and therapy resistance. Using AF10-rearranged leukemia as a prototypical example of the recurrently activated "stemness" network in AML, we screened for chromatin regulators that sustain its expression. We deployed a CRISPR-Cas9 screen with a bespoke domain-focused library and identified several novel chromatin-modifying complexes as regulators of the TALE domain transcription factor MEIS1, a key leukemia stem cell (LSC)-associated gene. CRISPR droplet sequencing revealed that many of these MEIS1 regulators coordinately controlled the transcription of several AML oncogenes. In particular, we identified a novel role for the Tudor-domain-containing chromatin reader protein SGF29 in the transcription of AML oncogenes. Furthermore, SGF29 deletion impaired leukemogenesis in models representative of multiple AML subtypes in multiple AML subtype models. Our studies reveal a novel role for SGF29 as a nononcogenic dependency in AML and identify the SGF29 Tudor domain as an attractive target for drug discovery.

Charting the transcriptomic landscape of primary and metastatic cancers in relation to their origin and target normal tissues.

Metastasis is a leading cause of cancer-related deaths, yet understanding how metastatic tumors adapt from their origin to target tissues is challenging. To address this, we assessed whether primary and metastatic tumors resemble their tissue of origin or target tissue in terms of gene expression. We analyzed gene expression profiles from various cancer types, including single-cell and bulk RNA-seq data, in both paired and unpaired primary and metastatic patient cohorts. We quantified the transcriptomic distances between tumor samples and their normal tissues, revealing that primary tumors are more similar to their tissue of origin, while metastases shift towards the target tissue. Pathway-level analysis highlighted critical transcriptomic changes during metastasis. Notably, primary cancers exhibited higher activity in cancer hallmarks, including Activating Invasion and Metastasis , compared to metastatic cancers. This comprehensive landscape analysis provides insight into how cancer tumors adapt to their metastatic environments, providing a transcriptome-wide view of the processes involved.

Temporal genomic analysis of melanoma rejection identifies regulators of tumor immune evasion.

Decreased intra-tumor heterogeneity (ITH) correlates with increased patient survival and immunotherapy response. However, even highly homogenous tumors may display variability in their aggressiveness, and how immunologic-factors impinge on their aggressiveness remains understudied. Here we studied the mechanisms responsible for the immune-escape of murine tumors with low ITH. We compared the temporal growth of homogeneous, genetically-similar single-cell clones that are rejected vs. those that are not-rejected after transplantation in-vivo using single-cell RNA sequencing and immunophenotyping. Non-rejected clones showed high infiltration of tumor-associated-macrophages (TAMs), lower T-cell infiltration, and increased T-cell exhaustion compared to rejected clones. Comparative analysis of rejection-associated gene expression programs, combined with in-vivo CRISPR knockout screens of candidate mediators, identified Mif (macrophage migration inhibitory factor) as a regulator of immune rejection. Mif knockout led to smaller tumors and reversed non-rejection-associated immune composition, particularly, leading to the reduction of immunosuppressive macrophage infiltration. Finally, we validated these results in melanoma patient data.

Elevated A-to-I RNA editing in COVID-19 infected individuals.

Given the current status of coronavirus disease 2019 (COVID-19) as a global pandemic, it is of high priority to gain a deeper understanding of the disease's development and how the virus impacts its host. Adenosine (A)-to-Inosine (I) RNA editing is a post-transcriptional modification, catalyzed by the ADAR family of enzymes, that can be considered part of the inherent cellular defense mechanism as it affects the innate immune response in a complex manner. It was previously reported that various viruses could interact with the host's ADAR enzymes, resulting in epigenetic changes both to the virus and the host. Here, we analyze RNA-seq of nasopharyngeal swab specimens as well as whole-blood samples of COVID-19 infected individuals and show a significant elevation in the global RNA editing activity in COVID-19 compared to healthy controls. We also detect specific coding sites that exhibit higher editing activity. We further show that the increment in editing activity during the disease is temporary and returns to baseline shortly after the symptomatic period. These significant epigenetic changes may contribute to the immune system response and affect adverse outcomes seen in post-viral cases.

Prediction of cancer treatment response from histopathology images through imputed transcriptomics.

Advances in artificial intelligence have paved the way for leveraging hematoxylin and eosin (H&E)-stained tumor slides for precision oncology. We present ENLIGHT-DeepPT, an approach for predicting response to multiple targeted and immunotherapies from H&E-slides. In difference from existing approaches that aim to predict treatment response directly from the slides, ENLIGHT-DeepPT is an indirect two-step approach consisting of (1) DeepPT, a new deep-learning framework that predicts genome-wide tumor mRNA expression from slides, and (2) ENLIGHT, which predicts response based on the DeepPT inferred expression values. DeepPT successfully predicts transcriptomics in all 16 TCGA cohorts tested and generalizes well to two independent datasets. Our key contribution is showing that ENLIGHT-DeepPT successfully predicts true responders in five independent patients' cohorts involving four different treatments spanning six cancer types with an overall odds ratio of 2.44, increasing the baseline response rate by 43.47% among predicted responders, without the need for any treatment data for training. Furthermore, its prediction accuracy on these datasets is comparable to a supervised approach predicting the response directly from the images, which needs to be trained and tested on the same cohort. ENLIGHT-DeepPT future application could provide clinicians with rapid treatment recommendations to an array of different therapies and importantly, may contribute to advancing precision oncology in developing countries.

Increased interleukin-6/C-reactive protein levels are associated with the upregulation of the adenosine pathway and serve as potential markers of therapeutic resistance to immune checkpoint inhibitor-based therapies in non-small cell lung cancer.

BACKGROUND: Systemic immune activation, hallmarked by C-reactive protein (CRP) and interleukin-6 (IL-6), can modulate antitumor immune responses. In this study, we evaluated the role of IL-6 and CRP in the stratification of patients with non-small cell lung cancer (NSCLC) treated with immune checkpoint inhibitors (ICIs). We also interrogated the underlying immunosuppressive mechanisms driven by the IL-6/CRP axis. METHODS: In cohort A (n=308), we estimated the association of baseline CRP with objective response rate (ORR), progression-free survival (PFS), and overall survival (OS) in patients with NSCLC treated with ICIs alone or with chemo-immunotherapy (Chemo-ICI). Baseline tumor bulk RNA sequencing (RNA-seq) of lung adenocarcinomas (LUADs) treated with pembrolizumab (cohort B, n=59) was used to evaluate differential expression of purine metabolism, as well as correlate IL-6 expression with PFS. CODEFACS approach was applied to deconvolve cohort B to characterize the tumor microenvironment by reconstructing the cell-type-specific transcriptome from bulk expression. Using the LUAD cohort from The Cancer Genome Atlas (TCGA) we explored the correlation between IL-6 expression and adenosine gene signatures. In a third cohort (cohort C, n=18), plasma concentrations of CRP, adenosine 2a receptor (A2aR), and IL-6 were measured using ELISA. RESULTS: In cohort A, 67.2% of patients had a baseline CRP≥10 mg/L (CRP-H). Patients with CRP-H achieved shorter OS (8.6 vs 14.8 months; p=0.006), shorter PFS (3.3 vs 6.6 months; p=0.013), and lower ORR (24.7% vs 46.3%; p=0.015). After adjusting for relevant clinical variables, CRP-H was confirmed as an independent predictor of increased risk of death (HR 1.51, 95% CI: 1.09 to 2.11) and lower probability of achieving disease response (OR 0.34, 95% CI: 0.13 to 0.89). In cohort B, RNA-seq analysis demonstrated higher IL-6 expression on tumor cells of non-responders, along with a shorter PFS (p<0.05) and enrichment of the purinergic pathway. Within the TCGA LUAD cohort, tumor IL-6 expression strongly correlated with the adenosine signature (R=0.65; p<2.2e-16). Plasma analysis in cohort C demonstrated that CRP-H patients had a greater median baseline level of A2aR (6.0 ng/mL vs 1.3 ng/mL; p=0.01). CONCLUSIONS: This study demonstrates CRP as a readily available blood-based prognostic biomarker in ICI-treated NSCLC. Additionally, we elucidate a potential link of the CRP/IL-6 axis with the immunosuppressive adenosine signature pathway that could drive inferior outcomes to ICIs in NSCLC and also offer novel therapeutic avenues.

Pan-Cancer Analysis of Patient Tumor Single-Cell Transcriptomes Identifies Promising Selective and Safe Chimeric Antigen Receptor Targets in Head and Neck Cancer.

Chimeric antigen receptor (CAR) T cell therapies have yielded transformative clinical successes for patients with blood tumors, but their full potential remains to be unleashed against solid tumors. One challenge is finding selective targets, which we define intuitively to be cell surface proteins that are expressed widely by cancer cells but minimally by healthy cells in the tumor microenvironment and other normal tissues. Analyzing patient tumor single-cell transcriptomics data, we first defined and quantified selectivity and safety scores of existing CAR targets for indications in which they are in clinical trials or approved. We then sought new candidate cell surface CAR targets that have better selectivity and safety scores than those currently being tested. Remarkably, in almost all cancer types, we could not find such better targets, testifying to the near optimality of the current target space. However, in human papillomavirus (HPV)-negative head and neck squamous cell carcinoma (HNSC), for which there is currently a dearth of existing CAR targets, we identified a total of twenty candidate novel CAR targets, five of which have both superior selectivity and safety scores. These newly identified cell surface targets lay a basis for future investigations that may lead to better CAR treatments in HNSC.

Expanding PROTACtable genome universe of E3 ligases.

Proteolysis-targeting chimera (PROTAC) and other targeted protein degradation (TPD) molecules that induce degradation by the ubiquitin-proteasome system (UPS) offer new opportunities to engage targets that remain challenging to be inhibited by conventional small molecules. One fundamental element in the degradation process is the E3 ligase. However, less than 2% amongst hundreds of E3 ligases in the human genome have been engaged in current studies in the TPD field, calling for the recruiting of additional ones to further enhance the therapeutic potential of TPD. To accelerate the development of PROTACs utilizing under-explored E3 ligases, we systematically characterize E3 ligases from seven different aspects, including chemical ligandability, expression patterns, protein-protein interactions (PPI), structure availability, functional essentiality, cellular location, and PPI interface by analyzing 30 large-scale data sets. Our analysis uncovers several E3 ligases as promising extant PROTACs. In total, combining confidence score, ligandability, expression pattern, and PPI, we identified 76 E3 ligases as PROTAC-interacting candidates. We develop a user-friendly and flexible web portal ( https://hanlaboratory.com/E3Atlas/ ) aimed at assisting researchers to rapidly identify E3 ligases with promising TPD activities against specifically desired targets, facilitating the development of these therapies in cancer and beyond.

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.