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.

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.

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.

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.

Accurate staging of colorectal cancer (CRC) with clinicopathological parameters is important for predicting prognosis and guiding treatment but provides no information about organ site of metastases. Patterns of genomic aberrations in primary colorectal tumors may reveal a chromosomal signature for organ specific metastases.

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.

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.

Personalized medicine for cancer patients requires a deep understanding of the underlying genetics that drive cancer and the subsequent identification of predictive biomarkers. To discover new genes and pathways contributing to oncogenesis and therapy resistance in HER2+ breast cancer, we performed Mouse Mammary Tumor Virus (MMTV)-induced insertional mutagenesis screens in ErbB2/cNeu-transgenic mouse models. The screens revealed 34 common integration sites (CIS) in mammary tumors of MMTV-infected mice, highlighting loci with multiple independent MMTV integrations in which potential oncogenes are activated, most of which had never been reported as MMTV CIS. The CIS most strongly associated with the ErbB2-transgenic genotype was the locus containing Eras (ES cell-expressed Ras), a constitutively active RAS-family GTPase. We show that upon expression, Eras acts as a potent oncogenic driver through hyperactivation of the PI3K/AKT pathway, in contrast to other RAS proteins that signal primarily via the MAPK/ERK pathway and require upstream activation or activating mutations to induce signaling. We additionally show that ERAS synergistically enhances HER2-induced tumorigenesis and, in this role, can functionally replace ERBB3/HER3 by acting as a more powerful activator of PI3K/AKT signaling. Although previously reported as pseudogene in humans, we observed ERAS RNA and protein expression in a substantial subset of human primary breast carcinomas. Importantly, we show that ERAS induces primary resistance to the widely used HER2-targeting drugs Trastuzumab (Herceptin) and Lapatinib (Tykerb/Tyverb) in vivo, and is involved in acquired resistance via selective upregulation during treatment in vitro, indicating that ERAS may serve as a novel clinical biomarker for PI3K/AKT pathway hyperactivation and HER2-targeted therapy resistance.

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.

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.

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.

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.

Previously, we developed breast cancer BRCA1-like and BRCA2-like copy-number profile shrunken centroid classifiers predictive for mutation status and response to therapy, targeting homologous recombination deficiency (HRD). Therefore, we investigated BRCA1- and BRCA2-like classification in ovarian cancer, aiming to acquire classifiers with similar properties as those in breast cancer.Experimental Design: We analyzed DNA copy-number profiles of germline BRCA1- and BRCA2-mutant ovarian cancers and control tumors and observed that existing breast cancer classifiers did not sufficiently predict mutation status. Hence, we trained new shrunken centroid classifiers on this set and validated them in the independent The Cancer Genome Atlas dataset. Subsequently, we assessed BRCA1/2-like classification and obtained germline and tumor mutation and methylation status of cancer predisposition genes, among them several involved in HR repair, of 300 ovarian cancer samples derived from the consecutive cohort trial AGO-TR1 (NCT02222883).

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.

Although platinum-based drugs are widely used chemotherapeutics for cancer treatment, the determinants of tumor cell responsiveness remain poorly understood. We show that the loss of subunits LRRC8A and LRRC8D of the heteromeric LRRC8 volume-regulated anion channels (VRACs) increased resistance to clinically relevant cisplatin/carboplatin concentrations. Under isotonic conditions, about 50% of cisplatin uptake depended on LRRC8A and LRRC8D, but neither on LRRC8C nor on LRRC8E. Cell swelling strongly enhanced LRRC8-dependent cisplatin uptake, bolstering the notion that cisplatin enters cells through VRAC. LRRC8A disruption also suppressed drug-induced apoptosis independently from drug uptake, possibly by impairing VRAC-dependent apoptotic cell volume decrease. Hence, by mediating cisplatin uptake and facilitating apoptosis, VRAC plays a dual role in the cellular drug response. Incorporation of the LRRC8D subunit into VRAC substantially increased its permeability for cisplatin and the cellular osmolyte taurine, indicating that LRRC8 proteins form the channel pore. Our work suggests that LRRC8D-containing VRACs are crucial for cell volume regulation by an important organic osmolyte and may influence cisplatin/carboplatin responsiveness of tumors.

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.

Oncogene-induced cellular senescence (OIS) is emerging as a potent cancer-protective response to oncogenic events, serving to eliminate early neoplastic cells from the proliferative pool. Using combined genetic and bioinformatic analysis, we find that OIS is linked specifically to the activation of an inflammatory transcriptome. Induced genes included the pleiotropic cytokine interleukin-6 (IL-6), which upon secretion by senescent cells acted mitogenically in a paracrine fashion. Unexpectedly, IL-6 was also required for the execution of OIS, but in a cell-autonomous mode. Its depletion caused the inflammatory network to collapse and abolished senescence entry and maintenance. Furthermore, we demonstrate that the transcription factor C/EBPbeta cooperates with IL-6 to amplify the activation of the inflammatory network, including IL-8. In human colon adenomas, IL-8 specifically colocalized with arrested, p16(INK4A)-positive epithelium. We propose a model in which the context-dependent cytostatic and promitogenic functions of specific interleukins contribute to connect senescence with an inflammatory phenotype and cancer.

The final event of the eukaryotic cell cycle is cytokinesis, when two new daughter cells are born. How the timing and execution of cytokinesis is controlled is poorly understood. Here, we show that downregulation of cyclin-dependent kinase (Cdk) activity, together with upregulation of its counteracting phosphatase Cdc14, controls each of the sequential steps of cytokinesis, including furrow ingression, membrane resolution and cell separation in budding yeast. We use phosphoproteome analysis of mitotic exit to identify Cdk targets that are dephosphorylated at the time of cytokinesis. We then apply a new and widely applicable tool to generate conditionally phosphorylated proteins to identify those whose dephosphorylation is required for cytokinesis. This approach identifies Aip1, Ede1 and Inn1 as cytokinetic regulators. Our results suggest that cytokinesis is coordinately controlled by the master cell cycle regulator Cdk together with its counteracting phosphatase and that it is executed by concerted dephosphorylation of Cdk targets involved in several cell biological processes.

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.