Current and future directions in network biology.

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

Centrality measures in psychological networks: A simulation study on identifying effective treatment targets.

The network theory of psychopathology suggests that symptoms in a disorder form a network and that identifying central symptoms within this network might be important for an effective and personalized treatment. However, recent evidence has been inconclusive. We analyzed contemporaneous idiographic networks of depression and anxiety symptoms. Two approaches were compared: a cascade-based attack where symptoms were deactivated in decreasing centrality order, and a normal attack where symptoms were deactivated based on original centrality estimates. Results showed that centrality measures significantly affected the attack's magnitude, particularly the number of components and average path length in both normal and cascade attacks. Degree centrality consistently had the highest impact on the network properties. This study emphasizes the importance of considering centrality measures when identifying treatment targets in psychological networks. Further research is needed to better understand the causal relationships and predictive capabilities of centrality measures in personalized treatments for mental disorders.

The clinical trials puzzle: How network effects limit drug discovery.

The depth of knowledge offered by post-genomic medicine has carried the promise of new drugs, and cures for multiple diseases. To explore the degree to which this capability has materialized, we extract meta-data from 356,403 clinical trials spanning four decades, aiming to offer mechanistic insights into the innovation practices in drug discovery. We find that convention dominates over innovation, as over 96% of the recorded trials focus on previously tested drug targets, and the tested drugs target only 12% of the human interactome. If current patterns persist, it would take 170 years to target all druggable proteins. We uncover two network-based fundamental mechanisms that currently limit target discovery: preferential attachment, leading to the repeated exploration of previously targeted proteins; and local network effects, limiting exploration to proteins interacting with highly explored proteins. We build on these insights to develop a quantitative network-based model to enhance drug discovery in clinical trials.

Noncoding RNAs improve the predictive power of network medicine.

Network medicine has improved the mechanistic understanding of disease, offering quantitative insights into disease mechanisms, comorbidities, and novel diagnostic tools and therapeutic treatments. Yet, most network-based approaches rely on a comprehensive map of protein-protein interactions (PPI), ignoring interactions mediated by noncoding RNAs (ncRNAs). Here, we systematically combine experimentally confirmed binding interactions mediated by ncRNA with PPI, constructing a comprehensive network of all physical interactions in the human cell. We find that the inclusion of ncRNA expands the number of genes in the interactome by 46% and the number of interactions by 107%, significantly enhancing our ability to identify disease modules. Indeed, we find that 132 diseases lacked a statistically significant disease module in the protein-based interactome but have a statistically significant disease module after inclusion of ncRNA-mediated interactions, making these diseases accessible to the tools of network medicine. We show that the inclusion of ncRNAs helps unveil disease-disease relationships that were not detectable before and expands our ability to predict comorbidity patterns between diseases. Taken together, we find that including noncoding interactions improves both the breath and the predictive accuracy of network medicine.

Improving the generalizability of protein-ligand binding predictions with AI-Bind.

Identifying novel drug-target interactions is a critical and rate-limiting step in drug discovery. While deep learning models have been proposed to accelerate the identification process, here we show that state-of-the-art models fail to generalize to novel (i.e., never-before-seen) structures. We unveil the mechanisms responsible for this shortcoming, demonstrating how models rely on shortcuts that leverage the topology of the protein-ligand bipartite network, rather than learning the node features. Here we introduce AI-Bind, a pipeline that combines network-based sampling strategies with unsupervised pre-training to improve binding predictions for novel proteins and ligands. We validate AI-Bind predictions via docking simulations and comparison with recent experimental evidence, and step up the process of interpreting machine learning prediction of protein-ligand binding by identifying potential active binding sites on the amino acid sequence. AI-Bind is a high-throughput approach to identify drug-target combinations with the potential of becoming a powerful tool in drug discovery.

Identification of potent inhibitors of SARS-CoV-2 infection by combined pharmacological evaluation and cellular network prioritization.

Pharmacologically active compounds with known biological targets were evaluated for inhibition of SARS-CoV-2 infection in cell and tissue models to help identify potent classes of active small molecules and to better understand host-virus interactions. We evaluated 6,710 clinical and preclinical compounds targeting 2,183 host proteins by immunocytofluorescence-based screening to identify SARS-CoV-2 infection inhibitors. Computationally integrating relationships between small molecule structure, dose-response antiviral activity, host target, and cell interactome produced cellular networks important for infection. This analysis revealed 389 small molecules with micromolar to low nanomolar activities, representing >12 scaffold classes and 813 host targets. Representatives were evaluated for mechanism of action in stable and primary human cell models with SARS-CoV-2 variants and MERS-CoV. One promising candidate, obatoclax, significantly reduced SARS-CoV-2 viral lung load in mice. Ultimately, this work establishes a rigorous approach for future pharmacological and computational identification of host factor dependencies and treatments for viral diseases.

The nosographic structure of posttraumatic stress symptoms across trauma types: An exploratory network analysis approach.

The nosographic structure of posttraumatic stress disorder (PTSD) remains unclear, and attempts to determine its symptomatic organization have been unsatisfactory. Several explanations have been suggested, and the impact of trauma type is receiving increasing attention. As little is known about the differential impact trauma type in the nosographic structure of PTSD, we explored the nosology of PTSD and the effect of trauma type on its symptomatic organization. We reanalyzed five cross-sectional psychopathological networks involving different trauma types, encompassing a broad range of traumatic events in veterans, war-related trauma in veterans, sexual abuse, terrorist attacks, and various traumatic events in refugees. The weighted topological overlap was used to estimate the networks and attribute weights to their links. Coexpression differential network analysis was used to identify the common and specific network structures of the connections across different trauma types and to determine the importance of symptoms across the networks. We found a set of symptoms with more common connections with other symptoms, suggesting that these might constitute the prototypical nosographic structure of PTSD. We also found a set of symptoms that had a high number of specific connections with other symptoms; these connections varied according to trauma type. The importance of symptoms across the common and specific networks was ascertained. The present findings offer new insights into the symptomatic organization of PTSD and support previous research on the impact of trauma type on the nosology of this disorder.

Identification of druggable host targets needed for SARS-CoV-2 infection by combined pharmacological evaluation and cellular network directed prioritization both in vitro and in vivo.

Identification of host factors contributing to replication of viruses and resulting disease progression remains a promising approach for development of new therapeutics. Here, we evaluated 6710 clinical and preclinical compounds targeting 2183 host proteins by immunocytofluorescence-based screening to identify SARS-CoV-2 infection inhibitors. Computationally integrating relationships between small molecule structure, dose-response antiviral activity, host target and cell interactome networking produced cellular networks important for infection. This analysis revealed 389 small molecules, >12 scaffold classes and 813 host targets with micromolar to low nanomolar activities. From these classes, representatives were extensively evaluated for mechanism of action in stable and primary human cell models, and additionally against Beta and Delta SARS-CoV-2 variants and MERS-CoV. One promising candidate, obatoclax, significantly reduced SARS-CoV-2 viral lung load in mice. Ultimately, this work establishes a rigorous approach for future pharmacological and computational identification of novel host factor dependencies and treatments for viral diseases.

Accelerated Evolution of Tissue-Specific Genes Mediates Divergence Amidst Gene Flow in European Green Lizards.

The European green lizards of the Lacerta viridis complex consist of two closely related species, L. viridis and Lacerta bilineata that split less than 7 million years ago in the presence of gene flow. Recently, a third lineage, referred to as the "Adriatic" was described within the L. viridis complex distributed from Slovenia to Greece. However, whether gene flow between the Adriatic lineage and L. viridis or L. bilineata has occurred and the evolutionary processes involved in their diversification are currently unknown. We hypothesized that divergence occurred in the presence of gene flow between multiple lineages and involved tissue-specific gene evolution. In this study, we sequenced the whole genome of an individual of the Adriatic lineage and tested for the presence of gene flow amongst L. viridis, L. bilineata, and Adriatic. Additionally, we sequenced transcriptomes from multiple tissues to understand tissue-specific effects. The species tree supports that the Adriatic lineage is a sister taxon to L. bilineata. We detected gene flow between the Adriatic lineage and L. viridis suggesting that the evolutionary history of the L. viridis complex is likely shaped by gene flow. Interestingly, we observed topological differences between the autosomal and Z-chromosome phylogenies with a few fast evolving genes on the Z-chromosome. Genes highly expressed in the ovaries and strongly co-expressed in the brain experienced accelerated evolution presumably contributing to establishing reproductive isolation in the L. viridis complex.

Network medicine framework for identifying drug-repurposing opportunities for COVID-19.

The COVID-19 pandemic has highlighted the need to quickly and reliably prioritize clinically approved compounds for their potential effectiveness for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections. Here, we deployed algorithms relying on artificial intelligence, network diffusion, and network proximity, tasking each of them to rank 6,340 drugs for their expected efficacy against SARS-CoV-2. To test the predictions, we used as ground truth 918 drugs experimentally screened in VeroE6 cells, as well as the list of drugs in clinical trials that capture the medical community's assessment of drugs with potential COVID-19 efficacy. We find that no single predictive algorithm offers consistently reliable outcomes across all datasets and metrics. This outcome prompted us to develop a multimodal technology that fuses the predictions of all algorithms, finding that a consensus among the different predictive methods consistently exceeds the performance of the best individual pipelines. We screened in human cells the top-ranked drugs, obtaining a 62% success rate, in contrast to the 0.8% hit rate of nonguided screenings. Of the six drugs that reduced viral infection, four could be directly repurposed to treat COVID-19, proposing novel treatments for COVID-19. We also found that 76 of the 77 drugs that successfully reduced viral infection do not bind the proteins targeted by SARS-CoV-2, indicating that these network drugs rely on network-based mechanisms that cannot be identified using docking-based strategies. These advances offer a methodological pathway to identify repurposable drugs for future pathogens and neglected diseases underserved by the costs and extended timeline of de novo drug development.

Whole transcriptomic network analysis using Co-expression Differential Network Analysis (CoDiNA).

Biological and medical sciences are increasingly acknowledging the significance of gene co-expression-networks for investigating complex-systems, phenotypes or diseases. Typically, complex phenotypes are investigated under varying conditions. While approaches for comparing nodes and links in two networks exist, almost no methods for the comparison of multiple networks are available and-to best of our knowledge-no comparative method allows for whole transcriptomic network analysis. However, it is the aim of many studies to compare networks of different conditions, for example, tissues, diseases, treatments, time points, or species. Here we present a method for the systematic comparison of an unlimited number of networks, with unlimited number of transcripts: Co-expression Differential Network Analysis (CoDiNA). In particular, CoDiNA detects links and nodes that are common, specific or different among the networks. We developed a statistical framework to normalize between these different categories of common or changed network links and nodes, resulting in a comprehensive network analysis method, more sophisticated than simply comparing the presence or absence of network nodes. Applying CoDiNA to a neurogenesis study we identified candidate genes involved in neuronal differentiation. We experimentally validated one candidate, demonstrating that its overexpression resulted in a significant disturbance in the underlying gene regulatory network of neurogenesis. Using clinical studies, we compared whole transcriptome co-expression networks from individuals with or without HIV and active tuberculosis (TB) and detected signature genes specific to HIV. Furthermore, analyzing multiple cancer transcription factor (TF) networks, we identified common and distinct features for particular cancer types. These CoDiNA applications demonstrate the successful detection of genes associated with specific phenotypes. Moreover, CoDiNA can also be used for comparing other types of undirected networks, for example, metabolic, protein-protein interaction, ecological and psychometric networks. CoDiNA is publicly available as an R package in CRAN (https://CRAN.R-project.org/package=CoDiNA).

Seasonal Dynamics of Pelagic Mycoplanktonic Communities: Interplay of Taxon Abundance, Temporal Occurrence, and Biotic Interactions.

Marine fungi are an important component of pelagic planktonic communities. However, it is not yet clear how individual fungal taxa are integrated in marine processes of the microbial loop and food webs. Most likely, biotic interactions play a major role in shaping the fungal community structure. Thus, the aim of our work was to identify possible biotic interactions of mycoplankton with phytoplankton and zooplankton groups and among fungi, and to investigate whether there is coherence between interactions and the dynamics, abundance and temporal occurrence of individual fungal OTUs. Marine surface water was sampled weekly over the course of 1 year, in the vicinity of the island of Helgoland in the German Bight (North Sea). The mycoplankton community was analyzed using 18S rRNA gene tag-sequencing and the identified dynamics were correlated to environmental data including phytoplankton, zooplankton, and abiotic factors. Finally, co-occurrence patterns of fungal taxa were detected with network analyses based on weighted topological overlaps (wTO). Of all abundant and persistent OTUs, 77% showed no biotic relations suggesting a saprotrophic lifestyle. Of all other fungal OTUs, nearly the half (44%) had at least one significant negative relationship, especially with zooplankton and other fungi, or to a lesser extent with phytoplankton. These findings suggest that mycoplankton OTUs are embedded into marine food web chains via highly complex and manifold relationships such as parasitism, predation, grazing, or allelopathy. Furthermore, about one third of all rare OTUs were part of a dense fungal co-occurrence network probably stabilizing the fungal community against environmental changes and acting as functional guilds or being involved in fungal cross-feeding. Placed in an ecological context, strong antagonistic relationships of the mycoplankton community with other components of the plankton suggest that: (i) there is a top-down control by fungi on zooplankton and phytoplankton; (ii) fungi serve as a food source for zooplankton and thereby transfer nutrients and organic material; (iii) the dynamics of fungi harmful to other plankton groups are controlled by antagonistic fungal taxa.

Network Medicine Framework for Identifying Drug Repurposing Opportunities for COVID-19.

The current pandemic has highlighted the need for methodologies that can quickly and reliably prioritize clinically approved compounds for their potential effectiveness for SARS-CoV-2 infections. In the past decade, network medicine has developed and validated multiple predictive algorithms for drug repurposing, exploiting the sub-cellular network-based relationship between a drug's targets and disease genes. Here, we deployed algorithms relying on artificial intelligence, network diffusion, and network proximity, tasking each of them to rank 6,340 drugs for their expected efficacy against SARS-CoV-2. To test the predictions, we used as ground truth 918 drugs that had been experimentally screened in VeroE6 cells, and the list of drugs under clinical trial, that capture the medical community's assessment of drugs with potential COVID-19 efficacy. We find that while most algorithms offer predictive power for these ground truth data, no single method offers consistently reliable outcomes across all datasets and metrics. This prompted us to develop a multimodal approach that fuses the predictions of all algorithms, showing that a consensus among the different predictive methods consistently exceeds the performance of the best individual pipelines. We find that 76 of the 77 drugs that successfully reduced viral infection do not bind the proteins targeted by SARS-CoV-2, indicating that these drugs rely on network-based actions that cannot be identified using docking-based strategies. These advances offer a methodological pathway to identify repurposable drugs for future pathogens and neglected diseases underserved by the costs and extended timeline of de novo drug development.

Construction, comparison and evolution of networks in life sciences and other disciplines.

Network approaches have become pervasive in many research fields. They allow for a more comprehensive understanding of complex relationships between entities as well as their group-level properties and dynamics. Many networks change over time, be it within seconds or millions of years, depending on the nature of the network. Our focus will be on comparative network analyses in life sciences, where deciphering temporal network changes is a core interest of molecular, ecological, neuropsychological and evolutionary biologists. Further, we will take a journey through different disciplines, such as social sciences, finance and computational gastronomy, to present commonalities and differences in how networks change and can be analysed. Finally, we envision how borrowing ideas from these disciplines could enrich the future of life science research.

Honey bee virus causes context-dependent changes in host social behavior.

Anthropogenic changes create evolutionarily novel environments that present opportunities for emerging diseases, potentially changing the balance between host and pathogen. Honey bees provide essential pollination services, but intensification and globalization of honey bee management has coincided with increased pathogen pressure, primarily due to a parasitic mite/virus complex. Here, we investigated how honey bee individual and group phenotypes are altered by a virus of concern, Israeli acute paralysis virus (IAPV). Using automated and manual behavioral monitoring of IAPV-inoculated individuals, we find evidence for pathogen manipulation of worker behavior by IAPV, and reveal that this effect depends on social context; that is, within versus between colony interactions. Experimental inoculation reduced social contacts between honey bee colony members, suggesting an adaptive host social immune response to diminish transmission. Parallel analyses with double-stranded RNA (dsRNA)-immunostimulated bees revealed these behaviors are part of a generalized social immune defensive response. Conversely, inoculated bees presented to groups of bees from other colonies experienced reduced aggression compared with dsRNA-immunostimulated bees, facilitating entry into susceptible colonies. This reduction was associated with a shift in cuticular hydrocarbons, the chemical signatures used by bees to discriminate colony members from intruders. These responses were specific to IAPV infection, suggestive of pathogen manipulation of the host. Emerging bee pathogens may thus shape host phenotypes to increase transmission, a strategy especially well-suited to the unnaturally high colony densities of modern apiculture. These findings demonstrate how anthropogenic changes could affect arms races between human-managed hosts and their pathogens to potentially affect global food security.

wTO: an R package for computing weighted topological overlap and a consensus network with integrated visualization tool.

BACKGROUND: Network analyses, such as of gene co-expression networks, metabolic networks and ecological networks have become a central approach for the systems-level study of biological data. Several software packages exist for generating and analyzing such networks, either from correlation scores or the absolute value of a transformed score called weighted topological overlap (wTO). However, since gene regulatory processes can up- or down-regulate genes, it is of great interest to explicitly consider both positive and negative correlations when constructing a gene co-expression network. RESULTS: Here, we present an R package for calculating the weighted topological overlap (wTO), that, in contrast to existing packages, explicitly addresses the sign of the wTO values, and is thus especially valuable for the analysis of gene regulatory networks. The package includes the calculation of p-values (raw and adjusted) for each pairwise gene score. Our package also allows the calculation of networks from time series (without replicates). Since networks from independent datasets (biological repeats or related studies) are not the same due to technical and biological noise in the data, we additionally, incorporated a novel method for calculating a consensus network (CN) from two or more networks into our R package. To graphically inspect the resulting networks, the R package contains a visualization tool, which allows for the direct network manipulation and access of node and link information. When testing the package on a standard laptop computer, we can conduct all calculations for systems of more than 20,000 genes in under two hours. We compare our new wTO package to state of art packages and demonstrate the application of the wTO and CN functions using 3 independently derived datasets from healthy human pre-frontal cortex samples. To showcase an example for the time series application we utilized a metagenomics data set. CONCLUSION: In this work, we developed a software package that allows the computation of wTO networks, CNs and a visualization tool in the R statistical environment. It is publicly available on CRAN repositories under the GPL -2 Open Source License ( https://cran.r-project.org/web/packages/wTO/ ).

Combined Experimental and System-Level Analyses Reveal the Complex Regulatory Network of miR-124 during Human Neurogenesis.

Non-coding RNAs regulate many biological processes including neurogenesis. The brain-enriched miR-124 has been assigned as a key player of neuronal differentiation via its complex but little understood regulation of thousands of annotated targets. To systematically chart its regulatory functions, we used CRISPR/Cas9 gene editing to disrupt all six miR-124 alleles in human induced pluripotent stem cells. Upon neuronal induction, miR-124-deleted cells underwent neurogenesis and became functional neurons, albeit with altered morphology and neurotransmitter specification. Using RNA-induced-silencing-complex precipitation, we identified 98 high-confidence miR-124 targets, of which some directly led to decreased viability. By performing advanced transcription-factor-network analysis, we identified indirect miR-124 effects on apoptosis, neuronal subtype differentiation, and the regulation of previously uncharacterized zinc finger transcription factors. Our data emphasize the need for combined experimental- and system-level analyses to comprehensively disentangle and reveal miRNA functions, including their involvement in the neurogenesis of diverse neuronal cell types found in the human brain.

Measuring hand hygiene compliance rates at hospital entrances.

BACKGROUND: Despite the importance of hand hygiene in the health care setting, there are no studies evaluating hand hygiene compliance at hospital entrances. METHODS: The study was prospectively performed over a 33-week period from March 30, 2014-November 15, 2014, to evaluate hand hygiene compliance in 2 hospital reception areas. We compared electronic handwash counters with the application of radiofrequency identification (GOJO SMARTLINK) (electronic observer) that counts each activation of alcohol gel dispensers to direct observation (human observer) via remote review of video surveillance. RESULTS: We found low hand hygiene compliance rates of 2.2% (99/4,412) and 1.7% (140/8,277), respectively, at reception areas A and D, detected by direct observation. Using the electronic observer, we measured rates of 17% (15,624/91,724) and 7.1% (51,605/730,357) at reception areas A and D, respectively. For the overall time period of simultaneous electronic and human observation, the human observer captured 1% of the hand hygiene episodes detected by the electronic observer. CONCLUSIONS: Our study showed very low hand hygiene compliance in hospital reception areas, and we found an electronic hand hygiene system to be a useful method to monitor hand hygiene compliance.