Comparing dynamic causal models of neurovascular coupling with fMRI and EEG/MEG

Contemporary neurocognitive models of memory: A descriptive comparative analysis

The great complexity involved in the study of memory has given rise to numerous hypotheses and models associated with various phenomena at different levels of analysis. This has allowed us to delve deeper in our knowledge about memory but has also made it difficult to synthesize and integrate data from different lines of research. In this context, this work presents a descriptive comparative analysis of contemporary models that address the structure and function of multiple memory systems. The main goal is to outline a panoramic view of the key elements that constitute these models in order to visualize both the current state of research and possible future directions. The elements that stand out from different levels of analysis are distributed neural networks, hierarchical organization, predictive coding, homeostasis, and evolutionary perspective.

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

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

Studying the Alzheimer's disease continuum using EEG and fMRI in single-modality and multi-modality settings

Alzheimer's disease (AD) is a biological, clinical continuum that covers the preclinical, prodromal, and clinical phases of the disease. Early diagnosis and identification of the stages of Alzheimer's disease (AD) are crucial in clinical practice. Ideally, biomarkers should reflect the underlying process (pathological or otherwise), be reproducible and non-invasive, and allow repeated measurements over time. However, the currently known biomarkers for AD are not suitable for differentiating the stages and predicting the trajectory of disease progression. Some objective parameters extracted using electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) are widely applied to diagnose the stages of the AD continuum. While electroencephalography (EEG) has a high temporal resolution, fMRI has a high spatial resolution. Combined EEG and fMRI (EEG-fMRI) can overcome single-modality drawbacks and obtain multi-dimensional information simultaneously, and it can help explore the hemodynamic changes associated with the neural oscillations that occur during information processing. This technique has been used in the cognitive field in recent years. This review focuses on the different techniques available for studying the AD continuum, including EEG and fMRI in single-modality and multi-modality settings, and the possible future directions of AD diagnosis using EEG-fMRI.

Low-intensity ultrasound stimulation modulates cortical neurovascular coupling in an attention deficit hyperactivity disorder rat model

Attention deficit hyperactivity disorder is accompanied by changes in cranial nerve function and cerebral blood flow (CBF). Low-intensity ultrasound stimulation can modulate brain neural activity in attention deficit hyperactivity disorder. However, to date, the modulatory effects of low-intensity ultrasound stimulation on CBF and neurovascular coupling in attention deficit hyperactivity disorder have not been reported. To address this question, Sprague-Dawley, Wistar-Kyoto, and spontaneously hypertensive (attention deficit hyperactivity disorder (ADHD) rat model) rats were divided into the control and low-intensity ultrasound stimulation (LIUS) groups. Cortical electrical stimulation was used to induce cortical excitability in different types of rats, and a penetrable laser speckle contrast imaging (LSCI) system and electrodes were used to evaluate the electrical stimulation-induced CBF, cortical excitability, and neurovascular coupling in free-moving rats. The CBF, cortical excitability, and neurovascular coupling (NVC) under cortical electrical stimulation in the attention deficit hyperactivity disorder rats were significantly different from those in the Sprague-Dawley and Wistar-Kyoto rats. We also found that low-intensity ultrasound stimulation significantly interfered with the cortical excitability and neurovascular coupling induced by cortical electrical stimulation in rats with attention deficit hyperactivity disorder. Our findings suggest that neurovascular coupling is a potential biomarker for attention deficit hyperactivity disorder. Furthermore, low-intensity ultrasound stimulation can improve abnormal brain function in attention deficit hyperactivity disorder and lay a research foundation for its application in the clinical treatment of attention deficit hyperactivity disorder.

Neuronal connected burst cascades bridge macroscale adaptive signatures across arousal states

The human brain displays a rich repertoire of states that emerge from the microscopic interactions of cortical and subcortical neurons. Difficulties inherent within large-scale simultaneous neuronal recording limit our ability to link biophysical processes at the microscale to emergent macroscopic brain states. Here we introduce a microscale biophysical network model of layer-5 pyramidal neurons that display graded coarse-sampled dynamics matching those observed in macroscale electrophysiological recordings from macaques and humans. We invert our model to identify the neuronal spike and burst dynamics that differentiate unconscious, dreaming, and awake arousal states and provide insights into their functional signatures. We further show that neuromodulatory arousal can mediate different modes of neuronal dynamics around a low-dimensional energy landscape, which in turn changes the response of the model to external stimuli. Our results highlight the promise of multiscale modelling to bridge theories of consciousness across spatiotemporal scales.

Multi-modal and multi-model interrogation of large-scale functional brain networks

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

Hi-Fi fMRI: High-resolution, fast-sampled and sub-second whole-brain functional MRI at 3T in humans

Functional magnetic resonance imaging (fMRI) is a methodological cornerstone of neuroscience. Most studies measure blood-oxygen-level-dependent (BOLD) signal using echo-planar imaging (EPI), Cartesian sampling, and image reconstruction with a one-to-one correspondence between the number of acquired volumes and reconstructed images. However, EPI schemes are subject to trade-offs between spatial and temporal resolutions. We overcome these limitations by measuring BOLD with a gradient recalled echo (GRE) with 3D radial-spiral phyllotaxis trajectory at a high sampling rate (28.24ms) on standard 3T field-strength. The framework enables the reconstruction of 3D signal time courses with whole-brain coverage at simultaneously higher spatial (1mm ³ ) and temporal (up to 250ms) resolutions, as compared to optimized EPI schemes. Additionally, artifacts are corrected before image reconstruction; the desired temporal resolution is chosen after scanning and without assumptions on the shape of the hemodynamic response. By showing activation in the calcarine sulcus of 20 participants performing an ON-OFF visual paradigm, we demonstrate the reliability of our method for cognitive neuroscience research.

Amortized Bayesian inference on generative dynamical network models of epilepsy using deep neural density estimators

Whole-brain modeling of epilepsy combines personalized anatomical data with dynamical models of abnormal activities to generate spatio-temporal seizure patterns as observed in brain imaging data. Such a parametric simulator is equipped with a stochastic generative process, which itself provides the basis for inference and prediction of the local and global brain dynamics affected by disorders. However, the calculation of likelihood function at whole-brain scale is often intractable. Thus, likelihood-free algorithms are required to efficiently estimate the parameters pertaining to the hypothetical areas, ideally including the uncertainty. In this study, we introduce the simulation-based inference for the virtual epileptic patient model (SBI-VEP), enabling us to amortize the approximate posterior of the generative process from a low-dimensional representation of whole-brain epileptic patterns. The state-of-the-art deep learning algorithms for conditional density estimation are used to readily retrieve the statistical relationships between parameters and observations through a sequence of invertible transformations. We show that the SBI-VEP is able to efficiently estimate the posterior distribution of parameters linked to the extent of the epileptogenic and propagation zones from sparse intracranial electroencephalography recordings. The presented Bayesian methodology can deal with non-linear latent dynamics and parameter degeneracy, paving the way for fast and reliable inference on brain disorders from neuroimaging modalities.

Dynamic causal modelling shows a prominent role of local inhibition in alpha power modulation in higher visual cortex

During resting-state EEG recordings, alpha activity is more prominent over the posterior cortex in eyes-closed (EC) conditions compared to eyes-open (EO). In this study, we characterized the difference in spectra between EO and EC conditions using dynamic causal modelling. Specifically, we investigated the role of intrinsic and extrinsic connectivity-within the visual cortex-in generating EC-EO alpha power differences over posterior electrodes. The primary visual cortex (V1) and the bilateral middle temporal visual areas (V5) were equipped with bidirectional extrinsic connections using a canonical microcircuit. The states of four intrinsically coupled subpopulations-within each occipital source-were also modelled. Using Bayesian model selection, we tested whether modulations of the intrinsic connections in V1, V5 or extrinsic connections (or a combination thereof) provided the best evidence for the data. In addition, using parametric empirical Bayes (PEB), we estimated group averages under the winning model. Bayesian model selection showed that the winning model contained both extrinsic connectivity modulations, as well as intrinsic connectivity modulations in all sources. The PEB analysis revealed increased extrinsic connectivity during EC. Overall, we found a reduction in the inhibitory intrinsic connections during EC. The results suggest that the intrinsic modulations in V5 played the most important role in producing EC-EO alpha differences, suggesting an intrinsic disinhibition in higher order visual cortex, during EC resting state.

Investigating the Influence of Different Umami Tastants on Brain Perception via Scalp Electroencephalogram

Three types of tastants are known as perceptually associated with umami taste: monosodium glutamate (MSG), disodium succinate (WSA), and disodium inosine monophosphate (IMP). While these tastants were confirmed to be perceptually similar in a sensory study, they could be discriminated (p < 0.05) by electroencephalogram (EEG) analysis on a time scale of 5-6 s. In comparison of the EEG responses of the participants, the brain could partly distinguish (p < 0.05) between different sensory intensities of MSG, WSA, or IMP. The EEG data indicated that the brain is partially sensitive to perceiving different sensory intensities (L, low; M, medium; and H, high) of the same umami stimuli; i.e., for MSG in μV²/Hz, L, 2.473 ± 0.181; M, 3.274 ± 0.181; and H, 3.202 ± 0.181. However, brain responses of perceptually equi-umami intensities could partially be discriminated, suggesting that the brain could partially discriminate (p < 0.05) MSG, WSA, and IMP, despite similar sensory intensities. Moreover, umami tastants were also found to significantly enhance (p < 0.05) the α wave activity, with the most responsive being at 10 Hz, particularly in the frontal and parietal and occipital regions of the brain (p < 0.001). This study shows the potential of EEG to investigate brain activity triggered by umami stimuli.

Auditory cortex modelled as a dynamical network of oscillators: understanding event-related fields and their adaptation

Adaptation, the reduction of neuronal responses by repetitive stimulation, is a ubiquitous feature of auditory cortex (AC). It is not clear what causes adaptation, but short-term synaptic depression (STSD) is a potential candidate for the underlying mechanism. In such a case, adaptation can be directly linked with the way AC produces context-sensitive responses such as mismatch negativity and stimulus-specific adaptation observed on the single-unit level. We examined this hypothesis via a computational model based on AC anatomy, which includes serially connected core, belt, and parabelt areas. The model replicates the event-related field (ERF) of the magnetoencephalogram as well as ERF adaptation. The model dynamics are described by excitatory and inhibitory state variables of cell populations, with the excitatory connections modulated by STSD. We analysed the system dynamics by linearising the firing rates and solving the STSD equation using time-scale separation. This allows for characterisation of AC dynamics as a superposition of damped harmonic oscillators, so-called normal modes. We show that repetition suppression of the N1m is due to a mixture of causes, with stimulus repetition modifying both the amplitudes and the frequencies of the normal modes. In this view, adaptation results from a complete reorganisation of AC dynamics rather than a reduction of activity in discrete sources. Further, both the network structure and the balance between excitation and inhibition contribute significantly to the rate with which AC recovers from adaptation. This lifetime of adaptation is longer in the belt and parabelt than in the core area, despite the time constants of STSD being spatially homogeneous. Finally, we critically evaluate the use of a single exponential function to describe recovery from adaptation.

Toward an integrative neurovascular framework for studying brain networks

Brain functional connectivity based on the measure of blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signals has become one of the most widely used measurements in human neuroimaging. However, the nature of the functional networks revealed by BOLD fMRI can be ambiguous, as highlighted by a recent series of experiments that have suggested that typical resting-state networks can be replicated from purely vascular or physiologically driven BOLD signals. After going through a brief review of the key concepts of brain network analysis, we explore how the vascular and neuronal systems interact to give rise to the brain functional networks measured with BOLD fMRI. This leads us to emphasize a view of the vascular network not only as a confounding element in fMRI but also as a functionally relevant system that is entangled with the neuronal network. To study the vascular and neuronal underpinnings of BOLD functional connectivity, we consider a combination of methodological avenues based on multiscale and multimodal optical imaging in mice, used in combination with computational models that allow the integration of vascular information to explain functional connectivity.

[Research progress and application of transfer entropy algorithm]

In recent years, exploring the physiological and pathological mechanisms of brain functional integration from the neural network level has become one of the focuses of neuroscience research. Due to the non-stationary and nonlinear characteristics of neural signals, its linear characteristics are not sufficient to fully explain the potential neurophysiological activity mechanism in the implementation of complex brain functions. In order to overcome the limitation that the linear algorithm cannot effectively analyze the nonlinear characteristics of signals, researchers proposed the transfer entropy (TE) algorithm. In recent years, with the introduction of the concept of brain functional network, TE has been continuously optimized as a powerful tool for nonlinear time series multivariate analysis. This paper first introduces the principle of TE algorithm and the research progress of related improved algorithms, discusses and compares their respective characteristics, and then summarizes the application of TE algorithm in the field of electrophysiological signal analysis. Finally, combined with the research progress in recent years, the existing problems of TE are discussed, and the future development direction is prospected.

Atypically high influence of subcortical activity on primary sensory regions in autism

•

The age-dependent decrease of subcortico-cortical connectivity is attenuated in ASD.

•

Primary sensory regions remain less segregated from subcortical activity in ASD.

•

This could underlie an excessive amount of sensory input relayed to the cortex.

Background

Hypersensitivity, stereotyped behaviors and attentional problems in autism spectrum disorder (ASD) are compatible with inefficient filtering of undesired or irrelevant sensory information at early stages of neural processing. This could stem from the persistent overconnectivity between primary sensory regions and deep brain nuclei in both children and adults with ASD – as reported by several previous studies – which could reflect a decreased or arrested maturation of brain connectivity. However, it has not yet been investigated whether this overconnectivity can be modelled as an excessive directional influence of subcortical brain activity on primary sensory cortical regions in ASD, with respect to age-matched typically developing (TD) individuals.

Methods

To this aim, we used dynamic causal modelling to estimate (1) the directional influence of subcortical activity on cortical processing and (2) the functional segregation of primary sensory cortical regions from subcortical activity in 166 participants with ASD and 193 TD participants from the Autism Brain Imaging Data Exchange (ABIDE). We then specifically tested the hypothesis that the age-related changes of these indicators of brain connectivity would differ between the two groups.

Results

We found that in TD participants age was significantly associated with decreased influence of subcortical activity on cortical processing, paralleled by an increased functional segregation of cortical sensory processing from subcortical activity. Instead these effects were highly reduced and mostly absent in ASD participants, suggesting a delayed or arrested development of the segregation between subcortical and cortical sensory processing in ASD.

Conclusion

This atypical configuration of subcortico-cortical connectivity in ASD can result in an excessive amount of unprocessed sensory input relayed to the cortex, which is likely to impact cognitive functioning in everyday situations where it is beneficial to limit the influence of basic sensory information on cognitive processing, such as activities requiring focused attention or social interactions.

Grey-box modeling and hypothesis testing of functional near-infrared spectroscopy-based cerebrovascular reactivity to anodal high-definition tDCS in healthy humans

Transcranial direct current stimulation (tDCS) has been shown to evoke hemodynamics response; however, the mechanisms have not been investigated systematically using systems biology approaches. Our study presents a grey-box linear model that was developed from a physiologically detailed multi-compartmental neurovascular unit model consisting of the vascular smooth muscle, perivascular space, synaptic space, and astrocyte glial cell. Then, model linearization was performed on the physiologically detailed nonlinear model to find appropriate complexity (Akaike information criterion) to fit functional near-infrared spectroscopy (fNIRS) based measure of blood volume changes, called cerebrovascular reactivity (CVR), to high-definition (HD) tDCS. The grey-box linear model was applied on the fNIRS-based CVR during the first 150 seconds of anodal HD-tDCS in eleven healthy humans. The grey-box linear models for each of the four nested pathways starting from tDCS scalp current density that perturbed synaptic potassium released from active neurons for Pathway 1, astrocytic transmembrane current for Pathway 2, perivascular potassium concentration for Pathway 3, and voltage-gated ion channel current on the smooth muscle cell for Pathway 4 were fitted to the total hemoglobin concentration (tHb) changes from optodes in the vicinity of 4x1 HD-tDCS electrodes as well as on the contralateral sensorimotor cortex. We found that the tDCS perturbation Pathway 3 presented the least mean square error (MSE, median <2.5%) and the lowest Akaike information criterion (AIC, median -1.726) from the individual grey-box linear model fitting at the targeted-region. Then, minimal realization transfer function with reduced-order approximations of the grey-box model pathways was fitted to the ensemble average tHb time series. Again, Pathway 3 with nine poles and two zeros (all free parameters), provided the best Goodness of Fit of 0.0078 for Chi-Square difference test of nested pathways. Therefore, our study provided a systems biology approach to investigate the initial transient hemodynamic response to tDCS based on fNIRS tHb data. Future studies need to investigate the steady-state responses, including steady-state oscillations found to be driven by calcium dynamics, where transcranial alternating current stimulation may provide frequency-dependent physiological entrainment for system identification. We postulate that such a mechanistic understanding from system identification of the hemodynamics response to transcranial electrical stimulation can facilitate adequate delivery of the current density to the neurovascular tissue under simultaneous portable imaging in various cerebrovascular diseases.

Adiabatic dynamic causal modelling

This technical note introduces adiabatic dynamic causal modelling, a method for inferring slow changes in biophysical parameters that control fluctuations of fast neuronal states. The application domain we have in mind is inferring slow changes in variables (e.g., extracellular ion concentrations or synaptic efficacy) that underlie phase transitions in brain activity (e.g., paroxysmal seizure activity). The scheme is efficient and yet retains a biophysical interpretation, in virtue of being based on established neural mass models that are equipped with a slow dynamic on the parameters (such as synaptic rate constants or effective connectivity). In brief, we use an adiabatic approximation to summarise fast fluctuations in hidden neuronal states (and their expression in sensors) in terms of their second order statistics; namely, their complex cross spectra. This allows one to specify and compare models of slowly changing parameters (using Bayesian model reduction) that generate a sequence of empirical cross spectra of electrophysiological recordings. Crucially, we use the slow fluctuations in the spectral power of neuronal activity as empirical priors on changes in synaptic parameters. This introduces a circular causality, in which synaptic parameters underwrite fast neuronal activity that, in turn, induces activity-dependent plasticity in synaptic parameters. In this foundational paper, we describe the underlying model, establish its face validity using simulations and provide an illustrative application to a chemoconvulsant animal model of seizure activity.

Ageing and the Ipsilateral M1 BOLD Response: A Connectivity Study

Young people exhibit a negative BOLD response in ipsilateral primary motor cortex (M1) when making unilateral movements, such as button presses. This negative BOLD response becomes more positive as people age. In this study, we investigated why this occurs, in terms of the underlying effective connectivity and haemodynamics. We applied dynamic causal modeling (DCM) to task fMRI data from 635 participants aged 18-88 from the Cam-CAN dataset, who performed a cued button pressing task with their right hand. We found that connectivity from contralateral supplementary motor area (SMA) and dorsal premotor cortex (PMd) to ipsilateral M1 became more positive with age, explaining 44% of the variability across people in ipsilateral M1 responses. In contrast, connectivity from contralateral M1 to ipsilateral M1 was weaker and did not correlate with individual differences in rM1 BOLD. Neurovascular and haemodynamic parameters in the model were not able to explain the age-related shift to positive BOLD. Our results add to a body of evidence implicating neural, rather than vascular factors as the predominant cause of negative BOLD-while emphasising the importance of inter-hemispheric connectivity. This study provides a foundation for investigating the clinical and lifestyle factors that determine the sign and amplitude of the M1 BOLD response in ageing, which could serve as a proxy for neural and vascular health, via the underlying neurovascular mechanisms.

The effects of age on resting-state BOLD signal variability is explained by cardiovascular and cerebrovascular factors

Accurate identification of brain function is necessary to understand neurocognitive aging, and thereby promote health and well-being. Many studies of neurocognitive aging have investigated brain function with the blood-oxygen level-dependent (BOLD) signal measured by functional magnetic resonance imaging. However, the BOLD signal is a composite of neural and vascular signals, which are differentially affected by aging. It is, therefore, essential to distinguish the age effects on vascular versus neural function. The BOLD signal variability at rest (known as resting state fluctuation amplitude, RSFA), is a safe, scalable, and robust means to calibrate vascular responsivity, as an alternative to breath-holding and hypercapnia. However, the use of RSFA for normalization of BOLD imaging assumes that age differences in RSFA reflecting only vascular factors, rather than age-related differences in neural function (activity) or neuronal loss (atrophy). Previous studies indicate that two vascular factors, cardiovascular health (CVH) and cerebrovascular function, are insufficient when used alone to fully explain age-related differences in RSFA. It remains possible that their joint consideration is required to fully capture age differences in RSFA. We tested the hypothesis that RSFA no longer varies with age after adjusting for a combination of cardiovascular and cerebrovascular measures. We also tested the hypothesis that RSFA variation with age is not associated with atrophy. We used data from the population-based, lifespan Cam-CAN cohort. After controlling for cardiovascular and cerebrovascular estimates alone, the residual variance in RSFA across individuals was significantly associated with age. However, when controlling for both cardiovascular and cerebrovascular estimates, the variance in RSFA was no longer associated with age. Grey matter volumes did not explain age differences in RSFA, after controlling for CVH. The results were consistent between voxel-level analysis and independent component analysis. Our findings indicate that cardiovascular and cerebrovascular signals are together sufficient predictors of age differences in RSFA. We suggest that RSFA can be used to separate vascular from neuronal factors, to characterize neurocognitive aging. We discuss the implications and make recommendations for the use of RSFA in the research of aging.

MRIES: A Matlab Toolbox for Mapping the Responses to Intracranial Electrical Stimulation

Propose

Directed cortical responses to intracranial electrical stimulation are a good standard for mapping inter-regional direct connectivity. Cortico-cortical evoked potential (CCEP), elicited by single pulse electrical stimulation (SPES), has been widely used to map the normal and abnormal brain effective network. However, automated processing of CCEP datasets and visualization of connectivity results remain challenging for researchers and clinicians. In this study, we develop a Matlab toolbox named MRIES (Mapping the Responses to Intracranial Electrical Stimulation) to automatically process CCEP data and visualize the connectivity results.

Method

The MRIES integrates the processing pipeline of the CCEP datasets and various methods for connectivity calculation based on low- and high-frequency signals with stimulation artifacts removed. The connectivity matrices are saved in different folders for visualization. Different visualization patterns (connectivity matrix, circle map, surface map, and volume map) are also integrated to the graphical user interface (GUI), which makes it easy to intuitively display and compare different connectivity measurements. Furthermore, one sample CCEP data set collected from eight epilepsy patients is used to validate the MRIES toolbox.

Result

We show the GUI and visualization functions of MRIES using one example CCEP data that has been described in a complete tutorial. We applied this toolbox to the sample CCEP data set to investigate the direct connectivity between the medial temporal lobe and the insular cortex. We find bidirectional connectivity between MTL and insular that are consistent with the findings of previous studies.

Conclusion

MRIES has a friendly GUI and integrates the full processing pipeline of CCEP data and various visualization methods. The MRIES toolbox, tutorial, and example data can be freely downloaded. As an open-source package, MRIES is expected to improve the reproducibility of CCEP findings and facilitate clinical translation.

Markov blankets in the brain

Recent characterisations of self-organising systems depend upon the presence of a 'Markov blanket': a statistical boundary that mediates the interactions between the inside and outside of a system. We leverage this idea to provide an analysis of partitions in neuronal systems. This is applicable to brain architectures at multiple scales, enabling partitions into single neurons, brain regions, and brain-wide networks. This treatment is based upon the canonical micro-circuitry used in empirical studies of effective connectivity, so as to speak directly to practical applications. The notion of effective connectivity depends upon the dynamic coupling between functional units, whose form recapitulates that of a Markov blanket at each level of analysis. The nuance afforded by partitioning neural systems in this way highlights certain limitations of 'modular' perspectives of brain function that only consider a single level of description.

Dynamic Causal Modeling for fMRI With Wilson-Cowan-Based Neuronal Equations

Dynamic causal modeling (DCM) is an analysis technique that has been successfully used to infer about directed connectivity between brain regions based on imaging data such as functional magnetic resonance imaging (fMRI). Most variants of DCM for fMRI rely on a simple bilinear differential equation for neural activation, making it difficult to interpret the results in terms of local neural dynamics. In this work, we introduce a modification to DCM for fMRI by replacing the bilinear equation with a non-linear Wilson-Cowan based equation and use Bayesian Model Comparison (BMC) to show that this modification improves the model evidences. Improved model evidence of the non-linear model is shown for our empirical data (imitation of facial expressions) and validated by synthetic data as well as an empirical test dataset (attention to visual motion) used in previous foundational papers. For our empirical data, we conduct the analysis for a group of 42 healthy participants who performed an imitation task, activating regions putatively containing the human mirror neuron system (MNS). In this regard, we build 540 models as one family for comparing the standard bilinear with the modified Wilson-Cowan models on the family-level. Using this modification, we can interpret the sigmoid transfer function as an averaged f-I curve of many neurons in a single region with a sigmoidal format. In this way, we can make a direct inference from the macroscopic model to detailed microscopic models. The new DCM variant shows superior model evidence on all tested data sets.