Sam Ng

Whole Genome Analysis Informs Breast Cancer Response to Aromatase Inhibition

To correlate the variable clinical features of oestrogen-receptor-positive breast cancer with somatic alterations, we studied pretreatment tumour biopsies accrued from patients in two studies of neoadjuvant aromatase inhibitor therapy by massively parallel sequencing and analysis. Eighteen significantly mutated genes were identified, including five genes (RUNX1, CBFB, MYH9, MLL3 and SF3B1) previously linked to haematopoietic disorders. Mutant MAP3K1 was associated with luminal A status, low-grade histology and low proliferation rates, whereas mutant TP53 was associated with the opposite pattern. Moreover, mutant GATA3 correlated with suppression of proliferation upon aromatase inhibitor treatment. Pathway analysis demonstrated that mutations in MAP2K4, a MAP3K1 substrate, produced similar perturbations as MAP3K1 loss. Distinct phenotypes in oestrogen-receptor-positive breast cancer are associated with specific patterns of somatic mutations that map into cellular pathways linked to tumour biology, but most recurrent mutations are relatively infrequent. Prospective clinical trials based on these findings will require comprehensive genome sequencing.

Subtype and pathway specific responses to anticancer compounds in breast cancer

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.

PARADIGM-SHIFT predicts the function of mutations in multiple cancers using pathway impact analysis

Motivation: A current challenge in understanding cancer processes is to pinpoint which mutations influence the onset and progression of disease. Toward this goal, we describe a method called PARADIGM-SHIFT that can predict whether a mutational event is neutral, gain-or loss-of-function in a tumor sample. The method uses a belief-propagation algorithm to infer gene activity from gene expression and copy number data in the context of a set of pathway interactions. Results: The method was found to be both sensitive and specific on a set of positive and negative controls for multiple cancers for which pathway information was available. Application to the Cancer Genome Atlas glioblastoma, ovarian and lung squamous cancer datasets revealed several novel mutations with predicted high impact including several genes mutated at low frequency suggesting the approach will be complementary to current approaches that rely on the prevalence of events to reach statistical significance. Availability: All source code is available at the github repository http:github.org/paradigmshift. Contact: jstuart@soe.ucsc.edu Supplementary information: Supplementary data are available at Bioinformatics online.

The UCSC Interaction Browser: multidimensional data views in pathway context.

High-throughput data sets such as genome-wide protein–protein interactions, protein–DNA interactions and gene expression data have been published for several model systems, especially for human cancer samples. The University of California, Santa Cruz (UCSC) Interaction Browser (http://sysbio.soe.ucsc.edu/nets) is an online tool for biologists to view high-throughput data sets simultaneously for the analysis of functional relationships between biological entities. Users can access several public interaction networks and functional genomics data sets through the portal as well as upload their own networks and data sets for analysis. Users can navigate through correlative relationships for focused sets of genes belonging to biological pathways using a standard web browser. Using a new visual modality called the CircleMap, multiple ‘omics’ data sets can be viewed simultaneously within the context of curated, predicted, directed and undirected regulatory interactions. The Interaction Browser provides an integrative viewing of biological networks based on the consensus of many observations about genes and their products, which may provide new insights about normal and disease processes not obvious from any isolated data set.

TNF-α modulates genome-wide redistribution of ΔNp63α/TAp73 and NF-κB cREL interactive binding on TP53 and AP-1 motifs to promote an oncogenic gene program in squamous cancer.

The Cancer Genome Atlas (TCGA) network study of 12 cancer types (PanCancer 12) revealed frequent mutation of TP53, and amplification and expression of related TP63 isoform ΔNp63 in squamous cancers. Further, aberrant expression of inflammatory genes and TP53/p63/p73 targets were detected in the PanCancer 12 project, reminiscent of gene programs comodulated by cREL/ΔNp63/TAp73 transcription factors we uncovered in head and neck squamous cell carcinomas (HNSCCs). However, how inflammatory gene signatures and cREL/p63/p73 targets are comodulated genome wide is unclear. Here, we examined how the inflammatory factor tumor necrosis factor-α (TNF-α) broadly modulates redistribution of cREL with ΔNp63α/TAp73 complexes and signatures genome wide in the HNSCC model UM-SCC46 using chromatin immunoprecipitation sequencing (ChIP-seq). TNF-α enhanced genome-wide co-occupancy of cREL with ΔNp63α on TP53/p63 sites, while unexpectedly promoting redistribution of TAp73 from TP53 to activator protein-1 (AP-1) sites. cREL, ΔNp63α and TAp73 binding and oligomerization on NF-κB-, TP53- or AP-1-specific sequences were independently validated by ChIP-qPCR (quantitative PCR), oligonucleotide-binding assays and analytical ultracentrifugation. Function of the binding activity was confirmed using TP53-, AP-1- and NF-κB-specific REs or p21, SERPINE1 and IL-6 promoter luciferase reporter activities. Concurrently, TNF-α regulated a broad gene network with cobinding activities for cREL, ΔNp63α and TAp73 observed upon array profiling and reverse transcription–PCR. Overlapping target gene signatures were observed in squamous cancer subsets and in inflamed skin of transgenic mice overexpressing ΔNp63α. Furthermore, multiple target genes identified in this study were linked to TP63 and TP73 activity and increased gene expression in large squamous cancer samples from PanCancer 12 TCGA by CircleMap. PARADIGM inferred pathway analysis revealed the network connection of TP63 and NF-κB complexes through an AP-1 hub, further supporting our findings. Thus, inflammatory cytokine TNF-α mediates genome-wide redistribution of the cREL/p63/p73, and AP-1 interactome, to diminish TAp73 tumor suppressor function and reciprocally activate NF-κB and AP-1 gene programs implicated in malignancy.

Abstract 4165: Differential pathway activation associated with domain-specific PIK3CA mutations

A recent study of the mutation landscape of >3000 cancers across 12 major cancer types from the Cancer Genome Atlas (TCGA) program revealed PIK3CA as the second most commonly mutated gene, occurring at >10% frequency in 8 types of cancer. Limited preclinical evidence suggested that mutations affecting PIK3CA catalytic vs. non-catalytic domains can produce different phenotypic consequences; however, whether domain-specific PIK3CA mutations results in distinct pathway consequences across multiple cancer types remain unclear. Thus, we used the PARADIGM algorithm, which integrates gene expression and copy number data into a superimposed pathway structure, to infer the activities of ∼13K pathway features and compared the signaling consequences associated with different domain-specific PIK3CA mutations within the TCGA Pan-Cancer dataset. Restricting to tumors harboring missense mutations in the coding region of a single PIK3CA domain resulted in 447 unique cases. PIK3CA mutations are distributed across the domains as follows: adaptor binding domain (ABD) = 23, Ras-binding domain (RBD) = 1, C2 = 50, helical = 199, and kinase = 174. Interestingly, the distribution of PIK3CA mutations among the domains is significantly different across cancer types (chi-square test p 50% of breast cancer PIK3CA mutations, while mutations in the helical domain predominate in head-and-neck and lung squamous carcinomas. Employing logistic regression and adjusting for cancer type, we identified 711 pathway features associated with kinase domain mutations (p Citation Format: Christina Yau, Stephen Benz, Charles Vaske, Sam Ng, Josh Stuart, Christopher C. Benz. Differential pathway activation associated with domain-specific PIK3CA mutations. [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 4165. doi:10.1158/1538-7445.AM2014-4165

TNF-α modulates genome-wide redistribution ofΔNp63α/TAp73 and NF-κB c-REL interactive binding on TP53and AP-1 motifs to promote an oncogenic gene program in squamouscancer

The Cancer Genome Atlas (TCGA) network study of 12 cancer types (PanCancer 12) revealed frequent mutation of TP53, and amplification and expression of related TP63 isoform ΔNp63 in squamous cancers. Further, aberrant expression of inflammatory genes and TP53/p63/p73 targets were detected in the PanCancer 12 project, reminiscent of gene programs comodulated by cREL/ΔNp63/TAp73 transcription factors we uncovered in head and neck squamous cell carcinomas (HNSCCs). However, how inflammatory gene signatures and cREL/p63/p73 targets are comodulated genome wide is unclear. Here, we examined how the inflammatory factor tumor necrosis factor-α (TNF-α) broadly modulates redistribution of cREL with ΔNp63α/TAp73 complexes and signatures genome wide in the HNSCC model UM-SCC46 using chromatin immunoprecipitation sequencing (ChIP-seq). TNF-α enhanced genome-wide co-occupancy of cREL with ΔNp63α on TP53/p63 sites, while unexpectedly promoting redistribution of TAp73 from TP53 to activator protein-1 (AP-1) sites. cREL, ΔNp63α and TAp73 binding and oligomerization on NF-κB-, TP53- or AP-1-specific sequences were independently validated by ChIP-qPCR (quantitative PCR), oligonucleotide-binding assays and analytical ultracentrifugation. Function of the binding activity was confirmed using TP53-, AP-1- and NF-κB-specific REs or p21, SERPINE1 and IL-6 promoter luciferase reporter activities. Concurrently, TNF-α regulated a broad gene network with cobinding activities for cREL, ΔNp63α and TAp73 observed upon array profiling and reverse transcription–PCR. Overlapping target gene signatures were observed in squamous cancer subsets and in inflamed skin of transgenic mice overexpressing ΔNp63α. Furthermore, multiple target genes identified in this study were linked to TP63 and TP73 activity and increased gene expression in large squamous cancer samples from PanCancer 12 TCGA by CircleMap. PARADIGM inferred pathway analysis revealed the network connection of TP63 and NF-κB complexes through an AP-1 hub, further supporting our findings. Thus, inflammatory cytokine TNF-α mediates genome-wide redistribution of the cREL/p63/p73, and AP-1 interactome, to diminish TAp73 tumor suppressor function and reciprocally activate NF-κB and AP-1 gene programs implicated in malignancy.

TNF-|[alpha]| modulates genome-wide redistribution of |[Delta]|Np63|[alpha]||[sol]|TAp73 and NF-|[kappa]|B cREL interactive binding on TP53 and AP-1 motifs to promote an oncogenic gene program in squamous cancer

The Cancer Genome Atlas (TCGA) network study of 12 cancer types (PanCancer 12) revealed frequent mutation of TP53, and amplification and expression of related TP63 isoform ΔNp63 in squamous cancers. Further, aberrant expression of inflammatory genes and TP53/p63/p73 targets were detected in the PanCancer 12 project, reminiscent of gene programs comodulated by cREL/ΔNp63/TAp73 transcription factors we uncovered in head and neck squamous cell carcinomas (HNSCCs). However, how inflammatory gene signatures and cREL/p63/p73 targets are comodulated genome wide is unclear. Here, we examined how the inflammatory factor tumor necrosis factor-α (TNF-α) broadly modulates redistribution of cREL with ΔNp63α/TAp73 complexes and signatures genome wide in the HNSCC model UM-SCC46 using chromatin immunoprecipitation sequencing (ChIP-seq). TNF-α enhanced genome-wide co-occupancy of cREL with ΔNp63α on TP53/p63 sites, while unexpectedly promoting redistribution of TAp73 from TP53 to activator protein-1 (AP-1) sites. cREL, ΔNp63α and TAp73 binding and oligomerization on NF-κB-, TP53- or AP-1-specific sequences were independently validated by ChIP-qPCR (quantitative PCR), oligonucleotide-binding assays and analytical ultracentrifugation. Function of the binding activity was confirmed using TP53-, AP-1- and NF-κB-specific REs or p21, SERPINE1 and IL-6 promoter luciferase reporter activities. Concurrently, TNF-α regulated a broad gene network with cobinding activities for cREL, ΔNp63α and TAp73 observed upon array profiling and reverse transcription–PCR. Overlapping target gene signatures were observed in squamous cancer subsets and in inflamed skin of transgenic mice overexpressing ΔNp63α. Furthermore, multiple target genes identified in this study were linked to TP63 and TP73 activity and increased gene expression in large squamous cancer samples from PanCancer 12 TCGA by CircleMap. PARADIGM inferred pathway analysis revealed the network connection of TP63 and NF-κB complexes through an AP-1 hub, further supporting our findings. Thus, inflammatory cytokine TNF-α mediates genome-wide redistribution of the cREL/p63/p73, and AP-1 interactome, to diminish TAp73 tumor suppressor function and reciprocally activate NF-κB and AP-1 gene programs implicated in malignancy.

Abstract PR02: Multiple Pathway Learning accurately predicts gene essentiality in the Cancer Cell Line Encyclopedia

We applied biologically-motivated feature transformations coupled with established machine learning methods to predict gene essentiality in CCLE cell line models. By leveraging additional large datasets, such as The Cancer Genome Atlas PanCancer12 data and MSigDB pathway definitions, we improved the robustness and biological interpretability of our models. We developed a multi-pathway learning (MPL) approach that associates a genetic pathway from MSigDB with a distinct kernel for use in a multiple kernel learning setting. We evaluated the performance of MPL compared to several other regression methods including random forests, kernel ridge regression, and elastic net linear models, We combined multiple approaches using an ensemble technique on the diverse set of predictors. We found that the best performing method was an ensemble combining MPL and random forest predictions. Both models utilized features derived from both gene expression and copy number data, the latter of which were filtered to those predicted as driver events in prior pan-cancer studies. The ensemble method was a joint winner in the recent DREAM 9 gene essentiality prediction challenge. MPL also demonstrated merit as a feature selector when used with other downstream methods. The ensemble performed best at predicting the essentiality of genes involved in cell cycle control (cyclins and cyclin-dependent kinases), protein degradation (proteasome complex), cell proliferation signaling (sonic hedgehog, Aurora-B, RAC1), apoptosis (RB1,TP53) and hypoxia response (VEGF, VHL). Many of the key genes in those pathways are known to be drivers of cancer progression, suggesting our method9s utility as a biomarker for detecting key tumorigenic events. The advantage of MPL is that mechanistically coherent gene sets are automatically selected as high scoring pathway kernels (HSPKs). We investigated whether the HSPKs identify cellular processes relevant to the loss of key genes. To do this, we inspected the HSPKs for a few of the most abundantly mutated genes in cancer. The MPL predictor for TP53 included the targets of this transcription factor as well as HSPKs involved in apoptosis, a cellular process regulated by TP53. The retinoblastoma gene (RB1) MPL predictor included RB1 targets as well as HSPKs involved in the regulation of histone deacetylase (HDAC) that interacts with RB1 to suppress DNA synthesis. These findings suggest trends in the MPL results could reveal a pathway-level view of the synthetic lethal architecture of cells. Such a map, that links patterns of pathway expression to potential genetic vulnerabilities, could provide an invaluable tool for exploring new avenues to target cancer cells. Citation Format: Vladislav Uzunangelov, Evan Paull, Sahil Chopra, Daniel Carlin, Adrian Bivol, Kyle Ellrott, Kiley Graim, Yulia Newton, Sam Ng, Artem Sokolov, Joshua Stuart. Multiple Pathway Learning accurately predicts gene essentiality in the Cancer Cell Line Encyclopedia. [abstract]. In: Proceedings of the AACR Special Conference on Computational and Systems Biology of Cancer; Feb 8-11 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 2):Abstract nr PR02.

Abstract B1-38: PARADIGM-SHIFT: Predicts the functional impact of ‘driver modules’ in multiple cancers using pathway impact analysis

The major mechanism by which cancer arises is through genomic alterations. These alterations can lead to changes in gene regulation, protein structure, and function. Individual tumors can contain hundreds to thousands of alterations. It is critical to distinguish alterations that have an important role defining the cancer – drivers – from alterations that are unimportant to the tumor – passengers. Driver genes can lead to significant changes in their pathways; however, alterations in a single gene may not be sufficient to explain all the pathway perturbations across patients. Additional alterations could range from DNA copy number changes, gene-fusions, or even lesser understood non-coding mutations. Identifying these ‘driver modules’ is essential for understanding cancer disease mechanisms, which can help guide treatment decisions as well as identify novel targets for treatment. We have developed a functional impact prediction method called PARADIGM-SHIFT based on integrated pathway analysis to discriminate loss-of-function, neutral, and gain-of-function alterations. Utilizing the set of regulatory interactions annotated for a given gene, we can detect a shift in the downstream effects of an altered gene compared to what is expected from its upstream influences. Additionally, since these shifts in pathway signal can be detected for all samples, PARADIGM-SHIFT can be used to identify additional genomic alterations that lead to similar changes to the altered pathway to form ‘driver modules.’ Application of our method to the TCGA Pan-Cancer cohort identifies many genes with significant alterations that lead to loss- and gain- of function. PARADIGM-SHIFT then identifies several additional genomic alterations, including non-coding mutations variants, which are significantly implicated in these pathway changes. This analysis offers insight into the mechanism of non-coding mutations orthogonal to sequence-based methods by interpreting the pathway consequences of these variants. Citation Format: Sam Ng, Christopher Benz, David Haussler, Joshua M. Stuart. PARADIGM-SHIFT: Predicts the functional impact of ‘driver modules’ in multiple cancers using pathway impact analysis. [abstract]. In: Proceedings of the AACR Special Conference on Computational and Systems Biology of Cancer; Feb 8-11 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 2):Abstract nr B1-38.