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

Reliability of energy landscape analysis of resting-state functional MRI data

Energy landscape analysis is a data-driven method to analyse multidimensional time series, including functional magnetic resonance imaging (fMRI) data. It has been shown to be a useful characterization of fMRI data in health and disease. It fits an Ising model to the data and captures the dynamics of the data as movement of a noisy ball constrained on the energy landscape derived from the estimated Ising model. In the present study, we examine test-retest reliability of the energy landscape analysis. To this end, we construct a permutation test that assesses whether or not indices characterizing the energy landscape are more consistent across different sets of scanning sessions from the same participant (i.e. within-participant reliability) than across different sets of sessions from different participants (i.e. between-participant reliability). We show that the energy landscape analysis has significantly higher within-participant than between-participant test-retest reliability with respect to four commonly used indices. We also show that a variational Bayesian method, which enables us to estimate energy landscapes tailored to each participant, displays comparable test-retest reliability to that using the conventional likelihood maximization method. The proposed methodology paves the way to perform individual-level energy landscape analysis for given data sets with a statistically controlled reliability.

Low-dimensional criticality embedded in high-dimensional awake brain dynamics

Whether cortical neurons operate in a strongly or weakly correlated dynamical regime determines fundamental information processing capabilities and has fueled decades of debate. We offer a resolution of this debate; we show that two important dynamical regimes, typically considered incompatible, can coexist in the same local cortical circuit by separating them into two different subspaces. In awake mouse motor cortex, we find a low-dimensional subspace with large fluctuations consistent with criticality-a dynamical regime with moderate correlations and multi-scale information capacity and transmission. Orthogonal to this critical subspace, we find a high-dimensional subspace containing a desynchronized dynamical regime, which may optimize input discrimination. The critical subspace is apparent only at long timescales, which explains discrepancies among some previous studies. Using a computational model, we show that the emergence of a low-dimensional critical subspace at large timescales agrees with established theory of critical dynamics. Our results suggest that the cortex leverages its high dimensionality to multiplex dynamical regimes across different subspaces.

Energy in functional brain states correlates with cognition in adolescent-onset schizophrenia and healthy persons

Adolescent-onset schizophrenia (AOS) is a relatively rare and under-studied form of schizophrenia with more severe cognitive impairments and poorer outcome compared to adult-onset schizophrenia. Several neuroimaging studies have reported alterations in regional activations that account for activity in individual regions (first-order model) and functional connectivity that reveals pairwise co-activations (second-order model) in AOS compared to controls. The pairwise maximum entropy model, also called the Ising model, can integrate both first-order and second-order terms to elucidate a comprehensive picture of neural dynamics and captures both individual and pairwise activity measures into a single quantity known as energy, which is inversely related to the probability of state occurrence. We applied the MEM framework to task functional MRI data collected on 23 AOS individuals in comparison with 53 healthy control subjects while performing the Penn Conditional Exclusion Test (PCET), which measures executive function that has been repeatedly shown to be more impaired in AOS compared to adult-onset schizophrenia. Accuracy of PCET performance was significantly reduced among AOS compared to controls as expected. Average cumulative energy achieved for a participant over the course of the fMRI negatively correlated with task performance, and the association was stronger than any first-order associations. The AOS subjects spent more time in higher energy states that represent lower probability of occurrence and were associated with impaired executive function and greater severity of psychopathology suggesting that the neural dynamics may be less efficient compared to controls who spent more time in lower energy states occurring with higher probability and hence are more stable and efficient. The energy landscapes in both conditions featured attractors that corresponded to two distinct subnetworks, namely fronto-temporal and parieto-motor. Attractor basins were larger in the controls than in AOS; moreover, fronto-temporal basin size was significantly correlated with cognitive performance in controls but not among the AOS. The single trial trajectories for the AOS group also showed higher variability in concordance with shallow attractor basins among AOS. These findings suggest that the neural dynamics of AOS features more frequent occurrence of less probable states with shallower attractors, which lack the relation to executive function associated with attractors in control subjects suggesting a diminished capacity of AOS to generate task-effective brain states.

Interpretable prediction of brain activity during conversations from multimodal behavioral signals

We present an analytical framework aimed at predicting the local brain activity in uncontrolled experimental conditions based on multimodal recordings of participants' behavior, and its application to a corpus of participants having conversations with another human or a conversational humanoid robot. The framework consists in extracting high-level features from the raw behavioral recordings and applying a dynamic prediction of binarized fMRI-recorded local brain activity using these behavioral features. The objective is to identify behavioral features required for this prediction, and their relative weights, depending on the brain area under investigation and the experimental condition. In order to validate our framework, we use a corpus of uncontrolled conversations of participants with a human or a robotic agent, focusing on brain regions involved in speech processing, and more generally in social interactions. The framework not only predicts local brain activity significantly better than random, it also quantifies the weights of behavioral features required for this prediction, depending on the brain area under investigation and on the nature of the conversational partner. In the left Superior Temporal Sulcus, perceived speech is the most important behavioral feature for predicting brain activity, regardless of the agent, while several features, which differ between the human and robot interlocutors, contribute to the prediction in regions involved in social cognition, such as the TemporoParietal Junction. This framework therefore allows us to study how multiple behavioral signals from different modalities are integrated in individual brain regions during complex social interactions.

Neural Geometrodynamics, Complexity, and Plasticity: A Psychedelics Perspective

We explore the intersection of neural dynamics and the effects of psychedelics in light of distinct timescales in a framework integrating concepts from dynamics, complexity, and plasticity. We call this framework neural geometrodynamics for its parallels with general relativity's description of the interplay of spacetime and matter. The geometry of trajectories within the dynamical landscape of "fast time" dynamics are shaped by the structure of a differential equation and its connectivity parameters, which themselves evolve over "slow time" driven by state-dependent and state-independent plasticity mechanisms. Finally, the adjustment of plasticity processes (metaplasticity) takes place in an "ultraslow" time scale. Psychedelics flatten the neural landscape, leading to heightened entropy and complexity of neural dynamics, as observed in neuroimaging and modeling studies linking increases in complexity with a disruption of functional integration. We highlight the relationship between criticality, the complexity of fast neural dynamics, and synaptic plasticity. Pathological, rigid, or "canalized" neural dynamics result in an ultrastable confined repertoire, allowing slower plastic changes to consolidate them further. However, under the influence of psychedelics, the destabilizing emergence of complex dynamics leads to a more fluid and adaptable neural state in a process that is amplified by the plasticity-enhancing effects of psychedelics. This shift manifests as an acute systemic increase of disorder and a possibly longer-lasting increase in complexity affecting both short-term dynamics and long-term plastic processes. Our framework offers a holistic perspective on the acute effects of these substances and their potential long-term impacts on neural structure and function.

Enhanced simulations of whole-brain dynamics using hybrid resting-state structural connectomes

The human brain, composed of billions of neurons and synaptic connections, is an intricate network coordinating a sophisticated balance of excitatory and inhibitory activities between brain regions. The dynamical balance between excitation and inhibition is vital for adjusting neural input/output relationships in cortical networks and regulating the dynamic range of their responses to stimuli. To infer this balance using connectomics, we recently introduced a computational framework based on the Ising model, which was first developed to explain phase transitions in ferromagnets, and proposed a novel hybrid resting-state structural connectome (rsSC). Here, we show that a generative model based on the Kuramoto phase oscillator can be used to simulate static and dynamic functional connectomes (FC) with rsSC as the coupling weight coefficients, such that the simulated FC aligns well with the observed FC when compared with that simulated traditional structural connectome.

Critical dynamics arise during structured information presentation within embodied in vitro neuronal networks

Understanding how brains process information is an incredibly difficult task. Amongst the metrics characterising information processing in the brain, observations of dynamic near-critical states have generated significant interest. However, theoretical and experimental limitations associated with human and animal models have precluded a definite answer about when and why neural criticality arises with links from attention, to cognition, and even to consciousness. To explore this topic, we used an in vitro neural network of cortical neurons that was trained to play a simplified game of 'Pong' to demonstrate Synthetic Biological Intelligence (SBI). We demonstrate that critical dynamics emerge when neural networks receive task-related structured sensory input, reorganizing the system to a near-critical state. Additionally, better task performance correlated with proximity to critical dynamics. However, criticality alone is insufficient for a neuronal network to demonstrate learning in the absence of additional information regarding the consequences of previous actions. These findings offer compelling support that neural criticality arises as a base feature of incoming structured information processing without the need for higher order cognition.

Low dimensional criticality embedded in high dimensional awake brain dynamics

Whether cortical neurons operate in a strongly or weakly correlated dynamical regime determines fundamental information processing capabilities and has fueled decades of debate. Here we offer a resolution of this debate; we show that two important dynamical regimes, typically considered incompatible, can coexist in the same local cortical circuit by separating them into two different subspaces. In awake mouse motor cortex, we find a low-dimensional subspace with large fluctuations consistent with criticality - a dynamical regime with moderate correlations and multi-scale information capacity and transmission. Orthogonal to this critical subspace, we find a high-dimensional subspace containing a desynchronized dynamical regime, which may optimize input discrimination. The critical subspace is apparent only at long timescales, which explains discrepancies among some previous studies. Using a computational model, we show that the emergence of a low-dimensional critical subspace at large timescale agrees with established theory of critical dynamics. Our results suggest that cortex leverages its high dimensionality to multiplex dynamical regimes across different subspaces.

Biases in inverse Ising estimates of near-critical behavior

Inverse Ising inference allows pairwise interactions of complex binary systems to be reconstructed from empirical correlations. Typical estimators used for this inference, such as pseudo-likelihood maximization (PLM), are biased. Using the Sherrington-Kirkpatrick model as a benchmark, we show that these biases are large in critical regimes close to phase boundaries, and they may alter the qualitative interpretation of the inferred model. In particular, we show that the small-sample bias causes models inferred through PLM to appear closer to criticality than one would expect from the data. Data-driven methods to correct this bias are explored and applied to a functional magnetic resonance imaging data set from neuroscience. Our results indicate that additional care should be taken when attributing criticality to real-world data sets.

Anxiety and depression: A top-down, bottom-up model of circuit function

A functional interplay of bottom-up and top-down processing allows an individual to appropriately respond to the dynamic environment around them. These processing modalities can be represented as attractor states using a dynamical systems model of the brain. The transition probability to move from one attractor state to another is dependent on the stability, depth, neuromodulatory tone, and tonic changes in plasticity. However, how does the relationship between these states change in disease states, such as anxiety or depression? We describe bottom-up and top-down processing from Marr's computational-algorithmic-implementation perspective to understand depressive and anxious disease states. We illustrate examples of bottom-up processing as basolateral amygdala signaling and projections and top-down processing as medial prefrontal cortex internal signaling and projections. Understanding these internal processing dynamics can help us better model the multifaceted elements of anxiety and depression.

Distance to Criticality Undergoes Critical Transition Before Epileptic Seizure Attacks

Epilepsy is a common neurological disorder characterized by recurring seizures, but its underlying mechanisms remain poorly understood. Despite extensive research, there are still gaps in our knowledge about the relationship between brain dynamics and seizures. In this study, our aim is to address these gaps by proposing a novel approach to assess the role of brain network dynamics in the onset of seizures. Specifically, we investigate the relationship between brain dynamics and seizures by tracking the distance to criticality. Our hypothesis is that this distance plays a crucial role in brain state changes and that seizures may be related to critical transitions of this distance. To test this hypothesis, we develop a method to measure the evolution of the brain network's distance to the critical dynamic systems (i.e., the distance to the tipping point, DTP) using dynamic network biomarker theory and random matrix theory. The results show that the DTP of the brain decreases significantly immediately after onset of an epileptic seizure, suggesting that the brain loses its well-defined quasi-critical state during seizures. We refer to this phenomenon as the "criticality of the criticality" (COC). Furthermore, we observe that DTP exhibits a shape transition before and after the onset of the seizures. This phenomenon suggests the possibility of early warning signal (EWS) identification in the dynamic sequence of DTP, which could be utilized for seizure prediction. Our results show that the Hurst exponent, skewness, kurtosis, autocorrelation, and variance of the DTP sequence are potential EWS features. This study advances our understanding of the relationship between brain dynamics and seizures and highlights the potential for using criticality-based measures to predict and prevent seizures.

Critical scaling of whole-brain resting-state dynamics

Scale invariance is a characteristic of neural activity. How this property emerges from neural interactions remains a fundamental question. Here, we studied the relation between scale-invariant brain dynamics and structural connectivity by analyzing human resting-state (rs-) fMRI signals, together with diffusion MRI (dMRI) connectivity and its approximation as an exponentially decaying function of the distance between brain regions. We analyzed the rs-fMRI dynamics using functional connectivity and a recently proposed phenomenological renormalization group (PRG) method that tracks the change of collective activity after successive coarse-graining at different scales. We found that brain dynamics display power-law correlations and power-law scaling as a function of PRG coarse-graining based on functional or structural connectivity. Moreover, we modeled the brain activity using a network of spins interacting through large-scale connectivity and presenting a phase transition between ordered and disordered phases. Within this simple model, we found that the observed scaling features were likely to emerge from critical dynamics and connections exponentially decaying with distance. In conclusion, our study tests the PRG method using large-scale brain activity and theoretical models and suggests that scaling of rs-fMRI activity relates to criticality.

A systematic approach to brain dynamics: cognitive evolution theory of consciousness

The brain integrates volition, cognition, and consciousness seamlessly over three hierarchical (scale-dependent) levels of neural activity for their emergence: a causal or 'hard' level, a computational (unconscious) or 'soft' level, and a phenomenal (conscious) or 'psyche' level respectively. The cognitive evolution theory (CET) is based on three general prerequisites: physicalism, dynamism, and emergentism, which entail five consequences about the nature of consciousness: discreteness, passivity, uniqueness, integrity, and graduation. CET starts from the assumption that brains should have primarily evolved as volitional subsystems of organisms, not as prediction machines. This emphasizes the dynamical nature of consciousness in terms of critical dynamics to account for metastability, avalanches, and self-organized criticality of brain processes, then coupling it with volition and cognition in a framework unified over the levels. Consciousness emerges near critical points, and unfolds as a discrete stream of momentary states, each volitionally driven from oldest subcortical arousal systems. The stream is the brain's way of making a difference via predictive (Bayesian) processing. Its objective observables could be complexity measures reflecting levels of consciousness and its dynamical coherency to reveal how much knowledge (information gain) the brain acquires over the stream. CET also proposes a quantitative classification of both disorders of consciousness and mental disorders within that unified framework.

Time-resolved correlation of distributed brain activity tracks E-I balance and accounts for diverse scale-free phenomena

Much of systems neuroscience posits the functional importance of brain activity patterns that lack natural scales of sizes, durations, or frequencies. The field has developed prominent, and sometimes competing, explanations for the nature of this scale-free activity. Here, we reconcile these explanations across species and modalities. First, we link estimates of excitation-inhibition (E-I) balance with time-resolved correlation of distributed brain activity. Second, we develop an unbiased method for sampling time series constrained by this time-resolved correlation. Third, we use this method to show that estimates of E-I balance account for diverse scale-free phenomena without need to attribute additional function or importance to these phenomena. Collectively, our results simplify existing explanations of scale-free brain activity and provide stringent tests on future theories that seek to transcend these explanations.

Enhanced simulations of whole-brain dynamics using hybrid resting-state structural connectomes

The human brain, composed of billions of neurons and synaptic connections, is an intricate network coordinating a sophisticated balance of excitatory and inhibitory activity between brain regions. The dynamical balance between excitation and inhibition is vital for adjusting neural input/output relationships in cortical networks and regulating the dynamic range of their responses to stimuli. To infer this balance using connectomics, we recently introduced a computational framework based on the Ising model, first developed to explain phase transitions in ferromagnets, and proposed a novel hybrid resting-state structural connectome (rsSC). Here, we show that a generative model based on the Kuramoto phase oscillator can be used to simulate static and dynamic functional connectomes (FC) with rsSC as the coupling weight coefficients, such that the simulated FC well aligns with the observed FC when compared to that simulated with traditional structural connectome. Simulations were performed using the open source framework The Virtual Brain on High Performance Computing infrastructure.

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.

Disrupted Excitation-Inhibition Balance in Cognitively Normal Individuals at Risk of Alzheimer's Disease

BACKGROUND

Sex differences impact Alzheimer's disease (AD) neuropathology, but cell-to-network level dysfunctions in the prodromal phase are unclear. Alterations in hippocampal excitation-inhibition balance (EIB) have recently been linked to early AD pathology.

OBJECTIVE

Examine how AD risk factors (age, APOEɛ4, amyloid-β) relate to hippocampal EIB in cognitively normal males and females using connectome-level measures.

METHODS

Individuals from the OASIS-3 cohort (age 42-95) were studied (N = 437), with a subset aged 65+ undergoing neuropsychological testing (N = 231).

RESULTS

In absence of AD risk factors (APOEɛ4/Aβ+), whole-brain EIB decreases with age more significantly in males than females (p = 0.021, β= -0.007). Regression modeling including APOEɛ4 allele carriers (Aβ-) yielded a significant positive AGE-by-APOE interaction in the right hippocampus for females only (p = 0.013, β= 0.014), persisting with inclusion of Aβ+ individuals (p = 0.012, β= 0.014). Partial correlation analyses of neuropsychological testing showed significant associations with EIB in females: positive correlations between right hippocampal EIB with categorical fluency and whole-brain EIB with the Trail Making Test (p < 0.05).

CONCLUSIONS

Sex differences in EIB emerge during normal aging and progresses differently with AD risk. Results suggest APOEɛ4 disrupts hippocampal balance more than amyloid in females. Increased excitation correlates positively with neuropsychological performance in the female group, suggesting a duality in terms of potential beneficial effects prior to cognitive impairment. This underscores the translational relevance of APOEɛ4 related hyperexcitation in females, potentially informing therapeutic targets or early interventions to mitigate AD progression in this vulnerable population.

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.

Recovery of neural dynamics criticality in personalized whole-brain models of stroke

The critical brain hypothesis states that biological neuronal networks, because of their structural and functional architecture, work near phase transitions for optimal response to internal and external inputs. Criticality thus provides optimal function and behavioral capabilities. We test this hypothesis by examining the influence of brain injury (strokes) on the criticality of neural dynamics estimated at the level of single participants using directly measured individual structural connectomes and whole-brain models. Lesions engender a sub-critical state that recovers over time in parallel with behavior. The improvement of criticality is associated with the re-modeling of specific white-matter connections. We show that personalized whole-brain dynamical models poised at criticality track neural dynamics, alteration post-stroke, and behavior at the level of single participants.

The authors investigate the influence of brain injury (strokes) on the criticality of neural dynamics using directly measured connectomes and whole-brain models. They show that lesions engender a sub-critical state that recovers over time in parallel with behavior.

Inferring excitation-inhibition dynamics using a maximum entropy model unifying brain structure and function

Neural activity coordinated across different scales from neuronal circuits to large-scale brain networks gives rise to complex cognitive functions. Bridging the gap between micro- and macroscale processes, we present a novel framework based on the maximum entropy model to infer a hybrid resting-state structural connectome, representing functional interactions constrained by structural connectivity. We demonstrate that the structurally informed network outperforms the unconstrained model in simulating brain dynamics, wherein by constraining the inference model with the network structure we may improve the estimation of pairwise BOLD signal interactions. Further, we simulate brain network dynamics using Monte Carlo simulations with the new hybrid connectome to probe connectome-level differences in excitation-inhibition balance between apolipoprotein E (APOE)-ε4 carriers and noncarriers. Our results reveal sex differences among APOE-ε4 carriers in functional dynamics at criticality; specifically, female carriers appear to exhibit a lower tolerance to network disruptions resulting from increased excitatory interactions. In sum, the new multimodal network explored here enables analysis of brain dynamics through the integration of structure and function, providing insight into the complex interactions underlying neural activity such as the balance of excitation and inhibition.

Physics Clues on the Mind Substrate and Attributes

The last decade has witnessed a remarkable progress in our understanding of the brain. This has mainly been based on the scrutiny and modeling of the transmission of activity among neurons across lively synapses. A main conclusion, thus far, is that essential features of the mind rely on collective phenomena that emerge from a willful interaction of many neurons that, mediating other cells, form a complex network whose details keep constantly adapting to their activity and surroundings. In parallel, theoretical and computational studies developed to understand many natural and artificial complex systems, which have truthfully explained their amazing emergent features and precise the role of the interaction dynamics and other conditions behind the different collective phenomena they happen to display. Focusing on promising ideas that arise when comparing these neurobiology and physics studies, the present perspective article shortly reviews such fascinating scenarios looking for clues about how high-level cognitive processes such as consciousness, intelligence, and identity can emerge. We, thus, show that basic concepts of physics, such as dynamical phases and non-equilibrium phase transitions, become quite relevant to the brain activity while determined by factors at the subcellular, cellular, and network levels. We also show how these transitions depend on details of the processing mechanism of stimuli in a noisy background and, most important, that one may detect them in familiar electroencephalogram (EEG) recordings. Thus, we associate the existence of such phases, which reveal a brain operating at (non-equilibrium) criticality, with the emergence of most interesting phenomena during memory tasks.

Network Analysis of Time Series: Novel Approaches to Network Neuroscience

In the last two decades, there has been an explosion of interest in modeling the brain as a network, where nodes correspond variously to brain regions or neurons, and edges correspond to structural or statistical dependencies between them. This kind of network construction, which preserves spatial, or structural, information while collapsing across time, has become broadly known as "network neuroscience." In this work, we provide an alternative application of network science to neural data: network-based analysis of non-linear time series and review applications of these methods to neural data. Instead of preserving spatial information and collapsing across time, network analysis of time series does the reverse: it collapses spatial information, instead preserving temporally extended dynamics, typically corresponding to evolution through some kind of phase/state-space. This allows researchers to infer a, possibly low-dimensional, "intrinsic manifold" from empirical brain data. We will discuss three methods of constructing networks from nonlinear time series, and how to interpret them in the context of neural data: recurrence networks, visibility networks, and ordinal partition networks. By capturing typically continuous, non-linear dynamics in the form of discrete networks, we show how techniques from network science, non-linear dynamics, and information theory can extract meaningful information distinct from what is normally accessible in standard network neuroscience approaches.

Avalanche criticality in individuals, fluid intelligence, and working memory

The critical brain hypothesis suggests that efficient neural computation can be achieved through critical brain dynamics. However, the relationship between human cognitive performance and scale‐free brain dynamics remains unclear. In this study, we investigated the whole‐brain avalanche activity and its individual variability in the human resting‐state functional magnetic resonance imaging (fMRI) data. We showed that though the group‐level analysis was inaccurate because of individual variability, the subject wise scale‐free avalanche activity was significantly associated with maximal synchronization entropy of their brain activity. Meanwhile, the complexity of functional connectivity, as well as structure–function coupling, is maximized in subjects with maximal synchronization entropy. We also observed order–disorder phase transitions in resting‐state brain dynamics and found that there were longer times spent in the subcritical regime. These results imply that large‐scale brain dynamics favor the slightly subcritical regime of phase transition. Finally, we showed evidence that the neural dynamics of human participants with higher fluid intelligence and working memory scores are closer to criticality. We identified brain regions whose critical dynamics showed significant positive correlations with fluid intelligence performance and found that these regions were located in the prefrontal cortex and inferior parietal cortex, which were believed to be important nodes of brain networks underlying human intelligence. Our results reveal the possible role that avalanche criticality plays in cognitive performance and provide a simple method to identify the critical point and map cortical states on a spectrum of neural dynamics, ranging from subcriticality to supercriticality.

The relationship between recent PTSD secondary to sexual assault, hippocampal volume and resting state functional connectivity in adolescent girls

Objective

Improved understanding of the time course of neural changes associated with adolescent PTSD would elucidate the development of the disorder and could inform approaches to treatment. We compared hippocampal volumes and resting state functional connectivity (RSFC) in adolescent girls with post-traumatic stress disorder (PTSD) secondary to sexual assault, within six months of onset and age- and gender-matched, non-trauma exposed healthy controls (HCs) in São Paulo, Brazil. We also examined the relationship between pre- and post-treatment PTSD symptoms and RSFC.

Method

We collected brain structure, RSFC, and PTSD symptoms in 30 adolescents with PTSD (mean age: 15.7 ± 1.04 years) and 21 HCs (mean age: 16.2 ± 1.21 years) at baseline. We collected repeated measures in 21 participants with PTSD following treatment; 9 participants dropped out. Hippocampal volume and RSFC from hippocampal and default mode network (DMN) seeds were compared between participants with PTSD and HCs. We examined associations between within-subject changes in RSFC and PTSD symptoms following treatment.

Results

No hippocampal volumetric differences between groups were found. Compared to HCs, adolescents with recent PTSD had reduced RSFC between hippocampus and the lateral parietal node of the DMN, encompassing the angular gyrus, peak coordinates: −38, −54, 16; 116 voxels; peak F_1,47 = 31.76; FDR corrected p = 0.038. Improvements in PTSD symptoms were associated with increased RSFC between hippocampus and part of the lateral parietal node of the DMN, peak coordinates: −38, −84, 38; 316 voxels; peak F_1,47 = 40.28; FDR corrected p < 0.001.

Conclusion

Adolescents with recent PTSD had reduced hippocampal-DMN RSFC, while no group differences in hippocampal volume were found, suggesting that hippocampal function, but not structure, is altered early in the course of PSTD. Following treatment, hippocampal-DMN RSFC increased with symptom improvement and may indicate an important neural mechanism related to successful PTSD treatment.

Empirical Bayes estimation of pairwise maximum entropy model for nonlinear brain state dynamics

The pairwise maximum entropy model (pMEM) has recently gained widespread attention to exploring the nonlinear characteristics of brain state dynamics observed in resting-state functional magnetic resonance imaging (rsfMRI). Despite its unique advantageous features, the practical application of pMEM for individuals is limited as it requires a much larger sample than conventional rsfMRI scans. Thus, this study proposes an empirical Bayes estimation of individual pMEM using the variational expectation-maximization algorithm (VEM-MEM). The performance of the VEM-MEM is evaluated for several simulation setups with various sample sizes and network sizes. Unlike conventional maximum likelihood estimation procedures, the VEM-MEM can reliably estimate the individual model parameters, even with small samples, by effectively incorporating the group information as the prior. As a test case, the individual rsfMRI of children with attention deficit hyperactivity disorder (ADHD) is analyzed compared to that of typically developed children using the default mode network, executive control network, and salient network, obtained from the Healthy Brain Network database. We found that the nonlinear dynamic properties uniquely established on the pMEM differ for each group. Furthermore, pMEM parameters are more sensitive to group differences and are better associated with the behavior scores of ADHD compared to the Pearson correlation-based functional connectivity. The simulation and experimental results suggest that the proposed method can reliably estimate the individual pMEM and characterize the dynamic properties of individuals by utilizing empirical information of the group brain state dynamics.

Dynamics of Long-Range Temporal Correlations in Broadband EEG During Different Motor Execution and Imagery Tasks

Brain activity is composed of oscillatory and broadband arrhythmic components; however, there is more focus on oscillatory sensorimotor rhythms to study movement, but temporal dynamics of broadband arrhythmic electroencephalography (EEG) remain unexplored. We have previously demonstrated that broadband arrhythmic EEG contains both short- and long-range temporal correlations that change significantly during movement. In this study, we build upon our previous work to gain a deeper understanding of these changes in the long-range temporal correlation (LRTC) in broadband EEG and contrast them with the well-known LRTC in alpha oscillation amplitude typically found in the literature. We investigate and validate changes in LRTCs during five different types of movements and motor imagery tasks using two independent EEG datasets recorded with two different paradigms-our finger tapping dataset with single self-initiated asynchronous finger taps and publicly available EEG dataset containing cued continuous movement and motor imagery of fists and feet. We quantified instantaneous changes in broadband LRTCs by detrended fluctuation analysis on single trial 2 s EEG sliding windows. The broadband LRTC increased significantly (p < 0.05) during all motor tasks as compared to the resting state. In contrast, the alpha oscillation LRTC, which had to be computed on longer stitched EEG segments, decreased significantly (p < 0.05) consistently with the literature. This suggests the complementarity of underlying fast and slow neuronal scale-free dynamics during movement and motor imagery. The single trial broadband LRTC gave high average binary classification accuracy in the range of 70.54±10.03% to 76.07±6.40% for all motor execution and imagery tasks and hence can be used in brain-computer interface (BCI). Thus, we demonstrate generalizability, robustness, and reproducibility of novel motor neural correlate, the single trial broadband LRTC, during different motor execution and imagery tasks in single asynchronous and cued continuous motor-BCI paradigms and its contrasting behavior with LRTC in alpha oscillation amplitude.

Bayesian estimation of maximum entropy model for individualized energy landscape analysis of brain state dynamics

The pairwise maximum entropy model (MEM) for resting state functional MRI (rsfMRI) has been used to generate energy landscape of brain states and to explore nonlinear brain state dynamics. Researches using MEM, however, has mostly been restricted to fixed‐effect group‐level analyses, using concatenated time series across individuals, due to the need for large samples in the parameter estimation of MEM. To mitigate the small sample problem in analyzing energy landscapes for individuals, we propose a Bayesian estimation of individual MEM using variational Bayes approximation (BMEM). We evaluated the performances of BMEM with respect to sample sizes and prior information using simulation. BMEM showed advantages over conventional maximum likelihood estimation in reliably estimating model parameters for individuals with small sample data, particularly utilizing the empirical priors derived from group data. We then analyzed individual rsfMRI of the Human Connectome Project to show the usefulness of MEM in differentiating individuals and in exploring neural correlates for human behavior. MEM and its energy landscape properties showed high subject specificity comparable to that of functional connectivity. Canonical correlation analysis identified canonical variables for MEM highly associated with cognitive scores. Inter‐individual variations of cognitive scores were also reflected in energy landscape properties such as energies, occupation times, and basin sizes at local minima. We conclude that BMEM provides an efficient method to characterize dynamic properties of individuals using energy landscape analysis of individual brain states.

Emulative, coherent, and causal dynamics between large-scale brain networks are neurobiomarkers of Accelerated Cognitive Ageing in epilepsy

Accelerated cognitive ageing (ACA) is an ageing co-morbidity in epilepsy that is diagnosed through the observation of an evident IQ decline of more than 1 standard deviation (15 points) around the age of 50 years old. To understand the mechanism of action of this pathology, we assessed brain dynamics with the use of resting-state fMRI data. In this paper, we present novel and promising methods to extract brain dynamics between large-scale resting-state networks: the emulative power, wavelet coherence, and granger causality between the networks were extracted in two resting-state sessions of 24 participants (10 ACA, 14 controls). We also calculated the widely used static functional connectivity to compare the methods. To find the best biomarkers of ACA, and have a better understanding of this epilepsy co-morbidity we compared the aforementioned between-network neurodynamics using classifiers and known machine learning algorithms; and assessed their performance. Results show that features based on the evolutionary game theory on networks approach, the emulative powers, are the best descriptors of the co-morbidity, using dynamics associated with the default mode and dorsal attention networks. With these dynamic markers, linear discriminant analysis could identify ACA patients at 82.9% accuracy. Using wavelet coherence features with decision-tree algorithm, and static functional connectivity features with support vector machine, ACA could be identified at 77.1% and 77.9% accuracy respectively. Granger causality fell short of being a relevant biomarker with best classifiers having an average accuracy of 67.9%. Combining the features based on the game theory, wavelet coherence, Granger-causality, and static functional connectivity- approaches increased the classification performance up to 90.0% average accuracy using support vector machine with a peak accuracy of 95.8%. The dynamics of the networks that lead to the best classifier performances are known to be challenged in elderly. Since our groups were age-matched, the results are in line with the idea of ACA patients having an accelerated cognitive decline. This classification pipeline is promising and could help to diagnose other neuropsychiatric disorders, and contribute to the field of psychoradiology.

The 2-D Cluster Variation Method: Topography Illustrations and Their Enthalpy Parameter Correlations

One of the biggest challenges in characterizing 2-D image topographies is finding a low-dimensional parameter set that can succinctly describe, not so much image patterns themselves, but the nature of these patterns. The 2-D cluster variation method (CVM), introduced by Kikuchi in 1951, can characterize very local image pattern distributions using configuration variables, identifying nearest-neighbor, next-nearest-neighbor, and triplet configurations. Using the 2-D CVM, we can characterize 2-D topographies using just two parameters; the activation enthalpy (ε0) and the interaction enthalpy (ε1). Two different initial topographies (“scale-free-like” and “extreme rich club-like”) were each computationally brought to a CVM free energy minimum, for the case where the activation enthalpy was zero and different values were used for the interaction enthalpy. The results are: (1) the computational configuration variable results differ significantly from the analytically-predicted values well before ε1 approaches the known divergence as ε1→0.881, (2) the range of potentially useful parameter values, favoring clustering of like-with-like units, is limited to the region where ε0<3 and ε1<0.25, and (3) the topographies in the systems that are brought to a free energy minimum show interesting visual features, such as extended “spider legs” connecting previously unconnected “islands,” and as well as evolution of “peninsulas” in what were previously solid masses.

Functional advantages of Lévy walks emerging near a critical point

A special class of random walks, so-called Lévy walks, has been observed in a variety of organisms ranging from cells, insects, fishes, and birds to mammals, including humans. Although their prevalence is considered to be a consequence of natural selection for higher search efficiency, some findings suggest that Lévy walks might also be epiphenomena that arise from interactions with the environment. Therefore, why they are common in biological movements remains an open question. Based on some evidence that Lévy walks are spontaneously generated in the brain and the fact that power-law distributions in Lévy walks can emerge at a critical point, we hypothesized that the advantages of Lévy walks might be enhanced by criticality. However, the functional advantages of Lévy walks are poorly understood. Here, we modeled nonlinear systems for the generation of locomotion and showed that Lévy walks emerging near a critical point had optimal dynamic ranges for coding information. This discovery suggested that Lévy walks could change movement trajectories based on the magnitude of environmental stimuli. We then showed that the high flexibility of Lévy walks enabled switching exploitation/exploration based on the nature of external cues. Finally, we analyzed the movement trajectories of freely moving Drosophila larvae and showed empirically that the Lévy walks may emerge near a critical point and have large dynamic range and high flexibility. Our results suggest that the commonly observed Lévy walks emerge near a critical point and could be explained on the basis of these functional advantages.

Recurrence quantification analysis of dynamic brain networks

Evidence suggests that brain network dynamics are a key determinant of brain function and dysfunction. Here we propose a new framework to assess the dynamics of brain networks based on recurrence analysis. Our framework uses recurrence plots and recurrence quantification analysis to characterize dynamic networks. For resting-state magnetoencephalographic dynamic functional networks (dFNs), we have found that functional networks recur more quickly in people with epilepsy than in healthy controls. This suggests that recurrence of dFNs may be used as a biomarker of epilepsy. For stereo electroencephalography data, we have found that dFNs involved in epileptic seizures emerge before seizure onset, and recurrence analysis allows us to detect seizures. We further observe distinct dFNs before and after seizures, which may inform neurostimulation strategies to prevent seizures. Our framework can also be used for understanding dFNs in healthy brain function and in other neurological disorders besides epilepsy.

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.