Oncogenic Signaling Pathways in The Cancer Genome Atlas.

Genetic alterations in signaling pathways that control cell-cycle progression, apoptosis, and cell growth are common hallmarks of cancer, but the extent, mechanisms, and co-occurrence of alterations in these pathways differ between individual tumors and tumor types. Using mutations, copy-number changes, mRNA expression, gene fusions and DNA methylation in 9,125 tumors profiled by The Cancer Genome Atlas (TCGA), we analyzed the mechanisms and patterns of somatic alterations in ten canonical pathways: cell cycle, Hippo, Myc, Notch, Nrf2, PI-3-Kinase/Akt, RTK-RAS, TGFß signaling, p53 and ß-catenin/Wnt. We charted the detailed landscape of pathway alterations in 33 cancer types, stratified into 64 subtypes, and identified patterns of co-occurrence and mutual exclusivity. Eighty-nine percent of tumors had at least one driver alteration in these pathways, and 57% percent of tumors had at least one alteration potentially targetable by currently available drugs. Thirty percent of tumors had multiple targetable alterations, indicating opportunities for combination therapy.

scPerturb: harmonized single-cell perturbation data.

Analysis across a growing number of single-cell perturbation datasets is hampered by poor data interoperability. To facilitate development and benchmarking of computational methods, we collect a set of 44 publicly available single-cell perturbation-response datasets with molecular readouts, including transcriptomics, proteomics and epigenomics. We apply uniform quality control pipelines and harmonize feature annotations. The resulting information resource, scPerturb, enables development and testing of computational methods, and facilitates comparison and integration across datasets. We describe energy statistics (E-statistics) for quantification of perturbation effects and significance testing, and demonstrate E-distance as a general distance measure between sets of single-cell expression profiles. We illustrate the application of E-statistics for quantifying similarity and efficacy of perturbations. The perturbation-response datasets and E-statistics computation software are publicly available at scperturb.org. This work provides an information resource for researchers working with single-cell perturbation data and recommendations for experimental design, including optimal cell counts and read depth.

LLGL2 rescues nutrient stress by promoting leucine uptake in ER+ breast cancer.

Drosophila Lgl and its mammalian homologues, LLGL1 and LLGL2, are scaffolding proteins that regulate the establishment of apical-basal polarity in epithelial cells1,2. Whereas Lgl functions as a tumour suppressor in Drosophila1, the roles of mammalian LLGL1 and LLGL2 in cancer are unclear. The majority (about 75%) of breast cancers express oestrogen receptors (ERs)3, and patients with these tumours receive endocrine treatment4. However, the development of resistance to endocrine therapy and metastatic progression are leading causes of death for patients with ER+ disease4. Here we report that, unlike LLGL1, LLGL2 is overexpressed in ER+ breast cancer and promotes cell proliferation under nutrient stress. LLGL2 regulates cell surface levels of a leucine transporter, SLC7A5, by forming a trimeric complex with SLC7A5 and a regulator of membrane fusion, YKT6, to promote leucine uptake and cell proliferation. The oestrogen receptor targets LLGL2 expression. Resistance to endocrine treatment in breast cancer cells was associated with SLC7A5- and LLGL2-dependent adaption to nutrient stress. SLC7A5 was necessary and sufficient to confer resistance to tamoxifen treatment, identifying SLC7A5 as a potential therapeutic target for overcoming resistance to endocrine treatments in breast cancer. Thus, LLGL2 functions as a promoter of tumour growth and not as a tumour suppressor in ER+ breast cancer. Beyond breast cancer, adaptation to nutrient stress is critically important5, and our findings identify an unexpected role for LLGL2 in this process.

Proteomic Dynamics of Breast Cancer Cell Lines Identifies Potential Therapeutic Protein Targets.

Treatment and relevant targets for breast cancer (BC) remain limited, especially for triple-negative BC (TNBC). We identified 6091 proteins of 76 human BC cell lines using data-independent acquisition (DIA). Integrating our proteomic findings with prior multi-omics datasets, we found that including proteomics data improved drug sensitivity predictions and provided insights into the mechanisms of action. We subsequently profiled the proteomic changes in nine cell lines (five TNBC and four non-TNBC) treated with EGFR/AKT/mTOR inhibitors. In TNBC, metabolism pathways were dysregulated after EGFR/mTOR inhibitor treatment, while RNA modification and cell cycle pathways were affected by AKT inhibitor. This systematic multi-omics and in-depth analysis of the proteome of BC cells can help prioritize potential therapeutic targets and provide insights into adaptive resistance in TNBC.

StonPy: a tool to parse and query collections of SBGN maps in a graph database.

SUMMARY: The systems biology graphical notation (SBGN) has become the de facto standard for the graphical representation of molecular maps. Having rapid and easy access to the content of large collections of maps is necessary to perform semantic or graph-based analysis of these resources. To this end, we propose StonPy, a new tool to store and query SBGN maps in a Neo4j graph database. StonPy notably includes a data model that takes into account all three SBGN languages and a completion module to automatically build valid SBGN maps from query results. StonPy is built as a library that can be integrated into other software and offers a command-line interface that allows users to easily perform all operations. AVAILABILITY AND IMPLEMENTATION: StonPy is implemented in Python 3 under a GPLv3 license. Its code and complete documentation are freely available from https://github.com/adrienrougny/stonpy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

COVIDpro: Database for Mining Protein Dysregulation in Patients with COVID-19.

The ongoing pandemic of the coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 still has limited treatment options. Our understanding of the molecular dysregulations that occur in response to infection remains incomplete. We developed a web application COVIDpro (https://www.guomics.com/covidPro/) that includes proteomics data obtained from 41 original studies conducted in 32 hospitals worldwide, involving 3077 patients and covering 19 types of clinical specimens, predominantly plasma and serum. The data set encompasses 53 protein expression matrices, comprising a total of 5434 samples and 14,403 unique proteins. We identified a panel of proteins that exhibit significant dysregulation, enabling the classification of COVID-19 patients into severe and non-severe disease categories. The proteomic signatures achieved promising results in distinguishing severe cases, with a mean area under the curve of 0.87 and accuracy of 0.80 across five independent test sets. COVIDpro serves as a valuable resource for testing hypotheses and exploring potential targets for novel treatments in COVID-19 patients.

CellMiner Cross-Database (CellMinerCDB) version 1.2: Exploration of patient-derived cancer cell line pharmacogenomics.

CellMiner Cross-Database (CellMinerCDB, discover.nci.nih.gov/cellminercdb) allows integration and analysis of molecular and pharmacological data within and across cancer cell line datasets from the National Cancer Institute (NCI), Broad Institute, Sanger/MGH and MD Anderson Cancer Center (MDACC). We present CellMinerCDB 1.2 with updates to datasets from NCI-60, Broad Cancer Cell Line Encyclopedia and Sanger/MGH, and the addition of new datasets, including NCI-ALMANAC drug combination, MDACC Cell Line Project proteomic, NCI-SCLC DNA copy number and methylation data, and Broad methylation, genetic dependency and metabolomic datasets. CellMinerCDB (v1.2) includes several improvements over the previously published version: (i) new and updated datasets; (ii) support for pattern comparisons and multivariate analyses across data sources; (iii) updated annotations with drug mechanism of action information and biologically relevant multigene signatures; (iv) analysis speedups via caching; (v) a new dataset download feature; (vi) improved visualization of subsets of multiple tissue types; (vii) breakdown of univariate associations by tissue type; and (viii) enhanced help information. The curation and common annotations (e.g. tissues of origin and identifiers) provided here across pharmacogenomic datasets increase the utility of the individual datasets to address multiple researcher question types, including data reproducibility, biomarker discovery and multivariate analysis of drug activity.

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.

The digital revolution in phenotyping.

Phenotypes have gained increased notoriety in the clinical and biological domain owing to their application in numerous areas such as the discovery of disease genes and drug targets, phylogenetics and pharmacogenomics. Phenotypes, defined as observable characteristics of organisms, can be seen as one of the bridges that lead to a translation of experimental findings into clinical applications and thereby support 'bench to bedside' efforts. However, to build this translational bridge, a common and universal understanding of phenotypes is required that goes beyond domain-specific definitions. To achieve this ambitious goal, a digital revolution is ongoing that enables the encoding of data in computer-readable formats and the data storage in specialized repositories, ready for integration, enabling translational research. While phenome research is an ongoing endeavor, the true potential hidden in the currently available data still needs to be unlocked, offering exciting opportunities for the forthcoming years. Here, we provide insights into the state-of-the-art in digital phenotyping, by means of representing, acquiring and analyzing phenotype data. In addition, we provide visions of this field for future research work that could enable better applications of phenotype data.

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.

rcellminer: exploring molecular profiles and drug response of the NCI-60 cell lines in R.

PURPOSE: The rcellminer R package provides a wide range of functionality to help R users access and explore molecular profiling and drug response data for the NCI-60. The package enables flexible programmatic access to CellMiner's unparalleled breadth of NCI-60 data, including gene and protein expression, copy number, whole exome mutations, as well as activity data for ∼21K compounds, with information on their structure, mechanism of action and repeat screens. Functions are available to easily visualize compound structures, activity patterns and molecular feature profiles. Additionally, embedded R Shiny applications allow interactive data exploration. AVAILABILITY AND IMPLEMENTATION: rcellminer is compatible with R 3.2 and above on Windows, Mac OS X and Linux. The package, documentation, tutorials and Shiny-based applications are available through Bioconductor (http://www.bioconductor.org/packages/rcellminer); ongoing updates will occur according to the Bioconductor release schedule with new CellMiner data. The package is free and open-source (LGPL 3). CONTACT: lunaa@cbio.mskcc.org or vinodh.rajapakse@nih.gov.

Predicted Role of NAD Utilization in the Control of Circadian Rhythms during DNA Damage Response.

The circadian clock is a set of regulatory steps that oscillate with a period of approximately 24 hours influencing many biological processes. These oscillations are robust to external stresses, and in the case of genotoxic stress (i.e. DNA damage), the circadian clock responds through phase shifting with primarily phase advancements. The effect of DNA damage on the circadian clock and the mechanism through which this effect operates remains to be thoroughly investigated. Here we build an in silico model to examine damage-induced circadian phase shifts by investigating a possible mechanism linking circadian rhythms to metabolism. The proposed model involves two DNA damage response proteins, SIRT1 and PARP1, that are each consumers of nicotinamide adenine dinucleotide (NAD), a metabolite involved in oxidation-reduction reactions and in ATP synthesis. This model builds on two key findings: 1) that SIRT1 (a protein deacetylase) is involved in both the positive (i.e. transcriptional activation) and negative (i.e. transcriptional repression) arms of the circadian regulation and 2) that PARP1 is a major consumer of NAD during the DNA damage response. In our simulations, we observe that increased PARP1 activity may be able to trigger SIRT1-induced circadian phase advancements by decreasing SIRT1 activity through competition for NAD supplies. We show how this competitive inhibition may operate through protein acetylation in conjunction with phosphorylation, consistent with reported observations. These findings suggest a possible mechanism through which multiple perturbations, each dominant during different points of the circadian cycle, may result in the phase advancement of the circadian clock seen during DNA damage.

Using drug response data to identify molecular effectors, and molecular "omic" data to identify candidate drugs in cancer.

The current convergence of molecular and pharmacological data provides unprecedented opportunities to gain insights into the relationships between the two types of data. Multiple forms of large-scale molecular data, including but not limited to gene and microRNA transcript expression, DNA somatic and germline variations from next-generation DNA and RNA sequencing, and DNA copy number from array comparative genomic hybridization are all potentially informative when one attempts to recognize the panoply of potentially influential events both for cancer progression and therapeutic outcome. Concurrently, there has also been a substantial expansion of the pharmacological data being accrued in a systematic fashion. For cancer cell lines, the National Cancer Institute cell line panel (NCI-60), the Cancer Cell Line Encyclopedia (CCLE), and the collaborative Genomics of Drug Sensitivity in Cancer (GDSC) databases all provide subsets of these forms of data. For the patient-derived data, The Cancer Genome Atlas (TCGA) provides analogous forms of genomic information along with treatment histories. Integration of these data in turn relies on the fields of statistics and statistical learning. Multiple algorithmic approaches may be chosen, depending on the data being considered, and the nature of the question being asked. Combining these algorithms with prior biological knowledge, the results of molecular biological studies, and the consideration of genes as pathways or functional groups provides both the challenge and the potential of the field. The ultimate goal is to provide a paradigm shift in the way that drugs are selected to provide a more targeted and efficacious outcome for the patient.

PathVisio-Faceted Search: an exploration tool for multi-dimensional navigation of large pathways.

PURPOSE: The PathVisio-Faceted Search plugin helps users explore and understand complex pathways by overlaying experimental data and data from webservices, such as Ensembl BioMart, onto diagrams drawn using formalized notations in PathVisio. The plugin then provides a filtering mechanism, known as a faceted search, to find and highlight diagram nodes (e.g. genes and proteins) of interest based on imported data. The tool additionally provides a flexible scripting mechanism to handle complex queries. AVAILABILITY: The PathVisio-Faceted Search plugin is compatible with PathVisio 3.0 and above. PathVisio is compatible with Windows, Mac OS X and Linux. The plugin, documentation, example diagrams and Groovy scripts are available at http://PathVisio.org/wiki/PathVisioFacetedSearchHelp. The plugin is free, open-source and licensed by the Apache 2.0 License.

drGAT: Attention-Guided Gene Assessment of Drug Response Utilizing a Drug-Cell-Gene Heterogeneous Network.

Drug development is a lengthy process with a high failure rate. Increasingly, machine learning is utilized to facilitate the drug development processes. These models aim to enhance our understanding of drug characteristics, including their activity in biological contexts. However, a major challenge in drug response (DR) prediction is model interpretability as it aids in the validation of findings. This is important in biomedicine, where models need to be understandable in comparison with established knowledge of drug interactions with proteins. drGAT, a graph deep learning model, leverages a heterogeneous graph composed of relationships between proteins, cell lines, and drugs. drGAT is designed with two objectives: DR prediction as a binary sensitivity prediction and elucidation of drug mechanism from attention coefficients. drGAT has demonstrated superior performance over existing models, achieving 78% accuracy (and precision), and 76% F1 score for 269 DNA-damaging compounds of the NCI60 drug response dataset. To assess the model's interpretability, we conducted a review of drug-gene co-occurrences in Pubmed abstracts in comparison to the top 5 genes with the highest attention coefficients for each drug. We also examined whether known relationships were retained in the model by inspecting the neighborhoods of topoisomerase-related drugs. For example, our model retained TOP1 as a highly weighted predictive feature for irinotecan and topotecan, in addition to other genes that could potentially be regulators of the drugs. Our method can be used to accurately predict sensitivity to drugs and may be useful in the identification of biomarkers relating to the treatment of cancer patients.

Sarcoma_CellminerCDB: A tool to interrogate the genomic and functional characteristics of a comprehensive collection of sarcoma cell lines.

Sarcomas are a diverse group of rare malignancies composed of multiple different clinical and molecular subtypes. Due to their rarity and heterogeneity, basic, translational, and clinical research in sarcoma has trailed behind that of other cancers. Outcomes for patients remain generally poor due to an incomplete understanding of disease biology and a lack of novel therapies. To address some of the limitations impeding preclinical sarcoma research, we have developed Sarcoma_CellMinerCDB, a publicly available interactive tool that merges publicly available sarcoma cell line data and newly generated omics data to create a comprehensive database of genomic, transcriptomic, methylomic, proteomic, metabolic, and pharmacologic data on 133 annotated sarcoma cell lines. The reproducibility, functionality, biological relevance, and therapeutic applications of Sarcoma_CellMinerCDB described herein are powerful tools to address and generate biological questions and test hypotheses for translational research. Sarcoma_CellMinerCDB (https://discover.nci.nih.gov/SarcomaCellMinerCDB) aims to contribute to advancing the preclinical study of sarcoma.

PathVisio-Validator: a rule-based validation plugin for graphical pathway notations.

PURPOSE: The PathVisio-Validator plugin aims to simplify the task of producing biological pathway diagrams that follow graphical standardized notations, such as Molecular Interaction Maps or the Systems Biology Graphical Notation. This plugin assists in the creation of pathway diagrams by ensuring correct usage of a notation, and thereby reducing ambiguity when diagrams are shared among biologists. Rulesets, needed in the validation process, can be generated for any graphical notation that a developer desires, using either Schematron or Groovy. The plugin also provides support for filtering validation results, validating on a subset of rules, and distinguishing errors and warnings.

Software support for SBGN maps: SBGN-ML and LibSBGN.

MOTIVATION: LibSBGN is a software library for reading, writing and manipulating Systems Biology Graphical Notation (SBGN) maps stored using the recently developed SBGN-ML file format. The library (available in C++ and Java) makes it easy for developers to add SBGN support to their tools, whereas the file format facilitates the exchange of maps between compatible software applications. The library also supports validation of maps, which simplifies the task of ensuring compliance with the detailed SBGN specifications. With this effort we hope to increase the adoption of SBGN in bioinformatics tools, ultimately enabling more researchers to visualize biological knowledge in a precise and unambiguous manner. AVAILABILITY AND IMPLEMENTATION: Milestone 2 was released in December 2011. Source code, example files and binaries are freely available under the terms of either the LGPL v2.1+ or Apache v2.0 open source licenses from http://libsbgn.sourceforge.net. CONTACT: sbgn-libsbgn@lists.sourceforge.net.

cyjShiny: A cytoscape.js R Shiny Widget for network visualization and analysis.

cyjShiny is an open-source R package that allows users to embed network visualization into Shiny apps and R Markdown documents. cyjShiny (https://github.com/cytoscape/cyjShiny) builds on the cytoscape.js Javascript graph library. Additionally, the package provides helper functions to convert common R data representations (e.g., data.frame) into forms compatible with cytoscape.js.

Mapping the functional interactions at the tumor-immune checkpoint interface.

The interactions between tumor intrinsic processes and immune checkpoints can mediate immune evasion by cancer cells and responses to immunotherapy. It is, however, challenging to identify functional interactions due to the prohibitively complex molecular landscape of the tumor-immune interfaces. We address this challenge with a statistical analysis framework, immuno-oncology gene interaction maps (ImogiMap). ImogiMap quantifies and statistically validates tumor-immune checkpoint interactions based on their co-associations with immune-associated phenotypes. The outcome is a catalog of tumor-immune checkpoint interaction maps for diverse immune-associated phenotypes. Applications of ImogiMap recapitulate the interaction of SERPINB9 and immune checkpoints with interferon gamma (IFNγ) expression. Our analyses suggest that CD86-CD70 and CD274-CD70 immunoregulatory interactions are significantly associated with IFNγ expression in uterine corpus endometrial carcinoma and basal-like breast cancer, respectively. The open-source ImogiMap software and user-friendly web application will enable future applications of ImogiMap. Such applications may guide the discovery of previously unknown tumor-immune interactions and immunotherapy targets.