Hydra: A mixture modeling framework for subtyping pediatric cancer cohorts using multimodal gene expression signatures.

Precision oncology has primarily relied on coding mutations as biomarkers of response to therapies. While transcriptome analysis can provide valuable information, incorporation into workflows has been difficult. For example, the relative rather than absolute gene expression level needs to be considered, requiring differential expression analysis across samples. However, expression programs related to the cell-of-origin and tumor microenvironment effects confound the search for cancer-specific expression changes. To address these challenges, we developed an unsupervised clustering approach for discovering differential pathway expression within cancer cohorts using gene expression measurements. The hydra approach uses a Dirichlet process mixture model to automatically detect multimodally distributed genes and expression signatures without the need for matched normal tissue. We demonstrate that the hydra approach is more sensitive than widely-used gene set enrichment approaches for detecting multimodal expression signatures. Application of the hydra analysis framework to small blue round cell tumors (including rhabdomyosarcoma, synovial sarcoma, neuroblastoma, Ewing sarcoma, and osteosarcoma) identified expression signatures associated with changes in the tumor microenvironment. The hydra approach also identified an association between ATRX deletions and elevated immune marker expression in high-risk neuroblastoma. Notably, hydra analysis of all small blue round cell tumors revealed similar subtypes, characterized by changes to infiltrating immune and stromal expression signatures

Chronic cisplatin treatment promotes enhanced damage repair and tumor progression in a mouse model of lung cancer

Chemotherapy resistance is a major obstacle in cancer treatment, yet the mechanisms of response to specific therapies have been largely unexplored in vivo. Employing genetic, genomic, and imaging approaches, we examined the dynamics of response to a mainstay chemotherapeutic, cisplatin, in multiple mouse models of human non-small-cell lung cancer (NSCLC). We show that lung tumors initially respond to cisplatin by sensing DNA damage, undergoing cell cycle arrest, and inducing apoptosis—leading to a significant reduction in tumor burden. Importantly, we demonstrate that this response does not depend on the tumor suppressor p53 or its transcriptional target, p21. Prolonged cisplatin treatment promotes the emergence of resistant tumors with enhanced repair capacity that are cross-resistant to platinum analogs, exhibit advanced histopathology, and possess an increased frequency of genomic alterations. Cisplatin-resistant tumors express elevated levels of multiple DNA damage repair and cell cycle arrest-related genes, including p53-inducible protein with a death domain (Pidd). We demonstrate a novel role for PIDD as a regulator of chemotherapy response in human lung tumor cells.National Institutes of Health (U.S.) (grant 5-UO1-CA84306)National Cancer Institute (U.S.) (CA034992

Multiple Oncogenic Pathway Signatures Show Coordinate Expression Patterns in Human Prostate Tumors

BACKGROUND: Gene transcription patterns associated with activation of oncogenes Myc, c-Src, beta-catenin, E2F3, H-Ras, HER2, EGFR, MEK, Raf, MAPK, Akt, and cyclin D1, as well as of the cell cycle and of androgen signaling have been generated in previous studies using experimental models. It was not clear whether genes in these "oncogenic signatures" would show coordinate expression patterns in human prostate tumors, particularly as most of the signatures were derived from cell types other than prostate. PRINCIPAL FINDINGS: The above oncogenic pathway signatures were examined in four different gene expression profile datasets of human prostate tumors (representing approximately 250 patients in all), using both Q1-Q2 and one-sided Fisher's exact enrichment analysis methods. A significant fraction (approximately 5%) of genes up-regulated experimentally by Myc, c-Src, HER2, Akt, or androgen were co-expressed in human tumors with the oncogene or biomarker corresponding to the pathway signature. Genes down-regulated experimentally, however, did not show anticipated patterns of anti-enrichment in the human tumors. CONCLUSIONS: Significant subsets of the genes in these experimentally-derived oncogenic signatures are relevant to the study of human prostate cancer. Both molecular biologists and clinical researchers could focus attention on the relatively small number of genes identified here as having coordinate patterns that arise from both the experimental system and the human disease system

Rank–rank hypergeometric overlap: identification of statistically significant overlap between gene-expression signatures

Comparing independent high-throughput gene-expression experiments can generate hypotheses about which gene-expression programs are shared between particular biological processes. Current techniques to compare expression profiles typically involve choosing a fixed differential expression threshold to summarize results, potentially reducing sensitivity to small but concordant changes. We present a threshold-free algorithm called Rank–rank Hypergeometric Overlap (RRHO). This algorithm steps through two gene lists ranked by the degree of differential expression observed in two profiling experiments, successively measuring the statistical significance of the number of overlapping genes. The output is a graphical map that shows the strength, pattern and bounds of correlation between two expression profiles. To demonstrate RRHO sensitivity and dynamic range, we identified shared expression networks in cancer microarray profiles driving tumor progression, stem cell properties and response to targeted kinase inhibition. We demonstrate how RRHO can be used to determine which model system or drug treatment best reflects a particular biological or disease response. The threshold-free and graphical aspects of RRHO complement other rank-based approaches such as Gene Set Enrichment Analysis (GSEA), for which RRHO is a 2D analog. Rank–rank overlap analysis is a sensitive, robust and web-accessible method for detecting and visualizing overlap trends between two complete, continuous gene-expression profiles. A web-based implementation of RRHO can be accessed at http://systems.crump.ucla.edu/rankrank/

Combined inhibition of BET family proteins and histone deacetylases as a potential epigenetics-based therapy for pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal human cancers and shows resistance to any therapeutic strategy used. Here we tested small-molecule inhibitors targeting chromatin regulators as possible therapeutic agents in PDAC. We show that JQ1, an inhibitor of the bromodomain and extraterminal (BET) family of proteins, suppresses PDAC development in mice by inhibiting both MYC activity and inflammatory signals. The histone deacetylase (HDAC) inhibitor SAHA synergizes with JQ1 to augment cell death and more potently suppress advanced PDAC. Finally, using a CRISPR-Cas9–based method for gene editing directly in the mouse adult pancreas, we show that de-repression of p57 (also known as KIP2 or CDKN1C) upon combined BET and HDAC inhibition is required for the induction of combination therapy–induced cell death in PDAC. SAHA is approved for human use, and molecules similar to JQ1 are being tested in clinical trials. Thus, these studies identify a promising epigenetic-based therapeutic strategy that may be rapidly implemented in fatal human tumors

Disruption of STAT3 signaling promotes KRAS induced lung tumorigenesis

STAT3 is considered to play an oncogenic role in several malignancies including lung cancer; consequently, targeting STAT3 is currently proposed as therapeutic intervention. Here we demonstrate that STAT3 plays an unexpected tumour-suppressive role in KRAS mutant lung adenocarcinoma (AC). Indeed, lung tissue-specific inactivation of Stat3 in mice results in increased KrasG12D-driven AC initiation and malignant progression leading to markedly reduced survival. Knockdown of STAT3 in xenografted human AC cells increases tumour growth. Clinically, low STAT3 expression levels correlate with poor survival and advanced malignancy in human lung AC patients with smoking history, which are prone to KRAS mutations. Consistently, KRAS mutant lung tumours exhibit reduced STAT3 levels. Mechanistically, we demonstrate that STAT3 controls NF-B-induced IL-8 expression by sequestering NF-B within the cytoplasm, thereby inhibiting IL-8-mediated myeloid tumour infiltration and tumour vascularization and hence tumour progression. These results elucidate a novel STAT3NF-BIL-8 axis in KRAS mutant AC with therapeutic and prognostic relevance.P 25599(VLID)183891

Babelomics: an integrative platform for the analysis of transcriptomics, proteomics and genomic data with advanced functional profiling

Babelomics is a response to the growing necessity of integrating and analyzing different types of genomic data in an environment that allows an easy functional interpretation of the results. Babelomics includes a complete suite of methods for the analysis of gene expression data that include normalization (covering most commercial platforms), pre-processing, differential gene expression (case-controls, multiclass, survival or continuous values), predictors, clustering; large-scale genotyping assays (case controls and TDTs, and allows population stratification analysis and correction). All these genomic data analysis facilities are integrated and connected to multiple options for the functional interpretation of the experiments. Different methods of functional enrichment or gene set enrichment can be used to understand the functional basis of the experiment analyzed. Many sources of biological information, which include functional (GO, KEGG, Biocarta, Reactome, etc.), regulatory (Transfac, Jaspar, ORegAnno, miRNAs, etc.), text-mining or protein–protein interaction modules can be used for this purpose. Finally a tool for the de novo functional annotation of sequences has been included in the system. This provides support for the functional analysis of non-model species. Mirrors of Babelomics or command line execution of their individual components are now possible. Babelomics is available at http://www.babelomics.org

SLEPR: A Sample-Level Enrichment-Based Pathway Ranking Method — Seeking Biological Themes through Pathway-Level Consistency

Analysis of microarray and other high throughput data often involves identification of genes consistently up or down-regulated across samples as the first step in extraction of biological meaning. This gene-level paradigm can be limited as a result of valid sample fluctuations and biological complexities. In this report, we describe a novel method, SLEPR, which eliminates this limitation by relying on pathway-level consistencies. Our method first selects the sample-level differentiated genes from each individual sample, capturing genes missed by other analysis methods, ascertains the enrichment levels of associated pathways from each of those lists, and then ranks annotated pathways based on the consistency of enrichment levels of individual samples from both sample classes. As a proof of concept, we have used this method to analyze three public microarray datasets with a direct comparison with the GSEA method, one of the most popular pathway-level analysis methods in the field. We found that our method was able to reproduce the earlier observations with significant improvements in depth of coverage for validated or expected biological themes, but also produced additional insights that make biological sense. This new method extends existing analyses approaches and facilitates integration of different types of HTP data

A Mouse Stromal Response to Tumor Invasion Predicts Prostate and Breast Cancer Patient Survival

Primary and metastatic tumor growth induces host tissue responses that are believed to support tumor progression. Understanding the molecular changes within the tumor microenvironment during tumor progression may therefore be relevant not only for discovering potential therapeutic targets, but also for identifying putative molecular signatures that may improve tumor classification and predict clinical outcome. To selectively address stromal gene expression changes during cancer progression, we performed cDNA microarray analysis of laser-microdissected stromal cells derived from prostate intraepithelial neoplasia (PIN) and invasive cancer in a multistage model of prostate carcinogenesis. Human orthologs of genes identified in the stromal reaction to tumor progression in this mouse model were observed to be expressed in several human cancers, and to cluster prostate and breast cancer patients into groups with statistically different clinical outcomes. Univariate Cox analysis showed that overexpression of these genes is associated with shorter survival and recurrence-free periods. Taken together, our observations provide evidence that the expression signature of the stromal response to tumor invasion in a mouse tumor model can be used to probe human cancer, and to provide a powerful prognostic indicator for some of the most frequent human malignancies