Search

Metastasis is the principal cause of cancer related deaths. Tumor invasion is essential for metastatic spread. However, determinants of invasion are poorly understood. We addressed this knowledge gap by leveraging a unique attribute of kidney cancer. Renal tumors invade into large vessels forming tumor thrombi (TT) that migrate extending sometimes into the heart. Over a decade, we prospectively enrolled 83 ethnically-diverse patients undergoing surgical resection for grossly invasive tumors at UT Southwestern Kidney Cancer Program. In this study, we perform comprehensive histological analyses, integrate multi-region genomic studies, generate in vivo models, and execute functional studies to define tumor invasion and metastatic competence. We find that invasion is not always associated with the most aggressive clone. Driven by immediate early genes, invasion appears to be an opportunistic trait attained by subclones with diverse oncogenomic status in geospatial proximity to vasculature. We show that not all invasive tumors metastasize and identify determinants of metastatic competency. TT associated with metastases are characterized by higher grade, mTOR activation and a particular immune contexture. Moreover, TT grade is a better predictor of metastasis than overall tumor grade, which may have implications for clinical practice., Tumour thrombi (TT) are intravascular extensions of renal cell carcinomas (RCC) which often lead to distant metastases. Here the authors examine the determinants of vascular invasion and metastasis in a unique cohort of RCC patients with TT using multi-region genomics and in vivo models.

By leveraging tumorgraft (patient-derived xenograft) RNA-sequencing data, we developed an empirical approach, DisHet, to dissect the tumor microenvironment (eTME). We found that 65% of previously defined immune signature genes are not abundantly expressed in renal cell carcinoma (RCC) and identified 610 novel immune/stromal transcripts. Using eTME, genomics, pathology, and medical record data involving >1,000 patients, we established an inflamed pan-RCC subtype (IS) enriched for regulatory T cells, natural killer cells, TH1 cells, neutrophils, macrophages, B cells, and CD8+ T cells. IS is enriched for aggressive RCCs, including BAP1-deficient clear-cell and type 2 papillary tumors. The IS subtype correlated with systemic manifestations of inflammation such as thrombocytosis and anemia, which are enigmatic predictors of poor prognosis. Furthermore, IS was a strong predictor of poor survival. Our analyses suggest that tumor cells drive the stromal immune response. These data provide a missing link between tumor cells, the TME, and systemic factors. Significance: We undertook a novel empirical approach to dissect the renal cell carcinoma TME by leveraging tumorgrafts. The dissection and downstream analyses uncovered missing links between tumor cells, the TME, systemic manifestations of inflammation, and poor prognosis. Cancer Discov; 8(9); 1142–55. ©2018 AACR. This article is highlighted in the In This Issue feature, p. 1047

A propensity for rewiring genetic and epigenetic regulatory networks, thus enabling sustained cell proliferation, suppression of apoptosis, and the ability to evade the immune system, is vital to cancer cell propagation. An increased understanding of how this is achieved is critical for identifying or improving therapeutic interventions. In this study, using acute myeloid leukemia (AML) human cell lines and a custom CRISPR/Cas9 screening platform, we identify the H3K9 methyltransferase SETDB1 as a novel, negative regulator of innate immunity. SETDB1 is overexpressed in many cancers, and loss of this gene in AML cells triggers desilencing of retrotransposable elements that leads to the production of double-stranded RNAs (dsRNAs). This is coincident with induction of a type I interferon response and apoptosis through the dsRNA-sensing pathway. Collectively, our findings establish a unique gene regulatory axis that cancer cells can exploit to circumvent the immune system.

The cellular origin of sporadic pancreatic neuroendocrine tumors (PNETs) is obscure. Dysregulated hormone production and the absence of obvious driver mutations mask tumor origins although they presumably arise from endocrine cell lineages such as glucagon-producing alpha cells or insulin-expressing beta cells. A reciprocal feedback circuit involving pancreatic glucagon and hepatic amino acids regulates alpha cell mass such that loss of hepatic glucagon receptor (GCGR) signaling drives alpha cell hyperproliferation and tumor formation. This genetically-defined system is an ideal model for parsing the mechanisms of PNET origin and progression. Utilizing a combination of genetically engineered Gcgr knockout mice and GCGR-inhibiting antibodies, we found that elevated plasma amino acids drive reversion to a proliferative SLC38A5+ embryonic progenitor-like status. Chronic activation of this mTOR-dependent state drives tumor initiation and single cell profiling reveals tumors originate from a subpopulation of alpha cells. Profiling RNA, tumor signaling, and immune content of biopsies from rare human PNETs found in individuals with germ-line GCGR loss-of-function variants showed active mTOR and prominent tumor cell-associated expression of amino acid transporter SLC7A8. This work, for the first time, identifies metabolic control of an adult-to-progenitor alpha cell identity switch and shows that chronic maintenance of progenitor identity in the adult organ drives proliferation and tumor initiation.

Background PIK3CA mutations are frequent in human breast cancer. Pik3caH1047R mutant expression in mouse mammary gland promotes tumorigenesis. TP53 mutations co-occur with PIK3CA mutations in human breast cancers. We previously generated a conditionally activatable Pik3caH1047R;MMTV-Cre mouse model and found a few malignant sarcomatoid (spindle cell) carcinomas that had acquired spontaneous dominant-negative Trp53 mutations. Methods A Pik3caH1047R;Trp53R270H;MMTV-Cre double mutant mouse breast cancer model was generated. Tumors were characterized by histology, marker analysis, transcriptional profiling, single-cell RNA-seq, and bioinformatics. Cell lines were developed from mutant tumors and used to identify and confirm genes involved in metastasis. Results We found Pik3caH1047R and Trp53R270H cooperate in driving oncogenesis in mammary glands leading to a shorter latency than either alone. Double mutant mice develop multiple histologically distinct mammary tumors, including adenocarcinoma and sarcomatoid (spindle cell) carcinoma. We found some tumors to be invasive and a few metastasized to the lung and/or the lymph node. Single-cell RNA-seq analysis of the tumors identified epithelial, stromal, myeloid, and T cell groups. Expression analysis of the metastatic tumors identified S100a4 as a top candidate gene associated with metastasis. Metastatic tumors contained a much higher percentage of epithelial–mesenchymal transition (EMT)-signature positive and S100a4-expressing cells. CRISPR/CAS9-mediated knockout of S100a4 in a metastatic tumor-derived cell line disrupted its metastatic potential indicating a role for S100a4 in metastasis. Conclusions Pik3caH1047R;Trp53R270H;MMTV-Cre mouse provides a preclinical model to mimic a subtype of human breast cancers that carry both PIK3CA and TP53 mutations. It also allows for understanding the cooperation between the two mutant genes in tumorigenesis. Our model also provides a system to study metastasis and develop therapeutic strategies for PIK3CA/TP53 double-positive cancers. S100a4 found involved in metastasis in this model can be a potential diagnostic and therapeutic target.

We report the unique pathogenesis and presentation of a rapidly progressive B-cell lymphoma in a 3-year-old female cynomolgus monkey on day 50 of a 13-week toxicity study. Clinical pathology evaluation revealed a marked leukocytosis with bicytopenia. A serum protein electrophoresis was consistent with monoclonal gammopathy. The architecture of the lymph node, spleen, and thymus were variably effaced by neoplastic cells, which also infiltrated other tissues. Immunohistochemistry of the affected tissues confirmed a predominant population of CD20+, CD79a+, CD3-, CD68-, and CD34-neoplastic cells. The full data best support a diagnosis of Stage V lymphoma. Nextgen sequencing and negative prestudy serology results suggested a recent infection by macaque lymphocryptovirus (mLCV) with a unique transcriptional profile comparable with a rarely observed direct LCV infection model. This infection model might be associated with a temporary lack of an LCV antigen-specific cytotoxic T-cell adaptive immune response. Consistent with the established mechanisms of LCV-related lymphoproliferation, MYC and BCL2L11 gene expression were increased and decreased, respectively. While there was no overt immunosuppression, immunophenotyping revealed the index animal had a relatively low NK cell count, which further decreased by >50% on day 24 of the study. In addition to the temporary lack of adaptive immunity, the low NK cell counts were suggestive of an impaired innate immunity to control the virally-transformed cells and the subsequent unchecked lymphoproliferation. To our knowledge, this is the first report of a Stage V lymphoma with a unique pathogenesis in an otherwise immunocompetent cynomolgus monkey.

Given the frequency with which MAP kinase signaling is dysregulated in cancer, much effort has been focused on inhibiting RAS signaling for therapeutic benefit. KSR1, a pseudokinase that interacts with RAF, is a potential target; it was originally cloned in screens for suppressors of constitutively active RAS, and its deletion prevents RAS-mediated transformation of mouse embryonic fibroblasts. In this work, we used a genetically engineered mouse model of pancreatic cancer to assess whether KSR1 deletion would influence tumor development in the setting of oncogenic RAS. We found that Ksr1-/- mice on this background had a modest but significant improvement in all-cause morbidity compared to Ksr1+/+ and Ksr1+/- cohorts. Ksr1-/- mice, however, still developed tumors, and precursor pancreatic intraepithelial neoplastic (PanIN) lesions were detected within a similar timeframe compared to Ksr1+/+ mice. No significant differences in pERK expression or in proliferation were noted. RNA sequencing also did not reveal any unique genetic signature in Ksr1-/- tumors. Further studies will be needed to determine whether and in what settings KSR inhibition may be clinically useful.

Colorectal cancer remains a major unmet medical need, prompting large-scale genomics efforts in the field to identify molecular drivers for which targeted therapies might be developed. We previously reported the identification of recurrent translocations in R-spondin genes present in a subset of colorectal tumours. Here we show that targeting RSPO3 in PTPRK-RSPO3-fusion-positive human tumour xenografts inhibits tumour growth and promotes differentiation. Notably, genes expressed in the stem-cell compartment of the intestine were among those most sensitive to anti-RSPO3 treatment. This observation, combined with functional assays, suggests that a stem-cell compartment drives PTPRK-RSPO3 colorectal tumour growth and indicates that the therapeutic targeting of stem-cell properties within tumours may be a clinically relevant approach for the treatment of colorectal tumours.

Identifying and understanding changes in cancer genomes is essential for the development of targeted therapeutics1. 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.

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.

The panel of 60 human cancer cell lines (the NCI-60) assembled by the National Cancer Institute for anticancer drug discovery is a widely used resource. We previously sequenced 24 cancer genes in those cell lines. Eleven of the genes were found to be mutated in three or more of the lines. Using a pharmacogenomic approach, we analyzed the relationship between drug activity and mutations in those 11 genes (APC, RB1, KRAS, NRAS, BRAF, PIK3CA, PTEN, STK11, MADH4, TP53, and CDKN2A). That analysis identified an association between mutation in BRAF and the antiproliferative potential of phenothiazine compounds. Phenothiazines have been used as antipsychotics and as adjunct antiemetics during cancer chemotherapy and more recently have been reported to have anticancer properties. However, to date, the anticancer mechanism of action of phenothiazines has not been elucidated. To follow up on the initial pharmacologic observations in the NCI-60 screen, we did pharmacologic experiments on 11 of the NCI-60 cell lines and, prospectively, on an additional 24 lines. The studies provide evidence that BRAF mutation (codon 600) in melanoma as opposed to RAS mutation is predictive of an increase in sensitivity to phenothiazines as determined by 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium, inner salt assay (Wilcoxon P = 0.007). That pattern of increased sensitivity to phenothiazines based on the presence of codon 600 BRAF mutation may be unique to melanomas, as we do not observe it in a panel of colorectal cancers. The findings reported here have potential implications for the use of phenothiazines in the treatment of V600E BRAF mutant melanoma. [Mol Cancer Ther 2008;7(6):1337–46]

During the type-setting of the final version of the article [1] some of the additional files were swapped. The correct files are republished in this Erratum.

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.

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., The genetic basis of gastric cancer, the fourth most common cancer worldwide, remains poorly understood. Here, the authors sequence and analyse the exomes and transcriptomes of primary gastric tumours and cell lines, and identify a ZAK kinase isoform that may have an oncogenic role in gastric cancer.

Malignant pleural mesothelioma (MPM) is an aggressive cancer arising from the mesothelial cells of the pleura. About 80% of mesothelioma cases are linked to asbestos exposure, while the remainder may be related to prior chest radiation, genetic predisposition or spontaneous occurrence. We analyzed transcriptomes (n = 211), whole exomes (n = 99), and targeted exomes (n = 103) from 216 malignant pleural mesothelioma (MPM) tumors. Four distinct molecular subtypes, sarcomatoid (S), epithelioid (E), biphasic-E, and biphasic-S were identified using RNA-seq data. Using exome analysis we found BAP1, NF2, TP53, SETD2, DDX3X, ULK2, RYR2, CFAP45, SETDB1, and DDX51 to be significantly mutated MPM genes. We identified recurrent mutations in several genes including splice complex factor SF3B1 (∼2% [4/216]) and TRAF7 (∼2% [5/216]). SF3B1 mutant samples showed a distinct splicing profile compared to wild-type tumors. Mutations in TRAF7 occurred primarily in the WD40 domain and except in one case were mutually exclusive with NF2 mutations. We found recurrent gene fusions and splice alterations to be frequent mechanisms for inactivation of NF2, BAP1 and SETD2. We analyzed the tumors for infiltrating immune cells and found that M2 macrophage to T-cell ratio was predictive of reduced overall survival. Further we found sarcomatoid subtype had a significantly higher level of PD-L1 expression suggesting that this group might be prioritized for checkpoint inhibitor therapy. Integrated analysis of alterations identified changes in Hippo, mTOR, Histone methylation, RNA helicase and p53 signaling pathways in MPMs. Overall our study significantly expands on the previous genomic studies and provides a comprehensive genomics profile of mesothelioma and identifies several previously unknown alterations and biological pathways that are potential targets for drug discovery and treatment. Incorporating genomic analysis for detection of actionable alterations as part of MPM patient care will help in developing rational individualized therapy. Citation Format: Somasekar Seshagiri, Eric Stawiski, Leonard Goldstein, Steffen Durinck, Zora Modrusan, Florian Gnad, Thong T. Nguyen, Bijay Jaiswal, Assunta De Rienzo, William Richards, Raphael Bueno. Comprehensive genomic analysis of malignant pleural mesothelioma identifies recurrent mutations, gene fusions and splicing alterations. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr LB-331.

International audience; Breast cancers are comprised of molecularly distinct subtypes that may respond differently to pathway-targeted therapies now under development. Collections of breast cancer cell lines mirror many of the molecular subtypes and pathways found in tumors, suggesting that treatment of cell lines with candidate therapeutic compounds can guide identification of associations between molecular subtypes, pathways, and drug response. In a test of 77 therapeutic compounds, nearly all drugs showed differential responses across these cell lines, and approximately one third showed subtype-, pathway-, and/or genomic aberration-specific responses. These observations suggest mechanisms of response and resistance and may inform efforts to develop molecular assays that predict clinical response.

Timely intervention for cancer requires knowledge of its earliest genetic aberrations. Sequencing of tumors and their metastases reveals numerous abnormalities occurring late in progression. A means to temporally order aberrations in a single cancer, rather than inferring them from serially acquired samples, would define changes preceding even clinically evident disease. We integrate DNA sequence and copy number information to reconstruct the order of abnormalities as individual tumors evolve for 2 separate cancer types. We detect vast, unreported expansion of simple mutations sharply demarcated by recombinative loss of the second copy of TP53 in cutaneous squamous cell carcinomas (cSCC) and serous ovarian adenocarcinomas, in the former surpassing 50 mutations per megabase. In cSCCs, we also report diverse secondary mutations in known and novel oncogenic pathways, illustrating how such expanded mutagenesis directly promotes malignant progression. These results reframe paradigms in which TP53 mutation is required later, to bypass senescence induced by driver oncogenes. Significance: Our approach reveals sequential ordering of oncogenic events in individual cancers, based on chromosomal rearrangements. Identifying the earliest abnormalities in cancer represents a critical step in timely diagnosis and deployment of targeted therapeutics. Cancer Discovery; 1(2); 137–43. © 2011 AACR. This article is highlighted in the In This Issue feature, p. 91

Purpose: microRNAs have been shown to be involved in different human cancers. We therefore have performed expression profiles on a panel of pediatric tumors to identify cancer-specific microRNAs. We also investigated if microRNAs are coregulated with their host gene. Experimental Design: We performed parallel microRNAs and mRNA expression profiling on 57 tumor xenografts and cell lines representing 10 different pediatric solid tumors using microarrays. For those microRNAs that map to their host mRNA, we calculated correlations between them. Results: We found that the majority of cancer types clustered together based on their global microRNA expression profiles by unsupervised hierarchical clustering. Fourteen microRNAs were significantly differentially expressed between rhabdomyosarcoma and neuroblastoma, and 8 of them were validated in independent patient tumor samples. Exploration of the expression of microRNAs in relationship with their host genes showed that the expression for 43 of 68 (63%) microRNAs located inside known coding genes was significantly correlated with that of their host genes. Among these 43 microRNAs, 5 of 7 microRNAs in the OncomiR-1 cluster correlated significantly with their host gene MIRHG1 (P < 0.01). In addition, high expression of MIRHG1 was significantly associated with high stage and MYCN amplification in neuroblastoma tumors, and the expression level of MIRHG1 could predict the outcome of neuroblastoma patients independently from the current neuroblastoma risk-stratification in two independent patient cohorts. Conclusion: Pediatric cancers express cancer-specific microRNAs. The high expression of the OncomiR-1 host gene MIRHG1 correlates with poor outcome for patients with neuroblastoma, indicating important oncogenic functions of this microRNA cluster in neuroblastoma biology. (Clin Cancer Res 2009;15(17):5560–8)

Loss of 1p36 heterozygosity commonly occurs with MYCN amplification in neuroblastoma tumors, and both are associated with an aggressive phenotype. Database searches identified five microRNAs that map to the commonly deleted region of 1p36 and we hypothesized that the loss of one or more of these microRNAs contributes to the malignant phenotype of MYCN-amplified tumors. By bioinformatic analysis, we identified that three out of the five microRNAs target MYCN and of these miR-34a caused the most significant suppression of cell growth through increased apoptosis and decreased DNA synthesis in neuroblastoma cell lines with MYCN amplification. Quantitative RT-PCR showed that neuroblastoma tumors with 1p36 loss expressed lower level of miR-34a than those with normal copies of 1p36. Furthermore, we demonstrated that MYCN is a direct target of miR-34a. Finally, using a series of mRNA expression profiling experiments, we identified other potential direct targets of miR-34a, and pathway analysis demonstrated that miR-34a suppresses cell-cycle genes and induces several neural-related genes. This study demonstrates one important regulatory role of miR-34a in cell growth and MYCN suppression in neuroblastoma.

The human epidermal growth factor receptor (HER) family of tyrosine kinases are deregulated in multiple cancers either through amplification, overexpression or mutation. ERBB3/HER3, the only ERBB family member with an impaired kinase domain, upon ligand binding heterodimerizes with ERBB2 to promote signaling. While amplification and overexpression of ERBB3 is observed in some cancers, occurrence and relevance of ERBB3 somatic mutations in oncogenesis is not established. Here we report the identification of ERBB3 somatic mutations in ∼11% of colon and gastric cancers. We found that the ERBB3 mutants together with ERBB2 promote oncogenic signaling and transformed colonic and breast epithelial cells in a ligand independent manner. Further, we found that multiple target therapeutics that acts on ERBB3, ERBB2 or their downstream signaling components are effective in blocking ERBB3-mutant mediated oncogenic signaling and disease progression in vivo. Identification of activating ERBB3 mutations along with the established contribution of ERBB3 to acquired resistance to EGFR/ERBB2-targeted drugs provides a rational to therapeutically target ERBB3. Citation Format: Bijay S. Jaiswal, Noelyn M. Kljavin, Eric Stawiski, Steffen Durinck, Subhra Chaudhuri, Charles Eigenbrot, Gabriele Schaefer, Frederic J. de Sauvage, Somasekar Seshagiri. Oncogenic ERBB3 mutations in human cancers. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 4428. doi:10.1158/1538-7445.AM2014-4428

Background First-generation molecular profiles for human breast cancers have enabled the identification of features that can predict therapeutic response; however, little is known about how the various data types can best be combined to yield optimal predictors. Collections of breast cancer cell lines mirror many aspects of breast cancer molecular pathobiology, and measurements of their omic and biological therapeutic responses are well-suited for development of strategies to identify the most predictive molecular feature sets. Results We used least squares-support vector machines and random forest algorithms to identify molecular features associated with responses of a collection of 70 breast cancer cell lines to 90 experimental or approved therapeutic agents. The datasets analyzed included measurements of copy number aberrations, mutations, gene and isoform expression, promoter methylation and protein expression. Transcriptional subtype contributed strongly to response predictors for 25% of compounds, and adding other molecular data types improved prediction for 65%. No single molecular dataset consistently out-performed the others, suggesting that therapeutic response is mediated at multiple levels in the genome. Response predictors were developed and applied to TCGA data, and were found to be present in subsets of those patient samples. Conclusions These results suggest that matching patients to treatments based on transcriptional subtype will improve response rates, and inclusion of additional features from other profiling data types may provide additional benefit. Further, we suggest a systems biology strategy for guiding clinical trials so that patient cohorts most likely to respond to new therapies may be more efficiently identified.