Resting-State Temporal Synchronization Networks Emerge from Connectivity Topology and Heterogeneity

Spatial patterns of coherent activity across different brain areas have been identified during the resting-state fluctuations of the brain. However, recent studies indicate that resting-state activity is not stationary, but shows complex temporal dynamics. We were interested in the spatiotemporal dynamics of the phase interactions among resting-state fMRI BOLD signals from human subjects. We found that the global phase synchrony of the BOLD signals evolves on a characteristic ultra-slow (<0.01Hz) time scale, and that its temporal variations reflect the transient formation and dissolution of multiple communities of synchronized brain regions. Synchronized communities reoccurred intermittently in time and across scanning sessions. We found that the synchronization communities relate to previously defined functional networks known to be engaged in sensory-motor or cognitive function, called resting-state networks (RSNs), including the default mode network, the somato-motor network, the visual network, the auditory network, the cognitive control networks, the self-referential network, and combinations of these and other RSNs. We studied the mechanism originating the observed spatiotemporal synchronization dynamics by using a network model of phase oscillators connected through the brain’s anatomical connectivity estimated using diffusion imaging human data. The model consistently approximates the temporal and spatial synchronization patterns of the empirical data, and reveals that multiple clusters that transiently synchronize and desynchronize emerge from the complex topology of anatomical connections, provided that oscillators are heterogeneous.     

Relative stability of core groups in pollination networks in a biodiversity hotspot over four years

Plants and their pollinators form pollination networks integral to the evolution and persistence of species in communities. Previous studies suggest that pollination network structure remains nested while network composition is highly dynamic. However, little is known about temporal variation in the structure and function of plant-pollinator networks, especially in species-rich communities where the strength of pollinator competition is predicted to be high. Here we quantify temporal variation of pollination networks over four consecutive years in an alpine meadow in the Hengduan Mountains biodiversity hotspot in China. We found that ranked positions and idiosyncratic temperatures of both plants and pollinators were more conservative between consecutive years than in non-consecutive years. Although network compositions exhibited high turnover, generalized core groups--decomposed by a k-core algorithm--were much more stable than peripheral groups. Given the high rate of turnover observed, we suggest that identical plants and pollinators that persist for at least two successive years sustain pollination services at the community level. Our data do not support theoretical predictions of a high proportion of specialized links within species-rich communities. Plants were relatively specialized, exhibiting less variability in pollinator composition at pollinator functional group level than at the species level. Both specialized and generalized plants experienced narrow variation in functional pollinator groups. The dynamic nature of pollination networks in the alpine meadow demonstrates the potential for networks to mitigate the effects of fluctuations in species composition in a high biodiversity area.     

Improving molecular cancer class discovery through sparse non-negative matrix factorization

MOTIVATION: Identifying different cancer classes or subclasses with similar morphological appearances presents a challenging problem and has important implication in cancer diagnosis and treatment. Clustering based on gene-expression data has been shown to be a powerful method in cancer class discovery. Non-negative matrix factorization is one such method and was shown to be advantageous over other clustering techniques, such as hierarchical clustering or self-organizing maps. In this paper, we investigate the benefit of explicitly enforcing sparseness in the factorization process. RESULTS: We report an improved unsupervised method for cancer classification by the use of gene-expression profile via sparse non-negative matrix factorization. We demonstrate the improvement by direct comparison with classic non-negative matrix factorization on the three well-studied datasets. In addition, we illustrate how to identify a small subset of co-expressed genes that may be directly involved in cancer.     

Linking dynamics of the inhibitory network to the input structure

Networks of inhibitory interneurons are found in many distinct classes of biological systems. Inhibitory interneurons govern the dynamics of principal cells and are likely to be critically involved in the coding of information. In this theoretical study, we describe the dynamics of a generic inhibitory network in terms of low-dimensional, simplified rate models. We study the relationship between the structure of external input applied to the network and the patterns of activity arising in response to that stimulation. We found that even a minimal inhibitory network can generate a great diversity of spatio-temporal patterning including complex bursting regimes with non-trivial ratios of burst firing. Despite the complexity of these dynamics, the network’s response patterns can be predicted from the rankings of the magnitudes of external inputs to the inhibitory neurons. This type of invariant dynamics is robust to noise and stable in densely connected networks with strong inhibitory coupling. Our study predicts that the response dynamics generated by an inhibitory network may provide critical insights about the temporal structure of the sensory input it receives.     

Dynamics and Control of Diseases in Networks with Community Structure

The dynamics of infectious diseases spread via direct person-to-person transmission (such as influenza, smallpox, HIV/AIDS, etc.) depends on the underlying host contact network. Human contact networks exhibit strong community structure. Understanding how such community structure affects epidemics may provide insights for preventing the spread of disease between communities by changing the structure of the contact network through pharmaceutical or non-pharmaceutical interventions. We use empirical and simulated networks to investigate the spread of disease in networks with community structure. We find that community structure has a major impact on disease dynamics, and we show that in networks with strong community structure, immunization interventions targeted at individuals bridging communities are more effective than those simply targeting highly connected individuals. Because the structure of relevant contact networks is generally not known, and vaccine supply is often limited, there is great need for efficient vaccination algorithms that do not require full knowledge of the network. We developed an algorithm that acts only on locally available network information and is able to quickly identify targets for successful immunization intervention. The algorithm generally outperforms existing algorithms when vaccine supply is limited, particularly in networks with strong community structure. Understanding the spread of infectious diseases and designing optimal control strategies is a major goal of public health. Social networks show marked patterns of community structure, and our results, based on empirical and simulated data, demonstrate that community structure strongly affects disease dynamics. These results have implications for the design of control strategies.     

Time Development in the Early History of Social Networks: Link Stabilization,Group Dynamics,and Segregation

Studies of the time development of empirical networks usually investigate late stages where lasting connections have already stabilized. Empirical data on early network history are rare but needed for a better understanding of how social network topology develops in real life. Studying students who are beginning their studies at a university with no or few prior connections to each other offers a unique opportunity to investigate the formation and early development of link patterns and community structure in social networks. During a nine week introductory physics course, first year physics students were asked to identify those with whom they communicated about problem solving in physics during the preceding week. We use these students'' self reports to produce time dependent student interaction networks. We investigate these networks to elucidate possible effects of different student attributes in early network formation. Changes in the weekly number of links show that while roughly half of all links change from week to week, students also reestablish a growing number of links as they progress through their first weeks of study. Using the Infomap community detection algorithm, we show that the networks exhibit community structure, and we use non-network student attributes, such as gender and end-of-course grade to characterize communities during their formation. Specifically, we develop a segregation measure and show that students structure themselves according to gender and pre-organized sections (in which students engage in problem solving and laboratory work), but not according to end-of-coure grade. Alluvial diagrams of consecutive weeks'' communities show that while student movement between groups are erratic in the beginnning of their studies, they stabilize somewhat towards the end of the course. Taken together, the analyses imply that student interaction networks stabilize quickly and that students establish collaborations based on who is immediately available to them and on observable personal characteristics.     

Sensitive pattern discovery with 'fuzzy' alignments of distantly related proteins

MOTIVATION: Evolutionary comparison leads to efficient functional characterisation of hypothetical proteins. Here, our goal is to map specific sequence patterns to putative functional classes. The evolutionary signal stands out most clearly in a maximally diverse set of homologues. This diversity, however, leads to a number of technical difficulties. The targeted patterns-as gleaned from structure comparisons-are too sparse for statistically significant signals of sequence similarity and accurate multiple sequence alignment. RESULTS: We address this problem by a fuzzy alignment model, which probabilistically assigns residues to structurally equivalent positions (attributes) of the proteins. We then apply multivariate analysis to the 'attributes x proteins' matrix. The dimensionality of the space is reduced using non-negative matrix factorization. The method is general, fully automatic and works without assumptions about pattern density, minimum support, explicit multiple alignments, phylogenetic trees, etc. We demonstrate the discovery of biologically meaningful patterns in an extremely diverse superfamily related to urease.     

Encoding time in neural dynamic regimes with distinct computational tradeoffs

Converging evidence suggests the brain encodes time in dynamic patterns of neural activity, including neural sequences, ramping activity, and complex dynamics. Most temporal tasks, however, require more than just encoding time, and can have distinct computational requirements including the need to exhibit temporal scaling, generalize to novel contexts, or robustness to noise. It is not known how neural circuits can encode time and satisfy distinct computational requirements, nor is it known whether similar patterns of neural activity at the population level can exhibit dramatically different computational or generalization properties. To begin to answer these questions, we trained RNNs on two timing tasks based on behavioral studies. The tasks had different input structures but required producing identically timed output patterns. Using a novel framework we quantified whether RNNs encoded two intervals using either of three different timing strategies: scaling, absolute, or stimulus-specific dynamics. We found that similar neural dynamic patterns at the level of single intervals, could exhibit fundamentally different properties, including, generalization, the connectivity structure of the trained networks, and the contribution of excitatory and inhibitory neurons. Critically, depending on the task structure RNNs were better suited for generalization or robustness to noise. Further analysis revealed different connection patterns underlying the different regimes. Our results predict that apparently similar neural dynamic patterns at the population level (e.g., neural sequences) can exhibit fundamentally different computational properties in regards to their ability to generalize to novel stimuli and their robustness to noise—and that these differences are associated with differences in network connectivity and distinct contributions of excitatory and inhibitory neurons. We also predict that the task structure used in different experimental studies accounts for some of the experimentally observed variability in how networks encode time.     

Anthropogenic disturbance and rainfall variation threaten the stability of plant–ant interactions in the Brazilian Caatinga

Climate change is projected to exacerbate the effects of anthropogenic disturbance, with negative impacts on ecosystem stability and functioning. We evaluate the additive and combined effects of chronic anthropogenic disturbance (CAD) and rainfall variation on the temporal stability of mutualistic EFN‐bearing plant–ant networks in a Caatinga dry forest. We evaluated whether changes in the stability of these interactions are driven by changes in the stability of the communities of partners involved and/or in ant behavior. We sampled EFN‐bearing plant–ant networks in sixteen 20 × 20 m plots distributed across CAD and rainfall gradients. The stability of EFN‐bearing plant and attendant–ant communities were measured as the inverse of temporal differences in their community structure and composition. We also computed the stability of EFN‐bearing plant–ant networks by measuring the inverse of temporal differences in network specialization metrics. We found that, in general, the structure and composition of plant and ant interacting communities were similarly stable along both environmental gradients. Only CAD and its interaction with rainfall affected the temporal stability of EFN‐bearing plant diversity, which declined as CAD increased, with a more pronounced relationship in wetter areas. However, variation in levels of CAD and, to a lesser extent, rainfall greatly modulated the stability of EFN‐bearing plant–ant network specialization. CAD reduced the stability of network generality (specialization at the ant level), an effect that was much stronger in wetter areas. Meanwhile, the stability in network vulnerability (specialization at the plant level) decreased with the increase of CAD and the decrease of rainfall levels. Finally, there was a trend of decreasing stability in specialization of the overall network with increasing CAD. Our results suggest that changes in the structure of interaction networks are mainly driven by a switch in ant behavior rather than by changes in the structure and composition of plant and ant communities between years.     

Link Community Detection Using Generative Model and Nonnegative Matrix Factorization

Discovery of communities in complex networks is a fundamental data analysis problem with applications in various domains. While most of the existing approaches have focused on discovering communities of nodes, recent studies have shown the advantages and uses of link community discovery in networks. Generative models provide a promising class of techniques for the identification of modular structures in networks, but most generative models mainly focus on the detection of node communities rather than link communities. In this work, we propose a generative model, which is based on the importance of each node when forming links in each community, to describe the structure of link communities. We proceed to fit the model parameters by taking it as an optimization problem, and solve it using nonnegative matrix factorization. Thereafter, in order to automatically determine the number of communities, we extend the above method by introducing a strategy of iterative bipartition. This extended method not only finds the number of communities all by itself, but also obtains high efficiency, and thus it is more suitable to deal with large and unexplored real networks. We test this approach on both synthetic benchmarks and real-world networks including an application on a large biological network, and compare it with two highly related methods. Results demonstrate the superior performance of our approach over competing methods for the detection of link communities.     

Effects of biotic and abiotic factors on the temporal dynamic of bat-fruit interactions

Mutualistic interactions between animals and plants vary over time and space based on the abundance of fruits or animals and seasonality. Little is known about this temporal dynamic and the influence of biotic and abiotic factors on the structure of interaction networks. We evaluated changes in the structure of network interactions between bats and fruits in relation to variations in rainfall. Our results suggest that fruit abundance is the main variable responsible for temporal changes in network attributes, such as network size, connectance, and number of interactions. In the same way, temperature positively affected the abundance of fruits and bats. An increase in temperature and alterations in rainfall patterns, due to human induced climate change, can cause changes in phenological patterns and fruit production, with negative consequences to biodiversity maintenance, ecological interactions, and ecosystem functioning.     

Seeing through the static: the temporal dimension of plant–animal mutualistic interactions

Most studies of plant–animal mutualistic networks have come from a temporally static perspective. This approach has revealed general patterns in network structure, but limits our ability to understand the ecological and evolutionary processes that shape these networks and to predict the consequences of natural and human‐driven disturbance on species interactions. We review the growing literature on temporal dynamics of plant–animal mutualistic networks including pollination, seed dispersal and ant defence mutualisms. We then discuss potential mechanisms underlying such variation in interactions, ranging from behavioural and physiological processes at the finest temporal scales to ecological and evolutionary processes at the broadest. We find that at the finest temporal scales (days, weeks, months) mutualistic interactions are highly dynamic, with considerable variation in network structure. At intermediate scales (years, decades), networks still exhibit high levels of temporal variation, but such variation appears to influence network properties only weakly. At the broadest temporal scales (many decades, centuries and beyond), continued shifts in interactions appear to reshape network structure, leading to dramatic community changes, including loss of species and function. Our review highlights the importance of considering the temporal dimension for understanding the ecology and evolution of complex webs of mutualistic interactions.     

Matrix models for quantifying competitive intransitivity from species abundance data

In a network of competing species, a competitive intransitivity occurs when the ranking of competitive abilities does not follow a linear hierarchy (A > B > C but C > A). A variety of mathematical models suggests that intransitive networks can prevent or slow down competitive exclusion and maintain biodiversity by enhancing species coexistence. However, it has been difficult to assess empirically the relative importance of intransitive competition because a large number of pairwise species competition experiments are needed to construct a competition matrix that is used to parameterize existing models. Here we introduce a statistical framework for evaluating the contribution of intransitivity to community structure using species abundance matrices that are commonly generated from replicated sampling of species assemblages. We provide metrics and analytical methods for using abundance matrices to estimate species competition and patch transition matrices by using reverse‐engineering and a colonization–competition model. These matrices provide complementary metrics to estimate the degree of intransitivity in the competition network of the sampled communities. Benchmark tests reveal that the proposed methods could successfully detect intransitive competition networks, even in the absence of direct measures of pairwise competitive strength. To illustrate the approach, we analyzed patterns of abundance and biomass of five species of necrophagous Diptera and eight species of their hymenopteran parasitoids that co‐occur in beech forests in Germany. We found evidence for a strong competitive hierarchy within communities of flies and parasitoids. However, for parasitoids, there was a tendency towards increasing intransitivity in higher weight classes, which represented larger resource patches. These tests provide novel methods for empirically estimating the degree of intransitivity in competitive networks from observational datasets. They can be applied to experimental measures of pairwise species interactions, as well as to spatio‐temporal samples of assemblages in homogenous environments or environmental gradients.     

Exploring the interaction patterns among taxa and environments from marine metagenomic data

The sequencing revolution driven by high-throughput technologies has generated a huge amount of marine microbial sequences which hide the interaction patterns among microbial species and environment factors. Exploring these patterns is helpful for exploiting the marine resources. In this paper, we use the complex network approach to mine and analyze the interaction patterns of marine taxa and environments in spring, summer, fall and winter seasons. With the 16S rRNA pyrosequencing data of 76 time point taken monthly over 6 years, we first use our MtHc clustering algorithm to generate the operational taxonomic units (OTUs). Then, employ the k-means method to divide 76 time point samples into four seasonal groups, and utilize mutual information (MI) to construct the four correlation networks among microbial species and environment factors. Finally, we adopt the symmetrical non-negative matrix factorization method to detect the interaction patterns, and analysis the relationship between marine species and environment factors. The results show that the four seasonal microbial interaction networks have the characters of complex networks, and interaction patterns are related with the seasonal variability; the same environmental factor influences different species in the four seasons; the four environmental factors of day length, photosynthetically active radiation, NO₂+NO₃ and silicate may have stronger influences on microbes than other environment factors.     

Top‐down network analysis characterizes hidden termite–termite interactions

The analysis of ecological networks is generally bottom‐up, where networks are established by observing interactions between individuals. Emergent network properties have been indicated to reflect the dominant mode of interactions in communities that might be mutualistic (e.g., pollination) or antagonistic (e.g., host–parasitoid communities). Many ecological communities, however, comprise species interactions that are difficult to observe directly. Here, we propose that a comparison of the emergent properties from detail‐rich reference communities with known modes of interaction can inform our understanding of detail‐sparse focal communities. With this top‐down approach, we consider patterns of coexistence between termite species that live as guests in mounds built by other host termite species as a case in point. Termite societies are extremely sensitive to perturbations, which precludes determining the nature of their interactions through direct observations. We perform a literature review to construct two networks representing termite mound cohabitation in a Brazilian savanna and in the tropical forest of Cameroon. We contrast the properties of these cohabitation networks with a total of 197 geographically diverse mutualistic plant–pollinator and antagonistic host–parasitoid networks. We analyze network properties for the networks, perform a principal components analysis (PCA), and compute the Mahalanobis distance of the termite networks to the cloud of mutualistic and antagonistic networks to assess the extent to which the termite networks overlap with the properties of the reference networks. Both termite networks overlap more closely with the mutualistic plant–pollinator communities than the antagonistic host–parasitoid communities, although the Brazilian community overlap with mutualistic communities is stronger. The analysis raises the hypothesis that termite–termite cohabitation networks may be overall mutualistic. More broadly, this work provides support for the argument that cryptic communities may be analyzed via comparison to well‐characterized communities.     

Temporal variation in bat-fruit interactions: Foraging strategies influence network structure over time

Mutualistic interactions, such as seed dispersal, are important for the maintenance of structure and stability of tropical communities. However, there is a lack of information about spatial and temporal variation in plant-animal interaction networks. Thus, our goal was to assess the effect of bat's foraging strategies on temporal variation in the structure and robustness of bat-fruit networks in both a dry and a rain tropical forest. We evaluated monthly variation in bat-fruit networks by using seven structure metrics: network size, average path length, nestedness, modularity, complementary specialization, normalized degree and betweenness centrality. Seed dispersal networks showed variations in size, species composition and modularity; did not present nested structures and their complementary specialization was high compared to other studies. Both networks presented short path lengths, and a constantly high robustness, despite their monthly variations. Sedentary bat species were recorded during all the study periods and occupied more central positions than nomadic species. We conclude that foraging strategies are important structuring factors that affect the dynamic of networks by determining the functional roles of frugivorous bats over time; thus sedentary bats are more important than nomadic species for the maintenance of the network structure, and their conservation is a must.     

Investigating the Temporal Patterns within and between Intrinsic Connectivity Networks under Eyes-Open and Eyes-Closed Resting States: A Dynamical Functional Connectivity Study Based on Phase Synchronization

The brain active patterns were organized differently under resting states of eyes open (EO) and eyes closed (EC). The altered voxel-wise and regional-wise resting state active patterns under EO/EC were found by static analysis. More importantly, dynamical spontaneous functional connectivity has been observed in the resting brain. To the best of our knowledge, the dynamical mechanisms of intrinsic connectivity networks (ICNs) under EO/EC remain largely unexplored. The goals of this paper were twofold: 1) investigating the dynamical intra-ICN and inter-ICN temporal patterns during resting state; 2) analyzing the altered dynamical temporal patterns of ICNs under EO/EC. To this end, a cohort of healthy subjects with scan conditions of EO/EC were recruited from 1000 Functional Connectomes Project. Through Hilbert transform, time-varying phase synchronization (PS) was applied to evaluate the inter-ICN synchrony. Meanwhile, time-varying amplitude was analyzed as dynamical intra-ICN temporal patterns. The results found six micro-states of inter-ICN synchrony. The medial visual network (MVN) showed decreased intra-ICN amplitude during EC relative to EO. The sensory-motor network (SMN) and auditory network (AN) exhibited enhanced intra-ICN amplitude during EC relative to EO. Altered inter-ICN PS was found between certain ICNs. Particularly, the SMN and AN exhibited enhanced PS to other ICNs during EC relative to EO. In addition, the intra-ICN amplitude might influence the inter-ICN synchrony. Moreover, default mode network (DMN) might play an important role in information processing during EO/EC. Together, the dynamical temporal patterns within and between ICNs were altered during different scan conditions of EO/EC. Overall, the dynamical intra-ICN and inter-ICN temporal patterns could benefit resting state fMRI-related research, and could be potential biomarkers for human functional connectome.     

No evidence of expertise-related changes in muscle synergies during rowing

The purpose of the present study was to determine whether expertise in rowing is driven by a specific structure in muscular coordination. We compared seven experienced rowers and eight untrained (i.e., inexperienced) subjects during rowing on an ergometer. Both surface electromyography activity and mechanical patterns (forces exerted at the handle and the foot-stretcher) were recorded during a high intensity rowing exercise. A non-negative matrix factorization was applied to 23 electromyographic patterns to differentiate muscle synergies. Results showed that expertise was not associated with different dimensionality in the electromyographic data and that three muscle synergies were sufficient to explain the majority of the variance accounted for (i.e., >90% of the total variance) in the two populations. The synergies extracted were similar in the two populations, with identical functional roles. While the temporal organization of the propulsive synergies was very similar, slight differences were found in the composition of the muscle synergies (muscle synergy vectors) between the two populations. The results suggests that rowing expertise would not require the development of novel muscle synergies but would imply intrinsic synergies already used in different behaviors. Performance in rowing is more probably linked to adjustments in the mechanical output of the muscle synergies rather than to differences in the shape and timing of their activations.     

Innate synchronous oscillations in freely-organized small neuronal circuits

Background

Information processing in neuronal networks relies on the network''s ability to generate temporal patterns of action potentials. Although the nature of neuronal network activity has been intensively investigated in the past several decades at the individual neuron level, the underlying principles of the collective network activity, such as the synchronization and coordination between neurons, are largely unknown. Here we focus on isolated neuronal clusters in culture and address the following simple, yet fundamental questions: What is the minimal number of cells needed to exhibit collective dynamics? What are the internal temporal characteristics of such dynamics and how do the temporal features of network activity alternate upon crossover from minimal networks to large networks?

Methodology/Principal Findings

We used network engineering techniques to induce self-organization of cultured networks into neuronal clusters of different sizes. We found that small clusters made of as few as 40 cells already exhibit spontaneous collective events characterized by innate synchronous network oscillations in the range of 25 to 100 Hz. The oscillation frequency of each network appeared to be independent of cluster size. The duration and rate of the network events scale with cluster size but converge to that of large uniform networks. Finally, the investigation of two coupled clusters revealed clear activity propagation with master/slave asymmetry.

Conclusions/Significance

The nature of the activity patterns observed in small networks, namely the consistent emergence of similar activity across networks of different size and morphology, suggests that neuronal clusters self-regulate their activity to sustain network bursts with internal oscillatory features. We therefore suggest that clusters of as few as tens of cells can serve as a minimal but sufficient functional network, capable of sustaining oscillatory activity. Interestingly, the frequencies of these oscillations are similar those observed in vivo.