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.

Drug mechanism-of-action discovery through the integration of pharmacological and CRISPR screens.

Low success rates during drug development are due, in part, to the difficulty of defining drug mechanism-of-action and molecular markers of therapeutic activity. Here, we integrated 199,219 drug sensitivity measurements for 397 unique anti-cancer drugs with genome-wide CRISPR loss-of-function screens in 484 cell lines to systematically investigate cellular drug mechanism-of-action. We observed an enrichment for positive associations between the profile of drug sensitivity and knockout of a drug's nominal target, and by leveraging protein-protein networks, we identified pathways underpinning drug sensitivity. This revealed an unappreciated positive association between mitochondrial E3 ubiquitin-protein ligase MARCH5 dependency and sensitivity to MCL1 inhibitors in breast cancer cell lines. We also estimated drug on-target and off-target activity, informing on specificity, potency and toxicity. Linking drug and gene dependency together with genomic data sets uncovered contexts in which molecular networks when perturbed mediate cancer cell loss-of-fitness and thereby provide independent and orthogonal evidence of biomarkers for drug development. This study illustrates how integrating cell line drug sensitivity with CRISPR loss-of-function screens can elucidate mechanism-of-action to advance drug development.

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.

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.

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.

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.

Complete genome sequence of the erythromycin-producing bacterium Saccharopolyspora erythraea NRRL23338.

Saccharopolyspora erythraea is used for the industrial-scale production of the antibiotic erythromycin A, derivatives of which play a vital role in medicine. The sequenced chromosome of this soil bacterium comprises 8,212,805 base pairs, predicted to encode 7,264 genes. It is circular, like those of the pathogenic actinomycetes Mycobacterium tuberculosis and Corynebacterium diphtheriae, but unlike the linear chromosomes of the model actinomycete Streptomyces coelicolor A3(2) and the closely related Streptomyces avermitilis. The S. erythraea genome contains at least 25 gene clusters for production of known or predicted secondary metabolites, at least 72 genes predicted to confer resistance to a range of common antibiotic classes and many sets of duplicated genes to support its saprophytic lifestyle. The availability of the genome sequence of S. erythraea will improve insight into its biology and facilitate rational development of strains to generate high-titer producers of clinically important antibiotics.

Insights into polyether biosynthesis from analysis of the nigericin biosynthetic gene cluster in Streptomyces sp. DSM4137.

Nigericin was among the first polyether ionophores to be discovered, but its biosynthesis remains obscure. The biosynthetic gene cluster for nigericin has been serendipitously cloned from Streptomyces sp. DSM4137, and deletion of this gene cluster abolished the production of both nigericin and the closely related metabolite abierixin. Detailed comparison of the nigericin biosynthetic genes with their counterparts in the biosynthetic clusters for other polyketides has prompted a significant revision of the proposed common pathway for polyether biosynthesis. In particular, we present evidence that in nigericin, nanchangmycin, and monensin, an unusual ketosynthase-like protein, KSX, transfers the initially formed linear polyketide chain to a discrete acyl carrier protein, ACPX, for oxidative cyclization. Consistent with this, deletion of either monACPX or monKSX from the monensin gene cluster effectively abolished monensin A biosynthesis.

The gene cluster for fluorometabolite biosynthesis in Streptomyces cattleya: a thioesterase confers resistance to fluoroacetyl-coenzyme A.

A genomic library of Streptomyces cattleya was screened to isolate a gene cluster encoding enzymes responsible for the production of fluorine-containing metabolites. In addition to the previously described fluorinase FlA which catalyzes the formation of 5'-fluoro-5'-deoxyadenosine from S-adenosylmethionine and fluoride, 11 other putative open reading frames have been identified. Three of the proteins encoded by these genes have been characterized. FlB was determined to be the second enzyme in the pathway, catalyzing the phosphorolytic cleavage of 5'-fluoro-5'-deoxyadenosine to produce 5-fluoro-5-deoxy-D-ribose-1-phosphate. The enzyme FlI was found to be an S-adenosylhomocysteine hydrolase, which may act to relieve S-adenosylhomocysteine inhibition of the fluorinase. Finally, flK encodes a thioesterase which catalyzes the selective breakdown of fluoroacetyl-CoA but not acetyl-CoA, suggesting that it provides the producing strain with a mechanism for resistance to fluoroacetate.

Mutation analysis of 24 known cancer genes in the NCI-60 cell line set.

The panel of 60 human cancer cell lines (the NCI-60) assembled by the National Cancer Institute for anticancer drug discovery is a widely used resource. The NCI-60 has been characterized pharmacologically and at the molecular level more extensively than any other set of cell lines. However, no systematic mutation analysis of genes causally implicated in oncogenesis has been reported. This study reports the sequence analysis of 24 known cancer genes in the NCI-60 and an assessment of 4 of the 24 genes for homozygous deletions. One hundred thirty-seven oncogenic mutations were identified in 14 (APC, BRAF, CDKN2, CTNNB1, HRAS, KRAS, NRAS, SMAD4, PIK3CA, PTEN, RB1, STK11, TP53, and VHL) of the 24 genes. All lines have at least one mutation among the cancer genes examined, with most lines (73%) having more than one. Identification of those cancer genes mutated in the NCI-60, in combination with pharmacologic and molecular profiles of the cells, will allow for more informed interpretation of anticancer agent screening and will enhance the use of the NCI-60 cell lines for molecularly targeted screens.

Biosynthetic gene cluster of the glycopeptide antibiotic teicoplanin: characterization of two glycosyltransferases and the key acyltransferase.

The gene cluster encoding biosynthesis of the clinically important glycopeptide antibiotic teicoplanin has been cloned from Actinoplanes teichomyceticus. Forty-nine putative open reading frames (ORFs) were identified within an 89 kbp genetic locus and assigned roles in teicoplanin biosynthesis, export, resistance, and regulation. Two ORFs, designated orfs 1 and 10*, showed significant homology to known glycosyltransferases. When heterologously expressed in Escherichia coli, these glycosyltransferases were shown to catalyze the transfer of UDP-(N-acetyl)-glucosamine onto, respectively, 3-chloro-beta-hydroxytyrosine-6 (3-Cl-6betaHty) and 4-hydroxyphenylglycine-4 (4Hpg) of the teicoplanin heptapeptide aglycone. The product of another ORF, orf11*, was demonstrated in vitro to transfer n-acetyl-, n-butyryl-, and n-octanoyl-groups from acyl-CoA donors either to a free UDP-aminosugar or to an aminosugar moiety in the teicoplanin pseudoaglycone, thus identifying Orf11* as the key acyltransferase in teicoplanin maturation. These findings should accelerate the combinatorial engineering of new and improved glycopeptide drugs.

Combinations of PARP Inhibitors with Temozolomide Drive PARP1 Trapping and Apoptosis in Ewing's Sarcoma.

Ewing's sarcoma is a malignant pediatric bone tumor with a poor prognosis for patients with metastatic or recurrent disease. Ewing's sarcoma cells are acutely hypersensitive to poly (ADP-ribose) polymerase (PARP) inhibition and this is being evaluated in clinical trials, although the mechanism of hypersensitivity has not been directly addressed. PARP inhibitors have efficacy in tumors with BRCA1/2 mutations, which confer deficiency in DNA double-strand break (DSB) repair by homologous recombination (HR). This drives dependence on PARP1/2 due to their function in DNA single-strand break (SSB) repair. PARP inhibitors are also cytotoxic through inhibiting PARP1/2 auto-PARylation, blocking PARP1/2 release from substrate DNA. Here, we show that PARP inhibitor sensitivity in Ewing's sarcoma cells is not through an apparent defect in DNA repair by HR, but through hypersensitivity to trapped PARP1-DNA complexes. This drives accumulation of DNA damage during replication, ultimately leading to apoptosis. We also show that the activity of PARP inhibitors is potentiated by temozolomide in Ewing's sarcoma cells and is associated with enhanced trapping of PARP1-DNA complexes. Furthermore, through mining of large-scale drug sensitivity datasets, we identify a subset of glioma, neuroblastoma and melanoma cell lines as hypersensitive to the combination of temozolomide and PARP inhibition, potentially identifying new avenues for therapeutic intervention. These data provide insights into the anti-cancer activity of PARP inhibitors with implications for the design of treatment for Ewing's sarcoma patients with PARP inhibitors.

The putative elaiophylin biosynthetic gene cluster in Streptomyces sp. DSM4137 is adjacent to genes encoding adenosylcobalamin-dependent methylmalonyl CoA mutase and to genes for synthesis of cobalamin.

A type I PKS gene probe obtained from RAPB of the rapamycin producer Streptomyces hygroscopicus, strongly hybridised to 92 out of 1120 cosmids from a genomic library of the elaiophylin-producing strain Streptomyces sp. DSM4137. Partial cosmid sequencing suggested the presence of 10 separate sequences encoding type I PKS genes. One entire DNA sequence was obtained and found exactly to match the gene organisation expected for the biosynthesis of the unusual macrodiolide polyketide elaiophylin. The putative elaiophylin gene cluster contains five large open-reading frames encoding typical modular polyketide synthases, which together catalyse the synthesis of the octaketide monomer of elaiophylin. Other genes were identified that would be required for provision of the ethylmalonate extender unit, for the synthesis and attachment of 2-deoxy-L-fucose and in regulation, or in export of the product. Immediately adjacent to the putative elaiophylin biosynthetic gene cluster is a 30-kbp region containing the gene for adenosylcobalamin-dependent methylmalonyl CoA mutase and also genes involved in the biosynthesis of the cobalamin cofactor. Analysis of the latter gene set confirms the view that cbiD of the anaerobic pathway and cobF in the aerobic pathway catalyse the same methylation of precorrin-5. The proximity of these genes to the putative elaiophylin gene cluster can best be rationalised if in this organism succinyl-CoA is a significant source of the methylmalonate units for complex polyketide biosynthesis.

Somatic mutations of the histone H3K27 demethylase gene UTX in human cancer.

Somatically acquired epigenetic changes are present in many cancers. Epigenetic regulation is maintained via post-translational modifications of core histones. Here, we describe inactivating somatic mutations in the histone lysine demethylase gene UTX, pointing to histone H3 lysine methylation deregulation in multiple tumor types. UTX reintroduction into cancer cells with inactivating UTX mutations resulted in slowing of proliferation and marked transcriptional changes. These data identify UTX as a new human cancer gene.

A systematic, large-scale resequencing screen of X-chromosome coding exons in mental retardation.

Large-scale systematic resequencing has been proposed as the key future strategy for the discovery of rare, disease-causing sequence variants across the spectrum of human complex disease. We have sequenced the coding exons of the X chromosome in 208 families with X-linked mental retardation (XLMR), the largest direct screen for constitutional disease-causing mutations thus far reported. The screen has discovered nine genes implicated in XLMR, including SYP, ZNF711 and CASK reported here, confirming the power of this strategy. The study has, however, also highlighted issues confronting whole-genome sequencing screens, including the observation that loss of function of 1% or more of X-chromosome genes is compatible with apparently normal existence.

X-linked protocadherin 19 mutations cause female-limited epilepsy and cognitive impairment.

Epilepsy and mental retardation limited to females (EFMR) is a disorder with an X-linked mode of inheritance and an unusual expression pattern. Disorders arising from mutations on the X chromosome are typically characterized by affected males and unaffected carrier females. In contrast, EFMR spares transmitting males and affects only carrier females. Aided by systematic resequencing of 737 X chromosome genes, we identified different protocadherin 19 (PCDH19) gene mutations in seven families with EFMR. Five mutations resulted in the introduction of a premature termination codon. Study of two of these demonstrated nonsense-mediated decay of PCDH19 mRNA. The two missense mutations were predicted to affect adhesiveness of PCDH19 through impaired calcium binding. PCDH19 is expressed in developing brains of human and mouse and is the first member of the cadherin superfamily to be directly implicated in epilepsy or mental retardation.

Mutations in CUL4B, which encodes a ubiquitin E3 ligase subunit, cause an X-linked mental retardation syndrome associated with aggressive outbursts, seizures, relative macrocephaly, central obesity, hypogonadism, pes cavus, and tremor.

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.