Structure-Dynamics Relationships in Bursting Neuronal Networks Revealed Using a Prediction Framework

The question of how the structure of a neuronal network affects its functionality has gained a lot of attention in neuroscience. However, the vast majority of the studies on structure-dynamics relationships consider few types of network structures and assess limited numbers of structural measures. In this in silico study, we employ a wide diversity of network topologies and search among many possibilities the aspects of structure that have the greatest effect on the network excitability. The network activity is simulated using two point-neuron models, where the neurons are activated by noisy fluctuation of the membrane potential and their connections are described by chemical synapse models, and statistics on the number and quality of the emergent network bursts are collected for each network type. We apply a prediction framework to the obtained data in order to find out the most relevant aspects of network structure. In this framework, predictors that use different sets of graph-theoretic measures are trained to estimate the activity properties, such as burst count or burst length, of the networks. The performances of these predictors are compared with each other. We show that the best performance in prediction of activity properties for networks with sharp in-degree distribution is obtained when the prediction is based on clustering coefficient. By contrast, for networks with broad in-degree distribution, the maximum eigenvalue of the connectivity graph gives the most accurate prediction. The results shown for small () networks hold with few exceptions when different neuron models, different choices of neuron population and different average degrees are applied. We confirm our conclusions using larger () networks as well. Our findings reveal the relevance of different aspects of network structure from the viewpoint of network excitability, and our integrative method could serve as a general framework for structure-dynamics studies in biosciences.     

Baboon (Papio anubis) social complexity--a network approach

Although many studies have analyzed the causes and consequences of social relationships, few studies have explicitly assessed how measures of social relationships are affected by the choice of behaviors used to quantify them. The use of many behaviors to measure social relationships in primates has long been advocated, but it was analytically difficult to implement this framework into primatological work. However, recent advances in social network analysis (SNA) now allow the comparison of multiple networks created from different behaviors. Here we use our database of baboon social behavior (Papio anubis, Gashaka Gumti National Park, Nigeria) to investigate (i) to what extent social networks created from different behaviors overlap, (ii) to what extent individuals occupy similar social positions in these networks and (iii) how sex affects social network position in this population of baboons. We used data on grooming, aggression, displacement, mounting and presenting, which were collected over a 15-month period. We calculated network parameters separately for each behavior. Networks based on displacement, mounting and presenting were very similar to each other, whereas grooming and aggression networks differed both from each other and from mounting, displacement and presenting networks. Overall, individual network positions were strongly affected by sex. Individuals central in one network tended to be central in most other networks as well, whereas other measures such as clustering coefficient were found to vary depending on the behavior analyzed. Thus, our results suggest that a baboon's social environment is best described by a multiplex network based on affiliative, aggressive and sexual behavior. Modern SNA provides a number of useful tools that will help us to better understand animals' social environment. We also discuss potential caveats related to their use.     

Functional coherence in domain interaction networks

MOTIVATION: Extracting functional information from protein-protein interactions (PPI) poses significant challenges arising from the noisy, incomplete, generic and static nature of data obtained from high-throughput screening. Typical proteins are composed of multiple domains, often regarded as their primary functional and structural units. Motivated by these considerations, domain-domain interactions (DDI) for network-based analyses have received significant recent attention. This article performs a formal comparative investigation of the relationship between functional coherence and topological proximity in PPI and DDI networks. Our investigation provides the necessary basis for continued and focused investigation of DDIs as abstractions for functional characterization and modularization of networks. RESULTS: We investigate the problem of assessing the functional coherence of two biomolecules (or segments thereof) in a formal framework. We establish essential attributes of admissible measures of functional coherence, and demonstrate that existing, well-accepted measures are ill-suited to comparative analyses involving different entities (i.e. domains versus proteins). We propose a statistically motivated functional similarity measure that takes into account functional specificity as well as the distribution of functional attributes across entity groups to assess functional similarity in a statistically meaningful and biologically interpretable manner. Results on diverse data, including high-throughput and computationally predicted PPIs, as well as structural and computationally inferred DDIs for different organisms show that: (i) the relationship between functional similarity and network proximity is captured in a much more (biologically) intuitive manner by our measure, compared to existing measures and (ii) network proximity and functional similarity are significantly more correlated in DDI networks than in PPI networks, and that structurally determined DDIs provide better functional relevance as compared to computationally inferred DDIs.     

Geometric robustness theory and biological networks

We provide a geometric framework for investigating the robustness of information flows over biological networks. We use information measures to quantify the impact of knockout perturbations on simple networks. Robustness has two components, a measure of the causal contribution of a node or nodes, and a measure of the change or exclusion dependence, of the network following node removal. Causality is measured as statistical contribution of a node to network function, wheras exclusion dependence measures a distance between unperturbed network and reconfigured network function. We explore the role that redundancy plays in increasing robustness, and how redundacy can be exploited through error-correcting codes implemented by networks. We provide examples of the robustness measure when applied to familiar boolean functions such as the AND, OR and XOR functions. We discuss the relationship between robustness measures and related measures of complexity and how robustness always implies a minimal level of complexity.     

Structural distance and evolutionary relationship of networks

Exploring common features and universal qualities shared by a particular class of networks in biological and other domains is one of the important aspects of evolutionary study. In an evolving system, evolutionary mechanism can cause functional changes that forces the system to adapt to new configurations of interaction pattern between the components of that system (e.g. gene duplication and mutation play a vital role for changing the connectivity structure in many biological networks. The evolutionary relation between two systems can be retraced by their structural differences). The eigenvalues of the normalized graph Laplacian not only capture the global properties of a network, but also local structures that are produced by graph evolutions (like motif duplication or joining). The spectrum of this operator carries many qualitative aspects of a graph. Given two networks of different sizes, we propose a method to quantify the topological distance between them based on the contrasting spectrum of normalized graph Laplacian. We find that network architectures are more similar within the same class compared to between classes. We also show that the evolutionary relationships can be retraced by the structural differences using our method. We analyze 43 metabolic networks from different species and mark the prominent separation of three groups: Bacteria, Archaea and Eukarya. This phenomenon is well captured in our findings that support the other cladistic results based on gene content and ribosomal RNA sequences. Our measure to quantify the structural distance between two networks is useful to elucidate evolutionary relationships.     

Discriminating Different Classes of Biological Networks by Analyzing the Graphs Spectra Distribution

The brain''s structural and functional systems, protein-protein interaction, and gene networks are examples of biological systems that share some features of complex networks, such as highly connected nodes, modularity, and small-world topology. Recent studies indicate that some pathologies present topological network alterations relative to norms seen in the general population. Therefore, methods to discriminate the processes that generate the different classes of networks (e.g., normal and disease) might be crucial for the diagnosis, prognosis, and treatment of the disease. It is known that several topological properties of a network (graph) can be described by the distribution of the spectrum of its adjacency matrix. Moreover, large networks generated by the same random process have the same spectrum distribution, allowing us to use it as a “fingerprint”. Based on this relationship, we introduce and propose the entropy of a graph spectrum to measure the “uncertainty” of a random graph and the Kullback-Leibler and Jensen-Shannon divergences between graph spectra to compare networks. We also introduce general methods for model selection and network model parameter estimation, as well as a statistical procedure to test the nullity of divergence between two classes of complex networks. Finally, we demonstrate the usefulness of the proposed methods by applying them to (1) protein-protein interaction networks of different species and (2) on networks derived from children diagnosed with Attention Deficit Hyperactivity Disorder (ADHD) and typically developing children. We conclude that scale-free networks best describe all the protein-protein interactions. Also, we show that our proposed measures succeeded in the identification of topological changes in the network while other commonly used measures (number of edges, clustering coefficient, average path length) failed.     

Algebraic connectivity may explain the evolution of gene regulatory networks

Gene expression is a result of the interplay between the structure, type, kinetics, and specificity of gene regulatory interactions, whose diversity gives rise to the variety of life forms. As the dynamic behavior of gene regulatory networks depends on their structure, here we attempt to determine structural reasons which, despite the similarities in global network properties, may explain the large differences in organismal complexity. We demonstrate that the algebraic connectivity, the smallest non-trivial eigenvalue of the Laplacian, of the directed gene regulatory networks decreases with the increase of organismal complexity, and may therefore explain the difference between the variety of analyzed regulatory networks. In addition, our results point out that, for the species considered in this study, evolution favours decreasing concentration of strategically positioned feed forward loops, so that the network as a whole can increase the specificity towards changing environments. Moreover, contrary to the existing results, we show that the average degree, the length of the longest cascade, and the average cascade length of gene regulatory networks cannot recover the evolutionary relationships between organisms. Whereas the dynamical properties of special subnetworks are relatively well understood, there is still limited knowledge about the evolutionary reasons for the already identified design principles pertaining to these special subnetworks, underlying the global quantitative features of gene regulatory networks of different organisms. The behavior of the algebraic connectivity, which we show valid on gene regulatory networks extracted from curated databases, can serve as an additional evolutionary principle of organism-specific regulatory networks.     

Beyond topology: coevolution of structure and flux in metabolic networks

Interactions between the structure of a metabolic network and its functional properties underlie its evolutionary diversification, but the mechanism by which such interactions arise remains elusive. Particularly unclear is whether metabolic fluxes that determine the concentrations of compounds produced by a metabolic network, are causally linked to a network's structure or emerge independently of it. A direct empirical study of populations where both structural and functional properties vary among individuals’ metabolic networks is required to establish whether changes in structure affect the distribution of metabolic flux. In a population of house finches (Haemorhous mexicanus), we reconstructed full carotenoid metabolic networks for 442 individuals and uncovered 11 structural variants of this network with different compounds and reactions. We examined the consequences of this structural diversity for the concentrations of plumage‐bound carotenoids produced by flux in these networks. We found that concentrations of metabolically derived, but not dietary carotenoids, depended on network structure. Flux was partitioned similarly among compounds in individuals of the same network structure: within each network, compound concentrations were closely correlated. The highest among‐individual variation in flux occurred in networks with the strongest among‐compound correlations, suggesting that changes in the magnitude, but not the distribution of flux, underlie individual differences in compound concentrations on a static network structure. These findings indicate that the distribution of flux in carotenoid metabolism closely follows network structure. Thus, evolutionary diversification and local adaptations in carotenoid metabolism may depend more on the gain or loss of enzymatic reactions than on changes in flux within a network structure.     

Quantifying and analyzing the network basis of genetic complexity

Genotype-to-phenotype maps exhibit complexity. This genetic complexity is mentioned frequently in the literature, but a consistent and quantitative definition is lacking. Here, we derive such a definition and investigate its consequences for model genetic systems. The definition equates genetic complexity with a surplus of genotypic diversity over phenotypic diversity. Applying this definition to ensembles of Boolean network models, we found that the in-degree distribution and the number of periodic attractors produced determine the relative complexity of different topology classes. We found evidence that networks that are difficult to control, or that exhibit a hierarchical structure, are genetically complex. We analyzed the complexity of the cell cycle network of Sacchoromyces cerevisiae and pinpointed genes and interactions that are most important for its high genetic complexity. The rigorous definition of genetic complexity is a tool for unraveling the structure and properties of genotype-to-phenotype maps by enabling the quantitative comparison of the relative complexities of different genetic systems. The definition also allows the identification of specific network elements and subnetworks that have the greatest effects on genetic complexity. Moreover, it suggests ways to engineer biological systems with desired genetic properties.     

Large-scale properties of clustered networks: implications for disease dynamics

We consider previously proposed procedures for generating clustered networks and investigate how these procedures lead to differences in network properties other than clustering. We interpret our findings in terms of the effect of the network structure on the disease outbreak threshold and disease dynamics. To generate null-model networks for comparison, we implement an assortativity-conserving rewiring algorithm that alters the level of clustering while causing minimal impact on other properties. We show that many theoretical network models used to generate networks with a particular property often lead to significant changes in network properties other than that of interest. For high levels of clustering, different procedures lead to networks that differ in degree heterogeneity and assortativity, and in broader scale measures such as ?(0) and the distribution of shortest path lengths. Hence, care must be taken when investigating the implications of network properties for disease transmission or other dynamic process that the network supports.     

Structural Properties and Complexity of a New Network Class: Collatz Step Graphs

In this paper, we introduce a biologically inspired model to generate complex networks. In contrast to many other construction procedures for growing networks introduced so far, our method generates networks from one-dimensional symbol sequences that are related to the so called Collatz problem from number theory. The major purpose of the present paper is, first, to derive a symbol sequence from the Collatz problem, we call the step sequence, and investigate its structural properties. Second, we introduce a construction procedure for growing networks that is based on these step sequences. Third, we investigate the structural properties of this new network class including their finite scaling and asymptotic behavior of their complexity, average shortest path lengths and clustering coefficients. Interestingly, in contrast to many other network models including the small-world network from Watts & Strogatz, we find that CS graphs become ‘smaller’ with an increasing size.     

A decade of seasonal dynamics and co-occurrences within freshwater bacterioplankton communities from eutrophic Lake Mendota,WI, USA

With an unprecedented decade-long time series from a temperate eutrophic lake, we analyzed bacterial and environmental co-occurrence networks to gain insight into seasonal dynamics at the community level. We found that (1) bacterial co-occurrence networks were non-random, (2) season explained the network complexity and (3) co-occurrence network complexity was negatively correlated with the underlying community diversity across different seasons. Network complexity was not related to the variance of associated environmental factors. Temperature and productivity may drive changes in diversity across seasons in temperate aquatic systems, much as they control diversity across latitude. While the implications of bacterioplankton network structure on ecosystem function are still largely unknown, network analysis, in conjunction with traditional multivariate techniques, continues to increase our understanding of bacterioplankton temporal dynamics.     

Structure–stability relationships in networks combining mutualistic and antagonistic interactions

The relationship between the structure of ecological networks and community stability has been studied for decades. Recent developments highlighted that this relationship depended on whether interactions were antagonistic or mutualistic. Different structures promoting stability in different types of ecological networks, i.e. mutualistic or antagonistic, have been pointed out. However, these findings come from studies considering mutualistic and antagonistic interactions separately whereas we know that species are part of both types of networks simultaneously. Understanding the relationship between network structure and community stability, when mutualistic and antagonistic interactions are merged in a single network, thus appears as the next challenge to improve our understanding of the dynamics of natural communities. Using a theoretical approach, we test whether the structural characteristics known to promote stability in networks made of a single interaction type still hold for network merging mutualistic and antagonistic interactions. We show that the effects of diversity and connectance remain unchanged. But the effects of nestedness and modularity are strongly weakened in networks combining mutualistic and antagonistic interactions. By challenging the stabilizing mechanisms proposed for networks with a single interaction type, our study calls for new measures of structure for networks that integrate the diversity of interaction.     

Patterns of genetic variation in populations of infectious agents

Background  

The analysis of genetic variation in populations of infectious agents may help us understand their epidemiology and evolution. Here we study a model for assessing the levels and patterns of genetic diversity in populations of infectious agents. The population is structured into many small subpopulations, which correspond to their hosts, that are connected according to a specific type of contact network. We considered different types of networks, including fully connected networks and scale free networks, which have been considered as a model that captures some properties of real contact networks. Infectious agents transmit between hosts, through migration, where they grow and mutate until elimination by the host immune system.     

Coemergence of regularity and complexity during neural network development

With the growing recognition that rhythmic and oscillatory patterns are widespread in the brain and play important roles in all aspects of the function of our nervous system, there has been a resurgence of interest in neuronal synchronized bursting activity. Here, we were interested in understanding the development of synchronized bursts as information-bearing neuronal activity patterns. For that, we have monitored the morphological organization and spontaneous activity of neuronal networks cultured on multielectrode-arrays during their self-executed evolvement from a mixture of dissociated cells into an active network. Complex collective network electrical activity evolved from sporadic firing patterns of the single neurons. On the system (network) level, the activity was marked by bursting events with interneuronal synchronization and nonarbitrary temporal ordering. We quantified these individual-to-collective activity transitions using newly-developed system level quantitative measures of time series regularity and complexity. We found that individual neuronal activity before synchronization was characterized by high regularity and low complexity. During neuronal wiring, there was a transient period of reorganization marked by low regularity, which then leads to coemergence of elevated regularity and functional (nonstochastic) complexity. We further investigated the morphology-activity interplay by modeling artificial neuronal networks with different topological organizations and connectivity schemes. The simulations support our experimental results by showing increased levels of complexity of neuronal activity patterns when neurons are wired up and organized in clusters (similar to mature real networks), as well as network-level activity regulation once collective activity forms.     

From social network (centralized vs. decentralized) to collective decision-making (unshared vs. shared consensus)

Relationships we have with our friends, family, or colleagues influence our personal decisions, as well as decisions we make together with others. As in human beings, despotism and egalitarian societies seem to also exist in animals. While studies have shown that social networks constrain many phenomena from amoebae to primates, we still do not know how consensus emerges from the properties of social networks in many biological systems. We created artificial social networks that represent the continuum from centralized to decentralized organization and used an agent-based model to make predictions about the patterns of consensus and collective movements we observed according to the social network. These theoretical results showed that different social networks and especially contrasted ones--star network vs. equal network--led to totally different patterns. Our model showed that, by moving from a centralized network to a decentralized one, the central individual seemed to lose its leadership in the collective movement's decisions. We, therefore, showed a link between the type of social network and the resulting consensus. By comparing our theoretical data with data on five groups of primates, we confirmed that this relationship between social network and consensus also appears to exist in animal societies.     

Modeling the Impact of Lesions in the Human Brain

Lesions of anatomical brain networks result in functional disturbances of brain systems and behavior which depend sensitively, often unpredictably, on the lesion site. The availability of whole-brain maps of structural connections within the human cerebrum and our increased understanding of the physiology and large-scale dynamics of cortical networks allow us to investigate the functional consequences of focal brain lesions in a computational model. We simulate the dynamic effects of lesions placed in different regions of the cerebral cortex by recording changes in the pattern of endogenous (“resting-state”) neural activity. We find that lesions produce specific patterns of altered functional connectivity among distant regions of cortex, often affecting both cortical hemispheres. The magnitude of these dynamic effects depends on the lesion location and is partly predicted by structural network properties of the lesion site. In the model, lesions along the cortical midline and in the vicinity of the temporo-parietal junction result in large and widely distributed changes in functional connectivity, while lesions of primary sensory or motor regions remain more localized. The model suggests that dynamic lesion effects can be predicted on the basis of specific network measures of structural brain networks and that these effects may be related to known behavioral and cognitive consequences of brain lesions.     

Brain rhythms reveal a hierarchical network organization

Recordings of ongoing neural activity with EEG and MEG exhibit oscillations of specific frequencies over a non-oscillatory background. The oscillations appear in the power spectrum as a collection of frequency bands that are evenly spaced on a logarithmic scale, thereby preventing mutual entrainment and cross-talk. Over the last few years, experimental, computational and theoretical studies have made substantial progress on our understanding of the biophysical mechanisms underlying the generation of network oscillations and their interactions, with emphasis on the role of neuronal synchronization. In this paper we ask a very different question. Rather than investigating how brain rhythms emerge, or whether they are necessary for neural function, we focus on what they tell us about functional brain connectivity. We hypothesized that if we were able to construct abstract networks, or "virtual brains", whose dynamics were similar to EEG/MEG recordings, those networks would share structural features among themselves, and also with real brains. Applying mathematical techniques for inverse problems, we have reverse-engineered network architectures that generate characteristic dynamics of actual brains, including spindles and sharp waves, which appear in the power spectrum as frequency bands superimposed on a non-oscillatory background dominated by low frequencies. We show that all reconstructed networks display similar topological features (e.g. structural motifs) and dynamics. We have also reverse-engineered putative diseased brains (epileptic and schizophrenic), in which the oscillatory activity is altered in different ways, as reported in clinical studies. These reconstructed networks show consistent alterations of functional connectivity and dynamics. In particular, we show that the complexity of the network, quantified as proposed by Tononi, Sporns and Edelman, is a good indicator of brain fitness, since virtual brains modeling diseased states display lower complexity than virtual brains modeling normal neural function. We finally discuss the implications of our results for the neurobiology of health and disease.     

Path lengths in protein-protein interaction networks and biological complexity

We investigated the biological significance of path lengths in 12 protein-protein interaction (PPI) networks. We put forward three predictions, based on the idea that biological complexity influences path lengths. First, at the network level, path lengths are generally longer in PPIs than in random networks. Second, this pattern is more pronounced in more complex organisms. Third, within a PPI network, path lengths of individual proteins are biologically significant. We found that in 11 of the 12 species, average path lengths in PPI networks are significantly longer than those in randomly rewired networks. The PPI network of the malaria parasite Plasmodium falciparum, however, does not exhibit deviation from rewired networks. Furthermore, eukaryotic PPIs exhibit significantly greater deviation from randomly rewired networks than prokaryotic PPIs. Thus our study highlights the potentially meaningful variation in path lengths of PPI networks. Moreover, node eccentricity, defined as the longest path from a protein to others, is significantly correlated with the levels of gene expression and dispensability in the yeast PPI network. We conclude that biological complexity influences both global and local properties of path lengths in PPI networks. Investigating variation of path lengths may provide new tools to analyze the evolution of functional modules in biological systems.