Dynamic information routing in complex networks

Neural interactions in the human frontal cortex dissociate reward and punishment learning

How human prefrontal and insular regions interact while maximizing rewards and minimizing punishments is unknown. Capitalizing on human intracranial recordings, we demonstrate that the functional specificity toward reward or punishment learning is better disentangled by interactions compared to local representations. Prefrontal and insular cortices display non-selective neural populations to rewards and punishments. Non-selective responses, however, give rise to context-specific interareal interactions. We identify a reward subsystem with redundant interactions between the orbitofrontal and ventromedial prefrontal cortices, with a driving role of the latter. In addition, we find a punishment subsystem with redundant interactions between the insular and dorsolateral cortices, with a driving role of the insula. Finally, switching between reward and punishment learning is mediated by synergistic interactions between the two subsystems. These results provide a unifying explanation of distributed cortical representations and interactions supporting reward and punishment learning.

Quantifying synergy and redundancy between networks

Understanding how different networks relate to each other is key for understanding complex systems. We introduce an intuitive yet powerful framework to disentangle different ways in which networks can be (dis)similar and complementary to each other. We decompose the shortest paths between nodes as uniquely contributed by one source network, or redundantly by either, or synergistically by both together. Our approach considers the networks' full topology, providing insights at multiple levels of resolution: from global statistics to individual paths. Our framework is widely applicable across scientific domains, from public transport to brain networks. In humans and 124 other species, we demonstrate the prevalence of unique contributions by long-range white-matter fibers in structural brain networks. Across species, efficient communication also relies on significantly greater synergy between long-range and short-range fibers than expected by chance. Our framework could find applications for designing network systems or evaluating existing ones.

Information decomposition and the informational architecture of the brain

To explain how the brain orchestrates information-processing for cognition, we must understand information itself. Importantly, information is not a monolithic entity. Information decomposition techniques provide a way to split information into its constituent elements: unique, redundant, and synergistic information. We review how disentangling synergistic and redundant interactions is redefining our understanding of integrative brain function and its neural organisation. To explain how the brain navigates the trade-offs between redundancy and synergy, we review converging evidence integrating the structural, molecular, and functional underpinnings of synergy and redundancy; their roles in cognition and computation; and how they might arise over evolution and development. Overall, disentangling synergistic and redundant information provides a guiding principle for understanding the informational architecture of the brain and cognition.

Decomposing Neural Circuit Function into Information Processing Primitives

It is challenging to measure how specific aspects of coordinated neural dynamics translate into operations of information processing and, ultimately, cognitive functions. An obstacle is that simple circuit mechanisms-such as self-sustained or propagating activity and nonlinear summation of inputs-do not directly give rise to high-level functions. Nevertheless, they already implement simple the information carried by neural activity. Here, we propose that distinct functions, such as stimulus representation, working memory, or selective attention, stem from different combinations and types of low-level manipulations of information or information processing primitives. To test this hypothesis, we combine approaches from information theory with simulations of multi-scale neural circuits involving interacting brain regions that emulate well-defined cognitive functions. Specifically, we track the information dynamics emergent from patterns of neural dynamics, using quantitative metrics to detect where and when information is actively buffered, transferred or nonlinearly merged, as possible modes of low-level processing (storage, transfer and modification). We find that neuronal subsets maintaining representations in working memory or performing attentional gain modulation are signaled by their boosted involvement in operations of information storage or modification, respectively. Thus, information dynamic metrics, beyond detecting which network units participate in cognitive processing, also promise to specify how and when they do it, that is, through which type of primitive computation, a capability that may be exploited for the analysis of experimental recordings.

Changes in functional connectivity preserve scale-free neuronal and behavioral dynamics

Does the brain optimize itself for storage and transmission of information, and if so, how? The critical brain hypothesis is based in statistical physics and posits that the brain self-tunes its dynamics to a critical point or regime to maximize the repertoire of neuronal responses. Yet, the robustness of this regime, especially with respect to changes in the functional connectivity, remains an unsolved fundamental challenge. Here, we show that both scale-free neuronal dynamics and self-similar features of behavioral dynamics persist following significant changes in functional connectivity. Specifically, we find that the psychedelic compound ibogaine that is associated with an altered state of consciousness fundamentally alters the functional connectivity in the retrosplenial cortex of mice. Yet, the scale-free statistics of movement and of neuronal avalanches among behaviorally related neurons remain largely unaltered. This indicates that the propagation of information within biological neural networks is robust to changes in functional organization of subpopulations of neurons, opening up a new perspective on how the adaptive nature of functional networks may lead to optimality of information transmission in the brain.

Perturbed Information Processing Complexity in Experimental Epilepsy

Comorbidities, such as cognitive deficits, which often accompany epilepsies, constitute a basal state, while seizures are rare and transient events. This suggests that neural dynamics, in particular those supporting cognitive function, are altered in a permanent manner in epilepsy. Here, we test the hypothesis that primitive processes of information processing at the core of cognitive function (i.e., storage and sharing of information) are altered in the hippocampus and the entorhinal cortex in experimental epilepsy in adult, male Wistar rats. We find that information storage and sharing are organized into substates across the stereotypic states of slow and theta oscillations in both epilepsy and control conditions. However, their internal composition and organization through time are disrupted in epilepsy, partially losing brain state selectivity compared with controls, and shifting toward a regimen of disorder. We propose that the alteration of information processing at this algorithmic level of computation, the theoretical intermediate level between structure and function, may be a mechanism behind the emergent and widespread comorbidities associated with epilepsy, and perhaps other disorders.SIGNIFICANCE STATEMENT Comorbidities, such as cognitive deficits, which often accompany epilepsies, constitute a basal state, while seizures are rare and transient events. This suggests that neural dynamics, in particular those supporting cognitive function, are altered in a permanent manner in epilepsy. Here, we show that basic processes of information processing at the core of cognitive function (i.e., storage and sharing of information) are altered in the hippocampus and the entorhinal cortex (two regions involved in memory processes) in experimental epilepsy. Such disruption of information processing at the algorithmic level itself could underlie the general performance impairments in epilepsy.

Spatially heterogeneous structure-function coupling in haemodynamic and electromagnetic brain networks

The relationship between structural and functional connectivity in the brain is a key question in connectomics. Here we quantify patterns of structure-function coupling across the neocortex, by comparing structural connectivity estimated using diffusion MRI with functional connectivity estimated using both neurophysiological (MEG-based) and haemodynamic (fMRI-based) recordings. We find that structure-function coupling is heterogeneous across brain regions and frequency bands. The link between structural and functional connectivity is generally stronger in multiple MEG frequency bands compared to resting state fMRI. Structure-function coupling is greater in slower and intermediate frequency bands compared to faster frequency bands. We also find that structure-function coupling systematically follows the archetypal sensorimotor-association hierarchy, as well as patterns of laminar differentiation, peaking in granular layer IV. Finally, structure-function coupling is better explained using structure-informed inter-regional communication metrics than using structural connectivity alone. Collectively, these results place neurophysiological and haemodynamic structure-function relationships in a common frame of reference and provide a starting point for a multi-modal understanding of structure-function coupling in the brain.

Co-evolving dynamics and topology in a coupled oscillator model of resting brain function

Dynamic models of ongoing BOLD fMRI brain dynamics and models of communication strategies have been two important approaches to understanding how brain network structure constrains function. However, dynamic models have yet to widely incorporate one of the most important insights from communication models: the brain may not use all of its connections in the same way or at the same time. Here we present a variation of a phase delayed Kuramoto coupled oscillator model that dynamically limits communication between nodes on each time step. An active subgraph of the empirically derived anatomical brain network is chosen in accordance with the local dynamic state on every time step, thus coupling dynamics and network structure in a novel way. We analyze this model with respect to its fit to empirical time-averaged functional connectivity, finding that, with the addition of only one parameter, it significantly outperforms standard Kuramoto models with phase delays. We also perform analyses on the novel time series of active edges it produces, demonstrating a slowly evolving topology moving through intermittent episodes of integration and segregation. We hope to demonstrate that the exploration of novel modeling mechanisms and the investigation of dynamics of networks in addition to dynamics on networks may advance our understanding of the relationship between brain structure and function.

Event-related variability is modulated by task and development

In carefully designed experimental paradigms, cognitive scientists interpret the mean event-related potentials (ERP) in terms of cognitive operations. However, the huge signal variability from one trial to the next, questions the representability of such mean events. We explored here whether this variability is an unwanted noise, or an informative part of the neural response. We took advantage of the rapid changes in the visual system during human infancy and analyzed the variability of visual responses to central and lateralized faces in 2-to 6-month-old infants compared to adults using high-density electroencephalography (EEG). We observed that neural trajectories of individual trials always remain very far from ERP components, only moderately bending their direction with a substantial temporal jitter across trials. However, single trial trajectories displayed characteristic patterns of acceleration and deceleration when approaching ERP components, as if they were under the active influence of steering forces causing transient attraction and stabilization. These dynamic events could only partly be accounted for by induced microstate transitions or phase reset phenomena. Importantly, these structured modulations of response variability, both between and within trials, had a rich sequential organization, which in infants, was modulated by the task difficulty and age. Our approaches to characterize Event Related Variability (ERV) expand on classic ERP analyses and provide the first evidence for the functional role of ongoing neural variability in human infants.

On quantification and maximization of information transfer in network dynamical systems

Information flow among nodes in a complex network describes the overall cause-effect relationships among the nodes and provides a better understanding of the contributions of these nodes individually or collectively towards the underlying network dynamics. Variations in network topologies result in varying information flows among nodes. We integrate theories from information science with control network theory into a framework that enables us to quantify and control the information flows among the nodes in a complex network. The framework explicates the relationships between the network topology and the functional patterns, such as the information transfers in biological networks, information rerouting in sensor nodes, and influence patterns in social networks. We show that by designing or re-configuring the network topology, we can optimize the information transfer function between two chosen nodes. As a proof of concept, we apply our proposed methods in the context of brain networks, where we reconfigure neural circuits to optimize excitation levels among the excitatory neurons.

The complex brain: connectivity, dynamics, information

Most would agree, the brain is complex. But, beyond metaphor, does the brain's complexity demand a paradigm shift in how we study its structure and function? I argue that complexity manifests in three domains - connectivity, dynamics, and information - and that unlocking their interactions will greatly advance our understanding of brain and cognition.

Modulations of local synchrony over time lead to resting-state functional connectivity in a parsimonious large-scale brain model

Biophysical models of large-scale brain activity are a fundamental tool for understanding the mechanisms underlying the patterns observed with neuroimaging. These models combine a macroscopic description of the within- and between-ensemble dynamics of neurons within a single architecture. A challenge for these models is accounting for modulations of within-ensemble synchrony over time. Such modulations in local synchrony are fundamental for modeling behavioral tasks and resting-state activity. Another challenge comes from the difficulty in parametrizing large scale brain models which hinders researching principles related with between-ensembles differences. Here we derive a parsimonious large scale brain model that can describe fluctuations of local synchrony. Crucially, we do not reduce within-ensemble dynamics to macroscopic variables first, instead we consider within and between-ensemble interactions similarly while preserving their physiological differences. The dynamics of within-ensemble synchrony can be tuned with a parameter which manipulates local connectivity strength. We simulated resting-state static and time-resolved functional connectivity of alpha band envelopes in models with identical and dissimilar local connectivities. We show that functional connectivity emerges when there are high fluctuations of local and global synchrony simultaneously (i.e. metastable dynamics). We also show that for most ensembles, leaning towards local asynchrony or synchrony correlates with the functional connectivity with other ensembles, with the exception of some regions belonging to the default-mode network.

Impact of basic network motifs on the collective response to perturbations

Many collective phenomena such as epidemic spreading and cascading failures in socioeconomic systems on networks are caused by perturbations of the dynamics. How perturbations propagate through networks, impact and disrupt their functions may depend on the network, the type and location of the perturbation as well as the spreading dynamics. Previous work has analyzed the retardation effects of the nodes along the propagation paths, suggesting a few transient propagation "scaling” regimes as a function of the nodes’ degree, but regardless of motifs such as triangles. Yet, empirical networks consist of motifs enabling the proper functioning of the system. Here, we show that basic motifs along the propagation path jointly determine the previously proposed scaling regimes of distance-limited propagation and degree-limited propagation, or even cease their existence. Our results suggest a radical departure from these scaling regimes and provide a deeper understanding of the interplay of self-dynamics, interaction dynamics, and topological properties.

Spreading processes and cascading failures on complex networks are often triggered by external perturbations. The authors uncover the impact of network motifs on the processes of perturbations propagation through networks, and networks’ response dynamics.

Functional control of oscillator networks

Oscillatory activity is ubiquitous in natural and engineered network systems. The interaction scheme underlying interdependent oscillatory components governs the emergence of network-wide patterns of synchrony that regulate and enable complex functions. Yet, understanding, and ultimately harnessing, the structure-function relationship in oscillator networks remains an outstanding challenge of modern science. Here, we address this challenge by presenting a principled method to prescribe exact and robust functional configurations from local network interactions through optimal tuning of the oscillators’ parameters. To quantify the behavioral synchrony between coupled oscillators, we introduce the notion of functional pattern, which encodes the pairwise relationships between the oscillators’ phases. Our procedure is computationally efficient and provably correct, accounts for constrained interaction types, and allows to concurrently assign multiple desired functional patterns. Further, we derive algebraic and graph-theoretic conditions to guarantee the feasibility and stability of target functional patterns. These conditions provide an interpretable mapping between the structural constraints and their functional implications in oscillator networks. As a proof of concept, we apply the proposed method to replicate empirically recorded functional relationships from cortical oscillations in a human brain, and to redistribute the active power flow in different models of electrical grids.

In network systems governed by oscillatory activity, such as brain networks or power grids, configurations of synchrony may define network functions. The authors introduce a control approach for the formation of desired synchrony patterns through optimal interventions on the network parameters.

A framework for macroscopic phase-resetting curves for generalised spiking neural networks

Brain rhythms emerge from synchronization among interconnected spiking neurons. Key properties of such rhythms can be gleaned from the phase-resetting curve (PRC). Inferring the PRC and developing a systematic phase reduction theory for large-scale brain rhythms remains an outstanding challenge. Here we present a theoretical framework and methodology to compute the PRC of generic spiking networks with emergent collective oscillations. We adopt a renewal approach where neurons are described by the time since their last action potential, a description that can reproduce the dynamical feature of many cell types. For a sufficiently large number of neurons, the network dynamics are well captured by a continuity equation known as the refractory density equation. We develop an adjoint method for this equation giving a semi-analytical expression of the infinitesimal PRC. We confirm the validity of our framework for specific examples of neural networks. Our theoretical framework can link key biological properties at the individual neuron scale and the macroscopic oscillatory network properties. Beyond spiking networks, the approach is applicable to a broad class of systems that can be described by renewal processes.

The formation of oscillatory neuronal assemblies at the network level has been hypothesized to be fundamental to many cognitive and motor functions. One prominent tool to understand the dynamics of oscillatory activity response to stimuli, and hence the neural code for which it is a substrate, is a nonlinear measure called Phase-Resetting Curve (PRC). At the network scale, the PRC defines the measure of how a given synaptic input perturbs the timing of next upcoming volley of spike assemblies: either advancing or delaying this timing. As a further application, one can use PRCs to make unambiguous predictions about whether communicating networks of neurons will phase-lock as it is often observed across the cortical areas and what would be this stable phase-configuration: synchronous, asynchronous or with asymmetric phase-shifts. The latter configuration also implies a preferential flow of information form the leading network to the follower, thereby giving causal signatures of directed functional connectivity. Because of the key position of the PRC in studying synchrony, information flow and entrainment to external forcing, it is crucial to move toward a theory that allows to compute the PRCs of network-wide oscillations not only for a restricted class of models, as has been done in the past, but to network descriptions that are generalized and can reflect flexibly single cell properties. In this manuscript, we tackle this issue by showing how the PRC for network oscillations can be computed using the adjoint systems of partial differential equations that define the dynamics of the neural activity density.

Oscillatory Resonance and Dynamic Manifolds in Cortical Networks with Time Delay and Multiple External Stimuli

Rhythmic oscillation is crucial for information transmission and neural communication among different brain areas. Stochastic resonance (SR) can evoke different patterns of neural oscillation. However, the characteristics of network resonance and underlying dynamical mechanisms are still unclear. In this paper, a biological model of cortical network is established and its dynamical response to external periodic stimulation is investigated. We explore the oscillatory resonance of excitatory and inhibitory populations in cortical network. It is found that the intrinsic parameters of neural populations determine the extent of resonant activity, indicating that the firing rate exhibits coherent oscillation when the frequency of external stimulation is close to intrinsic frequency of neural population. In addition, the nonlinear dynamics of cortical network in oscillatory resonance can be represented by helical manifolds in low-dimensional state space. The geometry of neural manifolds reveals the periodic dynamics and state transition in oscillatory resonance. Moreover, time delay in chemical synapses can induce multiple resonances, which appear intermittently at integer multiples of the period of input signal. The dynamical response of neural population achieves maximal periodically, due to the transition of network states induced by time delay. Furthermore, mean-field theory is applied to analyze theoretical dynamic of cortical networks with time delay and demonstrate the effective transmission of stimulation information via oscillatory resonance in the brain. Consequently, the obtained results contribute to the improvement of neuromodulation for neurological disease from the viewpoint of the neural basis.

Temperature elevations can induce switches to homoclinic action potentials that alter neural encoding and synchronization

Almost seventy years after the discovery of the mechanisms of action potential generation, some aspects of their computational consequences are still not fully understood. Based on mathematical modeling, we here explore a type of action potential dynamics – arising from a saddle-node homoclinic orbit bifurcation - that so far has received little attention. We show that this type of dynamics is to be expected by specific changes in common physiological parameters, like an elevation of temperature. Moreover, we demonstrate that it favours synchronization patterns in networks – a feature that becomes particularly prominent when system parameters change such that homoclinic spiking is induced. Supported by in-vitro hallmarks for homoclinic spikes in the rodent brain, we hypothesize that the prevalence of homoclinic spikes in the brain may be underestimated and provide a missing link between the impact of biophysical parameters on abrupt transitions between asynchronous and synchronous states of electrical activity in the brain.

The intrinsic dynamics of neurons, in particular the generation action potentials, can impact neural network states and processes of encoding information. The authors demonstrate how the elevation of temperature induces a type of action potential dynamics that favors synchronization patterns in neural networks.

Normalizing the brain connectome for communication through synchronization

Networks in neuroscience determine how brain function unfolds, and their perturbations lead to psychiatric disorders and brain disease. Brain networks are characterized by their connectomes, which comprise the totality of all connections, and are commonly described by graph theory. This approach is deeply rooted in a particle view of information processing, based on the quantification of informational bits such as firing rates. Oscillations and brain rhythms demand, however, a wave perspective of information processing based on synchronization. We extend traditional graph theory to a dual, particle-wave, perspective, integrate time delays due to finite transmission speeds, and derive a normalization of the connectome. When applied to the database of the Human Connectome Project, it explains the emergence of frequency-specific network cores including the visual and default mode networks. These findings are robust across human subjects (N = 100) and are a fundamental network property within the wave picture. The normalized connectome comprises the particle view in the limit of infinite transmission speeds and opens the applicability of graph theory to a wide range of novel network phenomena, including physiological and pathological brain rhythms. These two perspectives are orthogonal, but not incommensurable, when understood within the novel, here-proposed, generalized framework of structural connectivity.

Information evolution in complex networks

Many biological phenomena or social events critically depend on how information evolves in complex networks. However, a general theory to characterize information evolution is yet absent. Consequently, numerous unknowns remain about the mechanisms underlying information evolution. Among these unknowns, a fundamental problem, being a seeming paradox, lies in the coexistence of local randomness, manifested as the stochastic distortion of information content during individual-individual diffusion, and global regularity, illustrated by specific non-random patterns of information content on the network scale. Here, we attempt to formalize information evolution and explain the coexistence of randomness and regularity in complex networks. Applying network dynamics and information theory, we discover that a certain amount of information, determined by the selectivity of networks to the input information, frequently survives from random distortion. Other information will inevitably experience distortion or dissipation, whose speeds are shaped by the diversity of information selectivity in networks. The discovered laws exist irrespective of noise, but noise accounts for disturbing them. We further demonstrate the ubiquity of our discovered laws by analyzing the emergence of neural tuning properties in the primary visual and medial temporal cortices of animal brains and the emergence of extreme opinions in social networks.

Superconducting disordered neural networks for neuromorphic processing with fluxons

In superconductors, magnetic fields are quantized into discrete fluxons (flux quanta Φ₀), made of microscopic circulating supercurrents. We introduce a multiterminal synapse network comprising a disordered array of superconducting loops with Josephson junctions. The loops can trap fluxons defining memory, while the junctions allow their movement between loops. Dynamics of fluxons through such a disordered system through a complex reconfigurable energy landscape represents brain-like spiking information flow. In this work, we experimentally demonstrate a three-loop network using YBa₂Cu₃O_{7 - δ}-based superconducting loops and Josephson junctions, which exhibit stable memory configurations of trapped flux in loops that determine the rate of flow of fluxons through synaptic connections. The memory states are, in turn, affected by the applied input signals but can also be externally configured electrically through control current/feedback terminals. These results establish a previously unexplored, biologically similar architectural approach to neuromorphic computing that is scalable while dissipating energy of atto Joules/spike.

Mutual information resonances in delay-coupled limit cycle and quasi-cycle brain rhythms

We elucidate how coupling delays and noise impact phase and mutual information relationships between two stochastic brain rhythms. This impact depends on the dynamical regime of each PING-based rhythm, as well as on network heterogeneity and coupling asymmetry. The number of peaks at positive and negative time lags in the delayed mutual information between the two bi-directionally communicating rhythms defines our measure of flexibility of information sharing and reflects the number of ways in which the two networks can alternately lead one another. We identify two distinct mechanisms for the appearance of qualitatively similar flexible information sharing. The flexibility in the quasi-cycle regime arises from the coupling delay-induced bimodality of the phase difference distribution, and the related bimodal mutual information. It persists in the presence of asymmetric coupling and heterogeneity but is limited to two routes of information sharing. The second mechanism in noisy limit cycle regime is not induced by the delay. However, delay-coupling and heterogeneity enable communication routes at multiple time lags. Noise disrupts the shared compromise frequency, allowing the expression of individual network frequencies which leads to a slow beating pattern. Simulations of an envelope-phase description for delay-coupled quasi-cycles yield qualitatively similar properties as for the full system. Near the bifurcation from in-phase to out-of-phase behaviour, a single preferred phase difference can coexist with two information sharing routes; further, the phase laggard can be the mutual information leader, or vice versa. Overall, the coupling delay endows a two-rhythm system with an array of lead-lag relationships and mutual information resonances that exist in spite of the noise and across the Hopf bifurcation. These beg to be mapped out experimentally with the help of our predictions.

Two classes of functional connectivity in dynamical processes in networks

The relationship between network structure and dynamics is one of the most extensively investigated problems in the theory of complex systems of recent years. Understanding this relationship is of relevance to a range of disciplines-from neuroscience to geomorphology. A major strategy of investigating this relationship is the quantitative comparison of a representation of network architecture (structural connectivity, SC) with a (network) representation of the dynamics (functional connectivity, FC). Here, we show that one can distinguish two classes of functional connectivity-one based on simultaneous activity (co-activity) of nodes, the other based on sequential activity of nodes. We delineate these two classes in different categories of dynamical processes-excitations, regular and chaotic oscillators-and provide examples for SC/FC correlations of both classes in each of these models. We expand the theoretical view of the SC/FC relationships, with conceptual instances of the SC and the two classes of FC for various application scenarios in geomorphology, ecology, systems biology, neuroscience and socio-ecological systems. Seeing the organisation of dynamical processes in a network either as governed by co-activity or by sequential activity allows us to bring some order in the myriad of observations relating structure and function of complex networks.

Reconstructing network structures from partial measurements

The dynamics of systems of interacting agents is determined by the structure of their coupling network. The knowledge of the latter is, therefore, highly desirable, for instance, to develop efficient control schemes, to accurately predict the dynamics, or to better understand inter-agent processes. In many important and interesting situations, the network structure is not known, however, and previous investigations have shown how it may be inferred from complete measurement time series on each and every agent. These methods implicitly presuppose that, even though the network is not known, all its nodes are. Here, we investigate the different problem of inferring network structures within the observed/measured agents. For symmetrically coupled dynamical systems close to a stable equilibrium, we establish analytically and illustrate numerically that velocity signal correlators encode not only direct couplings, but also geodesic distances in the coupling network within the subset of measurable agents. When dynamical data are accessible for all agents, our method is furthermore algorithmically more efficient than the traditional ones because it does not rely on matrix inversion.

Stochastic facilitation in heteroclinic communication channels

Biological neural systems encode and transmit information as patterns of activity tracing complex trajectories in high-dimensional state spaces, inspiring alternative paradigms of information processing. Heteroclinic networks, naturally emerging in artificial neural systems, are networks of saddles in state space that provide a transparent approach to generate complex trajectories via controlled switches among interconnected saddles. External signals induce specific switching sequences, thus dynamically encoding inputs as trajectories. Recent works have focused either on computational aspects of heteroclinic networks, i.e., Heteroclinic Computing, or their stochastic properties under noise. Yet, how well such systems may transmit information remains an open question. Here, we investigate the information transmission properties of heteroclinic networks, studying them as communication channels. Choosing a tractable but representative system exhibiting a heteroclinic network, we investigate the mutual information rate (MIR) between input signals and the resulting sequences of states as the level of noise varies. Intriguingly, MIR does not decrease monotonically with increasing noise. Intermediate noise levels indeed maximize the information transmission capacity by promoting an increased yet controlled exploration of the underlying network of states. Complementing standard stochastic resonance, these results highlight the constructive effect of stochastic facilitation (i.e., noise-enhanced information transfer) on heteroclinic communication channels and possibly on more general dynamical systems exhibiting complex trajectories in state space.

Mechanisms of Flexible Information Sharing through Noisy Oscillations

Simple Summary

To properly interact with our environment, the brain must be able to identify external stimuli, process them, and make the right decisions all in a short time. This may involve several brain regions interacting together by sharing information birectionally via rhythmic activity. Such flexibility requires the functional connectivity between the areas to be dynamic, and a key question is the relevant parameter and operating regimes that make this possible in spite of fixed structural connectivity. Working towards this goal, we consider two coupled brain regions, each of which exhibits a noisy rhythm, a commonly observed type of neural activity. Such rhythms can be induced by the stochasticity in the neural circuitry, or be autonomously generated through nonlinearities and not necessitating noise. For these two types of rhythms, we computed the amount of information shared between the brain areas and the preferred direction(s) of sharing. We found that without the coupling delay, the flexibility needed by the brain to perform cognitive tasks requires the rhythms to be autogenerated rather than noise-induced. This is the case even with asymmetry or heterogeneity. This suggests that the importance of the dynamical regime has to be taken into account when modeling interacting neural rhythms from an information theoretical point of view.

Abstract

Brain areas must be able to interact and share information in a time-varying, dynamic manner on a fast timescale. Such flexibility in information sharing has been linked to the synchronization of rhythm phases between areas. One definition of flexibility is the number of local maxima in the delayed mutual information curve between two connected areas. However, the precise relationship between phase synchronization and information sharing is not clear, nor is the flexibility in the face of the fixed structural connectivity and noise. Here, we consider two coupled oscillatory excitatory-inhibitory networks connected through zero-delay excitatory connections, each of which mimics a rhythmic brain area. We numerically compute phase-locking and delayed mutual information between the phases of excitatory local field potential (LFPs) of the two networks, which measures the shared information and its direction. The flexibility in information sharing is shown to depend on the dynamical origin of oscillations, and its properties in different regimes are found to persist in the presence of asymmetry in the connectivity as well as system heterogeneity. For coupled noise-induced rhythms (quasi-cycles), phase synchronization is robust even in the presence of asymmetry and heterogeneity. However, they do not show flexibility, in contrast to noise-perturbed rhythms (noisy limit cycles), which are shown here to exhibit two local information maxima, i.e., flexibility. For quasi-cycles, phase difference and information measures for the envelope-phase dynamics obtained from previous analytical work using the Stochastic Averaging Method (SAM) are found to be in good qualitative agreement with those obtained from the original dynamics. The relation between phase synchronization and communication patterns is not trivial, particularly in the noisy limit cycle regime. There, complex patterns of information sharing can be observed for a single value of the phase difference. The mechanisms reported here can be extended to I-I networks since their phase synchronizations are similar. Our results set the stage for investigating information sharing between several connected noisy rhythms in neural and other complex biological networks.

Information Transmission in Delay-Coupled Neuronal Circuits in the Presence of a Relay Population

Synchronization between neuronal populations is hypothesized to play a crucial role in the communication between brain networks. The binding of features, or the association of computations occurring in spatially segregated areas, is supposed to take place when a stable synchronization between cortical areas occurs. While a direct cortico-cortical connection typically fails to support this mechanism, the participation of a third area, a relay element, mediating in the communication was proposed to overcome this limitation. Among the different structures that could play the role of coordination during the binding process, the thalamus is the best placed region to carry out this task. In this paper we study how information flows in a canonical motif that mimics a cortico-thalamo-cortical circuit composed by three mutually coupled neuronal populations (also called the V-motif). Through extensive numerical simulations, we found that the amount of information transferred between the oscillating neuronal populations is determined by the delay in their connections and the mismatch in their oscillation frequencies (detuning). While the transmission from a cortical population is mostly restricted to positive detuning, transmission from the relay (thalamic) population to the cortical populations is robust for a broad range of detuning values, including negative values, while permitting feedback communication from the cortex at high frequencies, thus supporting robust bottom up and top down interaction. In this case, a strong feedback transmission between the cortex to thalamus supports the possibility of robust bottom-up and top-down interactions in this motif. Interestingly, adding a cortico-cortical bidirectional connection to the V-motif (C-motif) expands the dynamics of the system with distinct operation modes. While overall transmission efficiency is decreased, new communication channels establish cortico-thalamo-cortical association loops. Switching between operation modes depends on the synaptic strength of the cortico-cortical connections. Our results support a role of the transthalamic V-motif in the binding of spatially segregated cortical computations, and suggest an important regulatory role of the direct cortico-cortical connection.

Virtual Connectomic Datasets in Alzheimer's Disease and Aging Using Whole-Brain Network Dynamics Modelling

Large neuroimaging datasets, including information about structural connectivity (SC) and functional connectivity (FC), play an increasingly important role in clinical research, where they guide the design of algorithms for automated stratification, diagnosis or prediction. A major obstacle is, however, the problem of missing features [e.g., lack of concurrent DTI SC and resting-state functional magnetic resonance imaging (rsfMRI) FC measurements for many of the subjects]. We propose here to address the missing connectivity features problem by introducing strategies based on computational whole-brain network modeling. Using two datasets, the Alzheimer’s Disease Neuroimaging Initiative (ADNI) dataset and a healthy aging dataset, for proof-of-concept, we demonstrate the feasibility of virtual data completion (i.e., inferring “virtual FC” from empirical SC or “virtual SC” from empirical FC), by using self-consistent simulations of linear and nonlinear brain network models. Furthermore, by performing machine learning classification (to separate age classes or control from patient subjects), we show that algorithms trained on virtual connectomes achieve discrimination performance comparable to when trained on actual empirical data; similarly, algorithms trained on virtual connectomes can be used to successfully classify novel empirical connectomes. Completion algorithms can be combined and reiterated to generate realistic surrogate connectivity matrices in arbitrarily large number, opening the way to the generation of virtual connectomic datasets with network connectivity information comparable to the one of the original data.

One-way dependent clusters and stability of cluster synchronization in directed networks

Cluster synchronization in networks of coupled oscillators is the subject of broad interest from the scientific community, with applications ranging from neural to social and animal networks and technological systems. Most of these networks are directed, with flows of information or energy that propagate unidirectionally from given nodes to other nodes. Nevertheless, most of the work on cluster synchronization has focused on undirected networks. Here we characterize cluster synchronization in general directed networks. Our first observation is that, in directed networks, a cluster A of nodes might be one-way dependent on another cluster B: in this case, A may remain synchronized provided that B is stable, but the opposite does not hold. The main contribution of this paper is a method to transform the cluster stability problem in an irreducible form. In this way, we decompose the original problem into subproblems of the lowest dimension, which allows us to immediately detect inter-dependencies among clusters. We apply our analysis to two examples of interest, a human network of violin players executing a musical piece for which directed interactions may be either activated or deactivated by the musicians, and a multilayer neural network with directed layer-to-layer connections.

Mechanisms of cluster formation in networks with directed links differ from those in undirected networks. Lodi et al. propose a method to compute interdependencies among clusters of nodes in directed networks. They show that clusters can be one-way dependent, as found in social and neural networks.

Paradoxical phase response of gamma rhythms facilitates their entrainment in heterogeneous networks

The synchronization of different γ-rhythms arising in different brain areas has been implicated in various cognitive functions. Here, we focus on the effect of the ubiquitous neuronal heterogeneity on the synchronization of ING (interneuronal network gamma) and PING (pyramidal-interneuronal network gamma) rhythms. The synchronization properties of rhythms depends on the response of their collective phase to external input. We therefore determine the macroscopic phase-response curve for finite-amplitude perturbations (fmPRC) of ING- and PING-rhythms in all-to-all coupled networks comprised of linear (IF) or quadratic (QIF) integrate-and-fire neurons. For the QIF networks we complement the direct simulations with the adjoint method to determine the infinitesimal macroscopic PRC (imPRC) within the exact mean-field theory. We show that the intrinsic neuronal heterogeneity can qualitatively modify the fmPRC and the imPRC. Both PRCs can be biphasic and change sign (type II), even though the phase-response curve for the individual neurons is strictly non-negative (type I). Thus, for ING rhythms, say, external inhibition to the inhibitory cells can, in fact, advance the collective oscillation of the network, even though the same inhibition would lead to a delay when applied to uncoupled neurons. This paradoxical advance arises when the external inhibition modifies the internal dynamics of the network by reducing the number of spikes of inhibitory neurons; the advance resulting from this disinhibition outweighs the immediate delay caused by the external inhibition. These results explain how intrinsic heterogeneity allows ING- and PING-rhythms to become synchronized with a periodic forcing or another rhythm for a wider range in the mismatch of their frequencies. Our results identify a potential function of neuronal heterogeneity in the synchronization of coupled γ-rhythms, which may play a role in neural information transfer via communication through coherence.

The interaction of a large number of oscillating units can lead to the emergence of a collective, macroscopic oscillation in which many units oscillate in near-unison or near-synchrony. This has been exploited technologically, e.g., to combine many coherently interacting, individual lasers to form a single powerful laser. Collective oscillations are also important in biology. For instance, the circadian rhythm of animals is controlled by the near-synchronous dynamics of a large number of individually oscillating cells. In animals and humans brain rhythms reflect the coherent dynamics of a large number of neurons and are surmised to play an important role in the communication between different brain areas. To be functionally relevant, these rhythms have to respond to external inputs and have to be able to synchronize with each other. We show that the ubiquitous heterogeneity in the properties of the individual neurons in a network can contribute to that ability. It can allow the external inputs to modify the internal network dynamics such that the network can follow these inputs over a wider range of frequencies. Paradoxically, while an external perturbation may delay individual neurons, their ensuing within-network interaction can overcompensate this delay, leading to an overall advance of the rhythm.

Long-term stability of avalanche scaling and integrative network organization in prefrontal and premotor cortex

Ongoing neuronal activity in the brain establishes functional networks that reflect normal and pathological brain function. Most estimates of these functional networks suffer from low spatiotemporal resolution and indirect measures of neuronal population activity, limiting the accuracy and reliability in their reconstruction over time. Here, we studied the stability of neuronal avalanche dynamics and corresponding reconstructed functional networks in the adult brain. Using chronically implanted high-density microelectrode arrays, the local field potential (LFP) of resting-state activity was recorded in prefrontal and premotor cortex of awake nonhuman primates. Avalanche dynamics revealed stable scaling exhibiting an inverted parabolic profile and collapse exponent of 2 in line with a critical branching process over many days and weeks. Functional networks were based on a Bayesian-derived estimator and demonstrated stable integrative properties characterized by nontrivial high neighborhood overlap between strongly connected nodes and robustness to weak-link pruning. Entropy-based mixing analysis revealed significant changes in strong link weights over weeks. The long-term stability in avalanche scaling and integrative network organization in the face of individual link weight changes should support the development of noninvasive biomarkers to characterize normal and abnormal brain states in the adult brain.

The brain is spontaneously active even in the absence of specific sensations or movements. This ongoing activity arises from the interactions among hundreds of thousands of neurons, which has been captured by a variety of distinct functional networks predictive of healthy and pathological brain functions. The global dynamical states and overarching network principles that underlie such ongoing brain activity are not well understood. Here we identify avalanche dynamics and “friendship” networks as two major features of ongoing activity in the frontal cortex of nonhuman primates. We demonstrate their stability over weeks in the face of local network changes. Deviation from avalanche dynamics and “friendship” organization might provide robust biomarkers to identify normal and pathological brain states.

A novel methodology to describe neuronal networks activity reveals spatiotemporal recruitment dynamics of synchronous bursting states

We propose a novel phase based analysis with the purpose of quantifying the periodic bursts of activity observed in various neuronal systems. The way bursts are intiated and propagate in a spatial network is still insufficiently characterized. In particular, we investigate here how these spatiotemporal dynamics depend on the mean connection length. We use a simplified description of a neuron's state as a time varying phase between firings. This leads to a definition of network bursts, that does not depend on the practitioner's individual judgment as the usage of subjective thresholds and time scales. This allows both an easy and objective characterization of the bursting dynamics, only depending on system's proper scales. Our approach thus ensures more reliable and reproducible measurements. We here use it to describe the spatiotemporal processes in networks of intrinsically oscillating neurons. The analysis rigorously reveals the role of the mean connectivity length in spatially embedded networks in determining the existence of "leader" neurons during burst initiation, a feature incompletely understood observed in several neuronal cultures experiments. The precise definition of a burst with our method allowed us to rigorously characterize the initiation dynamics of bursts and show how it depends on the mean connectivity length. Although presented with simulations, the methodology can be applied to other forms of neuronal spatiotemporal data. As shown in a preliminary study with MEA recordings, it is not limited to in silico modeling.

Coherence resonance in neuronal populations: Mean-field versus network model

The counterintuitive phenomenon of coherence resonance describes a nonmonotonic behavior of the regularity of noise-induced oscillations in the excitable regime, leading to an optimal response in terms of regularity of the excited oscillations for an intermediate noise intensity. We study this phenomenon in populations of FitzHugh-Nagumo (FHN) neurons with different coupling architectures. For networks of FHN systems in an excitable regime, coherence resonance has been previously analyzed numerically. Here we focus on an analytical approach studying the mean-field limits of the globally and locally coupled populations. The mean-field limit refers to an averaged behavior of a complex network as the number of elements goes to infinity. We apply the mean-field approach to the globally coupled FHN network. Further, we derive a mean-field limit approximating the locally coupled FHN network with low noise intensities. We study the effects of the coupling strength and noise intensity on coherence resonance for both the network and the mean-field models. We compare the results of the mean-field and network frameworks and find good agreement in the globally coupled case, where the correspondence between the two approaches is sufficiently good to capture the emergence of coherence resonance, as well as of anticoherence resonance.

Transmission delays and frequency detuning can regulate information flow between brain regions

Brain networks exhibit very variable and dynamical functional connectivity and flexible configurations of information exchange despite their overall fixed structure. Brain oscillations are hypothesized to underlie time-dependent functional connectivity by periodically changing the excitability of neural populations. In this paper, we investigate the role of the connection delay and the detuning between the natural frequencies of neural populations in the transmission of signals. Based on numerical simulations and analytical arguments, we show that the amount of information transfer between two oscillating neural populations could be determined by their connection delay and the mismatch in their oscillation frequencies. Our results highlight the role of the collective phase response curve of the oscillating neural populations for the efficacy of signal transmission and the quality of the information transfer in brain networks.

Dynamic core-periphery structure of information sharing networks in entorhinal cortex and hippocampus

Neural computation is associated with the emergence, reconfiguration, and dissolution of cell assemblies in the context of varying oscillatory states. Here, we describe the complex spatiotemporal dynamics of cell assemblies through temporal network formalism. We use a sliding window approach to extract sequences of networks of information sharing among single units in hippocampus and entorhinal cortex during anesthesia and study how global and node-wise functional connectivity properties evolve through time and as a function of changing global brain state (theta vs. slow-wave oscillations). First, we find that information sharing networks display, at any time, a core-periphery structure in which an integrated core of more tightly functionally interconnected units links to more loosely connected network leaves. However the units participating to the core or to the periphery substantially change across time windows, with units entering and leaving the core in a smooth way. Second, we find that discrete network states can be defined on top of this continuously ongoing liquid core-periphery reorganization. Switching between network states results in a more abrupt modification of the units belonging to the core and is only loosely linked to transitions between global oscillatory states. Third, we characterize different styles of temporal connectivity that cells can exhibit within each state of the sharing network. While inhibitory cells tend to be central, we show that, otherwise, anatomical localization only poorly influences the patterns of temporal connectivity of the different cells. Furthermore, cells can change temporal connectivity style when the network changes state. Altogether, these findings reveal that the sharing of information mediated by the intrinsic dynamics of hippocampal and entorhinal cortex cell assemblies have a rich spatiotemporal structure, which could not have been identified by more conventional time- or state-averaged analyses of functional connectivity.

It is generally thought that computations performed by local brain circuits rely on complex neural processes, associated with the flexible waxing and waning of cell assemblies, that is, an ensemble of cells firing in tight synchrony. Although cell assembly formation is inherently and unavoidably dynamical, it is still common to find studies in which essentially “static” approaches are used to characterize this process. In the present study, we adopt instead a temporal network approach. Avoiding usual time-averaging procedures, we reveal that hub neurons are not hardwired but that cells vary smoothly their degree of integration within the assembly core. Furthermore, our temporal network framework enables the definition of alternative possible styles of “hubness.” Some cells may share information with a multitude of other units but only in an intermittent manner, as “activists” in a flash mob. In contrast, some other cells may share information in a steadier manner, as resolute “lobbyists.” Finally, by avoiding averages over preimposed states, we show that within each global oscillatory state rich switching dynamics can take place between a repertoire of many available network states. We thus show that the temporal network framework provides a natural and effective language to rigorously describe the rich spatiotemporal patterns of information sharing instantiated by cell assembly evolution.

Time-evolving controllability of effective connectivity networks during seizure progression

Over one third of the estimated 3 million people with epilepsy in the United States are medication resistant. Responsive neurostimulation from chronically implanted electrodes provides a promising treatment alternative to resective surgery. However, determining optimal personalized stimulation parameters, including when and where to intervene to guarantee a positive patient outcome, is a major open challenge. Network neuroscience and control theory offer useful tools that may guide improvements in parameter selection for control of anomalous neural activity. Here we use a method to characterize dynamic controllability across consecutive effective connectivity (EC) networks based on regularized partial correlations between implanted electrodes during the onset, propagation, and termination regimes of 34 seizures. We estimate regularized partial correlation adjacency matrices from 1-s time windows of intracranial electrocorticography recordings using the Graphical Least Absolute Shrinkage and Selection Operator (GLASSO). Average and modal controllability metrics calculated from each resulting EC network track the time-varying controllability of the brain on an evolving landscape of conditionally dependent network interactions. We show that average controllability increases throughout a seizure and is negatively correlated with modal controllability throughout. Our results support the hypothesis that the energy required to drive the brain to a seizure-free state from an ictal state is smallest during seizure onset, yet we find that applying control energy at electrodes in the seizure onset zone may not always be energetically favorable. Our work suggests that a low-complexity model of time-evolving controllability may offer insights for developing and improving control strategies targeting seizure suppression.

Birth and destruction of collective oscillations in a network of two populations of coupled type 1 neurons

We study the macroscopic dynamics of large networks of excitable type 1 neurons composed of two populations interacting with disparate but symmetric intra- and inter-population coupling strengths. This nonuniform coupling scheme facilitates symmetric equilibria, where both populations display identical firing activity, characterized by either quiescent or spiking behavior, or asymmetric equilibria, where the firing activity of one population exhibits quiescent but the other exhibits spiking behavior. Oscillations in the firing rate are possible if neurons emit pulses with non-zero width but are otherwise quenched. Here, we explore how collective oscillations emerge for two statistically identical neuron populations in the limit of an infinite number of neurons. A detailed analysis reveals how collective oscillations are born and destroyed in various bifurcation scenarios and how they are organized around higher codimension bifurcation points. Since both symmetric and asymmetric equilibria display bistable behavior, a large configuration space with steady and oscillatory behavior is available. Switching between configurations of neural activity is relevant in functional processes such as working memory and the onset of collective oscillations in motor control.

Temporal complexity of fMRI is reproducible and correlates with higher order cognition

It has been hypothesized that resting state networks (RSNs), extracted from resting state functional magnetic resonance imaging (rsfMRI), likely display unique temporal complexity fingerprints, quantified by their multiscale entropy patterns (McDonough and Nashiro, 2014). This is a hypothesis with a potential capacity for developing digital biomarkers of normal brain function, as well as pathological brain dysfunction. Nevertheless, a limitation of McDonough and Nashiro (2014) was that rsfMRI data from only 20 healthy individuals was used for the analysis. To validate this hypothesis in a larger cohort, we used rsfMRI datasets of 987 healthy young adults from the Human Connectome Project (HCP), aged 22-35, each with four 14.4-min rsfMRI recordings and parcellated into 379 brain regions. We quantified multiscale entropy of rsfMRI time series averaged at different cortical and sub-cortical regions. We performed effect-size analysis on the data in 8 RSNs. Given that the morphology of multiscale entropy is affected by the choice of its tolerance parameter (r) and embedding dimension (m), we repeated the analyses at multiple values of r and m including the values used in McDonough and Nashiro (2014). Our results reinforced high temporal complexity in the default mode and frontoparietal networks. Lowest temporal complexity was observed in the subcortical areas and limbic system. We investigated the effect of temporal resolution (determined by the repetition time T_R) after downsampling of rsfMRI time series at two rates. At a low temporal resolution, we observed increased entropy and variance across datasets. Test-retest analysis showed that findings were likely reproducible across individuals over four rsfMRI runs, especially when the tolerance parameter r is equal to 0.5. The results confirmed that the relationship between functional brain connectivity strengths and rsfMRI temporal complexity changes over time scales. Finally, a non-random correlation was observed between temporal complexity of RSNs and fluid intelligence suggesting that complex dynamics of the human brain is an important attribute of high-level brain function.

Dynamic Functional Connectivity as a complex random walk: Definitions and the dFCwalk toolbox

•We have developed a framework to describe the dynamics of Functional Connectivity (dFC) estimated from brain activity time-series as a complex random walk in the space of possible functional networks. This conceptual and methodological framework considers dFC as a smooth reconfiguration process, combining "liquid" and "coordinated" aspects. Unlike other previous approaches, our method does not require the explicit extraction of discrete connectivity states.•In our previous work, we introduced several metrics for the quantitative characterization of the dFC random walk. First, dFC speed analyses extract the distribution of the time-resolved rate of reconfiguration of FC along time. These distributions have a clear peak (typical dFC speed, that can already serve as a biomarker) and fat tails (denoting deviations from Gaussianity that can be detected by suitable scaling analyses of FC network streams). Second, meta-connectivity (MC) analyses identify groups of functional links whose fluctuations co-vary in time and that define veritable dFC modules organized along specific dFC meta-hub controllers (differing from conventional FC modules and hubs). The decomposition of whole-brain dFC by MC allows performing dFC speed analyses separately for each of the detected dFC modules.•We present here blocks and pipelines for dFC random walk analyses that are made easily available through a dedicated MATLABⓇ toolbox (dFCwalk), openly downloadable. Although we applied such analyses mostly to fMRI resting state data, in principle our methods can be extended to any type of neural activity (from Local Field Potentials to EEG, MEG, fNIRS, etc.) or even non-neural time-series.

Relations between large-scale brain connectivity and effects of regional stimulation depend on collective dynamical state

At the macroscale, the brain operates as a network of interconnected neuronal populations, which display coordinated rhythmic dynamics that support interareal communication. Understanding how stimulation of different brain areas impacts such activity is important for gaining basic insights into brain function and for further developing therapeutic neurmodulation. However, the complexity of brain structure and dynamics hinders predictions regarding the downstream effects of focal stimulation. More specifically, little is known about how the collective oscillatory regime of brain network activity—in concert with network structure—affects the outcomes of perturbations. Here, we combine human connectome data and biophysical modeling to begin filling these gaps. By tuning parameters that control collective system dynamics, we identify distinct states of simulated brain activity and investigate how the distributed effects of stimulation manifest at different dynamical working points. When baseline oscillations are weak, the stimulated area exhibits enhanced power and frequency, and due to network interactions, activity in this excited frequency band propagates to nearby regions. Notably, beyond these linear effects, we further find that focal stimulation causes more distributed modifications to interareal coherence in a band containing regions’ baseline oscillation frequencies. Importantly, depending on the dynamical state of the system, these broadband effects can be better predicted by functional rather than structural connectivity, emphasizing a complex interplay between anatomical organization, dynamics, and response to perturbation. In contrast, when the network operates in a regime of strong regional oscillations, stimulation causes only slight shifts in power and frequency, and structural connectivity becomes most predictive of stimulation-induced changes in network activity patterns. In sum, this work builds upon and extends previous computational studies investigating the impacts of stimulation, and underscores the fact that both the stimulation site, and, crucially, the regime of brain network dynamics, can influence the network-wide responses to local perturbations.

Stimulation can be used to alter brain activity and is a therapeutic option for certain neurological conditions. However, predicting the distributed effects of local perturbations is difficult. Previous studies show that responses to stimulation depend on anatomical (or structural) coupling. In addition to structure, here we consider how stimulation effects also depend on the brain’s collective dynamical (or functional) state, arising from the coordination of rhythmic activity across large-scale networks. In a whole-brain computational model, we show that global responses to regional stimulation can indeed be contingent upon and differ across various dynamical working points. Notably, depending on the network’s oscillatory regime, stimulation can accelerate the activity of the stimulated site, and lead to widespread effects at both the new, excited frequency, as well as in a much broader frequency range including areas’ baseline frequencies. While structural connectivity is a good predictor of “excited band” changes, in some states “baseline band” effects can be better predicted by functional connectivity, which depends upon the system’s oscillatory regime. By integrating and extending past efforts, our results thus indicate that dynamical—in additional to structural—brain organization plays a role in governing how focal stimulation modulates interactions between distributed network elements.

Modular slowing of resting-state dynamic functional connectivity as a marker of cognitive dysfunction induced by sleep deprivation

Dynamic Functional Connectivity (dFC) in the resting state (rs) is considered as a correlate of cognitive processing. Describing dFC as a flow across morphing connectivity configurations, our notion of dFC speed quantifies the rate at which FC networks evolve in time. Here we probe the hypothesis that variations of rs dFC speed and cognitive performance are selectively interrelated within specific functional subnetworks. In particular, we focus on Sleep Deprivation (SD) as a reversible model of cognitive dysfunction. We found that whole-brain level (global) dFC speed significantly slows down after 24h of SD. However, the reduction in global dFC speed does not correlate with variations of cognitive performance in individual tasks, which are subtle and highly heterogeneous. On the contrary, we found strong correlations between performance variations in individual tasks -including Rapid Visual Processing (RVP, assessing sustained visual attention)- and dFC speed quantified at the level of functional sub-networks of interest. Providing a compromise between classic static FC (no time) and global dFC (no space), modular dFC speed analyses allow quantifying a different speed of dFC reconfiguration independently for sub-networks overseeing different tasks. Importantly, we found that RVP performance robustly correlates with the modular dFC speed of a characteristic frontoparietal module.

Dynamic Functional Connectivity between order and randomness and its evolution across the human adult lifespan

Functional Connectivity (FC) during resting-state or task conditions is not static but inherently dynamic. Yet, there is no consensus on whether fluctuations in FC may resemble isolated transitions between discrete FC states rather than continuous changes. This quarrel hampers advancing the study of dynamic FC. This is unfortunate as the structure of fluctuations in FC can certainly provide more information about developmental changes, aging, and progression of pathologies. We merge the two perspectives and consider dynamic FC as an ongoing network reconfiguration, including a stochastic exploration of the space of possible steady FC states. The statistical properties of this random walk deviate both from a purely "order-driven" dynamics, in which the mean FC is preserved, and from a purely "randomness-driven" scenario, in which fluctuations of FC remain uncorrelated over time. Instead, dynamic FC has a complex structure endowed with long-range sequential correlations that give rise to transient slowing and acceleration epochs in the continuous flow of reconfiguration. Our analysis for fMRI data in healthy elderly revealed that dynamic FC tends to slow down and becomes less complex as well as more random with increasing age. These effects appear to be strongly associated with age-related changes in behavioural and cognitive performance.

Cortical ignition dynamics is tightly linked to the core organisation of the human connectome

The capability of cortical regions to flexibly sustain an "ignited" state of activity has been discussed in relation to conscious perception or hierarchical information processing. Here, we investigate how the intrinsic propensity of different regions to get ignited is determined by the specific topological organisation of the structural connectome. More specifically, we simulated the resting-state dynamics of mean-field whole-brain models and assessed how dynamic multistability and ignition differ between a reference model embedding a realistic human connectome, and alternative models based on a variety of randomised connectome ensembles. We found that the strength of global excitation needed to first trigger ignition in a subset of regions is substantially smaller for the model embedding the empirical human connectome. Furthermore, when increasing the strength of excitation, the propagation of ignition outside of this initial core-which is able to self-sustain its high activity-is way more gradual than for any of the randomised connectomes, allowing for graded control of the number of ignited regions. We explain both these assets in terms of the exceptional weighted core-shell organisation of the empirical connectome, speculating that this topology of human structural connectivity may be attuned to support enhanced ignition dynamics.

Dynamic representations in networked neural systems

A group of neurons can generate patterns of activity that represent information about stimuli; subsequently, the group can transform and transmit activity patterns across synapses to spatially distributed areas. Recent studies in neuroscience have begun to independently address the two components of information processing: the representation of stimuli in neural activity and the transmission of information in networks that model neural interactions. Yet only recently are studies seeking to link these two types of approaches. Here we briefly review the two separate bodies of literature; we then review the recent strides made to address this gap. We continue with a discussion of how patterns of activity evolve from one representation to another, forming dynamic representations that unfold on the underlying network. Our goal is to offer a holistic framework for understanding and describing neural information representation and transmission while revealing exciting frontiers for future research.

Control of brain network dynamics across diverse scales of space and time

The human brain is composed of distinct regions that are each associated with particular functions and distinct propensities for the control of neural dynamics. However, the relation between these functions and control profiles is poorly understood, as is the variation in this relation across diverse scales of space and time. Here we probe the relation between control and dynamics in brain networks constructed from diffusion tensor imaging data in a large community sample of young adults. Specifically, we probe the control properties of each brain region and investigate their relationship with dynamics across various spatial scales using the Laplacian eigenspectrum. In addition, through analysis of regional modal controllability and partitioning of modes, we determine whether the associated dynamics are fast or slow, as well as whether they are alternating or monotone. We find that brain regions that facilitate the control of energetically easy transitions are associated with activity on short length scales and slow timescales. Conversely, brain regions that facilitate control of difficult transitions are associated with activity on long length scales and fast timescales. Built on linear dynamical models, our results offer parsimonious explanations for the activity propagation and network control profiles supported by regions of differing neuroanatomical structure.

Efficient communication over complex dynamical networks: The role of matrix non-normality

In both natural and engineered systems, communication often occurs dynamically over networks ranging from highly structured grids to largely disordered graphs. To use, or comprehend the use of, networks as efficient communication media requires understanding of how they propagate and transform information in the face of noise. Here, we develop a framework that enables us to examine how network structure, noise, and interference between consecutive packets jointly determine transmission performance in complex networks governed by linear dynamics. Mathematically, normal networks, which can be decomposed into separate low-dimensional information channels, suffer greatly from readout noise. Most details of their wiring have no impact on transmission quality. Non-normal networks, however, can largely cancel the effect of noise by transiently amplifying select input dimensions while ignoring others, resulting in higher net information throughput. Our theory could inform the design of new communication networks, as well as the optimal use of existing ones.

Determinants of Brain Rhythm Burst Statistics

Brain rhythms recorded in vivo, such as gamma oscillations, are notoriously variable both in amplitude and frequency. They are characterized by transient epochs of higher amplitude known as bursts. It has been suggested that, despite their short-life and random occurrence, bursts in gamma and other rhythms can efficiently contribute to working memory or communication tasks. Abnormalities in bursts have also been associated with e.g. motor and psychiatric disorders. It is thus crucial to understand how single cell and connectivity parameters influence burst statistics and the corresponding brain states. To address this problem, we consider a generic stochastic recurrent network of Pyramidal Interneuron Network Gamma (PING) type. Using the stochastic averaging method, we derive dynamics for the phase and envelope of the amplitude process, and find that they depend on only two meta-parameters that combine all the model parameters. This allows us to identify an optimal parameter regime of healthy variability with similar statistics to those seen in vivo; in this regime, oscillations and bursts are supported by synaptic noise. The probability density for the rhythm’s envelope as well as the mean burst duration are then derived using first passage time analysis. Our analysis enables us to link burst attributes, such as duration and frequency content, to system parameters. Our general approach can be extended to different frequency bands, network topologies and extra populations. It provides the much needed insight into the biophysical determinants of rhythm burst statistics, and into what needs to be changed to correct rhythms with pathological statistics.

Parallelism Strategies for Big Data Delayed Transfer Entropy Evaluation

Generated and collected data have been rising with the popularization of technologies such as Internet of Things, social media, and smartphone, leading big data term creation. One class of big data hidden information is causality. Among the tools to infer causal relationships, there is Delay Transfer Entropy (DTE); however, it has a high demanding processing power. Many approaches were proposed to overcome DTE performance issues such as GPU and FPGA implementations. Our study compared different parallel strategies to calculate DTE from big data series using a heterogeneous Beowulf cluster. Task Parallelism was significantly faster in comparison to Data Parallelism. With big data trend in sight, these results may enable bigger datasets analysis or better statistical evidence.

Fluctuation-induced distributed resonances in oscillatory networks

Across physics, biology, and engineering, the collective dynamics of oscillatory networks often evolve into self-organized operating states. How such networks respond to external fluctuating signals fundamentally underlies their function, yet is not well understood. Here, we present a theory of dynamic network response patterns and reveal how distributed resonance patterns emerge in oscillatory networks once the dynamics of the oscillatory units become more than one-dimensional. The network resonances are topology specific and emerge at an intermediate frequency content of the input signals, between global yet homogeneous responses at low frequencies and localized responses at high frequencies. Our analysis reveals why these patterns arise and where in the network they are most prominent. These results may thus provide general theoretical insights into how fluctuating signals induce response patterns in networked systems and simultaneously help to develop practical guiding principles for real-world network design and control.