Evaluation of Resting Spatio-Temporal Dynamics of a Neural Mass Model Using Resting fMRI Connectivity and EEG Microstates

Relationships and representations of brain structures, connectivity, dynamics and functions

The review explores the complex interplay between brain structures and their associated functions, presenting a diversity of hierarchical models that enhances our understanding of these relationships. Central to this approach are structure-function flow diagrams, which offer a visual representation of how specific neuroanatomical structures are linked to their functional roles. These diagrams are instrumental in mapping the intricate connections between different brain regions, providing a clearer understanding of how functions emerge from the underlying neural architecture. The study details innovative attempts to develop new functional hierarchies that integrate structural and functional data. These efforts leverage recent advancements in neuroimaging techniques such as fMRI, EEG, MEG, and PET, as well as computational models that simulate neural dynamics. By combining these approaches, the study seeks to create a more refined and dynamic hierarchy that can accommodate the brain's complexity, including its capacity for plasticity and adaptation. A significant focus is placed on the overlap of structures and functions within the brain. The manuscript acknowledges that many brain regions are multifunctional, contributing to different cognitive and behavioral processes depending on the context. This overlap highlights the need for a flexible, non-linear hierarchy that can capture the brain's intricate functional landscape. Moreover, the study examines the interdependence of these functions, emphasizing how the loss or impairment of one function can impact others. Another crucial aspect discussed is the brain's ability to compensate for functional deficits following neurological diseases or injuries. The investigation explores how the brain reorganizes itself, often through the recruitment of alternative neural pathways or the enhancement of existing ones, to maintain functionality despite structural damage. This compensatory mechanism underscores the brain's remarkable plasticity, demonstrating its ability to adapt and reconfigure itself in response to injury, thereby ensuring the continuation of essential functions. In conclusion, the study presents a system of brain functions that integrates structural, functional, and dynamic perspectives. It offers a robust framework for understanding how the brain's complex network of structures supports a wide range of cognitive and behavioral functions, with significant implications for both basic neuroscience and clinical applications.

Achieving Occam's razor: Deep learning for optimal model reduction

All fields of science depend on mathematical models. Occam's razor refers to the principle that good models should exclude parameters beyond those minimally required to describe the systems they represent. This is because redundancy can lead to incorrect estimates of model parameters from data, and thus inaccurate or ambiguous conclusions. Here, we show how deep learning can be powerfully leveraged to apply Occam's razor to model parameters. Our method, FixFit, uses a feedforward deep neural network with a bottleneck layer to characterize and predict the behavior of a given model from its input parameters. FixFit has three major benefits. First, it provides a metric to quantify the original model's degree of complexity. Second, it allows for the unique fitting of data. Third, it provides an unbiased way to discriminate between experimental hypotheses that add value versus those that do not. In three use cases, we demonstrate the broad applicability of this method across scientific domains. To validate the method using a known system, we apply FixFit to recover known composite parameters for the Kepler orbit model and a dynamic model of blood glucose regulation. In the latter, we demonstrate the ability to fit the latent parameters to real data. To illustrate how the method can be applied to less well-established fields, we use it to identify parameters for a multi-scale brain model and reduce the search space for viable candidate mechanisms.

Ketosis regulates K+ ion channels, strengthening brain-wide signaling disrupted by age

Aging is associated with impaired signaling between brain regions when measured using resting-state fMRI. This age-related destabilization and desynchronization of brain networks reverses itself when the brain switches from metabolizing glucose to ketones. Here, we probe the mechanistic basis for these effects. First, we confirmed their robustness across measurement modalities using two datasets acquired from resting-state EEG (Lifespan: standard diet, 20-80 years, N = 201; Metabolic: individually weight-dosed and calorically-matched glucose and ketone ester challenge,μ a g e = 26.9 ± 11.2 years , N = 36). Then, using a multiscale conductance-based neural mass model, we identified the unique set of mechanistic parameters consistent with our clinical data. Together, our results implicate potassium (K+) gradient dysregulation as a mechanism for age-related neural desynchronization and its reversal with ketosis, the latter finding of which is consistent with direct measurement of ion channels. As such, the approach facilitates the connection between macroscopic brain activity and cellular-level mechanisms.

Spatio-temporal dynamics of resting-state brain networks are associated with migraine disability

OBJECTIVE

The changes in resting-state functional networks and their correlations with clinical traits remain to be clarified in migraine. Here we aim to investigate the brain spatio-temporal dynamics of resting-state networks and their possible correlations with the clinical traits in migraine.

METHODS

Twenty Four migraine patients without aura and 26 healthy controls (HC) were enrolled. Each included subject underwent a resting-state EEG and echo planar imaging examination. The disability of migraine patients was evaluated by Migraine Disability Assessment (MIDAS). After data acquisition, EEG microstates (Ms) combining functional connectivity (FC) analysis based on Schafer 400-seven network atlas were performed. Then, the correlation between obtained parameters and clinical traits was investigated.

RESULTS

Compared with HC group, the brain temporal dynamics depicted by microstates showed significantly increased activity in functional networks involving MsB and decreased activity in functional networks involving MsD; The spatial dynamics were featured by decreased intra-network FC within the executive control network( ECN) and inter-network FC between dorsal attention network (DAN) and ECN (P < 0.05); Moreover, correlation analysis showed that the MIDAS score was positively correlated with the coverage and duration of MsC, and negatively correlated with the occurrence of MsA; The FC within default mode network (DMN), and the inter-FC of ECN- visual network (VN), ECN- limbic network, VN-limbic network was negatively correlated with MIDAS. However, the FC of DMN-ECN was positively correlated with MIDAS; Furthermore, significant interactions between the temporal and spatial dynamics were also obtained.

CONCLUSIONS

Our study confirmed the notion that altered spatio-temporal dynamics exist in migraine patients during resting-state. And the temporal dynamics, the spatial changes and the clinical traits such as migraine disability interact with each other. The spatio-temporal dynamics obtained from EEG microstate and fMRI FC analyses may be potential biomarkers for migraine and with a huge potential to change future clinical practice in migraine.

Ion gradient-driven bifurcations of a multi-scale neuronal model

Metabolic limitations within the brain frequently arise in the context of aging and disease. As the largest consumers of energy within the brain, ion pumps that maintain the neuronal membrane potential are the most affected when energy supply becomes limited. To characterize the effects of such limitations, we analyze the ion gradients present in a conductance-based (Morris-Lecar) neural mass model. We show the existence and locations of Neimark-Sacker and period-doubling bifurcations in the sodium, calcium, and potassium reversal potentials and demonstrate that these bifurcations form physiologically relevant bounds of ion gradient variability. Within these bounds, we show how depolarization of the gradients causes decreased neural activity. We also show that the depolarization of ion gradients decreases inter-regional coherence, causing a shift in the critical point at which the coupling occurs and thereby inducing loss of synchrony between regions. In this way, we show that the Larter-Breakspear model captures ion gradient variability present at the microscale level and propagates these changes to the macroscale effects such as those observed in human neuroimaging studies.

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.

Spatiotemporal Patterns of Adaptation-Induced Slow Oscillations in a Whole-Brain Model of Slow-Wave Sleep

During slow-wave sleep, the brain is in a self-organized regime in which slow oscillations (SOs) between up- and down-states travel across the cortex. While an isolated piece of cortex can produce SOs, the brain-wide propagation of these oscillations are thought to be mediated by the long-range axonal connections. We address the mechanism of how SOs emerge and recruit large parts of the brain using a whole-brain model constructed from empirical connectivity data in which SOs are induced independently in each brain area by a local adaptation mechanism. Using an evolutionary optimization approach, good fits to human resting-state fMRI data and sleep EEG data are found at values of the adaptation strength close to a bifurcation where the model produces a balance between local and global SOs with realistic spatiotemporal statistics. Local oscillations are more frequent, last shorter, and have a lower amplitude. Global oscillations spread as waves of silence across the undirected brain graph, traveling from anterior to posterior regions. These traveling waves are caused by heterogeneities in the brain network in which the connection strengths between brain areas determine which areas transition to a down-state first, and thus initiate traveling waves across the cortex. Our results demonstrate the utility of whole-brain models for explaining the origin of large-scale cortical oscillations and how they are shaped by the connectome.

Cortical traveling waves reflect state-dependent hierarchical sequencing of local regions in the human connectome network

Recent human studies using electrocorticography have demonstrated that alpha and theta band oscillations form traveling waves on the cortical surface. According to neural synchronization theories, the cortical traveling waves may group local cortical regions and sequence them by phase synchronization; however these contributions have not yet been assessed. This study aimed to evaluate the functional contributions of traveling waves using connectome-based network modeling. In the simulation, we observed stable traveling waves on the entire cortical surface wherein the topographical pattern of these phases was substantially correlated with the empirically obtained resting-state networks, and local radial waves also appeared within the size of the empirical networks (< 50 mm). Importantly, individual regions in the entire network were instantaneously sequenced by their internal frequencies, and regions with higher intrinsic frequency were seen in the earlier phases of the traveling waves. Based on the communication-through-coherence theory, this phase configuration produced a hierarchical organization of each region by unidirectional communication between the arbitrarily paired regions. In conclusion, cortical traveling waves reflect the intrinsic frequency-dependent hierarchical sequencing of local regions, global traveling waves sequence the set of large-scale cortical networks, and local traveling waves sequence local regions within individual cortical networks.

Resting state electroencephalography microstates in autism spectrum disorder: A mini-review

Atypical spatial organization and temporal characteristics, found via resting state electroencephalography (EEG) microstate analysis, have been associated with psychiatric disorders but these temporal and spatial parameters are less known in autism spectrum disorder (ASD). EEG microstates reflect a short time period of stable scalp potential topography. These canonical microstates (i.e., A, B, C, and D) and more are identified by their unique topographic map, mean duration, fraction of time covered, frequency of occurrence and global explained variance percentage; a measure of how well topographical maps represent EEG data. We reviewed the current literature for resting state microstate analysis in ASD and identified eight publications. This current review indicates there is significant alterations in microstate parameters in ASD populations as compared to typically developing (TD) populations. Microstate parameters were also found to change in relation to specific cognitive processes. However, as microstate parameters are found to be changed by cognitive states, the differently acquired data (e.g., eyes closed or open) resting state EEG are likely to produce disparate results. We also review the current understanding of EEG sources of microstates and the underlying brain networks.

neurolib: A Simulation Framework for Whole-Brain Neural Mass Modeling

neurolib is a computational framework for whole-brain modeling written in Python. It provides a set of neural mass models that represent the average activity of a brain region on a mesoscopic scale. In a whole-brain network model, brain regions are connected with each other based on biologically informed structural connectivity, i.e., the connectome of the brain. neurolib can load structural and functional datasets, set up a whole-brain model, manage its parameters, simulate it, and organize its outputs for later analysis. The activity of each brain region can be converted into a simulated BOLD signal in order to calibrate the model against empirical data from functional magnetic resonance imaging (fMRI). Extensive model analysis is made possible using a parameter exploration module, which allows one to characterize a model’s behavior as a function of changing parameters. An optimization module is provided for fitting models to multimodal empirical data using evolutionary algorithms. neurolib is designed to be extendable and allows for easy implementation of custom neural mass models, offering a versatile platform for computational neuroscientists for prototyping models, managing large numerical experiments, studying the structure–function relationship of brain networks, and for performing in-silico optimization of whole-brain models.

Tractography density affects whole-brain structural architecture and resting-state dynamical modeling

Dynamical modeling of the resting-state brain dynamics essentially relies on the empirical neuroimaging data utilized for the model derivation and validation. There is however still no standardized data processing for magnetic resonance imaging pipelines and the structural and functional connectomes involved in the models. In this study, we thus address how the parameters of diffusion-weighted data processing for structural connectivity (SC) can influence the validation results of the whole-brain mathematical models informed by SC. For this, we introduce a set of simulation conditions including the varying number of total streamlines of the whole-brain tractography (WBT) used for extraction of SC, cortical parcellations based on functional and anatomical brain properties and distinct model fitting modalities. The main objective of this study is to explore how the quality of the model validation can vary across the considered simulation conditions. We observed that the graph-theoretical network properties of structural connectome can be affected by varying tractography density and strongly relate to the model performance. We also found that the optimal number of the total streamlines of WBT can vary for different brain atlases. Consequently, we suggest a way how to improve the model performance based on the network properties and the optimal parameter configurations from multiple WBT conditions. Furthermore, the population of subjects can be stratified into subgroups with divergent behaviors induced by the varying WBT density such that different recommendations can be made with respect to the data processing for individual subjects and brain parcellations.

MEG current source reconstruction using a meta-analysis fMRI prior

Magnetoencephalography (MEG) offers a unique way to noninvasively investigate millisecond-order cortical activities by mapping sensor signals (magnetic fields outside the head) to cortical current sources using current source reconstruction methods. Current source reconstruction is defined as an ill-posed inverse problem, since the number of sensors is less than the number of current sources. One powerful approach to solving this problem is to use functional MRI (fMRI) data as a spatial constraint, although it boosts the cost of measurement and the burden on subjects. Here, we show how to use the meta-analysis fMRI data synthesized from thousands of papers instead of the individually recorded fMRI data. To mitigate the differences between the meta-analysis and individual data, the former are imported as prior information of the hierarchical Bayesian estimation. Using realistic simulations, we found out the performance of current source reconstruction using meta-analysis fMRI data to be better than that using low-quality individual fMRI data and conventional methods. By applying experimental data of a face recognition task, we qualitatively confirmed that group analysis results using the meta-analysis fMRI data showed a tendency similar to the results using the individual fMRI data. Our results indicate that the use of meta-analysis fMRI data improves current source reconstruction without additional measurement costs. We assume the proposed method would have greater effect for modalities with lower measurement costs, such as optically pumped magnetometers.