Pathway-specific model estimation for improved pathway annotation by network crosstalk

Disruption of the TP53 locus in osteosarcoma leads to TP53 promoter gene fusions and restoration of parts of the TP53 signalling pathway

TP53 is the most frequently mutated gene in human cancer. This gene shows not only loss-of-function mutations but also recurrent missense mutations with gain-of-function activity. We have studied the primary bone malignancy osteosarcoma, which harbours one of the most rearranged genomes of all cancers. This is odd since it primarily affects children and adolescents who have not lived the long life thought necessary to accumulate massive numbers of mutations. In osteosarcoma, TP53 is often disrupted by structural variants. Here, we show through combined whole-genome and transcriptome analyses of 148 osteosarcomas that TP53 structural variants commonly result in loss of coding parts of the gene while simultaneously preserving and relocating the promoter region. The transferred TP53 promoter region is fused to genes previously implicated in cancer development. Paradoxically, these erroneously upregulated genes are significantly associated with the TP53 signalling pathway itself. This suggests that while the classical tumour suppressor activities of TP53 are lost, certain parts of the TP53 signalling pathway that are necessary for cancer cell survival and proliferation are retained. In line with this, our data suggest that transposition of the TP53 promoter is an early event that allows for a new normal state of genome-wide rearrangements in osteosarcoma. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Benchmarking enrichment analysis methods with the disease pathway network

Enrichment analysis (EA) is a common approach to gain functional insights from genome-scale experiments. As a consequence, a large number of EA methods have been developed, yet it is unclear from previous studies which method is the best for a given dataset. The main issues with previous benchmarks include the complexity of correctly assigning true pathways to a test dataset, and lack of generality of the evaluation metrics, for which the rank of a single target pathway is commonly used. We here provide a generalized EA benchmark and apply it to the most widely used EA methods, representing all four categories of current approaches. The benchmark employs a new set of 82 curated gene expression datasets from DNA microarray and RNA-Seq experiments for 26 diseases, of which only 13 are cancers. In order to address the shortcomings of the single target pathway approach and to enhance the sensitivity evaluation, we present the Disease Pathway Network, in which related Kyoto Encyclopedia of Genes and Genomes pathways are linked. We introduce a novel approach to evaluate pathway EA by combining sensitivity and specificity to provide a balanced evaluation of EA methods. This approach identifies Network Enrichment Analysis methods as the overall top performers compared with overlap-based methods. By using randomized gene expression datasets, we explore the null hypothesis bias of each method, revealing that most of them produce skewed P-values.

Network-Based Methods and Their Applications in Drug Discovery

Drug discovery is time-consuming, expensive, and predominantly follows the "one drug → one target → one disease" paradigm. With the rapid development of systems biology and network pharmacology, a novel drug discovery paradigm, "multidrug → multitarget → multidisease", has emerged. This new holistic paradigm of drug discovery aligns well with the essence of networks, leading to the emergence of network-based methods in the field of drug discovery. In this Perspective, we initially introduce the concept and data sources of networks and highlight classical methodologies employed in network-based methods. Subsequently, we focus on the practical applications of network-based methods across various areas of drug discovery, such as target prediction, virtual screening, prediction of drug therapeutic effects or adverse drug events, and elucidation of molecular mechanisms. In addition, we provide representative web servers for researchers to use network-based methods in specific applications. Finally, we discuss several challenges of network-based methods and the directions for future development. In a word, network-based methods could serve as powerful tools to accelerate drug discovery.

FAK family proteins regulate in vivo breast cancer metastasis via distinct mechanisms

Breast cancer is the most commonly diagnosed malignancy and the major leading cause of tumor-related deaths in women. It is estimated that the majority of breast tumor-related deaths are a consequence of metastasis, to which no cure exists at present. The FAK family proteins Proline-rich tyrosine kinase (PYK2) and focal adhesion kinase (FAK) are highly expressed in breast cancer, but the exact cellular and signaling mechanisms by which they regulate in vivo tumor cell invasiveness and consequent metastatic dissemination are mostly unknown. Using a PYK2 and FAK knockdown xenograft model we show here, for the first time, that ablation of either PYK2 or FAK decreases primary tumor size and significantly reduces Tumor MicroEnvironment of Metastasis (TMEM) doorway activation, leading to decreased intravasation and reduced spontaneous lung metastasis. Intravital imaging analysis further demonstrates that PYK2, but not FAK, regulates a motility phenotype switch between focal adhesion-mediated fast motility and invadopodia-dependent, ECM-degradation associated slow motility within the primary tumor. Furthermore, we validate our in vivo and intravital imaging results with integrated transcriptomic and proteomic data analysis from xenograft knockdown tumors and reveal new and distinct pathways by which these two homologous kinases regulate breast tumor cell invasiveness and consequent metastatic dissemination. Our findings identify PYK2 and FAK as novel mediators of mammary tumor progression and metastasis and as candidate therapeutic targets for breast cancer metastasis.

Identifying cellular cancer mechanisms through pathway-driven data integration

MOTIVATION

Cancer is a genetic disease in which accumulated mutations of driver genes induce a functional reorganization of the cell by reprogramming cellular pathways. Current approaches identify cancer pathways as those most internally perturbed by gene expression changes. However, driver genes characteristically perform hub roles between pathways. Therefore, we hypothesize that cancer pathways should be identified by changes in their pathway-pathway relationships.

RESULTS

To learn an embedding space that captures the relationships between pathways in a healthy cell, we propose pathway-driven non-negative matrix tri-factorization. In this space, we determine condition-specific (i.e. diseased and healthy) embeddings of pathways and genes. Based on these embeddings, we define our 'NMTF centrality' to measure a pathway's or gene's functional importance, and our 'moving distance', to measure the change in its functional relationships. We combine both measures to predict 15 genes and pathways involved in four major cancers, predicting 60 gene-cancer associations in total, covering 28 unique genes. To further exploit driver genes' tendency to perform hub roles, we model our network data using graphlet adjacency, which considers nodes adjacent if their interaction patterns form specific shapes (e.g. paths or triangles). We find that the predicted genes rewire pathway-pathway interactions in the immune system and provide literary evidence that many are druggable (15/28) and implicated in the associated cancers (47/60). We predict six druggable cancer-specific drug targets.

AVAILABILITY AND IMPLEMENTATION

The code and data are available at: https://gitlab.bsc.es/swindels/pathway_driven_nmtf.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A systematic evaluation of the computational tools for ligand-receptor-based cell-cell interaction inference

Cell-cell interactions (CCIs) are essential for multicellular organisms to coordinate biological processes and functions. One classical type of CCI interaction is between secreted ligands and cell surface receptors, i.e. ligand-receptor (LR) interactions. With the recent development of single-cell technologies, a large amount of single-cell ribonucleic acid (RNA) sequencing (scRNA-Seq) data has become widely available. This data availability motivated the single-cell-resolution study of CCIs, particularly LR-based CCIs. Dozens of computational methods and tools have been developed to predict CCIs by identifying LR-based CCIs. Many of these tools have been theoretically reviewed. However, there is little study on current LR-based CCI prediction tools regarding their performance and running results on public scRNA-Seq datasets. In this work, to fill this gap, we tested and compared nine of the most recent computational tools for LR-based CCI prediction. We used 15 well-studied scRNA-Seq samples that correspond to approximately 100K single cells under different experimental conditions for testing and comparison. Besides briefing the methodology used in these nine tools, we summarized the similarities and differences of these tools in terms of both LR prediction and CCI inference between cell types. We provided insight into using these tools to make meaningful discoveries in understanding cell communications.

From Traditional Ethnopharmacology to Modern Natural Drug Discovery: A Methodology Discussion and Specific Examples

Ethnopharmacology, through the description of the beneficial effects of plants, has provided an early framework for the therapeutic use of natural compounds. Natural products, either in their native form or after crude extraction of their active ingredients, have long been used by different populations and explored as invaluable sources for drug design. The transition from traditional ethnopharmacology to drug discovery has followed a straightforward path, assisted by the evolution of isolation and characterization methods, the increase in computational power, and the development of specific chemoinformatic methods. The deriving extensive exploitation of the natural product chemical space has led to the discovery of novel compounds with pharmaceutical properties, although this was not followed by an analogous increase in novel drugs. In this work, we discuss the evolution of ideas and methods, from traditional ethnopharmacology to in silico drug discovery, applied to natural products. We point out that, in the past, the starting point was the plant itself, identified by sustained ethnopharmacological research, with the active compound deriving after extensive analysis and testing. In contrast, in recent years, the active substance has been pinpointed by computational methods (in silico docking and molecular dynamics, network pharmacology), followed by the identification of the plant(s) containing the active ingredient, identified by existing or putative ethnopharmacological information. We further stress the potential pitfalls of recent in silico methods and discuss the absolute need for in vitro and in vivo validation as an absolute requirement. Finally, we present our contribution to natural products’ drug discovery by discussing specific examples, applying the whole continuum of this rapidly evolving field. In detail, we report the isolation of novel antiviral compounds, based on natural products active against influenza and SARS-CoV-2 and novel substances active on a specific GPCR, OXER1.

Benefits and Challenges of Pre-clustered Network-Based Pathway Analysis

Functional analysis of gene sets derived from experiments is typically done by pathway annotation. Although many algorithms exist for analyzing the association between a gene set and a pathway, an issue which is generally ignored is that gene sets often represent multiple pathways. In such cases an association to a pathway is weakened by the presence of genes associated with other pathways. A way to counteract this is to cluster the gene set into more homogenous parts before performing pathway analysis on each module. We explored whether network-based pre-clustering of a query gene set can improve pathway analysis. The methods MCL, Infomap, and MGclus were used to cluster the gene set projected onto the FunCoup network. We characterized how well these methods are able to detect individual pathways in multi-pathway gene sets, and applied each of the clustering methods in combination with four pathway analysis methods: Gene Enrichment Analysis, BinoX, NEAT, and ANUBIX. Using benchmarks constructed from the KEGG pathway database we found that clustering can be beneficial by increasing the sensitivity of pathway analysis methods and by providing deeper insights of biological mechanisms related to the phenotype under study. However, keeping a high specificity is a challenge. For ANUBIX, clustering caused a minor loss of specificity, while for BinoX and NEAT it caused an unacceptable loss of specificity. GEA had very low sensitivity both before and after clustering. The choice of clustering method only had a minor effect on the results. We show examples of this approach and conclude that clustering can improve overall pathway annotation performance, but should only be used if the used enrichment method has a low false positive rate.

Network Crosstalk as a Basis for Drug Repurposing

The need for systematic drug repurposing has seen a steady increase over the past decade and may be particularly valuable to quickly remedy unexpected pandemics. The abundance of functional interaction data has allowed mapping of substantial parts of the human interactome modeled using functional association networks, favoring network-based drug repurposing. Network crosstalk-based approaches have never been tested for drug repurposing despite their success in the related and more mature field of pathway enrichment analysis. We have, therefore, evaluated the top performing crosstalk-based approaches for drug repurposing. Additionally, the volume of new interaction data as well as more sophisticated network integration approaches compelled us to construct a new benchmark for performance assessment of network-based drug repurposing tools, which we used to compare network crosstalk-based methods with a state-of-the-art technique. We find that network crosstalk-based drug repurposing is able to rival the state-of-the-art method and in some cases outperform it.

PathBIX-a web server for network-based pathway annotation with adaptive null models

Motivation

Pathway annotation is a vital tool for interpreting and giving meaning to experimental data in life sciences. Numerous tools exist for this task, where the most recent generation of pathway enrichment analysis tools, network-based methods, utilize biological networks to gain a richer source of information as a basis of the analysis than merely the gene content. Network-based methods use the network crosstalk between the query gene set and the genes in known pathways, and compare this to a null model of random expectation.

Results

We developed PathBIX, a novel web application for network-based pathway analysis, based on the recently published ANUBIX algorithm which has been shown to be more accurate than previous network-based methods. The PathBIX website performs pathway annotation for 21 species, and utilizes prefetched and preprocessed network data from FunCoup 5.0 networks and pathway data from three databases: KEGG, Reactome, and WikiPathways.

Availability

https://pathbix.sbc.su.se/.

Contact

erik.sonnhammer@scilifelab.se.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

FunCoup 5: Functional Association Networks in All Domains of Life, Supporting Directed Links and Tissue-Specificity

FunCoup (https://funcoup.sbc.su.se) is one of the most comprehensive functional association networks of genes/proteins available. Functional associations are inferred by integrating different types of evidence using a redundancy-weighted naïve Bayesian approach, combined with orthology transfer. FunCoup's high coverage comes from using eleven different types of evidence, and extensive transfer of information between species. Since the latest update of the database, the availability of source data has improved drastically, and user expectations on a tool for functional associations have grown. To meet these requirements, we have made a new release of FunCoup with updated source data and improved functionality. FunCoup 5 now includes 22 species from all domains of life, and the source data for evidences, gold standards, and genomes have been updated to the latest available versions. In this new release, directed regulatory links inferred from transcription factor binding can be visualized in the network viewer for the human interactome. Another new feature is the possibility to filter by genes expressed in a certain tissue in the network viewer. FunCoup 5 further includes the SARS-CoV-2 proteome, allowing users to visualize and analyze interactions between SARS-CoV-2 and human proteins in order to better understand COVID-19. This new release of FunCoup constitutes a major advance for the users, with updated sources, new species and improved functionality for analysis of the networks.

ViralLink: An integrated workflow to investigate the effect of SARS-CoV-2 on intracellular signalling and regulatory pathways

The SARS-CoV-2 pandemic of 2020 has mobilised scientists around the globe to research all aspects of the coronavirus virus and its infection. For fruitful and rapid investigation of viral pathomechanisms, a collaborative and interdisciplinary approach is required. Therefore, we have developed ViralLink: a systems biology workflow which reconstructs and analyses networks representing the effect of viruses on intracellular signalling. These networks trace the flow of signal from intracellular viral proteins through their human binding proteins and downstream signalling pathways, ending with transcription factors regulating genes differentially expressed upon viral exposure. In this way, the workflow provides a mechanistic insight from previously identified knowledge of virally infected cells. By default, the workflow is set up to analyse the intracellular effects of SARS-CoV-2, requiring only transcriptomics counts data as input from the user: thus, encouraging and enabling rapid multidisciplinary research. However, the wide-ranging applicability and modularity of the workflow facilitates customisation of viral context, a priori interactions and analysis methods. Through a case study of SARS-CoV-2 infected bronchial/tracheal epithelial cells, we evidence the functionality of the workflow and its ability to identify key pathways and proteins in the cellular response to infection. The application of ViralLink to different viral infections in a context specific manner using different available transcriptomics datasets will uncover key mechanisms in viral pathogenesis.