Tracking slow modulations in synaptic gain using dynamic causal modelling: validation in epilepsy

The right inferior frontal gyrus as pivotal node and effective regulator of the basal ganglia-thalamocortical response inhibition circuit

Background

The involvement of specific basal ganglia-thalamocortical circuits in response inhibition has been extensively mapped in animal models. However, the pivotal nodes and directed causal regulation within this inhibitory circuit in humans remains controversial.

Objective

The main aim of the present study was to determine the causal information flow and critical nodes in the basal ganglia-thalamocortical inhibitory circuits and also to examine whether these are modulated by biological factors (i.e. sex) and behavioral performance.

Methods

Here, we capitalize on the recent progress in robust and biologically plausible directed causal modeling (DCM-PEB) and a large response inhibition dataset (n = 250) acquired with concomitant functional magnetic resonance imaging to determine key nodes, their causal regulation and modulation via biological variables (sex) and inhibitory performance in the inhibitory circuit encompassing the right inferior frontal gyrus (rIFG), caudate nucleus (rCau), globus pallidum (rGP), and thalamus (rThal).

Results

The entire neural circuit exhibited high intrinsic connectivity and response inhibition critically increased causal projections from the rIFG to both rCau and rThal. Direct comparison further demonstrated that response inhibition induced an increasing rIFG inflow and increased the causal regulation of this region over the rCau and rThal. In addition, sex and performance influenced the functional architecture of the regulatory circuits such that women displayed increased rThal self-inhibition and decreased rThal to GP modulation, while better inhibitory performance was associated with stronger rThal to rIFG communication. Furthermore, control analyses did not reveal a similar key communication in a left lateralized model.

Conclusions

Together, these findings indicate a pivotal role of the rIFG as input and causal regulator of subcortical response inhibition nodes.

Computational Modeling of Oddball Sequence Processing Exposes Common and Differential Auditory Network Changes in First-Episode Schizophrenia-Spectrum Disorders and Schizophrenia

BACKGROUND AND HYPOTHESIS

Differences in sound relevance filtering in schizophrenia are proposed to represent a key index of biological changes in brain function in the illness. This study featured a computational modeling approach to test the hypothesis that processing differences might already be evident in first-episode, becoming more pronounced in the established illness.

STUDY DESIGN

Auditory event-related potentials to a typical oddball sequence (rare pitch deviations amongst regular sounds) were recorded from 90 persons with schizophrenia-spectrum disorders (40 first-episode schizophrenia-spectrum, 50 established illness) and age-matched healthy controls. The data were analyzed using dynamic causal modeling to identify the changes in effective connectivity that best explained group differences.

STUDY RESULTS

Group differences were linked to intrinsic (within brain region) connectivity changes. In activity-dependent measures these were restricted to the left auditory cortex in first-episode schizophrenia-spectrum but were more widespread in the established illness. Modeling suggested that both established illness and first-episode schizophrenia-spectrum groups expressed significantly lower inhibition of inhibitory interneuron activity and altered gain on superficial pyramidal cells with the data indicative of differences in both putative N-methyl-d-aspartate glutamate receptor activity-dependent plasticity and classic neuromodulation.

CONCLUSIONS

The study provides further support for the notion that examining the ability to alter responsiveness to structured sound sequences in schizophrenia and first-episode schizophrenia-spectrum could be informative to uncovering the nature and progression of changes in brain function during the illness. Furthermore, modeling suggested that limited differences present at first-episode schizophrenia-spectrum may become more expansive with illness progression.

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.

Delineating epileptogenic networks using brain imaging data and personalized modeling in drug-resistant epilepsy

Precise estimates of epileptogenic zone networks (EZNs) are crucial for planning intervention strategies to treat drug-resistant focal epilepsy. Here, we present the virtual epileptic patient (VEP), a workflow that uses personalized brain models and machine learning methods to estimate EZNs and to aid surgical strategies. The structural scaffold of the patient-specific whole-brain network model is constructed from anatomical T1 and diffusion-weighted magnetic resonance imaging. Each network node is equipped with a mathematical dynamical model to simulate seizure activity. Bayesian inference methods sample and optimize key parameters of the personalized model using functional stereoelectroencephalography recordings of patients' seizures. These key parameters together with their personalized model determine a given patient's EZN. Personalized models were further used to predict the outcome of surgical intervention using virtual surgeries. We evaluated the VEP workflow retrospectively using 53 patients with drug-resistant focal epilepsy. VEPs reproduced the clinically defined EZNs with a precision of 0.6, where the physical distance between epileptogenic regions identified by VEP and the clinically defined EZNs was small. Compared with the resected brain regions of 25 patients who underwent surgery, VEP showed lower false discovery rates in seizure-free patients (mean, 0.028) than in non-seizure-free patients (mean, 0.407). VEP is now being evaluated in an ongoing clinical trial (EPINOV) with an expected 356 prospective patients with epilepsy.

Task-Dependent Effective Connectivity of the Reward Network During Food Cue-Reactivity: A Dynamic Causal Modeling Investigation

Neural reactivity to food cues may play a central role in overeating and excess weight gain. Functional magnetic resonance imaging (fMRI) studies have implicated regions of the reward network in dysfunctional food cue-reactivity, but neural interactions underlying observed patterns of signal change remain poorly understood. Fifty overweight and obese participants with self-reported cue-induced food craving viewed food and neutral cues during fMRI scanning. Regions of the reward network with significantly greater food versus neutral cue-reactivity were used to specify plausible models of task-related neural interactions underlying the observed blood oxygenation level-dependent (BOLD) signal, and a bi-hemispheric winning model was identified in a dynamic causal modeling (DCM) framework. Neuro-behavioral correlations are investigated with group factor analysis (GFA) and Pearson's correlation tests. The ventral tegmental area (VTA), amygdalae, and orbitofrontal cortices (OFC) showed significant food cue-reactivity. DCM suggests these activations are produced by largely reciprocal dynamic signaling between these regions, with food cues causing regional disinhibition and an apparent shifting of activity to the right amygdala. Intrinsic self-inhibition in the VTA and right amygdala is negatively correlated with measures of food craving and hunger and right-amygdalar disinhibition by food cues is associated with the intensity of cue-induced food craving, but no robust cross-unit latent factors were identified between the neural group and behavioral or demographic variable groups. Our results suggest a rich array of dynamic signals drive reward network cue-reactivity, with the amygdalae mediating much of the dynamic signaling between the VTA and OFCs. Neuro-behavioral correlations suggest particularly crucial roles for the VTA, right amygdala, and the right OFC-amygdala connection but the more robust GFA identified no cross-unit factors, so these correlations should be interpreted with caution. This investigation provides novel insights into dynamic circuit mechanisms with etiologic relevance to obesity, suggesting pathways in biomarker development and intervention.

Investigation of two neural mass models for DCM-based effective connectivity inference in temporal epilepsy

BACKGROUND AND OBJECTIVE

Recently, spectral Dynamic Causal Modelling (DCM) has been used increasingly to infer effective connectivity from epileptic intracranial electroencephalographic (iEEG) signals. In this context, the Physiology-Based Model (PBM), a neural mass model, is used as a generative model. However, previous studies have highlighted out the inability of PBM to properly describe iEEG signals with specific power spectral densities (PSDs). More precisely, PSDs that have multiple peaks around β and γ rhythms (i.e. spectral characteristics at seizure onset) are concerned.

METHODS

To cope with this limitation, an alternative neural mass model, called the complete PBM (cPBM), is investigated. The spectral DCM and two recent variants are used to evaluate the relevance of cPBM over PBM.

RESULTS

The study is conducted on both simulated signals and real epileptic iEEG recordings. Our results confirm that, compared to PBM, cPBM shows (i) more ability to model the desired PSDs and (ii) lower numerical complexity whatever the method.

CONCLUSIONS

Thanks to its intrinsic and extrinsic connectivity parameters as well as the input coming into the fast inhibitory subpopulation, the cPBM provides a more expressive model of PSDs, leading to a better understanding of epileptic patterns and DCM-based effective connectivity inference.

State-Dependent Effective Connectivity in Resting-State fMRI

The human brain at rest exhibits intrinsic dynamics transitioning among the multiple metastable states of the inter-regional functional connectivity. Accordingly, the demand for exploring the state-specific functional connectivity increases for a deeper understanding of mental diseases. Functional connectivity, however, lacks information about the directed causal influences among the brain regions, called effective connectivity. This study presents the dynamic causal modeling (DCM) framework to explore the state-dependent effective connectivity using spectral DCM for the resting-state functional MRI (rsfMRI). We established the sequence of brain states using the hidden Markov model with the multivariate autoregressive coefficients of rsfMRI, summarizing the functional connectivity. We decomposed the state-dependent effective connectivity using a parametric empirical Bayes scheme that models the effective connectivity of consecutive windows with the time course of the discrete states as regressors. We showed the plausibility of the state-dependent effective connectivity analysis in a simulation setting. To test the clinical applicability, we applied the proposed method to characterize the state- and subtype-dependent effective connectivity of the default mode network in children with combined-type attention deficit hyperactivity disorder (ADHD-C) compared with age-matched, typically developed children (TDC). All 88 children were subtyped according to the occupation times (i.e., dwell times) of the three dominant functional connectivity states, independently of clinical diagnosis. The state-dependent effective connectivity differences between ADHD-C and TDC according to the subtypes and those between the subtypes of ADHD-C were expressed mainly in self-inhibition, magnifying the importance of excitation inhibition balance in the subtyping. These findings provide a clear motivation for decomposing the state-dependent dynamic effective connectivity and state-dependent analysis of the directed coupling in exploring mental diseases.

Identifying spatio-temporal seizure propagation patterns in epilepsy using Bayesian inference

Focal drug resistant epilepsy is a neurological disorder characterized by seizures caused by abnormal activity originating in one or more regions together called as epileptogenic zone. Treatment for such patients involves surgical resection of affected regions. Epileptogenic zone is typically identified using stereotactic EEG recordings from the electrodes implanted into the patient's brain. Identifying the epileptogenic zone is a challenging problem due to the spatial sparsity of electrode implantation. We propose a probabilistic hierarchical model of seizure propagation patterns, based on a phenomenological model of seizure dynamics called Epileptor. Using Bayesian inference, the Epileptor model is optimized to build patient specific virtual models that best fit to the log power of intracranial recordings. First, accuracy of the model predictions and identifiability of the model are investigated using synthetic data. Then, model predictions are evaluated against a retrospective patient cohort of 25 patients with varying surgical outcomes. In the patients who are seizure free after surgery, model predictions showed good match with the clinical hypothesis. In patients where surgery failed to achieve seizure freedom model predictions showed a strong mismatch. Our results demonstrate that proposed probabilistic model could be a valuable tool to aid the clinicians in identifying the seizure focus.

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.

Data-driven method to infer the seizure propagation patterns in an epileptic brain from intracranial electroencephalography

Surgical interventions in epileptic patients aimed at the removal of the epileptogenic zone have success rates at only 60-70%. This failure can be partly attributed to the insufficient spatial sampling by the implanted intracranial electrodes during the clinical evaluation, leading to an incomplete picture of spatio-temporal seizure organization in the regions that are not directly observed. Utilizing the partial observations of the seizure spreading through the brain network, complemented by the assumption that the epileptic seizures spread along the structural connections, we infer if and when are the unobserved regions recruited in the seizure. To this end we introduce a data-driven model of seizure recruitment and propagation across a weighted network, which we invert using the Bayesian inference framework. Using a leave-one-out cross-validation scheme on a cohort of 45 patients we demonstrate that the method can improve the predictions of the states of the unobserved regions compared to an empirical estimate that does not use the structural information, yet it is on the same level as the estimate that takes the structure into account. Furthermore, a comparison with the performed surgical resection and the surgery outcome indicates a link between the inferred excitable regions and the actual epileptogenic zone. The results emphasize the importance of the structural connectome in the large-scale spatio-temporal organization of epileptic seizures and introduce a novel way to integrate the patient-specific connectome and intracranial seizure recordings in a whole-brain computational model of seizure spread.

The Bayesian Virtual Epileptic Patient: A probabilistic framework designed to infer the spatial map of epileptogenicity in a personalized large-scale brain model of epilepsy spread

Despite the importance and frequent use of Bayesian frameworks in brain network modeling for parameter inference and model prediction, the advanced sampling algorithms implemented in probabilistic programming languages to overcome the inference difficulties have received relatively little attention in this context. In this technical note, we propose a probabilistic framework, namely the Bayesian Virtual Epileptic Patient (BVEP), which relies on the fusion of structural data of individuals to infer the spatial map of epileptogenicity in a personalized large-scale brain model of epilepsy spread. To invert the individualized whole-brain model employed in this study, we use the recently developed algorithms known as No-U-Turn Sampler (NUTS) as well as Automatic Differentiation Variational Inference (ADVI). Our results indicate that NUTS and ADVI accurately estimate the degree of epileptogenicity of brain regions, therefore, the hypothetical brain areas responsible for the seizure initiation and propagation, while the convergence diagnostics and posterior behavior analysis validate the reliability of the estimations. Moreover, we illustrate the efficiency of the transformed non-centered parameters in comparison to centered form of parameterization. The Bayesian framework used in this work proposes an appropriate patient-specific strategy for estimating the epileptogenicity of the brain regions to improve outcome after epilepsy surgery.

Dissociable neural effects of temporal expectations due to passage of time and contextual probability

The human brain is equipped with complex mechanisms to track the changing probability of events in time. While the passage of time itself usually leads to a mounting expectation, context can provide additional information about when events are likely to happen. In this study we dissociate these two sources of temporal expectation in terms of their neural correlates and underlying brain connectivity patterns. We analysed magnetoencephalographic (MEG) data acquired from N = 24 healthy participants listening to auditory stimuli. These stimuli could be presented at different temporal intervals but occurred most often at intermediate intervals, forming a contextual probability distribution. Evoked MEG response amplitude was sensitive to both passage of time (time elapsed since the cue) and contextual probability, albeit at different latencies: the effects of passage of time were observed earlier than the effects of context. The underlying sources of MEG activity were also different across the two types of temporal prediction: the effects of passage of time were localised to early auditory regions and superior temporal gyri, while context was additionally linked to activity in inferior parietal cortices. Finally, these differences were modelled using biophysical (dynamic causal) modelling: passage of time was explained in terms of widespread gain modulation and decreased prediction error signalling at lower levels of the hierarchy, while contextual expectation led to more localised gain modulation and decreased prediction error signalling at higher levels of the hierarchy. These results present a comprehensive account of how independent sources of temporal prediction may be differentially expressed in cortical circuits.

Structure learning in coupled dynamical systems and dynamic causal modelling

Identifying a coupled dynamical system out of many plausible candidates, each of which could serve as the underlying generator of some observed measurements, is a profoundly ill-posed problem that commonly arises when modelling real-world phenomena. In this review, we detail a set of statistical procedures for inferring the structure of nonlinear coupled dynamical systems (structure learning), which has proved useful in neuroscience research. A key focus here is the comparison of competing models of network architectures-and implicit coupling functions-in terms of their Bayesian model evidence. These methods are collectively referred to as dynamic causal modelling. We focus on a relatively new approach that is proving remarkably useful, namely Bayesian model reduction, which enables rapid evaluation and comparison of models that differ in their network architecture. We illustrate the usefulness of these techniques through modelling neurovascular coupling (cellular pathways linking neuronal and vascular systems), whose function is an active focus of research in neurobiology and the imaging of coupled neuronal systems. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

Dynamic causal modelling of fluctuating connectivity in resting-state EEG

Functional and effective connectivity are known to change systematically over time. These changes might be explained by several factors, including intrinsic fluctuations in activity-dependent neuronal coupling and contextual factors, like experimental condition and time. Furthermore, contextual effects may be subject-specific or conserved over subjects. To characterize fluctuations in effective connectivity, we used dynamic causal modelling (DCM) of cross spectral responses over 1- min of electroencephalogram (EEG) recordings during rest, divided into 1-sec windows. We focused on two intrinsic networks: the default mode and the saliency network. DCM was applied to estimate connectivity in each time-window for both networks. Fluctuations in DCM connectivity parameters were assessed using hierarchical parametric empirical Bayes (PEB). Within-subject, between-window effects were modelled with a second-level linear model with temporal basis functions as regressors. This procedure was conducted for every subject separately. Bayesian model reduction was then used to assess which (combination of) temporal basis functions best explain dynamic connectivity over windows. A third (between-subject) level model was used to infer which dynamic connectivity parameters are conserved over subjects. Our results indicate that connectivity fluctuations in the default mode network and to a lesser extent the saliency network comprised both subject-specific components and a common component. For both networks, connections to higher order regions appear to monotonically increase during the 1- min period. These results not only establish the predictive validity of dynamic connectivity estimates - in virtue of detecting systematic changes over subjects - they also suggest a network-specific dissociation in the relative contribution of fluctuations in connectivity that depend upon experimental context. We envisage these procedures could be useful for characterizing brain state transitions that may be explained by their cognitive or neuropathological underpinnings.

Disentangling causal webs in the brain using functional magnetic resonance imaging: A review of current approaches

In the past two decades, functional Magnetic Resonance Imaging (fMRI) has been used to relate neuronal network activity to cognitive processing and behavior. Recently this approach has been augmented by algorithms that allow us to infer causal links between component populations of neuronal networks. Multiple inference procedures have been proposed to approach this research question but so far, each method has limitations when it comes to establishing whole-brain connectivity patterns. In this paper, we discuss eight ways to infer causality in fMRI research: Bayesian Nets, Dynamical Causal Modelling, Granger Causality, Likelihood Ratios, Linear Non-Gaussian Acyclic Models, Patel's Tau, Structural Equation Modelling, and Transfer Entropy. We finish with formulating some recommendations for the future directions in this area.

Seizure pathways: A model-based investigation

We present the results of a model inversion algorithm for electrocorticography (ECoG) data recorded during epileptic seizures. The states and parameters of neural mass models were tracked during a total of over 3000 seizures from twelve patients with focal epilepsy. These models provide an estimate of the effective connectivity within intracortical circuits over the time course of seizures. Observing the dynamics of effective connectivity provides insight into mechanisms of seizures. Estimation of patients seizure dynamics revealed: 1) a highly stereotyped pattern of evolution for each patient, 2) distinct sub-groups of onset mechanisms amongst patients, and 3) different offset mechanisms for long and short seizures. Stereotypical dynamics suggest that, once initiated, seizures follow a deterministic path through the parameter space of a neural model. Furthermore, distinct sub-populations of patients were identified based on characteristic motifs in the dynamics at seizure onset. There were also distinct patterns between long and short duration seizures that were related to seizure offset. Understanding how these different patterns of seizure evolution arise may provide new insights into brain function and guide treatment for epilepsy, since specific therapies may have preferential effects on the various parameters that could potentially be individualized. Methods that unite computational models with data provide a powerful means to generate testable hypotheses for further experimental research. This work provides a demonstration that the hidden connectivity parameters of a neural mass model can be dynamically inferred from data. Our results underscore the power of theoretical models to inform epilepsy management. It is our hope that this work guides further efforts to apply computational models to clinical data.

Not All Predictions Are Equal: "What" and "When" Predictions Modulate Activity in Auditory Cortex through Different Mechanisms

Using predictions based on environmental regularities is fundamental for adaptive behavior. While it is widely accepted that predictions across different stimulus attributes (e.g., time and content) facilitate sensory processing, it is unknown whether predictions across these attributes rely on the same neural mechanism. Here, to elucidate the neural mechanisms of predictions, we combine invasive electrophysiological recordings (human electrocorticography in 4 females and 2 males) with computational modeling while manipulating predictions about content ("what") and time ("when"). We found that "when" predictions increased evoked activity over motor and prefrontal regions both at early (∼180 ms) and late (430-450 ms) latencies. "What" predictability, however, increased evoked activity only over prefrontal areas late in time (420-460 ms). Beyond these dissociable influences, we found that "what" and "when" predictability interactively modulated the amplitude of early (165 ms) evoked responses in the superior temporal gyrus. We modeled the observed neural responses using biophysically realistic neural mass models, to better understand whether "what" and "when" predictions tap into similar or different neurophysiological mechanisms. Our modeling results suggest that "what" and "when" predictability rely on complementary neural processes: "what" predictions increased short-term plasticity in auditory areas, whereas "when" predictability increased synaptic gain in motor areas. Thus, content and temporal predictions engage complementary neural mechanisms in different regions, suggesting domain-specific prediction signaling along the cortical hierarchy. Encoding predictions through different mechanisms may endow the brain with the flexibility to efficiently signal different sources of predictions, weight them by their reliability, and allow for their encoding without mutual interference.SIGNIFICANCE STATEMENT Predictions of different stimulus features facilitate sensory processing. However, it is unclear whether predictions of different attributes rely on similar or different neural mechanisms. By combining invasive electrophysiological recordings of cortical activity with experimental manipulations of participants' predictions about content and time of acoustic events, we found that the two types of predictions had dissociable influences on cortical activity, both in terms of the regions involved and the timing of the observed effects. Further, our biophysical modeling analysis suggests that predictability of content and time rely on complementary neural processes: short-term plasticity in auditory areas and synaptic gain in motor areas, respectively. This suggests that predictions of different features are encoded with complementary neural mechanisms in different brain regions.

Calcium imaging and dynamic causal modelling reveal brain-wide changes in effective connectivity and synaptic dynamics during epileptic seizures

Pathophysiological explanations of epilepsy typically focus on either the micro/mesoscale (e.g. excitation-inhibition imbalance), or on the macroscale (e.g. network architecture). Linking abnormalities across spatial scales remains difficult, partly because of technical limitations in measuring neuronal signatures concurrently at the scales involved. Here we use light sheet imaging of the larval zebrafish brain during acute epileptic seizure induced with pentylenetetrazole. Spectral changes of spontaneous neuronal activity during the seizure are then modelled using neural mass models, allowing Bayesian inference on changes in effective network connectivity and their underlying synaptic dynamics. This dynamic causal modelling of seizures in the zebrafish brain reveals concurrent changes in synaptic coupling at macro- and mesoscale. Fluctuations of both synaptic connection strength and their temporal dynamics are required to explain observed seizure patterns. These findings highlight distinct changes in local (intrinsic) and long-range (extrinsic) synaptic transmission dynamics as a possible seizure pathomechanism and illustrate how our Bayesian model inversion approach can be used to link existing neural mass models of seizure activity and novel experimental methods.

Network dynamics in the healthy and epileptic developing brain

Electroencephalography (EEG) allows recording of cortical activity at high temporal resolution. EEG recordings can be summarized along different dimensions using network-level quantitative measures, such as channel-to-channel correlation, or band power distributions across channels. These reveal network patterns that unfold over a range of different timescales and can be tracked dynamically. Here we describe the dynamics of network state transitions in EEG recordings of spontaneous brain activity in normally developing infants and infants with severe early infantile epileptic encephalopathies (n = 8, age: 1–8 months). We describe differences in measures of EEG dynamics derived from band power, and correlation-based summaries of network-wide brain activity. We further show that EEGs from different patient groups and controls may be distinguishable on a small set of the novel quantitative measures introduced here, which describe dynamic network state switching. Quantitative measures related to the sharpness of switching from one correlation pattern to another show the largest differences between groups. These findings reveal that the early epileptic encephalopathies are associated with characteristic dynamic features at the network level. Quantitative network-based analyses like the one presented here may in the future inform the clinical use of quantitative EEG for diagnosis.

The Cumulative Effects of Predictability on Synaptic Gain in the Auditory Processing Stream

Stimulus predictability can lead to substantial modulations of brain activity, such as shifts in sustained magnetic field amplitude, measured with magnetoencephalography (MEG). Here, we provide a mechanistic explanation of these effects using MEG data acquired from healthy human volunteers (N = 13, 7 female). In a source-level analysis of induced responses, we established the effects of orthogonal predictability manipulations of rapid tone-pip sequences (namely, sequence regularity and alphabet size) along the auditory processing stream. In auditory cortex, regular sequences with smaller alphabets induced greater gamma activity. Furthermore, sequence regularity shifted induced activity in frontal regions toward higher frequencies. To model these effects in terms of the underlying neurophysiology, we used dynamic causal modeling for cross-spectral density and estimated slow fluctuations in neural (postsynaptic) gain. Using the model-based parameters, we accurately explain the sensor-level sustained field amplitude, demonstrating that slow changes in synaptic efficacy, combined with sustained sensory input, can result in profound and sustained effects on neural responses to predictable sensory streams.

SIGNIFICANCE STATEMENT Brain activity can be strongly modulated by the predictability of stimuli it is currently processing. An example of such a modulation is a shift in sustained magnetic field amplitude, measured with magnetoencephalography. Here, we provide a mechanistic explanation of these effects. First, we establish the oscillatory neural correlates of independent predictability manipulations in hierarchically distinct areas of the auditory processing stream. Next, we use a biophysically realistic computational model to explain these effects in terms of the underlying neurophysiology. Finally, using the model-based parameters describing neural gain modulation, we can explain the previously unexplained effects observed at the sensor level. This demonstrates that slow modulations of synaptic gain can result in profound and sustained effects on neural activity.

A local adjustment strategy for the initialization of dynamic causal modelling to infer effective connectivity in brain epileptic structures

This paper addresses the question of effective connectivity in the human cerebral cortex in the context of epilepsy. Among model based approaches to infer brain connectivity, spectral Dynamic Causal Modelling is a conventional technique for which we propose an alternative to estimate cross spectral density. The proposed strategy we investigated tackles the sub-estimation of the free energy using the well-known variational Expectation-Maximization algorithm highly sensitive to the initialization of the parameters vector by a permanent local adjustment of the initialization process. The performance of the proposed strategy in terms of effective connectivity identification is assessed using simulated data generated by a neuronal mass model (simulating unidirectional and bidirectional flows) and real epileptic intracerebral Electroencephalographic signals. Results show the efficiency of proposed approach compared to the conventional Dynamic Causal Modelling and the one wherein a deterministic annealing scheme is employed.

Quantifying connectivity in a physiology based model using adaptive dynamic causal modelling

This paper proposes an Adaptive Dynamic Causal Modelling based approach to detect and quantify effective connectivity in human brain structures injured by epileptic activities. The identification of the parameters in the physiology based model subtended the Electroencephalographic observations is performed by improving the optimization step in the Expectation Maximization algorithm. Considering unidirectional flow propagation, we show the efficiency of our proposed approach compared to the conventional technique.

Dynamic causal modelling of seizure activity in a rat model

This paper presents a physiological account of seizure activity and its evolution over time using a rat model of induced epilepsy. We analyse spectral activity recorded in the hippocampi of three rats who received kainic acid injections in the right hippocampus. We use dynamic causal modelling of seizure activity and Bayesian model reduction to identify the key synaptic and connectivity parameters that underlie seizure onset. Using recent advances in hierarchical modelling (parametric empirical Bayes), we characterise seizure onset in terms of slow fluctuations in synaptic excitability of specific neuronal populations. Our results suggest differences in the pathophysiology - of seizure activity in the lesioned versus the non-lesioned hippocampus - with pronounced changes in excitation-inhibition balance and temporal summation on the lesioned side. In particular, our analyses suggest that marked reductions in the synaptic time constant of the deep pyramidal cells and the self-inhibition of inhibitory interneurons (in the lesioned hippocampus) are sufficient to explain changes in spectral activity. Although these synaptic changes are consistent over rats, the resulting electrophysiological phenotype can be quite diverse.

Suppressing epileptic activity in a neural mass model using a closed-loop proportional-integral controller

Closed-loop control is a promising deep brain stimulation (DBS) strategy that could be used to suppress high-amplitude epileptic activity. However, there are currently no analytical approaches to determine the stimulation parameters for effective and safe treatment protocols. Proportional-integral (PI) control is the most extensively used closed-loop control scheme in the field of control engineering because of its simple implementation and perfect performance. In this study, we took Jansen's neural mass model (NMM) as a test bed to develop a PI-type closed-loop controller for suppressing epileptic activity. A graphical stability analysis method was employed to determine the stabilizing region of the PI controller in the control parameter space, which provided a theoretical guideline for the choice of the PI control parameters. Furthermore, we established the relationship between the parameters of the PI controller and the parameters of the NMM in the form of a stabilizing region, which provided insights into the mechanisms that may suppress epileptic activity in the NMM. The simulation results demonstrated the validity and effectiveness of the proposed closed-loop PI control scheme.

New Perspectives on Spontaneous Brain Activity: Dynamic Networks and Energy Matter

Spontaneous brain activity has received increasing attention as demonstrated by the exponential rise in the number of published article on this topic over the last 30 years. Such “intrinsic” brain activity, generated in the absence of an explicit task, is frequently associated with resting-state or default-mode networks (DMN)s. The focus on characterizing spontaneous brain activity promises to shed new light on questions concerning the structural and functional architecture of the brain and how they are related to “mind”. However, many critical questions have yet to be addressed. In this review, we focus on a scarcely explored area, specifically the energetic requirements and constraints of spontaneous activity, taking into account both thermodynamical and informational perspectives. We argue that the “classical” definitions of spontaneous activity do not take into account an important feature, that is, the critical thermodynamic energetic differences between spontaneous and evoked brain activity. Spontaneous brain activity is associated with slower oscillations compared with evoked, task-related activity, hence it exhibits lower levels of enthalpy and “free-energy” (i.e., the energy that can be converted to do work), thus supporting noteworthy thermodynamic energetic differences between spontaneous and evoked brain activity. Increased spike frequency during evoked activity has a significant metabolic cost, consequently, brain functions traditionally associated with spontaneous activity, such as mind wandering, require less energy that other nervous activities. We also review recent empirical observations in neuroscience, in order to capture how spontaneous brain dynamics and mental function can be embedded in a non-linear dynamical framework, which considers nervous activity in terms of phase spaces, particle trajectories, random walks, attractors and/or paths at the edge of the chaos. This takes us from the thermodynamic free-energy, to the realm of “variational free-energy”, a theoretical construct pertaining to probability and information theory which allows explanation of unexplored features of spontaneous brain activity.

Eyes Open on Sleep and Wake: In Vivo to In Silico Neural Networks

Functional and effective connectivity of cortical areas are essential for normal brain function under different behavioral states. Appropriate cortical activity during sleep and wakefulness is ensured by the balanced activity of excitatory and inhibitory circuits. Ultimately, fast, millisecond cortical rhythmic oscillations shape cortical function in time and space. On a much longer time scale, brain function also depends on prior sleep-wake history and circadian processes. However, much remains to be established on how the brain operates at the neuronal level in humans during sleep and wakefulness. A key limitation of human neuroscience is the difficulty in isolating neuronal excitation/inhibition drive in vivo. Therefore, computational models are noninvasive approaches of choice to indirectly access hidden neuronal states. In this review, we present a physiologically driven in silico approach, Dynamic Causal Modelling (DCM), as a means to comprehend brain function under different experimental paradigms. Importantly, DCM has allowed for the understanding of how brain dynamics underscore brain plasticity, cognition, and different states of consciousness. In a broader perspective, noninvasive computational approaches, such as DCM, may help to puzzle out the spatial and temporal dynamics of human brain function at different behavioural states.

Estimating Directed Connectivity from Cortical Recordings and Reconstructed Sources

In cognitive neuroscience, electrical brain activity is most commonly recorded at the scalp. In order to infer the contributions and connectivity of underlying neuronal sources within the brain, it is necessary to reconstruct sensor data at the source level. Several approaches to this reconstruction have been developed, thereby solving the so-called implicit inverse problem Michel et al. (Clin Neurophysiol 115:2195–2222, 2004). However, a unifying premise against which to validate these source reconstructions is seldom available. The dataset provided in this work, in which brain activity is simultaneously recorded on the scalp (non-invasively) by electroencephalography (EEG) and on the cortex (invasively) by electrocorticography (ECoG), can be of a great help in this direction. These multimodal recordings were obtained from a macaque monkey under wakefulness and sedation. Our primary goal was to establish the connectivity architecture between two sources of interest (frontal and parietal), and to assess how their coupling changes over the conditions. We chose these sources because previous studies have shown that the connections between them are modified by anaesthesia Boly et al. (J Neurosci 32:7082–7090, 2012). Our secondary goal was to evaluate the consistency of the connectivity results when analyzing sources recorded from invasive data (128 implanted ECoG sources) and source activity reconstructed from scalp recordings (19 EEG sensors) at the same locations as the ECoG sources. We conclude that the directed connectivity in the frequency domain between cortical sources reconstructed from scalp EEG is qualitatively similar to the connectivity inferred directly from cortical recordings, using both data-driven (directed transfer function) and biologically grounded (dynamic causal modelling) methods. Furthermore, the connectivity changes identified were consistent with previous findings Boly et al. (J Neurosci 32:7082–7090, 2012). Our findings suggest that inferences about directed connectivity based upon non-invasive electrophysiological data have construct validity in relation to invasive recordings.

Dynamic causal modelling of electrographic seizure activity using Bayesian belief updating

Seizure activity in EEG recordings can persist for hours with seizure dynamics changing rapidly over time and space. To characterise the spatiotemporal evolution of seizure activity, large data sets often need to be analysed. Dynamic causal modelling (DCM) can be used to estimate the synaptic drivers of cortical dynamics during a seizure; however, the requisite (Bayesian) inversion procedure is computationally expensive. In this note, we describe a straightforward procedure, within the DCM framework, that provides efficient inversion of seizure activity measured with non-invasive and invasive physiological recordings; namely, EEG/ECoG. We describe the theoretical background behind a Bayesian belief updating scheme for DCM. The scheme is tested on simulated and empirical seizure activity (recorded both invasively and non-invasively) and compared with standard Bayesian inversion. We show that the Bayesian belief updating scheme provides similar estimates of time-varying synaptic parameters, compared to standard schemes, indicating no significant qualitative change in accuracy. The difference in variance explained was small (less than 5%). The updating method was substantially more efficient, taking approximately 5–10 min compared to approximately 1–2 h. Moreover, the setup of the model under the updating scheme allows for a clear specification of how neuronal variables fluctuate over separable timescales. This method now allows us to investigate the effect of fast (neuronal) activity on slow fluctuations in (synaptic) parameters, paving a way forward to understand how seizure activity is generated.

•

We describe a DCM procedure that provides efficient inversion of seizure activity.

•

Similar accuracy but substantially more efficient compared to standard DCM methods.

•

Physiological fluctuations over different timescales can be specified.

•

This scheme should contribute to understanding seizure activity using DCM.

Characterising seizures in anti-NMDA-receptor encephalitis with dynamic causal modelling

We characterised the pathophysiology of seizure onset in terms of slow fluctuations in synaptic efficacy using EEG in patients with anti-N-methyl-d-aspartate receptor (NMDA-R) encephalitis. EEG recordings were obtained from two female patients with anti-NMDA-R encephalitis with recurrent partial seizures (ages 19 and 31). Focal electrographic seizure activity was localised using an empirical Bayes beamformer. The spectral density of reconstructed source activity was then characterised with dynamic causal modelling (DCM). Eight models were compared for each patient, to evaluate the relative contribution of changes in intrinsic (excitatory and inhibitory) connectivity and endogenous afferent input. Bayesian model comparison established a role for changes in both excitatory and inhibitory connectivity during seizure activity (in addition to changes in the exogenous input). Seizures in both patients were associated with a sequence of changes in inhibitory and excitatory connectivity; a transient increase in inhibitory connectivity followed by a transient increase in excitatory connectivity and a final peak of excitatory–inhibitory balance at seizure offset. These systematic fluctuations in excitatory and inhibitory gain may be characteristic of (anti NMDA-R encephalitis) seizures. We present these results as a case study and replication to motivate analyses of larger patient cohorts, to see whether our findings generalise and further characterise the mechanisms of seizure activity in anti-NMDA-R encephalitis.

•

We characterised seizures in patient with anti-NMDA-R encephalitis using EEG.

•

Dynamic causal modelling was used to estimate causes of seizure activity.

•

Characteristic variation of excitatory–inhibitory balance during seizure activity.

•

This variation was seen for seizures within and between patients.