Network Diffusion Promotes the Integrative Analysis of Multiple Omics

Network dynamics and therapeutic aspects of mRNA and protein markers with the recurrence sites of pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is a deadly disease that typically manifests late patient presentation and poor outcomes. Furthermore, PDAC recurrence is a common challenge. Distinct patterns of PDAC recurrence have been associated with differential activation of immune pathway-related genes and specific inflammatory responses in their tumour microenvironment. However, the molecular associations between and within cellular components that underpin PDAC recurrence require further development, especially from a multi-omics integration perspective. In this study, we identified stable molecular associations across multiple PDAC recurrences and utilised integrative analytics to identify stable and novel associations via simultaneous feature selection. Spatial transcriptome and proteome datasets were used to perform univariate analysis, Spearman partial correlation analysis, and univariate analyses by Machine Learning methods, including regularised canonical correlation analysis and sparse partial least squares. Furthermore, networks were constructed for reported and new stable associations. Our findings revealed gene and protein associations across multiple PDAC recurrence groups, which can provide a better understanding of the multi-layer disease mechanisms that contribute to PDAC recurrence. These findings may help to provide novel association targets for clinical studies for constructing precision medicine and personalised surveillance tools for patients with PDAC recurrence.

Cuproptosis-Related Genes MTF1 and LIPT1 as Novel Prognostic Biomarker in Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is a life-threatening hematologic malignant disease with high morbidity and mortality in both adults and children. Cuproptosis, a novel mode of cell death, plays an important role in tumor development, but the functional mechanisms of cuproptosis-related genes (CRGs) in AML are unclear. The differential expression of CRGs between tumors such as AML and normal tissues in UCSC XENA, TCGA and GTEx was verified using R (version: 3.6.3). Lasso regression, Cox regression and Nomogram were used to screen for prognostic biomarkers of AML and to construct corresponding prognostic models. Kaplan-Meier analysis, ROC analysis, clinical correlation analysis, immune infiltration analysis and enrichment analysis were used to further investigate the correlation and functional mechanisms of CRGs with AML. The ceRNA regulatory network was used to identify the mRNA-miRNA-lncRNA regulatory axis. Cuproptosis-related genes LIPT1, MTF1, GLS and CDKN2A were highly expressed in AML, while FDX1, LIAS, DLD, DLAT, PDHA1, SLC31A1 and ATP7B were lowly expressed in AML. Lasso regression, Cox regression, Nomogram and calibration curve finally identified MTF1 and LIPT1 as two novel prognostic biomarkers of AML and constructed the corresponding prognostic models. In addition, all 12 CRGs had predictive power for AML, with MTF1, LIAS, SLC31A1 and CDKN2A showing more reliable results. Further analysis showed that ATP7B was closely associated with mutation types such as FLT3, NPM1, RAS and IDH1 R140 in AML, while the expression of MTF1, LIAS and ATP7B in AML was closely associated with immune infiltration. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA) revealed that biological functions such as metal ion transmembrane transporter activity, haptoglobin binding and oxygen carrier activity, pathways such as interferon alpha response, coagulation, UV response DN, apoptosis, hypoxia and heme metabolism all play a role in the development of AML. The ceRNA regulatory network revealed that 6 lncRNAs such as MALAT1, interfere with MTF1 expression through 6 miRNAs such as hsa-miR-32-5p, which in turn affect the development and progression of AML. In addition, APTO-253 has the potential to become an AML-targeted drug. The cuproptosis-related genes MTF1 and LIPT1 can be used as prognostic biomarkers in AML. A total of six lncRNAs, including MALAT1, are involved in the expression and regulation of MTF1 in AML through six miRNAs such as hsa-miR-32-5p.

Are trait-associated genes clustered together in a gene network?

Genome-wide association studies (GWAS) have provided an abundance of information about the genetic variants and their loci that are associated to complex traits and diseases. However, due to linkage disequilibrium (LD) and noncoding regions of loci, it remains a challenge to pinpoint the causal genes. Gene network-based approaches, paired with network diffusion methods, have been proposed to prioritize causal genes and to boost statistical power in GWAS based on the assumption that trait-associated genes are clustered in a gene network. Due to the difficulty in mapping trait-associated variants to genes in GWAS, this assumption has never been directly or rigorously tested empirically. On the other hand, whole exome sequencing (WES) data focuses on the protein-coding regions, directly identifying trait-associated genes. In this study, we tested the assumption by leveraging the recently available exome-based association statistics from the UK Biobank WES data along with two types of networks. We found that almost all trait-associated genes were significantly more proximal to each other than randomly selected genes within both networks. These results support the assumption that trait-associated genes are clustered in gene networks, which can be further leveraged to boost the power of GWAS such as by introducing less stringent p value thresholds.

Curation of causal interactions mediated by genes associated with autism accelerates the understanding of gene-phenotype relationships underlying neurodevelopmental disorders

Autism spectrum disorder (ASD) comprises a large group of neurodevelopmental conditions featuring, over a wide range of severity and combinations, a core set of manifestations (restricted sociality, stereotyped behavior and language impairment) alongside various comorbidities. Common and rare variants in several hundreds of genes and regulatory regions have been implicated in the molecular pathogenesis of ASD along a range of causation evidence strength. Despite significant progress in elucidating the impact of few paradigmatic individual loci, such sheer complexity in the genetic architecture underlying ASD as a whole has hampered the identification of convergent actionable hubs hypothesized to relay between the vastness of risk alleles and the core phenotypes. In turn this has limited the development of strategies that can revert or ameliorate this condition, calling for a systems-level approach to probe the cross-talk of cooperating genes in terms of causal interaction networks in order to make convergences experimentally tractable and reveal their clinical actionability. As a first step in this direction, we have captured from the scientific literature information on the causal links between the genes whose variants have been associated with ASD and the whole human proteome. This information has been annotated in a computer readable format in the SIGNOR database and is made freely available in the resource website. To link this information to cell functions and phenotypes, we have developed graph algorithms that estimate the functional distance of any protein in the SIGNOR causal interactome to phenotypes and pathways. The main novelty of our approach resides in the possibility to explore the mechanistic links connecting the suggested gene-phenotype relations.

Focal adhesion is associated with lithium response in bipolar disorder: evidence from a network-based multi-omics analysis

Lithium (Li) is one of the most effective drugs for treating bipolar disorder (BD), however, there is presently no way to predict response to guide treatment. The aim of this study is to identify functional genes and pathways that distinguish BD Li responders (LR) from BD Li non-responders (NR). An initial Pharmacogenomics of Bipolar Disorder study (PGBD) GWAS of lithium response did not provide any significant results. As a result, we then employed network-based integrative analysis of transcriptomic and genomic data. In transcriptomic study of iPSC-derived neurons, 41 significantly differentially expressed (DE) genes were identified in LR vs NR regardless of lithium exposure. In the PGBD, post-GWAS gene prioritization using the GWA-boosting (GWAB) approach identified 1119 candidate genes. Following DE-derived network propagation, there was a highly significant overlap of genes between the top 500- and top 2000-proximal gene networks and the GWAB gene list (Phypergeometric = 1.28E-09 and 4.10E-18, respectively). Functional enrichment analyses of the top 500 proximal network genes identified focal adhesion and the extracellular matrix (ECM) as the most significant functions. Our findings suggest that the difference between LR and NR was a much greater effect than that of lithium. The direct impact of dysregulation of focal adhesion on axon guidance and neuronal circuits could underpin mechanisms of response to lithium, as well as underlying BD. It also highlights the power of integrative multi-omics analysis of transcriptomic and genomic profiling to gain molecular insights into lithium response in BD.

A novel gene signature to diagnose MASLD in metabolically unhealthy obese individuals

Visceral adipose tissue (VAT) contributes to metabolic dysfunction-associated steatotic liver disease (MASLD), releasing lipogenic substrates and cytokines which promote inflammation. Metabolic healthy obese individuals (MHO) may shift towardsunhealthy ones (MUHO) who develop MASLD, although the mechanisms are still unexplained. Therefore, we aimed to identify dysfunctional pathways and transcriptomic signatures shared by liver and VAT and to outline novel obesity-related biomarkers which feature MASLD in MUHO subjects, at higher risk of progressive liver disease and extrahepatic comorbidities. We performed RNA-sequencing in 167 hepatic samples and in a subset of 79 matched VAT, stratified in MHO and MUHO. A validation analysis was performed in hepatic samples and primary adipocytes from 12 bariatric patients, by qRT-PCR and western blot. We identified a transcriptomic signature that discriminate MUHO vs MHO, including 498 deregulated genes in liver and 189 in VAT. According to pathway and network analyses, oxidative phosphorylation resulted the only significantly downregulated pathway in both tissues in MUHO subjects. Next, we highlighted 5 genes commonly deregulated in liver and VAT, encompassing C6, IGF1, OXA1L, NDUFB11 and KLHL5 and we built a tissue-related score by integrating their expressions. Accordingly to RNAseq data, serum levels of C6 and IGF1, which are the only secreted proteins among those included in the gene signature were downregulated in MUHO vs MHO. Finally, the expression pattern of this 5-genes was confirmed in hepatic and VAT samples. We firstly identified the liver and VAT transcriptional phenotype of MUHO and a gene signature associated with the presence of MASLD in these at risk individuals.

BraneMF: integration of biological networks for functional analysis of proteins

MOTIVATION

The cellular system of a living organism is composed of interacting bio-molecules that control cellular processes at multiple levels. Their correspondences are represented by tightly regulated molecular networks. The increase of omics technologies has favored the generation of large-scale disparate data and the consequent demand for simultaneously using molecular and functional interaction networks: gene co-expression, protein-protein interaction (PPI), genetic interaction and metabolic networks. They are rich sources of information at different molecular levels, and their effective integration is essential to understand cell functioning and their building blocks (proteins). Therefore, it is necessary to obtain informative representations of proteins and their proximity, that are not fully captured by features extracted directly from a single informational level. We propose BraneMF, a novel random walk-based matrix factorization method for learning node representation in a multilayer network, with application to omics data integration.

RESULTS

We test BraneMF with PPI networks of Saccharomyces cerevisiae, a well-studied yeast model organism. We demonstrate the applicability of the learned features for essential multi-omics inference tasks: clustering, function and PPI prediction. We compare it to the state-of-the-art integration methods for multilayer networks. BraneMF outperforms baseline methods by achieving high prediction scores for a variety of downstream tasks. The robustness of results is assessed by an extensive parameter sensitivity analysis.

AVAILABILITY AND IMPLEMENTATION

BraneMF's code is freely available at: https://github.com/Surabhivj/BraneMF, along with datasets, embeddings and result files.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Computational approaches for network-based integrative multi-omics analysis

Advances in omics technologies allow for holistic studies into biological systems. These studies rely on integrative data analysis techniques to obtain a comprehensive view of the dynamics of cellular processes, and molecular mechanisms. Network-based integrative approaches have revolutionized multi-omics analysis by providing the framework to represent interactions between multiple different omics-layers in a graph, which may faithfully reflect the molecular wiring in a cell. Here we review network-based multi-omics/multi-modal integrative analytical approaches. We classify these approaches according to the type of omics data supported, the methods and/or algorithms implemented, their node and/or edge weighting components, and their ability to identify key nodes and subnetworks. We show how these approaches can be used to identify biomarkers, disease subtypes, crosstalk, causality, and molecular drivers of physiological and pathological mechanisms. We provide insight into the most appropriate methods and tools for research questions as showcased around the aetiology and treatment of COVID-19 that can be informed by multi-omics data integration. We conclude with an overview of challenges associated with multi-omics network-based analysis, such as reproducibility, heterogeneity, (biological) interpretability of the results, and we highlight some future directions for network-based integration.

DrDimont: explainable drug response prediction from differential analysis of multi-omics networks

MOTIVATION

While it has been well established that drugs affect and help patients differently, personalized drug response predictions remain challenging. Solutions based on single omics measurements have been proposed, and networks provide means to incorporate molecular interactions into reasoning. However, how to integrate the wealth of information contained in multiple omics layers still poses a complex problem.

RESULTS

We present DrDimont, Drug response prediction from Differential analysis of multi-omics networks. It allows for comparative conclusions between two conditions and translates them into differential drug response predictions. DrDimont focuses on molecular interactions. It establishes condition-specific networks from correlation within an omics layer that are then reduced and combined into heterogeneous, multi-omics molecular networks. A novel semi-local, path-based integration step ensures integrative conclusions. Differential predictions are derived from comparing the condition-specific integrated networks. DrDimont's predictions are explainable, i.e. molecular differences that are the source of high differential drug scores can be retrieved. We predict differential drug response in breast cancer using transcriptomics, proteomics, phosphosite and metabolomics measurements and contrast estrogen receptor positive and receptor negative patients. DrDimont performs better than drug prediction based on differential protein expression or PageRank when evaluating it on ground truth data from cancer cell lines. We find proteomic and phosphosite layers to carry most information for distinguishing drug response.

AVAILABILITY AND IMPLEMENTATION

DrDimont is available on CRAN: https://cran.r-project.org/package=DrDimont.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Overview of methods for characterization and visualization of a protein–protein interaction network in a multi-omics integration context

At the heart of the cellular machinery through the regulation of cellular functions, protein–protein interactions (PPIs) have a significant role. PPIs can be analyzed with network approaches. Construction of a PPI network requires prediction of the interactions. All PPIs form a network. Different biases such as lack of data, recurrence of information, and false interactions make the network unstable. Integrated strategies allow solving these different challenges. These approaches have shown encouraging results for the understanding of molecular mechanisms, drug action mechanisms, and identification of target genes. In order to give more importance to an interaction, it is evaluated by different confidence scores. These scores allow the filtration of the network and thus facilitate the representation of the network, essential steps to the identification and understanding of molecular mechanisms. In this review, we will discuss the main computational methods for predicting PPI, including ones confirming an interaction as well as the integration of PPIs into a network, and we will discuss visualization of these complex data.

Construction and contextualization approaches for protein-protein interaction networks

Protein-protein interaction network (PPIN) analysis is a widely used method to study the contextual role of proteins of interest, to predict novel disease genes, disease or functional modules, and to identify novel drug targets. PPIN-based analysis uses both generic and context-specific networks. Multiple contextualization methodologies have been described, such as shortest-path algorithms, neighborhood-based methods, and diffusion/propagation algorithms. This review discusses these methods, provides intuitive representations of PPIN contextualization, and also examines how the quality of such context-specific networks could be improved by considering additional sources of evidence. As a heuristic, we observe that tasks such as identifying disease genes, drug targets, and protein complexes should consider local neighborhoods, while uncovering disease mechanisms and discovering disease-pathways would gain from diffusion-based construction.

In Silico Investigation of the Biological Implications of Complex DNA Damage with Emphasis in Cancer Radiotherapy through a Systems Biology Approach

Different types of DNA lesions forming in close vicinity, create clusters of damaged sites termed as “clustered/complex DNA damage” and they are considered to be a major challenge for DNA repair mechanisms resulting in significant repair delays and induction of genomic instability. Upon detection of DNA damage, the corresponding DNA damage response and repair (DDR/R) mechanisms are activated. The inability of cells to process clustered DNA lesions efficiently has a great impact on the normal function and survival of cells. If complex lesions are left unrepaired or misrepaired, they can lead to mutations and if persistent, they may lead to apoptotic cell death. In this in silico study, and through rigorous data mining, we have identified human genes that are activated upon complex DNA damage induction like in the case of ionizing radiation (IR) and beyond the standard DNA repair pathways, and are also involved in cancer pathways, by employing stringent bioinformatics and systems biology methodologies. Given that IR can cause repair resistant lesions within a short DNA segment (a few nm), thereby augmenting the hazardous and toxic effects of radiation, we also investigated the possible implication of the most biologically important of those genes in comorbid non-neoplastic diseases through network integration, as well as their potential for predicting survival in cancer patients.

Characterization and comparison of gene-centered human interactomes

The complex web of macromolecular interactions occurring within cells-the interactome-is the backbone of an increasing number of studies, but a clear consensus on the exact structure of this network is still lacking. Different genome-scale maps of human interactome have been obtained through several experimental techniques and functional analyses. Moreover, these maps can be enriched through literature-mining approaches, and different combinations of various 'source' databases have been used in the literature. It is therefore unclear to which extent the various interactomes yield similar results when used in the context of interactome-based approaches in network biology. We compared a comprehensive list of human interactomes on the basis of topology, protein complexes, molecular pathways, pathway cross-talk and disease gene prediction. In a general context of relevant heterogeneity, our study provides a series of qualitative and quantitative parameters that describe the state of the art of human interactomes and guidelines for selecting interactomes in future applications.

Identifying anti-TNF response biomarkers in ulcerative colitis using a diffusion-based signalling model

Motivation

Resistance to anti-TNF therapy in subgroups of ulcerative colitis (UC) patients is a major challenge and incurs significant treatment costs. Identification of patients at risk of nonresponse to anti-TNF is of major clinical importance. To date, no quantitative computational framework exists to develop a complex biomarker for the prognosis of UC treatment. Modelling patient-wise receptor to transcription factor (TF) network connectivity may enable personalized treatment.

Results

We present an approach for quantitative diffusion analysis between receptors and TFs using gene expression data. Key TFs were identified using pandaR. Network connectivities between immune-specific receptor-TF pairs were quantified using network diffusion in UC patients and controls. The patient-specific network could be considered a complex biomarker that separates anti-TNF treatment-resistant and responder patients both in the gene expression dataset used for model development and separate independent test datasets. The model was further validated in rheumatoid arthritis where it successfully discriminated resistant and responder patients to tocilizumab treatment. Our model may contribute to prognostic biomarkers that may identify treatment-resistant and responder subpopulations of UC patients.

Availability and implementation

Software is available at https://github.com/Amy3100/receptor2tfDiffusion.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Multi-layered network-based pathway activity inference using directed random walks: application to predicting clinical outcomes in urologic cancer

MOTIVATION

To better understand the molecular features of cancers, a comprehensive analysis using multiomics data has been conducted. Additionally, a pathway activity inference method has been developed to facilitate the integrative effects of multiple genes. In this respect, we have recently proposed a novel integrative pathway activity inference approach, iDRW, and demonstrated the effectiveness of the method with respect to dichotomizing two survival groups. However, there were several limitations, such as a lack of generality. In this study, we designed a directed gene-gene graph using pathway information by assigning interactions between genes in multiple layers of networks.

RESULTS

: As a proof-of-concept study, it was evaluated using three genomic profiles of urologic cancer patients. The proposed integrative approach achieved improved outcome prediction performances compared with a single genomic profile alone and other existing pathway activity inference methods. The integrative approach also identified common/cancer-specific candidate driver pathways as predictive prognostic features in urologic cancers. Furthermore, it provides better biological insights into the prioritized pathways and genes in an integrated view using a multi-layered gene-gene network. Our framework is not specifically designed for urologic cancers and can be generally applicable for various datasets.

AVAILABILITY

iDRW is implemented as the R software package. The source codes are available at https://github.com/sykim122/iDRW.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

MultiPaths: a python framework for analyzing multi-layer biological networks using diffusion algorithms

Summary

High-throughput screening yields vast amounts of biological data which can be highly challenging to interpret. In response, knowledge-driven approaches emerged as possible solutions to analyze large datasets by leveraging prior knowledge of biomolecular interactions represented in the form of biological networks. Nonetheless, given their size and complexity, their manual investigation quickly becomes impractical. Thus, computational approaches, such as diffusion algorithms, are often employed to interpret and contextualize the results of high-throughput experiments. Here, we present MultiPaths, a framework consisting of two independent Python packages for network analysis. While the first package, DiffuPy, comprises numerous commonly used diffusion algorithms applicable to any generic network, the second, DiffuPath, enables the application of these algorithms on multi-layer biological networks. To facilitate its usability, the framework includes a command line interface, reproducible examples and documentation. To demonstrate the framework, we conducted several diffusion experiments on three independent multi-omics datasets over disparate networks generated from pathway databases, thus, highlighting the ability of multi-layer networks to integrate multiple modalities. Finally, the results of these experiments demonstrate how the generation of harmonized networks from disparate databases can improve predictive performance with respect to individual resources.

Availability and implementation

DiffuPy and DiffuPath are publicly available under the Apache License 2.0 at https://github.com/multipaths.

Supplementary information

Supplementary data are available at Bioinformatics online.