Pathway Commons 2019 Update: integration, analysis and exploration of pathway data.

Pathway Commons (https://www.pathwaycommons.org) is an integrated resource of publicly available information about biological pathways including biochemical reactions, assembly of biomolecular complexes, transport and catalysis events and physical interactions involving proteins, DNA, RNA, and small molecules (e.g. metabolites and drug compounds). Data is collected from multiple providers in standard formats, including the Biological Pathway Exchange (BioPAX) language and the Proteomics Standards Initiative Molecular Interactions format, and then integrated. Pathway Commons provides biologists with (i) tools to search this comprehensive resource, (ii) a download site offering integrated bulk sets of pathway data (e.g. tables of interactions and gene sets), (iii) reusable software libraries for working with pathway information in several programming languages (Java, R, Python and Javascript) and (iv) a web service for programmatically querying the entire dataset. Visualization of pathways is supported using the Systems Biological Graphical Notation (SBGN). Pathway Commons currently contains data from 22 databases with 4794 detailed human biochemical processes (i.e. pathways) and ∼2.3 million interactions. To enhance the usability of this large resource for end-users, we develop and maintain interactive web applications and training materials that enable pathway exploration and advanced analysis.

Cytokit: a single-cell analysis toolkit for high dimensional fluorescent microscopy imaging.

BACKGROUND: Multiplexed in-situ fluorescent imaging offers several advantages over single-cell assays that do not preserve the spatial characteristics of biological samples. This spatial information, in addition to morphological properties and extensive intracellular or surface marker profiling, comprise promising avenues for rapid advancements in the understanding of disease progression and diagnosis. As protocols for conducting such imaging experiments continue to improve, it is the intent of this study to provide and validate software for processing the large quantity of associated data in kind. RESULTS: Cytokit offers (i) an end-to-end, GPU-accelerated image processing pipeline; (ii) efficient input/output (I/O) strategies for operations specific to high dimensional microscopy; and (iii) an interactive user interface for cross filtering of spatial, graphical, expression, and morphological cell properties within the 100+ GB image datasets common to multiplexed immunofluorescence. Image processing operations supported in Cytokit are generally sourced from existing deep learning models or are at least in part adapted from open source packages to run in a single or multi-GPU environment. The efficacy of these operations is demonstrated through several imaging experiments that pair Cytokit results with those from an independent but comparable assay. A further validation also demonstrates that previously published results can be reproduced from a publicly available multiplexed image dataset. CONCLUSION: Cytokit is a collection of open source tools for quantifying and analyzing properties of individual cells in large fluorescent microscopy datasets that are often, but not necessarily, generated from multiplexed antibody labeling protocols over many fields of view or time periods. This project is best suited to bioinformaticians or other technical users that wish to analyze such data in a batch-oriented, high-throughput setting. All source code, documentation, and data generated for this article are available under the Apache License 2.0 at https://github.com/hammerlab/cytokit .

Contribution of systemic and somatic factors to clinical response and resistance to PD-L1 blockade in urothelial cancer: An exploratory multi-omic analysis.

BACKGROUND: Inhibition of programmed death-ligand 1 (PD-L1) with atezolizumab can induce durable clinical benefit (DCB) in patients with metastatic urothelial cancers, including complete remissions in patients with chemotherapy refractory disease. Although mutation load and PD-L1 immune cell (IC) staining have been associated with response, they lack sufficient sensitivity and specificity for clinical use. Thus, there is a need to evaluate the peripheral blood immune environment and to conduct detailed analyses of mutation load, predicted neoantigens, and immune cellular infiltration in tumors to enhance our understanding of the biologic underpinnings of response and resistance. METHODS AND FINDINGS: The goals of this study were to (1) evaluate the association of mutation load and predicted neoantigen load with therapeutic benefit and (2) determine whether intratumoral and peripheral blood T cell receptor (TCR) clonality inform clinical outcomes in urothelial carcinoma treated with atezolizumab. We hypothesized that an elevated mutation load in combination with T cell clonal dominance among intratumoral lymphocytes prior to treatment or among peripheral T cells after treatment would be associated with effective tumor control upon treatment with anti-PD-L1 therapy. We performed whole exome sequencing (WES), RNA sequencing (RNA-seq), and T cell receptor sequencing (TCR-seq) of pretreatment tumor samples as well as TCR-seq of matched, serially collected peripheral blood, collected before and after treatment with atezolizumab. These parameters were assessed for correlation with DCB (defined as progression-free survival [PFS] >6 months), PFS, and overall survival (OS), both alone and in the context of clinical and intratumoral parameters known to be predictive of survival in this disease state. Patients with DCB displayed a higher proportion of tumor-infiltrating T lymphocytes (TIL) (n = 24, Mann-Whitney p = 0.047). Pretreatment peripheral blood TCR clonality below the median was associated with improved PFS (n = 29, log-rank p = 0.048) and OS (n = 29, log-rank p = 0.011). Patients with DCB also demonstrated more substantial expansion of tumor-associated TCR clones in the peripheral blood 3 weeks after starting treatment (n = 22, Mann-Whitney p = 0.022). The combination of high pretreatment peripheral blood TCR clonality with elevated PD-L1 IC staining in tumor tissue was strongly associated with poor clinical outcomes (n = 10, hazard ratio (HR) (mean) = 89.88, HR (median) = 23.41, 95% CI [2.43, 506.94], p(HR > 1) = 0.0014). Marked variations in mutation loads were seen with different somatic variant calling methodologies, which, in turn, impacted associations with clinical outcomes. Missense mutation load, predicted neoantigen load, and expressed neoantigen load did not demonstrate significant association with DCB (n = 25, Mann-Whitney p = 0.22, n = 25, Mann-Whitney p = 0.55, and n = 25, Mann-Whitney p = 0.29, respectively). Instead, we found evidence of time-varying effects of somatic mutation load on PFS in this cohort (n = 25, p = 0.044). A limitation of our study is its small sample size (n = 29), a subset of the patients treated on IMvigor 210 (NCT02108652). Given the number of exploratory analyses performed, we intend for these results to be hypothesis-generating. CONCLUSIONS: These results demonstrate the complex nature of immune response to checkpoint blockade and the compelling need for greater interrogation and data integration of both host and tumor factors. Incorporating these variables in prospective studies will facilitate identification and treatment of resistant patients.

PaxtoolsR: pathway analysis in R using Pathway Commons.

PURPOSE: PaxtoolsR package enables access to pathway data represented in the BioPAX format and made available through the Pathway Commons webservice for users of the R language to aid in advanced pathway analyses. Features include the extraction, merging and validation of pathway data represented in the BioPAX format. This package also provides novel pathway datasets and advanced querying features for R users through the Pathway Commons webservice allowing users to query, extract and retrieve data and integrate these data with local BioPAX datasets. AVAILABILITY AND IMPLEMENTATION: The PaxtoolsR package is compatible with versions of R 3.1.1 (and higher) on Windows, Mac OS X and Linux using Bioconductor 3.0 and is available through the Bioconductor R package repository along with source code and a tutorial vignette describing common tasks, such as data visualization and gene set enrichment analysis. Source code and documentation are at http://www.bioconductor.org/packages/paxtoolsr This plugin is free, open-source and licensed under the LGPL-3. CONTACT: paxtools@cbio.mskcc.org or lunaa@cbio.mskcc.org.

Pattern search in BioPAX models.

MOTIVATION: BioPAX is a standard language for representing complex cellular processes, including metabolic networks, signal transduction and gene regulation. Owing to the inherent complexity of a BioPAX model, searching for a specific type of subnetwork can be non-trivial and difficult. RESULTS: We developed an open source and extensible framework for defining and searching graph patterns in BioPAX models. We demonstrate its use with a sample pattern that captures directed signaling relations between proteins. We provide search results for the pattern obtained from the Pathway Commons database and compare these results with the current data in signaling databases SPIKE and SignaLink. Results show that a pattern search in public pathway data can identify a substantial amount of signaling relations that do not exist in signaling databases. AVAILABILITY: BioPAX-pattern software was developed in Java. Source code and documentation is freely available at http://code.google.com/p/biopax-pattern under Lesser GNU Public License.

Prediction of individualized therapeutic vulnerabilities in cancer from genomic profiles.

MOTIVATION: Somatic homozygous deletions of chromosomal regions in cancer, while not necessarily oncogenic, may lead to therapeutic vulnerabilities specific to cancer cells compared with normal cells. A recently reported example is the loss of one of the two isoenzymes in glioblastoma cancer cells such that the use of a specific inhibitor selectively inhibited growth of the cancer cells, which had become fully dependent on the second isoenzyme. We have now made use of the unprecedented conjunction of large-scale cancer genomics profiling of tumor samples in The Cancer Genome Atlas (TCGA) and of tumor-derived cell lines in the Cancer Cell Line Encyclopedia, as well as the availability of integrated pathway information systems, such as Pathway Commons, to systematically search for a comprehensive set of such epistatic vulnerabilities. RESULTS: Based on homozygous deletions affecting metabolic enzymes in 16 TCGA cancer studies and 972 cancer cell lines, we identified 4104 candidate metabolic vulnerabilities present in 1019 tumor samples and 482 cell lines. Up to 44% of these vulnerabilities can be targeted with at least one Food and Drug Administration-approved drug. We suggest focused experiments to test these vulnerabilities and clinical trials based on personalized genomic profiles of those that pass preclinical filters. We conclude that genomic profiling will in the future provide a promising basis for network pharmacology of epistatic vulnerabilities as a promising therapeutic strategy. AVAILABILITY AND IMPLEMENTATION: A web-based tool for exploring all vulnerabilities and their details is available at http://cbio.mskcc.org/cancergenomics/statius/ along with supplemental data files.

Integrating biological pathways and genomic profiles with ChiBE 2.

BACKGROUND: Dynamic visual exploration of detailed pathway information can help researchers digest and interpret complex mechanisms and genomic datasets. RESULTS: ChiBE is a free, open-source software tool for visualizing, querying, and analyzing human biological pathways in BioPAX format. The recently released version 2 can search for neighborhoods, paths between molecules, and common regulators/targets of molecules, on large integrated cellular networks in the Pathway Commons database as well as in local BioPAX models. Resulting networks can be automatically laid out for visualization using a graphically rich, process-centric notation. Profiling data from the cBioPortal for Cancer Genomics and expression data from the Gene Expression Omnibus can be overlaid on these networks. CONCLUSIONS: ChiBE's new capabilities are organized around a genomics-oriented workflow and offer a unique comprehensive pathway analysis solution for genomics researchers. The software is freely available at http://code.google.com/p/chibe.

PiHelper: an open source framework for drug-target and antibody-target data.

MOTIVATION: The interaction between drugs and their targets, often proteins, and between antibodies and their targets, is important for planning and analyzing investigational and therapeutic interventions in many biological systems. Although drug-target and antibody-target datasets are available in separate databases, they are not publicly available in an integrated bioinformatics resource. As medical therapeutics, especially in cancer, increasingly uses targeted drugs and measures their effects on biomolecular profiles, there is an unmet need for a user-friendly toolset that allows researchers to comprehensively and conveniently access and query information about drugs, antibodies and their targets. SUMMARY: The PiHelper framework integrates human drug-target and antibody-target associations from publicly available resources to help meet the needs of researchers in systems pharmacology, perturbation biology and proteomics. PiHelper has utilities to (i) import drug- and antibody-target information; (ii) search the associations either programmatically or through a web user interface (UI); (iii) visualize the data interactively in a network; and (iv) export relationships for use in publications or other analysis tools. AVAILABILITY: PiHelper is a free software under the GNU Lesser General Public License (LGPL) v3.0. Source code and documentation are at http://bit.ly/pihelper. We plan to coordinate contributions from the community by managing future releases.

Using biological pathway data with paxtools.

A rapidly growing corpus of formal, computable pathway information can be used to answer important biological questions including finding non-trivial connections between cellular processes, identifying significantly altered portions of the cellular network in a disease state and building predictive models that can be used for precision medicine. Due to its complexity and fragmented nature, however, working with pathway data is still difficult. We present Paxtools, a Java library that contains algorithms, software components and converters for biological pathways represented in the standard BioPAX language. Paxtools allows scientists to focus on their scientific problem by removing technical barriers to access and analyse pathway information. Paxtools can run on any platform that has a Java Runtime Environment and was tested on most modern operating systems. Paxtools is open source and is available under the Lesser GNU public license (LGPL), which allows users to freely use the code in their software systems with a requirement for attribution. Source code for the current release (4.2.0) can be found in Software S1. A detailed manual for obtaining and using Paxtools can be found in Protocol S1. The latest sources and release bundles can be obtained from biopax.org/paxtools.

Progressively Enhancing Stemness of Adoptively Transferred T Cells with PI3Kδ Blockade Improves Metabolism and Antitumor Immunity.

Generating stem-like memory T cells (TSCM) is a potential strategy to improve adoptive immunotherapy. Elucidating optimal ways to modulate signaling pathways that enrich TSCM properties could identify approaches to achieve this goal. We discovered herein that blocking the PI3Kδ pathway pharmaceutically to varying degrees can generate T cells with increasingly heightened stemness properties, based on the progressive enrichment of the transcription factors Tcf1 and Lef1. T cells with enhanced stemness features exhibited metabolic plasticity, marked by improved mitochondrial function and glucose uptake after tumor recognition. Conversely, T cells with low or medium stemness were less metabolically dynamic, vulnerable to antigen-induced cell death, and expressed more inhibitory checkpoint receptors. Only T-cell receptor-specific or chimeric antigen receptor (CAR)-specific T cells with high stemness persisted in vivo and mounted protective immunity to tumors. Likewise, the strongest level of PI3Kδ blockade in vitro generated human tumor-infiltrating lymphocytes and CAR T cells with elevated stemness properties, in turn bolstering their capacity to regress human solid tumors. The stemness level of T cells in vitro was important, ultimately impacting their efficacy in mice bearing three distinct solid tumors. Lef1 and Tcf1 sustained antitumor protection by donor high CD8+ TSCM or CD4+ Th17SCM, as deletion of either one compromised the therapeutic efficacy. Collectively, these findings highlight the importance of strategic modulation of PI3Kδ signaling in T cells to induce stemness and lasting protective responses to solid tumors. SIGNIFICANCE: Elevating T-cell stemness by progressively blocking PI3Kδ signaling during ex vivo manufacturing of adoptive cell therapies alters metabolic and functional properties to enhance antitumor immunity dependent on Tcf1 and Lef1.

Cancer-associated mutations in DICER1 RNase IIIa and IIIb domains exert similar effects on miRNA biogenesis.

Somatic mutations in the RNase IIIb domain of DICER1 arise in cancer and disrupt the cleavage of 5' pre-miRNA arms. Here, we characterize an unstudied, recurrent, mutation (S1344L) in the DICER1 RNase IIIa domain in tumors from The Cancer Genome Atlas (TCGA) project and MSK-IMPACT profiling. RNase IIIa/b hotspots are absent from most cancers, but are notably enriched in uterine cancers. Systematic analysis of TCGA small RNA datasets show that DICER1 RNase IIIa-S1344L tumors deplete 5p-miRNAs, analogous to RNase IIIb hotspot samples. Structural and evolutionary coupling analyses reveal constrained proximity of RNase IIIa-S1344 to the RNase IIIb catalytic site, rationalizing why mutation of this site phenocopies known hotspot alterations. Finally, examination of DICER1 hotspot endometrial tumors reveals derepression of specific miRNA target signatures. In summary, comprehensive analyses of DICER1 somatic mutations and small RNA data reveal a mechanistic aspect of pre-miRNA processing that manifests in specific cancer settings.

T cells genetically engineered to overcome death signaling enhance adoptive cancer immunotherapy.

Across clinical trials, T cell expansion and persistence following adoptive cell transfer (ACT) have correlated with superior patient outcomes. Herein, we undertook a pan-cancer analysis to identify actionable ligand-receptor pairs capable of compromising T cell durability following ACT. We discovered that FASLG, the gene encoding the apoptosis-inducing ligand FasL, is overexpressed within the majority of human tumor microenvironments (TMEs). Further, we uncovered that Fas, the receptor for FasL, is highly expressed on patient-derived T cells used for clinical ACT. We hypothesized that a cognate Fas-FasL interaction within the TME might limit both T cell persistence and antitumor efficacy. We discovered that genetic engineering of Fas variants impaired in the ability to bind FADD functioned as dominant negative receptors (DNRs), preventing FasL-induced apoptosis in Fas-competent T cells. T cells coengineered with a Fas DNR and either a T cell receptor or chimeric antigen receptor exhibited enhanced persistence following ACT, resulting in superior antitumor efficacy against established solid and hematologic cancers. Despite increased longevity, Fas DNR-engineered T cells did not undergo aberrant expansion or mediate autoimmunity. Thus, T cell-intrinsic disruption of Fas signaling through genetic engineering represents a potentially universal strategy to enhance ACT efficacy across a broad range of human malignancies.

Anti-CTLA-4 Activates Intratumoral NK Cells and Combined with IL15/IL15Rα Complexes Enhances Tumor Control.

Antibodies targeting CTLA-4 induce durable responses in some patients with melanoma and are being tested in a variety of human cancers. However, these therapies are ineffective for a majority of patients across tumor types. Further understanding the immune alterations induced by these therapies may enable the development of novel strategies to enhance tumor control and biomarkers to identify patients most likely to respond. In several murine models, including colon26, MC38, CT26, and B16 tumors cotreated with GVAX, anti-CTLA-4 efficacy depends on interactions between the Fc region of CTLA-4 antibodies and Fc receptors (FcR). Anti-CTLA-4 binding to FcRs has been linked to depletion of intratumoral T regulatory cells (Treg). In agreement with previous studies, we found that Tregs infiltrating CT26, B16-F1, and autochthonous Braf V600E Pten -/- melanoma tumors had higher expression of surface CTLA-4 (sCTLA-4) than other T-cell subsets, and anti-CTLA-4 treatment led to FcR-dependent depletion of Tregs infiltrating CT26 tumors. This Treg depletion coincided with activation and degranulation of intratumoral natural killer cells. Similarly, in non-small cell lung cancer (NSCLC) and melanoma patient-derived tumor tissue, Tregs had higher sCTLA-4 expression than other intratumoral T-cell subsets, and Tregs infiltrating NSCLC expressed more sCTLA-4 than circulating Tregs. Patients with cutaneous melanoma who benefited from ipilimumab, a mAb targeting CTLA-4, had higher intratumoral CD56 expression, compared with patients who received little to no benefit from this therapy. Furthermore, using the murine CT26 model we found that combination therapy with anti-CTLA-4 plus IL15/IL15Rα complexes enhanced tumor control compared with either monotherapy.

Coral: Clear and Customizable Visualization of Human Kinome Data.

Protein kinases represent one of the largest gene families in eukaryotes and play roles in a wide range of cell signaling processes and human diseases. Current tools for visualizing kinase data in the context of the human kinome superfamily are limited to encoding data through the addition of nodes to a low-resolution image of the kinome tree. We present Coral, a user-friendly interactive web application for visualizing both quantitative and qualitative data. Unlike previous tools, Coral can encode data in three features (node color, node size, and branch color), allows three modes of kinome visualization (the traditional kinome tree as well as radial and dynamic force networks), and generates high-resolution scalable vector graphics files suitable for publication without the need for refinement using graphics editing software. Due to its user-friendly, interactive, and highly customizable design, Coral is broadly applicable to high-throughput studies of the human kinome. The source code and web application are available at github.com/dphansti/CORAL and phanstiel-lab.med.unc.edu/Coral, respectively.

A Landscape of Metabolic Variation across Tumor Types.

Tumor metabolism is reorganized to support proliferation in the face of growth-related stress. Unlike the widespread profiling of changes to metabolic enzyme levels in cancer, comparatively less attention has been paid to the substrates/products of enzyme-catalyzed reactions, small-molecule metabolites. We developed an informatic pipeline to concurrently analyze metabolomics data from over 900 tissue samples spanning seven cancer types, revealing extensive heterogeneity in metabolic changes relative to normal tissue across cancers of different tissues of origin. Despite this heterogeneity, a number of metabolites were recurrently differentially abundant across many cancers, such as lactate and acyl-carnitine species. Through joint analysis of metabolomic data alongside clinical features of patient samples, we also identified a small number of metabolites, including several polyamines and kynurenine, which were associated with aggressive tumors across several tumor types. Our findings offer a glimpse onto common patterns of metabolic reprogramming across cancers, and the work serves as a large-scale resource accessible via a web application (http://www.sanderlab.org/pancanmet).

KMT2C Mutations in Diffuse-Type Gastric Adenocarcinoma Promote Epithelial-to-Mesenchymal Transition.

PURPOSE: Lauren diffuse-type gastric adenocarcinomas (DGAs) are generally genomically stable. We identified lysine (K)-specific methyltransferase 2C (KMT2C) as a frequently mutated gene and examined its role in DGA progression. EXPERIMENTAL DESIGN: We performed whole exome sequencing on tumor samples of 27 patients with DGA who underwent gastrectomy. Lysine (K)-specific methyltransferase 2C (KMT2C) was analyzed in DGA cell lines and in patient tumors. RESULTS: KMT2C was the most frequently mutated gene (11 of 27 tumors [41%]). KMT2C expression by immunohistochemistry in tumors from 135 patients with DGA undergoing gastrectomy inversely correlated with more advanced tumor stage (P = 0.023) and worse overall survival (P = 0.017). KMT2C shRNA knockdown in non-transformed HFE-145 gastric epithelial cells promoted epithelial-to-mesenchymal transition (EMT) as demonstrated by increased expression of EMT-related proteins N-cadherin and Slug. Migration and invasion in gastric epithelial cells following KMT2C knockdown increased by 47- to 88-fold. In the DGA cell lines MKN-45 and SNU-668, which have lost KMT2C expression, KMT2C re-expression decreased expression of EMT-related proteins, reduced cell migration by 52% to 60%, and reduced cell invasion by 50% to 74%. Flank xenografts derived from KMT2C-expressing DGA organoids, compared with wild-type organoids, grew more slowly and lost their infiltrative leading edge. EMT can lead to the acquisition of cancer stem cell (CSC) phenotypes. KMT2C re-expression in DGA cell lines reduced spheroid formation by 77% to 78% and reversed CSC resistance to chemotherapy via promotion of DNA damage and apoptosis. CONCLUSIONS: KMT2C is frequently mutated in certain populations with DGA. KMT2C loss in DGA promotes EMT and is associated with worse overall survival.

Computational Pipeline for the PGV-001 Neoantigen Vaccine Trial.

This paper describes the sequencing protocol and computational pipeline for the PGV-001 personalized vaccine trial. PGV-001 is a therapeutic peptide vaccine targeting neoantigens identified from patient tumor samples. Peptides are selected by a computational pipeline that identifies mutations from tumor/normal exome sequencing and ranks mutant sequences by a combination of predicted Class I MHC affinity and abundance estimated from tumor RNA. The personalized genomic vaccine (PGV) pipeline is modular and consists of independently usable tools and software libraries. We hope that the functionality of these tools may extend beyond the specifics of the PGV-001 trial and enable other research groups in their own neoantigen investigations.

Somatic Mutations and Neoepitope Homology in Melanomas Treated with CTLA-4 Blockade.

Immune checkpoint inhibitors are promising treatments for patients with a variety of malignancies. Toward understanding the determinants of response to immune checkpoint inhibitors, it was previously demonstrated that the presence of somatic mutations is associated with benefit from checkpoint inhibition. A hypothesis was posited that neoantigen homology to pathogens may in part explain the link between somatic mutations and response. To further examine this hypothesis, we reanalyzed cancer exome data obtained from our previously published study of 64 melanoma patients treated with CTLA-4 blockade and a new dataset of RNA-Seq data from 24 of these patients. We found that the ability to accurately predict patient benefit did not increase as the analysis narrowed from somatic mutation burden, to inclusion of only those mutations predicted to be MHC class I neoantigens, to only including those neoantigens that were expressed or that had homology to pathogens. The only association between somatic mutation burden and response was found when examining samples obtained prior to treatment. Neoantigen and expressed neoantigen burden were also associated with response, but neither was more predictive than somatic mutation burden. Neither the previously described tetrapeptide signature nor an updated method to evaluate neoepitope homology to pathogens was more predictive than mutation burden. Cancer Immunol Res; 5(1); 84-91. ©2016 AACR.