A Fast and Flexible Framework for Network-Assisted Genomic Association

State of the interactomes: an evaluation of molecular networks for generating biological insights

Advancements in genomic and proteomic technologies have powered the creation of large gene and protein networks ("interactomes") for understanding biological systems. However, the proliferation of interactomes complicates the selection of networks for specific applications. Here, we present a comprehensive evaluation of 45 current human interactomes, encompassing protein-protein interactions as well as gene regulatory, signaling, colocalization, and genetic interaction networks. Our analysis shows that large composite networks such as HumanNet, STRING, and FunCoup are most effective for identifying disease genes, while smaller networks such as DIP, Reactome, and SIGNOR demonstrate stronger performance in interaction prediction. Our study provides a benchmark for interactomes across diverse biological applications and clarifies factors that influence network performance. Furthermore, our evaluation pipeline paves the way for continued assessment of emerging and updated interaction networks in the future.

Analyses of GWAS signal using GRIN identify additional genes contributing to suicidal behavior

Genome-wide association studies (GWAS) identify genetic variants underlying complex traits but are limited by stringent genome-wide significance thresholds. We present GRIN (Gene set Refinement through Interacting Networks), which increases confidence in the expanded gene set by retaining genes strongly connected by biological networks when GWAS thresholds are relaxed. GRIN was validated on both simulated interrelated gene sets as well as multiple GWAS traits. From multiple GWAS summary statistics of suicide attempt, a complex phenotype, GRIN identified additional genes that replicated across independent cohorts and retained biologically interrelated genes despite a relaxed significance threshold. We present a conceptual model of how these retained genes interact through neurobiological pathways that may influence suicidal behavior, and identify existing drugs associated with these pathways that would not have been identified under traditional GWAS thresholds. We demonstrate GRIN's utility in boosting GWAS results by increasing the number of true positive genes identified from GWAS results.

Current and future directions in network biology

Summary

Network biology is an interdisciplinary field bridging computational and biological sciences that has proved pivotal in advancing the understanding of cellular functions and diseases across biological systems and scales. Although the field has been around for two decades, it remains nascent. It has witnessed rapid evolution, accompanied by emerging challenges. These stem from various factors, notably the growing complexity and volume of data together with the increased diversity of data types describing different tiers of biological organization. We discuss prevailing research directions in network biology, focusing on molecular/cellular networks but also on other biological network types such as biomedical knowledge graphs, patient similarity networks, brain networks, and social/contact networks relevant to disease spread. In more detail, we highlight areas of inference and comparison of biological networks, multimodal data integration and heterogeneous networks, higher-order network analysis, machine learning on networks, and network-based personalized medicine. Following the overview of recent breakthroughs across these five areas, we offer a perspective on future directions of network biology. Additionally, we discuss scientific communities, educational initiatives, and the importance of fostering diversity within the field. This article establishes a roadmap for an immediate and long-term vision for network biology.

Availability and implementation

Not applicable.

State of the Interactomes: an evaluation of molecular networks for generating biological insights

Advancements in genomic and proteomic technologies have powered the use of gene and protein networks ("interactomes") for understanding genotype-phenotype translation. However, the proliferation of interactomes complicates the selection of networks for specific applications. Here, we present a comprehensive evaluation of 46 current human interactomes, encompassing protein-protein interactions as well as gene regulatory, signaling, colocalization, and genetic interaction networks. Our analysis shows that large composite networks such as HumanNet, STRING, and FunCoup are most effective for identifying disease genes, while smaller networks such as DIP and SIGNOR demonstrate strong interaction prediction performance. These findings provide a benchmark for interactomes across diverse network biology applications and clarify factors that influence network performance. Furthermore, our evaluation pipeline paves the way for continued assessment of emerging and updated interaction networks in the future.

Lithium response in bipolar disorder is associated with focal adhesion and PI3K-Akt networks: a multi-omics replication study

Lithium is the gold standard treatment for bipolar disorder (BD). However, its mechanism of action is incompletely understood, and prediction of treatment outcomes is limited. In our previous multi-omics study of the Pharmacogenomics of Bipolar Disorder (PGBD) sample combining transcriptomic and genomic data, we found that focal adhesion, the extracellular matrix (ECM), and PI3K-Akt signaling networks were associated with response to lithium. In this study, we replicated the results of our previous study using network propagation methods in a genome-wide association study of an independent sample of 2039 patients from the International Consortium on Lithium Genetics (ConLiGen) study. We identified functional enrichment in focal adhesion and PI3K-Akt pathways, but we did not find an association with the ECM pathway. Our results suggest that deficits in the neuronal growth cone and PI3K-Akt signaling, but not in ECM proteins, may influence response to lithium in BD.

Lithium Response in Bipolar Disorder is Associated with Focal Adhesion and PI3K-Akt Networks: A Multi-omics Replication Study

Lithium is the gold standard treatment for bipolar disorder (BD). However, its mechanism of action is incompletely understood, and prediction of treatment outcomes is limited. In our previous multi-omics study of the Pharmacogenomics of Bipolar Disorder (PGBD) sample combining transcriptomic and genomic data, we found that focal adhesion, the extracellular matrix (ECM), and PI3K-Akt signaling networks were associated with response to lithium. In this study, we replicated the results of our previous study using network propagation methods in a genome-wide association study of an independent sample of 2,039 patients from the International Consortium on Lithium Genetics (ConLiGen) study. We identified functional enrichment in focal adhesion and PI3K-Akt pathways, but we did not find an association with the ECM pathway. Our results suggest that deficits in the neuronal growth cone and PI3K-Akt signaling, but not in ECM proteins, may influence response to lithium in BD.

Genome-wide association studies of human and rat BMI converge on synapse, epigenome, and hormone signaling networks

A vexing observation in genome-wide association studies (GWASs) is that parallel analyses in different species may not identify orthologous genes. Here, we demonstrate that cross-species translation of GWASs can be greatly improved by an analysis of co-localization within molecular networks. Using body mass index (BMI) as an example, we show that the genes associated with BMI in humans lack significant agreement with those identified in rats. However, the networks interconnecting these genes show substantial overlap, highlighting common mechanisms including synaptic signaling, epigenetic modification, and hormonal regulation. Genetic perturbations within these networks cause abnormal BMI phenotypes in mice, too, supporting their broad conservation across mammals. Other mechanisms appear species specific, including carbohydrate biosynthesis (humans) and glycerolipid metabolism (rodents). Finally, network co-localization also identifies cross-species convergence for height/body length. This study advances a general paradigm for determining whether and how phenotypes measured in model species recapitulate human biology.

Genome-wide association studies in plant pathosystems: success or failure?

Harnessing natural genetic variation is an established alternative to artificial genetic variation for investigating the molecular dialog between partners in plant pathosystems. Herein, we review the successes of genome-wide association studies (GWAS) in both plants and pathogens. While GWAS in plants confirmed that the genetic architecture of disease resistance is polygenic, dynamic during the infection kinetics, and dependent on the environment, GWAS shortened the time of identification of quantitative trait loci (QTLs) and revealed both complex epistatic networks and a genetic architecture dependent upon the geographical scale. A similar picture emerges from the few GWAS in pathogens. In addition, the ever-increasing number of functionally validated QTLs has revealed new molecular plant defense mechanisms and pathogenicity determinants. Finally, we propose recommendations to better decode the disease triangle.

Interaction-integrated linear mixed model reveals 3D-genetic basis underlying Autism

Genetic interactions play critical roles in genotype-phenotype associations. We developed a novel interaction-integrated linear mixed model (ILMM) that integrates a priori knowledge into linear mixed models. ILMM enables statistical integration of genetic interactions upfront and overcomes the problems of searching for combinations. To demonstrate its utility, with 3D genomic interactions (assessed by Hi-C experiments) as a priori, we applied ILMM to whole-genome sequencing data for Autism Spectrum Disorders (ASD) and brain transcriptome data, revealing the 3D-genetic basis of ASD and 3D-expression quantitative loci (3D-eQTLs) for brain tissues. Notably, we reported a potential mechanism involving distal regulation between FOXP2 and DNMT3A, conferring the risk of ASD.

Self-supervised graph representation learning integrates multiple molecular networks and decodes gene-disease relationships

Leveraging molecular networks to discover disease-relevant modules is a long-standing challenge. With the accumulation of interactomes, there is a pressing need for powerful computational approaches to handle the inevitable noise and context-specific nature of biological networks. Here, we introduce Graphene, a two-step self-supervised representation learning framework tailored to concisely integrate multiple molecular networks and adapted to gene functional analysis via downstream re-training. In practice, we first leverage GNN (graph neural network) pre-training techniques to obtain initial node embeddings followed by re-training Graphene using a graph attention architecture, achieving superior performance over competing methods for pathway gene recovery, disease gene reprioritization, and comorbidity prediction. Graphene successfully recapitulates tissue-specific gene expression across disease spectrum and demonstrates shared heritability of common mental disorders. Graphene can be updated with new interactomes or other omics features. Graphene holds promise to decipher gene function under network context and refine GWAS (genome-wide association study) hits and offers mechanistic insights via decoding diseases from genome to networks to phenotypes.

NetMix2: A Principled Network Propagation Algorithm for Identifying Altered Subnetworks

A standard paradigm in computational biology is to leverage interaction networks as prior knowledge in analyzing high-throughput biological data, where the data give a score for each vertex in the network. One classical approach is the identification of altered subnetworks, or subnetworks of the interaction network that have both outlier vertex scores and a defined network topology. One class of algorithms for identifying altered subnetworks search for high-scoring subnetworks in subnetwork families with simple topological constraints, such as connected subnetworks, and have sound statistical guarantees. A second class of algorithms employ network propagation-the smoothing of vertex scores over the network using a random walk or diffusion process-and utilize the global structure of the network. However, network propagation algorithms often rely on ad hoc heuristics that lack a rigorous statistical foundation. In this work, we unify the subnetwork family and network propagation approaches by deriving the propagation family, a subnetwork family that approximates the sets of vertices ranked highly by network propagation approaches. We introduce NetMix2, a principled algorithm for identifying altered subnetworks from a wide range of subnetwork families. When using the propagation family, NetMix2 combines the advantages of the subnetwork family and network propagation approaches. NetMix2 outperforms other methods, including network propagation on simulated data, pan-cancer somatic mutation data, and genome-wide association data from multiple human diseases.

Genetically personalised organ-specific metabolic models in health and disease

Understanding how genetic variants influence disease risk and complex traits (variant-to-function) is one of the major challenges in human genetics. Here we present a model-driven framework to leverage human genome-scale metabolic networks to define how genetic variants affect biochemical reaction fluxes across major human tissues, including skeletal muscle, adipose, liver, brain and heart. As proof of concept, we build personalised organ-specific metabolic flux models for 524,615 individuals of the INTERVAL and UK Biobank cohorts and perform a fluxome-wide association study (FWAS) to identify 4312 associations between personalised flux values and the concentration of metabolites in blood. Furthermore, we apply FWAS to identify 92 metabolic fluxes associated with the risk of developing coronary artery disease, many of which are linked to processes previously described to play in role in the disease. Our work demonstrates that genetically personalised metabolic models can elucidate the downstream effects of genetic variants on biochemical reactions involved in common human diseases.

mGWAS-Explorer: Linking SNPs, Genes, Metabolites, and Diseases for Functional Insights

Tens of thousands of single-nucleotide polymorphisms (SNPs) have been identified to be significantly associated with metabolite abundance in over 65 genome-wide association studies with metabolomics (mGWAS) to date. Obtaining mechanistic or functional insights from these associations for translational applications has become a key research area in the mGWAS community. Here, we introduce mGWAS-Explorer, a user-friendly web-based platform to help connect SNPs, metabolites, genes, and their known disease associations via powerful network visual analytics. The application of the mGWAS-Explorer was demonstrated using a COVID-19 and a type 2 diabetes case studies.

Associations of risk genes with onset age and plasma biomarkers of Alzheimer's disease: a large case-control study in mainland China

Most genetic studies concerning risk genes in Alzheimer's disease (AD) are from Caucasian populations, whereas the data remain limited in the Chinese population. In this study, we systematically explored the relationship between AD and risk genes in mainland China. We sequenced 33 risk genes previously reported to be associated with AD in a total of 3604 individuals in the mainland Chinese population. Common variant (MAF ≥ 0.01) based association analysis and gene-based (MAF < 0.01) association test were performed by PLINK 1.9 and Sequence Kernel Association Test-Optimal, respectively. Polygenic risk score (PRS) was calculated, and receiver operating characteristic curve (AUC) was computed. Plasma Aβ42, Aβ40, total tau (T-tau), and neurofilament light chain (NFL) were tested in a subgroup, and their associations with PRS were conducted using the Spearman correlation test. Six common variants varied significantly between AD patients and cognitively normal controls after the adjustment of age, gender, and APOE ε4 status, including variants in ABCA7 (n = 5) and APOE (n = 1). Among them, four variants were novel and two were reported previously. The AUC of PRS was 0.71. The high PRS was significantly associated with an earlier age at onset (P = 4.30 × 10-4). PRS was correlated with plasma Aβ42, Aβ42/Aβ40 ratio, T-tau, and NFL levels. Gene-based association test revealed that ABCA7 and UNC5C reached statistical significance. The common variants in APOE and ABCA7, as well as rare variants in ABCA7 and UNC5C, may contribute to the etiology of AD. Moreover, the PRS, to some extent, could predict the risk, onset age, and biological changes of AD.

Evaluation of Safety-Related Outcomes of One-Segment and More-Than-One-Segment High-Level Hepatectomy in Hepatocellular Carcinoma Based on the Japanese Board Certification System

Background

We evaluated the impact of the Japanese board certification system for expert surgeons (JBCSES) on complications and survival outcomes in hepatectomy for hepatocellular carcinoma.

Methods

The postoperative outcomes of 493 patients who underwent high-level liver surgery involving one-segment (OSeg) hepatectomy and more-than-one-segment (MOSeg) resection were compared before and after JBCSES establishment. After the establishment of the JBCSES, the patients’ postoperative outcomes were compared using propensity score matching (PSM) to determine the influence of expert surgeons.

Results

The establishment of the JBCSES was associated with a decrease in the overall postoperative complication rates after high-level liver surgery from 50.2 to 38.1% (P = 0.008) and a decrease in Clavien–Dindo class ≥ IIIb complications from 10.2 to 5.0% (P = 0.035). The 90-day mortality rate decreased from 5.1 to 0.7% (P = 0.003), and the 5-year survival rate increased from 51.4 to 63.9% (P = 0.009). Using PSM, a comparison of OSeg hepatectomies that involved expert surgeons (n = 48) and those that did not (n = 48) showed significantly lower intraoperative blood loss in surgeries involving an expert surgeon (mean, 340 vs. 473 mL; P = 0.033). There were no significant differences in complication rates or long-term prognosis between these groups. A comparison of MOSeg hepatectomies that involved expert surgeons (n = 26) and those that did not (n = 26) showed no significant difference in surgical factors, complications, or overall survival between the two groups.

Conclusions

After establishment of the JBCSES, postoperative complication rates and mortality rates decreased and survival rates increased following liver surgery. Expert surgeon participation significantly decreased intraoperative blood loss during OSeg hepatectomies.

NDEx: Accessing Network Models and Streamlining Network Biology Workflows

NDEx, the Network Data Exchange (https://www.ndexbio.org) is a web-based resource where users can find, store, share and publish network models of any type and size. NDEx is integrated with Cytoscape, the widely used desktop application for network analysis and visualization. NDEx and Cytoscape are the pillars of the Cytoscape Ecosystem, a diverse environment of resources, tools, applications and services for network biology workflows. In this article, we introduce researchers to NDEx and highlight how it can simplify common tasks in network biology workflows as well as streamline publication and access to). Finally, we show how NDEx can be used programmatically via Python with the 'ndex2' client library, and point readers to additional examples for other popular programming languages such as JavaScript and R. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Getting started with NDEx Basic Protocol 2: Using NDEx and Cytoscape in a publication-oriented workflow Basic Protocol 3: Manipulating networks in NDEx via Python.

Network and pathway expansion of genetic disease associations identifies successful drug targets

Genetic evidence of disease association has often been used as a basis for selecting of drug targets for complex common diseases. Likewise, the propagation of genetic evidence through gene or protein interaction networks has been shown to accurately infer novel disease associations at genes for which no direct genetic evidence can be observed. However, an empirical test of the utility of combining these approaches for drug discovery has been lacking. In this study, we examine genetic associations arising from an analysis of 648 UK Biobank GWAS and evaluate whether targets identified as proxies of direct genetic hits are enriched for successful drug targets, as measured by historical clinical trial data. We find that protein networks formed from specific functional linkages such as protein complexes and ligand–receptor pairs are suitable for even naïve guilt-by-association network propagation approaches. In addition, more sophisticated approaches applied to global protein–protein interaction networks and pathway databases, also successfully retrieve targets enriched for clinically successful drug targets. We conclude that network propagation of genetic evidence can be used for drug target identification.

FI-Net: Identification of Cancer Driver Genes by Using Functional Impact Prediction Neural Network

Identification of driver genes, whose mutations cause the development of tumors, is crucial for the improvement of cancer research and precision medicine. To overcome the problem that the traditional frequency-based methods cannot detect lowly recurrently mutated driver genes, researchers have focused on the functional impact of gene mutations and proposed the function-based methods. However, most of the function-based methods estimate the distribution of the null model through the non-parametric method, which is sensitive to sample size. Besides, such methods could probably lead to underselection or overselection results. In this study, we proposed a method to identify driver genes by using functional impact prediction neural network (FI-net). An artificial neural network as a parametric model was constructed to estimate the functional impact scores for genes, in which multi-omics features were used as the multivariate inputs. Then the estimation of the background distribution and the identification of driver genes were conducted in each cluster obtained by the hierarchical clustering algorithm. We applied FI-net and other 22 state-of-the-art methods to 31 datasets from The Cancer Genome Atlas project. According to the comprehensive evaluation criterion, FI-net was powerful among various datasets and outperformed the other methods in terms of the overlap fraction with Cancer Gene Census and Network of Cancer Genes database, and the consensus in predictions among methods. Furthermore, the results illustrated that FI-net can identify known and potential novel driver genes.

NetCore: a network propagation approach using node coreness

We present NetCore, a novel network propagation approach based on node coreness, for phenotype–genotype associations and module identification. NetCore addresses the node degree bias in PPI networks by using node coreness in the random walk with restart procedure, and achieves improved re-ranking of genes after propagation. Furthermore, NetCore implements a semi-supervised approach to identify phenotype-associated network modules, which anchors the identification of novel candidate genes at known genes associated with the phenotype. We evaluated NetCore on gene sets from 11 different GWAS traits and showed improved performance compared to the standard degree-based network propagation using cross-validation. Furthermore, we applied NetCore to identify disease genes and modules for Schizophrenia GWAS data and pan-cancer mutation data. We compared the novel approach to existing network propagation approaches and showed the benefits of using NetCore in comparison to those. We provide an easy-to-use implementation, together with a high confidence PPI network extracted from ConsensusPathDB, which can be applied to various types of genomics data in order to obtain a re-ranking of genes and functionally relevant network modules.

uKIN Combines New and Prior Information with Guided Network Propagation to Accurately Identify Disease Genes

Protein interaction networks provide a powerful framework for identifying genes causal for complex genetic diseases. Here, we introduce a general framework, uKIN, that uses prior knowledge of disease-associated genes to guide, within known protein-protein interaction networks, random walks that are initiated from newly identified candidate genes. In large-scale testing across 24 cancer types, we demonstrate that our network propagation approach for integrating both prior and new information not only better identifies cancer driver genes than using either source of information alone but also readily outperforms other state-of-the-art network-based approaches. We also apply our approach to genome-wide association data to identify genes functionally relevant for several complex diseases. Overall, our work suggests that guided network propagation approaches that utilize both prior and new data are a powerful means to identify disease genes. uKIN is freely available for download at: https://github.com/Singh-Lab/uKIN.

Strategies for Network GWAS Evaluated Using Classroom Crowd Science

Biological networks can substantially boost power to identify disease genes in genome-wide association studies. To explore different network GWAS methods, we challenged students of a UC San Diego graduate level bioinformatics course, Network Biology and Biomedicine, to explore and improve such algorithms during a four-week-long classroom competition. Here, we report the many creative solutions and share our experiences in conducting classroom crowd science as both a research and pedagogical tool.