Network inference in systems biology: recent developments, challenges, and applications

Camellia Japonica Radix modulates gut microbiota and 9(S)-HpODE-mediated ferroptosis to alleviate oxidative stress against MASLD

BACKGROUND

Camellia japonica radix (CJR), derived from the root of Camellia japonica L., has the potential to function as an herbal tea substitute for the prevention and intervention of metabolic dysfunction-associated steatotic liver disease (MASLD). It can provide systemic therapeutic benefits, boast a favorable safety profile, facilitate convenient consumption, and support long-term applicability. Despite its potential, research on CJR remains limited.

PURPOSE

The aim of this study aims is to elucidate the therapeutic mechanisms of CJR in MASLD, thereby providing evidence to support its clinical application.

METHODS

The therapeutic effects of CJR were evaluated using a water-supplementation model in MASLD mice. Integrated microbiome, transcriptome, proteome, and metabolome analyses were employed to comprehensively explore the mechanisms involved. A drug-target pull-down assay was performed to identify specific protein targets of small molecule metabolites in vitro. Fecal microbiota transplantation in antibiotic-treated ABX mice was conducted to confirm the critical role of gut microbiota and its metabolites. Furthermore, customized medicated feed supplemented with linoleic acid was used to explore the intervention effect of its metabolite, 9(S)-HpODE, as well as to evaluate its dietary intervention potential.

RESULTS

This present study explicitly elucidates the efficacy of CJR extract in alleviating hepatic inflammation and steatosis in a MASLD model mice, with its pharmacological mechanism associated with gut microbiota, linoleic acid metabolism, and GPX4-mediated ferroptosis. Notably, 9(S)-HpODE was discovered to be a key metabolite of linoleic acid, which could target both KEAP1 and SLC7A11, bidirectionally regulating GPX4-mediated ferroptosis, while acting as a signaling molecule at low doses to induce redox adaptation via oxidative preconditioning, thus ameliorating oxidative stress in MASLD.

CONCLUSION

Our findings indicate that both CJR and linoleic acid exhibit significant potential as dietary interventions for the management of MASLD, offering promising avenues for future research and clinical application.

The Evolution of Systems Biology and Systems Medicine: From Mechanistic Models to Uncertainty Quantification

Understanding interaction mechanisms within cells, tissues, and organisms is crucial for driving developments across biology and medicine. Mathematical modeling is an essential tool for simulating such biological systems. Building on experiments, mechanistic models are widely used to describe small-scale intracellular networks. The development of sequencing techniques and computational tools has recently enabled multiscale models. Combining such larger scale network modeling with mechanistic modeling provides us with an opportunity to reveal previously unknown disease mechanisms and pharmacological interventions. Here, we review systems biology models from mechanistic models to multiscale models that integrate multiple layers of cellular networks and discuss how they can be used to shed light on disease states and even wellness-related states. Additionally, we introduce several methods that increase the certainty and accuracy of model predictions. Thus, combining mechanistic models with emerging mathematical and computational techniques can provide us with increasingly powerful tools to understand disease states and inspire drug discoveries.

Unifying fragmented perspectives with additive deep learning for high-dimensional models from partial faceted datasets

Biological systems are complex networks where measurable functions emerge from interactions among thousands of components. Many studies aim to link biological function with molecular elements, yet quantifying their contributions simultaneously remains challenging, especially at the single-cell level. We propose a machine-learning approach that integrates faceted data subsets to reconstruct a complete view of the system using conditional distributions. We develop both polynomial regression and neural network models, validated with two examples: a mechanical spring network under external forces and an 8-dimensional biological network involving the senescence marker P53, using single-cell data. Our results demonstrate successful system reconstruction from partial datasets, with predictive accuracy improving as more variables are measured. This approach offers a systematic method to integrate fragmented experimental data, enabling unbiased and holistic modeling of complex biological functions.

Multi-omics in exploring the pathophysiology of diabetic retinopathy

Diabetic retinopathy (DR) is a leading global cause of vision impairment, with its prevalence increasing alongside the rising rates of diabetes mellitus (DM). Despite the retina's complex structure, the underlying pathology of DR remains incompletely understood. Single-cell RNA sequencing (scRNA-seq) and recent advancements in multi-omics analyses have revolutionized molecular profiling, enabling high-throughput analysis and comprehensive characterization of complex biological systems. This review highlights the significant contributions of scRNA-seq, in conjunction with other multi-omics technologies, to DR research. Integrated scRNA-seq and transcriptomic analyses have revealed novel insights into DR pathogenesis, including alternative transcription start site events, fluctuations in cell populations, altered gene expression profiles, and critical signaling pathways within retinal cells. Furthermore, by integrating scRNA-seq with genetic association studies and multi-omics analyses, researchers have identified novel biomarkers, susceptibility genes, and potential therapeutic targets for DR, emphasizing the importance of specific retinal cell types in disease progression. The integration of scRNA-seq with metabolomics has also been instrumental in identifying specific metabolites and dysregulated pathways associated with DR. It is highly conceivable that the continued synergy between scRNA-seq and other multi-omics approaches will accelerate the discovery of underlying mechanisms and the development of novel therapeutic interventions for DR.

Unravelling the genetic basis and regulation networks related to fibre quality improvement using chromosome segment substitution lines in cotton

The elucidation of genetic architecture and molecular regulatory networks underlying complex traits remains a significant challenge in life science, largely due to the substantial background effects that arise from epistasis and gene-environment interactions. The chromosome segment substitution line (CSSL) is an ideal material for genetic and molecular dissection of complex traits due to its near-isogenic properties; yet a comprehensive analysis, from the basic identification of substitution segments to advanced regulatory network, is still insufficient. Here, we developed two cotton CSSL populations on the Gossypium hirsutum background, representing wide adaptation and high lint yield, with introgression from G. barbadense, representing superior fibre quality. We sequenced 99 CSSLs that demonstrated significant differences from G. hirsutum in fibre, and characterized 836 dynamic fibre transcriptomes in three crucial developmental stages. We developed a workflow for precise resolution of chromosomal substitution segments; the genome sequencing revealed substitutions collectively representing 87.25% of the G. barbadense genome. Together, the genomic and transcriptomic survey identified 18 novel fibre-quality-related quantitative trait loci with high genetic contributions and the comprehensive landscape of fibre development regulation. Furthermore, analysis determined unique cis-expression patterns in CSSLs to be the driving force for fibre quality alteration; building upon this, the co-expression regulatory network revealed biological relationships among the noted pathways and accurately described the molecular interactions of GhHOX3, GhRDL1 and GhEXPA1 during fibre elongation, along with reliable predictions for their interactions with GhTBA8A5. Our study will enhance more strategic employment of CSSL in crop molecular biology and breeding programmes.

k- Strong Inference Algorithm: A Hybrid Information Theory Based Gene Network Inference Algorithm

Gene networks allow researchers to understand the underlying mechanisms between diseases and genes while reducing the need for wet lab experiments. Numerous gene network inference (GNI) algorithms have been presented in the literature to infer accurate gene networks. We proposed a hybrid GNI algorithm, k-Strong Inference Algorithm (ksia), to infer more reliable and robust gene networks from omics datasets. To increase reliability, ksia integrates Pearson correlation coefficient (PCC) and Spearman rank correlation coefficient (SCC) scores to determine mutual information scores between molecules to increase diversity of relation predictions. To infer a more robust gene network, ksia applies three different elimination steps to remove redundant and spurious relations between genes. The performance of ksia was evaluated on microbe microarrays database in the overlap analysis with other GNI algorithms, namely ARACNE, C3NET, CLR, and MRNET. Ksia inferred less number of relations due to its strict elimination steps. However, ksia generally performed better on Escherichia coli (E.coli) and Saccharomyces cerevisiae (yeast) gene expression datasets due to F- measure and precision values. The integration of association estimator scores and three elimination stages slightly increases the performance of ksia based gene networks. Users can access ksia R package and user manual of package via https://github.com/ozgurcingiz/ksia .

EIEPCF: accurate inference of functional gene regulatory networks by eliminating indirect effects from confounding factors

Reconstructing functional gene regulatory networks (GRNs) is a primary prerequisite for understanding pathogenic mechanisms and curing diseases in animals, and it also provides an important foundation for cultivating vegetable and fruit varieties that are resistant to diseases and corrosion in plants. Many computational methods have been developed to infer GRNs, but most of the regulatory relationships between genes obtained by these methods are biased. Eliminating indirect effects in GRNs remains a significant challenge for researchers. In this work, we propose a novel approach for inferring functional GRNs, named EIEPCF (eliminating indirect effects produced by confounding factors), which eliminates indirect effects caused by confounding factors. This method eliminates the influence of confounding factors on regulatory factors and target genes by measuring the similarity between their residuals. The validation results of the EIEPCF method on simulation studies, the gold-standard networks provided by the DREAM3 Challenge and the real gene networks of Escherichia coli demonstrate that it achieves significantly higher accuracy compared to other popular computational methods for inferring GRNs. As a case study, we utilized the EIEPCF method to reconstruct the cold-resistant specific GRN from gene expression data of cold-resistant in Arabidopsis thaliana. The source code and data are available at https://github.com/zhanglab-wbgcas/EIEPCF.

Exploring gene regulation and biological processes in insects: Insights from omics data using gene regulatory network models

Gene regulatory network (GRN) comprises complicated yet intertwined gene-regulator relationships. Understanding the GRN dynamics will unravel the complexity behind the observed gene expressions. Insect gene regulation is often complicated due to their complex life cycles and diverse ecological adaptations. The main interest of this review is to have an update on the current mathematical modelling methods of GRNs to explain insect science. Several popular GRN architecture models are discussed, together with examples of applications in insect science. In the last part of this review, each model is compared from different aspects, including network scalability, computation complexity, robustness to noise and biological relevancy.

Inference of gene regulatory networks based on directed graph convolutional networks

Inferring gene regulatory network (GRN) is one of the important challenges in systems biology, and many outstanding computational methods have been proposed; however there remains some challenges especially in real datasets. In this study, we propose Directed Graph Convolutional neural network-based method for GRN inference (DGCGRN). To better understand and process the directed graph structure data of GRN, a directed graph convolutional neural network is conducted which retains the structural information of the directed graph while also making full use of neighbor node features. The local augmentation strategy is adopted in graph neural network to solve the problem of poor prediction accuracy caused by a large number of low-degree nodes in GRN. In addition, for real data such as E.coli, sequence features are obtained by extracting hidden features using Bi-GRU and calculating the statistical physicochemical characteristics of gene sequence. At the training stage, a dynamic update strategy is used to convert the obtained edge prediction scores into edge weights to guide the subsequent training process of the model. The results on synthetic benchmark datasets and real datasets show that the prediction performance of DGCGRN is significantly better than existing models. Furthermore, the case studies on bladder uroepithelial carcinoma and lung cancer cells also illustrate the performance of the proposed model.

PathFinder: a novel graph transformer model to infer multi-cell intra- and inter-cellular signaling pathways and communications

Recently, large-scale scRNA-seq datasets have been generated to understand the complex signaling mechanisms within the microenvironment of Alzheimer's Disease (AD), which are critical for identifying novel therapeutic targets and precision medicine. However, the background signaling networks are highly complex and interactive. It remains challenging to infer the core intra- and inter-multi-cell signaling communication networks using scRNA-seq data. In this study, we introduced a novel graph transformer model, PathFinder, to infer multi-cell intra- and inter-cellular signaling pathways and communications among multi-cell types. Compared with existing models, the novel and unique design of PathFinder is based on the divide-and-conquer strategy. This model divides complex signaling networks into signaling paths, which are then scored and ranked using a novel graph transformer architecture to infer intra- and inter-cell signaling communications. We evaluated the performance of PathFinder using two scRNA-seq data cohorts. The first cohort is an APOE4 genotype-specific AD, and the second is a human cirrhosis cohort. The evaluation confirms the promising potential of using PathFinder as a general signaling network inference model.

Topological benchmarking of algorithms to infer gene regulatory networks from single-cell RNA-seq data

MOTIVATION

In recent years, many algorithms for inferring gene regulatory networks from single-cell transcriptomic data have been published. Several studies have evaluated their accuracy in estimating the presence of an interaction between pairs of genes. However, these benchmarking analyses do not quantify the algorithms' ability to capture structural properties of networks, which are fundamental, e.g., for studying the robustness of a gene network to external perturbations. Here, we devise a three-step benchmarking pipeline called STREAMLINE that quantifies the ability of algorithms to capture topological properties of networks and identify hubs.

RESULTS

To this aim, we use data simulated from different types of networks as well as experimental data from three different organisms. We apply our benchmarking pipeline to four inference algorithms and provide guidance on which algorithm should be used depending on the global network property of interest.

AVAILABILITY AND IMPLEMENTATION

STREAMLINE is available at https://github.com/ScialdoneLab/STREAMLINE. The data generated in this study are available at https://doi.org/10.5281/zenodo.10710444.

Organ-delimited gene regulatory networks provide high accuracy in candidate transcription factor selection across diverse processes

Organ-specific gene expression datasets that include hundreds to thousands of experiments allow the reconstruction of organ-level gene regulatory networks (GRNs). However, creating such datasets is greatly hampered by the requirements of extensive and tedious manual curation. Here, we trained a supervised classification model that can accurately classify the organ-of-origin for a plant transcriptome. This K-Nearest Neighbor-based multiclass classifier was used to create organ-specific gene expression datasets for the leaf, root, shoot, flower, and seed in Arabidopsis thaliana. A GRN inference approach was used to determine the: i. influential transcription factors (TFs) in each organ and, ii. most influential TFs for specific biological processes in that organ. These genome-wide, organ-delimited GRNs (OD-GRNs), recalled many known regulators of organ development and processes operating in those organs. Importantly, many previously unknown TF regulators were uncovered as potential regulators of these processes. As a proof-of-concept, we focused on experimentally validating the predicted TF regulators of lipid biosynthesis in seeds, an important food and biofuel trait. Of the top 20 predicted TFs, eight are known regulators of seed oil content, e.g., WRI1, LEC1, FUS3. Importantly, we validated our prediction of MybS2, TGA4, SPL12, AGL18, and DiV2 as regulators of seed lipid biosynthesis. We elucidated the molecular mechanism of MybS2 and show that it induces purple acid phosphatase family genes and lipid synthesis genes to enhance seed lipid content. This general approach has the potential to be extended to any species with sufficiently large gene expression datasets to find unique regulators of any trait-of-interest.

From biological data to oscillator models using SINDy

Periodic changes in the concentration or activity of different molecules regulate vital cellular processes such as cell division and circadian rhythms. Developing mathematical models is essential to better understand the mechanisms underlying these oscillations. Recent data-driven methods like SINDy have fundamentally changed model identification, yet their application to experimental biological data remains limited. This study investigates SINDy's constraints by directly applying it to biological oscillatory data. We identify insufficient resolution, noise, dimensionality, and limited prior knowledge as primary limitations. Using various generic oscillator models of different complexity and/or dimensionality, we systematically analyze these factors. We then propose a comprehensive guide for inferring models from biological data, addressing these challenges step by step. Our approach is validated using glycolytic oscillation data from yeast.

PathFinder: a novel graph transformer model to infer multi-cell intra- and inter-cellular signaling pathways and communications

Recently, large-scale scRNA-seq datasets have been generated to understand the complex and poorly understood signaling mechanisms within microenvironment of Alzheimer's Disease (AD), which are critical for identifying novel therapeutic targets and precision medicine. Though a set of targets have been identified, however, it remains a challenging to infer the core intra- and inter-multi-cell signaling communication networks using the scRNA-seq data, considering the complex and highly interactive background signaling network. Herein, we introduced a novel graph transformer model, PathFinder, to infer multi-cell intra- and inter-cellular signaling pathways and signaling communications among multi-cell types. Compared with existing models, the novel and unique design of PathFinder is based on the divide-and-conquer strategy, which divides the complex signaling networks into signaling paths, and then score and rank them using a novel graph transformer architecture to infer the intra- and inter-cell signaling communications. We evaluated PathFinder using scRNA-seq data of APOE4-genotype specific AD mice models and identified novel APOE4 altered intra- and inter-cell interaction networks among neurons, astrocytes, and microglia. PathFinder is a general signaling network inference model and can be applied to other omics data-driven signaling network inference.

Virus-host protein co-expression networks reveal temporal organization and strategies of viral infection

Viral replication is a complex dynamical process involving the global remodeling of the host cellular machinery across several stages. In this study, we provide a unified view of the virus-host interaction at the proteome level reconstructing protein co-expression networks from quantitative temporal data of four large DNA viruses. We take advantage of a formal framework, the theory of competing networks, to describe the viral infection as a dynamical system taking place on a network of networks where perturbations induced by viral proteins spread to hijack the host proteome for the virus benefit. Our methodology demonstrates how the viral replication cycle can be effectively examined as a complex interaction between protein networks, providing useful insights into the viral and host's temporal organization and strategies, key protein nodes targeted by the virus and dynamical bottlenecks during the course of the infection.

Graph embedding and geometric deep learning relevance to network biology and structural chemistry

Graphs are used as a model of complex relationships among data in biological science since the advent of systems biology in the early 2000. In particular, graph data analysis and graph data mining play an important role in biology interaction networks, where recent techniques of artificial intelligence, usually employed in other type of networks (e.g., social, citations, and trademark networks) aim to implement various data mining tasks including classification, clustering, recommendation, anomaly detection, and link prediction. The commitment and efforts of artificial intelligence research in network biology are motivated by the fact that machine learning techniques are often prohibitively computational demanding, low parallelizable, and ultimately inapplicable, since biological network of realistic size is a large system, which is characterised by a high density of interactions and often with a non-linear dynamics and a non-Euclidean latent geometry. Currently, graph embedding emerges as the new learning paradigm that shifts the tasks of building complex models for classification, clustering, and link prediction to learning an informative representation of the graph data in a vector space so that many graph mining and learning tasks can be more easily performed by employing efficient non-iterative traditional models (e.g., a linear support vector machine for the classification task). The great potential of graph embedding is the main reason of the flourishing of studies in this area and, in particular, the artificial intelligence learning techniques. In this mini review, we give a comprehensive summary of the main graph embedding algorithms in light of the recent burgeoning interest in geometric deep learning.

Studying stochastic systems biology of the cell with single-cell genomics data

Recent experimental developments in genome-wide RNA quantification hold considerable promise for systems biology. However, rigorously probing the biology of living cells requires a unified mathematical framework that accounts for single-molecule biological stochasticity in the context of technical variation associated with genomics assays. We review models for a variety of RNA transcription processes, as well as the encapsulation and library construction steps of microfluidics-based single-cell RNA sequencing, and present a framework to integrate these phenomena by the manipulation of generating functions. Finally, we use simulated scenarios and biological data to illustrate the implications and applications of the approach.

Uncovering Protein Networks in Cardiovascular Proteomics

Biological networks have been widely used in many different diseases to identify potential biomarkers and design drug targets. In the present review, we describe the main computational techniques for reconstructing and analyzing different types of protein networks and summarize the previous applications of such techniques in cardiovascular diseases. Existing tools are critically compared, discussing when each method is preferred such as the use of co-expression networks for functional annotation of protein clusters and the use of directed networks for inferring regulatory associations. Finally, we are presenting examples of reconstructing protein networks of different types (regulatory, co-expression, and protein-protein interaction networks). We demonstrate the necessity to reconstruct networks separately for each cardiovascular tissue type and disease entity and provide illustrative examples of the importance of taking into consideration relevant post-translational modifications. Finally, we demonstrate and discuss how the findings of protein networks could be interpreted using single-cell RNA-sequencing data.

From time-series transcriptomics to gene regulatory networks: A review on inference methods

Inference of gene regulatory networks has been an active area of research for around 20 years, leading to the development of sophisticated inference algorithms based on a variety of assumptions and approaches. With the ever increasing demand for more accurate and powerful models, the inference problem remains of broad scientific interest. The abstract representation of biological systems through gene regulatory networks represents a powerful method to study such systems, encoding different amounts and types of information. In this review, we summarize the different types of inference algorithms specifically based on time-series transcriptomics, giving an overview of the main applications of gene regulatory networks in computational biology. This review is intended to give an updated reference of regulatory networks inference tools to biologists and researchers new to the topic and guide them in selecting the appropriate inference method that best fits their questions, aims, and experimental data.

A general model-based causal inference method overcomes the curse of synchrony and indirect effect

To identify causation, model-free inference methods, such as Granger Causality, have been widely used due to their flexibility. However, they have difficulty distinguishing synchrony and indirect effects from direct causation, leading to false predictions. To overcome this, model-based inference methods that test the reproducibility of data with a specific mechanistic model to infer causality were developed. However, they can only be applied to systems described by a specific model, greatly limiting their applicability. Here, we address this limitation by deriving an easily testable condition for a general monotonic ODE model to reproduce time-series data. We built a user-friendly computational package, General ODE-Based Inference (GOBI), which is applicable to nearly any monotonic system with positive and negative regulations described by ODE. GOBI successfully inferred positive and negative regulations in various networks at both the molecular and population levels, unlike existing model-free methods. Thus, this accurate and broadly applicable inference method is a powerful tool for understanding complex dynamical systems.

Seiðr: Efficient calculation of robust ensemble gene networks

Gene regulatory and gene co-expression networks are powerful research tools for identifying biological signal within high-dimensional gene expression data. In recent years, research has focused on addressing shortcomings of these techniques with regard to the low signal-to-noise ratio, non-linear interactions and dataset dependent biases of published methods. Furthermore, it has been shown that aggregating networks from multiple methods provides improved results. Despite this, few useable and scalable software tools have been implemented to perform such best-practice analyses. Here, we present Seidr (stylized Seiðr), a software toolkit designed to assist scientists in gene regulatory and gene co-expression network inference. Seidr creates community networks to reduce algorithmic bias and utilizes noise corrected network backboning to prune noisy edges in the networks. Using benchmarks in real-world conditions across three eukaryotic model organisms, Saccharomyces cerevisiae, Drosophila melanogaster, and Arabidopsis thaliana, we show that individual algorithms are biased toward functional evidence for certain gene-gene interactions. We further demonstrate that the community network is less biased, providing robust performance across different standards and comparisons for the model organisms. Finally, we apply Seidr to a network of drought stress in Norway spruce (Picea abies (L.) H. Krast) as an example application in a non-model species. We demonstrate the use of a network inferred using Seidr for identifying key components, communities and suggesting gene function for non-annotated genes.

Studying stochastic systems biology of the cell with single-cell genomics data

Recent experimental developments in genome-wide RNA quantification hold considerable promise for systems biology. However, rigorously probing the biology of living cells requires a unified mathematical framework that accounts for single-molecule biological stochasticity in the context of technical variation associated with genomics assays. We review models for a variety of RNA transcription processes, as well as the encapsulation and library construction steps of microfluidics-based single-cell RNA sequencing, and present a framework to integrate these phenomena by the manipulation of generating functions. Finally, we use simulated scenarios and biological data to illustrate the implications and applications of the approach.

Benchmarking Gene Regulatory Network Inference Methods on Simulated and Experimental Data

Although the challenge of gene regulatory network inference has been studied for more than a decade, it is still unclear how well network inference methods work when applied to real data. Attempts to benchmark these methods on experimental data have yielded mixed results, in which sometimes even the best methods fail to outperform random guessing, and in other cases they perform reasonably well. So, one of the most valuable contributions one can currently make to the field of network inference is to benchmark methods on experimental data for which the true underlying network is already known, and report the results so that we can get a clearer picture of their efficacy. In this paper, we report results from the first, to our knowledge, benchmarking of network inference methods on single cell E. coli transcriptomic data. We report a moderate level of accuracy for the methods, better than random chance but still far from perfect. We also find that some methods that were quite strong and accurate on microarray and bulk RNA-seq data did not perform as well on the single cell data. Additionally, we benchmark a simple network inference method (Pearson correlation), on data generated through computer simulations in order to draw conclusions about general best practices in network inference studies. We predict that network inference would be more accurate using proteomic data rather than transcriptomic data, which could become relevant if high-throughput proteomic experimental methods are developed in the future. We also show through simulations that using a simplified model of gene expression that skips the mRNA step tends to substantially overestimate the accuracy of network inference methods, and advise against using this model for future in silico benchmarking studies.

EnsInfer: a simple ensemble approach to network inference outperforms any single method

This study evaluates both a variety of existing base causal inference methods and a variety of ensemble methods. We show that: (i) base network inference methods vary in their performance across different datasets, so a method that works poorly on one dataset may work well on another; (ii) a non-homogeneous ensemble method in the form of a Naive Bayes classifier leads overall to as good or better results than using the best single base method or any other ensemble method; (iii) for the best results, the ensemble method should integrate all methods that satisfy a statistical test of normality on training data. The resulting ensemble model EnsInfer easily integrates all kinds of RNA-seq data as well as new and existing inference methods. The paper categorizes and reviews state-of-the-art underlying methods, describes the EnsInfer ensemble approach in detail, and presents experimental results. The source code and data used will be made available to the community upon publication.

Inference of dynamic interaction networks: A comparison between Lotka-Volterra and multivariate autoregressive models

Networks are ubiquitous throughout biology, spanning the entire range from molecules to food webs and global environmental systems. Yet, despite substantial efforts by the scientific community, the inference of these networks from data still presents a problem that is unsolved in general. One frequent strategy of addressing the structure of networks is the assumption that the interactions among molecular or organismal populations are static and correlative. While often successful, these static methods are no panacea. They usually ignore the asymmetry of relationships between two species and inferences become more challenging if the network nodes represent dynamically changing quantities. Overcoming these challenges, two very different network inference approaches have been proposed in the literature: Lotka-Volterra (LV) models and Multivariate Autoregressive (MAR) models. These models are computational frameworks with different mathematical structures which, nevertheless, have both been proposed for the same purpose of inferring the interactions within coexisting population networks from observed time-series data. Here, we assess these dynamic network inference methods for the first time in a side-by-side comparison, using both synthetically generated and ecological datasets. Multivariate Autoregressive and Lotka-Volterra models are mathematically equivalent at the steady state, but the results of our comparison suggest that Lotka-Volterra models are generally superior in capturing the dynamics of networks with non-linear dynamics, whereas Multivariate Autoregressive models are better suited for analyses of networks of populations with process noise and close-to linear behavior. To the best of our knowledge, this is the first study comparing LV and MAR approaches. Both frameworks are valuable tools that address slightly different aspects of dynamic networks.

Network-based approaches for modeling disease regulation and progression

Molecular interaction networks lay the foundation for studying how biological functions are controlled by the complex interplay of genes and proteins. Investigating perturbed processes using biological networks has been instrumental in uncovering mechanisms that underlie complex disease phenotypes. Rapid advances in omics technologies have prompted the generation of high-throughput datasets, enabling large-scale, network-based analyses. Consequently, various modeling techniques, including network enrichment, differential network extraction, and network inference, have proven to be useful for gaining new mechanistic insights. We provide an overview of recent network-based methods and their core ideas to facilitate the discovery of disease modules or candidate mechanisms. Knowledge generated from these computational efforts will benefit biomedical research, especially drug development and precision medicine. We further discuss current challenges and provide perspectives in the field, highlighting the need for more integrative and dynamic network approaches to model disease development and progression.

DiffBrainNet: Differential analyses add new insights into the response to glucocorticoids at the level of genes, networks and brain regions

Genome-wide gene expression analyses are invaluable tools for studying biological and disease processes, allowing a hypothesis-free comparison of expression profiles. Traditionally, transcriptomic analysis has focused on gene-level effects found by differential expression. In recent years, network analysis has emerged as an important additional level of investigation, providing information on molecular connectivity, especially for diseases associated with a large number of linked effects of smaller magnitude, like neuropsychiatric disorders. Here, we describe how combined differential expression and prior-knowledge-based differential network analysis can be used to explore complex datasets. As an example, we analyze the transcriptional responses following administration of the glucocorticoid/stress receptor agonist dexamethasone in 8 mouse brain regions important for stress processing. By applying a combination of differential network- and expression-analyses, we find that these explain distinct but complementary biological mechanisms of the glucocorticoid responses. Additionally, network analysis identifies new differentially connected partners of risk genes and can be used to generate hypotheses on molecular pathways affected. With DiffBrainNet (http://diffbrainnet.psych.mpg.de), we provide an analysis framework and a publicly available resource for the study of the transcriptional landscape of the mouse brain which can identify molecular pathways important for basic functioning and response to glucocorticoids in a brain-region specific manner.

Gene regulatory network inference methodology for genomic and transcriptomic data acquired in genetically related heterozygote individuals

MOTIVATION

Inferring gene regulatory networks in non-independent genetically related panels is a methodological challenge. This hampers evolutionary and biological studies using heterozygote individuals such as in wild sunflower populations or cultivated hybrids.

RESULTS

First, we simulated 100 datasets of gene expressions and polymorphisms, displaying the same gene expression distributions, heterozygosities and heritabilities as in our dataset including 173 genes and 353 genotypes measured in sunflower hybrids. Secondly, we performed a meta-analysis based on six inference methods [least absolute shrinkage and selection operator (Lasso), Random Forests, Bayesian Networks, Markov Random Fields, Ordinary Least Square and fast inference of networks from directed regulation (Findr)] and selected the minimal density networks for better accuracy with 64 edges connecting 79 genes and 0.35 area under precision and recall (AUPR) score on average. We identified that triangles and mutual edges are prone to errors in the inferred networks. Applied on classical datasets without heterozygotes, our strategy produced a 0.65 AUPR score for one dataset of the DREAM5 Systems Genetics Challenge. Finally, we applied our method to an experimental dataset from sunflower hybrids. We successfully inferred a network composed of 105 genes connected by 106 putative regulations with a major connected component.

AVAILABILITY AND IMPLEMENTATION

Our inference methodology dedicated to genomic and transcriptomic data is available at https://forgemia.inra.fr/sunrise/inference_methods.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Review and assessment of Boolean approaches for inference of gene regulatory networks

Boolean descriptions of gene regulatory networks can provide an insight into interactions between genes. Boolean networks hold predictive power, are easy to understand, and can be used to simulate the observed networks in different scenarios. We review fundamental and state-of-the-art methods for inference of Boolean networks. We introduce a methodology for a straightforward evaluation of Boolean inference approaches based on the generation of evaluation datasets, application of selected inference methods, and evaluation of performance measures to guide the selection of the best method for a given inference problem. We demonstrate this procedure on inference methods REVEAL (REVerse Engineering ALgorithm), Best-Fit Extension, MIBNI (Mutual Information-based Boolean Network Inference), GABNI (Genetic Algorithm-based Boolean Network Inference) and ATEN (AND/OR Tree ENsemble algorithm), which infers Boolean descriptions of gene regulatory networks from discretised time series data. Boolean inference approaches tend to perform better in terms of dynamic accuracy, and slightly worse in terms of structural correctness. We believe that the proposed methodology and provided guidelines will help researchers to develop Boolean inference approaches with a good predictive capability while maintaining structural correctness and biological relevance.

Identification of T cell-related biomarkers for breast cancer based on weighted gene co-expression network analysis

Breast cancer is the most frequent malignancy worldwide, with immunotherapy and targeted therapy being key strategies to improving the prognosis. We downloaded mRNA expression dataset of breast cancer from The Cancer Genome Atlas (TCGA) database, and divided preprocessed genes into 12 modules based on gene expression profile by weighted gene co-expression network analysis (WGCNA). The StromalScore, ImmuneScore and ESTIMATEScore of samples were assessed. The Kaplan-Meier curve showed that ImmuneScore was notably correlated with breast cancer patient's prognosis. By analyzing the connectivity between module eigengenes and clinical traits, the gene module closely related to ImmuneScore was obtained. Further, through intramodular gene connectivity and protein-protein interaction network topology analysis of module genes, hub genes (HLA-E, HLA-DPB1 and HLA-DRB1) in immune-related module were screened out. Finally, bioinformatics analysis displayed that HLA-DPB1 and HLA-DRB1 were notably overexpressed and HLA-E was underexpressed in breast cancer tissues. TIMER database analysis showed that three hub gene levels were significantly correlated with infiltration levels of CD8+ T cells and CD4+ T cells. Meanwhile, Pearson correlation analysis revealed positive correlation between three hub genes and those of immune checkpoint genes (LAG3, PD-1, PD-L1). Additionally, prognosis could be effectively evaluated by HLA-DPB1 and HLA-DRB1 levels, and differentially activated signalling pathways between high- and low-expression groups of HLA-E and HLA-DPB1 were obtained by gene set enrichment analysis. To conclude, this study identified three T cell-related biomarkers for breast cancer based on TCGA-BRCA dataset, and the screened genes could provide references for breast cancer immunotherapy.

Data-driven dynamical modelling of a pathogen-infected plant gene regulatory network: A comparative analysis

Recent advances in synthetic biology have enabled the design of genetic feedback control circuits that could be implemented to build resilient plants against pathogen attacks. To facilitate the proper design of these genetic feedback control circuits, an accurate model that is able to capture the vital dynamical behaviour of the pathogen-infected plant is required. In this study, using a data-driven modelling approach, we develop and compare four dynamical models (i.e. linear, Michaelis-Menten with Hill coefficient (Hill Function), standard S-System and extended S-System) of a pathogen-infected plant gene regulatory network (GRN). These models are then assessed across several criteria, i.e. ease of identifying the type of gene regulation, the predictive capability, Akaike Information Criterion (AIC) and the robustness to parameter uncertainty to determine its viability of balancing between biological complexity and accuracy when modelling the pathogen-infected plant GRN. Using our defined ranking score, we obtain the following insights to the modelling of GRN. Our analyses show that despite commonly used and provide biological relevance, the Hill Function model ranks the lowest while the extended S-System model ranks highest in the overall comparison. Interestingly, the performance of the linear model is more consistent throughout the comparison, making it the preferred model for this pathogen-infected plant GRN when considering data-driven modelling approach.

Regulatory Modules of Metabolites and Protein Phosphorylation in Arabidopsis Genotypes With Altered Sucrose Allocation

Multi-omics data sets are increasingly being used for the interpretation of cellular processes in response to environmental cues. Especially, the posttranslational modification of proteins by phosphorylation is an important regulatory process affecting protein activity and/or localization, which, in turn, can have effects on metabolic processes and metabolite levels. Despite this importance, relationships between protein phosphorylation status and metabolite abundance remain largely underexplored. Here, we used a phosphoproteomics-metabolomics data set collected at the end of day and night in shoots and roots of Arabidopsis to propose regulatory relationships between protein phosphorylation and accumulation or allocation of metabolites. For this purpose, we introduced a novel, robust co-expression measure suited to the structure of our data sets, and we used this measure to construct metabolite-phosphopeptide networks. These networks were compared between wild type and plants with perturbations in key processes of sugar metabolism, namely, sucrose export (sweet11/12 mutant) and starch synthesis (pgm mutant). The phosphopeptide-metabolite network turned out to be highly sensitive to perturbations in sugar metabolism. Specifically, KING1, the regulatory subunit of SnRK1, was identified as a primary candidate connecting protein phosphorylation status with metabolism. We additionally identified strong changes in the fatty acid network of the sweet11/12 mutant, potentially resulting from a combination of fatty acid signaling and metabolic overflow reactions in response to high internal sucrose concentrations. Our results further suggest novel protein-metabolite relationships as candidates for future targeted research.

Inference of Molecular Regulatory Systems Using Statistical Path-Consistency Algorithm

One of the key challenges in systems biology and molecular sciences is how to infer regulatory relationships between genes and proteins using high-throughout omics datasets. Although a wide range of methods have been designed to reverse engineer the regulatory networks, recent studies show that the inferred network may depend on the variable order in the dataset. In this work, we develop a new algorithm, called the statistical path-consistency algorithm (SPCA), to solve the problem of the dependence of variable order. This method generates a number of different variable orders using random samples, and then infers a network by using the path-consistent algorithm based on each variable order. We propose measures to determine the edge weights using the corresponding edge weights in the inferred networks, and choose the edges with the largest weights as the putative regulations between genes or proteins. The developed method is rigorously assessed by the six benchmark networks in DREAM challenges, the mitogen-activated protein (MAP) kinase pathway, and a cancer-specific gene regulatory network. The inferred networks are compared with those obtained by using two up-to-date inference methods. The accuracy of the inferred networks shows that the developed method is effective for discovering molecular regulatory systems.

Emerging landscape of molecular interaction networks: Opportunities, challenges and prospects

Network biology finds application in interpreting molecular interaction networks and providing insightful inferences using graph theoretical analysis of biological systems. The integration of computational bio-modelling approaches with different hybrid network-based techniques provides additional information about the behaviour of complex systems. With increasing advances in high-throughput technologies in biological research, attempts have been made to incorporate this information into network structures, which has led to a continuous update of network biology approaches over time. The newly minted centrality measures accommodate the details of omics data and regulatory network structure information. The unification of graph network properties with classical mathematical and computational modelling approaches and technologically advanced approaches like machine-learning- and artificial intelligence-based algorithms leverages the potential application of these techniques. These computational advances prove beneficial and serve various applications such as essential gene prediction, identification of drug–disease interaction and gene prioritization. Hence, in this review, we have provided a comprehensive overview of the emerging landscape of molecular interaction networks using graph theoretical approaches. With the aim to provide information on the wide range of applications of network biology approaches in understanding the interaction and regulation of genes, proteins, enzymes and metabolites at different molecular levels, we have reviewed the methods that utilize network topological properties, emerging hybrid network-based approaches and applications that integrate machine learning techniques to analyse molecular interaction networks. Further, we have discussed the applications of these approaches in biomedical research with a note on future prospects.

Deep neural network prediction of genome-wide transcriptome signatures - beyond the Black-box

Prediction algorithms for protein or gene structures, including transcription factor binding from sequence information, have been transformative in understanding gene regulation. Here we ask whether human transcriptomic profiles can be predicted solely from the expression of transcription factors (TFs). We find that the expression of 1600 TFs can explain >95% of the variance in 25,000 genes. Using the light-up technique to inspect the trained NN, we find an over-representation of known TF-gene regulations. Furthermore, the learned prediction network has a hierarchical organization. A smaller set of around 125 core TFs could explain close to 80% of the variance. Interestingly, reducing the number of TFs below 500 induces a rapid decline in prediction performance. Next, we evaluated the prediction model using transcriptional data from 22 human diseases. The TFs were sufficient to predict the dysregulation of the target genes (rho = 0.61, P < 10⁻²¹⁶). By inspecting the model, key causative TFs could be extracted for subsequent validation using disease-associated genetic variants. We demonstrate a methodology for constructing an interpretable neural network predictor, where analyses of the predictors identified key TFs that were inducing transcriptional changes during disease.

Evaluating the cytotoxicity and pathogenicity of multi-walled carbon nanotube through weighted gene co-expression network analysis: a nanotoxicogenomics study

Background

Multi-walled carbon nanotube (MWCNT) is one of the most momentous carbonaceous nanoparticles which is widely used for various applications such as electronics, vehicles, and therapeutics. However, their possible toxicity and adverse effects convert them into a major health threat for humans and animals.

Results

In this study, we employed weighted gene co-expression network analysis (WGCNA) to identify the co-expressed gene groups and dysregulated pathways due to the MWCNT exposure. For this purpose, three weighted gene co-expression networks for the microarray gene expression profiles of the mouse after 1, 6, and 12-month post-exposure to MWCNT were constructed. The module-trait analysis specified the significant modules related to different doses (1, 10, 40, and 80 µg) of MWCNT. Afterward, common genes between co-regulated and differentially expressed genes were determined. The further pathway analysis highlighted the enrichment of genes including Actb, Ube2b, Psme3, Ezh2, Alas2, S100a10, Ypel5, Rhoa, Rac1, Ube2l6, Prdx2, Ctsb, Bnip3l, Gp6, Myh9, Ube2k, Mbnl1, Kbtbd8, Riok3, Itgb1, Rap1a, and Atp5h in immune-, inflammation-, and protein metabolism-related pathways.

Conclusions

This study discloses the genotoxicity and cytotoxicity effects of various doses of MWCNT which also affect the metabolism system. The identified genes can serve as potential biomarkers and therapeutic candidates. However, further studies should be performed to validate them in human cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01031-3.

Metabolic and Transcriptional Changes across Osteogenic Differentiation of Mesenchymal Stromal Cells

Mesenchymal stromal cells (MSCs) are multipotent post-natal stem cells with applications in tissue engineering and regenerative medicine. MSCs can differentiate into osteoblasts, chondrocytes, or adipocytes, with functional differences in cells during osteogenesis accompanied by metabolic changes. The temporal dynamics of these metabolic shifts have not yet been fully characterized and are suspected to be important for therapeutic applications such as osteogenesis optimization. Here, our goal was to characterize the metabolic shifts that occur during osteogenesis. We profiled five key extracellular metabolites longitudinally (glucose, lactate, glutamine, glutamate, and ammonia) from MSCs from four donors to classify osteogenic differentiation into three metabolic stages, defined by changes in the uptake and secretion rates of the metabolites in cell culture media. We used a combination of untargeted metabolomic analysis, targeted analysis of ¹³C-glucose labelled intracellular data, and RNA-sequencing data to reconstruct a gene regulatory network and further characterize cellular metabolism. The metabolic stages identified in this proof-of-concept study provide a framework for more detailed investigations aimed at identifying biomarkers of osteogenic differentiation and small molecule interventions to optimize MSC differentiation for clinical applications.

Rescaling Biology: Increasing Integration Across Biological Scales and Subdisciplines to Enhance Understanding and Prediction

The organization of the living world covers a vast range of spatiotemporal scales, from molecules to the biosphere, seconds to centuries. Biologists working within specialized subdisciplines tend to focus on different ranges of scales. Therefore, developing frameworks that enable testing questions and predictions of scaling require sufficient understanding of complex processes across biological subdisciplines and spatiotemporal scales. Frameworks that enable scaling across subdisciplines would ideally allow us to test hypotheses about the degree to which explicit integration across spatiotemporal scales is needed for predicting the outcome of biological processes. For instance, how does genomic variation within populations allow us to explain community structure? How do the dynamics of cellular metabolism translate to our understanding of whole-ecosystem metabolism? Do patterns and processes operate seamlessly across biological scales, or are there fundamental laws of biological scaling that limit our ability to make predictions from one scale to another? Similarly, can sub-organismal structures and processes be sufficiently understood in isolation of potential feedbacks from the population, community, or ecosystem levels? And can we infer the sub-organismal processes from data on the population, community, or ecosystem scale? Concerted efforts to develop more cross-disciplinary frameworks will open doors to a more fully integrated field of biology. In this paper we discuss how we might integrate across scales, specifically by 1. Identifying scales and boundaries, 2. Determining analogous units and processes across scales, 3. Developing frameworks to unite multiple scales, and 4. Extending frameworks to new empirical systems.

Inferring causality in biological oscillators

MOTIVATION

Fundamental to biological study is identifying regulatory interactions. The recent surge in time-series data collection in biology provides a unique opportunity to infer regulations computationally. However, when components oscillate, model-free inference methods, while easily implemented, struggle to distinguish periodic synchrony and causality. Alternatively, model-based methods test the reproducibility of time series given a specific model but require inefficient simulations and have limited applicability.

RESULTS

We develop an inference method based on a general model of molecular, neuronal and ecological oscillatory systems that merges the advantages of both model-based and model-free methods, namely accuracy, broad applicability and usability. Our method successfully infers the positive and negative regulations within various oscillatory networks, e.g. the repressilator and a network of cofactors at the pS2 promoter, outperforming popular inference methods.

AVAILABILITY AND IMPLEMENTATION

We provide a computational package, ION (Inferring Oscillatory Networks), that users can easily apply to noisy, oscillatory time series to uncover the mechanisms by which diverse systems generate oscillations. Accompanying MATLAB code under a BSD-style license and examples are available at https://github.com/Mathbiomed/ION. Additionally, the code is available under a CC-BY 4.0 License at https://doi.org/10.6084/m9.figshare.16431408.v1.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Systems biology for crop improvement

In recent years, generation of large-scale data from genome, transcriptome, proteome, metabolome, epigenome, and others, has become routine in several plant species. Most of these datasets in different crop species, however, were studied independently and as a result, full insight could not be gained on the molecular basis of complex traits and biological networks. A systems biology approach involving integration of multiple omics data, modeling, and prediction of the cellular functions is required to understand the flow of biological information that underlies complex traits. In this context, systems biology with multiomics data integration is crucial and allows a holistic understanding of the dynamic system with the different levels of biological organization interacting with external environment for a phenotypic expression. Here, we present recent progress made in the area of various omics studies-integrative and systems biology approaches with a special focus on application to crop improvement. We have also discussed the challenges and opportunities in multiomics data integration, modeling, and understanding of the biology of complex traits underpinning yield and stress tolerance in major cereals and legumes.

Integrated Inference of Asymmetric Protein Interaction Networks Using Dynamic Model and Individual Patient Proteomics Data

Recent advances in experimental biology studies have produced large amount of molecular activity data. In particular, individual patient data provide non-time series information for the molecular activities in disease conditions. The challenge is how to design effective algorithms to infer regulatory networks using the individual patient datasets and consequently address the issue of network symmetry. This work is aimed at developing an efficient pipeline to reverse-engineer regulatory networks based on the individual patient proteomic data. The first step uses the SCOUT algorithm to infer the pseudo-time trajectory of individual patients. Then the path-consistent method with part mutual information is used to construct a static network that contains the potential protein interactions. To address the issue of network symmetry in terms of undirected symmetric network, a dynamic model of ordinary differential equations is used to further remove false interactions to derive asymmetric networks. In this work a dataset from triple-negative breast cancer patients is used to develop a protein-protein interaction network with 15 proteins.

Computational analysis of fused co-expression networks for the identification of candidate cancer gene biomarkers

The complexity of cancer has always been a huge issue in understanding the source of this disease. However, by appreciating its complexity, we can shed some light on crucial gene associations across and in specific cancer types. In this study, we develop a general framework to infer relevant gene biomarkers and their gene-to-gene associations using multiple gene co-expression networks for each cancer type. Specifically, we infer computationally and biologically interesting communities of genes from kidney renal clear cell carcinoma, liver hepatocellular carcinoma, and prostate adenocarcinoma data sets of The Cancer Genome Atlas (TCGA) database. The gene communities are extracted through a data-driven pipeline and then evaluated through both functional analyses and literature findings. Furthermore, we provide a computational validation of their relevance for each cancer type by comparing the performance of normal/cancer classification for our identified gene sets and other gene signatures, including the typically-used differentially expressed genes. The hallmark of this study is its approach based on gene co-expression networks from different similarity measures: using a combination of multiple gene networks and then fusing normal and cancer networks for each cancer type, we can have better insights on the overall structure of the cancer-type-specific network.

Measuring topological descriptors of complex networks under uncertainty

Revealing the structural features of a complex system from the observed collective dynamics is a fundamental problem in network science. To compute the various topological descriptors commonly used to characterize the structure of a complex system (e.g., the degree, the clustering coefficient, etc.), it is usually necessary to completely reconstruct the network of relations between the subsystems. Several methods are available to detect the existence of interactions between the nodes of a network. By observing some physical quantities through time, the structural relationships are inferred using various discriminating statistics (e.g., correlations, mutual information, etc.). In this setting, the uncertainty about the existence of the edges is reflected in the uncertainty about the topological descriptors. In this study, we propose a methodological framework to evaluate this uncertainty, replacing the topological descriptors, even at the level of a single node, with appropriate probability distributions, eluding the reconstruction phase. Our theoretical framework agrees with the numerical experiments performed on a large set of synthetic and real-world networks. Our results provide a grounded framework for the analysis and the interpretation of widely used topological descriptors, such as degree centrality, clustering, and clusters, in scenarios in which the existence of network connectivity is statistically inferred or when the probabilities of existence π_{ij} of the edges are known. To this purpose, we also provide a simple and mathematically grounded process to transform the discriminating statistics into the probabilities π_{ij}.

A comprehensive overview and critical evaluation of gene regulatory network inference technologies

Gene regulatory network (GRN) is the important mechanism of maintaining life process, controlling biochemical reaction and regulating compound level, which plays an important role in various organisms and systems. Reconstructing GRN can help us to understand the molecular mechanism of organisms and to reveal the essential rules of a large number of biological processes and reactions in organisms. Various outstanding network reconstruction algorithms use specific assumptions that affect prediction accuracy, in order to deal with the uncertainty of processing. In order to study why a certain method is more suitable for specific research problem or experimental data, we conduct research from model-based, information-based and machine learning-based method classifications. There are obviously different types of computational tools that can be generated to distinguish GRNs. Furthermore, we discuss several classical, representative and latest methods in each category to analyze core ideas, general steps, characteristics, etc. We compare the performance of state-of-the-art GRN reconstruction technologies on simulated networks and real networks under different scaling conditions. Through standardized performance metrics and common benchmarks, we quantitatively evaluate the stability of various methods and the sensitivity of the same algorithm applying to different scaling networks. The aim of this study is to explore the most appropriate method for a specific GRN, which helps biologists and medical scientists in discovering potential drug targets and identifying cancer biomarkers.

Systems Biology Approaches to Understanding the Human Immune System

Systems biology is an approach to interrogate complex biological systems through large-scale quantification of numerous biomolecules. The immune system involves >1,500 genes/proteins in many interconnected pathways and processes, and a systems-level approach is critical in broadening our understanding of the immune response to vaccination. Changes in molecular pathways can be detected using high-throughput omics datasets (e.g., transcriptomics, proteomics, and metabolomics) by using methods such as pathway enrichment, network analysis, machine learning, etc. Importantly, integration of multiple omic datasets is becoming key to revealing novel biological insights. In this perspective article, we highlight the use of protein-protein interaction (PPI) networks as a multi-omics integration approach to unravel information flow and mechanisms during complex biological events, with a focus on the immune system. This involves a combination of tools, including: InnateDB, a database of curated interactions between genes and protein products involved in the innate immunity; NetworkAnalyst, a visualization and analysis platform for InnateDB interactions; and MetaBridge, a tool to integrate metabolite data into PPI networks. The application of these systems techniques is demonstrated for a variety of biological questions, including: the developmental trajectory of neonates during the first week of life, mechanisms in host-pathogen interaction, disease prognosis, biomarker discovery, and drug discovery and repurposing. Overall, systems biology analyses of omics data have been applied to a variety of immunology-related questions, and here we demonstrate the numerous ways in which PPI network analysis can be a powerful tool in contributing to our understanding of the immune system and the study of vaccines.