The Elongin BC Complex Negatively Regulates AXL and Marks a Differentiated Phenotype in Melanoma.

High expression of the receptor tyrosine kinase AXL is implicated in epithelial-to-mesenchymal transition, cancer progression, and therapy resistance. For example, AXL is abundant in BRAF mutant melanomas progressing on targeted BRAF/MEK inhibition. Therefore, AXL is thought to represent an attractive therapeutic target. This notwithstanding, little is known about the mechanisms governing expression of AXL. Here, we describe a FACS-based whole-genome-wide CRISPR-Cas9 screen to uncover regulators of AXL expression. We identified several genes, inactivation of which led to increased AXL expression. Most remarkable was the identification of five components that associate with the Elongin BC heterodimer. Elongin B/C engage in multiple protein-protein interactions, including the transcription factor complex subunit Elongin A, the von Hippel-Lindau (VHL) tumor suppressor protein, and members of the SOCS-box protein family. The screen identified ELOB, ELOC, SOCS5, UBE2F, and RNF7, each of which we demonstrate to serve as an inhibitor of AXL expression. Although the AXL promoter contains hypoxia response elements and Elongin B/C are found in the VHL complex, Elongin B/C unexpectedly regulate AXL independently of hypoxia. Instead, we demonstrate that the Elongin BC complex interacts with AXL through ELOB, and contributes to proteasomal AXL turnover. RNA-sequencing and IHC analyses of melanoma patient-derived xenografts and clinical samples revealed a negative association between Elongin B/C and dedifferentiation. Together, the Elongin BC complex regulates AXL and marks a differentiated melanoma phenotype. IMPLICATIONS: This study identifies the Elongin BC complex as a key regulator of AXL expression and marker of melanoma differentiation.

An adverse tumor-protective effect of IDO1 inhibition.

By restoring tryptophan, indoleamine 2,3-dioxygenase 1 (IDO1) inhibitors aim to reactivate anti-tumor T cells. However, a phase III trial assessing their clinical benefit failed, prompting us to revisit the role of IDO1 in tumor cells under T cell attack. We show here that IDO1 inhibition leads to an adverse protection of melanoma cells to T cell-derived interferon-gamma (IFNγ). RNA sequencing and ribosome profiling shows that IFNγ shuts down general protein translation, which is reversed by IDO1 inhibition. Impaired translation is accompanied by an amino acid deprivation-dependent stress response driving activating transcription factor-4 (ATF4)high/microphtalmia-associated transcription factor (MITF)low transcriptomic signatures, also in patient melanomas. Single-cell sequencing analysis reveals that MITF downregulation upon immune checkpoint blockade treatment predicts improved patient outcome. Conversely, MITF restoration in cultured melanoma cells causes T cell resistance. These results highlight the critical role of tryptophan and MITF in the melanoma response to T cell-derived IFNγ and uncover an unexpected negative consequence of IDO1 inhibition.

SERPINB9 is commonly amplified and high expression in cancer cells correlates with poor immune checkpoint blockade response.

Immunotherapies, in particular immune checkpoint blockade (ICB), have improved the clinical outcome of cancer patients, although many fail to mount a durable response. Several resistance mechanisms have been identified, but our understanding of the requirements for a robust ICB response is incomplete. We have engineered an MHC I/antigen: TCR-matched panel of human NSCLC cancer and T cells to identify tumor cell-intrinsic T cell resistance mechanisms. The top differentially expressed gene in resistant tumor cells was SERPINB9. This serine protease inhibitor of the effector T cell-derived molecule granzyme B prevents caspase-mediated tumor apoptosis. Concordantly, we show that genetic ablation of SERPINB9 reverts T cell resistance of NSCLC cell lines, whereas its overexpression reduces T cell sensitivity. SERPINB9 expression in NSCLC strongly correlates with a mesenchymal phenotype. We also find that SERPINB9 is commonly amplified in cancer, particularly melanoma in which it is indicative of poor prognosis. Single-cell RNA sequencing of ICB-treated melanomas revealed that SERPINB9 expression is elevated not only in cells from post- versus pre-treatment cancers, but also in ICB-refractory cancers. In NSCLC we commonly observed rare SERPINB9-positive cancer cells, possibly accounting for reservoirs of ICB-resistant cells. While underscoring SERPINB9 as a potential target to combat immunotherapy resistance, these results suggest its potential to serve as a prognostic and predictive biomarker.

RNF31 inhibition sensitizes tumors to bystander killing by innate and adaptive immune cells.

Tumor escape mechanisms for immunotherapy include deficiencies in antigen presentation, diminishing adaptive CD8+ T cell antitumor activity. Although innate natural killer (NK) cells are triggered by loss of MHC class I, their response is often inadequate. To increase tumor susceptibility to both innate and adaptive immune elimination, we performed parallel genome-wide CRISPR-Cas9 knockout screens under NK and CD8+ T cell pressure. We identify all components, RNF31, RBCK1, and SHARPIN, of the linear ubiquitination chain assembly complex (LUBAC). Genetic and pharmacologic ablation of RNF31, an E3 ubiquitin ligase, strongly sensitizes cancer cells to NK and CD8+ T cell killing. This occurs in a tumor necrosis factor (TNF)-dependent manner, causing loss of A20 and non-canonical IKK complexes from TNF receptor complex I. A small-molecule RNF31 inhibitor sensitizes colon carcinoma organoids to TNF and greatly enhances bystander killing of MHC antigen-deficient tumor cells. These results merit exploration of RNF31 inhibition as a clinical pharmacological opportunity for immunotherapy-refractory cancers.

Predictive Immune-Checkpoint Blockade Classifiers Identify Tumors Responding to Inhibition of PD-1 and/or CTLA-4.

PURPOSE: Combining anti-PD-1 + anti-CTLA-4 immune-checkpoint blockade (ICB) shows improved patient benefit, but it is associated with severe immune-related adverse events and exceedingly high cost. Therefore, there is a dire need to predict which patients respond to monotherapy and which require combination ICB treatment. EXPERIMENTAL DESIGN: In patient-derived melanoma xenografts (PDX), human tumor microenvironment (TME) cells were swiftly replaced by murine cells upon transplantation. Using our XenofilteR deconvolution algorithm we curated human tumor cell RNA reads, which were subsequently subtracted in silico from bulk (tumor cell + TME) patients' melanoma RNA. This produced a purely tumor cell-intrinsic signature ("InTumor") and a signature comprising tumor cell-extrinsic RNA reads ("ExTumor"). RESULTS: We show that whereas the InTumor signature predicts response to anti-PD-1, the ExTumor predicts anti-CTLA-4 benefit. In PDX, InTumorLO, but not InTumorHI, tumors are effectively eliminated by cytotoxic T cells. When used in conjunction, the InTumor and ExTumor signatures identify not only patients who have a substantially higher chance of responding to combination treatment than to either monotherapy, but also those who are likely to benefit little from anti-CTLA-4 on top of anti-PD-1. CONCLUSIONS: These signatures may be exploited to distinguish melanoma patients who need combination ICB blockade from those who likely benefit from either monotherapy.

Augmenting Immunotherapy Impact by Lowering Tumor TNF Cytotoxicity Threshold.

New opportunities are needed to increase immune checkpoint blockade (ICB) benefit. Whereas the interferon (IFN)γ pathway harbors both ICB resistance factors and therapeutic opportunities, this has not been systematically investigated for IFNγ-independent signaling routes. A genome-wide CRISPR/Cas9 screen to sensitize IFNγ receptor-deficient tumor cells to CD8 T cell elimination uncovered several hits mapping to the tumor necrosis factor (TNF) pathway. Clinically, we show that TNF antitumor activity is only limited in tumors at baseline and in ICB non-responders, correlating with its low abundance. Taking advantage of the genetic screen, we demonstrate that ablation of the top hit, TRAF2, lowers the TNF cytotoxicity threshold in tumors by redirecting TNF signaling to favor RIPK1-dependent apoptosis. TRAF2 loss greatly enhanced the therapeutic potential of pharmacologic inhibition of its interaction partner cIAP, another screen hit, thereby cooperating with ICB. Our results suggest that selective reduction of the TNF cytotoxicity threshold increases the susceptibility of tumors to immunotherapy.

Absence of PD-L1 expression on tumor cells in the context of an activated immune infiltrate may indicate impaired IFNγ signaling in non-small cell lung cancer.

BACKGROUND: In non-small cell lung cancer (NSCLC), PD-L1 expression on either tumor cells (TC) or both TC and tumor-infiltrating immune cells (IC) is currently the most used biomarker in cancer immunotherapy. However, the mechanisms involved in PD-L1 regulation are not fully understood. To provide better insight in these mechanisms, a multiangular analysis approach was used to combine protein and mRNA expression with several clinicopathological characteristics. PATIENTS AND METHODS: Archival tissues from 640 early stage, resected NSCLC patients were analyzed with immunohistochemistry for expression of PD-L1 and CD8 infiltration. In addition, mutational status and expression of a selection of immune genes involved in the PD-L1/PD-1 axis and T-cell response was determined. RESULTS: Tumors with high PD-L1 expression on TC or on IC represent two subsets of NSCLC with minimal overlap. We observed that PD-L1 expression on IC irrespective of expression on TC is a good marker for inflammation within tumors. In the tumors with the highest IC expression and absent TC expression an association with reduced IFNγ downstream signaling in tumor cells was observed. CONCLUSIONS: These results show that PD-L1 expression on TC and IC are both independent hallmarks of the inflamed phenotype in NSCLC, and TC-negative/IC-high tumors can also be categorized as inflamed. The lack of correlation between PD-L1 TC and IC expression in this subgroup may be caused by impaired IFNγ signaling in tumor cells. These findings may bring a better understanding of the tumor-immune system interaction and the clinical relevance of PD-L1 expression on IC irrespective of PD-L1 expression on TC.

XenofilteR: computational deconvolution of mouse and human reads in tumor xenograft sequence data.

BACKGROUND: Mouse xenografts from (patient-derived) tumors (PDX) or tumor cell lines are widely used as models to study various biological and preclinical aspects of cancer. However, analyses of their RNA and DNA profiles are challenging, because they comprise reads not only from the grafted human cancer but also from the murine host. The reads of murine origin result in false positives in mutation analysis of DNA samples and obscure gene expression levels when sequencing RNA. However, currently available algorithms are limited and improvements in accuracy and ease of use are necessary. RESULTS: We developed the R-package XenofilteR, which separates mouse from human sequence reads based on the edit-distance between a sequence read and reference genome. To assess the accuracy of XenofilteR, we generated sequence data by in silico mixing of mouse and human DNA sequence data. These analyses revealed that XenofilteR removes > 99.9% of sequence reads of mouse origin while retaining human sequences. This allowed for mutation analysis of xenograft samples with accurate variant allele frequencies, and retrieved all non-synonymous somatic tumor mutations. CONCLUSIONS: XenofilteR accurately dissects RNA and DNA sequences from mouse and human origin, thereby outperforming currently available tools. XenofilteR is open source and available at https://github.com/PeeperLab/XenofilteR .

Cooperative targeting of melanoma heterogeneity with an AXL antibody-drug conjugate and BRAF/MEK inhibitors.

Intratumor heterogeneity is a key factor contributing to therapeutic failure and, hence, cancer lethality. Heterogeneous tumors show partial therapy responses, allowing for the emergence of drug-resistant clones that often express high levels of the receptor tyrosine kinase AXL. In melanoma, AXL-high cells are resistant to MAPK pathway inhibitors, whereas AXL-low cells are sensitive to these inhibitors, rationalizing a differential therapeutic approach. We developed an antibody-drug conjugate, AXL-107-MMAE, comprising a human AXL antibody linked to the microtubule-disrupting agent monomethyl auristatin E. We found that AXL-107-MMAE, as a single agent, displayed potent in vivo anti-tumor activity in patient-derived xenografts, including melanoma, lung, pancreas and cervical cancer. By eliminating distinct populations in heterogeneous melanoma cell pools, AXL-107-MMAE and MAPK pathway inhibitors cooperatively inhibited tumor growth. Furthermore, by inducing AXL transcription, BRAF/MEK inhibitors potentiated the efficacy of AXL-107-MMAE. These findings provide proof of concept for the premise that rationalized combinatorial targeting of distinct populations in heterogeneous tumors may improve therapeutic effect, and merit clinical validation of AXL-107-MMAE in both treatment-naive and drug-resistant cancers in mono- or combination therapy.

Cancer drug addiction is relayed by an ERK2-dependent phenotype switch.

Observations from cultured cells, animal models and patients raise the possibility that the dependency of tumours on the therapeutic drugs to which they have acquired resistance represents a vulnerability with potential applications in cancer treatment. However, for this drug addiction trait to become of clinical interest, we must first define the mechanism that underlies it. We performed an unbiased CRISPR-Cas9 knockout screen on melanoma cells that were both resistant and addicted to inhibition of the serine/threonine-protein kinase BRAF, in order to functionally mine their genome for 'addiction genes'. Here we describe a signalling pathway comprising ERK2 kinase and JUNB and FRA1 transcription factors, disruption of which allowed addicted tumour cells to survive on treatment discontinuation. This occurred in both cultured cells and mice and was irrespective of the acquired drug resistance mechanism. In melanoma and lung cancer cells, death induced by drug withdrawal was preceded by a specific ERK2-dependent phenotype switch, alongside transcriptional reprogramming reminiscent of the epithelial-mesenchymal transition. In melanoma cells, this reprogramming caused the shutdown of microphthalmia-associated transcription factor (MITF), a lineage survival oncoprotein; restoring this protein reversed phenotype switching and prevented the lethality associated with drug addiction. In patients with melanoma that had progressed during treatment with a BRAF inhibitor, treatment cessation was followed by increased expression of the receptor tyrosine kinase AXL, which is associated with the phenotype switch. Drug discontinuation synergized with the melanoma chemotherapeutic agent dacarbazine by further suppressing MITF and its prosurvival target, B-cell lymphoma 2 (BCL-2), and by inducing DNA damage in cancer cells. Our results uncover a pathway that underpins drug addiction in cancer cells, which may help to guide the use of alternating therapeutic strategies for enhanced clinical responses in drug-resistant cancers.

Transcription Factor NFIB Is a Driver of Small Cell Lung Cancer Progression in Mice and Marks Metastatic Disease in Patients.

Small cell lung cancer (SCLC) is an aggressive neuroendocrine tumor, and no effective treatment is available to date. Mouse models of SCLC based on the inactivation of Rb1 and Trp53 show frequent amplifications of the Nfib and Mycl genes. Here, we report that, although overexpression of either transcription factor accelerates tumor growth, NFIB specifically promotes metastatic spread. High NFIB levels are associated with expansive growth of a poorly differentiated and almost exclusively E-cadherin (CDH1)-negative invasive tumor cell population. Consistent with the mouse data, we find that NFIB is overexpressed in almost all tested human metastatic high-grade neuroendocrine lung tumors, warranting further assessment of NFIB as a tumor progression marker in a clinical setting.

BRAF(V600E) Kinase Domain Duplication Identified in Therapy-Refractory Melanoma Patient-Derived Xenografts.

The therapeutic landscape of melanoma is improving rapidly. Targeted inhibitors show promising results, but drug resistance often limits durable clinical responses. There is a need for in vivo systems that allow for mechanistic drug resistance studies and (combinatorial) treatment optimization. Therefore, we established a large collection of patient-derived xenografts (PDXs), derived from BRAF(V600E), NRAS(Q61), or BRAF(WT)/NRAS(WT) melanoma metastases prior to treatment with BRAF inhibitor and after resistance had occurred. Taking advantage of PDXs as a limitless source, we screened tumor lysates for resistance mechanisms. We identified a BRAF(V600E) protein harboring a kinase domain duplication (BRAF(V600E/DK)) in ∼10% of the cases, both in PDXs and in an independent patient cohort. While BRAF(V600E/DK) depletion restored sensitivity to BRAF inhibition, a pan-RAF dimerization inhibitor effectively eliminated BRAF(V600E/DK)-expressing cells. These results illustrate the utility of this PDX platform and warrant clinical validation of BRAF dimerization inhibitors for this group of melanoma patients.

Neoadjuvant tamoxifen synchronizes ERα binding and gene expression profiles related to outcome and proliferation.

Estrogen receptor alpha (ERα)-positive breast cancers are frequently treated with tamoxifen, but resistance is common. It remains elusive how tamoxifen resistance occurs and predictive biomarkers for treatment outcome are needed. Because most biomarker discovery studies are performed using pre-treatment surgical resections, the effects of tamoxifen therapy directly on the tumor cell in vivo remain unexamined. In this study, we assessed DNA copy number, gene expression profiles and ERα/chromatin binding landscapes on breast tumor specimens, both before and after neoadjuvant tamoxifen treatment. We observed neoadjuvant tamoxifen treatment synchronized ERα/chromatin interactions and downstream gene expression, indicating that hormonal therapy reduces inter-tumor molecular variability. ERα-synchronized sites are associated with dynamic FOXA1 action at these sites, which is under control of growth factor signaling. Genes associated with tamoxifen-synchronized sites are capable of differentiating patients for tamoxifen benefit. Due to the direct effects of therapeutics on ERα behavior and transcriptional output, our study highlights the added value of biomarker discovery studies after neoadjuvant drug exposure.

Preserved genetic diversity in organoids cultured from biopsies of human colorectal cancer metastases.

Tumor organoids are 3D cultures of cancer cells. They can be derived from the tumor of each individual patient, thereby providing an attractive ex vivo assay to tailor treatment. Using patient-derived tumor organoids for this purpose requires that organoids derived from biopsies maintain the genetic diversity of the in vivo tumor. In this study tumor biopsies were obtained from 14 patients with metastatic colorectal cancer (i) to test the feasibility of organoid culture from metastatic biopsy specimens and (ii) to compare the genetic diversity of patient-derived tumor organoids and the original tumor biopsy. Genetic analysis was performed using SOLiD sequencing for 1,977 cancer-relevant genes. Copy number profiles were generated from sequencing data using CopywriteR. Here we demonstrate that organoid cultures can be established from tumor biopsies of patients with metastatic colorectal cancer with a success rate of 71%. Genetic analysis showed that organoids reflect the metastasis from which they were derived. Ninety percent of somatic mutations were shared between organoids and biopsies from the same patient, and the DNA copy number profiles of organoids and the corresponding original tumor show a correlation of 0.89. Most importantly, none of the mutations that were found exclusively in either the tumor or organoid culture are in driver genes or genes amenable for drug targeting. These findings support further exploration of patient-derived organoids as an ex vivo platform to personalize anticancer treatment.

Intra- and inter-tumor heterogeneity in a vemurafenib-resistant melanoma patient and derived xenografts.

The development of targeted inhibitors, like vemurafenib, has greatly improved the clinical outcome of BRAF(V600E) metastatic melanoma. However, resistance to such compounds represents a formidable problem. Using whole-exome sequencing and functional analyses, we have investigated the nature and pleiotropy of vemurafenib resistance in a melanoma patient carrying multiple drug-resistant metastases. Resistance was caused by a plethora of mechanisms, all of which reactivated the MAPK pathway. In addition to three independent amplifications and an aberrant form of BRAF(V600E), we identified a new activating insertion in MEK1. This MEK1(T55delins) (RT) mutation could be traced back to a fraction of the pre-treatment lesion and not only provided protection against vemurafenib but also promoted local invasion of transplanted melanomas. Analysis of patient-derived xenografts (PDX) from therapy-refractory metastases revealed that multiple resistance mechanisms were present within one metastasis. This heterogeneity, both inter- and intra-tumorally, caused an incomplete capture in the PDX of the resistance mechanisms observed in the patient. In conclusion, vemurafenib resistance in a single patient can be established through distinct events, which may be preexisting. Furthermore, our results indicate that PDX may not harbor the full genetic heterogeneity seen in the patient's melanoma.

Robust BRCA1-like classification of copy number profiles of samples repeated across different datasets and platforms.

Breast cancers with BRCA1 germline mutation have a characteristic DNA copy number (CN) pattern. We developed a test that assigns CN profiles to be 'BRCA1-like' or 'non-BRCA1-like', which refers to resembling a BRCA1-mutated tumor or resembling a tumor without a BRCA1 mutation, respectively. Approximately one third of the BRCA1-like breast cancers have a BRCA1 mutation, one third has hypermethylation of the BRCA1 promoter and one third has an unknown reason for being BRCA1-like. This classification is indicative of patients' response to high dose alkylating and platinum containing chemotherapy regimens, which targets the inability of BRCA1 deficient cells to repair DNA double strand breaks. We investigated whether this classification can be reliably obtained with next generation sequencing and copy number platforms other than the bacterial artificial chromosome (BAC) array Comparative Genomic Hybridization (aCGH) on which it was originally developed. We investigated samples from 230 breast cancer patients for which a CN profile had been generated on two to five platforms, comprising low coverage CN sequencing, CN extraction from targeted sequencing panels (CopywriteR), Affymetrix SNP6.0, 135K/720K oligonucleotide aCGH, Affymetrix Oncoscan FFPE (MIP) technology, 3K BAC and 32K BAC aCGH. Pairwise comparison of genomic position-mapped profiles from the original aCGH platform and other platforms revealed concordance. For most cases, biological differences between samples exceeded the differences between platforms within one sample. We observed the same classification across different platforms in over 80% of the patients and kappa values of at least 0.36. Differential classification could be attributed to CN profiles that were not strongly associated to one class. In conclusion, we have shown that the genomic regions that define our BRCA1-like classifier are robustly measured by different CN profiling technologies, providing the possibility to retro- and prospectively investigate BRCA1-like classification across a wide range of CN platforms.

CopywriteR: DNA copy number detection from off-target sequence data.

Current methods for detection of copy number variants (CNV) and aberrations (CNA) from targeted sequencing data are based on the depth of coverage of captured exons. Accurate CNA determination is complicated by uneven genomic distribution and non-uniform capture efficiency of targeted exons. Here we present CopywriteR, which eludes these problems by exploiting 'off-target' sequence reads. CopywriteR allows for extracting uniformly distributed copy number information, can be used without reference, and can be applied to sequencing data obtained from various techniques including chromatin immunoprecipitation and target enrichment on small gene panels. CopywriteR outperforms existing methods and constitutes a widely applicable alternative to available tools.

Nur1 dephosphorylation confers positive feedback to mitotic exit phosphatase activation in budding yeast.

Substrate dephosphorylation by the cyclin-dependent kinase (Cdk)-opposing phosphatase, Cdc14, is vital for many events during budding yeast mitotic exit. Cdc14 is sequestered in the nucleolus through inhibitory binding to Net1, from which it is released in anaphase following Net1 phosphorylation. Initial Net1 phosphorylation depends on Cdk itself, in conjunction with proteins of the Cdc14 Early Anaphase Release (FEAR) network. Later on, the Mitotic Exit Network (MEN) signaling cascade maintains Cdc14 release. An important unresolved question is how Cdc14 activity can increase in early anaphase, while Cdk activity, that is required for Net1 phosphorylation, decreases and the MEN is not yet active. Here we show that the nuclear rim protein Nur1 interacts with Net1 and, in its Cdk phosphorylated form, inhibits Cdc14 release. Nur1 is dephosphorylated by Cdc14 in early anaphase, relieving the inhibition and promoting further Cdc14 release. Nur1 dephosphorylation thus describes a positive feedback loop in Cdc14 phosphatase activation during mitotic exit, required for faithful chromosome segregation and completion of the cell division cycle.