Ising-like dynamics in large-scale functional brain networks

Diagnostically distinct resting state fMRI energy distributions: A subject-specific maximum entropy modeling study

Objective

Existing neuroimaging studies of psychotic and mood disorders have reported brain activation differences (first-order properties) and altered pairwise correlation-based functional connectivity (second-order properties). However, both approaches have certain limitations that can be overcome by integrating them in a pairwise maximum entropy model (MEM) that better represents a comprehensive picture of fMRI signal patterns and provides a system-wide summary measure called energy. This study examines the applicability of individual-level MEM for psychiatry and identifies image-derived model coefficients related to model parameters.

Method

MEMs are fit to resting state fMRI data from each individual with schizophrenia/schizoaffective disorder, bipolar disorder, and major depression (n=132) and demographically matched healthy controls (n=132) from the UK Biobank to different subsets of the default mode network (DMN) regions.

Results

The model satisfactorily explained observed brain energy state occurrence probabilities across all participants, and model parameters were significantly correlated with image-derived coefficients for all groups. Within clinical groups, averaged energy level distributions were higher in schizophrenia/schizoaffective disorder but lower in bipolar disorder compared to controls for both bilateral and unilateral DMN. Major depression energy distributions were higher compared to controls only in the right hemisphere DMN.

Conclusions

Diagnostically distinct energy states suggest that probability distributions of temporal changes in synchronously active nodes may underlie each diagnostic entity. Subject-specific MEMs allow for factoring in the individual variations compared to traditional group-level inferences, offering an improved measure of biologically meaningful correlates of brain activity that may have potential clinical utility.

Coarse graining and criticality in the human connectome

In the face of the stupefying complexity of the human brain, network analysis is a most useful tool that allows one to greatly simplify the problem, typically by approximating the billions of neurons making up the brain by means of a coarse-grained picture with a practicable number of nodes. But even such relatively small and coarse networks, such as the human connectome with its 100-1000 nodes, may present challenges for some computationally demanding analyses that are incapable of handling networks with more than a handful of nodes. With such applications in mind, we set out to study the extent to which dynamical behavior and critical phenomena in the brain may be preserved following a severe coarse-graining procedure. Thus we proceeded to further coarse grain the human connectome by taking a modularity-based approach, the goal being to produce a network of a relatively small number of modules. After finding that the qualitative dynamical behavior of the coarse-grained networks reflected that of the original networks, albeit to a less pronounced extent, we then formulated a hypothesis based on the coarse-grained networks in the context of criticality in the Wilson-Cowan and Ising models, and we verified the hypothesis, which connected a transition value of the former with the critical temperature of the latter, using the original networks. This preservation of dynamical and critical behavior following a severe coarse-graining procedure, in principle, allows for the drawing of similar qualitative conclusions by analyzing much smaller networks, which opens the door for studying the human connectome in contexts typically regarded as computationally intractable, such as Integrated Information Theory and quantum models of the human brain.

Stationary stable cross-correlation pattern and task specific deviations in unresponsive wakefulness syndrome as well as clinically healthy subjects

Brain dynamics is highly non-stationary, permanently subject to ever-changing external conditions and continuously monitoring and adjusting internal control mechanisms. Finding stationary structures in this system, as has been done recently, is therefore of great importance for understanding fundamental dynamic trade relationships. Here we analyse electroencephalographic recordings (EEG) of 13 subjects with unresponsive wakefulness syndrome (UWS) during rest and while being influenced by different acoustic stimuli. We compare the results with a control group under the same experimental conditions and with clinically healthy subjects during overnight sleep. The main objective of this study is to investigate whether a stationary correlation pattern is also present in the UWS group, and if so, to what extent this structure resembles the one found in healthy subjects. Furthermore, we extract transient dynamical features via specific deviations from the stationary interrelation pattern. We find that (i) the UWS group is more heterogeneous than the two groups of healthy subjects, (ii) also the EEGs of the UWS group contain a stationary cross-correlation pattern, although it is less pronounced and shows less similarity to that found for healthy subjects and (iii) deviations from the stationary pattern are notably larger for the UWS than for the two groups of healthy subjects. The results suggest that the nervous system of subjects with UWS receive external stimuli but show an overreaching reaction to them, which may disturb opportune information processing.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

Differences in Brain Network Topology Based on Alcohol Use History in Adolescents

Approximately 6 million youth aged 12 to 20 consume alcohol monthly in the United States. The effect of alcohol consumption in adolescence on behavior and cognition is heavily researched; however, little is known about how alcohol consumption in adolescence may alter brain function, leading to long-term developmental detriments. In order to investigate differences in brain connectivity associated with alcohol use in adolescents, brain networks were constructed using resting-state functional magnetic resonance imaging data collected by the National Consortium on Alcohol and NeuroDevelopment in Adolescence (NCANDA) from 698 youth (12-21 years; 117 hazardous drinkers and 581 no/low drinkers). Analyses assessed differences in brain network topology based on alcohol consumption in eight predefined brain networks, as well as in whole-brain connectivity. Within the central executive network (CEN), basal ganglia network (BGN), and sensorimotor network (SMN), no/low drinkers demonstrated stronger and more frequent connections between highly globally efficient nodes, with fewer and weaker connections between highly clustered nodes. Inverse results were observed within the dorsal attention network (DAN), visual network (VN), and frontotemporal network (FTN), with no/low drinkers demonstrating weaker connections between nodes with high efficiency and increased frequency of clustered nodes compared to hazardous drinkers. Cross-sectional results from this study show clear organizational differences between adolescents with no/low or hazardous alcohol use, suggesting that aberrant connectivity in these brain networks is associated with risky drinking behaviors.

Improvements in task performance after practice are associated with scale-free dynamics of brain activity

Although practicing a task generally benefits later performance on that same task, there are individual differences in practice effects. One avenue to model such differences comes from research showing that brain networks extract functional advantages from operating in the vicinity of criticality, a state in which brain network activity is more scale-free. We hypothesized that higher scale-free signal from fMRI data, measured with the Hurst exponent (H), indicates closer proximity to critical states. We tested whether individuals with higher H during repeated task performance would show greater practice effects. In Study 1, participants performed a dual-n-back task (DNB) twice during MRI (n = 56). In Study 2, we used two runs of n-back task (NBK) data from the Human Connectome Project sample (n = 599). In Study 3, participants performed a word completion task (CAST) across six runs (n = 44). In all three studies, multivariate analysis was used to test whether higher H was related to greater practice-related performance improvement. Supporting our hypothesis, we found patterns of higher H that reliably correlated with greater performance improvement across participants in all three studies. However, the predictive brain regions were distinct, suggesting that the specific spatial H↑ patterns are not task-general.

Investigating structural and functional aspects of the brain's criticality in stroke

This paper addresses the question of the brain's critical dynamics after an injury such as a stroke. It is hypothesized that the healthy brain operates near a phase transition (critical point), which provides optimal conditions for information transmission and responses to inputs. If structural damage could cause the critical point to disappear and thus make self-organized criticality unachievable, it would offer the theoretical explanation for the post-stroke impairment of brain function. In our contribution, however, we demonstrate using network models of the brain, that the dynamics remain critical even after a stroke. In cases where the average size of the second-largest cluster of active nodes, which is one of the commonly used indicators of criticality, shows an anomalous behavior, it results from the loss of integrity of the network, quantifiable within graph theory, and not from genuine non-critical dynamics. We propose a new simple model of an artificial stroke that explains this anomaly. The proposed interpretation of the results is confirmed by an analysis of real connectomes acquired from post-stroke patients and a control group. The results presented refer to neurobiological data; however, the conclusions reached apply to a broad class of complex systems that admit a critical state.

Global excitability and network structure in the human brain

We utilize a model of Wilson-Cowan oscillators to investigate structure-function relationships in the human brain by means of simulations of the spontaneous dynamics of brain networks generated through human connectome data. This allows us to establish relationships between the global excitability of such networks and global structural network quantities for connectomes of two different sizes for a number of individual subjects. We compare the qualitative behavior of such correlations between biological networks and shuffled networks, the latter generated by shuffling the pairwise connectivities of the former while preserving their distribution. Our results point towards a remarkable propensity of the brain to achieve a trade-off between low network wiring cost and strong functionality, and highlight the unique capacity of brain network topologies to exhibit a strong transition from an inactive state to a globally excited one.

LSD-induced increase of Ising temperature and algorithmic complexity of brain dynamics

A topic of growing interest in computational neuroscience is the discovery of fundamental principles underlying global dynamics and the self-organization of the brain. In particular, the notion that the brain operates near criticality has gained considerable support, and recent work has shown that the dynamics of different brain states may be modeled by pairwise maximum entropy Ising models at various distances from a phase transition, i.e., from criticality. Here we aim to characterize two brain states (psychedelics-induced and placebo) as captured by functional magnetic resonance imaging (fMRI), with features derived from the Ising spin model formalism (system temperature, critical point, susceptibility) and from algorithmic complexity. We hypothesized, along the lines of the entropic brain hypothesis, that psychedelics drive brain dynamics into a more disordered state at a higher Ising temperature and increased complexity. We analyze resting state blood-oxygen-level-dependent (BOLD) fMRI data collected in an earlier study from fifteen subjects in a control condition (placebo) and during ingestion of lysergic acid diethylamide (LSD). Working with the automated anatomical labeling (AAL) brain parcellation, we first create "archetype" Ising models representative of the entire dataset (global) and of the data in each condition. Remarkably, we find that such archetypes exhibit a strong correlation with an average structural connectome template obtained from dMRI (r = 0.6). We compare the archetypes from the two conditions and find that the Ising connectivity in the LSD condition is lower than in the placebo one, especially in homotopic links (interhemispheric connectivity), reflecting a significant decrease of homotopic functional connectivity in the LSD condition. The global archetype is then personalized for each individual and condition by adjusting the system temperature. The resulting temperatures are all near but above the critical point of the model in the paramagnetic (disordered) phase. The individualized Ising temperatures are higher in the LSD condition than in the placebo condition (p = 9 × 10-5). Next, we estimate the Lempel-Ziv-Welch (LZW) complexity of the binarized BOLD data and the synthetic data generated with the individualized model using the Metropolis algorithm for each participant and condition. The LZW complexity computed from experimental data reveals a weak statistical relationship with condition (p = 0.04 one-tailed Wilcoxon test) and none with Ising temperature (r(13) = 0.13, p = 0.65), presumably because of the limited length of the BOLD time series. Similarly, we explore complexity using the block decomposition method (BDM), a more advanced method for estimating algorithmic complexity. The BDM complexity of the experimental data displays a significant correlation with Ising temperature (r(13) = 0.56, p = 0.03) and a weak but significant correlation with condition (p = 0.04, one-tailed Wilcoxon test). This study suggests that the effects of LSD increase the complexity of brain dynamics by loosening interhemispheric connectivity-especially homotopic links. In agreement with earlier work using the Ising formalism with BOLD data, we find the brain state in the placebo condition is already above the critical point, with LSD resulting in a shift further away from criticality into a more disordered state.

Evidence of Critical Dynamics in Movements of Bees inside a Hive

Social insects such as honey bees exhibit complex behavioral patterns, and their distributed behavioral coordination enables decision-making at the colony level. It has, therefore, been proposed that a high-level description of their collective behavior might share commonalities with the dynamics of neural processes in brains. Here, we investigated this proposal by focusing on the possibility that brains are poised at the edge of a critical phase transition and that such a state is enabling increased computational power and adaptability. We applied mathematical tools developed in computational neuroscience to a dataset of bee movement trajectories that were recorded within the hive during the course of many days. We found that certain characteristics of the activity of the bee hive system are consistent with the Ising model when it operates at a critical temperature, and that the system's behavioral dynamics share features with the human brain in the resting state.

How critical is brain criticality?

Criticality is the singular state of complex systems poised at the brink of a phase transition between order and randomness. Such systems display remarkable information-processing capabilities, evoking the compelling hypothesis that the brain may itself be critical. This foundational idea is now drawing renewed interest thanks to high-density data and converging cross-disciplinary knowledge. Together, these lines of inquiry have shed light on the intimate link between criticality, computation, and cognition. Here, we review these emerging trends in criticality neuroscience, highlighting new data pertaining to the edge of chaos and near-criticality, and making a case for the distance to criticality as a useful metric for probing cognitive states and mental illness. This unfolding progress in the field contributes to establishing criticality theory as a powerful mechanistic framework for studying emergent function and its efficiency in both biological and artificial neural networks.

Astrocytes Learn to Detect and Signal Deviations from Critical Brain Dynamics

Astrocytes are nonneuronal brain cells that were recently shown to actively communicate with neurons and are implicated in memory, learning, and regulation of cognitive states. Interestingly, these information processing functions are also closely linked to the brain's ability to self-organize at a critical phase transition. Investigating the mechanistic link between astrocytes and critical brain dynamics remains beyond the reach of cellular experiments, but it becomes increasingly approachable through computational studies. We developed a biologically plausible computational model of astrocytes to analyze how astrocyte calcium waves can respond to changes in underlying network dynamics. Our results suggest that astrocytes detect synaptic activity and signal directional changes in neuronal network dynamics using the frequency of their calcium waves. We show that this function may be facilitated by receptor scaling plasticity by enabling astrocytes to learn the approximate information content of input synaptic activity. This resulted in a computationally simple, information-theoretic model, which we demonstrate replicating the signaling functionality of the biophysical astrocyte model with receptor scaling. Our findings provide several experimentally testable hypotheses that offer insight into the regulatory role of astrocytes in brain information processing.

Critical Spatial-Temporal Dynamics and Prominent Shape Collapse of Calcium Waves Observed in Human hNT Astrocytes in Vitro

Networks of neurons are typically studied in the field of Criticality. However, the study of astrocyte networks in the brain has been recently lauded to be of equal importance to that of the neural networks. To date criticality assessments have only been performed on networks astrocytes from healthy rats, and astrocytes from cultured dissociated resections of intractable epilepsy. This work, for the first time, presents studies of the critical dynamics and shape collapse of calcium waves observed in cultures of healthy human astrocyte networks in vitro, derived from the human hNT cell line. In this article, we demonstrate that avalanches of spontaneous calcium waves display strong critical dynamics, including power-laws in both the size and duration distributions. In addition, the temporal profiles of avalanches displayed self-similarity, leading to shape collapse of the temporal profiles. These findings are significant as they suggest that cultured networks of healthy human hNT astrocytes self-organize to a critical point, implying that healthy astrocytic networks operate at a critical point to process and transmit information. Furthermore, this work can serve as a point of reference to which other astrocyte criticality studies can be compared.

A regression framework for brain network distance metrics

Analyzing brain networks has long been a prominent research topic in neuroimaging. However, statistical methods to detect differences between these networks and relate them to phenotypic traits are still sorely needed. Our previous work developed a novel permutation testing framework to detect differences between two groups. Here we advance that work to allow both assessing differences by continuous phenotypes and controlling for confounding variables. To achieve this, we propose an innovative regression framework to relate distances (or similarities) between brain network features to functions of absolute differences in continuous covariates and indicators of difference for categorical variables. We explore several similarity metrics for comparing distances (or similarities) between connection matrices, and adapt several standard methods for estimation and inference within our framework: standard F test, F test with individual level effects (ILE), feasible generalized least squares (FGLS), and permutation. Via simulation studies, we assess all approaches for estimation and inference while comparing them with existing multivariate distance matrix regression (MDMR) methods. We then illustrate the utility of our framework by analyzing the relationship between fluid intelligence and brain network distances in Human Connectome Project (HCP) data.

Stationary EEG pattern relates to large-scale resting state networks - An EEG-fMRI study connecting brain networks across time-scales

Relating brain dynamics acting on time scales that differ by at least an order of magnitude is a fundamental issue in brain research. The same is true for the observation of stable dynamical structures in otherwise highly non-stationary signals. The present study addresses both problems by the analysis of simultaneous resting state EEG-fMRI recordings of 53 patients with epilepsy. Confirming previous findings, we observe a generic and temporally stable average correlation pattern in EEG recordings. We design a predictor for the General Linear Model describing fluctuations around the stationary EEG correlation pattern and detect resting state networks in fMRI data. The acquired statistical maps are contrasted to several surrogate tests and compared with maps derived by spatial Independent Component Analysis of the fMRI data. By means of the proposed EEG-predictor we observe core nodes of known fMRI resting state networks with high specificity in the default mode, the executive control and the salience network. Our results suggest that both, the stationary EEG pattern as well as resting state fMRI networks are different expressions of the same brain activity. This activity is interpreted as the dynamics on (or close to) a stable attractor in phase space that is necessary to maintain the brain in an efficient operational mode. We discuss that this interpretation is congruent with the theoretical framework of complex systems as well as with the brain's energy balance.

Brain Connectivity Studies on Structure-Function Relationships: A Short Survey with an Emphasis on Machine Learning

This short survey reviews the recent literature on the relationship between the brain structure and its functional dynamics. Imaging techniques such as diffusion tensor imaging (DTI) make it possible to reconstruct axonal fiber tracks and describe the structural connectivity (SC) between brain regions. By measuring fluctuations in neuronal activity, functional magnetic resonance imaging (fMRI) provides insights into the dynamics within this structural network. One key for a better understanding of brain mechanisms is to investigate how these fast dynamics emerge on a relatively stable structural backbone. So far, computational simulations and methods from graph theory have been mainly used for modeling this relationship. Machine learning techniques have already been established in neuroimaging for identifying functionally independent brain networks and classifying pathological brain states. This survey focuses on methods from machine learning, which contribute to our understanding of functional interactions between brain regions and their relation to the underlying anatomical substrate.

A model integrating multiple processes of synchronization and coherence for information instantiation within a cortical area

What is the form of dynamic, e.g., sensory, information in the mammalian cortex? Information in the cortex is modeled as a coherence map of a mixed chimera state of synchronous, phasic, and disordered minicolumns. The theoretical model is built on neurophysiological evidence. Complex spatiotemporal information is instantiated through a system of interacting biological processes that generate a synchronized cortical area, a coherent aperture. Minicolumn elements are grouped in macrocolumns in an array analogous to a phased-array radar, modeled as an aperture, a "hole through which radiant energy flows." Coherence maps in a cortical area transform inputs from multiple sources into outputs to multiple targets, while reducing complexity and entropy. Coherent apertures can assume extremely large numbers of different information states as coherence maps, which can be communicated among apertures with corresponding very large bandwidths. The coherent aperture model incorporates considerable reported research, integrating five conceptually and mathematically independent processes: 1) a damped Kuramoto network model, 2) a pumped area field potential, 3) the gating of nearly coincident spikes, 4) the coherence of activity across cortical lamina, and 5) complex information formed through functions in macrocolumns. Biological processes and their interactions are described in equations and a functional circuit such that the mathematical pieces can be assembled the same way the neurophysiological ones are. The model can be conceptually convolved over the specifics of local cortical areas within and across species. A coherent aperture becomes a node in a graph of cortical areas with a corresponding distribution of information.

A comparison of diffusion tractography techniques in simulating the generalized Ising model to predict the intrinsic activity of the brain

Diffusion tractography is a non-invasive technique that is being used to estimate the location and direction of white matter tracts in the brain. Identifying the characteristics of white matter plays an important role in research as well as in clinical practice that relies on finding the relationship between the structure and function of the brain. An Ising model implemented on a structural connectivity (SC) has proven to explain the spontaneous fluctuations in the brain at criticality using brain's structure depicted by white matter tracts. Since the SC is the only input of the model, identifying the tractography technique which provides a SC that delivers the highest prediction of the brain's intrinsic activity via the generalized Ising model (GIM) is essential. Hence an Ising model is simulated on SCs generated using two different acquisition schemes (single and multi-shell) and two different tractography approaches (deterministic and probabilistic) and analyzed at criticality across 69 healthy subjects. Results showed that by introducing the GIM, predictability of the empirical correlation matrix increases on average from 0.2 to 0.6 compared to the predictability using the empirical connectivity matrix directly. It is also observed that the SC generated using deterministic tractography without fractional anisotropy resulted in the highest correlation coefficient value of 0.65 between the simulated and empirical correlation matrices. Additionally, calculated dimensionalities per simulation illustrated that the dimensionality depends upon the method of tractography that has been used to extract the SC.

Modeling an auditory stimulated brain under altered states of consciousness using the generalized Ising model

Propofol is a short-acting medication that results in decreased levels of consciousness and is used for general anesthesia. Although it is the most commonly used anesthetic in the world, much remains unknown about the mechanisms by which it induces a loss of consciousness. Characterizing anesthesia-induced alterations to brain network activity might provide a powerful framework for understanding the neural mechanisms of unconsciousness. The aim of this work was to model brain activity in healthy brains during various stages of consciousness, as induced by propofol, in the auditory paradigm. We used the generalized Ising model (GIM) to fit the empirical fMRI data of healthy subjects while they listened to an audio clip from a movie. The external stimulus (audio clip) is believed to be at least partially driving a synchronization process of the brain activity and provides a similar conscious experience in different subjects. In order to observe the common synchronization among the subjects, a novel technique called the inter subject correlation (ISC) was implemented. We showed that the GIM-modified to incorporate the naturalistic external field-was able to fit the empirical task fMRI data in the awake state, in mild sedation, in deep sedation, and in recovery, at a temperature T* which is well above the critical temperature. To our knowledge this is the first study that captures human brain activity in response to real-life external stimuli at different levels of conscious awareness using mathematical modeling. This study might be helpful in the future to assess the level of consciousness of patients with disorders of consciousness and help in regaining their consciousness.

Why Brain Criticality Is Clinically Relevant: A Scoping Review

The past 25 years have seen a strong increase in the number of publications related to criticality in different areas of neuroscience. The potential of criticality to explain various brain properties, including optimal information processing, has made it an increasingly exciting area of investigation for neuroscientists. Recent reviews on this topic, sometimes termed brain criticality, make brief mention of clinical applications of these findings to several neurological disorders such as epilepsy, neurodegenerative disease, and neonatal hypoxia. Other clinicallyrelevant domains - including anesthesia, sleep medicine, developmental-behavioral pediatrics, and psychiatry - are seldom discussed in review papers of brain criticality. Thorough assessments of these application areas and their relevance for clinicians have also yet to be published. In this scoping review, studies of brain criticality involving human data of all ages are evaluated for their current and future clinical relevance. To make the results of these studies understandable to a more clinical audience, a review of the key concepts behind criticality (e.g., phase transitions, long-range temporal correlation, self-organized criticality, power laws, branching processes) precedes the discussion of human clinical studies. Open questions and forthcoming areas of investigation are also considered.

The Emergence of Integrated Information, Complexity, and 'Consciousness' at Criticality

Integrated Information Theory (IIT) posits that integrated information ( Φ ) represents the quantity of a conscious experience. Here, the generalized Ising model was used to calculate Φ as a function of temperature in toy models of fully connected neural networks. A Monte-Carlo simulation was run on 159 normalized, random, positively weighted networks analogous to small five-node excitatory neural network motifs. Integrated information generated by this sample of small Ising models was measured across model parameter spaces. It was observed that integrated information, as an order parameter, underwent a phase transition at the critical point in the model. This critical point was demarcated by the peak of the generalized susceptibility (or variance in configuration due to temperature) of integrated information. At this critical point, integrated information was maximally receptive and responsive to perturbations of its own states. The results of this study provide evidence that Φ can capture integrated information in an empirical dataset, and display critical behavior acting as an order parameter from the generalized Ising model.

Closer to critical resting-state neural dynamics in individuals with higher fluid intelligence

According to the critical brain hypothesis, the brain is considered to operate near criticality and realize efficient neural computations. Despite the prior theoretical and empirical evidence in favor of the hypothesis, no direct link has been provided between human cognitive performance and the neural criticality. Here we provide such a key link by analyzing resting-state dynamics of functional magnetic resonance imaging (fMRI) networks at a whole-brain level. We develop a data-driven analysis method, inspired from statistical physics theory of spin systems, to map out the whole-brain neural dynamics onto a phase diagram. Using this tool, we show evidence that neural dynamics of human participants with higher fluid intelligence quotient scores are closer to a critical state, i.e., the boundary between the paramagnetic phase and the spin-glass (SG) phase. The present results are consistent with the notion of "edge-of-chaos" neural computation.

Emergence of power laws in noncritical neuronal systems

Experimental and computational studies provide compelling evidence that neuronal systems are characterized by power-law distributions of neuronal avalanche sizes. This fact is interpreted as an indication that these systems are operating near criticality, and, in turn, typical properties of critical dynamical processes, such as optimal information transmission and stability, are attributed to neuronal systems. The purpose of this Rapid Communication is to show that the presence of power-law distributions for the size of neuronal avalanches is not a sufficient condition for the system to operate near criticality. Specifically, we consider a simplistic model of neuronal dynamics on networks and show that the degree distribution of the underlying neuronal network may trigger power-law distributions for neuronal avalanches even when the system is not in its critical regime. To certify and explain our findings we develop an analytical approach based on percolation theory and branching processes techniques.

Modeling functional resting-state brain networks through neural message passing on the human connectome

In this work, we propose a natural model for information flow in the brain through a neural message-passing dynamics on a structural network of macroscopic regions, such as the human connectome (HC). In our model, each brain region is assumed to have a binary behavior (active or not), the strengths of interactions among them are encoded in the anatomical connectivity matrix defined by the HC, and the dynamics of the system is defined by the Belief Propagation (BP) algorithm, working near the critical point of the network. We show that in the absence of direct external stimuli the BP algorithm converges to a spatial map of activations that is similar to the Default Mode Network (DMN) of the brain, which has been defined from the analysis of functional MRI data. Moreover, we use Susceptibility Propagation (SP) to compute the matrix of long-range correlations between the different regions and show that the modules defined by a clustering of this matrix resemble several Resting State Networks (RSN) determined experimentally. Both results suggest that the functional DMN and RSNs can be seen as simple consequences of the anatomical structure of the brain and a neural message-passing dynamics between macroscopic regions. With the new model, we explore predictions on how functional maps change when the anatomical brain network suffers structural alterations, like in Alzheimer's disease and in lesions of the Corpus Callosum. The implications and novel interpretations suggested by the model, as well as the role of criticality, are discussed.

Non-equilibrium critical dynamics of bursts in θ and δ rhythms as fundamental characteristic of sleep and wake micro-architecture

Origin and functions of intermittent transitions among sleep stages, including short awakenings and arousals, constitute a challenge to the current homeostatic framework for sleep regulation, focusing on factors modulating sleep over large time scales. Here we propose that the complex micro-architecture characterizing the sleep-wake cycle results from an underlying non-equilibrium critical dynamics, bridging collective behaviors across spatio-temporal scales. We investigate θ and δ wave dynamics in control rats and in rats with lesions of sleep-promoting neurons in the parafacial zone. We demonstrate that intermittent bursts in θ and δ rhythms exhibit a complex temporal organization, with long-range power-law correlations and a robust duality of power law (θ-bursts, active phase) and exponential-like (δ-bursts, quiescent phase) duration distributions, typical features of non-equilibrium systems self-organizing at criticality. Crucially, such temporal organization relates to anti-correlated coupling between θ- and δ-bursts, and is independent of the dominant physiologic state and lesions, a solid indication of a basic principle in sleep dynamics.

Avalanches and criticality in self-organized nanoscale networks

Current efforts to achieve neuromorphic computation are focused on highly organized architectures, such as integrated circuits and regular arrays of memristors, which lack the complex interconnectivity of the brain and so are unable to exhibit brain-like dynamics. New architectures are required, both to emulate the complexity of the brain and to achieve critical dynamics and consequent maximal computational performance. We show here that electrical signals from self-organized networks of nanoparticles exhibit brain-like spatiotemporal correlations and criticality when fabricated at a percolating phase transition. Specifically, the sizes and durations of avalanches of switching events are power law distributed, and the power law exponents satisfy rigorous criteria for criticality. These signals are therefore qualitatively and quantitatively similar to those measured in the cortex. Our self-organized networks provide a low-cost platform for computational approaches that rely on spatiotemporal correlations, such as reservoir computing, and are an important step toward creating neuromorphic device architectures.

The avalanche-like behaviour of large-scale haemodynamic activity from wakefulness to deep sleep

Increasing evidence suggests that responsiveness is associated with critical or near-critical cortical dynamics, which exhibit scale-free cascades of spatio-temporal activity. These cascades, or 'avalanches', have been detected at multiple scales, from in vitro and in vivo microcircuits to voltage imaging and brain-wide functional magnetic resonance imaging (fMRI) recordings. Criticality endows the cortex with certain information-processing capacities postulated as necessary for conscious wakefulness, yet it remains unknown how unresponsiveness impacts on the avalanche-like behaviour of large-scale human haemodynamic activity. We observed a scale-free hierarchy of co-activated connected clusters by applying a point-process transformation to fMRI data recorded during wakefulness and non-rapid eye movement (NREM) sleep. Maximum-likelihood estimates revealed a significant effect of sleep stage on the scaling parameters of the cluster size power-law distributions. Post hoc statistical tests showed that differences were maximal between wakefulness and N2 sleep. These results were robust against spatial coarse graining, fitting alternative statistical models and different point-process thresholds, and disappeared upon phase shuffling the fMRI time series. Evoked neural bistabilities preventing arousals during N2 sleep do not suffice to explain these differences, which point towards changes in the intrinsic dynamics of the brain that could be necessary to consolidate a state of deep unresponsiveness.

Critical excitation-inhibition balance in dense neural networks

The "edge of chaos" phase transition in artificial neural networks is of renewed interest in light of recent evidence for criticality in brain dynamics. Statistical mechanics traditionally studied this transition with connectivity k as the control parameter and an exactly balanced excitation-inhibition ratio. While critical connectivity has been found to be low in these model systems, typically around k=2, which is unrealistic for natural neural systems, a recent study utilizing the excitation-inhibition ratio as the control parameter found a new, nearly degree independent, critical point when connectivity is large. However, the new phase transition is accompanied by an unnaturally high level of activity in the network. Here we study random neural networks with the additional properties of (i) a high clustering coefficient and (ii) neurons that are solely either excitatory or inhibitory, a prominent property of natural neurons. As a result, we observe an additional critical point for networks with large connectivity, regardless of degree distribution, which exhibits low activity levels that compare well with neuronal brain networks.

Eigenvector Centrality Dynamics From Resting-State fMRI: Gender and Age Differences in Healthy Subjects

With the increasing use of functional brain network properties as markers of brain disorders, efficient visualization and evaluation methods have become essential. Eigenvector centrality mapping (ECM) of functional MRI (fMRI) data enables the representation of per-node graph theoretical measures as brain maps. This paper studies the use of centrality dynamics for measuring group differences in imaging studies. Imaging data were used from a publicly available imaging study, which included resting fMRI data. After warping the images to a standard space and masking cortical regions, ECM were computed in a sliding window. The dual regression method was used to identify dynamic centrality differences inside well-known resting-state networks between gender and age groups. Gender-related differences were found in the medial and lateral visual, motor, default mode, and executive control RSN, where male subjects had more consistent centrality variations within the network. Age-related differences between the youngest and oldest subjects, based on a median split, were found in the medial visual, executive control and left frontoparietal networks, where younger subjects had more consistent centrality variations within the network. Our findings show that centrality dynamics can be used to identify between-group functional brain network centrality differences, and that age and gender distributions studies need to be taken into account in functional imaging studies.

Retrieval of retrained and reconsolidated memories are associated with a distinct neural network

Consolidated memories can persist from a single day to years, and persistence is improved by retraining or retrieval-mediated plasticity. One retrieval-based way to strengthen memory is the reconsolidation process. Strengthening occurs simply by the presentation of specific cues associated with the original learning. This enhancement function has a fundamental role in the maintenance of memory relevance in animals everyday life. In the present study, we made a step forward in the identification of brain correlates imprinted by the reconsolidation process studying the long-term neural consequences when the strengthened memory is stable again. To reach such a goal, we compared the retention of paired-associate memories that went through retraining process or were labilizated-reconsolidated. Using functional magnetic resonance imaging (fMRI), we studied the specific areas activated during retrieval and analyzed the functional connectivity of the whole brain associated with the event-related design. We used Graph Theory tools to analyze the global features of the network. We show that reconsolidated memories imprint a more locally efficient network that is better at exchanging information, compared with memories that were retrained or untreated. For the first time, we report a method to elucidate the neural footprints associated with a relevant function of memory reconsolidation.

Homeostatic plasticity and emergence of functional networks in a whole-brain model at criticality

Understanding the relationship between large-scale structural and functional brain networks remains a crucial issue in modern neuroscience. Recently, there has been growing interest in investigating the role of homeostatic plasticity mechanisms, across different spatiotemporal scales, in regulating network activity and brain functioning against a wide range of environmental conditions and brain states (e.g., during learning, development, ageing, neurological diseases). In the present study, we investigate how the inclusion of homeostatic plasticity in a stochastic whole-brain model, implemented as a normalization of the incoming node's excitatory input, affects the macroscopic activity during rest and the formation of functional networks. Importantly, we address the structure-function relationship both at the group and individual-based levels. In this work, we show that normalization of the node's excitatory input improves the correspondence between simulated neural patterns of the model and various brain functional data. Indeed, we find that the best match is achieved when the model control parameter is in its critical value and that normalization minimizes both the variability of the critical points and neuronal activity patterns among subjects. Therefore, our results suggest that the inclusion of homeostatic principles lead to more realistic brain activity consistent with the hallmarks of criticality. Our theoretical framework open new perspectives in personalized brain modeling with potential applications to investigate the deviation from criticality due to structural lesions (e.g. stroke) or brain disorders.

A MATLAB toolbox for multivariate analysis of brain networks

Complex brain networks formed via structural and functional interactions among brain regions are believed to underlie information processing and cognitive function. A growing number of studies indicate that altered brain network topology is associated with physiological, behavioral, and cognitive abnormalities. Graph theory is showing promise as a method for evaluating and explaining brain networks. However, multivariate frameworks that provide statistical inferences about how such networks relate to covariates of interest, such as disease phenotypes, in different study populations are yet to be developed. We have developed a freely available MATLAB toolbox with a graphical user interface that bridges this important gap between brain network analyses and statistical inference. The modeling framework implemented in this toolbox utilizes a mixed-effects multivariate regression framework that allows assessing brain network differences between study populations as well as assessing the effects of covariates of interest such as age, disease phenotype, and risk factors on the density and strength of brain connections in global (i.e., whole-brain) and local (i.e., subnetworks) brain networks. Confounding variables, such as sex, are controlled for through the implemented framework. A variety of neuroimaging data such as fMRI, EEG, and DTI can be analyzed with this toolbox, which makes it useful for a wide range of studies examining the structure and function of brain networks. The toolbox uses SAS, R, or Python (depending on software availability) to perform the statistical modeling. We also provide a clustering-based data reduction method that helps with model convergence and substantially reduces modeling time for large data sets.

Robustness of functional networks at criticality against structural defects

The robustness of dynamical properties of neuronal networks against structural damages is a central problem in computational and experimental neuroscience. Research has shown that the cortical network of a healthy brain works near a critical state and, moreover, that functional neuronal networks often have scale-free and small-world properties. In this work, we study how the robustness of simple functional networks at criticality is affected by structural defects. In particular, we consider a two-dimensional Ising model at the critical temperature and investigate how its functional network changes with the increasing degree of structural defects. We show that the scale-free and small-world properties of the functional network at criticality are robust against large degrees of structural lesions while the system remains below the percolation limit. Although the Ising model is only a conceptual description of a two-state neuron, our research reveals fundamental robustness properties of functional networks derived from classical statistical mechanics models.

Role of Dimensionality in Predicting the Spontaneous Behavior of the Brain Using the Classical Ising Model and the Ising Model Implemented on a Structural Connectome

There is accumulating evidence that spontaneous fluctuations of the brain are sustained by a structural architecture of axonal fiber bundles. Various models have been used to investigate this structure–function relationship. In this work, we implemented the Ising model using the number of fibers between each pair of brain regions as input. The output of the Ising model simulations on a structural connectome was then compared with empirical functional connectivity data. A simpler two-dimensional classical Ising model was used as the baseline model for comparison purpose. Thermodynamic properties, such as the magnetic susceptibility and the specific heat, illustrated a phase transition from an ordered phase to a disordered phase at the critical temperature. Despite the differences between the two models, the lattice Ising model and the Ising model implemented on a structural connectome (the generalized Ising model) exhibited similar patterns of global properties. To study the behavior of the generalized Ising model around criticality, calculation of the dimensionality and critical exponents was performed for the first time, by introducing a new concept of distance based on structural connectivity. Same value inside the fitting error was found for the dimensionality in both models suggesting similar behavior of the models around criticality.

Frequency-Resolved Dynamic Functional Connectivity Reveals Scale-Stable Features of Connectivity-States

Investigating temporal variability of functional connectivity is an emerging field in connectomics. Entering dynamic functional connectivity by applying sliding window techniques on resting-state fMRI (rs-fMRI) time courses emerged from this topic. We introduce frequency-resolved dynamic functional connectivity (frdFC) by means of multivariate empirical mode decomposition (MEMD) followed up by filter-bank investigations. In general, we find that MEMD is capable of generating time courses to perform frdFC and we discover that the structure of connectivity-states is robust over frequency scales and even becomes more evident with decreasing frequency. This scale-stability varies with the number of extracted clusters when applying k-means. We find a scale-stability drop-off from k = 4 to k = 5 extracted connectivity-states, which is corroborated by null-models, simulations, theoretical considerations, filter-banks, and scale-adjusted windows. Our filter-bank studies show that filter design is more delicate in the rs-fMRI than in the simulated case. Besides offering a baseline for further frdFC research, we suggest and demonstrate the use of scale-stability as a possible quality criterion for connectivity-state and model selection. We present first evidence showing that connectivity-states are both a multivariate, and a multiscale phenomenon. A data repository of our frequency-resolved time-series is provided.

Adaptive scales of integration and response latencies in a critically-balanced model of the primary visual cortex

The primary visual cortex (V1) integrates information over scales in visual space, which have been shown to vary, in an input-dependent manner, as a function of contrast and other visual parameters. Which algorithms the brain uses to achieve this feat are largely unknown and an open problem in visual neuroscience. We demonstrate that a simple dynamical mechanism can account for this contrast-dependent scale of integration in visuotopic space as well as connect this property to two other stimulus-dependent features of V1: extents of lateral integration on the cortical surface and response latencies.

Altered avalanche dynamics in a developmental NMDAR hypofunction model of cognitive impairment

Disturbed activity patterns in cortical networks contribute to the pathophysiology of schizophrenia (SZ). Several lines of evidence implicate NMDA receptor hypofunction in SZ, and blocking NMDA receptor signaling during early neurodevelopment produces cognitive deficits in rodent models that resemble those seen in schizophrenic patients. However, the altered network dynamics underlying these cognitive impairments largely remain to be characterized, especially at the cellular level. Here, we use in vivo two-photon calcium imaging to describe pathological dynamics, occurring in parallel with cognitive dysfunction, in a developmental NMDA receptor hypofunction model. We observed increased synchrony and specific alterations in spatiotemporal activity propagation, which could be causally linked to a previously unidentified persistent bursting phenotype. This phenotype was rescued by acute treatment with the NMDA receptor co-agonist D-serine or the GABA_B receptor agonist baclofen, which similarly rescued working memory performance. It was not reproduced by optogenetic inhibition of fast-spiking interneurons. These results provide novel insight into network-level abnormalities mediating the cognitive impairment induced by NMDA receptor hypofunction.

Optimized connectome architecture for sensory-motor integration

The intricate connectivity patterns of neural circuits support a wide repertoire of communication processes and functional interactions. Here we systematically investigate how neural signaling is constrained by anatomical connectivity in the mesoscale Drosophila (fruit fly) brain network. We use a spreading model that describes how local perturbations, such as external stimuli, trigger global signaling cascades that spread through the network. Through a series of simple biological scenarios we demonstrate that anatomical embedding potentiates sensory-motor integration. We find that signal spreading is faster from nodes associated with sensory transduction (sensors) to nodes associated with motor output (effectors). Signal propagation was accelerated if sensor nodes were activated simultaneously, suggesting a topologically mediated synergy among sensors. In addition, the organization of the network increases the likelihood of convergence of multiple cascades towards effector nodes, thereby facilitating integration prior to motor output. Moreover, effector nodes tend to coactivate more frequently than other pairs of nodes, suggesting an anatomically enhanced coordination of motor output. Altogether, our results show that the organization of the mesoscale Drosophila connectome imparts privileged, behaviorally relevant communication patterns among sensors and effectors, shaping their capacity to collectively integrate information.

The complex network spanned by neurons and their axonal projections promotes a diverse set of functions. In the present report, we study how the topological organization of the fruit fly brain supports sensory-motor integration. Using a simple communication model, we demonstrate that the topology of this network allows efficient coordination among sensory and motor neurons. Our results suggest that brain network organization may profoundly shape the functional repertoire of this simple organism.

A multi-variate blind source separation algorithm

BACKGROUND AND OBJECTIVE

The study follows the proposal of decomposing a given data matrix into a product of independent spatial and temporal component matrices. A multi-variate decomposition approach is presented, based on an approximate diagonalization of a set of matrices computed using a latent space representation.

METHODS

The proposed methodology follows an algebraic approach, which is common to space, temporal or spatiotemporal blind source separation algorithms. More specifically, the algebraic approach relies on singular value decomposition techniques, which avoids computationally costly and numerically instable matrix inversion. The method is equally applicable to correlation matrices determined from second order correlations or by considering fourth order correlations.

RESULTS

The resulting algorithms are applied to fMRI data sets either to extract the underlying fMRI components or to extract connectivity maps from resting state fMRI data collected for a dynamic functional connectivity analysis. Intriguingly, our algorithm shows increased spatial specificity compared to common approaches, while temporal precision stays similar.

CONCLUSION

The study presents a novel spatiotemporal blind source separation algorithm, which is both robust and avoids parameters that are difficult to fine tune. Applied on experimental data sets, the new method yields highly confined and focused areas with least spatial extent in the retinotopy case, and similar results in the dynamic functional connectivity analyses compared to other blind source separation algorithms. Therefore, we conclude that our novel algorithm is highly competitive and yields results, which are superior or at least similar to existing approaches.

Criticality predicts maximum irregularity in recurrent networks of excitatory nodes

A rigorous understanding of brain dynamics and function requires a conceptual bridge between multiple levels of organization, including neural spiking and network-level population activity. Mounting evidence suggests that neural networks of cerebral cortex operate at a critical regime, which is defined as a transition point between two phases of short lasting and chaotic activity. However, despite the fact that criticality brings about certain functional advantages for information processing, its supporting evidence is still far from conclusive, as it has been mostly based on power law scaling of size and durations of cascades of activity. Moreover, to what degree such hypothesis could explain some fundamental features of neural activity is still largely unknown. One of the most prevalent features of cortical activity in vivo is known to be spike irregularity of spike trains, which is measured in terms of the coefficient of variation (CV) larger than one. Here, using a minimal computational model of excitatory nodes, we show that irregular spiking (CV > 1) naturally emerges in a recurrent network operating at criticality. More importantly, we show that even at the presence of other sources of spike irregularity, being at criticality maximizes the mean coefficient of variation of neurons, thereby maximizing their spike irregularity. Furthermore, we also show that such a maximized irregularity results in maximum correlation between neuronal firing rates and their corresponding spike irregularity (measured in terms of CV). On the one hand, using a model in the universality class of directed percolation, we propose new hallmarks of criticality at single-unit level, which could be applicable to any network of excitable nodes. On the other hand, given the controversy of the neural criticality hypothesis, we discuss the limitation of this approach to neural systems and to what degree they support the criticality hypothesis in real neural networks. Finally, we discuss the limitations of applying our results to real networks and to what degree they support the criticality hypothesis.

Energy landscape analysis of neuroimaging data

Computational neuroscience models have been used for understanding neural dynamics in the brain and how they may be altered when physiological or other conditions change. We review and develop a data-driven approach to neuroimaging data called the energy landscape analysis. The methods are rooted in statistical physics theory, in particular the Ising model, also known as the (pairwise) maximum entropy model and Boltzmann machine. The methods have been applied to fitting electrophysiological data in neuroscience for a decade, but their use in neuroimaging data is still in its infancy. We first review the methods and discuss some algorithms and technical aspects. Then, we apply the methods to functional magnetic resonance imaging data recorded from healthy individuals to inspect the relationship between the accuracy of fitting, the size of the brain system to be analysed and the data length.This article is part of the themed issue 'Mathematical methods in medicine: neuroscience, cardiology and pathology'.

The impacts of pesticide and nicotine exposures on functional brain networks in Latino immigrant workers

Latino immigrants that work on farms experience chronic exposures to potential neurotoxicants, such as pesticides, as part of their work. For tobacco farmworkers there is the additional risk of exposure to moderate to high doses of nicotine. Pesticide and nicotine exposures have been associated with neurological changes in the brain. Long-term exposure to cholinesterase-inhibiting pesticides, such as organophosphates and carbamates, and nicotine place this vulnerable population at risk for developing neurological dysfunction. In this study we examined whole-brain connectivity patterns and brain network properties of Latino immigrant workers. Comparisons were made between farmworkers and non-farmworkers using resting-state functional magnetic resonance imaging data and a mixed-effects modeling framework. We also evaluated how measures of pesticide and nicotine exposures contributed to the findings. Our results indicate that despite having the same functional connectivity density and strength, brain networks in farmworkers had more clustered and modular structures when compared to non-farmworkers. Our findings suggest increased functional specificity and decreased functional integration in farmworkers when compared to non-farmworkers. Cholinesterase activity was associated with population differences in community structure and the strength of brain network functional connections. Urinary cotinine, a marker of nicotine exposure, was associated with the differences in network community structure. Brain network differences between farmworkers and non-farmworkers, as well as pesticide and nicotine exposure effects on brain functional connections in this study, may illuminate underlying mechanisms that cause neurological implications in later life.

Ising model with conserved magnetization on the human connectome: Implications on the relation structure-function in wakefulness and anesthesia

Dynamical models implemented on the large scale architecture of the human brain may shed light on how a function arises from the underlying structure. This is the case notably for simple abstract models, such as the Ising model. We compare the spin correlations of the Ising model and the empirical functional brain correlations, both at the single link level and at the modular level, and show that their match increases at the modular level in anesthesia, in line with recent results and theories. Moreover, we show that at the peak of the specific heat (the critical state), the spin correlations are minimally shaped by the underlying structural network, explaining how the best match between the structure and function is obtained at the onset of criticality, as previously observed. These findings confirm that brain dynamics under anesthesia shows a departure from criticality and could open the way to novel perspectives when the conserved magnetization is interpreted in terms of a homeostatic principle imposed to neural activity.

Resting state connectivity patterns with near-infrared spectroscopy data of the whole head

Resting state cerebral dynamics has been a useful approach to explore the brain's functional organization. In this study, we employed graph theory to deeply investigate resting state functional connectivity (rs-FC) as measured by near-infrared spectroscopy (NIRS). Our results suggest that network parameters are very similar across time and subjects. We also identified the most frequent connections between brain regions and the main hubs that participate in the spontaneous activity of brain hemodynamics. Similar to previous findings, we verified that symmetrically located brain areas are highly connected. Overall, our results introduce new insights in NIRS-based functional connectivity at rest.

From phenotype to genotype in complex brain networks

Generative models are a popular instrument for illuminating the relationships between the hidden variables driving the growth of a complex network and its final topological characteristics, a process known as the "genotype to phenotype problem". However, the definition of a complete methodology encompassing all stages of the analysis, and in particular the validation of the final model, is still an open problem. We here discuss a framework that allows to quantitatively optimise and validate each step of the model creation process. It is based on the execution of a classification task, and on estimating the additional precision provided by the modelled genotype. This encompasses the three main steps of the model creation, namely the selection of topological features, the optimisation of the parameters of the generative model, and the validation of the obtained results. We provide a minimum requirement for a generative model to be useful, prescribing the function mapping genotype to phenotype to be non-monotonic; and we further show how a previously published model does not fulfil such condition, casting doubts on its fitness for the study of neurological disorders. The generality of such framework guarantees its applicability beyond neuroscience, like the emergence of social or technological networks.