Criticality in large-scale brain FMRI dynamics unveiled by a novel point process analysis

A novel method for sparse dynamic functional connectivity analysis from resting-state fMRI

BACKGROUND

There is growing interest in understanding the dynamic functional connectivity (DFC) between distributed brain regions. However, it remains challenging to reliably estimate the temporal dynamics from resting-state functional magnetic resonance imaging (rs-fMRI) due to the limitations of current methods.

NEW METHODS

We propose a new model called HDP-HSMM-BPCA for sparse DFC analysis of high-dimensional rs-fMRI data, which is a temporal extension of probabilistic principal component analysis using Bayesian nonparametric hidden semi-Markov model (HSMM). Specifically, we utilize a hierarchical Dirichlet process (HDP) prior to remove the parametric assumption of the HMM framework, overcoming the limitations of the standard HMM. An attractive superiority is its ability to automatically infer the state-specific latent space dimensionality within the Bayesian formulation.

RESULTS

The experiment results of synthetic data show that our model outperforms the competitive models with relatively higher estimation accuracy. In addition, the proposed framework is applied to real rs-fMRI data to explore sparse DFC patterns. The findings indicate that there is a time-varying underlying structure and sparse DFC patterns in high-dimensional rs-fMRI data.

COMPARISON WITH EXISTING METHODS

Compared with the existing DFC approaches based on HMM, our method overcomes the limitations of standard HMM. The observation model of HDP-HSMM-BPCA can discover the underlying temporal structure of rs-fMRI data. Furthermore, the relevant sparse DFC construction algorithm provides a scheme for estimating sparse DFC.

CONCLUSION

We describe a new computational framework for sparse DFC analysis to discover the underlying temporal structure of rs-fMRI data, which will facilitate the study of brain functional connectivity.

Stationary correlation pattern in highly non-stationary MEG recordings of healthy subjects and its relation to former EEG studies

In this study, we analyze magnetoencephalographic (MEG) recordings from 48 clinically healthy subjects obtained from the Human Connectome Project (HCP) while they performed a working memory task and a motor task. Our results reveal a well-developed, stable interrelation pattern that spans the entire scalp and is nearly universal, being almost task- and subject-independent. Additionally, we demonstrate that this pattern closely resembles a stationary correlation pattern (SCP) observed in EEG signals under various physiological and pathological conditions (the distribution of Pearson correlations are centered at about 0.75). Furthermore, we identify the most effective EEG reference for studying the brain's functional network derived from lag-zero cross-correlations. We contextualize these findings within the theory of complex dynamical systems operating near a critical point of a phase transition.

Subject-based assessment of large-scale integration dynamics in epileptic brain networks: insights from the intrinsic ignition framework

This study examined the dynamic properties of brain regions involved in the genesis and spread of seizures in 10 individuals diagnosed with pharmacoresistant focal epilepsy. The patients and 30 healthy controls underwent resting-state functional magnetic resonance imaging scans and the brain's functional network dynamics were analyzed using the intrinsic ignition framework. Comparative statistical analyses examined the differences in the integration and metastability measures in both groups in the whole brain and specific local brain regions. Invasive electroencephalography evaluations validated the findings of significant global and regional changes in the patient's brain network dynamics. There was a marked increase in global integration and metastability across the brain, reflecting substantial alterations in the overall connectivity and flexibility of the functional networks. Specific brain regions exhibited paradoxical dynamics within the seizure onset zone, with decreased intrinsic ignition and increased metastability. Increased intrinsic ignition was observed in remote brain regions, suggesting a reorganization of the brain network hubs and potential pathways for seizure propagation. Using the intrinsic ignition framework provided insights into dynamic alterations in the brain networks of patients with epilepsy. These have increased our understanding of the mechanisms underlying epileptic seizures and may guide the development of diagnostic biomarkers and targeted therapeutic interventions.

Criticality and universality in neuronal cultures during "up" and "down" states

The brain can be seen as a self-organized dynamical system that optimizes information processing and storage capabilities. This is supported by studies across scales, from small neuronal assemblies to the whole brain, where neuronal activity exhibits features typically associated with phase transitions in statistical physics. Such a critical state is characterized by the emergence of scale-free statistics as captured, for example, by the sizes and durations of activity avalanches corresponding to a cascading process of information flow. Another phenomenon observed during sleep, under anesthesia, and in in vitro cultures, is that cortical and hippocampal neuronal networks alternate between "up" and "down" states characterized by very distinct firing rates. Previous theoretical work has been able to relate these two concepts and proposed that only up states are critical whereas down states are subcritical, also indicating that the brain spontaneously transitions between the two. Using high-speed high-resolution calcium imaging recordings of neuronal cultures, we test this hypothesis here by analyzing the neuronal avalanche statistics in populations of thousands of neurons during "up" and "down" states separately. We find that both "up" and "down" states can exhibit scale-free behavior when taking into account their intrinsic time scales. In particular, the statistical signature of "down" states is indistinguishable from those observed previously in cultures without "up" states. We show that such behavior can not be explained by network models of non-conservative leaky integrate-and-fire neurons with short-term synaptic depression, even when realistic noise levels, spatial network embeddings, and heterogeneous populations are taken into account, which instead exhibits behavior consistent with previous theoretical models. Similar differences were also observed when taking into consideration finite-size scaling effects, suggesting that the intrinsic dynamics and self-organization mechanisms of these cultures might be more complex than previously thought. In particular, our findings point to the existence of different mechanisms of neuronal communication, with different time scales, acting during either high-activity or low-activity states, potentially requiring different plasticity mechanisms.

Stable individual differences from dynamic patterns of function: brain network flexibility predicts openness/intellect, intelligence, and psychoticism

A growing understanding of the nature of brain function has led to increased interest in interpreting the properties of large-scale brain networks. Methodological advances in network neuroscience provide means to decompose these networks into smaller functional communities and measure how they reconfigure over time as an index of their dynamic and flexible properties. Recent evidence has identified associations between flexibility and a variety of traits pertaining to complex cognition including creativity and working memory. The present study used measures of dynamic resting-state functional connectivity in data from the Human Connectome Project (n = 994) to test associations with Openness/Intellect, general intelligence, and psychoticism, three traits that involve flexible cognition. Using a machine-learning cross-validation approach, we identified reliable associations of intelligence with cohesive flexibility of parcels in large communities across the cortex, of psychoticism with disjoint flexibility, and of Openness/Intellect with overall flexibility among parcels in smaller communities. These findings are reasonably consistent with previous theories of the neural correlates of these traits and help to expand on previous associations of behavior with dynamic functional connectivity, in the context of broad personality dimensions.

Human brain state dynamics are highly reproducible and associated with neural and behavioral features

Neural activity and behavior vary within an individual (states) and between individuals (traits). However, the mapping of state-trait neural variation to behavior is not well understood. To address this gap, we quantify moment-to-moment changes in brain-wide co-activation patterns derived from resting-state functional magnetic resonance imaging. In healthy young adults, we identify reproducible spatiotemporal features of co-activation patterns at the single-subject level. We demonstrate that a joint analysis of state-trait neural variations and feature reduction reveal general motifs of individual differences, encompassing state-specific and general neural features that exhibit day-to-day variability. The principal neural variations co-vary with the principal variations of behavioral phenotypes, highlighting cognitive function, emotion regulation, alcohol and substance use. Person-specific probability of occupying a particular co-activation pattern is reproducible and associated with neural and behavioral features. This combined analysis of state-trait variations holds promise for developing reproducible neuroimaging markers of individual life functional outcome.

Role of the Locus Coeruleus Arousal Promoting Neurons in Maintaining Brain Criticality across the Sleep-Wake Cycle

Sleep control depends on a delicate interplay among brain regions. This generates a complex temporal architecture with numerous sleep-stage transitions and intermittent fluctuations to micro-states and brief arousals. These temporal dynamics exhibit hallmarks of criticality, suggesting that tuning to criticality is essential for spontaneous sleep-stage and arousal transitions. However, how the brain maintains criticality remains not understood. Here, we investigate θ- and δ-burst dynamics during the sleep-wake cycle of rats (Sprague-Dawley, adult male) with lesion in the wake-promoting locus coeruleus (LC). We show that, in control rats, θ- and δ-bursts exhibit power-law (θ-bursts, active phase) and exponential-like (δ-bursts, quiescent phase) duration distributions, as well as power-law long-range temporal correlations (LRTCs)-typical of non-equilibrium systems self-organizing at criticality. Furthermore, consecutive θ- and δ-bursts durations are characterized by anti-correlated coupling, indicating a new class of self-organized criticality that emerges from underlying feedback between neuronal populations and brain areas involved in generating arousals and sleep states. In contrast, we uncover that LC lesion leads to alteration of θ- and δ-burst critical features, with change in duration distributions and correlation properties, and increase in θ-δ coupling. Notably, these LC-lesion effects are opposite to those observed for lesions in the sleep-promoting ventrolateral preoptic (VLPO) nucleus. Our findings indicate that critical dynamics of θ- and δ-bursts arise from a balanced interplay of LC and VLPO, which maintains brain tuning to criticality across the sleep-wake cycle-a non-equilibrium behavior in sleep micro-architecture at short timescales that coexists with large-scale sleep-wake homeostasis.

Existence of multiple transitions of the critical state due to anesthetics

Scale-free statistics of coordinated neuronal activity, suggesting a universal operating mechanism across spatio-temporal scales, have been proposed as a necessary condition of healthy resting-state brain activity. Recent studies have focused on anesthetic agents to induce distinct neural states in which consciousness is altered to understand the importance of critical dynamics. However, variation in experimental techniques, species, and anesthetics, have made comparisons across studies difficult. Here we conduct a survey of several common anesthetics (isoflurane, pentobarbital, ketamine) at multiple dosages, using calcium wide-field optical imaging of the mouse cortex. We show that while low-dose anesthesia largely preserves scale-free statistics, surgical plane anesthesia induces multiple dynamical modes, most of which do not maintain critical avalanche dynamics. Our findings indicate multiple pathways away from default critical dynamics associated with quiet wakefulness, not only reflecting differences between these common anesthetics but also showing significant variations in individual responses. This is suggestive of a non-trivial relationship between criticality and the underlying state of the subject.

Imbalanced temporal states of cortical blood-oxygen-level-dependent signal variability during rest in episodic migraine

BACKGROUND

Migraine has been associated with functional brain changes including altered connectivity and activity both during and between headache attacks. Recent studies established that the variability of the blood-oxygen-level-dependent (BOLD) signal is an important attribute of brain activity, which has so far been understudied in migraine. In this study, we investigate how time-varying measures of BOLD variability change interictally in episodic migraine patients.

METHODS

Two independent resting state functional MRI datasets acquired on 3T (discovery cohort) and 1.5T MRI scanners (replication cohort) including 99 episodic migraine patients (n3T = 42, n1.5T=57) and 78 healthy controls (n3T = 46, n1.5T=32) were analyzed in this cross-sectional study. A framework using time-varying measures of BOLD variability was applied to derive BOLD variability states. Descriptors of BOLD variability states such as dwell time and fractional occupancy were calculated, then compared between migraine patients and healthy controls using Mann-Whitney U-tests. Spearman's rank correlation was calculated to test associations with clinical parameters.

RESULTS

Resting-state activity was characterized by states of high and low BOLD signal variability. Migraine patients in the discovery cohort spent more time in the low variability state (mean dwell time: p = 0.014, median dwell time: p = 0.022, maximum dwell time: p = 0.013, fractional occupancy: p = 0.013) and less time in the high variability state (mean dwell time: p = 0.021, median dwell time: p = 0.021, maximum dwell time: p = 0.025, fractional occupancy: p = 0.013). Higher uptime of the low variability state was associated with greater disability as measured by MIDAS scores (maximum dwell time: R = 0.45, p = 0.007; fractional occupancy: R = 0.36, p = 0.035). Similar results were observed in the replication cohort.

CONCLUSION

Episodic migraine patients spend more time in a state of low BOLD variability during rest in headache-free periods, which is associated with greater disability. BOLD variability states show potential as a replicable functional imaging marker in episodic migraine.

Differences in spatiotemporal brain network dynamics of Montessori and traditionally schooled students

Across development, experience has a strong impact on the way we think and adapt. School experience affects academic and social-emotional outcomes, yet whether differences in pedagogical experience modulate underlying brain network development is still unknown. In this study, we compared the brain network dynamics of students with different pedagogical backgrounds. Specifically, we characterized the diversity and stability of brain activity at rest by combining both resting-state fMRI and diffusion-weighted structural imaging data of 87 4-18 years old students experiencing either the Montessori pedagogy (i.e., student-led, trial-and-error pedagogy) or the traditional pedagogy (i.e., teacher-led, test-based pedagogy). Our results revealed spatiotemporal brain dynamics differences between students as a function of schooling experience at the whole-brain level. Students from Montessori schools showed overall higher functional integration (higher system diversity) and neural stability (lower spatiotemporal diversity) compared to traditionally schooled students. Higher integration was explained mainly through the cerebellar (CBL) functional network. In contrast, higher temporal stability was observed in the ventral attention, dorsal attention, somatomotor, frontoparietal, and CBL functional networks. This study suggests a form of experience-dependent dynamic functional connectivity plasticity, in learning-related networks.

Time-varying phase synchronization of resting-state functional magnetic resonance imaging reveals a shift toward self-referential processes during sustained pain

ABSTRACT

Growing evidence has suggested that time-varying functional connectivity between different brain regions might underlie the dynamic experience of pain. This study used a novel, data-driven framework to characterize the dynamic interactions of large-scale brain networks during sustained pain by estimating recurrent patterns of phase-synchronization. Resting-state functional magnetic resonance imaging signals were collected from 50 healthy participants before (once) and after (twice) the onset of sustained pain that was induced by topical application of capsaicin cream. We first decoded the instantaneous phase of neural activity and then applied leading eigenvector dynamic analysis on the time-varying phase-synchronization. We identified 3 recurrent brain states that show distinctive phase-synchronization. The presence of state 1, characterized by phase-synchronization between the default mode network and auditory, visual, and sensorimotor networks, together with transitions towards this brain state, increased during sustained pain. These changes can account for the perceived pain intensity and reported unpleasantness induced by capsaicin application. In contrast, state 3, characterized by phase-synchronization between the cognitive control network and sensory networks, decreased after the onset of sustained pain. These results are indicative of a shift toward internally directed self-referential processes (state 1) and away from externally directed cognitive control processes (state 3) during sustained pain.

Altered spreading of fast aperiodic brain waves relates to disease duration in Amyotrophic Lateral Sclerosis

OBJECTIVE

To test the hypothesis that patients affected by Amyotrophic Lateral Sclerosis (ALS) show an altered spatio-temporal spreading of neuronal avalanches in the brain, and that this may related to the clinical picture.

METHODS

We obtained the source-reconstructed magnetoencephalography (MEG) signals from thirty-six ALS patients and forty-two healthy controls. Then, we used the construct of the avalanche transition matrix (ATM) and the corresponding network parameter nodal strength to quantify the changes in each region, since this parameter provides key information about which brain regions are mostly involved in the spreading avalanches.

RESULTS

ALS patients presented higher values of the nodal strength in both cortical and sub-cortical brain areas. This parameter correlated directly with disease duration.

CONCLUSIONS

In this work, we provide a deeper characterization of neuronal avalanches propagation in ALS, describing their spatio-temporal trajectories and identifying the brain regions most likely to be involved in the process. This makes it possible to recognize the brain areas that take part in the pathogenic mechanisms of ALS. Furthermore, the nodal strength of the involved regions correlates directly with disease duration.

SIGNIFICANCE

Our results corroborate the clinical relevance of aperiodic, fast large-scale brain activity as a biomarker of microscopic changes induced by neurophysiological processes.

Preconfigured architecture of the developing mouse brain

In the adult brain, structural and functional parameters, such as synaptic sizes and neuronal firing rates, follow right-skewed and heavy-tailed distributions. While this organization is thought to have significant implications, its development is still largely unknown. Here, we address this knowledge gap by investigating a large-scale dataset recorded from the prefrontal cortex and the olfactory bulb of mice aged 4-60 postnatal days. We show that firing rates and spike train interactions have a largely stable distribution shape throughout the first 60 postnatal days and that the prefrontal cortex displays a functional small-world architecture. Moreover, early brain activity exhibits an oligarchical organization, where high-firing neurons have hub-like properties. In a neural network model, we show that analogously right-skewed and heavy-tailed synaptic parameters are instrumental to consistently recapitulate the experimental data. Thus, functional and structural parameters in the developing brain are already extremely distributed, suggesting that this organization is preconfigured and not experience dependent.

Estimation of static and dynamic functional connectivity in resting-state fMRI using zero-frequency resonator

Resting-state functional magnetic resonance imaging (rs-fMRI) is increasingly being used to infer the functional organization of the brain. Blood oxygen level-dependent (BOLD) features related to spontaneous neuronal activity, are yet to be clearly understood. Prior studies have hypothesized that rs-fMRI is spontaneous event-related and these events convey crucial information about the neuronal activity in estimating resting state functional connectivity (FC). Attempts have been made to extract these temporal events using a predetermined threshold. However, the thresholding methods in addition to being very sensitive to noise, may consider redundant events or exclude the low-valued inflection points. Here, we extract the event-related temporal onsets from the rs-fMRI time courses using a zero-frequency resonator (ZFR). The ZFR reflects the transient behavior of the BOLD events at its output. The conditional rate (CR) of the BOLD events occurring in a time course with respect to a seed time course is used to derive static FC. The temporal activity around the estimated events called high signal-to-noise ratio (SNR) segments are also obtained in the rs-fMRI time course and are then used to compute static and dynamic FCs during rest. Coactivation pattern (CAP) is the dynamic FC obtained using the high SNR segments driven by the ZFR. The static FC demonstrates that the ZFR-based CR distinguishes the coactivation and non-coactivation scores well in the distribution. CAP analysis demonstrated the stable and longer dwell time dominant resting state functional networks with high SNR segments driven by the ZFR. Static and dynamic FC analysis underpins that the ZFR-driven temporal onsets of BOLD events derive reliable and consistent FCs in the resting brain using a subset of the time points.

Coactivation pattern analysis reveals altered whole-brain functional transient dynamics in autism spectrum disorder

Recent studies on autism spectrum disorder (ASD) have identified recurring states dominated by similar coactivation pattern (CAP) and revealed associations between dysfunction in seed-based large-scale brain networks and clinical symptoms. However, the presence of abnormalities in moment-to-moment whole-brain dynamics in ASD remains uncertain. In this study, we employed seed-free CAP analysis to identify transient brain activity configurations and investigate dynamic abnormalities in ASD. We utilized a substantial multisite resting-state fMRI dataset consisting of 354 individuals with ASD and 446 healthy controls (HCs, from HC groups and 2). CAP were generated from a subgroup of all HC subjects (HC group 1) through temporal K-means clustering, identifying four CAPs. These four CAPs exhibited either the activation or inhibition of the default mode network (DMN) and were grouped into two pairs with opposing spatial CAPs. CAPs for HC group 2 and ASD were identified by their spatial similarity to those for HC group 1. Compared with individuals in HC group 2, those with ASD spent more time in CAPs involving the ventral attention network but less time in CAPs related to executive control and the dorsal attention network. Support vector machine analysis demonstrated that the aberrant dynamic characteristics of CAPs achieved an accuracy of 74.87% in multisite classification. In addition, we used whole-brain dynamics to predict symptom severity in ASD. Our findings revealed whole-brain dynamic functional abnormalities in ASD from a single transient perspective, emphasizing the importance of the DMN in abnormal dynamic functional activity in ASD and suggesting that temporally dynamic techniques offer novel insights into time-varying neural processes.

The flattening of spacetime hierarchy of the N,N-dimethyltryptamine brain state is characterized by harmonic decomposition of spacetime (HADES) framework

The human brain is a complex system, whose activity exhibits flexible and continuous reorganization across space and time. The decomposition of whole-brain recordings into harmonic modes has revealed a repertoire of gradient-like activity patterns associated with distinct brain functions. However, the way these activity patterns are expressed over time with their changes in various brain states remains unclear. Here, we investigate healthy participants taking the serotonergic psychedelic N,N-dimethyltryptamine (DMT) with the Harmonic Decomposition of Spacetime (HADES) framework that can characterize how different harmonic modes defined in space are expressed over time. HADES demonstrates significant decreases in contributions across most low-frequency harmonic modes in the DMT-induced brain state. When normalizing the contributions by condition (DMT and non-DMT), we detect a decrease specifically in the second functional harmonic, which represents the uni- to transmodal functional hierarchy of the brain, supporting the leading hypothesis that functional hierarchy is changed in psychedelics. Moreover, HADES' dynamic spacetime measures of fractional occupancy, life time and latent space provide a precise description of the significant changes of the spacetime hierarchical organization of brain activity in the psychedelic state.

Associations of Brain Entropy Estimated by Resting State fMRI With Physiological Indices, Body Mass Index, and Cognition

BACKGROUND

In recent years, resting-state fMRI (rsfMRI)-based brain entropy (BEN) has gained increasing interest as a tool to characterize brain activity. While previous studies indicate that BEN is correlated with cognition, it remains unclear whether BEN is influenced by other factors that typically affect brain activity measured by fMRI.

PURPOSE

To investigate the relationship between BEN and physiological indices, including respiratory rate (RR), heart rate (HR), systolic blood pressure (s-BP), and body mass index (BMI), and to investigate whether and to what extent the relationship between BEN and cognition is influenced by physiological variables.

STUDY TYPE

Retrospective.

SUBJECTS

One thousand two hundred six healthy subjects (mean age: 28.83 ± 3.69 years; 550 male) with rsfMRI datasets selected from the Human Connectome Project (HCP).

FIELD STRENGTH/SEQUENCE

Multiband echo planar imaging (EPI) sequence at 3.0 Tesla.

ASSESSMENT

Neurocognitive, physical health (RR, HR, s-BP, BMI), and rsfMRI data were retrieved from the HCP datasets. Neurocognition was measured through the total cognition composite (TCC) score provided by HCP. BEN maps were calculated from rsfMRI data.

STATISTICAL TESTS

Multiple regression models, pheight-family wise error (FWE) < 0.05 and pcluster-FWE < 0.05 were considered statistically significant.

RESULTS

BEN was negatively associated with RR (T-thresholds ranging from 4.75 to 4.8; r-threshold = |0.15|) and positively associated with s-BP and BMI (T-thresholds ranging from 4.75 to 4.8; r-threshold = |0.15|) in areas overlapping with the default mode network. After controlling the physiological effects, BEN still showed regional associations with TCC, including negative associations (T-thresholds = 3.09; r-threshold = |0.1|) in the fronto-parietal cortex and positive associations (T-thresholds = 3.09; r-threshold = |0.1|) in the sensorimotor system (motor network and the limbic system).

DATA CONCLUSIONS

RR negatively affects rsfMRI-derived BEN, while s-BP and BMI positively affect BEN. The positive associations between BEN and cognition in the motor network and the limbic system might indicate a facilitation of information processing in the sensorimotor system.

EVIDENCE LEVEL

3 TECHNICAL EFFICACY: Stage 3.

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.

Laying the foundations for a theory of consciousness: the significance of critical brain dynamics for the formation of conscious states

Empirical evidence indicates that conscious states, distinguished by the presence of phenomenal qualities, are closely linked to synchronized neural activity patterns whose dynamical characteristics can be attributed to self-organized criticality and phase transitions. These findings imply that insight into the mechanism by which the brain controls phase transitions will provide a deeper understanding of the fundamental mechanism by which the brain manages to transcend the threshold of consciousness. This article aims to show that the initiation of phase transitions and the formation of synchronized activity patterns is due to the coupling of the brain to the zero-point field (ZPF), which plays a central role in quantum electrodynamics (QED). The ZPF stands for the presence of ubiquitous vacuum fluctuations of the electromagnetic field, represented by a spectrum of normal modes. With reference to QED-based model calculations, the details of the coupling mechanism are revealed, suggesting that critical brain dynamics is governed by the resonant interaction of the ZPF with the most abundant neurotransmitter glutamate. The pyramidal neurons in the cortical microcolumns turn out to be ideally suited to control this interaction. A direct consequence of resonant glutamate-ZPF coupling is the amplification of specific ZPF modes, which leads us to conclude that the ZPF is the key to the understanding of consciousness and that the distinctive feature of neurophysiological processes associated with conscious experience consists in modulating the ZPF. Postulating that the ZPF is an inherently sentient field and assuming that the spectrum of phenomenal qualities is represented by the normal modes of the ZPF, the significance of resonant glutamate-ZPF interaction for the formation of conscious states becomes apparent in that the amplification of specific ZPF modes is inextricably linked with the excitation of specific phenomenal qualities. This theory of consciousness, according to which phenomenal states arise through resonant amplification of zero-point modes, is given the acronym TRAZE. An experimental setup is specified that can be used to test a corollary of the theory, namely, the prediction that normally occurring conscious perceptions are absent under experimental conditions in which resonant glutamate-ZPF coupling is disrupted.

Emergence of complex oscillatory dynamics in the neuronal networks with long activity time of inhibitory synapses

The brain displays complex dynamics, including collective oscillations, and extensive research has been conducted to understand their generation. However, our understanding of how biological constraints influence these oscillations is incomplete. This study investigates the essential properties of neuronal networks needed to generate oscillations resembling those in the brain. A simple discrete-time model of interconnected excitable elements is developed, capable of closely resembling the complex oscillations observed in biological neural networks. In the model, synaptic connections remain active for a duration exceeding individual neuron activity. We show that the inhibitory synapses must exhibit longer activity than excitatory synapses to produce a diverse range of the dynamical states, including biologically plausible oscillations. Upon meeting this condition, the transition between different dynamical states can be controlled by external stochastic input to the neurons. The study provides a comprehensive explanation for the emergence of distinct dynamical states in neural networks based on specific parameters.

Early altered directionality of resting brain network state transitions in the TgF344-AD rat model of Alzheimer's disease

Introduction

Alzheimer's disease (AD) is a progressive neurodegenerative disease resulting in memory loss and cognitive decline. Synaptic dysfunction is an early hallmark of the disease whose effects on whole-brain functional architecture can be identified using resting-state functional MRI (rsfMRI). Insights into mechanisms of early, whole-brain network alterations can help our understanding of the functional impact of AD's pathophysiology.

Methods

Here, we obtained rsfMRI data in the TgF344-AD rat model at the pre- and early-plaque stages. This model recapitulates the major pathological and behavioral hallmarks of AD. We used co-activation pattern (CAP) analysis to investigate if and how the dynamic organization of intrinsic brain functional networks states, undetectable by earlier methods, is altered at these early stages.

Results

We identified and characterized six intrinsic brain states as CAPs, their spatial and temporal features, and the transitions between the different states. At the pre-plaque stage, the TgF344-AD rats showed reduced co-activation of hub regions in the CAPs corresponding to the default mode-like and lateral cortical network. Default mode-like network activity segregated into two distinct brain states, with one state characterized by high co-activation of the basal forebrain. This basal forebrain co-activation was reduced in TgF344-AD animals mainly at the pre-plaque stage. Brain state transition probabilities were altered at the pre-plaque stage between states involving the default mode-like network, lateral cortical network, and basal forebrain regions. Additionally, while the directionality preference in the network-state transitions observed in the wild-type animals at the pre-plaque stage had diminished at the early-plaque stage, TgF344-AD animals continued to show directionality preference at both stages.

Discussion

Our study enhances the understanding of intrinsic brain state dynamics and how they are impacted at the early stages of AD, providing a nuanced characterization of the early, functional impact of the disease's neurodegenerative process.

Whole-Brain Dynamics Disruptions in the Progression of Alzheimer's Disease: Understanding the Influence of Amyloid-Beta and Tau

INTRODUCTION

Alzheimer's disease (AD) affects brain structure and function along its evolution, but brain network dynamic changes remain largely unknown.

METHODS

To understand how AD shapes brain activity, we investigated the spatiotemporal dynamics and resting state functional networks using the intrinsic ignition framework, which characterizes how an area transmits neuronal activity to others, resulting in different degrees of integration. Healthy participants, MCI, and AD patients were scanned using resting state fMRI. Mixed effects models were used to assess the impact of ABeta and tau, at the regional and whole-brain levels.

RESULTS

Dynamic complexity is progressively reduced, with Healthy participants showing higher metastability (i.e., a more complex dynamical regime over time) than observed in the other stages, while AD subjects showed the lowest.

DISCUSSION

Our study provides further insight into how AD modulates brain network dynamics along its evolution, progressively disrupting the whole-brain and resting state network dynamics.

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.

Complexity Measures for EEG Microstate Sequences: Concepts and Algorithms

EEG microstate sequence analysis quantifies properties of ongoing brain electrical activity which is known to exhibit complex dynamics across many time scales. In this report we review recent developments in quantifying microstate sequence complexity, we classify these approaches with regard to different complexity concepts, and we evaluate excess entropy as a yet unexplored quantity in microstate research. We determined the quantities entropy rate, excess entropy, Lempel-Ziv complexity (LZC), and Hurst exponents on Potts model data, a discrete statistical mechanics model with a temperature-controlled phase transition. We then applied the same techniques to EEG microstate sequences from wakefulness and non-REM sleep stages and used first-order Markov surrogate data to determine which time scales contributed to the different complexity measures. We demonstrate that entropy rate and LZC measure the Kolmogorov complexity (randomness) of microstate sequences, whereas excess entropy and Hurst exponents describe statistical complexity which attains its maximum at intermediate levels of randomness. We confirmed the equivalence of entropy rate and LZC when the LZ-76 algorithm is used, a result previously reported for neural spike train analysis (Amigó et al., Neural Comput 16:717-736, https://doi.org/10.1162/089976604322860677 , 2004). Surrogate data analyses prove that entropy-based quantities and LZC focus on short-range temporal correlations, whereas Hurst exponents include short and long time scales. Sleep data analysis reveals that deeper sleep stages are accompanied by a decrease in Kolmogorov complexity and an increase in statistical complexity. Microstate jump sequences, where duplicate states have been removed, show higher randomness, lower statistical complexity, and no long-range correlations. Regarding the practical use of these methods, we suggest that LZC can be used as an efficient entropy rate estimator that avoids the estimation of joint entropies, whereas entropy rate estimation via joint entropies has the advantage of providing excess entropy as the second parameter of the same linear fit. We conclude that metrics of statistical complexity are a useful addition to microstate analysis and address a complexity concept that is not yet covered by existing microstate algorithms while being actively explored in other areas of brain research.

A recurrent Hopfield network for estimating meso-scale effective connectivity in MEG

The architecture of communication within the brain, represented by the human connectome, has gained a paramount role in the neuroscience community. Several features of this communication, e.g., the frequency content, spatial topology, and temporal dynamics are currently well established. However, identifying generative models providing the underlying patterns of inhibition/excitation is very challenging. To address this issue, we present a novel generative model to estimate large-scale effective connectivity from MEG. The dynamic evolution of this model is determined by a recurrent Hopfield neural network with asymmetric connections, and thus denoted Recurrent Hopfield Mass Model (RHoMM). Since RHoMM must be applied to binary neurons, it is suitable for analyzing Band Limited Power (BLP) dynamics following a binarization process. We trained RHoMM to predict the MEG dynamics through a gradient descent minimization and we validated it in two steps. First, we showed a significant agreement between the similarity of the effective connectivity patterns and that of the interregional BLP correlation, demonstrating RHoMM's ability to capture individual variability of BLP dynamics. Second, we showed that the simulated BLP correlation connectomes, obtained from RHoMM evolutions of BLP, preserved some important topological features, e.g, the centrality of the real data, assuring the reliability of RHoMM. Compared to other biophysical models, RHoMM is based on recurrent Hopfield neural networks, thus, it has the advantage of being data-driven, less demanding in terms of hyperparameters and scalable to encompass large-scale system interactions. These features are promising for investigating the dynamics of inhibition/excitation at different spatial scales.

Excitation-Inhibition Balance, Neural Criticality, and Activities in Neuronal Circuits

Neural activities in local circuits exhibit complex and multilevel dynamic features. Individual neurons spike irregularly, which is believed to originate from receiving balanced amounts of excitatory and inhibitory inputs, known as the excitation-inhibition balance. The spatial-temporal cascades of clustered neuronal spikes occur in variable sizes and durations, manifested as neural avalanches with scale-free features. These may be explained by the neural criticality hypothesis, which posits that neural systems operate around the transition between distinct dynamic states. Here, we summarize the experimental evidence for and the underlying theory of excitation-inhibition balance and neural criticality. Furthermore, we review recent studies of excitatory-inhibitory networks with synaptic kinetics as a simple solution to reconcile these two apparently distinct theories in a single circuit model. This provides a more unified understanding of multilevel neural activities in local circuits, from spontaneous to stimulus-response dynamics.

Complexity of STG signals and linguistic rhythm: a methodological study for EEG data

The superior temporal and the Heschl's gyri of the human brain play a fundamental role in speech processing. Neurons synchronize their activity to the amplitude envelope of the speech signal to extract acoustic and linguistic features, a process known as neural tracking/entrainment. Electroencephalography has been extensively used in language-related research due to its high temporal resolution and reduced cost, but it does not allow for a precise source localization. Motivated by the lack of a unified methodology for the interpretation of source reconstructed signals, we propose a method based on modularity and signal complexity. The procedure was tested on data from an experiment in which we investigated the impact of native language on tracking to linguistic rhythms in two groups: English natives and Spanish natives. In the experiment, we found no effect of native language but an effect of language rhythm. Here, we compare source projected signals in the auditory areas of both hemispheres for the different conditions using nonparametric permutation tests, modularity, and a dynamical complexity measure. We found increasing values of complexity for decreased regularity in the stimuli, giving us the possibility to conclude that languages with less complex rhythms are easier to track by the auditory cortex.

Connectome-based reservoir computing with the conn2res toolbox

The connection patterns of neural circuits form a complex network. How signaling in these circuits manifests as complex cognition and adaptive behaviour remains the central question in neuroscience. Concomitant advances in connectomics and artificial intelligence open fundamentally new opportunities to understand how connection patterns shape computational capacity in biological brain networks. Reservoir computing is a versatile paradigm that uses high-dimensional, nonlinear dynamical systems to perform computations and approximate cognitive functions. Here we present conn2res: an open-source Python toolbox for implementing biological neural networks as artificial neural networks. conn2res is modular, allowing arbitrary network architecture and dynamics to be imposed. The toolbox allows researchers to input connectomes reconstructed using multiple techniques, from tract tracing to noninvasive diffusion imaging, and to impose multiple dynamical systems, from spiking neurons to memristive dynamics. The versatility of the conn2res toolbox allows us to ask new questions at the confluence of neuroscience and artificial intelligence. By reconceptualizing function as computation, conn2res sets the stage for a more mechanistic understanding of structure-function relationships in brain networks.

Glioma-Induced Alterations in Excitatory Neurons are Reversed by mTOR Inhibition

Gliomas are highly aggressive brain tumors characterized by poor prognosis and composed of diffusely infiltrating tumor cells that intermingle with non-neoplastic cells in the tumor microenvironment, including neurons. Neurons are increasingly appreciated as important reactive components of the glioma microenvironment, due to their role in causing hallmark glioma symptoms, such as cognitive deficits and seizures, as well as their potential ability to drive glioma progression. Separately, mTOR signaling has been shown to have pleiotropic effects in the brain tumor microenvironment, including regulation of neuronal hyperexcitability. However, the local cellular-level effects of mTOR inhibition on glioma-induced neuronal alterations are not well understood. Here we employed neuron-specific profiling of ribosome-bound mRNA via 'RiboTag,' morphometric analysis of dendritic spines, and in vivo calcium imaging, along with pharmacological mTOR inhibition to investigate the impact of glioma burden and mTOR inhibition on these neuronal alterations. The RiboTag analysis of tumor-associated excitatory neurons showed a downregulation of transcripts encoding excitatory and inhibitory postsynaptic proteins and dendritic spine development, and an upregulation of transcripts encoding cytoskeletal proteins involved in dendritic spine turnover. Light and electron microscopy of tumor-associated excitatory neurons demonstrated marked decreases in dendritic spine density. In vivo two-photon calcium imaging in tumor-associated excitatory neurons revealed progressive alterations in neuronal activity, both at the population and single-neuron level, throughout tumor growth. This in vivo calcium imaging also revealed altered stimulus-evoked somatic calcium events, with changes in event rate, size, and temporal alignment to stimulus, which was most pronounced in neurons with high-tumor burden. A single acute dose of AZD8055, a combined mTORC1/2 inhibitor, reversed the glioma-induced alterations on the excitatory neurons, including the alterations in ribosome-bound transcripts, dendritic spine density, and stimulus evoked responses seen by calcium imaging. These results point to mTOR-driven pathological plasticity in neurons at the infiltrative margin of glioma - manifested by alterations in ribosome-bound mRNA, dendritic spine density, and stimulus-evoked neuronal activity. Collectively, our work identifies the pathological changes that tumor-associated excitatory neurons experience as both hyperlocal and reversible under the influence of mTOR inhibition, providing a foundation for developing therapies targeting neuronal signaling in glioma.

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.

Mapping and comparing fMRI connectivity networks across species

Technical advances in neuroimaging, notably in fMRI, have allowed distributed patterns of functional connectivity to be mapped in the human brain with increasing spatiotemporal resolution. Recent years have seen a growing interest in extending this approach to rodents and non-human primates to understand the mechanism of fMRI connectivity and complement human investigations of the functional connectome. Here, we discuss current challenges and opportunities of fMRI connectivity mapping across species. We underscore the critical importance of physiologically decoding neuroimaging measures of brain (dys)connectivity via multiscale mechanistic investigations in animals. We next highlight a set of general principles governing the organization of mammalian connectivity networks across species. These include the presence of evolutionarily conserved network systems, a dominant cortical axis of functional connectivity, and a common repertoire of topographically conserved fMRI spatiotemporal modes. We finally describe emerging approaches allowing comparisons and extrapolations of fMRI connectivity findings across species. As neuroscientists gain access to increasingly sophisticated perturbational, computational and recording tools, cross-species fMRI offers novel opportunities to investigate the large-scale organization of the mammalian brain in health and disease.

Large-scale momentary brain co-activation patterns are associated with hyperalgesia and mediate focal neurochemistry and cross-network functional connectivity in fibromyalgia

ABSTRACT

Fibromyalgia has been characterized by augmented cross-network functional communication between the brain's sensorimotor, default mode, and attentional (salience/ventral and dorsal) networks. However, the underlying mechanisms of these aberrant communication patterns are unknown. In this study, we sought to understand large-scale topographic patterns at instantaneous timepoints, known as co-activation patterns (CAPs). We found that a sustained pressure pain challenge temporally modulated the occurrence of CAPs. Using proton magnetic resonance spectroscopy, we found that greater basal excitatory over inhibitory neurotransmitter levels within the anterior insula orchestrated higher cross-network connectivity between the anterior insula and the default mode network through lower occurrence of a CAP encompassing the attentional networks during sustained pain. Moreover, we found that hyperalgesia in fibromyalgia was mediated through increased occurrence of a CAP encompassing the sensorimotor network during sustained pain. In conclusion, this study elucidates the role of momentary large-scale topographic brain patterns in shaping noxious information in patients with fibromyalgia, while laying the groundwork for using precise spatiotemporal dynamics of the brain for the potential development of therapeutics.

Topological changes of fast large-scale brain dynamics in mild cognitive impairment predict early memory impairment: a resting-state, source reconstructed, magnetoencephalography study

Functional connectivity has been used as a framework to investigate widespread brain interactions underlying cognitive deficits in mild cognitive impairment (MCI). However, many functional connectivity metrics focus on the average of the periodic activities, disregarding the aperiodic bursts of activity (i.e., the neuronal avalanches) characterizing the large-scale dynamic activities of the brain. Here, we apply the recently described avalanche transition matrix framework to source-reconstructed magnetoencephalography signals in a cohort of 32 MCI patients and 32 healthy controls to describe the spatio-temporal features of neuronal avalanches and explore their topological properties. Our results showed that MCI patients showed a more centralized network (as assessed by higher values of the degree divergence and leaf fraction) as compared to healthy controls. Furthermore, we found that the degree divergence (in the theta band) was predictive of hippocampal memory impairment. These findings highlight the role of the changes of aperiodic bursts in clinical conditions and may contribute to a more thorough phenotypical assessment of patients.

A thalamocortical substrate for integrated information via critical synchronous bursting

Understanding the neurobiological mechanisms underlying consciousness remains a significant challenge. Recent evidence suggests that the coupling between distal-apical and basal-somatic dendrites in thick-tufted layer 5 pyramidal neurons (L5PN), regulated by the nonspecific-projecting thalamus, is crucial for consciousness. Yet, it is uncertain whether this thalamocortical mechanism can support emergent signatures of consciousness, such as integrated information. To address this question, we constructed a biophysical network of dual-compartment thick-tufted L5PN, with dendrosomatic coupling controlled by thalamic inputs. Our findings demonstrate that integrated information is maximized when nonspecific thalamic inputs drive the system into a regime of time-varying synchronous bursting. Here, the system exhibits variable spiking dynamics with broad pairwise correlations, supporting the enhanced integrated information. Further, the observed peak in integrated information aligns with criticality signatures and empirically observed layer 5 pyramidal bursting rates. These results suggest that the thalamocortical core of the mammalian brain may be evolutionarily configured to optimize effective information processing, providing a potential neuronal mechanism that integrates microscale theories with macroscale signatures of consciousness.

Living on the edge: network neuroscience beyond nodes

Network neuroscience has emphasized the connectional properties of neural elements - cells, populations, and regions. This has come at the expense of the anatomical and functional connections that link these elements to one another. A new perspective - namely one that emphasizes 'edges' - may prove fruitful in addressing outstanding questions in network neuroscience. We highlight one recently proposed 'edge-centric' method and review its current applications, merits, and limitations. We also seek to establish conceptual and mathematical links between this method and previously proposed approaches in the network science and neuroimaging literature. We conclude by presenting several avenues for future work to extend and refine existing edge-centric analysis.

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

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

Beyond pulsed inhibition: Alpha oscillations modulate attenuation and amplification of neural activity in the awake resting state

Alpha oscillations are a distinctive feature of the awake resting state of the human brain. However, their functional role in resting-state neuronal dynamics remains poorly understood. Here we show that, during resting wakefulness, alpha oscillations drive an alternation of attenuation and amplification bouts in neural activity. Our analysis indicates that inhibition is activated in pulses that last for a single alpha cycle and gradually suppress neural activity, while excitation is successively enhanced over a few alpha cycles to amplify neural activity. Furthermore, we show that long-term alpha amplitude fluctuations-the "waxing and waning" phenomenon-are an attenuation-amplification mechanism described by a power-law decay of the activity rate in the "waning" phase. Importantly, we do not observe such dynamics during non-rapid eye movement (NREM) sleep with marginal alpha oscillations. The results suggest that alpha oscillations modulate neural activity not only through pulses of inhibition (pulsed inhibition hypothesis) but also by timely enhancement of excitation (or disinhibition).

Frequency dependent whole-brain coactivation patterns analysis in Alzheimer's disease

Background

The brain in resting state has complex dynamic properties and shows frequency dependent characteristics. The frequency-dependent whole-brain dynamic changes of resting state across the scans have been ignored in Alzheimer's disease (AD).

Objective

Coactivation pattern (CAP) analysis can identify different brain states. This paper aimed to investigate the dynamic characteristics of frequency dependent whole-brain CAPs in AD.

Methods

We utilized a multiband CAP approach to model the state space and study brain dynamics in both AD and NC. The correlation between the dynamic characteristics and the subjects' clinical index was further analyzed.

Results

The results showed similar CAP patterns at different frequency bands, but the occurrence of patterns was different. In addition, CAPs associated with the default mode network (DMN) and the ventral/dorsal visual network (dorsal/ventral VN) were altered significantly between the AD and NC groups. This study also found the correlation between the altered dynamic characteristics of frequency dependent CAPs and the patients' clinical Mini-Mental State Examination assessment scale scores.

Conclusion

This study revealed that while similar CAP spatial patterns appear in different frequency bands, their dynamic characteristics in subbands vary. In addition, delineating subbands was more helpful in distinguishing AD from NC in terms of CAP.

Criticality of neuronal avalanches in human sleep and their relationship with sleep macro- and micro-architecture

Sleep plays a key role in preserving brain function, keeping brain networks in a state that ensures optimal computation. Empirical evidence indicates that this state is consistent with criticality, where scale-free neuronal avalanches emerge. However, the connection between sleep architecture and brain tuning to criticality remains poorly understood. Here, we characterize the critical behavior of avalanches and study their relationship with sleep macro- and micro-architectures, in particular, the cyclic alternating pattern (CAP). We show that avalanches exhibit robust scaling behaviors, with exponents obeying scaling relations consistent with the mean-field directed percolation universality class. We demonstrate that avalanche dynamics is modulated by the NREM-REM cycles and that, within NREM sleep, avalanche occurrence correlates with CAP activation phases-indicating a potential link between CAP and brain tuning to criticality. The results open new perspectives on the collective dynamics underlying CAP function, and on the relationship between sleep architecture, avalanches, and self-organization to criticality.

Hierarchical organization of spontaneous co-fluctuations in densely sampled individuals using fMRI

Edge time series decompose functional connectivity into its framewise contributions. Previous studies have focused on characterizing the properties of high-amplitude frames (time points when the global co-fluctuation amplitude takes on its largest value), including their cluster structure. Less is known about middle- and low-amplitude co-fluctuations (peaks in co-fluctuation time series but of lower amplitude). Here, we directly address those questions, using data from two dense-sampling studies: the MyConnectome project and Midnight Scan Club. We develop a hierarchical clustering algorithm to group peak co-fluctuations of all magnitudes into nested and multiscale clusters based on their pairwise concordance. At a coarse scale, we find evidence of three large clusters that, collectively, engage virtually all canonical brain systems. At finer scales, however, each cluster is dissolved, giving way to increasingly refined patterns of co-fluctuations involving specific sets of brain systems. We also find an increase in global co-fluctuation magnitude with hierarchical scale. Finally, we comment on the amount of data needed to estimate co-fluctuation pattern clusters and implications for brain-behavior studies. Collectively, the findings reported here fill several gaps in current knowledge concerning the heterogeneity and richness of co-fluctuation patterns as estimated with edge time series while providing some practical guidance for future studies.

Establishing brain states in neuroimaging data

The definition of a brain state remains elusive, with varying interpretations across different sub-fields of neuroscience-from the level of wakefulness in anaesthesia, to activity of individual neurons, voltage in EEG, and blood flow in fMRI. This lack of consensus presents a significant challenge to the development of accurate models of neural dynamics. However, at the foundation of dynamical systems theory lies a definition of what constitutes the 'state' of a system-i.e., a specification of the system's future. Here, we propose to adopt this definition to establish brain states in neuroimaging timeseries by applying Dynamic Causal Modelling (DCM) to low-dimensional embedding of resting and task condition fMRI data. We find that ~90% of subjects in resting conditions are better described by first-order models, whereas ~55% of subjects in task conditions are better described by second-order models. Our work calls into question the status quo of using first-order equations almost exclusively within computational neuroscience and provides a new way of establishing brain states, as well as their associated phase space representations, in neuroimaging datasets.

Improved dynamic functional connectivity estimation with an alternating hidden Markov model

Dynamic functional connectivity (DFC) analysis has been widely applied to functional magnetic resonance imaging (fMRI) data to reveal the time-varying functional interactions between brain regions. Although the sliding window (SW) method is popular for DFC analysis, the selection of window length is hard, and the temporal resolution is limited by the window length. The hidden Markov model (HMM) without the limitation of window length has been proven to be able to estimate time-varying brain states from fMRI data. However, HMM tends to be overfitted in DFC analysis of fMRI data because of the high spatial dimension and the limited sample size of fMRI data. In this study, we proposed an alternating HMM (aHMM) method that used the functional connectivity estimation of SW to initialize the covariance matrix of HMM and adopted an alternating HMM procedure to reduce the number of parameters during each optimization. The simulated and real fMRI resting data from the Human Connectome Projects showed that aHMM produced better robustness to noise, parameter number and sample size in DFC estimation than SW and HMM. For the real fMRI resting data of cerebral small vessel disease (CSVD), results of aHMM revealed that amnesia and mild cognitive impairment (aMCI) caused the CSVD with aMCI (CSVD-aMCI) group tended to spend more time on the brain state with overall weak connections and less time on the state with overall strong connections than the CSVD-controls. Moreover, CSVD-aMCI showed significantly lower connectivity amplitude and higher connectivity fluctuation than CSVD-control. In contrast, HMM did not detect intergroup differences of the connectivity amplitude and fluctuations and SW did not detect intergroup differences of connectivity fluctuations and fraction of time. The results further indicated that aHMM outperformed HMM and SW in detecting inter-group differences of temporal properties of DFC and connectivity fluctuations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-022-09874-3.

Scale-free correlations in the dynamics of a small (N∼10000) cortical network

The advent of novel optogenetics technology allows the recording of brain activity with a resolution never seen before. The characterization of these very large data sets offers new challenges as well as unique theory-testing opportunities. Here we discuss whether the spatial and temporal correlations of the collective activity of thousands of neurons are tangled as predicted by the theory of critical phenomena. The analysis shows that both the correlation length ξ and the correlation time τ scale as predicted as a function of the system size. With some peculiarities that we discuss, the analysis uncovers evidence consistent with the view that the large-scale brain cortical dynamics corresponds to critical phenomena.

The neurophysiologic landscape of the sleep onset: a systematic review

Background

The sleep onset process is an ill-defined complex process of transition from wakefulness to sleep, characterized by progressive modifications at the subjective, behavioural, cognitive, and physiological levels. To this date, there is no international consensus which could aid a principled characterisation of this process for clinical research purposes. The current review aims to systemise the current knowledge about the underlying mechanisms of the natural heterogeneity of this process.

Methods

In this systematic review, studies investigating the process of the sleep onset from 1970 to 2022 were identified using electronic database searches of PsychINFO, MEDLINE, and Embase.

Results

A total of 139 studies were included; 110 studies in healthy participants and 29 studies in participants with sleep disorders. Overall, there is a limited consensus across a body of research about what distinct biomarkers of the sleep onset constitute. Only sparse data exists on the physiology, neurophysiology and behavioural mechanisms of the sleep onset, with majority of studies concentrating on the non-rapid eye movement stage 2 (NREM 2) as a potentially better defined and a more reliable time point that separates sleep from the wake, on the sleep wake continuum.

Conclusions

The neurophysiologic landscape of sleep onset bears a complex pattern associated with a multitude of behavioural and physiological markers and remains poorly understood. The methodological variation and a heterogenous definition of the wake-sleep transition in various studies to date is understandable, given that sleep onset is a process that has fluctuating and ill-defined boundaries. Nonetheless, the principled characterisation of the sleep onset process is needed which will allow for a greater conceptualisation of the mechanisms underlying this process, further influencing the efficacy of current treatments for sleep disorders.

Leading basic modes of spontaneous activity drive individual functional connectivity organization in the resting human brain

Spontaneous activity of the human brain provides a window to explore intrinsic principles of functional organization. However, most studies have focused on interregional functional connectivity. The principles underlying rich repertoires of instantaneous activity remain largely unknown. We apply a recently proposed eigen-microstate analysis to three resting-state functional MRI datasets to identify basic modes that represent fundamental activity patterns that coexist over time. We identify five leading basic modes that dominate activity fluctuations. Each mode exhibits a distinct functional system-dependent coactivation pattern and corresponds to specific cognitive profiles. In particular, the spatial pattern of the first leading basis mode shows the separation of activity between the default-mode and primary and attention regions. Based on theoretical modelling, we further reconstruct individual functional connectivity as the weighted superposition of coactivation patterns corresponding to these leading basic modes. Moreover, these leading basic modes capture sleep deprivation-induced changes in brain activity and interregional connectivity, primarily involving the default-mode and task-positive regions. Our findings reveal a dominant set of basic modes of spontaneous activity that reflect multiplexed interregional coordination and drive conventional functional connectivity, furthering the understanding of the functional significance of spontaneous brain activity.

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.

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.

The confound of hemodynamic response function variability in human resting-state functional MRI studies

Functional magnetic resonance imaging (fMRI) is an indirect measure of neural activity with the hemodynamic response function (HRF) coupling it with unmeasured neural activity. The HRF, modulated by several non-neural factors, is variable across brain regions, individuals and populations. Yet, a majority of human resting-state fMRI connectivity studies continue to assume a non-variable HRF. In this article, with supportive prior evidence, we argue that HRF variability cannot be ignored as it substantially confounds within-subject connectivity estimates and between-subjects connectivity group differences. We also discuss its clinical relevance with connectivity impairments confounded by HRF aberrations in several disorders. We present limited data on HRF differences between women and men, which resulted in a 15.4% median error in functional connectivity estimates in a group-level comparison. We also discuss the implications of HRF variability for fMRI studies in the spinal cord. There is a need for more dialogue within the community on the HRF confound, and we hope that our article is a catalyst in the process.

Intermediately synchronised brain states optimise trade-off between subject specificity and predictive capacity

Functional connectivity (FC) refers to the statistical dependencies between activity of distinct brain areas. To study temporal fluctuations in FC within the duration of a functional magnetic resonance imaging (fMRI) scanning session, researchers have proposed the computation of an edge time series (ETS) and their derivatives. Evidence suggests that FC is driven by a few time points of high-amplitude co-fluctuation (HACF) in the ETS, which may also contribute disproportionately to interindividual differences. However, it remains unclear to what degree different time points actually contribute to brain-behaviour associations. Here, we systematically evaluate this question by assessing the predictive utility of FC estimates at different levels of co-fluctuation using machine learning (ML) approaches. We demonstrate that time points of lower and intermediate co-fluctuation levels provide overall highest subject specificity as well as highest predictive capacity of individual-level phenotypes.

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.