MotivationVenn and Euler diagrams are a popular yet inadequate solution for quantitative visualization of set intersections. A scalable alternative to Venn and Euler diagrams for visualizing intersecting sets and their properties is needed.ResultsWe developed UpSetR, an open source R package that employs a scalable matrix-based visualization to show intersections of sets, their size, and other properties.Availability and implementationUpSetR is available at https://github.com/hms-dbmi/UpSetR/ and released under the MIT License. A Shiny app is available at https://gehlenborglab.shinyapps.io/upsetr/.Supplementary information Supplementary data are available at Bioinformatics online.

Immune-checkpoint blockade (ICB) has demonstrated efficacy in many tumor types, but predictors of responsiveness to anti-PD1 ICB are incompletely characterized. In this study, we analyzed a clinically annotated cohort of patients with melanoma (n = 144) treated with anti-PD1 ICB, with whole-exome and whole-transcriptome sequencing of pre-treatment tumors. We found that tumor mutational burden as a predictor of response was confounded by melanoma subtype, whereas multiple novel genomic and transcriptomic features predicted selective response, including features associated with MHC-I and MHC-II antigen presentation. Furthermore, previous anti-CTLA4 ICB exposure was associated with different predictors of response compared to tumors that were naive to ICB, suggesting selective immune effects of previous exposure to anti-CTLA4 ICB. Finally, we developed parsimonious models integrating clinical, genomic and transcriptomic features to predict intrinsic resistance to anti-PD1 ICB in individual tumors, with validation in smaller independent cohorts limited by the availability of comprehensive data. Broadly, we present a framework to discover predictive features and build models of ICB therapeutic response.

Chromatin regulators play a broad role in regulating gene expression and, when gone awry, can lead to cancer. Here, we demonstrate that ablation of the histone demethylase LSD1 in cancer cells increases repetitive element expression, including endogenous retroviral elements (ERVs), and decreases expression of RNA-induced silencing complex (RISC) components. Significantly, this leads to double-stranded RNA (dsRNA) stress and activation of type 1 interferon, which stimulates anti-tumor T cell immunity and restrains tumor growth. Furthermore, LSD1 depletion enhances tumor immunogenicity and T cell infiltration in poorly immunogenic tumors and elicits significant responses of checkpoint blockade-refractory mouse melanoma to anti-PD-1 therapy. Consistently, TCGA data analysis shows an inverse correlation between LSD1 expression and CD8 T cell infiltration in various human cancers. Our study identifies LSD1 as a potent inhibitor of anti-tumor immunity and responsiveness to immunotherapy and suggests LSD1 inhibition combined with PD-(L)1 blockade as a novel cancer treatment strategy.

Venn and Euler diagrams are a popular yet inadequate solution for quantitative visualization of set intersections. A scalable alternative to Venn and Euler diagrams for visualizing intersecting sets and their properties is needed. We developed UpSetR, an open source R package that employs a scalable matrix-based visualization to show intersections of sets, their size, and other properties. UpSetR is available at https://cran.r-project.org/package=UpSetR and released under the MIT License. A Shiny app is available at https://gehlenborglab.shinyapps.io/upsetr/.

Immune checkpoint blockade (ICB) therapies, which potentiate the body’s natural immune response against tumor cells, have shown immense promise in the treatment of various cancers. Currently, tumor mutational burden (TMB) and programmed death ligand 1 (PD-L1) expression are the primary biomarkers evaluated for clinical management of cancer patients across histologies. However, the wide range of responses has demonstrated that the specific molecular and genetic characteristics of each patient’s tumor and immune system must be considered to maximize treatment efficacy. Here, we review the various biological pathways and emerging biomarkers implicated in response to PD-(L)1 and cytotoxic T lymphocyte-associated antigen 4 (CTLA-4) therapies, including oncogenic signaling pathways, human leukocyte antigen (HLA) variability, mutation and neoantigen burden, microbiome composition, endogenous retroviruses (ERV), and deficiencies in chromatin remodeling and DNA damage repair (DDR) machinery. We also discuss several mechanisms that have been observed to confer resistance to ICB, such as loss of phosphatase and tensin homolog (PTEN), loss of major histocompatibility complex (MHC) I/II expression, and activation of the indoleamine 2,3-dioxygenase 1 (IDO1) and transforming growth factor beta (TGFβ) pathways. Clinical trials testing the combination of PD-(L)1 or CTLA-4 blockade with molecular mediators of these pathways are becoming more common and may hold promise for improving treatment efficacy and response. Ultimately, some of the genes and molecular mechanisms highlighted in this review may serve as novel biological targets or therapeutic vulnerabilities to improve clinical outcomes in patients.

Immune checkpoint blockade (ICB) therapy revolutionized cancer treatment, but many patients with impaired MHC-I expression remain refractory. Here, we combined FACS-based genome-wide CRISPR screens with a data-mining approach to identify drugs that can upregulate MHC-I without inducing PD-L1. CRISPR screening identified TRAF3, a suppressor of the NFκB pathway, as a negative regulator of MHC-I but not PD-L1. The Traf3-knockout gene expression signature is associated with better survival in ICB-naïve patients with cancer and better ICB response. We then screened for drugs with similar transcriptional effects as this signature and identified Second Mitochondria-derived Activator of Caspase (SMAC) mimetics. We experimentally validated that the SMAC mimetic birinapant upregulates MHC-I, sensitizes cancer cells to T cell–dependent killing, and adds to ICB efficacy. Our findings provide preclinical rationale for treating tumors expressing low MHC-I expression with SMAC mimetics to enhance sensitivity to immunotherapy. The approach used in this study can be generalized to identify other drugs that enhance immunotherapy efficacy. Significance: MHC-I loss or downregulation in cancer cells is a major mechanism of resistance to T cell–based immunotherapies. Our study reveals that birinapant may be used for patients with low baseline MHC-I to enhance ICB response. This represents promising immunotherapy opportunities given the biosafety profile of birinapant from multiple clinical trials. This article is highlighted in the In This Issue feature, p. 1307

We performed harmonized molecular and clinical analysis on 1,048 melanomas and discovered markedly different global genomic properties among subtypes ( BRAF, (N)RAS, NF1 , Triple Wild-Type), subtype-specific preferences for secondary driver genes, and active mutational processes previously unreported in melanoma. Secondary driver genes significantly enriched in specific subtypes reflected preferential dysregulation of additional pathways, such as induction of TGF-β signaling in BRAF melanomas and inactivation of the SWI/SNF complex in (N)RAS melanomas, and select co-mutation patterns coordinated selective response to immune checkpoint blockade. We also defined the mutational landscape of Triple Wild-Type melanomas and identified enrichment of DNA repair defect signatures in this subtype, which were associated with transcriptional downregulation of key DNA repair genes and may revive previously discarded or currently unconsidered therapeutic modalities for genomically stratified melanoma patient subsets. Broadly, harmonized meta-analysis of melanoma whole-exomes identified distinct molecular drivers that may point to multiple opportunities for biological and therapeutic investigation.

PURPOSE Next-generation sequencing (NGS) of tumor and germline DNA is foundational for precision oncology, with rapidly expanding diagnostic, prognostic, and therapeutic implications. Although few question the importance of NGS in modern oncology care, the process of gathering primary molecular data, integrating it into electronic health records, and optimally using it as part of a clinical workflow remains far from seamless. Numerous challenges persist around data standards and interoperability, and clinicians frequently face difficulties in managing the growing amount of genomic knowledge required to care for patients and keep up to date. METHODS This review provides a descriptive analysis of genomic data workflows for NGS data in clinical oncology and issues that arise from the inconsistent use of standards for sharing data across systems. Potential solutions are described. RESULTS NGS technology, especially for somatic genomics, is well established and widely used in routine patient care, quality measurement, and research. Available genomic knowledge bases play an evolving role in patient management but lack harmonization with one another. Questions about their provenance and timeliness of updating remain. Potentially useful standards for sharing genomic data, such as HL7 FHIR and mCODE, remain primarily in the research and/or development stage. Nonetheless, their impact will likely be seen as uptake increases across care settings and laboratories. The specific use case of ASCO CancerLinQ, as a clinicogenomic database, is discussed. CONCLUSION Because the electronic health records of today seem ill suited for managing genomic data, other solutions are required, including universal data standards and applications that use application programming interfaces, along with a commitment on the part of sequencing laboratories to consistently provide structured genomic data for clinical use.