Reproducible pharmacogenomic profiling of cancer cell line panels.

The use of large-scale genomic and drug response screening of cancer cell lines depends crucially on the reproducibility of results. Here we consider two previously published screens, plus a later critique of these studies. Using independent data, we show that consistency is achievable, and provide a systematic description of the best laboratory and analysis practices for future studies.

Germline genetic polymorphisms influence tumor gene expression and immune cell infiltration.

Cancer immunotherapy has emerged as an effective therapy in a variety of cancers. However, a key challenge in the field is that only a subset of patients who receive immunotherapy exhibit durable response. It has been hypothesized that host genetics influences the inherent immune profiles of patients and may underlie their differential response to immunotherapy. Herein, we systematically determined the association of common germline genetic variants with gene expression and immune cell infiltration of the tumor. We identified 64,094 expression quantitative trait loci (eQTLs) that associated with 18,210 genes (eGenes) across 24 human cancers. Overall, eGenes were enriched for their being involved in immune processes, suggesting that expression of immune genes can be shaped by hereditary genetic variants. We identified the endoplasmic reticulum aminopeptidase 2 (ERAP2) gene as a pan-cancer type eGene whose expression levels stratified overall survival in a subset of patients with bladder cancer receiving anti-PD-L1 (atezolizumab) therapy. Finally, we identified 103 gene signature QTLs (gsQTLs) that were associated with predicted immune cell abundance within the tumor microenvironment. Our findings highlight the impact of germline SNPs on cancer-immune phenotypes and response to therapy; and these analyses provide a resource for integration of germline genetics as a component of personalized cancer immunotherapy.

Recurrent R-spondin fusions in colon cancer.

Identifying and understanding changes in cancer genomes is essential for the development of targeted therapeutics. Here we analyse systematically more than 70 pairs of primary human colon tumours by applying next-generation sequencing to characterize their exomes, transcriptomes and copy-number alterations. We have identified 36,303 protein-altering somatic changes that include several new recurrent mutations in the Wnt pathway gene TCF7L2, chromatin-remodelling genes such as TET2 and TET3 and receptor tyrosine kinases including ERBB3. Our analysis for significantly mutated cancer genes identified 23 candidates, including the cell cycle checkpoint kinase ATM. Copy-number and RNA-seq data analysis identified amplifications and corresponding overexpression of IGF2 in a subset of colon tumours. Furthermore, using RNA-seq data we identified multiple fusion transcripts including recurrent gene fusions involving R-spondin family members RSPO2 and RSPO3 that together occur in 10% of colon tumours. The RSPO fusions were mutually exclusive with APC mutations, indicating that they probably have a role in the activation of Wnt signalling and tumorigenesis. Consistent with this we show that the RSPO fusion proteins were capable of potentiating Wnt signalling. The R-spondin gene fusions and several other gene mutations identified in this study provide new potential opportunities for therapeutic intervention in colon cancer.

Application of Mass Spectrometry Profiling to Establish Brusatol as an Inhibitor of Global Protein Synthesis.

The KEAP1/Nrf2 pathway senses and responds to changes in intracellular oxidative stress. Mutations that result in constitutive activation of Nrf2 are present in several human tumors, especially non-small cell lung cancer. Therefore, compounds that inhibit Nrf2 activity might be beneficial in treating patients whose tumors show activation of this pathway. Recent reports suggest that the natural product brusatol can potently and selectively inhibit Nrf2 activity, resulting in cell cytotoxicity, and can be effectively combined with chemotherapeutic agents. Here, we analyzed the effects of brusatol on the cellular proteome in the KEAP1 mutant non-small cell lung cancer cell line A549. Brusatol was found to rapidly and potently decrease the expression of the majority of detected proteins, including Nrf2. The most dramatically decreased proteins are those that display a short half-life, like Nrf2. This effect was confirmed by restricting the analysis to newly synthesized proteins using a labeled methionine analogue. Moreover, brusatol increased the expression of multiple components of the ribosome, suggesting that it regulates the function of this macromolecular complex. Finally, we show that brusatol induces its potent cellular cytotoxicity effects on multiple cancer cell lines in a manner independent of KEAP1/Nrf2 activity and with a profile similar to the protein translation inhibitor silvestrol. In conclusion, our data show that the activity of brusatol is not restricted to Nrf2 but, rather, functions as a global protein synthesis inhibitor.

The mutation spectrum revealed by paired genome sequences from a lung cancer patient.

Lung cancer is the leading cause of cancer-related mortality worldwide, with non-small-cell lung carcinomas in smokers being the predominant form of the disease. Although previous studies have identified important common somatic mutations in lung cancers, they have primarily focused on a limited set of genes and have thus provided a constrained view of the mutational spectrum. Recent cancer sequencing efforts have used next-generation sequencing technologies to provide a genome-wide view of mutations in leukaemia, breast cancer and cancer cell lines. Here we present the complete sequences of a primary lung tumour (60x coverage) and adjacent normal tissue (46x). Comparing the two genomes, we identify a wide variety of somatic variations, including >50,000 high-confidence single nucleotide variants. We validated 530 somatic single nucleotide variants in this tumour, including one in the KRAS proto-oncogene and 391 others in coding regions, as well as 43 large-scale structural variations. These constitute a large set of new somatic mutations and yield an estimated 17.7 per megabase genome-wide somatic mutation rate. Notably, we observe a distinct pattern of selection against mutations within expressed genes compared to non-expressed genes and in promoter regions up to 5 kilobases upstream of all protein-coding genes. Furthermore, we observe a higher rate of amino acid-changing mutations in kinase genes. We present a comprehensive view of somatic alterations in a single lung tumour, and provide the first evidence, to our knowledge, of distinct selective pressures present within the tumour environment.

Diverse somatic mutation patterns and pathway alterations in human cancers.

The systematic characterization of somatic mutations in cancer genomes is essential for understanding the disease and for developing targeted therapeutics. Here we report the identification of 2,576 somatic mutations across approximately 1,800 megabases of DNA representing 1,507 coding genes from 441 tumours comprising breast, lung, ovarian and prostate cancer types and subtypes. We found that mutation rates and the sets of mutated genes varied substantially across tumour types and subtypes. Statistical analysis identified 77 significantly mutated genes including protein kinases, G-protein-coupled receptors such as GRM8, BAI3, AGTRL1 (also called APLNR) and LPHN3, and other druggable targets. Integrated analysis of somatic mutations and copy number alterations identified another 35 significantly altered genes including GNAS, indicating an expanded role for galpha subunits in multiple cancer types. Furthermore, our experimental analyses demonstrate the functional roles of mutant GNAO1 (a Galpha subunit) and mutant MAP2K4 (a member of the JNK signalling pathway) in oncogenesis. Our study provides an overview of the mutational spectra across major human cancers and identifies several potential therapeutic targets.

Small molecule inhibition of group I p21-activated kinases in breast cancer induces apoptosis and potentiates the activity of microtubule stabilizing agents.

INTRODUCTION: Breast cancer, the most common cause of cancer-related deaths worldwide among women, is a molecularly and clinically heterogeneous disease. Extensive genetic and epigenetic profiling of breast tumors has recently revealed novel putative driver genes, including p21-activated kinase (PAK)1. PAK1 is a serine/threonine kinase downstream of small GTP-binding proteins, Rac1 and Cdc42, and is an integral component of growth factor signaling networks and cellular functions fundamental to tumorigenesis. METHODS: PAK1 dysregulation (copy number gain, mRNA and protein expression) was evaluated in two cohorts of breast cancer tissues (n=980 and 1,108). A novel small molecule inhibitor, FRAX1036, and RNA interference were used to examine PAK1 loss of function and combination with docetaxel in vitro. Mechanism of action for the therapeutic combination, both cellular and molecular, was assessed via time-lapse microscopy and immunoblotting. RESULTS: We demonstrate that focal genomic amplification and overexpression of PAK1 are associated with poor clinical outcome in the luminal subtype of breast cancer (P=1.29×10(-4) and P=0.015, respectively). Given the role for PAK1 in regulating cytoskeletal organization, we hypothesized that combination of PAK1 inhibition with taxane treatment could be combined to further interfere with microtubule dynamics and cell survival. Consistent with this, administration of docetaxel with either a novel small molecule inhibitor of group I PAKs, FRAX1036, or PAK1 small interfering RNA oligonucleotides dramatically altered signaling to cytoskeletal-associated proteins, such as stathmin, and induced microtubule disorganization and cellular apoptosis. Live-cell imaging revealed that the duration of mitotic arrest mediated by docetaxel was significantly reduced in the presence of FRAX1036, and this was associated with increased kinetics of apoptosis. CONCLUSIONS: Taken together, these findings further support PAK1 as a potential target in breast cancer and suggest combination with taxanes as a viable strategy to increase anti-tumor efficacy.

The effects of hepatitis B virus integration into the genomes of hepatocellular carcinoma patients.

Hepatitis B virus (HBV) infection is a leading risk factor for hepatocellular carcinoma (HCC). HBV integration into the host genome has been reported, but its scale, impact and contribution to HCC development is not clear. Here, we sequenced the tumor and nontumor genomes (>80× coverage) and transcriptomes of four HCC patients and identified 255 HBV integration sites. Increased sequencing to 240× coverage revealed a proportionally higher number of integration sites. Clonal expansion of HBV-integrated hepatocytes was found specifically in tumor samples. We observe a diverse collection of genomic perturbations near viral integration sites, including direct gene disruption, viral promoter-driven human transcription, viral-human transcript fusion, and DNA copy number alteration. Thus, we report the most comprehensive characterization of HBV integration in hepatocellular carcinoma patients. Such widespread random viral integration will likely increase carcinogenic opportunities in HBV-infected individuals.

Genome and transcriptome sequencing of lung cancers reveal diverse mutational and splicing events.

Lung cancer is a highly heterogeneous disease in terms of both underlying genetic lesions and response to therapeutic treatments. We performed deep whole-genome sequencing and transcriptome sequencing on 19 lung cancer cell lines and three lung tumor/normal pairs. Overall, our data show that cell line models exhibit similar mutation spectra to human tumor samples. Smoker and never-smoker cancer samples exhibit distinguishable patterns of mutations. A number of epigenetic regulators, including KDM6A, ASH1L, SMARCA4, and ATAD2, are frequently altered by mutations or copy number changes. A systematic survey of splice-site mutations identified 106 splice site mutations associated with cancer specific aberrant splicing, including mutations in several known cancer-related genes. RAC1b, an isoform of the RAC1 GTPase that includes one additional exon, was found to be preferentially up-regulated in lung cancer. We further show that its expression is significantly associated with sensitivity to a MAP2K (MEK) inhibitor PD-0325901. Taken together, these data present a comprehensive genomic landscape of a large number of lung cancer samples and further demonstrate that cancer-specific alternative splicing is a widespread phenomenon that has potential utility as therapeutic biomarkers. The detailed characterizations of the lung cancer cell lines also provide genomic context to the vast amount of experimental data gathered for these lines over the decades, and represent highly valuable resources for cancer biology.

PAK1 mediates pancreatic cancer cell migration and resistance to MET inhibition.

Pancreatic adenocarcinoma (PDAC) is a major unmet medical need and a deeper understanding of molecular drivers is needed to advance therapeutic options for patients. We report here that p21-activated kinase 1 (PAK1) is a central node in PDAC cells downstream of multiple growth factor signalling pathways, including hepatocyte growth factor (HGF) and MET receptor tyrosine kinase. PAK1 inhibition blocks signalling to cytoskeletal effectors and tumour cell motility driven by HGF/MET. MET antagonists, such as onartuzumab and crizotinib, are currently in clinical development. Given that even highly effective therapies have resistance mechanisms, we show that combination with PAK1 inhibition overcomes potential resistance mechanisms mediated either by activation of parallel growth factor pathways or by direct amplification of PAK1. Inhibition of PAK1 attenuated in vivo tumour growth and metastasis in a model of pancreatic adenocarcinoma. In human tissues, PAK1 is highly expressed in a proportion of PDACs (33% IHC score 2 or 3; n = 304) and its expression is significantly associated with MET positivity (p < 0.0001) and linked to a widespread metastatic pattern in patients (p = 0.067). Taken together, our results provide evidence for a functional role of MET/PAK1 signalling in pancreatic adenocarcinoma and support further characterization of therapeutic inhibitors in this indication.

Targeting p21-activated kinase 1 (PAK1) to induce apoptosis of tumor cells.

p21-activated kinases (PAKs) are serine/threonine protein kinases that serve as important mediators of Rac and Cdc42 GTPase function as well as pathways required for Ras-driven tumorigenesis. PAK1 has been implicated in signaling by growth factor receptors and morphogenetic processes that control cell polarity, invasion, and actin cytoskeleton organization. To better understand the role of PAK1 in tumorigenesis, PAK1 genomic copy number and expression were determined for a large panel of breast, lung, and head and neck tumors. PAK1 genomic amplification at 11q13 was prevalent in luminal breast cancer, and PAK1 protein expression was associated with lymph node metastasis. Breast cancer cells with PAK1 genomic amplification rapidly underwent apoptosis after inhibition of this kinase. Strong nuclear and cytoplasmic PAK1 expression was also prevalent in squamous nonsmall cell lung carcinomas (NSCLCs), and selective PAK1 inhibition was associated with delayed cell-cycle progression in vitro and in vivo. NSCLC cells were profiled using a library of pathway-targeted small-molecule inhibitors, and several synergistic combination therapies, including combination with antagonists of inhibitor of apoptosis proteins, were revealed for PAK1. Dual inhibition of PAK1 and X chromosome-linked inhibitor of apoptosis efficiently increased effector caspase activation and apoptosis of NSCLC cells. Together, our results provide evidence for dysregulation of PAK1 in breast and squamous NSCLCs and a role for PAK1 in cellular survival and proliferation in these indications.

Large-scale integrative network-based analysis identifies common pathways disrupted by copy number alterations across cancers.

BACKGROUND: Many large-scale studies analyzed high-throughput genomic data to identify altered pathways essential to the development and progression of specific types of cancer. However, no previous study has been extended to provide a comprehensive analysis of pathways disrupted by copy number alterations across different human cancers. Towards this goal, we propose a network-based method to integrate copy number alteration data with human protein-protein interaction networks and pathway databases to identify pathways that are commonly disrupted in many different types of cancer. RESULTS: We applied our approach to a data set of 2,172 cancer patients across 16 different types of cancers, and discovered a set of commonly disrupted pathways, which are likely essential for tumor formation in majority of the cancers. We also identified pathways that are only disrupted in specific cancer types, providing molecular markers for different human cancers. Analysis with independent microarray gene expression datasets confirms that the commonly disrupted pathways can be used to identify patient subgroups with significantly different survival outcomes. We also provide a network view of disrupted pathways to explain how copy number alterations affect pathways that regulate cell growth, cycle, and differentiation for tumorigenesis. CONCLUSIONS: In this work, we demonstrated that the network-based integrative analysis can help to identify pathways disrupted by copy number alterations across 16 types of human cancers, which are not readily identifiable by conventional overrepresentation-based and other pathway-based methods. All the results and source code are available at http://compbio.cs.umn.edu/NetPathID/.

Genetic alterations and oncogenic pathways associated with breast cancer subtypes.

Breast cancers can be divided into subtypes with important implications for prognosis and treatment. We set out to characterize the genetic alterations observed in different breast cancer subtypes and to identify specific candidate genes and pathways associated with subtype biology. mRNA expression levels of estrogen receptor, progesterone receptor, and HER2 were shown to predict marker status determined by immunohistochemistry and to be effective at assigning samples to subtypes. HER2(+) cancers were shown to have the greatest frequency of high-level amplification (independent of the ERBB2 amplicon itself), but triple-negative cancers had the highest overall frequencies of copy gain. Triple-negative cancers also were shown to have more frequent loss of phosphatase and tensin homologue and mutation of RB1, which may contribute to genomic instability. We identified and validated seven regions of copy number alteration associated with different subtypes, and used integrative bioinformatics analysis to identify candidate oncogenes and tumor suppressors, including ERBB2, GRB7, MYST2, PPM1D, CCND1, HDAC2, FOXA1, and RASA1. We tested the candidate oncogene MYST2 and showed that it enhances the anchorage-independent growth of breast cancer cells. The genome-wide and region-specific differences between subtypes suggest the differential activation of oncogenic pathways.

Super-enhancer acquisition drives oncogene expression in triple negative breast cancer.

Triple Negative Breast Cancer (TNBC) is a heterogeneous disease lacking known molecular drivers and effective targeted therapies. Cytotoxic chemotherapy remains the mainstay of treatment for TNBCs, which have significantly poorer survival rates compared to other breast cancer subtypes. In addition to changes within the coding genome, aberrant enhancer activity is a well-established contributor to tumorigenesis. Here we use H3K27Ac chromatin immunoprecipitation followed by sequencing (ChIP-Seq) to map the active cis-regulatory landscape in TNBC. We identify distinct disease subtypes associated with specific enhancer activity, and over 2,500 unique superenhancers acquired by tumor cells but absent from normal breast tissue. To identify potential actionable disease drivers, we probed the dependency on genes that associate with tumor-specific enhancers by CRISPR screening. In this way we identify a number of tumor-specific dependencies, including a previously uncharacterized dependency on the TGFß pseudo-receptor BAMBI.

High-resolution genomic and expression analyses of copy number alterations in breast tumors.

Analysis of recurrent DNA amplification can lead to the identification of cancer driver genes, but this process is often hampered by the low resolution of existing copy number analysis platforms. Fifty-one breast tumors were profiled for copy number alterations (CNAs) with the high-resolution Affymetrix 500K SNP array. These tumors were also expression-profiled and surveyed for mutations in selected genes commonly mutated in breast cancer (TP53, CDKN2A, ERBB2, KRAS, PIK3CA, PTEN). Combined analysis of common CNAs and mutations revealed putative associations between features. Analysis of both the prevalence and amplitude of CNAs defined regions of recurrent alteration. Compared with previous array comparative genomic hybridization studies, our analysis provided boundaries for frequently altered regions that were approximately one-fourth the size, greatly reducing the number of potential alteration-driving genes. Expression data from matched tumor samples were used to further interrogate the functional relevance of genes located in recurrent amplicons. Although our data support the importance of some known driver genes such as ERBB2, refined amplicon boundaries at other locations, such as 8p11-12 and 11q13.5-q14.2, greatly reduce the number of potential driver genes and indicate alternatives to commonly suggested driver genes in some cases. For example, the previously reported recurrent amplification at 17q23.2 is reduced to a 249 kb minimal region containing the putative driver RPS6KB1 as well as the putative oncogenic microRNA mir-21. High-resolution copy number analysis provides refined insight into many breast cancer amplicons and their relationships to gene expression, point mutations and breast cancer subtype classifications. This article contains Supplementary Material available at http://www.interscience.wiley.com/jpages/1045-2257/suppmat.

Distinguishing cancer-associated missense mutations from common polymorphisms.

Missense variants are commonly identified in genomic sequence but only a small fraction directly contribute to oncogenesis. The ability to distinguish those missense changes that contribute to cancer progression from those that do not is a difficult problem usually only accomplished through functional in vivo analyses. Using two computational algorithms, Sorting Intolerant from Tolerant (SIFT) and the Pfam-based LogR.E-value method, we have identified features that distinguish cancer-associated missense mutations from other classes of missense change. Our data reveal that cancer mutants behave similarly to Mendelian disease mutations, but are clearly distinct from either complex disease mutations or common single-nucleotide polymorphisms. We show that both activating and inactivating oncogenic mutations are predicted to be deleterious, although activating changes are likely to increase protein activity. Using the Gene Ontology and data from the SIFT and LogR.E-value metrics, a classifier was built that predicts cancer-associated missense mutations with a very low false-positive rate. The classifier does remarkably well in a number of different experiments designed to distinguish polymorphisms from true cancer-associated mutations. We also show that recurrently observed mutations are much more likely to be predicted to be cancer-associated than rare mutations, suggesting that our classifier will be useful in distinguishing causal from passenger mutations. In addition, from an expressed sequence tag-based screen, we identified a previously unknown germ line change (P1104A) in tumor tissues that is predicted to disrupt the function of the TYK2 protein. The data presented here show that this novel bioinformatics approach to classifying cancer-associated variants is robust and can be used for large-scale analyses.

Multiple-gene targeting and mismatch tolerance can confound analysis of genome-wide pooled CRISPR screens.

BACKGROUND: Genome-wide loss-of-function screens using the CRISPR/Cas9 system allow the efficient discovery of cancer cell vulnerabilities. While several studies have focused on correcting for DNA cleavage toxicity biases associated with copy number alterations, the effects of sgRNAs co-targeting multiple genomic loci in CRISPR screens have not been discussed. RESULTS: In this work, we analyze CRISPR essentiality screen data from 391 cancer cell lines to characterize biases induced by multi-target sgRNAs. We investigate two types of multi-targets: on-targets predicted through perfect sequence complementarity and off-targets predicted through sequence complementarity with up to two nucleotide mismatches. We find that the number of on-targets and off-targets both increase sgRNA activity in a cell line-specific manner and that existing additive models of gene knockout effects fail at capturing genetic interactions that may occur between co-targeted genes. We use synthetic lethality between paralog genes to show that genetic interactions can introduce biases in essentiality scores estimated from multi-target sgRNAs. We further show that single-mismatch tolerant sgRNAs can confound the analysis of gene essentiality and lead to incorrect co-essentiality functional networks. Lastly, we also find that single nucleotide polymorphisms located in protospacer regions can impair on-target activity as a result of mismatch tolerance. CONCLUSION: We show the impact of multi-target effects on estimating cancer cell dependencies and the impact of off-target effects caused by mismatch tolerance in sgRNA-DNA binding.

Transposable element expression in tumors is associated with immune infiltration and increased antigenicity.

Profound global loss of DNA methylation is a hallmark of many cancers. One potential consequence of this is the reactivation of transposable elements (TEs) which could stimulate the immune system via cell-intrinsic antiviral responses. Here, we develop REdiscoverTE, a computational method for quantifying genome-wide TE expression in RNA sequencing data. Using The Cancer Genome Atlas database, we observe increased expression of over 400 TE subfamilies, of which 262 appear to result from a proximal loss of DNA methylation. The most recurrent TEs are among the evolutionarily youngest in the genome, predominantly expressed from intergenic loci, and associated with antiviral or DNA damage responses. Treatment of glioblastoma cells with a demethylation agent results in both increased TE expression and de novo presentation of TE-derived peptides on MHC class I molecules. Therapeutic reactivation of tumor-specific TEs may synergize with immunotherapy by inducing inflammation and the display of potentially immunogenic neoantigens.

Enhancer Activity Requires CBP/P300 Bromodomain-Dependent Histone H3K27 Acetylation.

Acetylation of histone H3 at lysine 27 is a well-defined marker of enhancer activity. However, the functional impact of this modification at enhancers is poorly understood. Here, we use a chemical genetics approach to acutely block the function of the cAMP response element binding protein (CREB) binding protein (CBP)/P300 bromodomain in models of hematological malignancies and describe a consequent loss of H3K27Ac specifically from enhancers, despite the continued presence of CBP/P300 at chromatin. Using this approach to dissect the role of H3K27Ac at enhancers, we identify a critical role for this modification in the production of enhancer RNAs and transcription of enhancer-regulated gene networks.

A transcriptional MAPK Pathway Activity Score (MPAS) is a clinically relevant biomarker in multiple cancer types.

KRAS- and BRAF-mutant tumors are often dependent on MAPK signaling for proliferation and survival and thus sensitive to MAPK pathway inhibitors. However, clinical studies have shown that MEK inhibitors are not uniformly effective in these cancers indicating that mutational status of these oncogenes does not accurately capture MAPK pathway activity. A number of transcripts are regulated by this pathway and are recurrently identified in genome-based MAPK transcriptional signatures. To test whether the transcriptional output of only 10 of these targets could quantify MAPK pathway activity with potential predictive or prognostic clinical utility, we created a MAPK Pathway Activity Score (MPAS) derived from aggregated gene expression. In vitro, MPAS predicted sensitivity to MAPK inhibitors in multiple cell lines, comparable to or better than larger genome-based statistical models. Bridging in vitro studies and clinical samples, median MPAS from a given tumor type correlated with cobimetinib (MEK inhibitor) sensitivity of cancer cell lines originating from the same tissue type. Retrospective analyses of clinical datasets showed that MPAS was associated with the sensitivity of melanomas to vemurafenib (HR: 0.596) and negatively prognostic of overall or progression-free survival in both adjuvant and metastatic CRC (HR: 1.5 and 1.4), adrenal cancer (HR: 1.7), and HER2+ breast cancer (HR: 1.6). MPAS thus demonstrates potential clinical utility that warrants further exploration.