Integrating biological knowledge into variable selection : an empirical Bayes approach with an application in cancer biology

Background: An important question in the analysis of biochemical data is that of identifying subsets of molecular variables that may jointly influence a biological response. Statistical variable selection methods have been widely used for this purpose. In many settings, it may be important to incorporate ancillary biological information concerning the variables of interest. Pathway and network maps are one example of a source of such information. However, although ancillary information is increasingly available, it is not always clear how it should be used nor how it should be weighted in relation to primary data. Results: We put forward an approach in which biological knowledge is incorporated using informative prior distributions over variable subsets, with prior information selected and weighted in an automated, objective manner using an empirical Bayes formulation. We employ continuous, linear models with interaction terms and exploit biochemically-motivated sparsity constraints to permit exact inference. We show an example of priors for pathway- and network-based information and illustrate our proposed method on both synthetic response data and by an application to cancer drug response data. Comparisons are also made to alternative Bayesian and frequentist penalised-likelihood methods for incorporating network-based information. Conclusions: The empirical Bayes method proposed here can aid prior elicitation for Bayesian variable selection studies and help to guard against mis-specification of priors. Empirical Bayes, together with the proposed pathway-based priors, results in an approach with a competitive variable selection performance. In addition, the overall procedure is fast, deterministic, and has very few user-set parameters, yet is capable of capturing interplay between molecular players. The approach presented is general and readily applicable in any setting with multiple sources of biological prior knowledge

Genome co-amplification upregulates a mitotic gene network activity that predicts outcome and response to mitotic protein inhibitors in breast cancer.

BACKGROUND: High mitotic activity is associated with the genesis and progression of many cancers. Small molecule inhibitors of mitotic apparatus proteins are now being developed and evaluated clinically as anticancer agents. With clinical trials of several of these experimental compounds underway, it is important to understand the molecular mechanisms that determine high mitotic activity, identify tumor subtypes that carry molecular aberrations that confer high mitotic activity, and to develop molecular markers that distinguish which tumors will be most responsive to mitotic apparatus inhibitors. METHODS: We identified a coordinately regulated mitotic apparatus network by analyzing gene expression profiles for 53 malignant and non-malignant human breast cancer cell lines and two separate primary breast tumor datasets. We defined the mitotic network activity index (MNAI) as the sum of the transcriptional levels of the 54 coordinately regulated mitotic apparatus genes. The effect of those genes on cell growth was evaluated by small interfering RNA (siRNA). RESULTS: High MNAI was enriched in basal-like breast tumors and was associated with reduced survival duration and preferential sensitivity to inhibitors of the mitotic apparatus proteins, polo-like kinase, centromere associated protein E and aurora kinase designated GSK462364, GSK923295 and GSK1070916, respectively. Co-amplification of regions of chromosomes 8q24, 10p15-p12, 12p13, and 17q24-q25 was associated with the transcriptional upregulation of this network of 54 mitotic apparatus genes, and we identify transcription factors that localize to these regions and putatively regulate mitotic activity. Knockdown of the mitotic network by siRNA identified 22 genes that might be considered as additional therapeutic targets for this clinically relevant patient subgroup. CONCLUSIONS: We define a molecular signature which may guide therapeutic approaches for tumors with high mitotic network activity

Kinome rewiring reveals AURKA limits PI3K-pathway inhibitor efficacy in breast cancer.

Dysregulation of the PI3K-AKT-mTOR signaling network is a prominent feature of breast cancers. However, clinical responses to drugs targeting this pathway have been modest, possibly because of dynamic changes in cellular signaling that drive resistance and limit drug efficacy. Using a quantitative chemoproteomics approach, we mapped kinome dynamics in response to inhibitors of this pathway and identified signaling changes that correlate with drug sensitivity. Maintenance of AURKA after drug treatment was associated with resistance in breast cancer models. Incomplete inhibition of AURKA was a common source of therapy failure, and combinations of PI3K, AKT or mTOR inhibitors with the AURKA inhibitor MLN8237 were highly synergistic and durably suppressed mTOR signaling, resulting in apoptosis and tumor regression in vivo. This signaling map identifies survival factors whose presence limits the efficacy of targeted therapies and reveals new drug combinations that may unlock the full potential of PI3K-AKT-mTOR pathway inhibitors in breast cancer

The expression level of HJURP has an independent prognostic impact and predicts the sensitivity to radiotherapy in breast cancer

INTRODUCTION. HJURP (Holliday Junction Recognition Protein) is a newly discovered gene reported to function at centromeres and to interact with CENPA. However its role in tumor development remains largely unknown. The goal of this study was to investigate the clinical significance of HJURP in breast cancer and its correlation with radiotherapeutic outcome. METHODS. We measured HJURP expression level in human breast cancer cell lines and primary breast cancers by Western blot and/or by Affymetrix Microarray; and determined its associations with clinical variables using standard statistical methods. Validation was performed with the use of published microarray data. We assessed cell growth and apoptosis of breast cancer cells after radiation using high-content image analysis. RESULTS. HJURP was expressed at higher level in breast cancer than in normal breast tissue. HJURP mRNA levels were significantly associated with estrogen receptor (ER), progesterone receptor (PR), Scarff-Bloom-Richardson (SBR) grade, age and Ki67 proliferation indices, but not with pathologic stage, ERBB2, tumor size, or lymph node status. Higher HJURP mRNA levels significantly decreased disease-free and overall survival. HJURP mRNA levels predicted the prognosis better than Ki67 proliferation indices. In a multivariate Cox proportional-hazard regression, including clinical variables as covariates, HJURP mRNA levels remained an independent prognostic factor for disease-free and overall survival. In addition HJURP mRNA levels were an independent prognostic factor over molecular subtypes (normal like, luminal, Erbb2 and basal). Poor clinical outcomes among patients with high HJURP expression were validated in five additional breast cancer cohorts. Furthermore, the patients with high HJURP levels were much more sensitive to radiotherapy. In vitro studies in breast cancer cell lines showed that cells with high HJURP levels were more sensitive to radiation treatment and had a higher rate of apoptosis than those with low levels. Knock down of HJURP in human breast cancer cells using shRNA reduced the sensitivity to radiation treatment. HJURP mRNA levels were significantly correlated with CENPA mRNA levels. CONCLUSIONS. HJURP mRNA level is a prognostic factor for disease-free and overall survival in patients with breast cancer and is a predictive biomarker for sensitivity to radiotherapy.National Institutes of Health, National Cancer Institute (R01 CA116481, P50 CA 5820, P30 CA 82103, U54 CA 112970); Office of Science; U.S. Department of Energy Office of Science, Office of Biological & Environmental Research (DE-AC02-05CH11231

Basal-Subtype and MEK-Pl3K Feedback Signaling Determine Susceptibility of Breast Cancer Cells to MEK Inhibition

Specific inhibitors of MEK have been developed that efficiently inhibit the oncogenic RAF-MEK-ERK pathway. We employed a systems-based approach to identify breast cancer subtypes particularly susceptible to MEK inhibitors and to understand molecular mechanisms conferring resistance to such compounds. Basal-type breast cancer cells were found to be particularly susceptible to growth-inhibition by small-molecule MEK inhibitors. Activation of the PI3 kinase pathway in response to MEK inhibition through a negative MEK-EGFR-PI3 kinase feedback loop was found to limit efficacy. Interruption of this feedback mechanism by targeting MEK and PI3 kinase produced synergistic effects, including induction of apoptosis and, in some cell lines, cell cycle arrest and protection from apoptosis induced by proapoptotic agents. These findings enhance our understanding of the interconnectivity of oncogenic signal transduction circuits and have implications for the design of future clinical trials of MEK inhibitors in breast cancer by guiding patient selection and suggesting rational combination therapies

A systems analysis of the chemosensitivity of breast cancer cells to the polyamine analogue PG-11047

<p>Abstract</p> <p>Background</p> <p>Polyamines regulate important cellular functions and polyamine dysregulation frequently occurs in cancer. The objective of this study was to use a systems approach to study the relative effects of PG-11047, a polyamine analogue, across breast cancer cells derived from different patients and to identify genetic markers associated with differential cytotoxicity.</p> <p>Methods</p> <p>A panel of 48 breast cell lines that mirror many transcriptional and genomic features present in primary human breast tumours were used to study the antiproliferative activity of PG-11047. Sensitive cell lines were further examined for cell cycle distribution and apoptotic response. Cell line responses, quantified by the GI₅₀(dose required for 50% relative growth inhibition) were correlated with the omic profiles of the cell lines to identify markers that predict response and cellular functions associated with drug sensitivity.</p> <p>Results</p> <p>The concentrations of PG-11047 needed to inhibit growth of members of the panel of breast cell lines varied over a wide range, with basal-like cell lines being inhibited at lower concentrations than the luminal cell lines. Sensitive cell lines showed a significant decrease in S phase fraction at doses that produced little apoptosis. Correlation of the GI₅₀values with the omic profiles of the cell lines identified genomic, transcriptional and proteomic variables associated with response.</p> <p>Conclusions</p> <p>A 13-gene transcriptional marker set was developed as a predictor of response to PG-11047 that warrants clinical evaluation. Analyses of the pathways, networks and genes associated with response to PG-11047 suggest that response may be influenced by interferon signalling and differential inhibition of aspects of motility and epithelial to mesenchymal transition.</p> <p>See the related commentary by Benes and Settleman: <url>http://www.biomedcentral.com/1741-7015/7/78</url></p

Enzalutamide response in a panel of prostate cancer cell lines reveals a role for glucocorticoid receptor in enzalutamide resistant disease.

Representative in vitro model systems that accurately model response to therapy and allow the identification of new targets are important for improving our treatment of prostate cancer. Here we describe molecular characterization and drug testing in a panel of 20 prostate cancer cell lines. The cell lines cluster into distinct subsets based on RNA expression, which is largely driven by functional Androgen Receptor (AR) expression. KLK3, the AR-responsive gene that encodes prostate specific antigen, shows the greatest variability in expression across the cell line panel. Other common prostate cancer associated genes such as TMPRSS2 and ERG show similar expression patterns. Copy number analysis demonstrates that many of the most commonly gained (including regions containing TERC and MYC) and lost regions (including regions containing TP53 and PTEN) that were identified in patient samples by the TCGA are mirrored in the prostate cancer cell lines. Assessment of response to the anti-androgen enzalutamide shows a distinct separation of responders and non-responders, predominantly related to status of wild-type AR. Surprisingly, several AR-null lines responded to enzalutamide. These AR-null, enzalutamide-responsive cells were characterized by high levels of expression of glucocorticoid receptor (GR) encoded by NR3C1. Treatment of these cells with the anti-GR agent mifepristone showed that they were more sensitive to this drug than enzalutamide, as were several of the enzalutamide non-responsive lines. This is consistent with several recent reports that suggest that GR expression is an alternative signaling mechanism that can bypass AR blockade. This study reinforces the utility of large cell line panels for the study of cancer and identifies several cell lines that represent ideal models to study AR-null cells that have upregulated GR to sustain growth

Targeting KRAS4A splicing through the RBM39/DCAF15 pathway inhibits cancer stem cells.

The commonly mutated human KRAS oncogene encodes two distinct KRAS4A and KRAS4B proteins generated by differential splicing. We demonstrate here that coordinated regulation of both isoforms through control of splicing is essential for development of Kras mutant tumors. The minor KRAS4A isoform is enriched in cancer stem-like cells, where it responds to hypoxia, while the major KRAS4B is induced by ER stress. KRAS4A splicing is controlled by the DCAF15/RBM39 pathway, and deletion of KRAS4A or pharmacological inhibition of RBM39 using Indisulam leads to inhibition of cancer stem cells. Our data identify existing clinical drugs that target KRAS4A splicing, and suggest that levels of the minor KRAS4A isoform in human tumors can be a biomarker of sensitivity to some existing cancer therapeutics