Dynamical Mean Field approximation of a canonical cortical model for studying inter-population synchrony

The goal of this paper is twofold. We propose and explore a model to study the synchronization among populations in the canonical model of the neocortex proposed previously by (R.J. Douglas, K.A.C. Martin, A functional microcircuit for cat visual cortex. J.Physiol. 440(1991) 735–769). For this, a model describing N synapses of each m-population (m = 1, 2,3) is proposed. Each synapse is described by a system of 2 stochastic differential equations (SDEs). Then, by using the dynamical mean field approximation (DMA) (H. Hasegawa, Dynamical mean-field theory of spiking neuron ensembles: Response to a single spike with independent noises, Phys. Rev. E. (2003)1-19.) the system of several SDEs is reduced to 12 ordinary differential equations for the means and the second-order moments of global variables. The connectivity among populations is obtained by summarizing in the canonical model the detailed information from a quantitative description of the circuits formed in cat area 17 given in (T.Binzegger, R.J. Douglas, K.A. Martin, A Quantitative Map of the Circuit of Cat Primary Visual Cortex, J. Neurosci. 24 (2004) 8441- 8453). In the framework of the used DMA we propose a measure for inter-population synchronization. Simulations are carried out for exploring how inter-population synchrony is related to the variation of firing frequency of each population. Our results suggest that superficial pyramidal clusters appear to have a predominant influence on the synchronization process among pyramidal populations as well as put forward the active role of inhibition in the rest of the synchronizations between populations

Bayesian Comparison of Neurovascular Coupling Models Using EEG-fMRI

Functional magnetic resonance imaging (fMRI), with blood oxygenation level-dependent (BOLD) contrast, is a widely used technique for studying the human brain. However, it is an indirect measure of underlying neuronal activity and the processes that link this activity to BOLD signals are still a topic of much debate. In order to relate findings from fMRI research to other measures of neuronal activity it is vital to understand the underlying neurovascular coupling mechanism. Currently, there is no consensus on the relative roles of synaptic and spiking activity in the generation of the BOLD response. Here we designed a modelling framework to investigate different neurovascular coupling mechanisms. We use Electroencephalographic (EEG) and fMRI data from a visual stimulation task together with biophysically informed mathematical models describing how neuronal activity generates the BOLD signals. These models allow us to non-invasively infer the degree of local synaptic and spiking activity in the healthy human brain. In addition, we use Bayesian model comparison to decide between neurovascular coupling mechanisms. We show that the BOLD signal is dependent upon both the synaptic and spiking activity but that the relative contributions of these two inputs are dependent upon the underlying neuronal firing rate. When the underlying neuronal firing is low then the BOLD response is best explained by synaptic activity. However, when the neuronal firing rate is high then both synaptic and spiking activity are required to explain the BOLD signal

MEG and fMRI Fusion for Non-Linear Estimation of Neural and BOLD Signal Changes

The combined analysis of magnetoencephalography (MEG)/electroencephalography and functional magnetic resonance imaging (fMRI) measurements can lead to improvement in the description of the dynamical and spatial properties of brain activity. In this paper we empirically demonstrate this improvement using simulated and recorded task related MEG and fMRI activity. Neural activity estimates were derived using a dynamic Bayesian network with continuous real valued parameters by means of a sequential Monte Carlo technique. In synthetic data, we show that MEG and fMRI fusion improves estimation of the indirectly observed neural activity and smooths tracking of the blood oxygenation level dependent (BOLD) response. In recordings of task related neural activity the combination of MEG and fMRI produces a result with greater signal-to-noise ratio, that confirms the expectation arising from the nature of the experiment. The highly non-linear model of the BOLD response poses a difficult inference problem for neural activity estimation; computational requirements are also high due to the time and space complexity. We show that joint analysis of the data improves the system's behavior by stabilizing the differential equations system and by requiring fewer computational resources

Neurovascular coupling: insights from multi-modal dynamic causal modelling of fMRI and MEG

This technical note presents a framework for investigating the underlying mechanisms of neurovascular coupling in the human brain using multi-modal magnetoencephalography (MEG) and functional magnetic resonance (fMRI) neuroimaging data. This amounts to estimating the evidence for several biologically informed models of neurovascular coupling using variational Bayesian methods and selecting the most plausible explanation using Bayesian model comparison. First, fMRI data is used to localise active neuronal sources. The coordinates of neuronal sources are then used as priors in the specification of a DCM for MEG, in order to estimate the underlying generators of the electrophysiological responses. The ensuing estimates of neuronal parameters are used to generate neuronal drive functions, which model the pre or post synaptic responses to each experimental condition in the fMRI paradigm. These functions form the input to a model of neurovascular coupling, the parameters of which are estimated from the fMRI data. This establishes a Bayesian fusion technique that characterises the BOLD response - asking, for example, whether instantaneous or delayed pre or post synaptic signals mediate haemodynamic responses. Bayesian model comparison is used to identify the most plausible hypotheses about the causes of the multimodal data. We illustrate this procedure by comparing a set of models of a single-subject auditory fMRI and MEG dataset. Our exemplar analysis suggests that the origin of the BOLD signal is mediated instantaneously by intrinsic neuronal dynamics and that neurovascular coupling mechanisms are region-specific. The code and example dataset associated with this technical note are available through the statistical parametric mapping (SPM) software package

Simultaneous EEG-fMRI Analysis with Application to Detection of Seizure Signal Sources

Non

Dynamic models of brain imaging data and their Bayesian inversion

This work is about understanding the dynamics of neuronal systems, in particular with respect to brain connectivity. It addresses complex neuronal systems by looking at neuronal interactions and their causal relations. These systems are characterized using a generic approach to dynamical system analysis of brain signals - dynamic causal modelling (DCM). DCM is a technique for inferring directed connectivity among brain regions, which distinguishes between a neuronal and an observation level. DCM is a natural extension of the convolution models used in the standard analysis of neuroimaging data. This thesis develops biologically constrained and plausible models, informed by anatomic and physiological principles. Within this framework, it uses mathematical formalisms of neural mass, mean-field and ensemble dynamic causal models as generative models for observed neuronal activity. These models allow for the evaluation of intrinsic neuronal connections and high-order statistics of neuronal states, using Bayesian estimation and inference. Critically it employs Bayesian model selection (BMS) to discover the best among several equally plausible models. In the first part of this thesis, a two-state DCM for functional magnetic resonance imaging (fMRI) is described, where each region can model selective changes in both extrinsic and intrinsic connectivity. The second part is concerned with how the sigmoid activation function of neural-mass models (NMM) can be understood in terms of the variance or dispersion of neuronal states. The third part presents a mean-field model (MFM) for neuronal dynamics as observed with magneto- and electroencephalographic data (M/EEG). In the final part, the MFM is used as a generative model in a DCM for M/EEG and compared to the NMM using Bayesian model selection

Quantifying the time course of visual object processing using ERPs: it's time to up the game

Hundreds of studies have investigated the early ERPs to faces and objects using scalp and intracranial recordings. The vast majority of these studies have used uncontrolled stimuli, inappropriate designs, peak measurements, poor figures, and poor inferential and descriptive group statistics. These problems, together with a tendency to discuss any effect p < 0.05 rather than to report effect sizes, have led to a research field very much qualitative in nature, despite its quantitative inspirations, and in which predictions do not go beyond condition A > condition B. Here we describe the main limitations of face and object ERP research and suggest alternative strategies to move forward. The problems plague intracranial and surface ERP studies, but also studies using more advanced techniques – e.g., source space analyses and measurements of network dynamics, as well as many behavioral, fMRI, TMS, and LFP studies. In essence, it is time to stop amassing binary results and start using single-trial analyses to build models of visual perception

The Importance of Cerebellar Connectivity on Simulated Brain Dynamics

The brain shows a complex multiscale organization that prevents a direct understanding of how structure, function and dynamics are correlated. To date, advances in neural modeling offer a unique opportunity for simulating global brain dynamics by embedding empirical data on different scales in a mathematical framework. The Virtual Brain (TVB) is an advanced data-driven model allowing to simulate brain dynamics starting from individual subjects’ structural and functional connectivity obtained, for example, from magnetic resonance imaging (MRI). The use of TVB has been limited so far to cerebral connectivity but here, for the first time, we have introduced cerebellar nodes and interconnecting tracts to demonstrate the impact of cerebro-cerebellar loops on brain dynamics. Indeed, the matching between the empirical and simulated functional connectome was significantly improved when including the cerebro-cerebellar loops. This positive result should be considered as a first step, since issues remain open about the best strategy to reconstruct effective structural connectivity and the nature of the neural mass or mean-field models generating local activity in the nodes. For example, signal processing is known to differ remarkably between cortical and cerebellar microcircuits. Tackling these challenges is expected to further improve the predictive power of functional brain activity simulations, using TVB or other similar tools, in explaining not just global brain dynamics but also the role of cerebellum in determining brain states in physiological conditions and in the numerous pathologies affecting the cerebro-cerebellar loop

A biophysical model of dynamic balancing of excitation and inhibition in fast oscillatory large-scale networks

Over long timescales, neuronal dynamics can be robust to quite large perturbations, such as changes in white matter connectivity and grey matter structure through processes including learning, aging, development and certain disease processes. One possible explanation is that robust dynamics are facilitated by homeostatic mechanisms that can dynamically rebalance brain networks. In this study, we simulate a cortical brain network using the Wilson-Cowan neural mass model with conduction delays and noise, and use inhibitory synaptic plasticity (ISP) to dynamically achieve a spatially local balance between excitation and inhibition. Using MEG data from 55 subjects we find that ISP enables us to simultaneously achieve high correlation with multiple measures of functional connectivity, including amplitude envelope correlation and phase locking. Further, we find that ISP successfully achieves local E/I balance, and can consistently predict the functional connectivity computed from real MEG data, for a much wider range of model parameters than is possible with a model without ISP