Multilayer networks: aspects, implementations, and application in biomedicine

Multilayer biological networks to upscale marine research to global change-smart management and sustainable resource use

Human activities are having a massive negative impact on biodiversity and ecological processes worldwide. The rate and magnitude of ecological transformations induced by climate change, habitat destruction, overexploitation and pollution are now so substantial that a sixth mass extinction event is currently underway. The biodiversity crisis of the Anthropocene urges scientists to put forward a transformative vision to promote the conservation of biodiversity, and thus indirectly the preservation of ecosystem functions. Here, we identify pressing issues in global change biology research and propose an integrative framework based on multilayer biological networks as a tool to support conservation actions and marine risk assessments in multi-stressor scenarios. Multilayer networks can integrate different levels of environmental and biotic complexity, enabling us to combine information on molecular, physiological and behaviour responses, species interactions and biotic communities. The ultimate aim of this framework is to link human-induced environmental changes to species physiology, fitness, biogeography and ecosystem impacts across vast seascapes and time frames, to help guide solutions to address biodiversity loss and ecological tipping points. Further, we also define our current ability to adopt a widespread use of multilayer networks within ecology, evolution and conservation by providing examples of case-studies. We also assess which approaches are ready to be transferred and which ones require further development before use. We conclude that multilayer biological networks will be crucial to inform (using reliable multi-levels integrative indicators) stakeholders and support their decision-making concerning the sustainable use of resources and marine conservation.

The Molecular Landscape of Premature Aging Diseases Defined by Multilayer Network Exploration

Premature Aging (PA) diseases are rare genetic disorders that mimic some aspects of physiological aging at an early age. Various causative genes of PA diseases have been identified in recent years, providing insights into some dysfunctional cellular processes. However, the identification of PA genes also revealed significant genetic heterogeneity and highlighted the gaps in this understanding of PA-associated molecular mechanisms. Furthermore, many patients remain undiagnosed. Overall, the current lack of knowledge about PA diseases hinders the development of effective diagnosis and therapies and poses significant challenges to improving patient care. Here, a network-based approach to systematically unravel the cellular functions disrupted in PA diseases is presented. Leveraging a network community identification algorithm, it is delved into a vast multilayer network of biological interactions to extract the communities of 67 PA diseases from their 132 associated genes. It is found that these communities can be grouped into six distinct clusters, each reflecting specific cellular functions: DNA repair, cell cycle, transcription regulation, inflammation, cell communication, and vesicle-mediated transport. That these clusters collectively represent the landscape of the molecular mechanisms that are perturbed in PA diseases, providing a framework for better understanding their pathogenesis is proposed. Intriguingly, most clusters also exhibited a significant enrichment in genes associated with physiological aging, suggesting a potential overlap between the molecular underpinnings of PA diseases and natural aging.

Dataset on global trade networks of COVID-19 medical products

This paper presents a comprehensive dataset on the global trade dynamics of COVID-19-related medical products for the years 2019 and 2020. The dataset, derived from the BACI database, focuses on eight distinct product categories identified by six-digit codes. The trade flow data for 224 countries is structured as a multilevel network, with countries as nodes and product categories as layers. Directed edges represent trading activities, and edge weights are determined by the difference in exported values between 2019 and 2020. The dataset is provided in an edges-and-nodes format. Additionally, the associated R script transforms the data into the MuxViz R package format, facilitating further analysis and visualization of the dataset. The dataset is valuable for researchers in the field of foreign trade or medical products, and for decision-makers in these fields, whether at company or national level.

The spatial distribution of coupling between tau and neurodegeneration in amyloid-β positive mild cognitive impairment

Synergies between amyloid-β (Aβ), tau, and neurodegeneration persist along the Alzheimer's disease (AD) continuum. This study aimed to evaluate the extent of spatial coupling between tau and neurodegeneration (atrophy) and its relation to Aβ positivity in mild cognitive impairment (MCI). Data from 409 participants were included (95 cognitively normal controls, 158 Aβ positive (Aβ+) MCI, and 156 Aβ negative (Aβ-) MCI). Florbetapir PET, Flortaucipir PET, and structural MRI were used as biomarkers for Aβ, tau and atrophy, respectively. Individual correlation matrices for tau load and atrophy were used to layer a multilayer network, with separate layers for tau and atrophy. A measure of coupling between corresponding regions of interest (ROIs) in the tau and atrophy layers was computed, as a function of Aβ positivity. Fewer than 25% of the ROIs across the brain showed heightened coupling between tau and atrophy in Aβ+ , relative to Aβ- MCI. Coupling strengths in the right rostral middle frontal and right paracentral gyri, in particular, mediated the association between Aβ burden and cognition in this sample.

Reconstructing the post-glacial spread of the sand fly Phlebotomus mascittii Grassi, 1908 (Diptera: Psychodidae) in Europe

Phlebotomine sand flies (Diptera: Phlebotominae) are the principal vectors of Leishmania spp. (Kinetoplastida: Trypanosomatidae). In Central Europe, Phlebotomus mascittii is the predominant species, but largely understudied. To better understand factors driving its current distribution, we infer patterns of genetic diversity by testing for signals of population expansion based on two mitochondrial genes and model current and past climate and habitat suitability for seven post-glacial maximum periods, taking 19 climatic variables into account. Consequently, we elucidate their connections by environmental-geographical network analysis. Most analyzed populations share a main haplotype tracing back to a single glacial maximum refuge area on the Mediterranean coasts of South France, which is supported by network analysis. The rapid range expansion of Ph. mascittii likely started in the early mid-Holocene epoch until today and its spread possibly followed two routes. The first one was through northern France to Germany and then Belgium, and the second across the Ligurian coast through present-day Slovenia to Austria, toward the northern Balkans. Here we present a combined approach to reveal glacial refugia and post-glacial spread of Ph. mascittii and observed discrepancies between the modelled and the current known distribution might reveal yet overlooked populations and potential further spread.

Multilayer network alignment based on topological assessment via embeddings

BACKGROUND

Network graphs allow modelling the real world objects in terms of interactions. In a multilayer network, the interactions are distributed over layers (i.e., intralayer and interlayer edges). Network alignment (NA) is a methodology that allows mapping nodes between two or multiple given networks, by preserving topologically similar regions. For instance, NA can be applied to transfer knowledge from one biological species to another. In this paper, we present DANTEml, a software tool for the Pairwise Global NA (PGNA) of multilayer networks, based on topological assessment. It builds its own similarity matrix by processing the node embeddings computed from two multilayer networks of interest, to evaluate their topological similarities. The proposed solution can be used via a user-friendly command line interface, also having a built-in guided mode (step-by-step) for defining input parameters.

RESULTS

We investigated the performance of DANTEml based on (i) performance evaluation on synthetic multilayer networks, (ii) statistical assessment of the resulting alignments, and (iii) alignment of real multilayer networks. DANTEml over performed a method that does not consider the distribution of nodes and edges over multiple layers by 1193.62%, and a method for temporal NA by 25.88%; we also performed the statistical assessment, which corroborates the significance of its own node mappings. In addition, we tested the proposed solution by using a real multilayer network in presence of several levels of noise, in accordance with the same outcome pursued for the NA on our dataset of synthetic networks. In this case, the improvement is even more evident: +4008.75% and +111.72%, compared to a method that does not consider the distribution of nodes and edges over multiple layers and a method for temporal NA, respectively.

CONCLUSIONS

DANTEml is a software tool for the PGNA of multilayer networks based on topological assessment, that is able to provide effective alignments both on synthetic and real multi layer networks, of which node mappings can be validated statistically. Our experimentation reported a high degree of reliability and effectiveness for the proposed solution.

Gene communities in co-expression networks across different tissues

With the recent availability of tissue-specific gene expression data, e.g., provided by the GTEx Consortium, there is interest in comparing gene co-expression patterns across tissues. One promising approach to this problem is to use a multilayer network analysis framework and perform multilayer community detection. Communities in gene co-expression networks reveal groups of genes similarly expressed across individuals, potentially involved in related biological processes responding to specific environmental stimuli or sharing common regulatory variations. We construct a multilayer network in which each of the four layers is an exocrine gland tissue-specific gene co-expression network. We develop methods for multilayer community detection with correlation matrix input and an appropriate null model. Our correlation matrix input method identifies five groups of genes that are similarly co-expressed in multiple tissues (a community that spans multiple layers, which we call a generalist community) and two groups of genes that are co-expressed in just one tissue (a community that lies primarily within just one layer, which we call a specialist community). We further found gene co-expression communities where the genes physically cluster across the genome significantly more than expected by chance (on chromosomes 1 and 11). This clustering hints at underlying regulatory elements determining similar expression patterns across individuals and cell types. We suggest that KRTAP3-1, KRTAP3-3, and KRTAP3-5 share regulatory elements in skin and pancreas. Furthermore, we find that CELA3A and CELA3B share associated expression quantitative trait loci in the pancreas. The results indicate that our multilayer community detection method for correlation matrix input extracts biologically interesting communities of genes.

Network medicine for patients' stratification: From single-layer to multi-omics

Precision medicine research increasingly relies on the integrated analysis of multiple types of omics. In the era of big data, the large availability of different health-related information represents a great, but at the same time untapped, chance with a potentially fundamental role in the prevention, diagnosis and prognosis of diseases. Computational methods are needed to combine this data to create a comprehensive view of a given disease. Network science can model biomedical data in terms of relationships among molecular players of different nature and has been successfully proposed as a new paradigm for studying human diseases. Patient stratification is an open challenge aimed at identifying subtypes with different disease manifestations, severity, and expected survival time. Several stratification approaches based on high-throughput gene expression measurements have been successfully applied. However, few attempts have been proposed to exploit the integration of various genotypic and phenotypic data to discover novel sub-types or improve the detection of known groupings. This article is categorized under: Cancer > Biomedical Engineering Cancer > Computational Models Cancer > Genetics/Genomics/Epigenetics.

PLEX.I: a tool to discover features in multiplex networks that reflect clinical variation

Molecular profiling technologies, such as RNA sequencing, offer new opportunities to better discover and understand the molecular networks involved in complex biological processes. Clinically important variations of diseases, or responses to treatment, are often reflected, or even caused, by the dysregulation of molecular interaction networks specific to particular network regions. In this work, we propose the R package PLEX.I, that allows quantifying and testing variation in the direct neighborhood of a given node between networks corresponding to different conditions or states. We illustrate PLEX.I in two applications in which we discover variation that is associated with different responses to tamoxifen treatment and to sex-specific responses to bacterial stimuli. In the first case, PLEX.I analysis identifies two known pathways i) that have already been implicated in the same context as the tamoxifen mechanism of action, and ii) that would have not have been identified using classical differential gene expression analysis.

An Exploratory Application of Multilayer Networks and Pathway Analysis in Pharmacogenomics

Over the years, network analysis has become a promising strategy for analysing complex system, i.e., systems composed of a large number of interacting elements. In particular, multilayer networks have emerged as a powerful framework for modelling and analysing complex systems with multiple types of interactions. Network analysis can be applied to pharmacogenomics to gain insights into the interactions between genes, drugs, and diseases. By integrating network analysis techniques with pharmacogenomic data, the goal consists of uncovering complex relationships and identifying key genes to use in pathway enrichment analysis to figure out biological pathways involved in drug response and adverse reactions. In this study, we modelled omics, disease, and drug data together through multilayer network representation. Then, we mined the multilayer network with a community detection algorithm to obtain the top communities. After that, we used the identified list of genes from the communities to perform pathway enrichment analysis (PEA) to figure out the biological function affected by the selected genes. The results show that the genes forming the top community have multiple roles through different pathways.

Systems Biology and Cytokines Potential Role in Lung Cancer Immunotherapy Targeting Autophagic Axis

Lung cancer accounts for the highest number of deaths among men and women worldwide. Although extensive therapies, either alone or in conjunction with some specific drugs, continue to be the principal regimen for evolving lung cancer, significant improvements are still needed to understand the inherent biology behind progressive inflammation and its detection. Unfortunately, despite every advancement in its treatment, lung cancer patients display different growth mechanisms and continue to die at significant rates. Autophagy, which is a physiological defense mechanism, serves to meet the energy demands of nutrient-deprived cancer cells and sustain the tumor cells under stressed conditions. In contrast, autophagy is believed to play a dual role during different stages of tumorigenesis. During early stages, it acts as a tumor suppressor, degrading oncogenic proteins; however, during later stages, autophagy supports tumor cell survival by minimizing stress in the tumor microenvironment. The pivotal role of the IL6-IL17-IL23 signaling axis has been observed to trigger autophagic events in lung cancer patients. Since the obvious roles of autophagy are a result of different immune signaling cascades, systems biology can be an effective tool to understand these interconnections and enhance cancer treatment and immunotherapy. In this review, we focus on how systems biology can be exploited to target autophagic processes that resolve inflammatory responses and contribute to better treatment in carcinogenesis.

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.

Capturing the dynamics of microbial interactions through individual-specific networks

Longitudinal analysis of multivariate individual-specific microbiome profiles over time or across conditions remains dauntin. Most statistical tools and methods that are available to study microbiomes are based on cross-sectional data. Over the past few years, several attempts have been made to model the dynamics of bacterial species over time or across conditions. However, the field needs novel views on handling microbial interactions in temporal analyses. This study proposes a novel data analysis framework, MNDA, that combines representation learning and individual-specific microbial co-occurrence networks to uncover taxon neighborhood dynamics. As a use case, we consider a cohort of newborns with microbiomes available at 6 and 9 months after birth, and extraneous data available on the mode of delivery and diet changes between the considered time points. Our results show that prediction models for these extraneous outcomes based on an MNDA measure of local neighborhood dynamics for each taxon outperform traditional prediction models solely based on individual-specific microbial abundances. Furthermore, our results show that unsupervised similarity analysis of newborns in the study, again using the notion of a taxon's dynamic neighborhood derived from time-matched individual-specific microbial networks, can reveal different subpopulations of individuals, compared to standard microbiome-based clustering, with potential relevance to clinical practice. This study highlights the complementarity of microbial interactions and abundances in downstream analyses and opens new avenues to personalized prediction or stratified medicine with temporal microbiome data.

MultiCens: Multilayer network centrality measures to uncover molecular mediators of tissue-tissue communication

With the evolution of multicellularity, communication among cells in different tissues and organs became pivotal to life. Molecular basis of such communication has long been studied, but genome-wide screens for genes and other biomolecules mediating tissue-tissue signaling are lacking. To systematically identify inter-tissue mediators, we present a novel computational approach MultiCens (Multilayer/Multi-tissue network Centrality measures). Unlike single-layer network methods, MultiCens can distinguish within- vs. across-layer connectivity to quantify the "influence" of any gene in a tissue on a query set of genes of interest in another tissue. MultiCens enjoys theoretical guarantees on convergence and decomposability, and performs well on synthetic benchmarks. On human multi-tissue datasets, MultiCens predicts known and novel genes linked to hormones. MultiCens further reveals shifts in gene network architecture among four brain regions in Alzheimer's disease. MultiCens-prioritized hypotheses from these two diverse applications, and potential future ones like "Multi-tissue-expanded Gene Ontology" analysis, can enable whole-body yet molecular-level systems investigations in humans.

Network modeling helps to tackle the complexity of drug-disease systems

From the (patho)physiological point of view, diseases can be considered as emergent properties of living systems stemming from the complexity of these systems. Complex systems display some typical features, including the presence of emergent behavior and the organization in successive hierarchic levels. Drug treatments increase this complexity scenario, and from some years the use of network models has been introduced to describe drug-disease systems and to make predictions about them with regard to several aspects related to drug discovery. Here, we review some recent examples thereof with the aim to illustrate how network science tools can be very effective in addressing both tasks. We will examine the use of bipartite networks that lead to the important concept of "disease module", as well as the introduction of more articulated models, like multi-scale and multiplex networks, able to describe disease systems at increasing levels of organization. Examples of predictive models will then be discussed, considering both those that exploit approaches purely based on graph theory and those that integrate machine learning methods. A short account of both kinds of methodological applications will be provided. Finally, the point will be made on the present situation of modeling complex drug-disease systems highlighting some open issues. This article is categorized under: Neurological Diseases > Computational Models Infectious Diseases > Computational Models Cardiovascular Diseases > Computational Models.

konnect2prot: a web application to explore the protein properties in a functional protein-protein interaction network

MOTIVATION

The regulation of proteins governs the biological processes and functions and, therefore, the organisms' phenotype. So there is an unmet need for a systematic tool for identifying the proteins that play a crucial role in information processing in a protein-protein interaction (PPI) network. However, the current protein databases and web servers still lag behind to provide an end-to-end pipeline that can leverage the topological understanding of a context-specific PPI network to identify the influential spreaders. Addressing this, we developed a web application, 'konnect2prot' (k2p), which can generate context-specific directional PPI network from the input proteins and detect their biological and topological importance in the network.

RESULTS

We pooled together a large amount of ontological knowledge, parsed it down into a functional network, and gained insight into the molecular underpinnings of the disease development by creating a one-stop junction for PPI data. k2p contains both local and global information about a protein, such as protein class, disease mutations, ligands and PDB structure, enriched processes and pathways, multi-disease interactome and hubs and bottlenecks in the directional network. It also identifies spreaders in the network and maps them to disease hallmarks to determine whether they can affect the disease state or not.

AVAILABILITY AND IMPLEMENTATION

konnect2prot is freely accessible using the link https://konnect2prot.thsti.in. The code repository is https://github.com/samrat-lab/k2p_bioinfo-2022.

Metabolite annotation from knowns to unknowns through knowledge-guided multi-layer metabolic networking

Liquid chromatography - mass spectrometry (LC-MS) based untargeted metabolomics allows to measure both known and unknown metabolites in the metabolome. However, unknown metabolite annotation is a major challenge in untargeted metabolomics. Here, we develop an approach, namely, knowledge-guided multi-layer network (KGMN), to enable global metabolite annotation from knowns to unknowns in untargeted metabolomics. The KGMN approach integrates three-layer networks, including knowledge-based metabolic reaction network, knowledge-guided MS/MS similarity network, and global peak correlation network. To demonstrate the principle, we apply KGMN in an in vitro enzymatic reaction system and different biological samples, with ~100-300 putative unknowns annotated in each data set. Among them, >80% unknown metabolites are corroborated with in silico MS/MS tools. Finally, we validate 5 metabolites that are absent in common MS/MS libraries through repository mining and synthesis of chemical standards. Together, the KGMN approach enables efficient unknown annotations, and substantially advances the discovery of recurrent unknown metabolites for common biological samples from model organisms, towards deciphering dark matter in untargeted metabolomics.

Design and Implementation of a New Local Alignment Algorithm for Multilayer Networks

Network alignment (NA) is a popular research field that aims to develop algorithms for comparing networks. Applications of network alignment span many fields, from biology to social network analysis. NA comes in two forms: global network alignment (GNA), which aims to find a global similarity, and LNA, which aims to find local regions of similarity. Recently, there has been an increasing interest in introducing complex network models such as multilayer networks. Multilayer networks are common in many application scenarios, such as modelling of relations among people in a social network or representing the interplay of different molecules in a cell or different cells in the brain. Consequently, the need to introduce algorithms for the comparison of such multilayer networks, i.e., local network alignment, arises. Existing algorithms for LNA do not perform well on multilayer networks since they cannot consider inter-layer edges. Thus, we propose local alignment of multilayer networks (MultiLoAl), a novel algorithm for the local alignment of multilayer networks. We define the local alignment of multilayer networks and propose a heuristic for solving it. We present an extensive assessment indicating the strength of the algorithm. Furthermore, we implemented a synthetic multilayer network generator to build the data for the algorithm's evaluation.

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

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

From Isles of Königsberg to Islets of Langerhans: Examining the Function of the Endocrine Pancreas Through Network Science

Islets of Langerhans are multicellular microorgans located in the pancreas that play a central role in whole-body energy homeostasis. Through secretion of insulin and other hormones they regulate postprandial storage and interprandial usage of energy-rich nutrients. In these clusters of hormone-secreting endocrine cells, intricate cell-cell communication is essential for proper function. Electrical coupling between the insulin-secreting beta cells through gap junctions composed of connexin36 is particularly important, as it provides the required, most important, basis for coordinated responses of the beta cell population. The increasing evidence that gap-junctional communication and its modulation are vital to well-regulated secretion of insulin has stimulated immense interest in how subpopulations of heterogeneous beta cells are functionally arranged throughout the islets and how they mediate intercellular signals. In the last decade, several novel techniques have been proposed to assess cooperation between cells in islets, including the prosperous combination of multicellular imaging and network science. In the present contribution, we review recent advances related to the application of complex network approaches to uncover the functional connectivity patterns among cells within the islets. We first provide an accessible introduction to the basic principles of network theory, enumerating the measures characterizing the intercellular interactions and quantifying the functional integration and segregation of a multicellular system. Then we describe methodological approaches to construct functional beta cell networks, point out possible pitfalls, and specify the functional implications of beta cell network examinations. We continue by highlighting the recent findings obtained through advanced multicellular imaging techniques supported by network-based analyses, giving special emphasis to the current developments in both mouse and human islets, as well as outlining challenges offered by the multilayer network formalism in exploring the collective activity of islet cell populations. Finally, we emphasize that the combination of these imaging techniques and network-based analyses does not only represent an innovative concept that can be used to describe and interpret the physiology of islets, but also provides fertile ground for delineating normal from pathological function and for quantifying the changes in islet communication networks associated with the development of diabetes mellitus.

Networks and Graphs Discovery in Metabolomics Data Analysis and Interpretation

Both targeted and untargeted mass spectrometry-based metabolomics approaches are used to understand the metabolic processes taking place in various organisms, from prokaryotes, plants, fungi to animals and humans. Untargeted approaches allow to detect as many metabolites as possible at once, identify unexpected metabolic changes, and characterize novel metabolites in biological samples. However, the identification of metabolites and the biological interpretation of such large and complex datasets remain challenging. One approach to address these challenges is considering that metabolites are connected through informative relationships. Such relationships can be formalized as networks, where the nodes correspond to the metabolites or features (when there is no or only partial identification), and edges connect nodes if the corresponding metabolites are related. Several networks can be built from a single dataset (or a list of metabolites), where each network represents different relationships, such as statistical (correlated metabolites), biochemical (known or putative substrates and products of reactions), or chemical (structural similarities, ontological relations). Once these networks are built, they can subsequently be mined using algorithms from network (or graph) theory to gain insights into metabolism. For instance, we can connect metabolites based on prior knowledge on enzymatic reactions, then provide suggestions for potential metabolite identifications, or detect clusters of co-regulated metabolites. In this review, we first aim at settling a nomenclature and formalism to avoid confusion when referring to different networks used in the field of metabolomics. Then, we present the state of the art of network-based methods for mass spectrometry-based metabolomics data analysis, as well as future developments expected in this area. We cover the use of networks applications using biochemical reactions, mass spectrometry features, chemical structural similarities, and correlations between metabolites. We also describe the application of knowledge networks such as metabolic reaction networks. Finally, we discuss the possibility of combining different networks to analyze and interpret them simultaneously.

Brain Network Organization Following Post-Stroke Neurorehabilitation

Brain network analysis can offer useful information to guide the rehabilitation of post-stroke patients. We applied functional network connection models based on multiplex-multilayer network analysis (MMN) to explore functional network connectivity changes induced by robot-aided gait training (RAGT) using the Ekso, a wearable exoskeleton, and compared it to conventional overground gait training (COGT) in chronic stroke patients. We extracted the coreness of individual nodes at multiple locations in the brain from EEG recordings obtained before and after gait training in a resting state. We found that patients provided with RAGT achieved a greater motor function recovery than those receiving COGT. This difference in clinical outcome was paralleled by greater changes in connectivity patterns among different brain areas central to motor programming and execution, as well as a recruitment of other areas beyond the sensorimotor cortices and at multiple frequency ranges, contemporarily. The magnitude of these changes correlated with motor function recovery chances. Our data suggest that the use of RAGT as an add-on treatment to COGT may provide post-stroke patients with a greater modification of the functional brain network impairment following a stroke. This might have potential clinical implications if confirmed in large clinical trials.

A semantic framework supporting multilayer networks analysis for rare diseases

Understanding the role played by genetic variations in diseases, exploring genomic variants and discovering disease-associated loci are among the most pressing challenges of genomic medicine. A huge and ever-increasing amount of information is available to researchers to address these challenges. Unfortunately, it is stored in fragmented ontologies and databases, which use heterogeneous formats and poorly integrated schemas. To overcome these limitations, we propose a linked data approach, based on the formalism of multilayer networks, able to integrate and harmonize biomedical information from multiple sources into a single dense network covering different aspects on Neuroendocrine Neoplasms (NENs). The proposed integration schema consists of three interconnected layers representing, respectively, information on the disease, on the affected genes, on the related biological processes and molecular functions. An easy-to-use client-server application was also developed to browse and search for information on the model supporting multilayer network analysis.

Network-Based Approaches Reveal Potential Therapeutic Targets for Host-Directed Antileishmanial Therapy Driving Drug Repurposing

Leishmania parasites are the causal agent of leishmaniasis, an endemic disease in more than 90 countries worldwide. Over the years, traditional approaches focused on the parasite when developing treatments against leishmaniasis. Despite numerous attempts, there is not yet a universal treatment, and those available have allowed for the appearance of resistance. Here, we propose and follow a host-directed approach that aims to overcome the current lack of treatment. Our approach identifies potential therapeutic targets in the host cell and proposes known drug interactions aiming to improve the immune response and to block the host machinery necessary for the survival of the parasite. We started analyzing transcription factor regulatory networks of macrophages infected with Leishmania major. Next, based on the regulatory dynamics of the infection and available gene expression profiles, we selected potential therapeutic target proteins. The function of these proteins was then analyzed following a multilayered network scheme in which we combined information on metabolic pathways with known drugs that have a direct connection with the activity carried out by these proteins. Using our approach, we were able to identify five host protein-coding gene products that are potential therapeutic targets for treating leishmaniasis. Moreover, from the 11 drugs known to interact with the function performed by these proteins, 3 have already been tested against this parasite, verifying in this way our novel methodology. More importantly, the remaining eight drugs previously employed to treat other diseases, remain as promising yet-untested antileishmanial therapies. IMPORTANCE This work opens a new path to fight parasites by targeting host molecular functions by repurposing available and approved drugs. We created a novel approach to identify key proteins involved in any biological process by combining gene regulatory networks and expression profiles. Once proteins have been selected, our approach employs a multilayered network methodology that relates proteins to functions to drugs that alter these functions. By applying our novel approach to macrophages during the Leishmania infection process, we both validated our work and found eight drugs already approved for use in humans that to the best of our knowledge were never employed to treat leishmaniasis, rendering our work as a new tool in the box available to the scientific community fighting parasites.

Network Biology Approaches to Achieve Precision Medicine in Inflammatory Bowel Disease

Inflammatory bowel disease (IBD) is a chronic immune-mediated condition arising due to complex interactions between multiple genetic and environmental factors. Despite recent advances, the pathogenesis of the condition is not fully understood and patients still experience suboptimal clinical outcomes. Over the past few years, investigators are increasingly capturing multi-omics data from patient cohorts to better characterise the disease. However, reaching clinically translatable endpoints from these complex multi-omics datasets is an arduous task. Network biology, a branch of systems biology that utilises mathematical graph theory to represent, integrate and analyse biological data through networks, will be key to addressing this challenge. In this narrative review, we provide an overview of various types of network biology approaches that have been utilised in IBD including protein-protein interaction networks, metabolic networks, gene regulatory networks and gene co-expression networks. We also include examples of multi-layered networks that have combined various network types to gain deeper insights into IBD pathogenesis. Finally, we discuss the need to incorporate other data sources including metabolomic, histopathological, and high-quality clinical meta-data. Together with more robust network data integration and analysis frameworks, such efforts have the potential to realise the key goal of precision medicine in IBD.

Multilayer networks: An untapped tool for understanding bilingual neurocognition

Cross-linguistic similarity is a term so broad and multi-faceted that it is not easily defined. The degree of overlap between languages is known to affect lexical competition during online processing and production, and its relevance for second language acquisition has also been established. Nevertheless, determining what makes two languages similar (or not) increases in complexity when multiple levels of the language hierarchy (e.g., phonology, syntax) are considered. How can we feasibly account for the patterns of convergence and divergence at each level of representation, as well as the interactions between them? The growing field of network science brings new methodologies to bear on this longstanding question. Below, we summarize current network science approaches to modeling language structure and discuss implications for understanding various linguistic processes. Critically, we stress the particular value of multilayer techniques, unique and powerful in their ability to simultaneously accommodate an array of node-to-node (or word-to-word) relationships.

Comparative Analyses of Gene Co-expression Networks: Implementations and Applications in the Study of Evolution

Similarities and differences in the associations of biological entities among species can provide us with a better understanding of evolutionary relationships. Often the evolution of new phenotypes results from changes to interactions in pre-existing biological networks and comparing networks across species can identify evidence of conservation or adaptation. Gene co-expression networks (GCNs), constructed from high-throughput gene expression data, can be used to understand evolution and the rise of new phenotypes. The increasing abundance of gene expression data makes GCNs a valuable tool for the study of evolution in non-model organisms. In this paper, we cover motivations for why comparing these networks across species can be valuable for the study of evolution. We also review techniques for comparing GCNs in the context of evolution, including local and global methods of graph alignment. While some protein-protein interaction (PPI) bioinformatic methods can be used to compare co-expression networks, they often disregard highly relevant properties, including the existence of continuous and negative values for edge weights. Also, the lack of comparative datasets in non-model organisms has hindered the study of evolution using PPI networks. We also discuss limitations and challenges associated with cross-species comparison using GCNs, and provide suggestions for utilizing co-expression network alignments as an indispensable tool for evolutionary studies going forward.

NIDM: network impulsive dynamics on multiplex biological network for disease-gene prediction

The prediction of genes related to diseases is important to the study of the diseases due to high cost and time consumption of biological experiments. Network propagation is a popular strategy for disease-gene prediction. However, existing methods focus on the stable solution of dynamics while ignoring the useful information hidden in the dynamical process, and it is still a challenge to make use of multiple types of physical/functional relationships between proteins/genes to effectively predict disease-related genes. Therefore, we proposed a framework of network impulsive dynamics on multiplex biological network (NIDM) to predict disease-related genes, along with four variants of NIDM models and four kinds of impulsive dynamical signatures (IDSs). NIDM is to identify disease-related genes by mining the dynamical responses of nodes to impulsive signals being exerted at specific nodes. By a series of experimental evaluations in various types of biological networks, we confirmed the advantage of multiplex network and the important roles of functional associations in disease-gene prediction, demonstrated superior performance of NIDM compared with four types of network-based algorithms and then gave the effective recommendations of NIDM models and IDS signatures. To facilitate the prioritization and analysis of (candidate) genes associated to specific diseases, we developed a user-friendly web server, which provides three kinds of filtering patterns for genes, network visualization, enrichment analysis and a wealth of external links (http://bioinformatics.csu.edu.cn/DGP/NID.jsp). NIDM is a protocol for disease-gene prediction integrating different types of biological networks, which may become a very useful computational tool for the study of disease-related genes.