Predicting perturbation patterns from the topology of biological networks

Topological signatures in regulatory network enable phenotypic heterogeneity in small cell lung cancer

Phenotypic (non-genetic) heterogeneity has significant implications for the development and evolution of organs, organisms, and populations. Recent observations in multiple cancers have unraveled the role of phenotypic heterogeneity in driving metastasis and therapy recalcitrance. However, the origins of such phenotypic heterogeneity are poorly understood in most cancers. Here, we investigate a regulatory network underlying phenotypic heterogeneity in small cell lung cancer, a devastating disease with no molecular targeted therapy. Discrete and continuous dynamical simulations of this network reveal its multistable behavior that can explain co-existence of four experimentally observed phenotypes. Analysis of the network topology uncovers that multistability emerges from two teams of players that mutually inhibit each other, but members of a team activate one another, forming a 'toggle switch' between the two teams. Deciphering these topological signatures in cancer-related regulatory networks can unravel their 'latent' design principles and offer a rational approach to characterize phenotypic heterogeneity in a tumor.

Percolation of heterogeneous flows uncovers the bottlenecks of infrastructure networks

Whether it be the passengers’ mobility demand in transportation systems, or the consumers’ energy demand in power grids, the primary purpose of many infrastructure networks is to best serve this flow demand. In reality, the volume of flow demand fluctuates unevenly across complex networks while simultaneously being hindered by some form of congestion or overload. Nevertheless, there is little known about how the heterogeneity of flow demand influences the network flow dynamics under congestion. To explore this, we introduce a percolation-based network analysis framework underpinned by flow heterogeneity. Thereby, we theoretically identify bottleneck links with guaranteed decisive impact on how flows are passed through the network. The effectiveness of the framework is demonstrated on large-scale real transportation networks, where mitigating the congestion on a small fraction of the links identified as bottlenecks results in a significant network improvement.

Infrastructure networks are characterized by fluctuations of flow demand between different points and temporal congestion or overload on flow pathways. Hamedmoghadam et al. identify congestion bottlenecks in networks relevant to communication, transportation, water supply, and power distribution.

ComHub: Community predictions of hubs in gene regulatory networks

BACKGROUND

Hub transcription factors, regulating many target genes in gene regulatory networks (GRNs), play important roles as disease regulators and potential drug targets. However, while numerous methods have been developed to predict individual regulator-gene interactions from gene expression data, few methods focus on inferring these hubs.

RESULTS

We have developed ComHub, a tool to predict hubs in GRNs. ComHub makes a community prediction of hubs by averaging over predictions by a compendium of network inference methods. Benchmarking ComHub against the DREAM5 challenge data and two independent gene expression datasets showed a robust performance of ComHub over all datasets.

CONCLUSIONS

In contrast to other evaluated methods, ComHub consistently scored among the top performing methods on data from different sources. Lastly, we implemented ComHub to work with both predefined networks and to perform stand-alone network inference, which will make the method generally applicable.

Generating Ensembles of Gene Regulatory Networks to Assess Robustness of Disease Modules

The use of biological networks such as protein-protein interaction and transcriptional regulatory networks is becoming an integral part of genomics research. However, these networks are not static, and during phenotypic transitions like disease onset, they can acquire new "communities" (or highly interacting groups) of genes that carry out cellular processes. Disease communities can be detected by maximizing a modularity-based score, but since biological systems and network inference algorithms are inherently noisy, it remains a challenge to determine whether these changes represent real cellular responses or whether they appeared by random chance. Here, we introduce Constrained Random Alteration of Network Edges (CRANE), a method for randomizing networks with fixed node strengths. CRANE can be used to generate a null distribution of gene regulatory networks that can in turn be used to rank the most significant changes in candidate disease communities. Compared to other approaches, such as consensus clustering or commonly used generative models, CRANE emulates biologically realistic networks and recovers simulated disease modules with higher accuracy. When applied to breast and ovarian cancer networks, CRANE improves the identification of cancer-relevant GO terms while reducing the signal from non-specific housekeeping processes.

Exploring Biocontrol Agents From Microbial Keystone Taxa Associated to Suppressive Soil: A New Attempt for a Biocontrol Strategy

Recent studies have observed differing microbiomes between disease-suppressive and disease-conducive soils. However, it remains unclear whether the microbial keystone taxa in suppressive soil are critical for the suppression of diseases. Bacterial wilt is a common soil-borne disease caused by Ralstonia solanacearum that affects tobacco plants. In this study, two contrasting tobacco fields with bacterial wilt disease incidences of 0% (disease suppressive) and 100% (disease conducive) were observed. Through amplicon sequencing, as expected, a high abundance of Ralstonia was found in the disease-conducive soil, while large amounts of potential beneficial bacteria were found in the disease-suppressive soil. In the fungal community, an abundance of the Fusarium genus, which contains species that cause Fusarium wilt, showed a positive correlation (p < 0.001) with the abundance of Ralstonia. Network analysis revealed that the healthy plants had more complex bacterial networks than the diseased plants. A total of 9 and 13 bacterial keystone taxa were identified from the disease-suppressive soil and healthy root, respectively. Accumulated abundance of these bacterial keystones showed a negative correlation (p < 0.001) with the abundance of Ralstonia. To complement network analysis, culturable strains were isolated, and three species belonging to Pseudomonas showed high 16S rRNA gene similarity (98.4-100%) with keystone taxa. These strains displayed strong inhibition on pathogens and reduced the incidence of bacterial wilt disease in greenhouse condition. This study highlighted the importance of keystone species in the protection of crops against pathogen infection and proposed an approach to obtain beneficial bacteria through identifying keystone species, avoiding large-scale bacterial isolation and cultivation.

DReSS: a method to quantitatively describe the influence of structural perturbations on state spaces of genetic regulatory networks

Structures of genetic regulatory networks are not fixed. These structural perturbations can cause changes to the reachability of systems' state spaces. As system structures are related to genotypes and state spaces are related to phenotypes, it is important to study the relationship between structures and state spaces. However, there is still no method can quantitively describe the reachability differences of two state spaces caused by structural perturbations. Therefore, Difference in Reachability between State Spaces (DReSS) is proposed. DReSS index family can quantitively describe differences of reachability, attractor sets between two state spaces and can help find the key structure in a system, which may influence system's state space significantly. First, basic properties of DReSS including non-negativity, symmetry and subadditivity are proved. Then, typical examples are shown to explain the meaning of DReSS and the differences between DReSS and traditional graph distance. Finally, differences of DReSS distribution between real biological regulatory networks and random networks are compared. Results show most structural perturbations in biological networks tend to affect reachability inside and between attractor basins rather than to affect attractor set itself when compared with random networks, which illustrates that most genotype differences tend to influence the proportion of different phenotypes and only a few ones can create new phenotypes. DReSS can provide researchers with a new insight to study the relation between genotypes and phenotypes.

Cellular reprogramming: Mathematics meets medicine

Generating needed cell types using cellular reprogramming is a promising strategy for restoring tissue function in injury or disease. A common method for reprogramming is addition of one or more transcription factors that confer a new function or identity. Advancements in transcription factor selection and delivery have culminated in successful grafting of autologous reprogrammed cells, an early demonstration of their clinical utility. Though cellular reprogramming has been successful in a number of settings, identification of appropriate transcription factors for a particular transformation has been challenging. Computational methods enable more sophisticated prediction of relevant transcription factors for reprogramming by leveraging gene expression data of initial and target cell types, and are built on mathematical frameworks ranging from information theory to control theory. This review highlights the utility and impact of these mathematical frameworks in the field of cellular reprogramming. This article is categorized under: Reproductive System Diseases > Reproductive System Diseases>Genetics/Genomics/Epigenetics Reproductive System Diseases > Reproductive System Diseases>Stem Cells and Development Reproductive System Diseases > Reproductive System Diseases>Computational Models.

Drug2ways: Reasoning over causal paths in biological networks for drug discovery

Elucidating the causal mechanisms responsible for disease can reveal potential therapeutic targets for pharmacological intervention and, accordingly, guide drug repositioning and discovery. In essence, the topology of a network can reveal the impact a drug candidate may have on a given biological state, leading the way for enhanced disease characterization and the design of advanced therapies. Network-based approaches, in particular, are highly suited for these purposes as they hold the capacity to identify the molecular mechanisms underlying disease. Here, we present drug2ways, a novel methodology that leverages multimodal causal networks for predicting drug candidates. Drug2ways implements an efficient algorithm which reasons over causal paths in large-scale biological networks to propose drug candidates for a given disease. We validate our approach using clinical trial information and demonstrate how drug2ways can be used for multiple applications to identify: i) single-target drug candidates, ii) candidates with polypharmacological properties that can optimize multiple targets, and iii) candidates for combination therapy. Finally, we make drug2ways available to the scientific community as a Python package that enables conducting these applications on multiple standard network formats.

Topological Determinants of Perturbation Spreading in Networks

Spreading phenomena essentially underlie the dynamics of various natural and technological networked systems, yet how spatiotemporal propagation patterns emerge from such networks remains largely unknown. Here we propose a novel approach that reveals universal features determining the spreading dynamics in diffusively coupled networks and disentangles them from factors that are system specific. In particular, we first analytically identify a purely topological factor encoding the interaction structure and strength, and second, numerically estimate a master function characterizing the universal scaling of the perturbation arrival times across topologically different networks. The proposed approach thereby provides intuitive insights into complex propagation patterns as well as accurate predictions for the perturbation arrival times. The approach readily generalizes to a wide range of networked systems with diffusive couplings and may contribute to assess the risks of transient influences of ubiquitous perturbations in real-world systems.

The Rhodamine Isothiocyanate Analogue as a Quorum Sensing Inhibitor Has the Potential to Control Microbially-Induced Biofouling

Quorum sensing inhibitors (QSIs) have been proven to be an innovative approach to interfering with biofilm formation, since this process is regulated by QS signals. However, most studies have focused on single-species biofilm formation, whereas studies of the effects of signal interference on the development of multispecies biofilm, especially in the natural environment, are still lacking. Here we develop and evaluate the anti-biofilm capability of a new QSI (rhodamine isothiocyanate analogue, RIA) in natural seawater. During the experiment, biofilm characteristics, microbial communities/functions and network interactions were monitored at 36, 80, and 180 h, respectively. The results showed that the biomass and 3D structure of the biofilm were significantly different in the presence of the QSI. The expression of genes involved in extracellular polysaccharide synthesis was also downregulated in the QSI-treated group. Dramatic differences in microbial composition, β-diversity and functions between the RIA-treated group and the control group were also observed, especially in the early stage of biofilm development. Furthermore, co-occurrence model analysis showed that RIA reduced the complexity of the microbial network. This study demonstrates that rhodamine isothiocyanate analogue is an efficient QS inhibitor and has potential applications in controlling biofouling caused by multispecies biofilm, especially in the early stage of biofouling formation.

Fractional diffusion on the human proteome as an alternative to the multi-organ damage of SARS-CoV-2

The coronavirus 2019 (COVID-19) respiratory disease is caused by the novel coronavirus SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), which uses the enzyme ACE2 to enter human cells. This disease is characterized by important damage at a multi-organ level, partially due to the abundant expression of ACE2 in practically all human tissues. However, not every organ in which ACE2 is abundant is affected by SARS-CoV-2, which suggests the existence of other multi-organ routes for transmitting the perturbations produced by the virus. We consider here diffusive processes through the protein-protein interaction (PPI) network of proteins targeted by SARS-CoV-2 as an alternative route. We found a subdiffusive regime that allows the propagation of virus perturbations through the PPI network at a significant rate. By following the main subdiffusive routes across the PPI network, we identify proteins mainly expressed in the heart, cerebral cortex, thymus, testis, lymph node, kidney, among others of the organs reported to be affected by COVID-19.

Sensitivity of collective outcomes identifies pivotal components

A social system is susceptible to perturbation when its collective properties depend sensitively on a few pivotal components. Using the information geometry of minimal models from statistical physics, we develop an approach to identify pivotal components to which coarse-grained, or aggregate, properties are sensitive. As an example, we introduce our approach on a reduced toy model with a median voter who always votes in the majority. The sensitivity of majority-minority divisions to changing voter behaviour pinpoints the unique role of the median. More generally, the sensitivity identifies pivotal components that precisely determine collective outcomes generated by a complex network of interactions. Using perturbations to target pivotal components in the models, we analyse datasets from political voting, finance and Twitter. Across these systems, we find remarkable variety, from systems dominated by a median-like component to those whose components behave more equally. In the context of political institutions such as courts or legislatures, our methodology can help describe how changes in voters map to new collective voting outcomes. For economic indices, differing system response reflects varying fiscal conditions across time. Thus, our information-geometric approach provides a principled, quantitative framework that may help assess the robustness of collective outcomes to targeted perturbation and compare social institutions, or even biological networks, with one another and across time.

Predicting miRNA-based disease-disease relationships through network diffusion on multi-omics biological data

With critical roles in regulating gene expression, miRNAs are strongly implicated in the pathophysiology of many complex diseases. Experimental methods to determine disease related miRNAs are time consuming and costly. Computationally predicting miRNA-disease associations has potential applications in finding miRNA therapeutic pathways and in understanding the role of miRNAs in disease-disease relationships. In this study, we propose the MiRNA-disease Association Prediction (MAP) method, an in-silico method to predict and prioritize miRNA-disease associations. The MAP method applies a network diffusion approach, starting from the known disease genes in a heterogenous network constructed from miRNA-gene associations, protein-protein interactions, and gene-disease associations. Validation using experimental data on miRNA-disease associations demonstrated superior performance to two current state-of-the-art methods, with areas under the ROC curve all over 0.8 for four types of cancer. MAP is successfully applied to predict differential miRNA expression in four cancer types. Most strikingly, disease-disease relationships in terms of shared miRNAs revealed hidden disease subtyping comparable to that of previous work on shared genes between diseases, with applications for multi-omics characterization of disease relationships.

The exposome and health: Where chemistry meets biology

Despite extensive evidence showing that exposure to specific chemicals can lead to disease, current research approaches and regulatory policies fail to address the chemical complexity of our world. To safeguard current and future generations from the increasing number of chemicals polluting our environment, a systematic and agnostic approach is needed. The "exposome" concept strives to capture the diversity and range of exposures to synthetic chemicals, dietary constituents, psychosocial stressors, and physical factors, as well as their corresponding biological responses. Technological advances such as high-resolution mass spectrometry and network science have allowed us to take the first steps toward a comprehensive assessment of the exposome. Given the increased recognition of the dominant role that nongenetic factors play in disease, an effort to characterize the exposome at a scale comparable to that of the human genome is warranted.

Opportunistic bacteria use quorum sensing to disturb coral symbiotic communities and mediate the occurrence of coral bleaching

Coral associated microorganisms, especially some opportunistic pathogens can utilize quorum-sensing (QS) signals to affect population structure and host health. However, direct evidence about the link between coral bleaching and dysbiotic microbiomes under QS regulation was lacking. Here, using 11 opportunistic bacteria and their QS products (AHLs, acyl-homoserine-lactones), we exposed Pocillopora damicornis to three different treatments: test groups (A and B: mixture of AHLs-producing bacteria and cocktail of AHLs signals respectively); control groups (C and D: group A and B with furanone added respectively); and a blank control (group E: only seawater) for 21 days. The results showed that remarkable bleaching phenomenon was observed in groups A and B. The operational taxonomic units-sequencing analysis shown that the bacterial network interactions and communities composition were significantly changed, becoming especially enhanced in the relative abundances of Vibrio, Edwardsiella, Enterobacter, Pseudomonas, and Aeromonas. Interestingly, the control groups (C and D) were found to have a limited influence upon host microbial composition and reduced bleaching susceptibility of P. damicornis. These results indicate bleaching's initiation and progression may be caused by opportunistic bacteria of resident microbes in a process under regulation by AHLs. These findings add a new dimension to our understanding of the complexity of bleaching mechanisms from a chemoecological perspective.

LiPLike: towards gene regulatory network predictions of high certainty

MOTIVATION

High correlation in expression between regulatory elements is a persistent obstacle for the reverse-engineering of gene regulatory networks. If two potential regulators have matching expression patterns, it becomes challenging to differentiate between them, thus increasing the risk of false positive identifications.

RESULTS

To allow for gene regulation predictions of high confidence, we propose a novel method, the Linear Profile Likelihood (LiPLike), that assumes a regression model and iteratively searches for interactions that cannot be replaced by a linear combination of other predictors. To compare the performance of LiPLike with other available inference methods, we benchmarked LiPLike using three independent datasets from the Dialogue on Reverse Engineering Assessment and Methods 5 (DREAM5) network inference challenge. We found that LiPLike could be used to stratify predictions of other inference tools, and when applied to the predictions of DREAM5 participants, we observed an average improvement in accuracy of >140% compared to individual methods. Furthermore, LiPLike was able to independently predict networks better than all DREAM5 participants when applied to biological data. When predicting the Escherichia coli network, LiPLike had an accuracy of 0.38 for the top-ranked 100 interactions, whereas the corresponding DREAM5 consensus model yielded an accuracy of 0.11.

AVAILABILITY AND IMPLEMENTATION

We made LiPLike available to the community as a Python toolbox, available at https://gitlab.com/Gustafsson-lab/liplike. We believe that LiPLike will be used for high confidence predictions in studies where individual model interactions are of high importance, and to remove false positive predictions made by other state-of-the-art gene-gene regulation prediction tools.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Charting the Fragmented Landscape of Drug Synergy

Even as the clinical impact of drug combinations continues to accelerate, no consensus on how to quantify drug synergy has emerged. Rather, surveying the landscape of drug synergy reveals the persistence of historical fissures regarding the appropriate domains of conflicting synergy models - fissures impacting all aspects of combination therapy discovery and deployment. Herein we chronicle the impact of these divisions on: (i) the design, interpretation, and reproducibility of high-throughput combination screens; (ii) the performance of algorithms to predict synergistic mixtures; and (iii) the search for higher-order synergistic interactions. Further progress in each of these subfields hinges on reaching a consensus regarding the long-standing rifts in the field.

Introducing Serine as Cardiovascular Disease Biomarker Candidate via Pathway Analysis

Background:

The rate of death due to cardiovascular disease (CVD) is growing. Investigations about CVD that leading to introduce varieties of metabolites is available. The monitoring of these metabolites to find effective ones in the future of clinic applications is the main aim of this study.

Materials and Methods:

Numbers of 34 metabolites for the CVD are extracted from literature and designated for interaction determinations by MetScape V 3.1.3. The compound-reaction-enzyme-gene network was constructed and the pathways were analyzed. Based on the presence of metabolites in the pathways the critical compounds were determined.

Results:

Pathway analysis revealed 18 disturbed pathways related to the CVD. glycerophospholipid metabolism pathway including 27 compounds is related to the 9 queried metabolites. L-Serine which was communed between 5 pathways and also was presented in the largest pathway was identified as the critical compound.

Conclusion:

It can be concluded that L-Serine is a proper biomarker candidate for CVD diagnosis and also patients follow up approaches.

A Guide to Conquer the Biological Network Era Using Graph Theory

Networks are one of the most common ways to represent biological systems as complex sets of binary interactions or relations between different bioentities. In this article, we discuss the basic graph theory concepts and the various graph types, as well as the available data structures for storing and reading graphs. In addition, we describe several network properties and we highlight some of the widely used network topological features. We briefly mention the network patterns, motifs and models, and we further comment on the types of biological and biomedical networks along with their corresponding computer- and human-readable file formats. Finally, we discuss a variety of algorithms and metrics for network analyses regarding graph drawing, clustering, visualization, link prediction, perturbation, and network alignment as well as the current state-of-the-art tools. We expect this review to reach a very broad spectrum of readers varying from experts to beginners while encouraging them to enhance the field further.

The exposome and health: Where chemistry meets biology

Despite extensive evidence showing that exposure to specific chemicals can lead to disease, current research approaches and regulatory policies fail to address the chemical complexity of our world. To safeguard current and future generations from the increasing number of chemicals polluting our environment, a systematic and agnostic approach is needed. The "exposome" concept strives to capture the diversity and range of exposures to synthetic chemicals, dietary constituents, psychosocial stressors, and physical factors, as well as their corresponding biological responses. Technological advances such as high-resolution mass spectrometry and network science have allowed us to take the first steps toward a comprehensive assessment of the exposome. Given the increased recognition of the dominant role that nongenetic factors play in disease, an effort to characterize the exposome at a scale comparable to that of the human genome is warranted.

Analysis of Laser Therapy Effects on Squamous Cell Carcinoma Patients: A System Biology Study

Introduction: The Mechanism of laser therapy and also its safety are 2 important features of the application of different types of lasers in medicine. This study aims to investigate the critically affected genes after the treatment of squamous cell carcinoma patients. Methods: The gene expression profiles of 4 squamous cell carcinoma patients that were treated via chemoradiotherapy (CRT) plus the laser and 3 similar patients without laser exposure from Gene Expression Omnibus (GEO) were downloaded and were screened to find critical genes via network analysis. The STRING database, Cytoscape software, and the Clue GO plug-in of Cytoscape software were used. Results: The genes HSX70 and NCC27 were determined as neighbors and HSPA1B, CLIC1, RAB13, PPIF, and LCE3D as hub genes. The over-expression of LCE3D was interpreted as the side effect of laser therapy. Apoptosis and the cell cycle were the dominant biological processes regulated by the HSP molecules in the laser-treated patients. Conclusion: The laser affected the main biological processes and simultaneously issued side effects.

Understanding allergic multimorbidity within the non-eosinophilic interactome

BACKGROUND

The mechanisms explaining multimorbidity between asthma, dermatitis and rhinitis (allergic multimorbidity) are not well known. We investigated these mechanisms and their specificity in distinct cell types by means of an interactome-based analysis of expression data.

METHODS

Genes associated to the diseases were identified using data mining approaches, and their multimorbidity mechanisms in distinct cell types were characterized by means of an in silico analysis of the topology of the human interactome.

RESULTS

We characterized specific pathomechanisms for multimorbidities between asthma, dermatitis and rhinitis for distinct emergent non-eosinophilic cell types. We observed differential roles for cytokine signaling, TLR-mediated signaling and metabolic pathways for multimorbidities across distinct cell types. Furthermore, we also identified individual genes potentially associated to multimorbidity mechanisms.

CONCLUSIONS

Our results support the existence of differentiated multimorbidity mechanisms between asthma, dermatitis and rhinitis at cell type level, as well as mechanisms common to distinct cell types. These results will help understanding the biology underlying allergic multimorbidity, assisting in the design of new clinical studies.

SourceSet: A graphical model approach to identify primary genes in perturbed biological pathways

Topological gene-set analysis has emerged as a powerful means for omic data interpretation. Although numerous methods for identifying dysregulated genes have been proposed, few of them aim to distinguish genes that are the real source of perturbation from those that merely respond to the signal dysregulation. Here, we propose a new method, called SourceSet, able to distinguish between the primary and the secondary dysregulation within a Gaussian graphical model context. The proposed method compares gene expression profiles in the control and in the perturbed condition and detects the differences in both the mean and the covariance parameters with a series of likelihood ratio tests. The resulting evidence is used to infer the primary and the secondary set, i.e. the genes responsible for the primary dysregulation, and the genes affected by the perturbation through network propagation. The proposed method demonstrates high specificity and sensitivity in different simulated scenarios and on several real biological case studies. In order to fit into the more traditional pathway analysis framework, SourceSet R package also extends the analysis from a single to multiple pathways and provides several graphical outputs, including Cytoscape visualization to browse the results.

The rapid increase in omic studies has created a need to understand the biological implications of their results. Gene-set analysis has emerged as a powerful means for gaining such understanding, evolving in the last decade from the classical enrichment analysis to the more powerful topological approaches. Although numerous methods for identifying dysregulated genes have been proposed, few of them aim to distinguish genes that are the real source of perturbation from those that merely respond to the signal dysregulation. This distinction is crucial for network medicine, where the prioritization of the effect of biological perturbations may help in the molecular understanding of drug treatments and diseases. Here we propose a new method, called SourceSet, able to distinguish between primary and secondary dysregulation within a graphical model context, demonstrating a high specificity and sensitivity in different simulated scenarios and on real biological case studies.

Signal flow control of complex signaling networks

Complex disease such as cancer is often caused by genetic mutations that eventually alter the signal flow in the intra-cellular signaling network and result in different cell fate. Therefore, it is crucial to identify control targets that can most effectively block such unwanted signal flow. For this purpose, systems biological analysis provides a useful framework, but mathematical modeling of complicated signaling networks requires massive time-series measurements of signaling protein activity levels for accurate estimation of kinetic parameter values or regulatory logics. Here, we present a novel method, called SFC (Signal Flow Control), for identifying control targets without the information of kinetic parameter values or regulatory logics. Our method requires only the structural information of a signaling network and is based on the topological estimation of signal flow through the network. SFC will be particularly useful for a large-scale signaling network to which parameter estimation or inference of regulatory logics is no longer applicable in practice. The identified control targets have significant implication in drug development as they can be putative drug targets.

Genetic interaction networks mediate individual statin drug response in Saccharomyces cerevisiae

Eukaryotic genetic interaction networks (GINs) are extensively described in the Saccharomyces cerevisiae S288C model using deletion libraries, yet being limited to this one genetic background, not informative to individual drug response. Here we created deletion libraries in three additional genetic backgrounds. Statin response was probed with five queries against four genetic backgrounds. The 20 resultant GINs representing drug-gene and gene-gene interactions were not conserved by functional enrichment, hierarchical clustering, and topology-based community partitioning. An unfolded protein response (UPR) community exhibited genetic background variation including different betweenness genes that were network bottlenecks, and we experimentally validated this UPR community via measurements of the UPR that were differentially activated and regulated in statin-resistant strains relative to the statin-sensitive S288C background. These network analyses by topology and function provide insight into the complexity of drug response influenced by genetic background.

LCN2-interacting proteins and their expression patterns in brain tumors

Lipocalin 2 (LCN2) is a member of the lipocalin family. Elevated expression of LCN2 has been observed in many human tumors, suggesting it might be a potential biomarker and/or therapeutic target in malignancies. In this study, we aimed to explore LCN2 interacting proteins through bioinformatics, as well as their biological functions. Protein-protein interaction networks (PPIN) were constructed using LCN2 and its interacting proteins as the core node. These PPINs were scale free biological networks in which LCN2 and its interacting proteins could connect or cross-talk with at least one partner protein. Both functional and KEGG pathway enrichment analyses identified the known and potential biological functions of the PPIN, such as cell migration and cancer-related pathways. Expression levels of the PPIN proteins, as well as their expression correlations, in five types of brain tumor, were analyzed and integrated into the PPIN to illustrate a dynamic change. A significant correlation was found between the survival time of glioblastoma patients and the expression level of 10 genes (LCN2, MMP9, MMP2, PDE4DIP, L2HGDH, HNRNPA1, DDX31, LOXL2, FAM60A and RNF25). Taken together, our results suggest that LCN2 and its interacting proteins are mostly differentially expressed and have a distinguishing co-expression pattern. They might promote proliferation and migration via cell migration signaling and cancer-related pathways. LCN2 and its interacting proteins might be potential biomarkers in glioblastoma.

Deciphering the multicomponent synergy mechanisms of SiNiSan prescription on irritable bowel syndrome using a bioinformatics/network topology based strategy

BACKGROUND

SiNiSan (SNS) is a traditional Chinese medicine (TCM) prescription that has been widely used in the clinical treatment of irritable bowel syndrome (IBS). However, the underlying active substances and molecular mechanisms remain obscure.

PURPOSE

A bioinformatics/topology based strategy was proposed for identification of the drug targets, therapeutic agents and molecular mechanisms of SiNiSan against irritable bowel syndrome.

MATERIALS AND METHODS

In this work, a bioinformatics/network topology based strategy was employed by integrating ADME filtering, text mining, bioinformatics, network topology, Venn analysis and molecular docking to uncover systematically the multicomponent synergy mechanisms. In vivo experimental validation was executed in a Visceral Hypersensitivity (VHS) rat model.

RESULTS

76 protein targets and 109 active components of SNS were identified. Bioinformatics analysis revealed that 116 disease pathways associated with IBS therapy could be classified into the 19 statistically enriched functional sub-groups. The multi-functional co-synergism of SNS against IBS were predicted, including inflammatory reaction regulation, oxidative-stress depression regulation and hormone and immune regulation. The multi-component synergetic effects were also revealed on the herbal combination of SNS. The hub-bottleneck genes of the protein networks including PTGS2, CALM2, NOS2, SLC6A3 and MAOB, MAOA, CREB1 could become potential drug targets and Paeoniflorin, Naringin, Glycyrrhizic acid may be candidate agents. Experimental results showed that the potential mechanisms of SiNiSan treatment involved in the suppression of activation of Dopaminergic synapse and Amphetamine addiction signaling pathways, which are congruent with the prediction by the systematic approach.

CONCLUSION

The integrative investigation based on bioinformatics/network topology strategy may elaborate the multicomponent synergy mechanisms of SNS against IBS and provide the way out to develop new combination medicines for IBS.

Scale-variant topological information for characterizing the structure of complex networks

The structure of real-world networks is usually difficult to characterize owing to the variation of topological scales, the nondyadic complex interactions, and the fluctuations in the network. We aim to address these problems by introducing a general framework using a method based on topological data analysis. By considering the diffusion process at a single specified timescale in a network, we map the network nodes to a finite set of points that contains the topological information of the network at a single scale. Subsequently, we study the shape of these point sets over variable timescales that provide scale-variant topological information, to understand the varying topological scales and the complex interactions in the network. We conduct experiments on synthetic and real-world data to demonstrate the effectiveness of the proposed framework in identifying network models, classifying real-world networks, and detecting transition points in time-evolving networks. Overall, our study presents a unified analysis that can be applied to more complex network structures, as in the case of multilayer and multiplex networks.

Towards perturbation prediction of biological networks using deep learning

The mapping of the physical interactions between biochemical entities enables quantitative analysis of dynamic biological living systems. While developing a precise dynamical model on biological entity interaction is still challenging due to the limitation of kinetic parameter detection of the underlying biological system. This challenge promotes the needs of topology-based models to predict biochemical perturbation patterns. Pure topology-based model, however, is limited on the scale and heterogeneity of biological networks. Here we propose a learning based model that adopts graph convolutional networks to learn the implicit perturbation pattern factors and thus enhance the perturbation pattern prediction on the basic topology model. Our experimental studies on 87 biological models show an average of 73% accuracy on perturbation pattern prediction and outperforms the best topology-based model by 7%, indicating that the graph-driven neural network model is robust and beneficial for accurate prediction of the perturbation spread modeling and giving an inspiration of the implementation of the deep neural networks on biological network modeling.

Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots

Root-associated microbes play a key role in plant performance and productivity, making them important players in agroecosystems. So far, very few studies have assessed the impact of different farming systems on the root microbiota and it is still unclear whether agricultural intensification influences the structure and complexity of microbial communities. We investigated the impact of conventional, no-till, and organic farming on wheat root fungal communities using PacBio SMRT sequencing on samples collected from 60 farmlands in Switzerland. Organic farming harbored a much more complex fungal network with significantly higher connectivity than conventional and no-till farming systems. The abundance of keystone taxa was the highest under organic farming where agricultural intensification was the lowest. We also found a strong negative association (R2 = 0.366; P < 0.0001) between agricultural intensification and root fungal network connectivity. The occurrence of keystone taxa was best explained by soil phosphorus levels, bulk density, pH, and mycorrhizal colonization. The majority of keystone taxa are known to form arbuscular mycorrhizal associations with plants and belong to the orders Glomerales, Paraglomerales, and Diversisporales. Supporting this, the abundance of mycorrhizal fungi in roots and soils was also significantly higher under organic farming. To our knowledge, this is the first study to report mycorrhizal keystone taxa for agroecosystems, and we demonstrate that agricultural intensification reduces network complexity and the abundance of keystone taxa in the root microbiome.

The role that choice of model plays in predictions for epilepsy surgery

Mathematical modelling has been widely used to predict the effects of perturbations to brain networks. An important example is epilepsy surgery, where the perturbation in question is the removal of brain tissue in order to render the patient free of seizures. Different dynamical models have been proposed to represent transitions to ictal states in this context. However, our choice of which mathematical model to use to address this question relies on making assumptions regarding the mechanism that defines the transition from background to the seizure state. Since these mechanisms are unknown, it is important to understand how predictions from alternative dynamical descriptions compare. Herein we evaluate to what extent three different dynamical models provide consistent predictions for the effect of removing nodes from networks. We show that for small, directed, connected networks the three considered models provide consistent predictions. For larger networks, predictions are shown to be less consistent. However consistency is higher in networks that have sufficiently large differences in ictogenicity between nodes. We further demonstrate that heterogeneity in ictogenicity across nodes correlates with variability in the number of connections for each node.

Network Medicine in Pathobiology

The past decade has witnessed exponential growth in the generation of high-throughput human data across almost all known dimensions of biological systems. The discipline of network medicine has rapidly evolved in parallel, providing an unbiased, comprehensive biological framework through which to interrogate and integrate systematically these large-scale, multi-omic data to enhance our understanding of disease mechanisms and to design drugs that reflect a deep knowledge of molecular pathobiology. In this review, we discuss the key principles of network medicine and the human disease network and explore the latest applications of network medicine in this multi-omic era. We also highlight the current conceptual and technological challenges, which serve as exciting opportunities by which to improve and expand the network-based applications beyond the artificial boundaries of the current state of human pathobiology.

Analysis of Topological Parameters of Complex Disease Genes Reveals the Importance of Location in a Biomolecular Network

Network biology and medicine provide unprecedented opportunities and challenges for deciphering disease mechanisms from integrative viewpoints. The disease genes and their products perform their dysfunctions via physical and biochemical interactions in the form of a molecular network. The topological parameters of these disease genes in the interactome are of prominent interest to the understanding of their functionality from a systematic perspective. In this work, we provide a systems biology analysis of the topological features of complex disease genes in an integrated biomolecular network. Firstly, we identify the characteristics of four network parameters in the ten most frequently studied disease genes and identify several specific patterns of their topologies. Then, we confirm our findings in the other disease genes of three complex disorders (i.e., Alzheimer’s disease, diabetes mellitus, and hepatocellular carcinoma). The results reveal that the disease genes tend to have a higher betweenness centrality, a smaller average shortest path length, and a smaller clustering coefficient when compared to normal genes, whereas they have no significant degree prominence. The features highlight the importance of gene location in the integrated functional linkages.

Simplicial closure and higher-order link prediction

Networks provide a powerful formalism for modeling complex systems by using a model of pairwise interactions. But much of the structure within these systems involves interactions that take place among more than two nodes at once-for example, communication within a group rather than person to person, collaboration among a team rather than a pair of coauthors, or biological interaction between a set of molecules rather than just two. Such higher-order interactions are ubiquitous, but their empirical study has received limited attention, and little is known about possible organizational principles of such structures. Here we study the temporal evolution of 19 datasets with explicit accounting for higher-order interactions. We show that there is a rich variety of structure in our datasets but datasets from the same system types have consistent patterns of higher-order structure. Furthermore, we find that tie strength and edge density are competing positive indicators of higher-order organization, and these trends are consistent across interactions involving differing numbers of nodes. To systematically further the study of theories for such higher-order structures, we propose higher-order link prediction as a benchmark problem to assess models and algorithms that predict higher-order structure. We find a fundamental difference from traditional pairwise link prediction, with a greater role for local rather than long-range information in predicting the appearance of new interactions.