BRAF (v-raf murine sarcoma viral oncogene homolog B) inhibitors elicit a transient anti-tumor response in ∼ 80% of BRAF(V600)-mutant melanoma patients that almost uniformly precedes the emergence of resistance. Here we used a mouse model of melanoma in which melanocyte-specific expression of Braf(V618E) (analogous to the human BRAF(V600E) mutation) led to the development of skin hyperpigmentation and nevi, as well as melanoma formation with incomplete penetrance. Sleeping Beauty insertional mutagenesis in this model led to accelerated and fully penetrant melanomagenesis and synchronous tumor formation. Treatment of Braf(V618E) transposon mice with the BRAF inhibitor PLX4720 resulted in tumor regression followed by relapse. Analysis of transposon insertions identified eight genes including Braf, Mitf, and ERas (ES-cell expressed Ras) as candidate resistance genes. Expression of ERAS in human melanoma cell lines conferred resistance to PLX4720 and induced hyperphosphorylation of AKT (v-akt murine thymoma viral oncogene homolog 1), a phenotype reverted by combinatorial treatment with PLX4720 and the AKT inhibitor MK2206. We show that ERAS expression elicits a prosurvival signal associated with phosphorylation/inactivation of BAD, and that the resistance of hepatocyte growth factor-treated human melanoma cells to PLX4720 can be reverted by treatment with the BAD-like BH3 mimetic ABT-737. Thus, we define a role for the AKT/BAD pathway in resistance to BRAF inhibition and illustrate an in vivo approach for finding drug resistance genes.

In a systematic sequencing screen of the coding exons of the X chromosome in 250 families with X-linked mental retardation (XLMR), we identified two nonsense mutations and one consensus splice-site mutation in the AP1S2 gene on Xp22 in three families. Affected individuals in these families showed mild-to-profound mental retardation. Other features included hypotonia early in life and delay in walking. AP1S2 encodes an adaptin protein that constitutes part of the adaptor protein complex found at the cytoplasmic face of coated vesicles located at the Golgi complex. The complex mediates the recruitment of clathrin to the vesicle membrane. Aberrant endocytic processing through disruption of adaptor protein complexes is likely to result from the AP1S2 mutations identified in the three XLMR-affected families, and such defects may plausibly cause abnormal synaptic development and function. AP1S2 is the first reported XLMR gene that encodes a protein directly involved in the assembly of endocytic vesicles.

Transmissible cancers are clonal lineages that spread through populations via contagious cancer cells. Although rare in nature, two facial tumor clones affect Tasmanian devils. Here we perform comparative genetic and functional characterization of these lineages. The two cancers have similar patterns of mutation and show no evidence of exposure to exogenous mutagens or viruses. Genes encoding PDGF receptors have copy number gains and are present on extrachromosomal double minutes. Drug screening indicates causative roles for receptor tyrosine kinases and sensitivity to inhibitors of DNA repair. Y chromosome loss from a male clone infecting a female host suggests immunoediting. These results imply that Tasmanian devils may have inherent susceptibility to transmissible cancers and present a suite of therapeutic compounds for use in conservation.

Combinations of anti-cancer drugs can overcome resistance and provide new treatments1,2. The number of possible drug combinations vastly exceeds what could be tested clinically. Efforts to systematically identify active combinations and the tissues and molecular contexts in which they are most effective could accelerate the development of combination treatments. Here we evaluate the potency and efficacy of 2,025 clinically relevant two-drug combinations, generating a dataset encompassing 125 molecularly characterized breast, colorectal and pancreatic cancer cell lines. We show that synergy between drugs is rare and highly context-dependent, and that combinations of targeted agents are most likely to be synergistic. We incorporate multi-omic molecular features to identify combination biomarkers and specify synergistic drug combinations and their active contexts, including in basal-like breast cancer, and microsatellite-stable or KRAS-mutant colon cancer. Our results show that irinotecan and CHEK1 inhibition have synergistic effects in microsatellite-stable or KRAS-TP53 double-mutant colon cancer cells, leading to apoptosis and suppression of tumour xenograft growth. This study identifies clinically relevant effective drug combinations in distinct molecular subpopulations and is a resource to guide rational efforts to develop combinatorial drug treatments.

We have identified one frameshift mutation, one splice-site mutation, and two missense mutations in highly conserved residues in ZDHHC9 at Xq26.1 in 4 of 250 families with X-linked mental retardation (XLMR). In three of the families, the mental retardation phenotype is associated with a Marfanoid habitus, although none of the affected individuals meets the Ghent criteria for Marfan syndrome. ZDHHC9 is a palmitoyltransferase that catalyzes the posttranslational modification of NRAS and HRAS. The degree of palmitoylation determines the temporal and spatial location of these proteins in the plasma membrane and Golgi complex. The finding of mutations in ZDHHC9 suggests that alterations in the concentrations and cellular distribution of target proteins are sufficient to cause disease. This is the first XLMR gene to be reported that encodes a posttranslational modification enzyme, palmitoyltransferase. Furthermore, now that the first palmitoyltransferase that causes mental retardation has been identified, defects in other palmitoylation transferases become good candidates for causing other mental retardation syndromes.

Drug resistance is an almost inevitable consequence of cancer therapy and ultimately proves fatal for the majority of patients. In many cases, this is the consequence of specific gene mutations that have the potential to be targeted to resensitize the tumor. The ability to uniformly saturate the genome with point mutations without chromosome or nucleotide sequence context bias would open the door to identify all putative drug resistance mutations in cancer models. Here, we describe such a method for elucidating drug resistance mechanisms using genome-wide chemical mutagenesis allied to next-generation sequencing. We show that chemically mutagenizing the genome of cancer cells dramatically increases the number of drug-resistant clones and allows the detection of both known and novel drug resistance mutations. We used an efficient computational process that allows for the rapid identification of involved pathways and druggable targets. Such a priori knowledge would greatly empower serial monitoring strategies for drug resistance in the clinic as well as the development of trials for drug-resistant patients.

Multivalent second-generation TRAIL-R2 agonists are currently in late preclinical development and early clinical trials. Herein, we use a representative second-generation agent, MEDI3039, to address two major clinical challenges facing these agents: lack of predictive biomarkers to enable patient selection and emergence of resistance. Genome-wide CRISPR knockout screens were notable for the lack of resistance mechanisms beyond the canonical TRAIL-R2 pathway (caspase-8, FADD, BID) as well as p53 and BAX in TP53 wild-type models, whereas a CRISPR activatory screen identified cell death inhibitors MCL-1 and BCL-XL as mechanisms to suppress MEDI3039-induced cell death. High-throughput drug screening failed to identify genomic alterations associated with response to MEDI3039; however, transcriptomics analysis revealed striking association between MEDI3039 sensitivity and expression of core components of the extrinsic apoptotic pathway, most notably its main apoptotic effector caspase-8 in solid tumor cell lines. Further analyses of colorectal cell lines and patient-derived xenografts identified caspase-8 expression ratio to its endogenous regulator FLIP(L) as predictive of sensitivity to MEDI3039 in several major solid tumor types and a further subset indicated by caspase-8:MCL-1 ratio. Subsequent MEDI3039 combination screening of TRAIL-R2, caspase-8, FADD, and BID knockout models with 60 compounds with varying mechanisms of action identified two inhibitor of apoptosis proteins (IAP) that exhibited strong synergy with MEDI3039 that could reverse resistance only in BID-deleted models. In summary, we identify the ratios of caspase-8:FLIP(L) and caspase-8:MCL-1 as potential predictive biomarkers for second-generation TRAIL-R2 agonists and loss of key effectors such as FADD and caspase-8 as likely drivers of clinical resistance in solid tumors.

Interferon-γ (IFN-γ) signaling mediates host responses to infection, inflammation and anti-tumor immunity. Mutations in the IFN-γ signaling pathway cause immunological disorders, hematological malignancies, and resistance to immune checkpoint blockade (ICB) in cancer; however, the function of most clinically observed variants remains unknown. Here, we systematically investigate the genetic determinants of IFN-γ response in colorectal cancer cells using CRISPR-Cas9 screens and base editing mutagenesis. Deep mutagenesis of JAK1 with cytidine and adenine base editors, combined with pathway-wide screens, reveal loss-of-function and gain-of-function mutations, including causal variants in hematological malignancies and mutations detected in patients refractory to ICB. We functionally validate variants of uncertain significance in primary tumor organoids, where engineering missense mutations in JAK1 enhanced or reduced sensitivity to autologous tumor-reactive T cells. We identify more than 300 predicted missense mutations altering IFN-γ pathway activity, generating a valuable resource for interpreting gene variant function.

Malignant mesothelioma (MM) is poorly responsive to systemic cytotoxic chemotherapy and invariably fatal. Here we describe a screen of 94 drugs in 15 exome-sequenced MM lines and the discovery of a subset defined by loss of function of the nuclear deubiquitinase BRCA associated protein-1 (BAP1) that demonstrate heightened sensitivity to TRAIL (tumour necrosis factor-related apoptosis-inducing ligand). This association is observed across human early passage MM cultures, mouse xenografts and human tumour explants. We demonstrate that BAP1 deubiquitinase activity and its association with ASXL1 to form the Polycomb repressive deubiquitinase complex (PR-DUB) impacts TRAIL sensitivity implicating transcriptional modulation as an underlying mechanism. Death receptor agonists are well-tolerated anti-cancer agents demonstrating limited therapeutic benefit in trials without a targeting biomarker. We identify BAP1 loss-of-function mutations, which are frequent in MM, as a potential genomic stratification tool for TRAIL sensitivity with immediate and actionable therapeutic implications.

To develop an international consensus on managing penile cancer patients during the COVID-19 acute waves. A major concern for patients with penile cancer during the coronavirus disease 2019 (COVID-19) pandemic is how the enforced safety measures will affect their disease management. Delays in diagnosis and treatment initiation may have an impact on the extent of the primary lesion as well as the cancer-specific survival because of the development and progression of inguinal lymph node metastases.

Given the immune system's importance for cancer surveillance and treatment, we have investigated how it may be affected by SARS-CoV-2 infection of cancer patients. Across some heterogeneity in tumor type, stage, and treatment, virus-exposed solid cancer patients display a dominant impact of SARS-CoV-2, apparent from the resemblance of their immune signatures to those for COVID-19+ non-cancer patients. This is not the case for hematological malignancies, with virus-exposed patients collectively displaying heterogeneous humoral responses, an exhausted T cell phenotype and a high prevalence of prolonged virus shedding. Furthermore, while recovered solid cancer patients' immunophenotypes resemble those of non-virus-exposed cancer patients, recovered hematological cancer patients display distinct, lingering immunological legacies. Thus, while solid cancer patients, including those with advanced disease, seem no more at risk of SARS-CoV-2-associated immune dysregulation than the general population, hematological cancer patients show complex immunological consequences of SARS-CoV-2 exposure that might usefully inform their care.

Organoid cell culture methodologies are enabling the generation of cell models from healthy and diseased tissue. Patient-derived cancer organoids that recapitulate the genetic and histopathological diversity of patient tumours are being systematically generated, providing an opportunity to investigate new cancer biology and therapeutic approaches. The use of organoid cultures for many applications, including genetic and chemical perturbation screens, is limited due to the technical demands and cost associated with their handling and propagation. Here we report and benchmark a suspension culture technique for cancer organoids which allows for the expansion of models to tens of millions of cells with increased efficiency in comparison to standard organoid culturing protocols. Using whole-genome DNA and RNA sequencing analyses, as well as medium-throughput drug sensitivity testing and genome-wide CRISPR-Cas9 screening, we demonstrate that cancer organoids grown as a suspension culture are genetically and phenotypically similar to their counterparts grown in standard conditions. This culture technique simplifies organoid cell culture and extends the range of organoid applications, including for routine use in large-scale perturbation screens.

A major challenge in cancer genetics is to determine which low-frequency somatic mutations are drivers of tumorigenesis. Here we interrogate the genomes of 7,651 diverse human cancers and find inactivating mutations in the homeodomain transcription factor gene CUX1 (cut-like homeobox 1) in ~1-5% of various tumors. Meta-analysis of CUX1 mutational status in 2,519 cases of myeloid malignancies reveals disruptive mutations associated with poor survival, highlighting the clinical significance of CUX1 loss. In parallel, we validate CUX1 as a bona fide tumor suppressor using mouse transposon-mediated insertional mutagenesis and Drosophila cancer models. We demonstrate that CUX1 deficiency activates phosphoinositide 3-kinase (PI3K) signaling through direct transcriptional downregulation of the PI3K inhibitor PIK3IP1 (phosphoinositide-3-kinase interacting protein 1), leading to increased tumor growth and susceptibility to PI3K-AKT inhibition. Thus, our complementary approaches identify CUX1 as a pan-driver of tumorigenesis and uncover a potential strategy for treating CUX1-mutant tumors.

Clear cell renal cell carcinoma (ccRCC) is the most common form of adult kidney cancer, characterized by the presence of inactivating mutations in the VHL gene in most cases, and by infrequent somatic mutations in known cancer genes. To determine further the genetics of ccRCC, we have sequenced 101 cases through 3,544 protein-coding genes. Here we report the identification of inactivating mutations in two genes encoding enzymes involved in histone modification-SETD2, a histone H3 lysine 36 methyltransferase, and JARID1C (also known as KDM5C), a histone H3 lysine 4 demethylase-as well as mutations in the histone H3 lysine 27 demethylase, UTX (KMD6A), that we recently reported. The results highlight the role of mutations in components of the chromatin modification machinery in human cancer. Furthermore, NF2 mutations were found in non-VHL mutated ccRCC, and several other probable cancer genes were identified. These results indicate that substantial genetic heterogeneity exists in a cancer type dominated by mutations in a single gene, and that systematic screens will be key to fully determining the somatic genetic architecture of cancer.

Castration-resistant prostate cancer (CRPC) is the major cause of death from prostate cancer. Biomarkers to improve early detection and prediction of CRPC especially using non-invasive liquid biopsies could improve outcomes. Therefore, we investigated the plasma exosomal miRNAs associated with CRPC and their potential for development into non-invasive early detection biomarkers for resistance to treatment. RNA-sequencing, which generated approximately five million reads per patient, was performed to identify differentially expressed plasma exosomal miRNAs in 24 treatment-naive prostate cancer and 24 CRPC patients. RT-qPCR was used to confirm the differential expressions of six exosomal miRNAs, miR-423-3p, miR-320a, miR-99a-5p, miR-320d, miR-320b, and miR-150-5p (p = 7.3 × 10-8, 0.0020, 0.018, 0.0028, 0.0013, and 0.0058, respectively) firstly in a validation cohort of 108 treatment-naive prostate cancer and 42 CRPC patients. The most significant differentially expressed miRNA, miR-423-3p, was shown to be associated with CRPC with area under the ROC curve (AUC) = 0.784. Combining miR-423-3p with prostate-specific antigen (PSA) enhanced the prediction of CRPC (AUC = 0.908). A separate research center validation with 30 treatment-naive and 30 CRPC patients also confirmed the differential expression of miR-423-3p (p = 0.016). Finally, plasma exosomal miR-423-3p expression in CRPC patients was compared to 36 non-CRPC patients under androgen depletion therapy, which showed significantly higher expression in CRPC than treated non-CRPC patients (p < 0.0001) with AUC = 0.879 to predict CRPC with no difference between treatment-naive and treated non-CRPC patients. Therefore, our findings demonstrate that a number of plasma exosomal miRNAs are associated with CRPC and miR-423-3p may serve as a biomarker for early detection/prediction of castration-resistance.

The proteome provides unique insights into disease biology beyond the genome and transcriptome. A lack of large proteomic datasets has restricted the identification of new cancer biomarkers. Here, proteomes of 949 cancer cell lines across 28 tissue types are analyzed by mass spectrometry. Deploying a workflow to quantify 8,498 proteins, these data capture evidence of cell-type and post-transcriptional modifications. Integrating multi-omics, drug response, and CRISPR-Cas9 gene essentiality screens with a deep learning-based pipeline reveals thousands of protein biomarkers of cancer vulnerabilities that are not significant at the transcript level. The power of the proteome to predict drug response is very similar to that of the transcriptome. Further, random downsampling to only 1,500 proteins has limited impact on predictive power, consistent with protein networks being highly connected and co-regulated. This pan-cancer proteomic map (ProCan-DepMapSanger) is a comprehensive resource available at https://cellmodelpassports.sanger.ac.uk.

Nasopharyngeal carcinoma (NPC), a cancer that is etiologically associated with the Epstein-Barr virus (EBV), is endemic in Southern China and Southeast Asia. The scarcity of representative NPC cell lines owing to the frequent loss of EBV episomes following prolonged propagation and compromised authenticity of previous models underscores the critical need for new EBV-positive NPC models. Herein, we describe the establishment of a new EBV-positive NPC cell line, designated NPC268 from a primary non-keratinizing, differentiated NPC tissue. NPC268 can undergo productive lytic reactivation of EBV and is highly tumorigenic in immunodeficient mice. Whole-genome sequencing revealed close similarities with the tissue of origin, including large chromosomal rearrangements, while whole-genome bisulfite sequencing and RNA sequencing demonstrated a hypomethylated genome and enrichment in immune-related pathways, respectively. Drug screening of NPC268 together with six other NPC cell lines using 339 compounds, representing the largest high-throughput drug testing in NPC, revealed biomarkers associated with specific drug classes. NPC268 represents the first and only available EBV-positive non-keratinizing differentiated NPC model, and extensive genomic, methylomic, transcriptomic, and drug response data should facilitate research in EBV and NPC, where current models are limited.

Systematic studies of cancer genomes have provided unprecedented insights into the molecular nature of cancer. Using this information to guide the development and application of therapies in the clinic is challenging. Here, we report how cancer-driven alterations identified in 11,289 tumors from 29 tissues (integrating somatic mutations, copy number alterations, DNA methylation, and gene expression) can be mapped onto 1,001 molecularly annotated human cancer cell lines and correlated with sensitivity to 265 drugs. We find that cell lines faithfully recapitulate oncogenic alterations identified in tumors, find that many of these associate with drug sensitivity/resistance, and highlight the importance of tissue lineage in mediating drug response. Logic-based modeling uncovers combinations of alterations that sensitize to drugs, while machine learning demonstrates the relative importance of different data types in predicting drug response. Our analysis and datasets are rich resources to link genotypes with cellular phenotypes and to identify therapeutic options for selected cancer sub-populations.

Clinical responses to anticancer therapies are often restricted to a subset of patients. In some cases, mutated cancer genes are potent biomarkers for responses to targeted agents. Here, to uncover new biomarkers of sensitivity and resistance to cancer therapeutics, we screened a panel of several hundred cancer cell lines--which represent much of the tissue-type and genetic diversity of human cancers--with 130 drugs under clinical and preclinical investigation. In aggregate, we found that mutated cancer genes were associated with cellular response to most currently available cancer drugs. Classic oncogene addiction paradigms were modified by additional tissue-specific or expression biomarkers, and some frequently mutated genes were associated with sensitivity to a broad range of therapeutic agents. Unexpected relationships were revealed, including the marked sensitivity of Ewing's sarcoma cells harbouring the EWS (also known as EWSR1)-FLI1 gene translocation to poly(ADP-ribose) polymerase (PARP) inhibitors. By linking drug activity to the functional complexity of cancer genomes, systematic pharmacogenomic profiling in cancer cell lines provides a powerful biomarker discovery platform to guide rational cancer therapeutic strategies.

We have identified three truncating, two splice-site, and three missense variants at conserved amino acids in the CUL4B gene on Xq24 in 8 of 250 families with X-linked mental retardation (XLMR). During affected subjects' adolescence, a syndrome emerged with delayed puberty, hypogonadism, relative macrocephaly, moderate short stature, central obesity, unprovoked aggressive outbursts, fine intention tremor, pes cavus, and abnormalities of the toes. This syndrome was first described by Cazebas et al., in a family that was included in our study and that carried a CUL4B missense variant. CUL4B is a ubiquitin E3 ligase subunit implicated in the regulation of several biological processes, and CUL4B is the first XLMR gene that encodes an E3 ubiquitin ligase. The relatively high frequency of CUL4B mutations in this series indicates that it is one of the most commonly mutated genes underlying XLMR and suggests that its introduction into clinical diagnostics should be a high priority.