Search

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., Competing Interests: I have read the journal's policy and the authors of this manuscript have the following competing interests: Ryan Raisner, Peter M Haverty, & Karen E Gascoigne are employees of Genentech and shareholders of Roche. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

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.

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.

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., Competing Interests: Conflict of interest statement: All authors are employees of Genentech, Inc., (Copyright © 2018 the Author(s). Published by PNAS.)

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., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: All authors are employees and shareholders of Roche/Genentech.

[This corrects the article DOI: 10.1371/journal.pone.0185862.].

Resistance invariably develops to antiandrogen therapies used to treat newly diagnosed prostate cancers, but effective treatments for castration-resistant disease remain elusive. Here, we report that the transcriptional coactivator CBP/p300 is required to maintain the growth of castration-resistant prostate cancer. To exploit this vulnerability, we developed a novel small-molecule inhibitor of the CBP/p300 bromodomain that blocks prostate cancer growth in vitro and in vivo Molecular dissection of the consequences of drug treatment revealed a critical role for CBP/p300 in histone acetylation required for the transcriptional activity of the androgen receptor and its target gene expression. Our findings offer a preclinical proof of concept for small-molecule therapies to target the CBP/p300 bromodomain as a strategy to treat castration-resistant prostate cancer. Cancer Res; 77(20); 5564-75. ©2017 AACR ., (©2017 American Association for Cancer Research.)

Mitogen-activated protein kinase (MAPK) pathway dysregulation is implicated in >30% of all cancers, rationalizing the development of RAF, MEK and ERK inhibitors. While BRAF and MEK inhibitors improve BRAF mutant melanoma patient outcomes, these inhibitors had limited success in other MAPK dysregulated tumors, with insufficient pathway suppression and likely pathway reactivation. In this study we show that inhibition of either MEK or ERK alone only transiently inhibits the MAPK pathway due to feedback reactivation. Simultaneous targeting of both MEK and ERK nodes results in deeper and more durable suppression of MAPK signaling that is not achievable with any dose of single agent, in tumors where feedback reactivation occurs. Strikingly, combined MEK and ERK inhibition is synergistic in RAS mutant models but only additive in BRAF mutant models where the RAF complex is dissociated from RAS and thus feedback productivity is disabled. We discovered that pathway reactivation in RAS mutant models occurs at the level of CRAF with combination treatment resulting in a markedly more active pool of CRAF. However, distinct from single node targeting, combining MEK and ERK inhibitor treatment effectively blocks the downstream signaling as assessed by transcriptional signatures and phospho-p90RSK. Importantly, these data reveal that MAPK pathway inhibitors whose activity is attenuated due to feedback reactivation can be rescued with sufficient inhibition by using a combination of MEK and ERK inhibitors. The MEK and ERK combination significantly suppresses MAPK pathway output and tumor growth in vivo to a greater extent than the maximum tolerated doses of single agents, and results in improved anti-tumor activity in multiple xenografts as well as in two Kras mutant genetically engineered mouse (GEM) models. Collectively, these data demonstrate that combined MEK and ERK inhibition is functionally unique, yielding greater than additive anti-tumor effects and elucidates a highly effective combination strategy in MAPK-dependent cancer, such as KRAS mutant tumors.

The role of essential amino acids in metabolic reprogramming of cancer cells is now well established, whereas the role of non-essential amino acids (NEAAs) in malignancy remains less clear. Here, we have identified an important role for the NEAA proline in the tumorigenic potential of a subset of cancer cells. By profiling a large panel of cancer cell lines, we observed that proline consumption and expression of proline biosynthesis enzymes were well correlated with clonogenic and tumorigenic potential. Moreover, proline starvation or inhibition of proline biosynthesis enzymes impaired clonogenic/tumorigenic potential. Cancer cells exhibiting dependency on exogenous proline displayed hyperactivation of the mTORC1-mediated 4EBP1 signaling axis, as well as unresolved ER stress. Exogenous proline alleviated ER stress and promoted cellular homeostasis and clonogenicity. Increased dependence on proline may therefore define a specific vulnerability in some cancers that can be exploited by proline depletion., (Copyright © 2016 Elsevier Inc. All rights reserved.)

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.

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., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

Purpose: Chemotherapies are limited by a narrow therapeutic index resulting in suboptimal exposure of the tumor to the drug and acquired tumor resistance. One approach to overcome this is through antibody-drug conjugates (ADC) that facilitate greater potency via target-specific delivery of highly potent cytotoxic agents., Experimental Design: In this study, we used a bioinformatics approach to identify the lymphocyte antigen 6 complex locus E (LY6E), an IFN-inducible glycosylphosphatidylinositol (GPI)-linked cell membrane protein as a promising ADC target. We developed a monoclonal anti-LY6E antibody and characterized in situ LY6E expression in over 750 cancer specimens and normal tissues. Target-dependent anti-LY6E ADC killing was investigated both in vitro and in vivo using patient-derived xenograft models., Results: Using in silico approaches, we found that LY6E was significantly overexpressed and amplified in a wide array of different human solid tumors. IHC analysis revealed high LY6E protein expression in a number of tumor types, such as breast, lung, gastric, ovarian, pancreatic, kidney and head/neck carcinomas. Characterization of the endocytic pathways for LY6E revealed that the LY6E-specific antibody is internalized into cells leading to lysosomal accumulation. Consistent with this, a LY6E-specific ADC inhibited in vitro cell proliferation and produced durable tumor regression in vivo in clinically relevant LY6E-expressing xenograft models., Conclusions: Our results identify LY6E as a highly promising molecular ADC target for a variety of solid tumor types with current unmet medical need., (©2015 American Association for Cancer Research.)

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.

Smoothened (SMO) inhibitors are under clinical investigation for the treatment of several cancers. Vismodegib is approved for the treatment of locally advanced and metastatic basal cell carcinoma (BCC). Most BCC patients experience significant clinical benefit on vismodegib, but some develop resistance. Genomic analysis of tumor biopsies revealed that vismodegib resistance is associated with Hedgehog (Hh) pathway reactivation, predominantly through mutation of the drug target SMO and to a lesser extent through concurrent copy number changes in SUFU and GLI2. SMO mutations either directly impaired drug binding or activated SMO to varying levels. Furthermore, we found evidence for intra-tumor heterogeneity, suggesting that a combination of therapies targeting components at multiple levels of the Hh pathway is required to overcome resistance., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Tumor-derived cell lines have served as vital models to advance our understanding of oncogene function and therapeutic responses. Although substantial effort has been made to define the genomic constitution of cancer cell line panels, the transcriptome remains understudied. Here we describe RNA sequencing and single-nucleotide polymorphism (SNP) array analysis of 675 human cancer cell lines. We report comprehensive analyses of transcriptome features including gene expression, mutations, gene fusions and expression of non-human sequences. Of the 2,200 gene fusions catalogued, 1,435 consist of genes not previously found in fusions, providing many leads for further investigation. We combine multiple genome and transcriptome features in a pathway-based approach to enhance prediction of response to targeted therapeutics. Our results provide a valuable resource for studies that use cancer cell lines.

Cell lines are the foundation for much of the fundamental research into the mechanisms underlying normal biologic processes and disease mechanisms. It is estimated that 15%-35% of human cell lines are misidentified or contaminated, resulting in a huge waste of resources and publication of false or misleading data. Here we evaluate a panel of 96 single-nucleotide polymorphism (SNP) assays utilizing Fluidigm microfluidics technology for authentication and sex determination of human cell lines. The SNPtrace Panel was tested on 907 human cell lines. Pairwise comparison of these data show the SNPtrace Panel discriminated among identical, related and unrelated pairs of samples with a high degree of confidence, equivalent to short tandem repeat (STR) profiling. We also compared annotated sex calls with those determined by the SNPtrace Panel, STR and Illumina SNP arrays, revealing a high number of male samples are identified as female due to loss of the Y chromosome. Finally we assessed the sensitivity of the SNPtrace Panel to detect intra-human cross-contamination, resulting in detection of as little as 2% contaminating cell population. In conclusion, this study has generated a database of SNP fingerprints for 907 cell lines used in biomedical research and provides a reliable, fast, and economic alternative to STR profiling which can be applied to any human cell line or tissue sample.

To further understand the molecular distinctions between kidney cancer subtypes, we analyzed exome, transcriptome and copy number alteration data from 167 primary human tumors that included renal oncocytomas and non-clear cell renal cell carcinomas (nccRCCs), consisting of papillary (pRCC), chromophobe (chRCC) and translocation (tRCC) subtypes. We identified ten significantly mutated genes in pRCC, including MET, NF2, SLC5A3, PNKD and CPQ. MET mutations occurred in 15% (10/65) of pRCC samples and included previously unreported recurrent activating mutations. In chRCC, we found TP53, PTEN, FAAH2, PDHB, PDXDC1 and ZNF765 to be significantly mutated. Gene expression analysis identified a five-gene set that enabled the molecular classification of chRCC, renal oncocytoma and pRCC. Using RNA sequencing, we identified previously unreported gene fusions, including ACTG1-MITF fusion. Ectopic expression of the ACTG1-MITF fusion led to cellular transformation and induced the expression of downstream target genes. Finally, we observed upregulation of the anti-apoptotic factor BIRC7 in MiTF-high RCC tumors, suggesting a potential therapeutic role for BIRC7 inhibitors.

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., (Copyright © 2014 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.)

Allele-specific gene expression, ASE, is an important aspect of gene regulation. We developed a novel method MBASED, meta-analysis based allele-specific expression detection for ASE detection using RNA-seq data that aggregates information across multiple single nucleotide variation loci to obtain a gene-level measure of ASE, even when prior phasing information is unavailable. MBASED is capable of one-sample and two-sample analyses and performs well in simulations. We applied MBASED to a panel of cancer cell lines and paired tumor-normal tissue samples, and observed extensive ASE in cancer, but not normal, samples, mainly driven by genomic copy number alterations.

Cancer genomes maintain a complex array of somatic alterations required for maintenance and progression of the disease, posing a challenge to identify driver genes among this genetic disorder. Toward this end, we mapped regions of recurrent amplification in a large collection (n=392) of primary human cancers and selected 620 genes whose expression is elevated in tumors. An RNAi loss-of-function screen targeting these genes across a panel of 32 cancer cell lines identified potential driver genes. Subsequent functional assays identified SHMT2, a key enzyme in the serine/glycine synthesis pathway, as necessary for tumor cell survival but insufficient for transformation. The 26S proteasomal subunit, PSMB4, was identified as the first proteasomal subunit with oncogenic properties promoting cancer cell survival and tumor growth in vivo. Elevated expression of SHMT2 and PSMB4 was found to be associated with poor prognosis in human cancer, supporting the development of molecular therapies targeting these genes or components of their pathways., (©2014 American Association for Cancer Research.)

Gastric cancer is the second leading cause of worldwide cancer mortality, yet the underlying genomic alterations remain poorly understood. Here we perform exome and transcriptome sequencing and SNP array assays to characterize 51 primary gastric tumours and 32 cell lines. Meta-analysis of exome data and previously published data sets reveals 24 significantly mutated genes in microsatellite stable (MSS) tumours and 16 in microsatellite instable (MSI) tumours. Over half the patients in our collection could potentially benefit from targeted therapies. We identify 55 splice site mutations accompanied by aberrant splicing products, in addition to mutation-independent differential isoform usage in tumours. ZAK kinase isoform TV1 is preferentially upregulated in gastric tumours and cell lines relative to normal samples. This pattern is also observed in colorectal, bladder and breast cancers. Overexpression of this particular isoform activates multiple cancer-related transcription factor reporters, while depletion of ZAK in gastric cell lines inhibits proliferation. These results reveal the spectrum of genomic and transcriptomic alterations in gastric cancer, and identify isoform-specific oncogenic properties of ZAK.

Introduction: Rheumatoid arthritis (RA) is a complex and clinically heterogeneous autoimmune disease. Currently, the relationship between pathogenic molecular drivers of disease in RA and therapeutic response is poorly understood., Methods: We analyzed synovial tissue samples from two RA cohorts of 49 and 20 patients using a combination of global gene expression, histologic and cellular analyses, and analysis of gene expression data from two further publicly available RA cohorts. To identify candidate serum biomarkers that correspond to differential synovial biology and clinical response to targeted therapies, we performed pre-treatment biomarker analysis compared with therapeutic outcome at week 24 in serum samples from 198 patients from the ADACTA (ADalimumab ACTemrA) phase 4 trial of tocilizumab (anti-IL-6R) monotherapy versus adalimumab (anti-TNFα) monotherapy., Results: We documented evidence for four major phenotypes of RA synovium - lymphoid, myeloid, low inflammatory, and fibroid - each with distinct underlying gene expression signatures. We observed that baseline synovial myeloid, but not lymphoid, gene signature expression was higher in patients with good compared with poor European league against rheumatism (EULAR) clinical response to anti-TNFα therapy at week 16 (P =0.011). We observed that high baseline serum soluble intercellular adhesion molecule 1 (sICAM1), associated with the myeloid phenotype, and high serum C-X-C motif chemokine 13 (CXCL13), associated with the lymphoid phenotype, had differential relationships with clinical response to anti-TNFα compared with anti-IL6R treatment. sICAM1-high/CXCL13-low patients showed the highest week 24 American College of Rheumatology (ACR) 50 response rate to anti-TNFα treatment as compared with sICAM1-low/CXCL13-high patients (42% versus 13%, respectively, P =0.05) while anti-IL-6R patients showed the opposite relationship with these biomarker subgroups (ACR50 20% versus 69%, P =0.004)., Conclusions: These data demonstrate that underlying molecular and cellular heterogeneity in RA impacts clinical outcome to therapies targeting different biological pathways, with patients with the myeloid phenotype exhibiting the most robust response to anti-TNFα. These data suggest a path to identify and validate serum biomarkers that predict response to targeted therapies in rheumatoid arthritis and possibly other autoimmune diseases., Trial Registration: ClinicalTrials.gov NCT01119859

Cancer cells derived from different stages of tumor progression may exhibit distinct biological properties, as exemplified by the paired lung cancer cell lines H1993 and H2073. While H1993 was derived from chemo-naive metastasized tumor, H2073 originated from the chemo-resistant primary tumor from the same patient and exhibits strikingly different drug response profile. To understand the underlying genetic and epigenetic bases for their biological properties, we investigated these cells using a wide range of large-scale methods including whole genome sequencing, RNA sequencing, SNP array, DNA methylation array, and de novo genome assembly. We conducted an integrative analysis of both cell lines to distinguish between potential driver and passenger alterations. Although many genes are mutated in these cell lines, the combination of DNA- and RNA-based variant information strongly implicates a small number of genes including TP53 and STK11 as likely drivers. Likewise, we found a diverse set of genes differentially expressed between these cell lines, but only a fraction can be attributed to changes in DNA copy number or methylation. This set included the ABC transporter ABCC4, implicated in drug resistance, and the metastasis associated MET oncogene. While the rich data content allowed us to reduce the space of hypotheses that could explain most of the observed biological properties, we also caution there is a lack of statistical power and inherent limitations in such single patient case studies.

Purpose: We sought to identify predictive biomarkers for a novel nicotinamide phosphoribosyltransferase (NAMPT) inhibitor., Experimental Design: We use a NAMPT inhibitor, GNE-617, to evaluate nicotinic acid rescue status in a panel of more than 400 cancer cell lines. Using correlative analysis and RNA interference (RNAi), we identify a specific biomarker for nicotinic acid rescue status. We next determine the mechanism of regulation of expression of the biomarker. Finally, we develop immunohistochemical (IHC) and DNA methylation assays and evaluate cancer tissue for prevalence of the biomarker across indications., Results: Nicotinate phosphoribosyltransferase (NAPRT1) is necessary for nicotinic acid rescue and its expression is the major determinant of rescue status. We demonstrate that NAPRT1 promoter methylation accounts for NAPRT1 deficiency in cancer cells, and NAPRT1 methylation is predictive of rescue status in cancer cell lines. Bisulfite next-generation sequencing mapping of the NAPRT1 promoter identified tumor-specific sites of NAPRT1 DNA methylation and enabled the development of a quantitative methylation-specific PCR (QMSP) assay suitable for use on archival formalin-fixed paraffin-embedded tumor tissue., Conclusions: Tumor-specific promoter hypermethylation of NAPRT1 inactivates one of two NAD salvage pathways, resulting in synthetic lethality with the coadministration of a NAMPT inhibitor. NAPRT1 expression is lost due to promoter hypermethylation in most cancer types evaluated at frequencies ranging from 5% to 65%. NAPRT1-specific immunohistochemical or DNA methylation assays can be used on archival formalin paraffin-embedded cancer tissue to identify patients likely to benefit from coadministration of a Nampt inhibitor and nicotinic acid., (©2013 AACR.)

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/.

Background: Although remarkable clinical response rates in melanoma have been observed using vemurafenib or dabrafenib in patients with tumors carrying oncogenic mutations in BRAF, a substantial unmet medical need remains for the subset of patients with wild-type BRAF tumors., Methods: To investigate the role of p21-activated kinases (PAKs) in melanoma, we determined PAK1 genomic copy number and protein expression for a panel of human melanoma tissues. PAK1 was inhibited in vitro and in vivo using RNA interference or PF-3758309 inhibitor treatment in a panel of melanoma cell lines with known BRAF and RAS (rat sarcoma) genotype to better understand its role in melanoma cell proliferation and migration. Tumorigenesis was assessed in vivo in female NCR nude mice and analyzed with cubic spline regression and area under the curve analyses. All statistical tests were two-sided., Results: Strong cytoplasmic PAK1 protein expression was prevalent in melanomas (27%) and negatively associated with activating mutation of the BRAF oncogene (P < .001). Focal copy number gain of PAK1 at 11q13 was also observed in 9% of melanomas (n = 87; copy number ≥ 2.5) and was mutually exclusive with BRAF mutation (P < .005). Selective PAK1 inhibition attenuated signaling through mitogen-activated protein kinase (MAPK) as well as cytoskeleton-regulating pathways to modulate the proliferation and migration of BRAF wild-type melanoma cells. Treatment of BRAF wild-type melanomas with PF-3758309 PAK inhibitor decreased tumor growth for SK-MEL23 and 537MEL xenografts (91% and 63% inhibition, respectively; P < .001) and MAPK pathway activation in vivo., Conclusions: Taken together, our results provide evidence for a functional role of PAK1 in BRAF wild-type melanoma and therapeutic use of PAK inhibitors in this indication.

Oncogenic mutations in PIK3CA, which encodes the phosphoinositide-3-kinase (PI3K) catalytic subunit p110α, occur in ∼25% of human breast cancers. In this study, we report the development of a knock-in mouse model for breast cancer where the endogenous Pik3ca allele was modified to allow tissue-specific conditional expression of a frequently found Pik3ca(H1047R) (Pik3ca(e20H1047R)) mutant allele. We found that activation of the latent Pik3ca(H1047R) allele resulted in breast tumors with multiple histological types. Whole-exome analysis of the Pik3ca(H1047R)-driven mammary tumors identified multiple mutations, including Trp53 mutations that appeared spontaneously during the development of adenocarinoma and spindle cell tumors. Further, we used this model to test the efficacy of GDC-0941, a PI3K inhibitor, in clinical development, and showed that the tumors respond to PI3K inhibition.

Purpose: Class 1 phosphatidylinositol 3-kinase (PI3K) plays a major role in cell proliferation and survival in a wide variety of human cancers. Here, we investigated biomarker strategies for PI3K pathway inhibitors in non-small-cell lung cancer (NSCLC)., Experimental Design: Molecular profiling for candidate PI3K predictive biomarkers was conducted on a collection of NSCLC tumor samples. Assays included comparative genomic hybridization, reverse-transcription polymerase chain reaction gene expression, mutation detection for PIK3CA and other oncogenes, PTEN immunohistochemistry, and FISH for PIK3CA copy number. In addition, a panel of NSCLC cell lines characterized for alterations in the PI3K pathway was screened with PI3K and dual PI3K/mTOR inhibitors to assess the preclinical predictive value of candidate biomarkers., Results: PIK3CA amplification was detected in 37% of squamous tumors and 5% of adenocarcinomas, whereas PIK3CA mutations were found in 9% of squamous and 0% of adenocarcinomas. Total loss of PTEN immunostaining was found in 21% of squamous tumors and 4% of adenocarcinomas. Cell lines harboring pathway alterations (receptor tyrosine kinase activation, PI3K mutation or amplification, and PTEN loss) were exquisitely sensitive to the PI3K inhibitor GDC-0941. A dual PI3K/mTOR inhibitor had broader activity across the cell line panel and in tumor xenografts. The combination of GDC-0941 with paclitaxel, erlotinib, or a mitogen-activated protein-extracellular signal-regulated kinase inhibitor had greater effects on cell viability than PI3K inhibition alone., Conclusions: Candidate biomarkers for PI3K inhibitors have predictive value in preclinical models and show histology-specific alterations in primary tumors, suggesting that distinct biomarker strategies may be required in squamous compared with nonsquamous NSCLC patient populations., (©2012 AACR.)

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.

Small-cell lung cancer (SCLC) is an exceptionally aggressive disease with poor prognosis. Here, we obtained exome, transcriptome and copy-number alteration data from approximately 53 samples consisting of 36 primary human SCLC and normal tissue pairs and 17 matched SCLC and lymphoblastoid cell lines. We also obtained data for 4 primary tumors and 23 SCLC cell lines. We identified 22 significantly mutated genes in SCLC, including genes encoding kinases, G protein-coupled receptors and chromatin-modifying proteins. We found that several members of the SOX family of genes were mutated in SCLC. We also found SOX2 amplification in ∼27% of the samples. Suppression of SOX2 using shRNAs blocked proliferation of SOX2-amplified SCLC lines. RNA sequencing identified multiple fusion transcripts and a recurrent RLF-MYCL1 fusion. Silencing of MYCL1 in SCLC cell lines that had the RLF-MYCL1 fusion decreased cell proliferation. These data provide an in-depth view of the spectrum of genomic alterations in SCLC and identify several potential targets for therapeutic intervention.

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.

The RAS/RAF/MEK pathway is activated in more than 30% of human cancers, most commonly via mutation in the K-ras oncogene and also via mutations in BRAF. Several allosteric mitogen-activated protein/extracellular signal-regulated kinase (MEK) inhibitors, aimed at treating tumors with RAS/RAF pathway alterations, are in clinical development. However, acquired resistance to these inhibitors has been documented both in preclinical and clinical samples. To identify strategies to overcome this resistance, we have derived three independent MEK inhibitor-resistant cell lines. Resistance to allosteric MEK inhibitors in these cell lines was consistently linked to acquired mutations in the allosteric binding pocket of MEK. In one cell line, concurrent amplification of mutant K-ras was observed in conjunction with MEK allosteric pocket mutations. Clonal analysis showed that both resistance mechanisms occur in the same cell and contribute to enhanced resistance. Importantly, in all cases the MEK-resistant cell lines retained their addiction to the mitogen-activated protein kinase (MAPK) pathway, as evidenced by their sensitivity to a selective inhibitor of the ERK1/2 kinases. These data suggest that tumors with acquired MEK inhibitor resistance remain dependent on the MAPK pathway and are therefore sensitive to inhibitors that act downstream of the mutated MEK target. Importantly, we show that dual inhibition of MEK and ERK by small molecule inhibitors was synergistic and acted to both inhibit the emergence of resistance, as well as to overcome acquired resistance to MEK inhibitors. Therefore, our data provide a rationale for cotargeting multiple nodes within the MAPK signaling cascade in K-ras mutant tumors to maximize therapeutic benefit for patients., (©2012 AACR)

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.

The PTEN/PI3K pathway is commonly mutated in cancer and therefore represents an attractive target for therapeutic intervention. To investigate the primary phenotypes mediated by increased pathway signaling in a clean, patient-relevant context, an activating PIK3CA mutation (H1047R) was knocked-in to an endogenous allele of the MCF10A non-tumorigenic human breast epithelial cell line. Introduction of an endogenously mutated PIK3CA allele resulted in a marked epithelial-mesenchymal transition (EMT) and invasive phenotype, compared to isogenic wild-type cells. The invasive phenotype was linked to enhanced PIP(3) production via a S6K-IRS positive feedback mechanism. Moreover, potent and selective inhibitors of PI3K were highly effective in reversing this phenotype, which is optimally revealed in 3-dimensional cell culture. In contrast, inhibition of Akt or mTOR exacerbated the invasive phenotype. Our results suggest that invasion is a core phenotype mediated by increased PTEN/PI3K pathway activity and that therapeutic agents targeting different nodes of the PI3K pathway may have dramatic differences in their ability to reverse or promote cancer metastasis.

The HER2 oncogene is overexpressed or amplified in 20% of breast cancers. HER2-positive cancer historically portends a poor prognosis, but the HER2-targeted therapy trastuzumab mitigates this otherwise ominous distinction. Nevertheless, some patients suffer disease recurrence despite trastuzumab, and metastatic disease remains largely incurable due to innate and acquired resistance. Thus, understanding trastuzumab resistance remains an unmet medical need. Through RNA interference screening, we discovered that knockdown of the serine/threonine phosphatase PPM1H confers trastuzumab resistance via reduction in protein levels of the tumor suppressor p27. PPM1H dephosphorylates p27 at threonine 187, thus removing a signal for proteasomal degradation. We further determined that patients whose tumors express low levels of PPM1H trend towards worse clinical outcome on trastuzumab. Identifying PPM1H as a novel p27 phosphatase reveals new insight into how cancer cells destabilize a well-recognized tumor suppressor. Furthermore, low PPM1H expression may identify a subset of HER2-positive tumors that are harder to treat., (©2011 AACR.)

The p21-activated kinases (PAKs) are central players in growth factor signaling networks and morphogenetic processes that control proliferation, cell polarity, invasion and actin cytoskeleton organization. This raises the possibility that interfering with PAK activity may produce significant anti-tumor activity. In this perspective, we summarize recent data concerning the contribution of the PAK family member, PAK1, in growth factor signaling and tumorigenesis. We further discuss mechanisms by which inhibition of PAK1 can arrest tumor growth and promote cell apoptosis, and the types of cancers in which PAK1 inhibition may hold promise.

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.

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.

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.

Purpose: Oncogenic activation of the phosphatidylinositol 3-kinase (PI3K) signaling pathway is prevalent in breast cancer and has been associated with resistance to HER2 inhibitors in the clinic. We therefore investigated the combinatorial activity of GDC-0941, a novel class I PI3K inhibitor, with standard-of-care therapies for HER2-amplified breast cancer., Experimental Design: Three-dimensional laminin-rich extracellular matrix cultures of human breast cancer cells were utilized to provide a physiologically relevant approach to analyze the efficacy and molecular mechanism of combination therapies ex vivo. Combination studies were done using GDC-0941 with trastuzumab (Herceptin), pertuzumab, lapatinib (Tykerb), and docetaxel, the principal therapeutic agents that are either approved or being evaluated for treatment of early HER2-positive breast cancer., Results: Significant GDC-0941 activity (EC(50) <1 micromol/L) was observed for >70% of breast cancer cell lines that were examined in three-dimensional laminin-rich extracellular matrix culture. Differential responsiveness to GDC-0941 as a single agent was observed for luminal breast cancer cells upon stimulation with the HER3 ligand, heregulin. Combined treatment of GDC-0941, trastuzumab, and pertuzumab resulted in growth inhibition, altered acinar morphology, and suppression of AKT mitogen-activated protein kinase (MAPK) / extracellular signed-regulated kinase (ERK) kinase and MEK effector signaling pathways for HER2-amplified cells in both normal and heregulin-supplemented media. The GDC-0941 and lapatinib combination further showed that inhibition of HER2 activity was essential for maximum combinatorial efficacy. PI3K inhibition also rendered HER2-amplified BT-474M1 cells and tumor xenografts more sensitive to docetaxel., Conclusions: GDC-0941 is efficacious in preclinical models of breast cancer. The addition of GDC-0941 to HER2-directed treatment could augment clinical benefit in breast cancer patients.

Background: DNA copy number alterations are frequently observed in ovarian cancer, but it remains a challenge to identify the most relevant alterations and the specific causal genes in those regions., Methods: We obtained high-resolution 500K SNP array data for 52 ovarian tumors and identified the most statistically significant minimal genomic regions with the most prevalent and highest-level copy number alterations (recurrent CNAs). Within a region of recurrent CNA, comparison of expression levels in tumors with a given CNA to tumors lacking that CNA and to whole normal ovary samples was used to select genes with CNA-specific expression patterns. A public expression array data set of laser capture micro-dissected (LCM) non-malignant fallopian tube epithelia and LCM ovarian serous adenocarcinoma was used to evaluate the effect of cell-type mixture biases., Results: Fourteen recurrent deletions were detected on chromosomes 4, 6, 9, 12, 13, 15, 16, 17, 18, 22 and most prevalently on X and 8. Copy number and expression data suggest several apoptosis mediators as candidate drivers of the 8p deletions. Sixteen recurrent gains were identified on chromosomes 1, 2, 3, 5, 8, 10, 12, 15, 17, 19, and 20, with the most prevalent gains localized to 8q and 3q. Within the 8q amplicon, PVT1, but not MYC, was strongly over-expressed relative to tumors lacking this CNA and showed over-expression relative to normal ovary. Likewise, the cell polarity regulators PRKCI and ECT2 were identified as putative drivers of two distinct amplicons on 3q. Co-occurrence analyses suggested potential synergistic or antagonistic relationships between recurrent CNAs. Genes within regions of recurrent CNA showed an enrichment of Cancer Census genes, particularly when filtered for CNA-specific expression., Conclusion: These analyses provide detailed views of ovarian cancer genomic changes and highlight the benefits of using multiple reference sample types for the evaluation of CNA-specific expression changes.

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.

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.

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.

The adenosine A2A receptor (A2AR) is highly expressed in the striatum, where it modulates motor and emotional behaviors. We used both microarray and bioinformatics analyses to compare gene expression profiles by genetic and pharmacological inactivation of A2AR and inferred an A2AR-controlled transcription network in the mouse striatum. A comparison between vehicle (VEH)-treated A2AR knockout (KO) mice (A2AR KO-VEH) and wild-type (WT) mice (WT-VEH) revealed 36 upregulated genes that were partially mimicked by treatment with SCH-58261 (SCH; an A2AR antagonist) and 54 downregulated genes that were not mimicked by SCH treatment. We validated the A2AR as a specific drug target for SCH by comparing A2AR KO-SCH and A2AR KO-VEH groups. The unique downregulation effect of A2AR KO was confirmed by comparing A2AR KO-SCH with WT-SCH gene groups. The distinct striatal gene expression profiles induced by A2AR KO and SCH should provide clues to the molecular mechanisms underlying the different phenotypes observed after genetic and pharmacological inactivation of A2AR. Furthermore, bioinformatics analysis discovered that Egr-2 binding sites were statistically overrepresented in the proximal promoters of A2AR KO-affected genes relative to the unaffected genes. This finding was further substantiated by the demonstration that the Egr-2 mRNA level increased in the striatum of both A2AR KO and SCH-treated mice and that striatal Egr-2 binding activity in the promoters of two A2AR KO-affected genes was enhanced in A2AR KO mice as assayed by chromatin immunoprecipitation. Taken together, these results strongly support the existence of an Egr-2-directed transcriptional regulatory network controlled by striatal A2ARs.

Motivation: DNA microarrays have revolutionized biological research, but their reliability and accuracy have not been extensively evaluated. Thorough testing of microarrays through comparison to dissimilar gene expression methods is necessary in order to determine their accuracy., Results: We have systematically compared three global gene expression methods on all available histologically normal samples from five human organ types. The data included 25 Affymetrix high-density oligonucleotide array experiments, 23 expressed sequence tag based expression (EBE) experiments and 5 SAGE experiments. The reported gene-by-gene expression patterns showed a wide range of correlations between pairs of methods. This level of agreement was sufficient for accurate clustering of datasets from the same tissue and dissimilar methods, but highlights the need for thorough validation of individual gene expression measurements by alternate, non-global methods. Furthermore, analyses of mRNA abundance distributions indicate limitations in the EBE and SAGE methods at both high- and low-expression levels.

Summary: CisML is an XML-based format for sequence motif detection software. This proposed standard is applicable to many types of sequence motif detection programs. It is intended to facilitate the integration of data and the comparison of results from different software packages, and to simplify the development of downstream tools. XSL stylesheets are provided for easy generation of text, html and graphical reports from CisML-formatted data., Availability: http://zlab.bu.edu/CisML/ SUPPLEMENTARY INFORMATION: Example CisML-formatted data and XSL stylesheets for report generation are available along with the sample output.

Detecting overrepresented known transcription factor binding motifs in a set of promoter sequences of co-regulated genes has become an important approach to deciphering transcriptional regulatory mechanisms. In this paper, we present an interactive web server, MotifViz, for three motif discovery programs, Clover, Rover and Motifish, covering most available flavors of algorithms for achieving this goal. For comparison, we have also implemented the simple motif-matching program Possum. MotifViz provides uniform and intuitive input and output formats for all four programs. It can be accessed at http://biowulf.bu.edu/MotifViz.