Noise during rest enables the exploration of the brain's dynamic repertoire

Neuroimaging markers of aberrant brain activity and treatment response in schizophrenia patients based on brain complexity

The complexity of brain activity reflects its ability to process information, adapt to environmental changes, and transition between states. However, it remains unclear how schizophrenia (SZ) affects brain activity complexity, particularly its dynamic changes. This study aimed to investigate the abnormal patterns of brain activity complexity in SZ, their relationship with cognitive deficits, and the impact of antipsychotic medication. Forty-four drug-naive first-episode (DNFE) SZ patients and thirty demographically matched healthy controls (HC) were included. Functional MRI-based sliding window analysis was utilized for the first time to calculate weighted permutation entropy to characterize complex patterns of brain activity in SZ patients before and after 12 weeks of risperidone treatment. Results revealed reduced complexity in the caudate, putamen, and pallidum at baseline in SZ patients compared to HC, with reduced complexity in the left caudate positively correlated with Continuous Performance Test (CPT) and Category Fluency Test scores. After treatment, the complexity of the left caudate increased. Regions with abnormal complexity showed decreased functional connectivity, with complexity positively correlated with connectivity strength. We observed that the dynamic complexity of the brain exhibited the characteristic of spontaneous, recurring "complexity drop", potentially reflecting transient state transitions in the resting brain. Compared to HC, patients exhibited reduced scope, intensity, and duration of complexity drop, all of which improved after treatment. Reduced duration was negatively correlated with CPT scores and positively with clinical symptoms. The results suggest that abnormalities in brain activity complexity and its dynamic changes may underlie cognitive deficits and clinical symptoms in SZ patients. Antipsychotic treatment partially restores these abnormalities, highlighting their potential as indicators of treatment efficacy and biomarkers for personalized therapy.

Resting-state neural dynamics changes in older adults with post-COVID syndrome and the modulatory effect of cognitive training and sex

Post-COVID syndrome manifests with numerous neurological and cognitive symptoms, the precise origins of which are still not fully understood. As females and older adults are more susceptible to developing this condition, our study aimed to investigate how post-COVID syndrome alters intrinsic brain dynamics in older adults and whether biological sex and cognitive training might modulate these effects, with a specific focus on older females. The participants, aged between 60 and 75 years, were divided into three experimental groups: healthy old female, post-COVID old female and post-COVID old male. They underwent an adaptive task-switching training protocol. We analysed multiscale entropy and spectral power density of resting-state EEG data collected before and after the training to assess neural signal complexity and oscillatory power, respectively. We found no difference between post-COVID females and males before training, indicating that post-COVID similarly affected both sexes. However, cognitive training was effective only in post-COVID females and not in males, by modulating local neural processing capacity. This improvement was further evidenced by comparing healthy and post-COVID females, wherein the latter group showed increased finer timescale entropy (1-30 ms) and higher frequency band power (11-40 Hz) before training, but these differences disappeared following cognitive training. Our results suggest that in older adults with post-COVID syndrome, there is a pronounced shift from more global to local neural processing, potentially contributing to accelerated neural aging in this condition. However, cognitive training seems to offer a promising intervention method for modulating these changes in brain dynamics, especially among females.

Reliability of variability and complexity measures for task and task-free BOLD fMRI

Brain activity continuously fluctuates over time, even if the brain is in controlled (e.g., experimentally induced) states. Recent years have seen an increasing interest in understanding the complexity of these temporal variations, for example with respect to developmental changes in brain function or between-person differences in healthy and clinical populations. However, the psychometric reliability of brain signal variability and complexity measures-which is an important precondition for robust individual differences as well as longitudinal research-is not yet sufficiently studied. We examined reliability (split-half correlations) and test-retest correlations for task-free (resting-state) BOLD fMRI as well as split-half correlations for seven functional task data sets from the Human Connectome Project to evaluate their reliability. We observed good to excellent split-half reliability for temporal variability measures derived from rest and task fMRI activation time series (standard deviation, mean absolute successive difference, mean squared successive difference), and moderate test-retest correlations for the same variability measures under rest conditions. Brain signal complexity estimates (several entropy and dimensionality measures) showed moderate to good reliabilities under both, rest and task activation conditions. We calculated the same measures also for time-resolved (dynamic) functional connectivity time series and observed moderate to good reliabilities for variability measures, but poor reliabilities for complexity measures derived from functional connectivity time series. Global (i.e., mean across cortical regions) measures tended to show higher reliability than region-specific variability or complexity estimates. Larger subcortical regions showed similar reliability as cortical regions, but small regions showed lower reliability, especially for complexity measures. Lastly, we also show that reliability scores are only minorly dependent on differences in scan length and replicate our results across different parcellation and denoising strategies. These results suggest that the variability and complexity of BOLD activation time series are robust measures well-suited for individual differences research. Temporal variability of global functional connectivity over time provides an important novel approach to robustly quantifying the dynamics of brain function. PRACTITIONER POINTS: Variability and complexity measures of BOLD activation show good split-half reliability and moderate test-retest reliability. Measures of variability of global functional connectivity over time can robustly quantify neural dynamics. Length of fMRI data has only a minor effect on reliability.

Inducing a meditative state by artificial perturbations: A mechanistic understanding of brain dynamics underlying meditation

Contemplative neuroscience has increasingly explored meditation using neuroimaging. However, the brain mechanisms underlying meditation remain elusive. Here, we implemented a mechanistic framework to explore the spatiotemporal dynamics of expert meditators during meditation and rest, and controls during rest. We first applied a model-free approach by defining a probabilistic metastable substate (PMS) space for each condition, consisting of different probabilities of occurrence from a repertoire of dynamic patterns. Moreover, we implemented a model-based approach by adjusting the PMS of each condition to a whole-brain model, which enabled us to explore in silico perturbations to transition from resting-state to meditation and vice versa. Consequently, we assessed the sensitivity of different brain areas regarding their perturbability and their mechanistic local-global effects. Overall, our work reveals distinct whole-brain dynamics in meditation compared to rest, and how transitions can be induced with localized artificial perturbations. It motivates future work regarding meditation as a practice in health and as a potential therapy for brain disorders.

Disrupted topological properties of structural brain networks present a glutamatergic neuropathophysiology in people with narcolepsy

STUDY OBJECTIVES

Growing evidences have documented various abnormalities of the white matter bundles in people with narcolepsy. We sought to evaluate topological properties of brain structural networks, and their association with symptoms and neuropathophysiological features in people with narcolepsy.

METHODS

Diffusion tensor imaging was conducted for people with narcolepsy (n = 30) and matched healthy controls as well as symptoms assessment. Structural connectivity for each participant was generated to analyze global and regional topological properties and their correlations with narcoleptic features. Further human brain transcriptome was extracted and spatially registered for connectivity vulnerability. Genetic functional enrichment analysis was performed and further clarified using in vivo emission computed tomography data.

RESULTS

A wide and dramatic decrease in structural connectivities was observed in people with narcolepsy, with descending network degree and global efficiency. These metrics were not only correlated with sleep latency and awakening features, but also reflected alterations of sleep macrostructure in people with narcolepsy. Network-based statistics identified a small hyperenhanced subnetwork of cingulate gyrus that was closely related to rapid eye movement sleep behavior disorder (RBD) in narcolepsy. Further imaging genetics analysis suggested glutamatergic signatures were responsible for the preferential vulnerability of connectivity alterations in people with narcolepsy, while additional PET/SPECT data verified that structural alteration was significantly correlated with metabotropic glutamate receptor 5 (mGlutR5) and N-methyl-D-aspartate receptor (NMDA).

CONCLUSIONS

People with narcolepsy endured a remarkable decrease in the structural architecture, which was not only closely related to narcolepsy symptoms but also glutamatergic signatures.

In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Cholinergic modulation supports dynamic switching of resting state networks through selective DMN suppression

Brain activity during the resting state is widely used to examine brain organization, cognition and alterations in disease states. While it is known that neuromodulation and the state of alertness impact resting-state activity, neural mechanisms behind such modulation of resting-state activity are unknown. In this work, we used a computational model to demonstrate that change in excitability and recurrent connections, due to cholinergic modulation, impacts resting-state activity. The results of such modulation in the model match closely with experimental work on direct cholinergic modulation of Default Mode Network (DMN) in rodents. We further extended our study to the human connectome derived from diffusion-weighted MRI. In human resting-state simulations, an increase in cholinergic input resulted in a brain-wide reduction of functional connectivity. Furthermore, selective cholinergic modulation of DMN closely captured experimentally observed transitions between the baseline resting state and states with suppressed DMN fluctuations associated with attention to external tasks. Our study thus provides insight into potential neural mechanisms for the effects of cholinergic neuromodulation on resting-state activity and its dynamics.

Virtual brain twins: from basic neuroscience to clinical use

Virtual brain twins are personalized, generative and adaptive brain models based on data from an individual's brain for scientific and clinical use. After a description of the key elements of virtual brain twins, we present the standard model for personalized whole-brain network models. The personalization is accomplished using a subject's brain imaging data by three means: (1) assemble cortical and subcortical areas in the subject-specific brain space; (2) directly map connectivity into the brain models, which can be generalized to other parameters; and (3) estimate relevant parameters through model inversion, typically using probabilistic machine learning. We present the use of personalized whole-brain network models in healthy ageing and five clinical diseases: epilepsy, Alzheimer's disease, multiple sclerosis, Parkinson's disease and psychiatric disorders. Specifically, we introduce spatial masks for relevant parameters and demonstrate their use based on the physiological and pathophysiological hypotheses. Finally, we pinpoint the key challenges and future directions.

Investigating the Effect of Contextual Cueing with Face Stimuli on Electrophysiological Measures in Younger and Older Adults

Extracting repeated patterns from our surroundings plays a crucial role in contextualizing information, making predictions, and guiding our behavior implicitly. Previous research showed that contextual cueing enhances visual search performance in younger adults. In this study, we investigated whether contextual cueing could also improve older adults' performance and whether age-related differences in the neural processes underlying implicit contextual learning could be detected. Twenty-four younger and 25 older participants performed a visual search task with contextual cueing. Contextual information was generated using repeated face configurations alongside random new configurations. We measured RT difference between new and repeated configurations; ERPs to uncover the neural processes underlying contextual cueing for early (N2pc), intermediate (P3b), and late (r-LRP) processes; and multiscale entropy and spectral power density analyses to examine neural dynamics. Both younger and older adults showed similar contextual cueing benefits in their visual search efficiency at the behavioral level. In addition, they showed similar patterns regarding contextual information processing: Repeated face configurations evoked decreased finer timescale entropy (1-20 msec) and higher frequency band power (13-30 Hz) compared with new configurations. However, we detected age-related differences in ERPs: Younger, but not older adults, had larger N2pc and P3b components for repeated compared with new configurations. These results suggest that contextual cueing remains intact with aging. Although attention- and target-evaluation-related ERPs differed between the age groups, the neural dynamics of contextual learning were preserved with aging, as both age groups increasingly utilized more globally grouped representations for repeated face configurations during the learning process.

In-vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here we non-invasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically-plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the GABA-agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 years old) and Asian (7.2 to 7.9 years old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Modeling brain network flexibility in networks of coupled oscillators: a feasibility study

Modeling the functionality of the human brain is a major goal in neuroscience for which many powerful methodologies have been developed over the last decade. The impact of working memory and the associated brain regions on the brain dynamics is of particular interest due to their connection with many functions and malfunctions in the brain. In this context, the concept of brain flexibility has been developed for the characterization of brain functionality. We discuss emergence of brain flexibility that is commonly measured by the identification of changes in the cluster structure of co-active brain regions. We provide evidence that brain flexibility can be modeled by a system of coupled FitzHugh-Nagumo oscillators where the network structure is obtained from human brain Diffusion Tensor Imaging (DTI). Additionally, we propose a straightforward and computationally efficient alternative macroscopic measure, which is derived from the Pearson distance of functional brain matrices. This metric exhibits similarities to the established patterns of brain template flexibility that have been observed in prior investigations. Furthermore, we explore the significance of the brain's network structure and the strength of connections between network nodes or brain regions associated with working memory in the observation of patterns in networks flexibility. This work enriches our understanding of the interplay between the structure and function of dynamic brain networks and proposes a modeling strategy to study brain flexibility.

The Hopf whole-brain model and its linear approximation

Whole-brain models have proven to be useful to understand the emergence of collective activity among neural populations or brain regions. These models combine connectivity matrices, or connectomes, with local node dynamics, noise, and, eventually, transmission delays. Multiple choices for the local dynamics have been proposed. Among them, nonlinear oscillators corresponding to a supercritical Hopf bifurcation have been used to link brain connectivity and collective phase and amplitude dynamics in different brain states. Here, we studied the linear fluctuations of this model to estimate its stationary statistics, i.e., the instantaneous and lagged covariances and the power spectral densities. This linear approximation-that holds in the case of heterogeneous parameters and time-delays-allows analytical estimation of the statistics and it can be used for fast parameter explorations to study changes in brain state, changes in brain activity due to alterations in structural connectivity, and modulations of parameter due to non-equilibrium dynamics.

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.

Dynamical noise can enhance high-order statistical structure in complex systems

Recent research has provided a wealth of evidence highlighting the pivotal role of high-order interdependencies in supporting the information-processing capabilities of distributed complex systems. These findings may suggest that high-order interdependencies constitute a powerful resource that is, however, challenging to harness and can be readily disrupted. In this paper, we contest this perspective by demonstrating that high-order interdependencies can not only exhibit robustness to stochastic perturbations, but can in fact be enhanced by them. Using elementary cellular automata as a general testbed, our results unveil the capacity of dynamical noise to enhance the statistical regularities between agents and, intriguingly, even alter the prevailing character of their interdependencies. Furthermore, our results show that these effects are related to the high-order structure of the local rules, which affect the system's susceptibility to noise and characteristic time scales. These results deepen our understanding of how high-order interdependencies may spontaneously emerge within distributed systems interacting with stochastic environments, thus providing an initial step toward elucidating their origin and function in complex systems like the human brain.

Spatiotemporal dedifferentiation of the global brain signal topography along the adult lifespan

Age-related variations in many regions and/or networks of the human brain have been uncovered using resting-state functional magnetic resonance imaging. However, these findings did not account for the dynamical effect the brain's global activity (global signal [GS]) causes on local characteristics, which is measured by GS topography. To address this gap, we tested GS topography including its correlation with age using a large-scale cross-sectional adult lifespan dataset (n = 492). Both GS topography and its variation with age showed frequency-specific patterns, reflecting the spatiotemporal characteristics of the dynamic change of GS topography with age. A general trend toward dedifferentiation of GS topography with age was observed in both spatial (i.e., less differences of GS between different regions) and temporal (i.e., less differences of GS between different frequencies) dimensions. Further, methodological control analyses suggested that although most age-related dedifferentiation effects remained across different preprocessing strategies, some were triggered by neuro-vascular coupling and physiological noises. Together, these results provide the first evidence for age-related effects on global brain activity and its topographic-dynamic representation in terms of spatiotemporal dedifferentiation.

Violations of the fluctuation-dissipation theorem reveal distinct nonequilibrium dynamics of brain states

The brain is a nonequilibrium system whose dynamics change in different brain states, such as wakefulness and deep sleep. Thermodynamics provides the tools for revealing these nonequilibrium dynamics. We used violations of the fluctuation-dissipation theorem to describe the hierarchy of nonequilibrium dynamics associated with different brain states. Together with a whole-brain model fitted to empirical human neuroimaging data, and deriving the appropriate analytical expressions, we were able to capture the deviation from equilibrium in different brain states that arises from asymmetric interactions and hierarchical organization.

The influence of inter-regional delays in generating large-scale brain networks of phase synchronization

Large-scale networks of phase synchronization are considered to regulate the communication between brain regions fundamental to cognitive function, but the mapping to their structural substrates, i.e., the structure-function relationship, remains poorly understood. Biophysical Network Models (BNMs) have demonstrated the influences of local oscillatory activity and inter-regional anatomical connections in generating alpha-band (8-12 Hz) networks of phase synchronization observed with Electroencephalography (EEG) and Magnetoencephalography (MEG). Yet, the influence of inter-regional conduction delays remains unknown. In this study, we compared a BNM with standard "distance-dependent delays", which assumes constant conduction velocity, to BNMs with delays specified by two alternative methods accounting for spatially varying conduction velocities, "isochronous delays" and "mixed delays". We followed the Approximate Bayesian Computation (ABC) workflow, i) specifying neurophysiologically informed prior distributions of BNM parameters, ii) verifying the suitability of the prior distributions with Prior Predictive Checks, iii) fitting each of the three BNMs to alpha-band MEG resting-state data (N = 75) with Bayesian optimization for Likelihood-Free Inference (BOLFI), and iv) choosing between the fitted BNMs with ABC model comparison on a separate MEG dataset (N = 30). Prior Predictive Checks revealed the range of dynamics generated by each of the BNMs to encompass those seen in the MEG data, suggesting the suitability of the prior distributions. Fitting the models to MEG data yielded reliable posterior distributions of the parameters of each of the BNMs. Finally, model comparison revealed the BNM with "distance-dependent delays", as the most probable to describe the generation of alpha-band networks of phase synchronization seen in MEG. These findings suggest that distance-dependent delays might contribute to the neocortical architecture of human alpha-band networks of phase synchronization. Hence, our study illuminates the role of inter-regional delays in generating the large-scale networks of phase synchronization that might subserve the communication between regions vital to cognition.

Cross-attractor modeling of resting-state functional connectivity in psychiatric disorders

Resting-state functional connectivity (RSFC) is altered across various psychiatric disorders. Brain network modeling (BNM) has the potential to reveal the neurobiological underpinnings of such abnormalities by dynamically modeling the structure-function relationship and examining biologically relevant parameters after fitting the models with real data. Although innovative BNM approaches have been developed, two main issues need to be further addressed. First, previous BNM approaches are primarily limited to simulating noise-driven dynamics near a chosen attractor (or a stable brain state). An alternative approach is to examine multi(or cross)-attractor dynamics, which can be used to better capture non-stationarity and switching between states in the resting brain. Second, previous BNM work is limited to characterizing one disorder at a time. Given the large degree of co-morbidity across psychiatric disorders, comparing BNMs across disorders might provide a novel avenue to generate insights regarding the dynamical features that are common across (vs. specific to) disorders. Here, we address these issues by (1) examining the layout of the attractor repertoire over the entire multi-attractor landscape using a recently developed cross-attractor BNM approach; and (2) characterizing and comparing multiple disorders (schizophrenia, bipolar, and ADHD) with healthy controls using an openly available and moderately large multimodal dataset from the UCLA Consortium for Neuropsychiatric Phenomics. Both global and local differences were observed across disorders. Specifically, the global coupling between regions was significantly decreased in schizophrenia patients relative to healthy controls. At the same time, the ratio between local excitation and inhibition was significantly higher in the schizophrenia group than the ADHD group. In line with these results, the schizophrenia group had the lowest switching costs (energy gaps) across groups for several networks including the default mode network. Paired comparison also showed that schizophrenia patients had significantly lower energy gaps than healthy controls for the somatomotor and visual networks. Overall, this study provides preliminary evidence supporting transdiagnostic multi-attractor BNM approaches to better understand psychiatric disorders' pathophysiology.

Multiscale entropy of ADHD children during resting state condition

This present study aims to investigate neural mechanisms underlying ADHD compared to healthy children through the analysis of the complexity and the variability of the EEG brain signal using multiscale entropy (MSE), EEG signal standard deviation (SDs), as well as the mean, standard deviation (SDp) and coefficient of variation (CV) of absolute spectral power (PSD). For this purpose, a sample of children diagnosed with attention-deficit/hyperactivity disorder (ADHD) between 6 and 17 years old were selected based on the number of trials and diagnostic agreement, 32 for the open-eyes (OE) experimental condition and 25 children for the close-eyes (CE) experimental condition. Healthy control subjects were age- and gender-matched with the ADHD group. The MSE and SDs of resting-state EEG activity were calculated on 34 time scales using a coarse-grained procedure. In addition, the PSD was averaged in delta, theta, alpha, and beta frequency bands, and its mean, SDp, and CV were calculated. The results show that the MSE changes with age during development, increases as the number of scales increases and has a higher amplitude in controls than in ADHD. The absolute PSD results show CV differences between subjects in low and beta frequency bands, with higher variability values in the ADHD group. All these results suggest an increased EEG variability and reduced complexity in ADHD compared to controls.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-022-09869-0.

The impact of regional heterogeneity in whole-brain dynamics in the presence of oscillations

Large variability exists across brain regions in health and disease, considering their cellular and molecular composition, connectivity, and function. Large-scale whole-brain models comprising coupled brain regions provide insights into the underlying dynamics that shape complex patterns of spontaneous brain activity. In particular, biophysically grounded mean-field whole-brain models in the asynchronous regime were used to demonstrate the dynamical consequences of including regional variability. Nevertheless, the role of heterogeneities when brain dynamics are supported by synchronous oscillating state, which is a ubiquitous phenomenon in brain, remains poorly understood. Here, we implemented two models capable of presenting oscillatory behavior with different levels of abstraction: a phenomenological Stuart-Landau model and an exact mean-field model. The fit of these models informed by structural- to functional-weighted MRI signal (T1w/T2w) allowed us to explore the implication of the inclusion of heterogeneities for modeling resting-state fMRI recordings from healthy participants. We found that disease-specific regional functional heterogeneity imposed dynamical consequences within the oscillatory regime in fMRI recordings from neurodegeneration with specific impacts on brain atrophy/structure (Alzheimer's patients). Overall, we found that models with oscillations perform better when structural and functional regional heterogeneities are considered, showing that phenomenological and biophysical models behave similarly at the brink of the Hopf bifurcation.

Nonlinear changes in delayed functional network topology in Alzheimer's disease: relationship with amyloid and tau pathology

BACKGROUND

Alzheimer's disease is a neurodegenerative disorder associated with the abnormal deposition of pathological processes, such as amyloid-ß and tau, which produces nonlinear changes in the functional connectivity patterns between different brain regions across the Alzheimer's disease continuum. However, the mechanisms underlying these nonlinear changes remain largely unknown. Here, we address this question using a novel method based on temporal or delayed correlations and calculate new whole-brain functional networks to tackle these mechanisms.

METHODS

To assess our method, we evaluated 166 individuals from the ADNI database, including amyloid-beta negative and positive cognitively normal subjects, patients with mild cognitive impairment, and patients with Alzheimer's disease dementia. We used the clustering coefficient and the global efficiency to measure the functional network topology and assessed their relationship with amyloid and tau pathology measured by positron emission tomography, as well as cognitive performance using tests measuring memory, executive function, attention, and global cognition.

RESULTS

Our study found nonlinear changes in the global efficiency, but not in the clustering coefficient, showing that the nonlinear changes in functional connectivity are due to an altered ability of brain regions to communicate with each other through direct paths. These changes in global efficiency were most prominent in early disease stages. However, later stages of Alzheimer's disease were associated with widespread network disruptions characterized by changes in both network measures. The temporal delays required for the detection of these changes varied across the Alzheimer's disease continuum, with shorter delays necessary to detect changes in early stages and longer delays necessary to detect changes in late stages. Both global efficiency and clustering coefficient showed quadratic associations with pathological amyloid and tau burden as well as cognitive decline.

CONCLUSIONS

This study suggests that global efficiency is a more sensitive indicator of network changes in Alzheimer's disease when compared to clustering coefficient. Both network properties were associated with pathology and cognitive performance, demonstrating their relevance in clinical settings. Our findings provide an insight into the mechanisms underlying nonlinear changes in functional network organization in Alzheimer's disease, suggesting that it is the lack of direct connections that drives these functional changes.

Hierarchical fluctuation shapes a dynamic flow linked to states of consciousness

Consciousness arises from the spatiotemporal neural dynamics, however, its relationship with neural flexibility and regional specialization remains elusive. We identified a consciousness-related signature marked by shifting spontaneous fluctuations along a unimodal-transmodal cortical axis. This simple signature is sensitive to altered states of consciousness in single individuals, exhibiting abnormal elevation under psychedelics and in psychosis. The hierarchical dynamic reflects brain state changes in global integration and connectome diversity under task-free conditions. Quasi-periodic pattern detection revealed that hierarchical heterogeneity as spatiotemporally propagating waves linking to arousal. A similar pattern can be observed in macaque electrocorticography. Furthermore, the spatial distribution of principal cortical gradient preferentially recapitulated the genetic transcription levels of the histaminergic system and that of the functional connectome mapping of the tuberomammillary nucleus, which promotes wakefulness. Combining behavioral, neuroimaging, electrophysiological, and transcriptomic evidence, we propose that global consciousness is supported by efficient hierarchical processing constrained along a low-dimensional macroscale gradient.

Non-invasive mapping of epileptogenic networks predicts surgical outcome

Epilepsy is increasingly considered a disorder of brain networks. Studying these networks with functional connectivity can help identify hubs that facilitate the spread of epileptiform activity. Surgical resection of these hubs may lead patients who suffer from drug-resistant epilepsy to seizure freedom. Here, we aim to map non-invasively epileptogenic networks, through the virtual implantation of sensors estimated with electric and magnetic source imaging, in patients with drug-resistant epilepsy. We hypothesize that highly connected hubs identified non-invasively with source imaging can predict the epileptogenic zone and the surgical outcome better than spikes localized with conventional source localization methods (dipoles). We retrospectively analysed simultaneous high-density electroencephalography (EEG) and magnetoencephalography data recorded from 37 children and young adults with drug-resistant epilepsy who underwent neurosurgery. Using source imaging, we estimated virtual sensors at locations where intracranial EEG contacts were placed. On data with and without spikes, we computed undirected functional connectivity between sensors/contacts using amplitude envelope correlation and phase locking value for physiologically relevant frequency bands. From each functional connectivity matrix, we generated an undirected network containing the strongest connections within sensors/contacts using the minimum spanning tree. For each sensor/contact, we computed graph centrality measures. We compared functional connectivity and their derived graph centrality of sensors/contacts inside resection for good (n = 22, ILAE I) and poor (n = 15, ILAE II-VI) outcome patients, tested their ability to predict the epileptogenic zone in good-outcome patients, examined the association between highly connected hubs removal and surgical outcome and performed leave-one-out cross-validation to support their prognostic value. We also compared the predictive values of functional connectivity with those of dipoles. Finally, we tested the reliability of virtual sensor measures via Spearman's correlation with intracranial EEG at population- and patient-level. We observed higher functional connectivity inside than outside resection (P < 0.05, Wilcoxon signed-rank test) for good-outcome patients, on data with and without spikes across different bands for intracranial EEG and electric/magnetic source imaging and few differences for poor-outcome patients. These functional connectivity measures were predictive of both the epileptogenic zone and outcome (positive and negative predictive values ≥55%, validated using leave-one-out cross-validation) outperforming dipoles on spikes. Significant correlations were found between source imaging and intracranial EEG measures (0.4 ≤ rho ≤ 0.9, P < 0.05). Our findings suggest that virtual implantation of sensors through source imaging can non-invasively identify highly connected hubs in patients with drug-resistant epilepsy, even in the absence of frank epileptiform activity. Surgical resection of these hubs predicts outcome better than dipoles.

Noise improves the association between effects of local stimulation and structural degree of brain networks

Stimulation to local areas remarkably affects brain activity patterns, which can be exploited to investigate neural bases of cognitive function and modify pathological brain statuses. There has been growing interest in exploring the fundamental action mechanisms of local stimulation. Nevertheless, how noise amplitude, an essential element in neural dynamics, influences stimulation-induced brain states remains unknown. Here, we systematically examine the effects of local stimulation by using a large-scale biophysical model under different combinations of noise amplitudes and stimulation sites. We demonstrate that noise amplitude nonlinearly and heterogeneously tunes the stimulation effects from both regional and network perspectives. Furthermore, by incorporating the role of the anatomical network, we show that the peak frequencies of unstimulated areas at different stimulation sites averaged across noise amplitudes are highly positively related to structural connectivity. Crucially, the association between the overall changes in functional connectivity as well as the alterations in the constraints imposed by structural connectivity with the structural degree of stimulation sites is nonmonotonically influenced by the noise amplitude, with the association increasing in specific noise amplitude ranges. Moreover, the impacts of local stimulation of cognitive systems depend on the complex interplay between the noise amplitude and average structural degree. Overall, this work provides theoretical insights into how noise amplitude and network structure jointly modulate brain dynamics during stimulation and introduces possibilities for better predicting and controlling stimulation outcomes.

Activity-State Entropy: A Novel Brain Entropy Measure Based on Spatial Patterns of Activity

BACKGROUND

Brain entropy is a measure of the complexity of brain activity that has been linked to various cognitive abilities. The measure is based on Shannon Entropy, a measure from Information Theory that quantifies the information capacity of a system from the probability distribution of its states. Most fMRI studies measure brain entropy at the voxel level as time-series entropy and assume that entropic time-series indicate complex large-scale spatiotemporal patterns of activity. New Method We developed a novel measure of brain entropy called Activity-State Entropy. The method quantifies entropy based on underlying patterns of coactivation identified using Principal Components Analysis. These patterns, termed eigenactivity states, combine in time-varying proportions.

RESULTS

We showed that Activity-State Entropy is sensitive to the complexity of the spatiotemporal patterns of activity in simulated fMRI data. We then applied this measure to real resting-state fMRI data and found that the eigenactivity states that explained the most variance in the data were comprised of large clusters of coactivating voxels, including clusters within Default Mode Network regions. More entropic brains were increasingly influenced by eigenactivity states comprised of smaller and more sparsely distributed clusters. Comparison to Existing Methods We compared Activity-State Entropy to Sample Entropy and Dispersion Entropy, two time-series entropy measures commonly used in neuroimaging research, and found all three measures were positively correlated.

CONCLUSIONS

Activity-State Entropy provides a measure of the spatiotemporal complexity of brain activity that complements time-series based measures of brain entropy.

Stability of sensorimotor network sculpts the dynamic repertoire of resting state over lifespan

Temporally stable patterns of neural coordination among distributed brain regions are crucial for survival. Recently, many studies highlight association between healthy aging and modifications in organization of functional brain networks, across various time-scales. Nonetheless, quantitative characterization of temporal stability of functional brain networks across healthy aging remains unexplored. This study introduces a data-driven unsupervised approach to capture high-dimensional dynamic functional connectivity (dFC) via low-dimensional patterns and subsequent estimation of temporal stability using quantitative metrics. Healthy aging related changes in temporal stability of dFC were characterized across resting-state, movie-viewing, and sensorimotor tasks (SMT) on a large (n = 645) healthy aging dataset (18-88 years). Prominent results reveal that (1) whole-brain temporal dynamics of dFC movie-watching task is closer to resting-state than to SMT with an overall trend of highest temporal stability observed during SMT followed by movie-watching and resting-state, invariant across lifespan aging, (2) in both tasks conditions stability of neurocognitive networks in young adults is higher than older adults, and (3) temporal stability of whole brain resting-state follows a U-shaped curve along lifespan-a pattern shared by sensorimotor network stability indicating their deeper relationship. Overall, the results can be applied generally for studying cohorts of neurological disorders using neuroimaging tools.

White-matter degradation and dynamical compensation support age-related functional alterations in human brain

Structural connectivity of the brain at different ages is analyzed using diffusion-weighted magnetic resonance imaging (MRI) data. The largest decrease of streamlines is found in frontal regions and for long inter-hemispheric links. The average length of the tracts also decreases, but the clustering is unaffected. From functional MRI we identify age-related changes of dynamic functional connectivity (dFC) and spatial covariation features of functional connectivity (FC) links captured by metaconnectivity. They indicate more stable dFC, but wider range and variance of MC, whereas static features of FC did not show any significant differences with age. We implement individual connectivity in whole-brain models and test several hypotheses for the mechanisms of operation among underlying neural system. We demonstrate that age-related functional fingerprints are only supported if the model accounts for: (i) compensation of the individual brains for the overall loss of structural connectivity and (ii) decrease of propagation velocity due to the loss of myelination. We also show that with these 2 conditions, it is sufficient to decompose the time-delays as bimodal distribution that only distinguishes between intra- and inter-hemispheric delays, and that the same working point also captures the static FC the best, and produces the largest variability at slow time-scales.

Dynamic sensitivity analysis: Defining personalised strategies to drive brain state transitions via whole brain modelling

Traditionally, in neuroimaging, model-free analyses are used to find significant differences between brain states via signal detection theory. Depending on the a priori assumptions about the underlying data, different spatio-temporal features can be analysed. Alternatively, model-based techniques infer features from the data and compare significance from model parameters. However, to assess transitions from one brain state to another remains a challenge in current paradigms. Here, we introduce a "Dynamic Sensitivity Analysis" framework that quantifies transitions between brain states in terms of stimulation ability to rebalance spatio-temporal brain activity towards a target state such as healthy brain dynamics. In practice, it means building a whole-brain model fitted to the spatio-temporal description of brain dynamics, and applying systematic stimulations in-silico to assess the optimal strategy to drive brain dynamics towards a target state. Further, we show how Dynamic Sensitivity Analysis extends to various brain stimulation paradigms, ultimately contributing to improving the efficacy of personalised clinical interventions.

First few seconds for flow: A comprehensive proposal of the neurobiology and neurodynamics of state onset

Flow is a cognitive state that manifests when there is complete attentional absorption while performing a task. Flow occurs when certain internal as well as external conditions are present, including intense concentration, a sense of control, feedback, and a balance between the challenge of the task and the relevant skillset. Phenomenologically, flow is accompanied by a loss of self-consciousness, seamless integration of action and awareness, and acute changes in time perception. Research has begun to uncover some of the neurophysiological correlates of flow, as well as some of the state's neuromodulatory processes. We comprehensively review this work and consider the neurodynamics of the onset of the state, considering large-scale brain networks, as well as dopaminergic, noradrenergic, and endocannabinoid systems. To accomplish this, we outline an evidence-based hypothetical situation, and consider the flow state in a broader context including other profound alterations in consciousness, such as the psychedelic state and the state of traumatic stress that can induce PTSD. We present a broad theoretical framework which may motivate future testable hypotheses.

Whole-Brain Propagation Delays in Multiple Sclerosis, a Combined Tractography-Magnetoencephalography Study

Two structurally connected brain regions are more likely to interact, with the lengths of the structural bundles, their widths, myelination, and the topology of the structural connectome influencing the timing of the interactions. We introduce an in vivo approach for measuring functional delays across the whole brain in humans (of either sex) using magneto/electroencephalography (MEG/EEG) and integrating them with the structural bundles. The resulting topochronic map of the functional delays/velocities shows that larger bundles have faster velocities. We estimated the topochronic map in multiple sclerosis patients, who have damaged myelin sheaths, and controls, demonstrating greater delays in patients across the network and that structurally lesioned tracts were slowed down more than unaffected ones. We provide a novel framework for estimating functional transmission delays in vivo at the single-subject and single-tract level.SIGNIFICANCE STATEMENT This article provides a straightforward way to estimate patient-specific delays and conduction velocities in the CNS, at the individual level, in healthy and diseased subjects. To do so, it uses a principled way to merge magnetoencephalography (MEG)/electroencephalography (EEG) and tractography.

Strength-dependent perturbation of whole-brain model working in different regimes reveals the role of fluctuations in brain dynamics

Despite decades of research, there is still a lack of understanding of the role and generating mechanisms of the ubiquitous fluctuations and oscillations found in recordings of brain dynamics. Here, we used whole-brain computational models capable of presenting different dynamical regimes to reproduce empirical data's turbulence level. We showed that the model's fluctuations regime fitted to turbulence more faithfully reproduces the empirical functional connectivity compared to oscillatory and noise regimes. By applying global and local strength-dependent perturbations and subsequently measuring the responsiveness of the model, we revealed each regime's computational capacity demonstrating that brain dynamics is shifted towards fluctuations to provide much-needed flexibility. Importantly, fluctuation regime stimulation in a brain region within a given resting state network modulates that network, aligned with previous empirical and computational studies. Furthermore, this framework generates specific, testable empirical predictions for human stimulation studies using strength-dependent rather than constant perturbation. Overall, the whole-brain models fitted to the level of empirical turbulence together with functional connectivity unveil that the fluctuation regime best captures empirical data, and the strength-dependent perturbative framework demonstrates how this regime provides maximal flexibility to the human brain.

Data-driven discovery of canonical large-scale brain dynamics

Human behavior and cognitive function correlate with complex patterns of spatio-temporal brain dynamics, which can be simulated using computational models with different degrees of biophysical realism. We used a data-driven optimization algorithm to determine and classify the types of local dynamics that enable the reproduction of different observables derived from functional magnetic resonance recordings. The phase space analysis of the resulting equations revealed a predominance of stable spiral attractors, which optimized the similarity to the empirical data in terms of the synchronization, metastability, and functional connectivity dynamics. For stable limit cycles, departures from harmonic oscillations improved the fit in terms of functional connectivity dynamics. Eigenvalue analyses showed that proximity to a bifurcation improved the accuracy of the simulation for wakefulness, whereas deep sleep was associated with increased stability. Our results provide testable predictions that constrain the landscape of suitable biophysical models, while supporting noise-driven dynamics close to a bifurcation as a canonical mechanism underlying the complex fluctuations that characterize endogenous brain activity.

Resting-state BOLD signal variability is associated with individual differences in metacontrol

Numerous studies demonstrate that moment-to-moment neural variability is behaviorally relevant and beneficial for tasks and behaviors requiring cognitive flexibility. However, it remains unclear whether the positive effect of neural variability also holds for cognitive persistence. Moreover, different brain variability measures have been used in previous studies, yet comparisons between them are lacking. In the current study, we examined the association between resting-state BOLD signal variability and two metacontrol policies (i.e., persistence vs. flexibility). Brain variability was estimated from resting-state fMRI (rsfMRI) data using two different approaches (i.e., Standard Deviation (SD), and Mean Square Successive Difference (MSSD)) and metacontrol biases were assessed by three metacontrol-sensitive tasks. Results showed that brain variability measured by SD and MSSD was highly positively related. Critically, higher variability measured by MSSD in the attention network, parietal and frontal network, frontal and ACC network, parietal and motor network, and higher variability measured by SD in the parietal and motor network, parietal and frontal network were associated with reduced persistence (or greater flexibility) of metacontrol (i.e., larger Stroop effect or worse RAT performance). These results show that the beneficial effect of brain signal variability on cognitive control depends on the metacontrol states involved. Our study highlights the importance of temporal variability of rsfMRI activity in understanding the neural underpinnings of cognitive control.

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

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

Targeting disrupted rich-club network organization with neuroplasticity-based computerized cognitive remediation in major depressive disorder patients

Disrupted rich-club organization has been extensively studied in major depressive disorder (MDD) patients. Although data indicate that neuroplasticity-based computerized cognitive remediation (nCCR) can accelerate clinical responses in MDD patients, the mechanisms underlying its antidepressant efficacy are unknown. In this study, all MDD patients underwent two (baseline and week 4) neuropsychological assessments and DTI imaging. Additionally, 17 MDD patients did nCCR for 30 hours spread across 4 weeks. Rich-club organization was calculated with a graph-theoretical approach, and SC-FC coupling was explored. After 4 weeks of treatment, the number of rich-club connections, global efficiency, and SC-FC coupling strength increased significantly and were negatively associated with TMT-B scores. The effects of nCCR on disrupted rich-club organization may partly underlie its efficacy in improving the executive function of patients with MDD. Effects of nCCR on disrupted rich-club organization may partly underlie its efficacy in improving the executive function of patients with MDD.

It's about time: Linking dynamical systems with human neuroimaging to understand the brain

Most human neuroscience research to date has focused on statistical approaches that describe stationary patterns of localized neural activity or blood flow. While these patterns are often interpreted in light of dynamic, information-processing concepts, the static, local, and inferential nature of the statistical approach makes it challenging to directly link neuroimaging results to plausible underlying neural mechanisms. Here, we argue that dynamical systems theory provides the crucial mechanistic framework for characterizing both the brain's time-varying quality and its partial stability in the face of perturbations, and hence, that this perspective can have a profound impact on the interpretation of human neuroimaging results and their relationship with behavior. After briefly reviewing some key terminology, we identify three key ways in which neuroimaging analyses can embrace a dynamical systems perspective: by shifting from a local to a more global perspective, by focusing on dynamics instead of static snapshots of neural activity, and by embracing modeling approaches that map neural dynamics using "forward" models. Through this approach, we envisage ample opportunities for neuroimaging researchers to enrich their understanding of the dynamic neural mechanisms that support a wide array of brain functions, both in health and in the setting of psychopathology.

Cross-attractor repertoire provides new perspective on structure-function relationship in the brain

The brain exhibits complex intrinsic dynamics, i.e., spontaneously arising activity patterns without any external inputs or tasks. Such intrinsic dynamics and their alteration are thought to play crucial roles in typical as well as atypical cognitive functioning. Linking the ever-changing intrinsic dynamics to the rather static anatomy is a challenging endeavor. Dynamical systems models are important tools for understanding how structure and function are linked in the brain. Here, we provide a novel modeling framework to examine how functional connectivity depends on structural connectivity in the brain. Existing modeling frameworks typically focus on noise-driven (or stochastic) dynamics near a single attractor. Complementing existing approaches, we examine deterministic features of the distribution of attractors, in particular, how regional states are correlated across all attractors - cross-attractor coordination. We found that cross-attractor coordination between brain regions better predicts human functional connectivity than noise-driven single-attractor dynamics. Importantly, cross-attractor coordination better accounts for the nonlinear dependency of functional connectivity on structural connectivity. Our findings suggest that functional connectivity patterns in the brain may reflect transitions between attractors, which impose an energy cost. The framework may be used to predict transitions and energy costs associated with experimental or clinical interventions.

Fully Bayesian estimation of virtual brain parameters with self-tuning Hamiltonian Monte Carlo

Virtual brain models are data-driven patient-specific brain models integrating individual brain imaging data with neural mass modeling in a single computational framework, capable of autonomously generating brain activity and its associated brain imaging signals. Along the example of epilepsy, we develop an efficient and accurate Bayesian methodology estimating the parameters linked to the extent of the epileptogenic zone. State-of-the-art advances in Bayesian inference using Hamiltonian Monte Carlo (HMC) algorithms have remained elusive for large-scale differential-equations based models due to their slow convergence. We propose appropriate priors and a novel reparameterization to facilitate efficient exploration of the posterior distribution in terms of computational time and convergence diagnostics. The methodology is illustrated for in-silico dataset and then, applied to infer the personalized model parameters based on the empirical stereotactic electroencephalography (SEEG) recordings of retrospective patients. This improved methodology may pave the way to render HMC methods sufficiently easy and efficient to use, thus applicable in personalized medicine.

From Brain-Body Function to Conscious Interactions

In this review, we discuss empirical results inspiring the introduction of a formal mathematical multilayer model for the biological neuroscience of conscious experience. First, we motivate the discussion through evidence regarding the dynamic brain. Second, we review different brain-body couplings associated with conscious experience and its potential role in driving brain dynamics. Third, we introduce the machinery of multilayer networks to account for several types of interactions in brain-body systems. Then, a multilayer structure consists of two main generalizations: a formal semantic to study biological systems, and an integrative account for several signatures and models of consciousness. Finally, under this framework, we define composition of layers to account for entangled features of brain-body systems related to conscious experience. As such, a multilayer mathematical framework is highly integrative and thus may be more complete than other models. In this short review, we discuss a variety of empirical results inspiring the introduction of a formal mathematical multilayer model for the biological neuroscience of conscious experience.

Reliability and subject specificity of personalized whole-brain dynamical models

Dynamical whole-brain models were developed to link structural (SC) and functional connectivity (FC) together into one framework. Nowadays, they are used to investigate the dynamical regimes of the brain and how these relate to behavioral, clinical and demographic traits. However, there is no comprehensive investigation on how reliable and subject specific the modeling results are given the variability of the empirical FC. In this study, we show that the parameters of these models can be fitted with a "poor" to "good" reliability depending on the exact implementation of the modeling paradigm. We find, as a general rule of thumb, that enhanced model personalization leads to increasingly reliable model parameters. In addition, we observe no clear effect of the model complexity evaluated by separately sampling results for linear, phase oscillator and neural mass network models. In fact, the most complex neural mass model often yields modeling results with "poor" reliability comparable to the simple linear model, but demonstrates an enhanced subject specificity of the model similarity maps. Subsequently, we show that the FC simulated by these models can outperform the empirical FC in terms of both reliability and subject specificity. For the structure-function relationship, simulated FC of individual subjects may be identified from the correlations with the empirical SC with an accuracy up to 70%, but not vice versa for non-linear models. We sample all our findings for 8 distinct brain parcellations and 6 modeling conditions and show that the parcellation-induced effect is much more pronounced for the modeling results than for the empirical data. In sum, this study provides an exploratory account on the reliability and subject specificity of dynamical whole-brain models and may be relevant for their further development and application. In particular, our findings suggest that the application of the dynamical whole-brain modeling should be tightly connected with an estimate of the reliability of the results.

Human brain structural connectivity matrices-ready for modelling

The human brain represents a complex computational system, the function and structure of which may be measured using various neuroimaging techniques focusing on separate properties of the brain tissue and activity. We capture the organization of white matter fibers acquired by diffusion-weighted imaging using probabilistic diffusion tractography. By segmenting the results of tractography into larger anatomical units, it is possible to draw inferences about the structural relationships between these parts of the system. This pipeline results in a structural connectivity matrix, which contains an estimate of connection strength among all regions. However, raw data processing is complex, computationally intensive, and requires expert quality control, which may be discouraging for researchers with less experience in the field. We thus provide brain structural connectivity matrices in a form ready for modelling and analysis and thus usable by a wide community of scientists. The presented dataset contains brain structural connectivity matrices together with the underlying raw diffusion and structural data, as well as basic demographic data of 88 healthy subjects.

Measurement(s)	Diffusion Weighted Imaging
Technology Type(s)	Magnetic Resonance Imaging
Sample Characteristic - Organism	Homo sapiens
Sample Characteristic - Location	Czech Republic

Normalizing the brain connectome for communication through synchronization

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

Unifying turbulent dynamics framework distinguishes different brain states

Significant advances have been made by identifying the levels of synchrony of the underlying dynamics of a given brain state. This research has demonstrated that non-conscious dynamics tend to be more synchronous than in conscious states, which are more asynchronous. Here we go beyond this dichotomy to demonstrate that different brain states are underpinned by dissociable spatiotemporal dynamics. We investigated human neuroimaging data from different brain states (resting state, meditation, deep sleep and disorders of consciousness after coma). The model-free approach was based on Kuramoto's turbulence framework using coupled oscillators. This was extended by a measure of the information cascade across spatial scales. Complementarily, the model-based approach used exhaustive in silico perturbations of whole-brain models fitted to these measures. This allowed studying of the information encoding capabilities in given brain states. Overall, this framework demonstrates that elements from turbulence theory provide excellent tools for describing and differentiating between brain states.

The modulatory effect of adaptive task-switching training on resting-state neural network dynamics in younger and older adults

With increasing life expectancy and active aging, it becomes crucial to investigate methods which could compensate for generally detected cognitive aging processes. A promising candidate is adaptive cognitive training, during which task difficulty is adjusted to the participants' performance level to enhance the training and potential transfer effects. Measuring intrinsic brain activity is suitable for detecting possible distributed training-effects since resting-state dynamics are linked to the brain's functional flexibility and the effectiveness of different cognitive processes. Therefore, we investigated if adaptive task-switching training could modulate resting-state neural dynamics in younger (18-25 years) and older (60-75 years) adults (79 people altogether). We examined spectral power density on resting-state EEG data for measuring oscillatory activity, and multiscale entropy for detecting intrinsic neural complexity. Decreased coarse timescale entropy and lower frequency band power as well as increased fine timescale entropy and higher frequency band power revealed a shift from more global to local information processing with aging before training. However, cognitive training modulated these age-group differences, as coarse timescale entropy and lower frequency band power increased from pre- to post-training in the old-training group. Overall, our results suggest that cognitive training can modulate neural dynamics even when measured outside of the trained task.

Unsupervised classification reveals consistency and degeneracy in neural network patterns of emotion

In the present study, we used an unsupervised classification algorithm to reveal both consistency and degeneracy in neural network connectivity during anger and anxiety. Degeneracy refers to the ability of different biological pathways to produce the same outcomes. Previous research is suggestive of degeneracy in emotion, but little research has explicitly examined whether degenerate functional connectivity patterns exist for emotion categories such as anger and anxiety. Twenty-four subjects underwent functional magnetic resonance imaging (fMRI) while listening to unpleasant music and self-generating experiences of anger and anxiety. A data-driven model building algorithm with unsupervised classification (subgrouping Group Iterative Multiple Model Estimation) identified patterns of connectivity among 11 intrinsic networks that were associated with anger vs anxiety. As predicted, degenerate functional connectivity patterns existed within these overarching consistent patterns. Degenerate patterns were not attributable to differences in emotional experience or other individual-level factors. These findings are consistent with the constructionist account that emotions emerge from flexible functional neuronal assemblies and that emotion categories such as anger and anxiety each describe populations of highly variable instances.

Towards an efficient validation of dynamical whole-brain models

Simulating the resting-state brain dynamics via mathematical whole-brain models requires an optimal selection of parameters, which determine the model’s capability to replicate empirical data. Since the parameter optimization via a grid search (GS) becomes unfeasible for high-dimensional models, we evaluate several alternative approaches to maximize the correspondence between simulated and empirical functional connectivity. A dense GS serves as a benchmark to assess the performance of four optimization schemes: Nelder-Mead Algorithm (NMA), Particle Swarm Optimization (PSO), Covariance Matrix Adaptation Evolution Strategy (CMAES) and Bayesian Optimization (BO). To compare them, we employ an ensemble of coupled phase oscillators built upon individual empirical structural connectivity of 105 healthy subjects. We determine optimal model parameters from two- and three-dimensional parameter spaces and show that the overall fitting quality of the tested methods can compete with the GS. There are, however, marked differences in the required computational resources and stability properties, which we also investigate before proposing CMAES and BO as efficient alternatives to a high-dimensional GS. For the three-dimensional case, these methods generated similar results as the GS, but within less than 6% of the computation time. Our results contribute to an efficient validation of models for personalized simulations of brain dynamics.

Dynamic primitives of brain network interaction Special Issue "Advances in Mapping the Connectome"

What dynamic processes underly functional brain networks? Functional connectivity (FC) and functional connectivity dynamics (FCD) are used to represent the patterns and dynamics of functional brain networks. FC(D) is related to the synchrony of brain activity: when brain areas oscillate in a coordinated manner this yields a high correlation between their signal time series. To explain the processes underlying FC(D) we review how synchronized oscillations emerge from coupled neural populations in brain network models (BNMs). From detailed spiking networks to more abstract population models, there is strong support for the idea that the brain operates near critical instabilities that give rise to multistable or metastable dynamics that in turn lead to the intermittently synchronized slow oscillations underlying FC(D). We explore further consequences from these fundamental mechanisms and how they fit with reality. We conclude by highlighting the need for integrative brain models that connect separate mechanisms across levels of description and spatiotemporal scales and link them with cognitive function.

Directed Brain Connectivity Identifies Widespread Functional Network Abnormalities in Parkinson's Disease

Parkinson's disease (PD) is a neurodegenerative disorder characterized by topological abnormalities in large-scale functional brain networks, which are commonly analyzed using undirected correlations in the activation signals between brain regions. This approach assumes simultaneous activation of brain regions, despite previous evidence showing that brain activation entails causality, with signals being typically generated in one region and then propagated to other ones. To address this limitation, here, we developed a new method to assess whole-brain directed functional connectivity in participants with PD and healthy controls using antisymmetric delayed correlations, which capture better this underlying causality. Our results show that whole-brain directed connectivity, computed on functional magnetic resonance imaging data, identifies widespread differences in the functional networks of PD participants compared with controls, in contrast to undirected methods. These differences are characterized by increased global efficiency, clustering, and transitivity combined with lower modularity. Moreover, directed connectivity patterns in the precuneus, thalamus, and cerebellum were associated with motor, executive, and memory deficits in PD participants. Altogether, these findings suggest that directional brain connectivity is more sensitive to functional network differences occurring in PD compared with standard methods, opening new opportunities for brain connectivity analysis and development of new markers to track PD progression.