Assessing whether next-generation DNA sequencing (NGS) can be used to screen prostate cancer for multiple gene alterations in men routinely diagnosed with this disease and/or who are entered into clinical trials. Previous studies are limited and have reported only low success rates.

Cancer cells possess aberrant proteomes that can arise by the disruption of genes involved in physiological protein degradation. Here we demonstrate the presence of promoter CpG island hypermethylation-linked inactivation of DERL3 (Derlin-3), a key gene in the endoplasmic reticulum-associated protein degradation pathway, in human tumours. The restoration of in vitro and in vivo DERL3 activity highlights the tumour suppressor features of the gene. Using the stable isotopic labelling of amino acids in cell culture workflow for differential proteome analysis, we identify SLC2A1 (glucose transporter 1, GLUT1) as a downstream target of DERL3. Most importantly, SLC2A1 overexpression mediated by DERL3 epigenetic loss contributes to the Warburg effect in the studied cells and pinpoints a subset of human tumours with greater vulnerability to drugs targeting glycolysis.

Ionizing radiation is a potent carcinogen, inducing cancer through DNA damage. The signatures of mutations arising in human tissues following in vivo exposure to ionizing radiation have not been documented. Here, we searched for signatures of ionizing radiation in 12 radiation-associated second malignancies of different tumour types. Two signatures of somatic mutation characterize ionizing radiation exposure irrespective of tumour type. Compared with 319 radiation-naive tumours, radiation-associated tumours carry a median extra 201 deletions genome-wide, sized 1-100 base pairs often with microhomology at the junction. Unlike deletions of radiation-naive tumours, these show no variation in density across the genome or correlation with sequence context, replication timing or chromatin structure. Furthermore, we observe a significant increase in balanced inversions in radiation-associated tumours. Both small deletions and inversions generate driver mutations. Thus, ionizing radiation generates distinctive mutational signatures that explain its carcinogenic potential.

Predicting the response of a specific cancer to a therapy is a major goal in modern oncology that should ultimately lead to a personalised treatment. High-throughput screenings of potentially active compounds against a panel of genomically heterogeneous cancer cell lines have unveiled multiple relationships between genomic alterations and drug responses. Various computational approaches have been proposed to predict sensitivity based on genomic features, while others have used the chemical properties of the drugs to ascertain their effect. In an effort to integrate these complementary approaches, we developed machine learning models to predict the response of cancer cell lines to drug treatment, quantified through IC₅₀ values, based on both the genomic features of the cell lines and the chemical properties of the considered drugs. Models predicted IC₅₀ values in a 8-fold cross-validation and an independent blind test with coefficient of determination R² of 0.72 and 0.64 respectively. Furthermore, models were able to predict with comparable accuracy (R² of 0.61) IC50s of cell lines from a tissue not used in the training stage. Our in silico models can be used to optimise the experimental design of drug-cell screenings by estimating a large proportion of missing IC₅₀ values rather than experimentally measuring them. The implications of our results go beyond virtual drug screening design: potentially thousands of drugs could be probed in silico to systematically test their potential efficacy as anti-tumour agents based on their structure, thus providing a computational framework to identify new drug repositioning opportunities as well as ultimately be useful for personalized medicine by linking the genomic traits of patients to drug sensitivity.

The tumor suppressor BRCA2 plays a key role in genome integrity by promoting replication-fork stability and homologous recombination (HR) DNA repair. Here we report that human cancer cells lacking BRCA2 rely on the Fanconi anemia protein FANCD2 to limit replication-fork progression and genomic instability. Our results identify a new role of FANCD2 in limiting constitutive replication stress in BRCA2-deficient cells, thereby affecting cell survival and treatment responses.

Human cancer cell lines that can be propagated and manipulated in culture have proven to be excellent models for studying many aspects of gene function in cancer. In addition, they can provide a powerful system for assessing the molecular determinants of sensitivity to anticancer drugs. They have also been used in recent studies to identify genomic alterations and gene expression patterns that provide important insights into the genetic features that distinguish the properties of tumor cells associated with similar histologies. We have established a large repository of human tumor cell lines (>1000) corresponding to a wide variety of tumor types, and we have developed a methodology for profiling the collection for sensitivity to putative anticancer compounds. The rationale for examining tumor cell lines on this relatively large scale reflects accumulating evidence indicating that there is substantial genetic heterogeneity among human tumor cells-even those derived from tumors of similar histologies. Thus, to develop an accurate picture of the molecular determinants of tumorigenesis and response to therapy, it is essential to study the nature of such heterogeneity in a relatively large sample set. Here, we describe the methodologies used to conduct such screens and we describe a "proof-of-concept" screen using the EGFR kinase inhibitor, erlotinib (Tarceva), with a panel of lung cancer lines to demonstrate a correlation between EGFR mutations and drug sensitivity.

Improving early detection of colorectal cancer (CRC) is a key public health priority as adenomas and stage I cancer can be treated with minimally invasive procedures. Population screening strategies based on detection of occult blood in the feces have contributed to enhance detection rates of localized disease, but new approaches based on genetic analyses able to increase specificity and sensitivity could provide additional advantages compared to current screening methodologies. Recently, circulating cell-free DNA (cfDNA) has received much attention as a cancer biomarker for its ability to monitor the progression of advanced disease, predict tumor recurrence and reflect the complex genetic heterogeneity of cancers. Here, we tested whether analysis of cfDNA is a viable tool to enhance detection of colon adenomas. To address this, we assessed a cohort of patients with adenomas and healthy controls using droplet digital PCR (ddPCR) and mutation-specific assays targeted to trunk mutations. Additionally, we performed multiregional, targeted next-generation sequencing (NGS) of adenomas and unmasked extensive heterogeneity, affecting known drivers such as APC, KRAS and mismatch repair (MMR) genes. However, tumor-related mutations were undetectable in patients' plasma. Finally, we employed a preclinical mouse model of Apc-driven intestinal adenomas and confirmed the inability to identify tumor-related alterations via cfDNA, despite the enhanced disease burden displayed by this experimental cancer model. Therefore, we conclude that benign colon lesions display extensive genetic heterogeneity, that they are not prone to release DNA into the circulation and are unlikely to be reliably detected with liquid biopsies, at least with the current technologies.

Cancer stem cells (CSCs) correlate with recurrence, metastasis and poor survival in clinical studies. Encouraging results from clinical trials of CSC inhibitors have further validated CSCs as therapeutic targets. ONC201 is a first-in-class small molecule imipridone in Phase I/II clinical trials for advanced cancer. We have previously shown that ONC201 targets self-renewing, chemotherapy-resistant colorectal CSCs via Akt/ERK inhibition and DR5/TRAIL induction. In this study, we demonstrate that the anti-CSC effects of ONC201 involve early changes in stem cell-related gene expression prior to tumor cell death induction. A targeted network analysis of gene expression profiles in colorectal cancer cells revealed that ONC201 downregulates stem cell pathways such as Wnt signaling and modulates genes (ID1, ID2, ID3 and ALDH7A1) known to regulate self-renewal in colorectal, prostate cancer and glioblastoma. ONC201-mediated changes in CSC-related gene expression were validated at the RNA and protein level for each tumor type. Accordingly, we observed inhibition of self-renewal and CSC markers in prostate cancer cell lines and patient-derived glioblastoma cells upon ONC201 treatment. Interestingly, ONC201-mediated CSC depletion does not occur in colorectal cancer cells with acquired resistance to ONC201. Finally, we observed that basal expression of CSC-related genes (ID1, CD44, HES7 and TCF3) significantly correlate with ONC201 efficacy in >1000 cancer cell lines and combining the expression of multiple genes leads to a stronger overall prediction. These proof-of-concept studies provide a rationale for testing CSC expression at the RNA and protein level as a predictive and pharmacodynamic biomarker of ONC201 response in ongoing clinical studies.

Osteosarcoma is a primary malignancy of bone that affects children and adults. Here, we present the largest sequencing study of osteosarcoma to date, comprising 112 childhood and adult tumours encompassing all major histological subtypes. A key finding of our study is the identification of mutations in insulin-like growth factor (IGF) signalling genes in 8/112 (7%) of cases. We validate this observation using fluorescence in situ hybridization (FISH) in an additional 87 osteosarcomas, with IGF1 receptor (IGF1R) amplification observed in 14% of tumours. These findings may inform patient selection in future trials of IGF1R inhibitors in osteosarcoma. Analysing patterns of mutation, we identify distinct rearrangement profiles including a process characterized by chromothripsis and amplification. This process operates recurrently at discrete genomic regions and generates driver mutations. It may represent an age-independent mutational mechanism that contributes to the development of osteosarcoma in children and adults alike.

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.

Single-cell technologies have described heterogeneity across tissues, but the spatial distribution and forces that drive single-cell phenotypes have not been well defined. Combining single-cell RNA and protein analytics in studying the role of stromal cancer-associated fibroblasts (CAFs) in modulating heterogeneity in pancreatic cancer (pancreatic ductal adenocarcinoma [PDAC]) model systems, we have identified significant single-cell population shifts toward invasive epithelial-to-mesenchymal transition (EMT) and proliferative (PRO) phenotypes linked with mitogen-activated protein kinase (MAPK) and signal transducer and activator of transcription 3 (STAT3) signaling. Using high-content digital imaging of RNA in situ hybridization in 195 PDAC tumors, we quantified these EMT and PRO subpopulations in 319,626 individual cancer cells that can be classified within the context of distinct tumor gland "units." Tumor gland typing provided an additional layer of intratumoral heterogeneity that was associated with differences in stromal abundance and clinical outcomes. This demonstrates the impact of the stroma in shaping tumor architecture by altering inherent patterns of tumor glands in human PDAC.

Evolved resistance to tyrosine kinase inhibitor (TKI)-targeted therapies remains a major clinical challenge. In epidermal growth factor receptor (EGFR) mutant non-small-cell lung cancer (NSCLC), failure of EGFR TKIs can result from both genetic and epigenetic mechanisms of acquired drug resistance. Widespread reports of histologic and gene expression changes consistent with an epithelial-to-mesenchymal transition (EMT) have been associated with initially surviving drug-tolerant persister cells, which can seed bona fide genetic mechanisms of resistance to EGFR TKIs. While therapeutic approaches targeting fully resistant cells, such as those harboring an EGFRT790M mutation, have been developed, a clinical strategy for preventing the emergence of persister cells remains elusive. Using mesenchymal cell lines derived from biopsies of patients who progressed on EGFR TKI as surrogates for persister populations, we performed whole-genome CRISPR screening and identified fibroblast growth factor receptor 1 (FGFR1) as the top target promoting survival of mesenchymal EGFR mutant cancers. Although numerous previous reports of FGFR signaling contributing to EGFR TKI resistance in vitro exist, the data have not yet been sufficiently compelling to instigate a clinical trial testing this hypothesis, nor has the role of FGFR in promoting the survival of persister cells been elucidated. In this study, we find that combining EGFR and FGFR inhibitors inhibited the survival and expansion of EGFR mutant drug-tolerant cells over long time periods, preventing the development of fully resistant cancers in multiple vitro models and in vivo. These results suggest that dual EGFR and FGFR blockade may be a promising clinical strategy for both preventing and overcoming EMT-associated acquired drug resistance and provide motivation for the clinical study of combined EGFR and FGFR inhibition in EGFR-mutated NSCLCs.

Metastatic colorectal cancer is a leading cause of cancer death. However, current therapy options are limited to chemotherapy, with the addition of anti-EGFR antibodies for patients with RAS wild-type tumours. Novel drug targets, or drug combinations that induce a synergistic response, would be of great benefit to patients. The identification of genes that are essential for cell survival can be undertaken using functional genomics screens. Furthermore, performing such screens in the presence of a targeted agent would allow the identification of combinations that result in a synthetic lethal interaction. Here, we present a dataset containing the results of a large scale RNAi screen (815 genes) to detect essential genes as well as synergistic combinations with targeted therapeutic agents using a panel of 27 colorectal cancer cell lines. These data identify genes that are essential for colorectal cancer cell survival as well as synthetic lethal treatment combinations using novel computational approaches. Moreover, this dataset could be utilised in combination with genomic profiling to identify predictive biomarkers of response.

Patients with seemingly the same tumour can respond very differently to treatment. There are strong, well-established effects of somatic mutations on drug efficacy, but there is at-most anecdotal evidence of a germline component to drug response. Here, we report a systematic survey of how inherited germline variants affect drug susceptibility in cancer cell lines. We develop a joint analysis approach that leverages both germline and somatic variants, before applying it to screening data from 993 cell lines and 265 drugs. Surprisingly, we find that the germline contribution to variation in drug susceptibility can be as large or larger than effects due to somatic mutations. Several of the associations identified have a direct relationship to the drug target. Finally, using 17-AAG response as an example, we show how germline effects in combination with transcriptomic data can be leveraged for improved patient stratification and to identify new markers for drug sensitivity.

Many gene fusions are reported in tumours and for most their role remains unknown. As fusions are used for diagnostic and prognostic purposes, and are targets for treatment, it is crucial to assess their function in cancer. To systematically investigate the role of fusions in tumour cell fitness, we utilized RNA-sequencing data from 1011 human cancer cell lines to functionally link 8354 fusion events with genomic data, sensitivity to >350 anti-cancer drugs and CRISPR-Cas9 loss-of-fitness effects. Established clinically-relevant fusions were identified. Overall, detection of functional fusions was rare, including those involving cancer driver genes, suggesting that many fusions are dispensable for tumour fitness. Therapeutically actionable fusions involving RAF1, BRD4 and ROS1 were verified in new histologies. In addition, recurrent YAP1-MAML2 fusions were identified as activators of Hippo-pathway signaling in multiple cancer types. Our approach discriminates functional fusions, identifying new drivers of carcinogenesis and fusions that could have clinical implications.

Mitochondria are typically essential for the viability of eukaryotic cells, and utilize oxygen and nutrients (e.g. glucose) to perform key metabolic functions that maintain energetic homeostasis and support proliferation. Here we provide a comprehensive functional annotation of mitochondrial genes that are essential for the viability of a large panel (625) of tumour cell lines. We perform genome-wide CRISPR/Cas9 deletion screening in normoxia-glucose, hypoxia-glucose and normoxia-galactose conditions, and identify both unique and overlapping genes whose loss influences tumour cell viability under these different metabolic conditions. We discover that loss of certain oxidative phosphorylation (OXPHOS) genes (e.g. SDHC) improves tumour cell growth in hypoxia-glucose, but reduces growth in normoxia, indicating a metabolic switch in OXPHOS gene function. Moreover, compared to normoxia-glucose, loss of genes involved in energy-consuming processes that are energetically demanding, such as translation and actin polymerization, improve cell viability under both hypoxia-glucose and normoxia-galactose. Collectively, our study defines mitochondrial gene essentiality in tumour cells, highlighting that essentiality is dependent on the metabolic environment, and identifies routes for regulating tumour cell viability in hypoxia.

Cancer-associated fibroblasts (CAFs) are highly heterogeneous. With the lack of a comprehensive understanding of CAFs' functional distinctions, it remains unclear how cancer treatments could be personalized based on CAFs in a patient's tumor. We have established a living biobank of CAFs derived from biopsies of patients' non-small lung cancer (NSCLC) that encompasses a broad molecular spectrum of CAFs in clinical NSCLC. By functionally interrogating CAF heterogeneity using the same therapeutics received by patients, we identify three functional subtypes: (1) robustly protective of cancers and highly expressing HGF and FGF7; (2) moderately protective of cancers and highly expressing FGF7; and (3) those providing minimal protection. These functional differences among CAFs are governed by their intrinsic TGF-β signaling, which suppresses HGF and FGF7 expression. This CAF functional classification correlates with patients' clinical response to targeted therapies and also associates with the tumor immune microenvironment, therefore providing an avenue to guide personalized treatment.

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.

Osimertinib is an EGFR tyrosine kinase inhibitor (TKI) with proven clinical efficacy; however, acquired resistance presents an obstacle to curing EGFR-driven disease. Recent studies have shown that drug-tolerant persister cells (DTP) have a distinct transcriptional profile that may confer specific vulnerabilities. By definition these cells avoid apoptosis, yet little is known about how their survival is regulated. We found that paradoxically, the proapoptotic gene BIM was upregulated in osimertinib DTPs, and cotreatment with BH3 mimetics could trigger DTP cell death. Furthermore, cIAP proteins, antiapoptotic members of the extrinsic pathway, were significantly elevated in DTPs. cIAP antagonists could block DTP formation as an up-front combination, and could eliminate preformed DTPs. Critically, when treated at the time of maximal osimertinib response, cIAP or MCL1 inhibitor treatment could significantly attenuate the regrowth of EGFRm cell line mouse xenografts. Finally, we show that apoptosis can be maximized in cell lines with acquired osimertinib resistance by combining BH3 or SMAC mimetics with agents that target the resistance driver in these models. Taken together, these data suggest novel therapeutic strategies at the point of minimal residual disease or full osimertinib resistance for patients in this critical area of unmet need.

Recent clinical data have linked KRAS/TP53 comutation (mut) to resistance to radiotherapy (RT), but supporting laboratory in vivo evidence is lacking. In addition, the ability of different radiation doses, with/without epidermal growth factor receptor (EGFR)-directed treatment, to achieve local tumor control as a function of KRAS status is unknown. Here, we assessed clonogenic radiation survival of a panel of annotated lung cancer cell lines. KRASmut/TP53mut was associated with the highest radioresistance in nonisogenic and isogenic comparisons. To validate these findings, isogenic TP53mut NCI-H1703 models, KRASmut or wild-type (wt), were grown as heterotopic xenografts in nude mice. A clinical RT schedule of 30 fractions over 6 weeks was employed. The dose that controlled 50% of tumors (TCD50 ) was calculated. The TCD50 for KRASwt/TP53mut xenografts was 43.1 Gy whereas KRASmut/TP53mut tumors required a 1.9-fold higher TCD50 of 81.4 Gy. The EGFR inhibitor erlotinib radiosensitized KRASmut but not KRASwt cells and xenografts. The TCD50 associated with adding erlotinib to RT was 58.8 Gy for KRASmut, that is, a ~1.4-fold dose enhancement. However, the EGFR antibody cetuximab did not have a radiosensitizing effect. In conclusion, we demonstrate for the first time that KRASmut in a TP53mut background confers radioresistance when studying a clinical RT schedule and local control rather than tumor growth delay. Despite the known unresponsiveness of KRASmut tumors to EGFR inhibitors, erlotinib radiosensitized KRASmut tumors. Our data highlight KRAS/TP53 comutation as a candidate biomarker of radioresistance that can be at least partially reversed by dose escalation or the addition of a targeted agent.