A Landscape of Pharmacogenomic Interactions in Cancer.

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.

Effective drug combinations in breast, colon and pancreatic cancer cells.

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.

Cell Model Passports-a hub for clinical, genetic and functional datasets of preclinical cancer models.

In vitro cancer cell cultures are facile experimental models used widely for research and drug development. Many cancer cell lines are available and efforts are ongoing to derive new models representing the histopathological and molecular diversity of tumours. Cell models have been generated by multiple laboratories over decades and consequently their annotation is incomplete and inconsistent. Furthermore, the relationships between many patient-matched and derivative cell lines have been lost, and accessing information and datasets is time-consuming and difficult. Here, we describe the Cell Model Passports database; cellmodelpassports.sanger.ac.uk, which provides details of cell model relationships, patient and clinical information, as well as access to associated genetic and functional datasets. The Passports database currently contains curated details and standardized annotation for >1200 cell models, including cancer organoid cultures. The Passports will be updated with newly derived cell models and datasets as they are generated. Users can navigate the database via tissue, cancer-type, genetic feature and data availability to select a model most suitable for specific applications. A flexible REST-API provides programmatic data access and exploration. The Cell Model Passports are a valuable tool enabling access to high-dimensional genomic and phenotypic cancer cell model datasets empowering diverse research applications.

Genome-wide chemical mutagenesis screens allow unbiased saturation of the cancer genome and identification of drug resistance mutations.

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.

Systematic identification of genomic markers of drug sensitivity in cancer cells.

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.

Systematic sequencing of renal carcinoma reveals inactivation of histone modifying genes.

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.

Mutations in UPF3B, a member of the nonsense-mediated mRNA decay complex, cause syndromic and nonsyndromic mental retardation.

Nonsense-mediated mRNA decay (NMD) is of universal biological significance. It has emerged as an important global RNA, DNA and translation regulatory pathway. By systematically sequencing 737 genes (annotated in the Vertebrate Genome Annotation database) on the human X chromosome in 250 families with X-linked mental retardation, we identified mutations in the UPF3 regulator of nonsense transcripts homolog B (yeast) (UPF3B) leading to protein truncations in three families: two with the Lujan-Fryns phenotype and one with the FG phenotype. We also identified a missense mutation in another family with nonsyndromic mental retardation. Three mutations lead to the introduction of a premature termination codon and subsequent NMD of mutant UPF3B mRNA. Protein blot analysis using lymphoblastoid cell lines from affected individuals showed an absence of the UPF3B protein in two families. The UPF3B protein is an important component of the NMD surveillance machinery. Our results directly implicate abnormalities of NMD in human disease and suggest at least partial redundancy of NMD pathways.

High TNF and NF-ĸB pathway dependency are associated with AZD5582 sensitivity in OSCC via CASP8-dependent apoptosis.

Drug repurposing can accelerate the development of novel therapeutic strategies. Oral squamous cell carcinoma (OSCC) is prevalent in Asians and is associated with a poor prognosis. Here, we integrate newly generated sensitivity data for 339 anticancer drugs with genome-wide CRISPR‒Cas9 gene essentiality data for 21 cell lines mostly derived from Asian OSCC patients to identify drug repurposing candidates. We observe potent activity of AZD5582 (an antagonist of the inhibitor of apoptosis (IAP) family of proteins) in a subset of OSCC cells, which was associated with dependency on NF-ĸB pathway genes (RNF31, MAP3K7 and IKBKG). We confirm the on-target specificity and potency of AZD5582 in vitro and its ability to suppress cell growth in xenograft models. Furthermore, we identify tumor necrosis factor (TNF) as a key mediator of AZD5582 sensitivity in OSCC. CASP8-dependent apoptotic and CASP8-independent necroptotic cellular programs mediate AZD5582-induced cell death. In summary, through the systematic integration of pharmacological and CRISPR data, we identified a subset of OSCC with potent sensitivity to AZD5582 mediated through the NF-κB and TNF signalling pathways.

A pan-cancer screen identifies drug combination benefit in cancer cell lines at the individual and population level.

Combining drugs can enhance their clinical efficacy, but the number of possible combinations and inter-tumor heterogeneity make identifying effective combinations challenging, while existing approaches often overlook clinically relevant activity. We screen one of the largest cell line panels (N = 757) with 51 clinically relevant combinations and identify responses at the level of individual cell lines and tissue populations. We establish three response classes to model cellular effects beyond monotherapy: synergy, Bliss additivity, and independent drug action (IDA). Synergy is rare (11% of responses) and frequently efficacious (>50% viability reduction), whereas Bliss and IDA are more frequent but less frequently efficacious. We introduce "efficacious combination benefit" (ECB) to describe high-efficacy responses classified as either synergy, Bliss, or IDA. We identify ECB biomarkers in vitro and show that ECB predicts response in patient-derived xenografts better than synergy alone. Our work here provides a valuable resource and framework for preclinical evaluation and the development of combination treatments.

Establishment and Characterization of an Epstein-Barr Virus-positive Cell Line from a Non-keratinizing Differentiated Primary Nasopharyngeal Carcinoma.

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. SIGNIFICANCE: NPC268 is the first and only EBV-positive cell line derived from a primary non-keratinizing, differentiated nasopharyngeal carcinoma, an understudied but important subtype in Southeast Asian countries. This model adds to the limited number of authentic EBV-positive lines globally that will facilitate mechanistic studies and drug development for NPC.

Large-scale Pan-cancer Cell Line Screening Identifies Actionable and Effective Drug Combinations.

Oncology drug combinations can improve therapeutic responses and increase treatment options for patients. The number of possible combinations is vast and responses can be context-specific. Systematic screens can identify clinically relevant, actionable combinations in defined patient subtypes. We present data for 109 anticancer drug combinations from AstraZeneca's oncology small molecule portfolio screened in 755 pan-cancer cell lines. Combinations were screened in a 7 × 7 concentration matrix, with more than 4 million measurements of sensitivity, producing an exceptionally data-rich resource. We implement a new approach using combination Emax (viability effect) and highest single agent (HSA) to assess combination benefit. We designed a clinical translatability workflow to identify combinations with clearly defined patient populations, rationale for tolerability based on tumor type and combination-specific "emergent" biomarkers, and exposures relevant to clinical doses. We describe three actionable combinations in defined cancer types, confirmed in vitro and in vivo, with a focus on hematologic cancers and apoptotic targets. SIGNIFICANCE: We present the largest cancer drug combination screen published to date with 7 × 7 concentration response matrices for 109 combinations in more than 750 cell lines, complemented by multi-omics predictors of response and identification of "emergent" combination biomarkers. We prioritize hits to optimize clinical translatability, and experimentally validate novel combination hypotheses. This article is featured in Selected Articles from This Issue, p. 695.

Novel WRN Helicase Inhibitors Selectively Target Microsatellite-Unstable Cancer Cells.

Microsatellite-unstable (MSI) cancers require WRN helicase to resolve replication stress due to expanded DNA (TA)n dinucleotide repeats. WRN is a promising synthetic lethal target for MSI tumors, and WRN inhibitors are in development. In this study, we used CRISPR-Cas9 base editing to map WRN residues critical for MSI cells, validating the helicase domain as the primary drug target. Fragment-based screening led to the development of potent and highly selective WRN helicase covalent inhibitors. These compounds selectively suppressed MSI model growth in vitro and in vivo by mimicking WRN loss, inducing DNA double-strand breaks at expanded TA repeats and DNA damage. Assessment of biomarkers in preclinical models linked TA-repeat expansions and mismatch repair alterations to compound activity. Efficacy was confirmed in immunotherapy-resistant organoids and patient-derived xenograft models. The discovery of potent, selective covalent WRN inhibitors provides proof of concept for synthetic lethal targeting of WRN in MSI cancer and tools to dissect WRN biology. Significance: We report the discovery and characterization of potent, selective WRN helicase inhibitors for MSI cancer treatment, with biomarker analysis and evaluation of efficacy in vivo and in immunotherapy-refractory preclinical models. These findings pave the way to translate WRN inhibition into MSI cancer therapies and provide tools to investigate WRN biology. See related commentary by Wainberg, p. 1369.

Patterns of somatic mutation in human cancer genomes.

Cancers arise owing to mutations in a subset of genes that confer growth advantage. The availability of the human genome sequence led us to propose that systematic resequencing of cancer genomes for mutations would lead to the discovery of many additional cancer genes. Here we report more than 1,000 somatic mutations found in 274 megabases (Mb) of DNA corresponding to the coding exons of 518 protein kinase genes in 210 diverse human cancers. There was substantial variation in the number and pattern of mutations in individual cancers reflecting different exposures, DNA repair defects and cellular origins. Most somatic mutations are likely to be 'passengers' that do not contribute to oncogenesis. However, there was evidence for 'driver' mutations contributing to the development of the cancers studied in approximately 120 genes. Systematic sequencing of cancer genomes therefore reveals the evolutionary diversity of cancers and implicates a larger repertoire of cancer genes than previously anticipated.

Base editing screens map mutations affecting interferon-γ signaling in cancer.

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.

A suspension technique for efficient large-scale cancer organoid culturing and perturbation screens.

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.

Functional Genomic Identification of Predictors of Sensitivity and Mechanisms of Resistance to Multivalent Second-Generation TRAIL-R2 Agonists.

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.

Pan-cancer proteomic map of 949 human cell lines.

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.

A hypermutation phenotype and somatic MSH6 mutations in recurrent human malignant gliomas after alkylator chemotherapy.

Malignant gliomas have a very poor prognosis. The current standard of care for these cancers consists of extended adjuvant treatment with the alkylating agent temozolomide after surgical resection and radiotherapy. Although a statistically significant increase in survival has been reported with this regimen, nearly all gliomas recur and become insensitive to further treatment with this class of agents. We sequenced 500 kb of genomic DNA corresponding to the kinase domains of 518 protein kinases in each of nine gliomas. Large numbers of somatic mutations were observed in two gliomas recurrent after alkylating agent treatment. The pattern of mutations in these cases showed strong similarity to that induced by alkylating agents in experimental systems. Further investigation revealed inactivating somatic mutations of the mismatch repair gene MSH6 in each case. We propose that inactivating somatic mutations of MSH6 confer resistance to alkylating agents in gliomas in vivo and concurrently unleash accelerated mutagenesis in resistant clones as a consequence of continued exposure to alkylating agents in the presence of defective mismatch repair. The evidence therefore suggests that when MSH6 is inactivated in gliomas, alkylating agents convert from induction of tumor cell death to promotion of neoplastic progression. These observations highlight the potential of large scale sequencing for revealing and elucidating mutagenic processes operative in individual human cancers.

Loss of functional BAP1 augments sensitivity to TRAIL in cancer cells.

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.

The Origins and Vulnerabilities of Two Transmissible Cancers in Tasmanian Devils.

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.