The functional consequences of mutualistic network architecture

The architecture and properties of many complex networks play a significant role in the functioning of the systems they describe. Recently, complex network theory has been applied to ecological entities, like food webs or mutualistic plant-animal interactions. Unfortunately, we still lack an accurate view of the relationship between the architecture and functioning of ecological networks. In this study we explore this link by building individual-based pollination networks from eight Erysimum mediohispanicum (Brassicaceae) populations. In these individual-based networks, each individual plant in a population was considered a node, and was connected by means of undirected links to conspecifics sharing pollinators. The architecture of these unipartite networks was described by means of nestedness, connectivity and transitivity. Network functioning was estimated by quantifying the performance of the population described by each network as the number of per-capita juvenile plants produced per population. We found a consistent relationship between the topology of the networks and their functioning, since variation across populations in the average per-capita production of juvenile plants was positively and significantly related with network nestedness, connectivity and clustering. Subtle changes in the composition of diverse pollinator assemblages can drive major consequences for plant population performance and local persistence through modifications in the structure of the inter-plant pollination networks.     

New moves in motor control

Motor behaviour results from information processing across multiple neural networks acting at all levels from initial selection of the behaviour to its final generation. Understanding how motor behaviour is produced requires identifying the constituent neurons of these networks, their cellular properties, and their pattern of synaptic connectivity. Neural networks have been traditionally studied with neurophysiological and neuroanatomical approaches. These approaches have been highly successful in particularly suitable 'model' preparations, typically ones in which the numbers of neurons in the networks were relatively small, neural network composition was unvarying across individual animals, and the preparations continued to produce fictive motor patterns in?vitro. However, analysing networks without these characteristics, and analysing the complete ensemble of networks that cooperatively generate behaviours, is difficult with these approaches. Recently developed molecular and neurogenetic tools provide additional avenues for analysing motor networks by allowing individual or groups of neurons within networks to be manipulated in novel ways and allowing experiments to be performed not only in?vitro but also in?vivo. We review here some of the new insights into motor network function that these advances have provided and indicate how these advances might bridge gaps in our understanding of motor control. To these ends, we first review motor neural network organisation highlighting cross-phylum principles. We then use prominent examples from the field to show how neurogenetic approaches can complement classical physiological studies, and identify additional areas where these approaches could be advantageously applied.     

Thematic review series: systems biology approaches to metabolic and cardiovascular disorders. Multi-organ whole-genome measurements and reverse engineering to uncover gene networks underlying complex traits

Together with computational analysis and modeling, the development of whole-genome measurement technologies holds the potential to fundamentally change research on complex disorders such as coronary artery disease. With these tools, the stage has been set to reveal the full repertoire of biological components (genes, proteins, and metabolites) in complex diseases and their interplay in modules and networks. Here we review how network identification based on reverse engineering, as applied to whole-genome datasets from simpler organisms, is now being adapted to more complex settings such as datasets from human cell lines and organs in relation to physiological and pathological states. Our focus is on the use of a systems biological approach to identify gene networks in coronary atherosclerosis. We also address how gene networks will probably play a key role in the development of early diagnostics and treatments for complex disorders in the coming era of individualized medicine.     

Automatic network fingerprinting through single-node motifs

Complex networks have been characterised by their specific connectivity patterns (network motifs), but their building blocks can also be identified and described by node-motifs-a combination of local network features. One technique to identify single node-motifs has been presented by Costa et al. (L. D. F. Costa, F. A. Rodrigues, C. C. Hilgetag, and M. Kaiser, Europhys. Lett., 87, 1, 2009). Here, we first suggest improvements to the method including how its parameters can be determined automatically. Such automatic routines make high-throughput studies of many networks feasible. Second, the new routines are validated in different network-series. Third, we provide an example of how the method can be used to analyse network time-series. In conclusion, we provide a robust method for systematically discovering and classifying characteristic nodes of a network. In contrast to classical motif analysis, our approach can identify individual components (here: nodes) that are specific to a network. Such special nodes, as hubs before, might be found to play critical roles in real-world networks.     

Driving and driven architectures of directed small-world human brain functional networks

Recently, increasing attention has been focused on the investigation of the human brain connectome that describes the patterns of structural and functional connectivity networks of the human brain. Many studies of the human connectome have demonstrated that the brain network follows a small-world topology with an intrinsically cohesive modular structure and includes several network hubs in the medial parietal regions. However, most of these studies have only focused on undirected connections between regions in which the directions of information flow are not taken into account. How the brain regions causally influence each other and how the directed network of human brain is topologically organized remain largely unknown. Here, we applied linear multivariate Granger causality analysis (GCA) and graph theoretical approaches to a resting-state functional MRI dataset with a large cohort of young healthy participants (n = 86) to explore connectivity patterns of the population-based whole-brain functional directed network. This directed brain network exhibited prominent small-world properties, which obviously improved previous results of functional MRI studies showing weak small-world properties in the directed brain networks in terms of a kernel-based GCA and individual analysis. This brain network also showed significant modular structures associated with 5 well known subsystems: fronto-parietal, visual, paralimbic/limbic, subcortical and primary systems. Importantly, we identified several driving hubs predominantly located in the components of the attentional network (e.g., the inferior frontal gyrus, supplementary motor area, insula and fusiform gyrus) and several driven hubs predominantly located in the components of the default mode network (e.g., the precuneus, posterior cingulate gyrus, medial prefrontal cortex and inferior parietal lobule). Further split-half analyses indicated that our results were highly reproducible between two independent subgroups. The current study demonstrated the directions of spontaneous information flow and causal influences in the directed brain networks, thus providing new insights into our understanding of human brain functional connectome.     

Multiscale Systems Analysis of Root Growth and Development: Modeling Beyond the Network and Cellular Scales

Over recent decades, we have gained detailed knowledge of many processes involved in root growth and development. However, with this knowledge come increasing complexity and an increasing need for mechanistic modeling to understand how those individual processes interact. One major challenge is in relating genotypes to phenotypes, requiring us to move beyond the network and cellular scales, to use multiscale modeling to predict emergent dynamics at the tissue and organ levels. In this review, we highlight recent developments in multiscale modeling, illustrating how these are generating new mechanistic insights into the regulation of root growth and development. We consider how these models are motivating new biological data analysis and explore directions for future research. This modeling progress will be crucial as we move from a qualitative to an increasingly quantitative understanding of root biology, generating predictive tools that accelerate the development of improved crop varieties.     

Protein interaction networks in medicine and disease

The physical interaction of proteins is subject to intense investigation that has revealed that proteins are assembled into large densely connected networks. In this review, we will examine how signaling pathways can be combined to form higher order protein interaction networks. By using network graph theory, these interaction networks can be further analyzed for global organization, which has revealed unique aspects of the relationships between protein networks and complex biological phenotypes. Moreover, several studies have shown that the structure and dynamics of protein networks are disturbed in complex diseases such as cancer progression. These relationships suggest a novel paradigm for treatment of complex multigenic disease where the protein interaction network is the target of therapy more so than individual molecules within the network.     

Aspects of randomness in neural graph structures

In the past two decades, significant advances have been made in understanding the structural and functional properties of biological networks, via graph-theoretic analysis. In general, most graph-theoretic studies are conducted in the presence of serious uncertainties, such as major undersampling of the experimental data. In the specific case of neural systems, however, a few moderately robust experimental reconstructions have been reported, and these have long served as fundamental prototypes for studying connectivity patterns in the nervous system. In this paper, we provide a comparative analysis of these “historical” graphs, both in their directed (original) and symmetrized (a common preprocessing step) forms, and provide a set of measures that can be consistently applied across graphs (directed or undirected, with or without self-loops). We focus on simple structural characterizations of network connectivity and find that in many measures, the networks studied are captured by simple random graph models. In a few key measures, however, we observe a marked departure from the random graph prediction. Our results suggest that the mechanism of graph formation in the networks studied is not well captured by existing abstract graph models in their first- and second-order connectivity.     

Information processing in social insect networks

Investigating local-scale interactions within a network makes it possible to test hypotheses about the mechanisms of global network connectivity and to ask whether there are general rules underlying network function across systems. Here we use motif analysis to determine whether the interactions within social insect colonies resemble the patterns exhibited by other animal associations or if they exhibit characteristics of biological regulatory systems. Colonies exhibit a predominance of feed-forward interaction motifs, in contrast to the densely interconnected clique patterns that characterize human interaction and animal social networks. The regulatory motif signature supports the hypothesis that social insect colonies are shaped by selection for network patterns that integrate colony functionality at the group rather than individual level, and demonstrates the utility of this approach for analysis of selection effects on complex systems across biological levels of organization.     

Fine-Scale Dissection of Functional Protein Network Organization by Statistical Network Analysis

Revealing organizational principles of biological networks is an important goal of systems biology. In this study, we sought to analyze the dynamic organizational principles within the protein interaction network by studying the characteristics of individual neighborhoods of proteins within the network based on their gene expression as well as protein-protein interaction patterns. By clustering proteins into distinct groups based on their neighborhood gene expression characteristics, we identify several significant trends in the dynamic organization of the protein interaction network. We show that proteins with distinct neighborhood gene expression characteristics are positioned in specific localities in the protein interaction network thereby playing specific roles in the dynamic network connectivity. Remarkably, our analysis reveals a neighborhood characteristic that corresponds to the most centrally located group of proteins within the network. Further, we show that the connectivity pattern displayed by this group is consistent with the notion of “rich club connectivity” in complex networks. Importantly, our findings are largely reproducible in networks constructed using independent and different datasets.     

Proteins encoded in genomic regions associated with immune-mediated disease physically interact and suggest underlying biology

Genome-wide association studies (GWAS) have defined over 150 genomic regions unequivocally containing variation predisposing to immune-mediated disease. Inferring disease biology from these observations, however, hinges on our ability to discover the molecular processes being perturbed by these risk variants. It has previously been observed that different genes harboring causal mutations for the same Mendelian disease often physically interact. We sought to evaluate the degree to which this is true of genes within strongly associated loci in complex disease. Using sets of loci defined in rheumatoid arthritis (RA) and Crohn's disease (CD) GWAS, we build protein-protein interaction (PPI) networks for genes within associated loci and find abundant physical interactions between protein products of associated genes. We apply multiple permutation approaches to show that these networks are more densely connected than chance expectation. To confirm biological relevance, we show that the components of the networks tend to be expressed in similar tissues relevant to the phenotypes in question, suggesting the network indicates common underlying processes perturbed by risk loci. Furthermore, we show that the RA and CD networks have predictive power by demonstrating that proteins in these networks, not encoded in the confirmed list of disease associated loci, are significantly enriched for association to the phenotypes in question in extended GWAS analysis. Finally, we test our method in 3 non-immune traits to assess its applicability to complex traits in general. We find that genes in loci associated to height and lipid levels assemble into significantly connected networks but did not detect excess connectivity among Type 2 Diabetes (T2D) loci beyond chance. Taken together, our results constitute evidence that, for many of the complex diseases studied here, common genetic associations implicate regions encoding proteins that physically interact in a preferential manner, in line with observations in Mendelian disease.     

Raccoon Social Networks and the Potential for Disease Transmission

Raccoons are an important vector of rabies and other pathogens. The degree to which these pathogens can spread through a raccoon population should be closely linked to association rates between individual raccoons. Most studies of raccoon sociality have found patterns consistent with low levels of social connectivity within populations, thus the likelihood of direct pathogen transmission between raccoons is theoretically low. We used proximity detecting collars and social network metrics to calculate the degree of social connectivity in an urban raccoon population for purposes of estimating potential pathogen spread. In contrast to previous assumptions, raccoon social association networks were highly connected, and all individuals were connected to one large social network during 15 out of 18 months of study. However, these metrics may overestimate the potential for a pathogen to spread through a population, as many of the social connections were based on relatively short contact periods. To more closely reflect varying probabilities of pathogen spread, we censored the raccoon social networks based on the total amount of time spent in close proximity between two individuals per month. As this time criteria for censoring the social networks increased from one to thirty minutes, corresponding measures of network connectivity declined. These findings demonstrate that raccoon populations are much more tightly connected than would have been predicted based on previous studies, but also point out that additional research is needed to calculate more precise transmission probabilities by infected individuals, and determine how disease infection changes normal social behaviors.     

A modified Ising model of Barabási–Albert network with gene-type spins

The central question of systems biology is to understand how individual components of a biological system such as genes or proteins cooperate in emerging phenotypes resulting in the evolution of diseases. As living cells are open systems in quasi-steady state type equilibrium in continuous exchange with their environment, computational techniques that have been successfully applied in statistical thermodynamics to describe phase transitions may provide new insights to the emerging behavior of biological systems. Here we systematically evaluate the translation of computational techniques from solid-state physics to network models that closely resemble biological networks and develop specific translational rules to tackle problems unique to living systems. We focus on logic models exhibiting only two states in each network node. Motivated by the apparent asymmetry between biological states where an entity exhibits boolean states i.e. is active or inactive, we present an adaptation of symmetric Ising model towards an asymmetric one fitting to living systems here referred to as the modified Ising model with gene-type spins. We analyze phase transitions by Monte Carlo simulations and propose a mean-field solution of a modified Ising model of a network type that closely resembles a real-world network, the Barabási–Albert model of scale-free networks. We show that asymmetric Ising models show similarities to symmetric Ising models with the external field and undergoes a discontinuous phase transition of the first-order and exhibits hysteresis. The simulation setup presented herein can be directly used for any biological network connectivity dataset and is also applicable for other networks that exhibit similar states of activity. The method proposed here is a general statistical method to deal with non-linear large scale models arising in the context of biological systems and is scalable to any network size.

     

Identification of potential targets in biological signalling systems through network perturbation analysis

The network-based representation and analysis of biological systems contributes to a greater understanding of their structures and functions at different levels of complexity. These techniques can also be used to identify potential novel therapeutic targets based on the characterisation of vulnerable or highly influential network components. There is a need to investigate methods for estimating the impact of molecular perturbations. The prediction of high-impact or critical targets can aid in the identification of novel strategies for controlling the level of activation of specific, therapeutically relevant genes or proteins. Here, we report a new computational strategy for the analysis of the vulnerability of cellular signalling networks based on the quantitative assessment of the impact of large-scale, dynamic perturbations. To show the usefulness of this methodology, two complex signalling networks were analysed: the caspase-3 and the adenosine-regulated calcium signalling systems. This allowed us to estimate and rank the perturbation impact of the components defining these networks. Testable hypotheses about how these targets could modify the dynamic operation of the systems are provided. In the case of the caspase-3 system, the predictions and rankings were in line with results obtained from previous experimental validations of computational predictions generated by a relatively more computationally complex technique. In the case of the adenosine-regulated calcium system, we offer new testable predictions on the potential effect of different targets on the control of calcium flux. Unlike previous methods, the proposed approach provides perturbation-specific scores for each network component. The proposed perturbation assessment methodology may be applied to other systems to gain a deeper understanding of their dynamic operation and to assist the discovery of new therapeutic targets and strategies.     

Graph theoretical analysis of complex networks in the brain

Since the discovery of small-world and scale-free networks the study of complex systems from a network perspective has taken an enormous flight. In recent years many important properties of complex networks have been delineated. In particular, significant progress has been made in understanding the relationship between the structural properties of networks and the nature of dynamics taking place on these networks. For instance, the 'synchronizability' of complex networks of coupled oscillators can be determined by graph spectral analysis. These developments in the theory of complex networks have inspired new applications in the field of neuroscience. Graph analysis has been used in the study of models of neural networks, anatomical connectivity, and functional connectivity based upon fMRI, EEG and MEG. These studies suggest that the human brain can be modelled as a complex network, and may have a small-world structure both at the level of anatomical as well as functional connectivity. This small-world structure is hypothesized to reflect an optimal situation associated with rapid synchronization and information transfer, minimal wiring costs, as well as a balance between local processing and global integration. The topological structure of functional networks is probably restrained by genetic and anatomical factors, but can be modified during tasks. There is also increasing evidence that various types of brain disease such as Alzheimer's disease, schizophrenia, brain tumours and epilepsy may be associated with deviations of the functional network topology from the optimal small-world pattern.     

Failure tolerance of motif structure in biological networks

Complex networks serve as generic models for many biological systems that have been shown to share a number of common structural properties such as power-law degree distribution and small-worldness. Real-world networks are composed of building blocks called motifs that are indeed specific subgraphs of (usually) small number of nodes. Network motifs are important in the functionality of complex networks, and the role of some motifs such as feed-forward loop in many biological networks has been heavily studied. On the other hand, many biological networks have shown some degrees of robustness in terms of their efficiency and connectedness against failures in their components. In this paper we investigated how random and systematic failures in the edges of biological networks influenced their motif structure. We considered two biological networks, namely, protein structure network and human brain functional network. Furthermore, we considered random failures as well as systematic failures based on different strategies for choosing candidate edges for removal. Failure in the edges tipping to high degree nodes had the most destructive role in the motif structure of the networks by decreasing their significance level, while removing edges that were connected to nodes with high values of betweenness centrality had the least effect on the significance profiles. In some cases, the latter caused increase in the significance levels of the motifs.