Consistent spectral predictors for dynamic causal models of steady-state responses

Reliability of dynamic causal modelling of resting-state magnetoencephalography

This study assesses the reliability of resting-state dynamic causal modelling (DCM) of magnetoencephalography (MEG) under conductance-based canonical microcircuit models, in terms of both posterior parameter estimates and model evidence. We use resting-state MEG data from two sessions, acquired 2 weeks apart, from a cohort with high between-subject variance arising from Alzheimer's disease. Our focus is not on the effect of disease, but on the reliability of the methods (as within-subject between-session agreement), which is crucial for future studies of disease progression and drug intervention. To assess the reliability of first-level DCMs, we compare model evidence associated with the covariance among subject-specific free energies (i.e., the 'quality' of the models) with versus without interclass correlations. We then used parametric empirical Bayes (PEB) to investigate the differences between the inferred DCM parameter probability distributions at the between subject level. Specifically, we examined the evidence for or against parameter differences (i) within-subject, within-session, and between-epochs; (ii) within-subject between-session; and (iii) within-site between-subjects, accommodating the conditional dependency among parameter estimates. We show that for data acquired close in time, and under similar circumstances, more than 95% of inferred DCM parameters are unlikely to differ, speaking to mutual predictability over sessions. Using PEB, we show a reciprocal relationship between a conventional definition of 'reliability' and the conditional dependency among inferred model parameters. Our analyses confirm the reliability and reproducibility of the conductance-based DCMs for resting-state neurophysiological data. In this respect, the implicit generative modelling is suitable for interventional and longitudinal studies of neurological and psychiatric disorders.

Linking fast and slow: The case for generative models

A pervasive challenge in neuroscience is testing whether neuronal connectivity changes over time due to specific causes, such as stimuli, events, or clinical interventions. Recent hardware innovations and falling data storage costs enable longer, more naturalistic neuronal recordings. The implicit opportunity for understanding the self-organised brain calls for new analysis methods that link temporal scales: from the order of milliseconds over which neuronal dynamics evolve, to the order of minutes, days, or even years over which experimental observations unfold. This review article demonstrates how hierarchical generative models and Bayesian inference help to characterise neuronal activity across different time scales. Crucially, these methods go beyond describing statistical associations among observations and enable inference about underlying mechanisms. We offer an overview of fundamental concepts in state-space modeling and suggest a taxonomy for these methods. Additionally, we introduce key mathematical principles that underscore a separation of temporal scales, such as the slaving principle, and review Bayesian methods that are being used to test hypotheses about the brain with multiscale data. We hope that this review will serve as a useful primer for experimental and computational neuroscientists on the state of the art and current directions of travel in the complex systems modelling literature.

Effective connectivity of working memory performance: a DCM study of MEG data

Visual working memory (WM) engages several nodes of a large-scale network that includes frontal, parietal, and visual regions; however, little is understood about how these regions interact to support WM behavior. In particular, it is unclear whether network dynamics during WM maintenance primarily represent feedforward or feedback connections. This question has important implications for current debates about the relative roles of frontoparietal and visual regions in WM maintenance. In the current study, we investigated the network activity supporting WM using MEG data acquired while healthy subjects performed a multi-item delayed estimation WM task. We used computational modeling of behavior to discriminate correct responses (high accuracy trials) from two different types of incorrect responses (low accuracy and swap trials), and dynamic causal modeling of MEG data to measure effective connectivity. We observed behaviorally dependent changes in effective connectivity in a brain network comprising frontoparietal and early visual areas. In comparison with high accuracy trials, frontoparietal and frontooccipital networks showed disrupted signals depending on type of behavioral error. Low accuracy trials showed disrupted feedback signals during early portions of WM maintenance and disrupted feedforward signals during later portions of maintenance delay, while swap errors showed disrupted feedback signals during the whole delay period. These results support a distributed model of WM that emphasizes the role of visual regions in WM storage and where changes in large scale network configurations can have important consequences for memory-guided behavior.

Clinical and electrophysiological correlates of hopelessness in the context of suicide risk✰

Hopelessness is a key risk factor for suicide. This analysis explored whether hopelessness indicates a recent suicide crisis state and is linked with magnetoencephalography (MEG) oscillatory power and effective connectivity differences. Change in hopelessness ratings and effective connectivity post-ketamine were also evaluated in a subsample of high-risk individuals to evaluate correlates of dynamic changes over time. Participants (66F;44 M;1 transgender) included individuals with suicide crisis in the last two weeks (High Risk (HR), n = 14), those with past suicide attempt but no recent suicide ideation (SI) (Low Risk (LR), n = 37), clinical controls (CC, n = 33), and healthy volunteers at minimal risk (MinR, n = 27). MEG oscillatory power and clinical hopelessness ratings (via the Beck Hopelessness Scale (BHS)) were evaluated across groups. Dynamic casual modeling (DCM) evaluated connectivity within and between the anterior insula (AI) and anterior cingulate cortex (ACC). A subsample of HR individuals who received ketamine (n = 10) were evaluated at Day 1 post-infusion. The HR group reported the highest levels of hopelessness, even when adjusting for SI. MEG results linked hopelessness with reduced activity across frequency bands in salience network regions, with no group or group-by-interaction effects. Using DCM, the HR group had reduced intrinsic drive from granular Layer IV stellate cells to superficial pyramidal cells in the ACC and AI. In the pilot HR study, reduced hopelessness was linked with increased drive for this same connection post-ketamine. Hopelessness is a possible proxy for suicide risk. Electrophysiological targets for hopelessness include widespread reductions in salience network activity, particularly in the ACC and AI.

Neurochemistry-enriched dynamic causal models of magnetoencephalography, using magnetic resonance spectroscopy

We present a hierarchical empirical Bayesian framework for testing hypotheses about neurotransmitters' concertation as empirical prior for synaptic physiology using ultra-high field magnetic resonance spectroscopy (7T-MRS) and magnetoencephalography data (MEG). A first level dynamic causal modelling of cortical microcircuits is used to infer the connectivity parameters of a generative model of individuals' neurophysiological observations. At the second level, individuals' 7T-MRS estimates of regional neurotransmitter concentration supply empirical priors on synaptic connectivity. We compare the group-wise evidence for alternative empirical priors, defined by monotonic functions of spectroscopic estimates, on subsets of synaptic connections. For efficiency and reproducibility, we used Bayesian model reduction (BMR), parametric empirical Bayes and variational Bayesian inversion. In particular, we used Bayesian model reduction to compare alternative model evidence of how spectroscopic neurotransmitter measures inform estimates of synaptic connectivity. This identifies the subset of synaptic connections that are influenced by individual differences in neurotransmitter levels, as measured by 7T-MRS. We demonstrate the method using resting-state MEG (i.e., task-free recording) and 7T-MRS data from healthy adults. Our results confirm the hypotheses that GABA concentration influences local recurrent inhibitory intrinsic connectivity in deep and superficial cortical layers, while glutamate influences the excitatory connections between superficial and deep layers and connections from superficial to inhibitory interneurons. Using within-subject split-sampling of the MEG dataset (i.e., validation by means of a held-out dataset), we show that model comparison for hypothesis testing can be highly reliable. The method is suitable for applications with magnetoencephalography or electroencephalography, and is well-suited to reveal the mechanisms of neurological and psychiatric disorders, including responses to psychopharmacological interventions.

Neuropharmacological computational analysis of longitudinal electroencephalograms in clozapine-treated patients with schizophrenia using hierarchical dynamic causal modeling

The hierarchical characteristics of the brain are prominent in the pharmacological treatment of psychiatric diseases, primarily targeting cellular receptors that extend upward to intrinsic connectivity within a region, interregional connectivity, and, consequently, clinical observations such as an electroencephalogram (EEG). To understand the long-term effects of neuropharmacological intervention on neurobiological properties at different hierarchical levels, we explored long-term changes in neurobiological parameters of an N-methyl-D-aspartate canonical microcircuit model (CMM-NMDA) in the default mode network (DMN) and auditory hallucination network (AHN) using dynamic causal modeling of longitudinal EEG in clozapine-treated patients with schizophrenia. The neurobiological properties of the CMM-NMDA model associated with symptom improvement in schizophrenia were found across hierarchical levels, from a reduced membrane capacity of the deep pyramidal cell and intrinsic connectivity with the inhibitory population in DMN and intrinsic and extrinsic connectivity in AHN. The medication duration mainly affects the intrinsic connectivity and NMDA time constant in DMN. Virtual perturbation analysis specified the contribution of each parameter to the cross-spectral density (CSD) of the EEG, particularly intrinsic connectivity and membrane capacitances for CSD frequency shifts and progression. It further reveals that excitatory and inhibitory connectivity complements frequency-specific CSD changes, notably the alpha frequency band in DMN. Positive and negative synergistic interactions exist between neurobiological properties primarily within the same region in patients treated with clozapine. The current study shows how computational neuropharmacology helps explore the multiscale link between neurobiological properties and clinical observations and understand the long-term mechanism of neuropharmacological intervention reflected in clinical EEG.

Conductance-based dynamic causal modeling: A mathematical review of its application to cross-power spectral densities

Dynamic Causal Modeling (DCM) is a Bayesian framework for inferring on hidden (latent) neuronal states, based on measurements of brain activity. Since its introduction in 2003 for functional magnetic resonance imaging data, DCM has been extended to electrophysiological data, and several variants have been developed. Their biophysically motivated formulations make these models promising candidates for providing a mechanistic understanding of human brain dynamics, both in health and disease. However, due to their complexity and reliance on concepts from several fields, fully understanding the mathematical and conceptual basis behind certain variants of DCM can be challenging. At the same time, a solid theoretical knowledge of the models is crucial to avoid pitfalls in the application of these models and interpretation of their results. In this paper, we focus on one of the most advanced formulations of DCM, i.e. conductance-based DCM for cross-spectral densities, whose components are described across multiple technical papers. The aim of the present article is to provide an accessible exposition of the mathematical background, together with an illustration of the model's behavior. To this end, we include step-by-step derivations of the model equations, point to important aspects in the software implementation of those models, and use simulations to provide an intuitive understanding of the type of responses that can be generated and the role that specific parameters play in the model. Furthermore, all code utilized for our simulations is made publicly available alongside the manuscript to allow readers an easy hands-on experience with conductance-based DCM.

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.

In Vivo Assay of Cortical Microcircuitry in Frontotemporal Dementia: A Platform for Experimental Medicine Studies

The analysis of neural circuits can provide crucial insights into the mechanisms of neurodegeneration and dementias, and offer potential quantitative biological tools to assess novel therapeutics. Here we use behavioral variant frontotemporal dementia (bvFTD) as a model disease. We demonstrate that inversion of canonical microcircuit models to noninvasive human magnetoencephalography, using dynamic causal modeling, can identify the regional- and laminar-specificity of bvFTD pathophysiology, and their parameters can accurately differentiate patients from matched healthy controls. Using such models, we show that changes in local coupling in frontotemporal dementia underlie the failure to adequately establish sensory predictions, leading to altered prediction error responses in a cortical information-processing hierarchy. Using machine learning, this model-based approach provided greater case-control classification accuracy than conventional evoked cortical responses. We suggest that this approach provides an in vivo platform for testing mechanistic hypotheses about disease progression and pharmacotherapeutics.

A Connectome-Based, Corticothalamic Model of State- and Stimulation-Dependent Modulation of Rhythmic Neural Activity and Connectivity

Rhythmic activity in the brain fluctuates with behaviour and cognitive state, through a combination of coexisting and interacting frequencies. At large spatial scales such as those studied in human M/EEG, measured oscillatory dynamics are believed to arise primarily from a combination of cortical (intracolumnar) and corticothalamic rhythmogenic mechanisms. Whilst considerable progress has been made in characterizing these two types of neural circuit separately, relatively little work has been done that attempts to unify them into a single consistent picture. This is the aim of the present paper. We present and examine a whole-brain, connectome-based neural mass model with detailed long-range cortico-cortical connectivity and strong, recurrent corticothalamic circuitry. This system reproduces a variety of known features of human M/EEG recordings, including spectral peaks at canonical frequencies, and functional connectivity structure that is shaped by the underlying anatomical connectivity. Importantly, our model is able to capture state- (e.g., idling/active) dependent fluctuations in oscillatory activity and the coexistence of multiple oscillatory phenomena, as well as frequency-specific modulation of functional connectivity. We find that increasing the level of sensory drive to the thalamus triggers a suppression of the dominant low frequency rhythms generated by corticothalamic loops, and subsequent disinhibition of higher frequency endogenous rhythmic behaviour of intracolumnar microcircuits. These combine to yield simultaneous decreases in lower frequency and increases in higher frequency components of the M/EEG power spectrum during states of high sensory or cognitive drive. Building on this, we also explored the effect of pulsatile brain stimulation on ongoing oscillatory activity, and evaluated the impact of coexistent frequencies and state-dependent fluctuations on the response of cortical networks. Our results provide new insight into the role played by cortical and corticothalamic circuits in shaping intrinsic brain rhythms, and suggest new directions for brain stimulation therapies aimed at state-and frequency-specific control of oscillatory brain activity.

Generative modelling of the thalamo-cortical circuit mechanisms underlying the neurophysiological effects of ketamine

Cortical recordings of task-induced oscillations following subanaesthetic ketamine administration demonstrate alterations in amplitude, including increases at high-frequencies (gamma) and reductions at low frequencies (theta, alpha). To investigate the population-level interactions underlying these changes, we implemented a thalamo-cortical model (TCM) capable of recapitulating broadband spectral responses. Compared with an existing cortex-only 4-population model, Bayesian Model Selection preferred the TCM. The model was able to accurately and significantly recapitulate ketamine-induced reductions in alpha amplitude and increases in gamma amplitude. Parameter analysis revealed no change in receptor time-constants but significant increases in select synaptic connectivity with ketamine. Significantly increased connections included both AMPA and NMDA mediated connections from layer 2/3 superficial pyramidal cells to inhibitory interneurons and both GABAA and NMDA mediated within-population gain control of layer 5 pyramidal cells. These results support the use of extended generative models for explaining oscillatory data and provide in silico support for ketamine's ability to alter local coupling mediated by NMDA, AMPA and GABA-A.

VEP atlas: An anatomic and functional human brain atlas dedicated to epilepsy patients

BACKGROUND

Several automated parcellation atlases of the human brain have been developed over the past decades, based on various criteria, and have been applied in basic and clinical research.

NEW METHOD

Here we present the Virtual Epileptic Patient (VEP) atlas that offers a new automated brain region parcellation and labeling, which has been developed for the specific use in the domains of epileptology and functional neurosurgery and is able to apply at individual patient's level.

RESULTS

It comprises 162 brain regions, including 73 cortical and 8 subcortical regions per hemisphere. We demonstrate the successful application of the VEP atlas in a cohort of 50 retrospective patients. The structural organization is complemented by the functional variation of stereotactic intracerebral EEG (SEEG) signal data features establishing brain region-specific 3d-maps.

COMPARISON WITH EXISTING METHODS

The VEP atlas integrates both anatomical and functional definitions in the same atlas, adapted to applications for epilepsy patients and individualizable.

CONCLUSION

The covariation of structural and functional organization is the basis for current efforts of patient-specific large-scale brain network modeling exploiting virtual brain technologies for the identification of the epileptogenic regions in an ongoing prospective clinical trial EPINOV.

Lateralized memory circuit dropout in Alzheimer’s disease patients

Altered connectivity within neuronal networks is often observed in Alzheimer’s disease. However, delineating pro-cognitive compensatory changes from pathological network decline relies on characterizing network and task effects together. In this study, we interrogated the dynamics of occipito-temporo-frontal brain networks responsible for implicit and explicit memory processes using high-density EEG and dynamic causal modelling. We examined source-localized network activity from patients with Alzheimer’s disease (n = 21) and healthy controls (n = 21), while they performed both visual recognition (explicit memory) and implicit priming tasks. Parametric empirical Bayes analyses identified significant reductions in temporo-frontal connectivity and in subcortical visual input in patients, specifically in the left hemisphere during the recognition task. There was also slowing in frontal left hemisphere signal transmission during the implicit priming task, with significantly more distinct dropout in connectivity during the recognition task, suggesting that these network drop-out effects are affected by task difficulty. Furthermore, during the implicit memory task, increased right frontal activity was correlated with improved task performance in patients only, suggesting that right-hemisphere compensatory mechanisms may be employed to mitigate left-lateralized network dropout in Alzheimer’s disease. Taken together, these findings suggest that Alzheimer’s disease is associated with lateralized memory circuit dropout and potential compensation from the right hemisphere, at least for simpler memory tasks.

Changes in the Effective Connectivity of the Social Brain When Making Inferences About Close Others vs. the Self

Previous research showed that the ability to make inferences about our own and other’s mental states rely on common brain pathways; particularly in the case of close relationships (e.g., romantic relationships). Despite the evidence for shared neural representations of self and others, less is known about the distributed processing within these common neural networks, particularly whether there are specific patterns of internode communication when focusing on other vs. self. This study aimed to characterize context-sensitive coupling among social brain regions involved in self and other understanding. Participants underwent an fMRI while watching emotional video vignettes of their romantic partner and elaborated on their partner’s (other-condition) or on their own experience (self-condition). We used dynamic causal modeling (DCM) to quantify the associated changes in effective connectivity (EC) in a network of brain regions involved in social cognition including the temporoparietal junction (TPJ), the posterior cingulate (PCC)/precuneus and middle temporal gyrus (MTG). DCM revealed that: the PCC plays a central coordination role within this network, the bilateral MTG receives driving inputs from other nodes suggesting that social information is first processed in language comprehension regions; the right TPJ evidenced a selective increase in its sensitivity when focusing on the other’s experience, relative to focusing on oneself.

Magnetoencephalographic Correlates of Suicidal Ideation in Major Depression

BACKGROUND

Defining the neurobiological underpinnings of suicidal ideation (SI) is crucial to improving our understanding of suicide. This study used magnetoencephalographic gamma power as a surrogate marker for population-level excitation-inhibition balance to explore the underlying neurobiology of SI and depression. In addition, effects of pharmacological intervention with ketamine, which has been shown to rapidly reduce SI and depression, were assessed.

METHODS

Data were obtained from 29 drug-free patients with major depressive disorder who participated in an experiment comparing subanesthetic ketamine (0.5 mg/kg) with a placebo saline infusion. Magnetoencephalographic recordings were collected at baseline and after ketamine and placebo infusions. During scanning, patients rested with their eyes closed. SI and depression were assessed, and a linear mixed-effects model was used to identify brain regions where gamma power and both SI and depression were associated. Two regions of the salience network (anterior insula, anterior cingulate) were then probed using dynamic causal modeling to test for ketamine effects.

RESULTS

Clinically, patients showed significantly reduced SI and depression after ketamine administration. In addition, distinct regions in the anterior insula were found to be associated with SI compared with depression. In modeling of insula-anterior cingulate connectivity, ketamine lowered the membrane capacitance for superficial pyramidal cells. Finally, connectivity between the insula and anterior cingulate was associated with improvements in depression symptoms.

CONCLUSIONS

These findings suggest that the anterior insula plays a key role in SI, perhaps via its role in salience detection. In addition, transient changes in superficial pyramidal cell membrane capacitance and subsequent increases in cortical excitability might be a mechanism through which ketamine improves SI.

Oscillatory, Computational, and Behavioral Evidence for Impaired GABAergic Inhibition in Schizophrenia

The dysconnection hypothesis of schizophrenia (SZ) proposes that psychosis is best understood in terms of aberrant connectivity. Specifically, it suggests that dysconnectivity arises through aberrant synaptic modulation associated with deficits in GABAergic inhibition, excitation-inhibition balance and disturbances of high-frequency oscillations. Using a computational model combined with a graded-difficulty visual orientation discrimination paradigm, we demonstrate that, in SZ, perceptual performance is determined by the balance of excitation-inhibition in superficial cortical layers. Twenty-eight individuals with a DSM-IV diagnosis of SZ, and 30 age- and gender-matched healthy controls participated in a psychophysics orientation discrimination task, a visual grating magnetoencephalography (MEG) recording, and a magnetic resonance spectroscopy (MRS) scan for GABA. Using a neurophysiologically informed model, we quantified group differences in GABA, gamma measures, and the predictive validity of model parameters for orientation discrimination in the SZ group. MEG visual gamma frequency was reduced in SZ, with lower peak frequency associated with more severe negative symptoms. Orientation discrimination performance was impaired in SZ. Dynamic causal modeling of the MEG data showed that local synaptic connections were reduced in SZ and local inhibition correlated negatively with the severity of negative symptoms. The effective connectivity between inhibitory interneurons and superficial pyramidal cells predicted orientation discrimination performance within the SZ group; consistent with graded, behaviorally relevant, disease-related changes in local GABAergic connections. Occipital GABA levels were significantly reduced in SZ but did not predict behavioral performance or oscillatory measures. These findings endorse the importance, and behavioral relevance, of GABAergic synaptic disconnection in schizophrenia that underwrites excitation-inhibition balance.

Dynamic causal modelling of phase-amplitude interactions

Models of coupled phase oscillators are used to describe a wide variety of phenomena in neuroimaging. These models typically rest on the premise that oscillator dynamics do not evolve beyond their respective limit cycles, and hence that interactions can be described purely in terms of phase differences. Whilst mathematically convenient, the restrictive nature of phase-only models can limit their explanatory power. We therefore propose a generalisation of dynamic causal modelling that incorporates both phase and amplitude. This allows for the separate quantifications of phase and amplitude contributions to the connectivity between neural regions. We show, using model-generated data and simulations of coupled pendula, that phase-amplitude models can describe strongly coupled systems more effectively than their phase-only counterparts. We relate our findings to four metrics commonly used in neuroimaging: the Kuramoto order parameter, cross-correlation, phase-lag index, and spectral entropy. We find that, with the exception of spectral entropy, the phase-amplitude model is able to capture all metrics more effectively than the phase-only model. We then demonstrate, using local field potential recordings in rodents and functional magnetic resonance imaging in macaque monkeys, that amplitudes in oscillator models play an important role in describing neural dynamics in anaesthetised brain states.

Neuronal network models of epileptogenesis

Epilepsy is a chronic neurological condition, following some trigger, transforming a normal brain to one that produces recurrent unprovoked seizures. In the search for the mechanisms that best explain the epileptogenic process, there is a growing body of evidence suggesting that the epilepsies are network level disorders. In this review, we briefly describe the concept of neuronal networks and highlight 2 methods used to analyse such networks. The first method, graph theory, is used to describe general characteristics of a network to facilitate comparison between normal and abnormal networks. The second, dynamic causal modelling, is useful in the analysis of the pathways of seizure spread. We concluded that the end results of the epileptogenic process are best understood as abnormalities of neuronal circuitry and not simply as molecular or cellular abnormalities. The network approach promises to generate new understanding and more targeted treatment of epilepsy.

Neurophysiological effects of continuous cortical stimulation in epilepsy - Spike and spontaneous ECoG activity

OBJECTIVE

The effect of continuous subthreshold cortical stimulation (CSCS) over the seizure onset zone (SOZ) in epilepsy was analyzed to delineate the affected physiological processes.

METHOD

ECoG data was recorded over SOZ and adjacent regions in patients (n = 7) with refractory-epilepsy. Data was reviewed before and during 2 Hz cortical electrical stimulation. Group differences were estimated using ANOVA and correlation with Pearson's r.

RESULTS

CSCS reduced background ECoG power at SOZ (p < 0.05), increased spectral coherence (p < 0.05) and reduced spike rate (p < 0.01) over all recorded sites. Spectral power and coherence (p < 0.01) correlated with spike rate at SOZ but not with each other at any location. Spike morphology correlated with spike-rate over all recorded sites (p < 0.0001) and with spectral power and coherence at SOZ (p < 0.01).

CONCLUSION

This study shows changes in cortical electrophysiology during CSCS over the SOZ where spike rate reduction correlated with two independent electrophysiological parameters, background power and coherence. These results suggest the possibility of a causal relationship between spectral power, coherence and interictal spikes which may be related to seizure rate.

SIGNIFICANCE

Improved understanding of the effect of electrical stimulation on epileptic tissue could suggest improvements in stimulation paradigms to reduce seizure frequency.

NMDA-receptor antibodies alter cortical microcircuit dynamics

Recently, autoantibodies against NMDA receptors (NMDARs) were identified as a major cause of autoimmune encephalitis. They cause abnormalities in brain function often associated with significant changes in patients’ brain dynamics. Here we use computational modeling to identify how NMDAR dysfunction causes abnormalities in brain dynamics using patient EEGs and local field potential recordings in a mouse model of NMDAR-Ab encephalitis. NMDAR autoantibodies cause a specific shift in excitatory coupling within cortical circuits that places the circuits closer to pathological transitions between dynamic brain states. Because of the proximity to these phase transitions, otherwise benign fluctuations in neuronal coupling cause abnormal EEG responses in the presence of the antibodies. Our modeling results thus explain fluctuating abnormalities in brain dynamics observed in patients.

NMDA-receptor antibodies (NMDAR-Abs) cause an autoimmune encephalitis with a diverse range of EEG abnormalities. NMDAR-Abs are believed to disrupt receptor function, but how blocking this excitatory synaptic receptor can lead to paroxysmal EEG abnormalities—or even seizures—is poorly understood. Here we show that NMDAR-Abs change intrinsic cortical connections and neuronal population dynamics to alter the spectral composition of spontaneous EEG activity and predispose brain dynamics to paroxysmal abnormalities. Based on local field potential recordings in a mouse model, we first validate a dynamic causal model of NMDAR-Ab effects on cortical microcircuitry. Using this model, we then identify the key synaptic parameters that best explain EEG paroxysms in pediatric patients with NMDAR-Ab encephalitis. Finally, we use the mouse model to show that NMDAR-Ab–related changes render microcircuitry critically susceptible to overt EEG paroxysms when these key parameters are changed, even though the same parameter fluctuations are tolerated in the in silico model of the control condition. These findings offer mechanistic insights into circuit-level dysfunction induced by NMDAR-Ab.

Generic dynamic causal modelling: An illustrative application to Parkinson's disease

We present a technical development in the dynamic causal modelling of electrophysiological responses that combines qualitatively different neural mass models within a single network. This affords the option to couple various cortical and subcortical nodes that differ in their form and dynamics. Moreover, it enables users to implement new neural mass models in a straightforward and standardized way. This generic framework hence supports flexibility and facilitates the exploration of increasingly plausible models. We illustrate this by coupling a basal ganglia-thalamus model to a (previously validated) cortical model developed specifically for motor cortex. The ensuing DCM is used to infer pathways that contribute to the suppression of beta oscillations induced by dopaminergic medication in patients with Parkinson's disease. Experimental recordings were obtained from deep brain stimulation electrodes (implanted in the subthalamic nucleus) and simultaneous magnetoencephalography. In line with previous studies, our results indicate a reduction of synaptic efficacy within the circuit between the subthalamic nucleus and external pallidum, as well as reduced efficacy in connections of the hyperdirect and indirect pathway leading to this circuit. This work forms the foundation for a range of modelling studies of the synaptic mechanisms (and pathophysiology) underlying event-related potentials and cross-spectral densities.

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.

Ion channels in EEG: isolating channel dysfunction in NMDA receptor antibody encephalitis

See Roberts and Breakspear (doi:10.1093/brain/awy136) for a scientific commentary on this article.Neurological and psychiatric practice frequently lack diagnostic probes that can assess mechanisms of neuronal communication non-invasively in humans. In N-methyl-d-aspartate (NMDA) receptor antibody encephalitis, functional molecular assays are particularly important given the presence of NMDA antibodies in healthy populations, the multifarious symptomology and the lack of radiological signs. Recent advances in biophysical modelling techniques suggest that inferring cellular-level properties of neural circuits from macroscopic measures of brain activity is possible. Here, we estimated receptor function from EEG in patients with NMDA receptor antibody encephalitis (n = 29) as well as from encephalopathic and neurological patient controls (n = 36). We show that the autoimmune patients exhibit distinct fronto-parietal network changes from which ion channel estimates can be obtained using a microcircuit model. Specifically, a dynamic causal model of EEG data applied to spontaneous brain responses identifies a selective deficit in signalling at NMDA receptors in patients with NMDA receptor antibody encephalitis but not at other ionotropic receptors. Moreover, though these changes are observed across brain regions, these effects predominate at the NMDA receptors of excitatory neurons rather than at inhibitory interneurons. Given that EEG is a ubiquitously available clinical method, our findings suggest a unique re-purposing of EEG data as an assay of brain network dysfunction at the molecular level.

Adaptive optimal input design and parametric estimation of nonlinear dynamical systems: application to neuronal modeling

OBJECTIVE

Many physical models of biological processes including neural systems are characterized by parametric nonlinear dynamical relations between driving inputs, internal states, and measured outputs of the process. Fitting such models using experimental data (data assimilation) is a challenging task since the physical process often operates in a noisy, possibly non-stationary environment; moreover, conducting multiple experiments under controlled and repeatable conditions can be impractical, time consuming or costly. The accuracy of model identification, therefore, is dictated principally by the quality and dynamic richness of collected data over single or few experimental sessions. Accordingly, it is highly desirable to design efficient experiments that, by exciting the physical process with smart inputs, yields fast convergence and increased accuracy of the model.

APPROACH

We herein introduce an adaptive framework in which optimal input design is integrated with square root cubature Kalman filters (OID-SCKF) to develop an online estimation procedure that first, converges significantly quicker, thereby permitting model fitting over shorter time windows, and second, enhances model accuracy when only few process outputs are accessible. The methodology is demonstrated on common nonlinear models and on a four-area neural mass model with noisy and limited measurements. Estimation quality (speed and accuracy) is benchmarked against high-performance SCKF-based methods that commonly employ dynamically rich informed inputs for accurate model identification.

MAIN RESULTS

For all the tested models, simulated single-trial and ensemble averages showed that OID-SCKF exhibited (i) faster convergence of parameter estimates and (ii) lower dependence on inter-trial noise variability with gains up to around 1000 ms in speed and 81% increase in variability for the neural mass models. In terms of accuracy, OID-SCKF estimation was superior, and exhibited considerably less variability across experiments, in identifying model parameters of (a) systems with challenging model inversion dynamics and (b) systems with fewer measurable outputs that directly relate to the underlying processes.

SIGNIFICANCE

Fast and accurate identification therefore carries particular promise for modeling of transient (short-lived) neuronal network dynamics using a spatially under-sampled set of noisy measurements, as is commonly encountered in neural engineering applications.

Generative models for clinical applications in computational psychiatry

Despite the success of modern neuroimaging techniques in furthering our understanding of cognitive and pathophysiological processes, translation of these advances into clinically relevant tools has been virtually absent until now. Neuromodeling represents a powerful framework for overcoming this translational deadlock, and the development of computational models to solve clinical problems has become a major scientific goal over the last decade, as reflected by the emergence of clinically oriented neuromodeling fields like Computational Psychiatry, Computational Neurology, and Computational Psychosomatics. Generative models of brain physiology and connectivity in the human brain play a key role in this endeavor, striving for computational assays that can be applied to neuroimaging data from individual patients for differential diagnosis and treatment prediction. In this review, we focus on dynamic causal modeling (DCM) and its use for Computational Psychiatry. DCM is a widely used generative modeling framework for functional magnetic resonance imaging (fMRI) and magneto-/electroencephalography (M/EEG) data. This article reviews the basic concepts of DCM, revisits examples where it has proven valuable for addressing clinically relevant questions, and critically discusses methodological challenges and recent methodological advances. We conclude this review with a more general discussion of the promises and pitfalls of generative models in Computational Psychiatry and highlight the path that lies ahead of us. This article is categorized under: Neuroscience > Computation Neuroscience > Clinical Neuroscience.

Hybrid Cubature Kalman filtering for identifying nonlinear models from sampled recording: Estimation of neuronal dynamics

Kalman filtering methods have long been regarded as efficient adaptive Bayesian techniques for estimating hidden states in models of linear dynamical systems under Gaussian uncertainty. Recent advents of the Cubature Kalman filter (CKF) have extended this efficient estimation property to nonlinear systems, and also to hybrid nonlinear problems where by the processes are continuous and the observations are discrete (continuous-discrete CD-CKF). Employing CKF techniques, therefore, carries high promise for modeling many biological phenomena where the underlying processes exhibit inherently nonlinear, continuous, and noisy dynamics and the associated measurements are uncertain and time-sampled. This paper investigates the performance of cubature filtering (CKF and CD-CKF) in two flagship problems arising in the field of neuroscience upon relating brain functionality to aggregate neurophysiological recordings: (i) estimation of the firing dynamics and the neural circuit model parameters from electric potentials (EP) observations, and (ii) estimation of the hemodynamic model parameters and the underlying neural drive from BOLD (fMRI) signals. First, in simulated neural circuit models, estimation accuracy was investigated under varying levels of observation noise (SNR), process noise structures, and observation sampling intervals (dt). When compared to the CKF, the CD-CKF consistently exhibited better accuracy for a given SNR, sharp accuracy increase with higher SNR, and persistent error reduction with smaller dt. Remarkably, CD-CKF accuracy shows only a mild deterioration for non-Gaussian process noise, specifically with Poisson noise, a commonly assumed form of background fluctuations in neuronal systems. Second, in simulated hemodynamic models, parametric estimates were consistently improved under CD-CKF. Critically, time-localization of the underlying neural drive, a determinant factor in fMRI-based functional connectivity studies, was significantly more accurate under CD-CKF. In conclusion, and with the CKF recently benchmarked against other advanced Bayesian techniques, the CD-CKF framework could provide significant gains in robustness and accuracy when estimating a variety of biological phenomena models where the underlying process dynamics unfold at time scales faster than those seen in collected measurements.

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.

Multistability in Large Scale Models of Brain Activity

Noise driven exploration of a brain network’s dynamic repertoire has been hypothesized to be causally involved in cognitive function, aging and neurodegeneration. The dynamic repertoire crucially depends on the network’s capacity to store patterns, as well as their stability. Here we systematically explore the capacity of networks derived from human connectomes to store attractor states, as well as various network mechanisms to control the brain’s dynamic repertoire. Using a deterministic graded response Hopfield model with connectome-based interactions, we reconstruct the system’s attractor space through a uniform sampling of the initial conditions. Large fixed-point attractor sets are obtained in the low temperature condition, with a bigger number of attractors than ever reported so far. Different variants of the initial model, including (i) a uniform activation threshold or (ii) a global negative feedback, produce a similarly robust multistability in a limited parameter range. A numerical analysis of the distribution of the attractors identifies spatially-segregated components, with a centro-medial core and several well-delineated regional patches. Those different modes share similarity with the fMRI independent components observed in the “resting state” condition. We demonstrate non-stationary behavior in noise-driven generalizations of the models, with different meta-stable attractors visited along the same time course. Only the model with a global dynamic density control is found to display robust and long-lasting non-stationarity with no tendency toward either overactivity or extinction. The best fit with empirical signals is observed at the edge of multistability, a parameter region that also corresponds to the highest entropy of the attractors.

Recent developments in non-invasive brain imaging allow reconstructing axonal tracts in the human brain and building realistic network models of the human brain. These models resemble brain systems in their network character and allow deciphering how different regions share signals and process information. Inspired by the metastable dynamics of the spin glass model in statistical physics, we systematically explore the brain network’s capacity to process information and investigate novel avenues how to enhance it. In particular, we study how the brain activates and switches between different functional networks across time. Such non-stationary behavior has been observed in human brain imaging data and hypothesized to be linked to information processsing. To shed light on the conditions under which large-scale brain network models exhibit such dynamics, we characterize the principal network patterns and confront them with modular structures observed both in graph theoretical analysis and resting-state functional Magnetic Resonance Imaging (rs-fMRI).

Evidence that Subanesthetic Doses of Ketamine Cause Sustained Disruptions of NMDA and AMPA-Mediated Frontoparietal Connectivity in Humans

Following the discovery of the antidepressant properties of ketamine, there has been a recent resurgence in the interest in this NMDA receptor antagonist. Although detailed animal models of the molecular mechanisms underlying ketamine's effects have emerged, there are few MEG/EEG studies examining the acute subanesthetic effects of ketamine infusion in man. We recorded 275 channel MEG in two experiments (n = 25 human males) examining the effects of subanesthetic ketamine infusion. MEG power spectra revealed a rich set of significant oscillatory changes compared with placebo sessions, including decreases in occipital, parietal, and anterior cingulate alpha power, increases in medial frontal theta power, and increases in parietal and cingulate cortex high gamma power. Each of these spectral effects demonstrated their own set of temporal dynamics. Dynamic causal modeling of frontoparietal connectivity changes with ketamine indicated a decrease in NMDA and AMPA-mediated frontal-to-parietal connectivity. AMPA-mediated connectivity changes were sustained for up to 50 min after ketamine infusion had ceased, by which time perceptual distortions were absent. The results also indicated a decrease in gain of parietal pyramidal cells, which was correlated with participants' self-reports of blissful state. Based on these results, we suggest that the antidepressant effects of ketamine may depend on its ability to change the balance of frontoparietal connectivity patterns.

SIGNIFICANCE STATEMENT

In this paper, we found that subanesthetic doses of ketamine, similar to those used in antidepressant studies, increase anterior theta and gamma power but decrease posterior theta, delta, and alpha power, as revealed by magnetoencephalographic recordings. Dynamic causal modeling of frontoparietal connectivity changes with ketamine indicated a decrease in NMDA and AMPA-mediated frontal-to-parietal connectivity. AMPA-mediated connectivity changes were sustained for up to 50 min after ketamine infusion had ceased, by which time perceptual distortions were absent. The results also indicated a decrease in gain of parietal pyramidal cells, which was correlated with participants' self-reports of blissful state. The alterations in frontoparietal connectivity patterns we observe here may be important in generating the antidepressant response to ketamine.

Profiling neuronal ion channelopathies with non-invasive brain imaging and dynamic causal models: Case studies of single gene mutations

Clinical assessments of brain function rely upon visual inspection of electroencephalographic waveform abnormalities in tandem with functional magnetic resonance imaging. However, no current technology proffers in vivo assessments of activity at synapses, receptors and ion-channels, the basis of neuronal communication. Using dynamic causal modeling we compared electrophysiological responses from two patients with distinct monogenic ion channelopathies and a large cohort of healthy controls to demonstrate the feasibility of assaying synaptic-level channel communication non-invasively. Synaptic channel abnormality was identified in both patients (100% sensitivity) with assay specificity above 89%, furnishing estimates of neurotransmitter and voltage-gated ion throughput of sodium, calcium, chloride and potassium. This performance indicates a potential novel application as an adjunct for clinical assessments in neurological and psychiatric settings. More broadly, these findings indicate that biophysical models of synaptic channels can be estimated non-invasively, having important implications for advancing human neuroimaging to the level of non-invasive ion channel assays.

•

Dynamic causal modeling (DCM) for M/EEG includes ion channel parameter estimates.

•

Parameter estimates from patients with monogenic ion channelopathies were compared.

•

Synaptic channel abnormality was identified in patients, with specificity above 89%.

•

DCM could serve as a platform for non-invasively assaying brain molecular dynamics.

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.

Losing control under ketamine: suppressed cortico-hippocampal drive following acute ketamine in rats

Systemic doses of the psychotomimetic ketamine alter the spectral characteristics of hippocampal and prefrontal cortical network activity. Using dynamic causal modeling (DCM) of cross-spectral densities, we quantify the putative synaptic mechanisms underlying ketamine effects in terms of changes in directed, effective connectivity between dorsal hippocampus and medial prefrontal (dCA1-mPFC) cortex of freely moving rats. We parameterize dose-dependent changes in spectral signatures of dCA1-mPFC local field potential recordings, using neural mass models of glutamatergic and GABAergic circuits. Optimizing DCMs of theta and gamma frequency range responses, model comparisons suggest that both enhanced gamma and depressed theta power result from a reduction in top-down connectivity from mPFC to the hippocampus, mediated by postsynaptic NMDA receptors (NMDARs). This is accompanied by an alteration in the bottom-up pathway from dCA1 to mPFC, which exhibits a distinct asymmetry: here, feed-forward drive at AMPA receptors increases in the presence of decreased NMDAR-mediated inputs. Setting these findings in the context of predictive coding suggests that NMDAR antagonism by ketamine in recurrent hierarchical networks may result in the failure of top-down connections from higher cortical regions to signal predictions to lower regions in the hierarchy, which consequently fail to respond consistently to errors. Given that NMDAR dysfunction has a central role in pathophysiological theories of schizophrenia and that theta and gamma rhythm abnormalities are evident in schizophrenic patients, the approach followed here may furnish a framework for the study of aberrant hierarchical message passing (of prediction errors) in schizophrenia-and the false perceptual inferences that ensue.

The brain ages optimally to model its environment: evidence from sensory learning over the adult lifespan

The aging brain shows a progressive loss of neuropil, which is accompanied by subtle changes in neuronal plasticity, sensory learning and memory. Neurophysiologically, aging attenuates evoked responses--including the mismatch negativity (MMN). This is accompanied by a shift in cortical responsivity from sensory (posterior) regions to executive (anterior) regions, which has been interpreted as a compensatory response for cognitive decline. Theoretical neurobiology offers a simpler explanation for all of these effects--from a Bayesian perspective, as the brain is progressively optimized to model its world, its complexity will decrease. A corollary of this complexity reduction is an attenuation of Bayesian updating or sensory learning. Here we confirmed this hypothesis using magnetoencephalographic recordings of the mismatch negativity elicited in a large cohort of human subjects, in their third to ninth decade. Employing dynamic causal modeling to assay the synaptic mechanisms underlying these non-invasive recordings, we found a selective age-related attenuation of synaptic connectivity changes that underpin rapid sensory learning. In contrast, baseline synaptic connectivity strengths were consistently strong over the decades. Our findings suggest that the lifetime accrual of sensory experience optimizes functional brain architectures to enable efficient and generalizable predictions of the world.

Network diffusion accurately models the relationship between structural and functional brain connectivity networks

The relationship between anatomic connectivity of large-scale brain networks and their functional connectivity is of immense importance and an area of active research. Previous attempts have required complex simulations which model the dynamics of each cortical region, and explore the coupling between regions as derived by anatomic connections. While much insight is gained from these non-linear simulations, they can be computationally taxing tools for predicting functional from anatomic connectivities. Little attention has been paid to linear models. Here we show that a properly designed linear model appears to be superior to previous non-linear approaches in capturing the brain's long-range second order correlation structure that governs the relationship between anatomic and functional connectivities. We derive a linear network of brain dynamics based on graph diffusion, whereby the diffusing quantity undergoes a random walk on a graph. We test our model using subjects who underwent diffusion MRI and resting state fMRI. The network diffusion model applied to the structural networks largely predicts the correlation structures derived from their fMRI data, to a greater extent than other approaches. The utility of the proposed approach is that it can routinely be used to infer functional correlation from anatomic connectivity. And since it is linear, anatomic connectivity can also be inferred from functional data. The success of our model confirms the linearity of ensemble average signals in the brain, and implies that their long-range correlation structure may percolate within the brain via purely mechanistic processes enacted on its structural connectivity pathways.

Modeling ketamine effects on synaptic plasticity during the mismatch negativity

This paper presents a model-based investigation of mechanisms underlying the reduction of mismatch negativity (MMN) amplitudes under the NMDA-receptor antagonist ketamine. We applied dynamic causal modeling and Bayesian model selection to data from a recent ketamine study of the roving MMN paradigm, using a cross-over, double-blind, placebo-controlled design. Our modeling was guided by a predictive coding framework that unifies contemporary "adaptation" and "model adjustment" MMN theories. Comparing a series of dynamic causal models that allowed for different expressions of neuronal adaptation and synaptic plasticity, we obtained 3 major results: 1) We replicated previous results that both adaptation and short-term plasticity are necessary to explain MMN generation per se; 2) we found significant ketamine effects on synaptic plasticity, but not adaptation, and a selective ketamine effect on the forward connection from left primary auditory cortex to superior temporal gyrus; 3) this model-based estimate of ketamine effects on synaptic plasticity correlated significantly with ratings of ketamine-induced impairments in cognition and control. Our modeling approach thus suggests a concrete mechanism for ketamine effects on MMN that correlates with drug-induced psychopathology. More generally, this demonstrates the potential of modeling for inferring on synaptic physiology, and its pharmacological modulation, from electroencephalography data.

Neural masses and fields in dynamic causal modeling

Dynamic causal modeling (DCM) provides a framework for the analysis of effective connectivity among neuronal subpopulations that subtend invasive (electrocorticograms and local field potentials) and non-invasive (electroencephalography and magnetoencephalography) electrophysiological responses. This paper reviews the suite of neuronal population models including neural masses, fields and conductance-based models that are used in DCM. These models are expressed in terms of sets of differential equations that allow one to model the synaptic underpinnings of connectivity. We describe early developments using neural mass models, where convolution-based dynamics are used to generate responses in laminar-specific populations of excitatory and inhibitory cells. We show that these models, though resting on only two simple transforms, can recapitulate the characteristics of both evoked and spectral responses observed empirically. Using an identical neuronal architecture, we show that a set of conductance based models-that consider the dynamics of specific ion-channels-present a richer space of responses; owing to non-linear interactions between conductances and membrane potentials. We propose that conductance-based models may be more appropriate when spectra present with multiple resonances. Finally, we outline a third class of models, where each neuronal subpopulation is treated as a field; in other words, as a manifold on the cortical surface. By explicitly accounting for the spatial propagation of cortical activity through partial differential equations (PDEs), we show that the topology of connectivity-through local lateral interactions among cortical layers-may be inferred, even in the absence of spatially resolved data. We also show that these models allow for a detailed analysis of structure-function relationships in the cortex. Our review highlights the relationship among these models and how the hypothesis asked of empirical data suggests an appropriate model class.

Combined neurostimulation and neuroimaging in cognitive neuroscience: past, present, and future

Modern neurostimulation approaches in humans provide controlled inputs into the operations of cortical regions, with highly specific behavioral consequences. This enables causal structure–function inferences, and in combination with neuroimaging, has provided novel insights into the basic mechanisms of action of neurostimulation on distributed networks. For example, more recent work has established the capacity of transcranial magnetic stimulation (TMS) to probe causal interregional influences, and their interaction with cognitive state changes. Combinations of neurostimulation and neuroimaging now face the challenge of integrating the known physiological effects of neurostimulation with theoretical and biological models of cognition, for example, when theoretical stalemates between opposing cognitive theories need to be resolved. This will be driven by novel developments, including biologically informed computational network analyses for predicting the impact of neurostimulation on brain networks, as well as novel neuroimaging and neurostimulation techniques. Such future developments may offer an expanded set of tools with which to investigate structure–function relationships, and to formulate and reconceptualize testable hypotheses about complex neural network interactions and their causal roles in cognition.

DCM for complex-valued data: cross-spectra, coherence and phase-delays

This note describes an extension of Bayesian model inversion procedures for the Dynamic Causal Modeling (DCM) of complex-valued data. Modeling complex data can be particularly useful in the analysis of multivariate ergodic (stationary) time-series. We illustrate this with a generalization of DCM for steady-state responses that models both the real and imaginary parts of sample cross-spectra. DCM allows one to infer underlying biophysical parameters generating data (like synaptic time constants, connection strengths and conduction delays). Because transfer functions and complex cross-spectra can be generated from these parameters, one can also describe the implicit system architecture in terms of conventional (linear systems) measures; like coherence, phase-delay or cross-correlation functions. Crucially, these measures can be derived in both sensor and source-space. In other words, one can examine the cross-correlation or phase-delay functions between hidden neuronal sources using non-invasive data and relate these functions to synaptic parameters and neuronal conduction delays. We illustrate these points using local field potential recordings from the subthalamic nucleus and globus pallidus, with a special focus on the relationship between conduction delays and the ensuing phase relationships and cross-correlation time lags between population activities.

Dynamic causal models and physiological inference: a validation study using isoflurane anaesthesia in rodents

Generative models of neuroimaging and electrophysiological data present new opportunities for accessing hidden or latent brain states. Dynamic causal modeling (DCM) uses Bayesian model inversion and selection to infer the synaptic mechanisms underlying empirically observed brain responses. DCM for electrophysiological data, in particular, aims to estimate the relative strength of synaptic transmission at different cell types and via specific neurotransmitters. Here, we report a DCM validation study concerning inference on excitatory and inhibitory synaptic transmission, using different doses of a volatile anaesthetic agent (isoflurane) to parametrically modify excitatory and inhibitory synaptic processing while recording local field potentials (LFPs) from primary auditory cortex (A1) and the posterior auditory field (PAF) in the auditory belt region in rodents. We test whether DCM can infer, from the LFP measurements, the expected drug-induced changes in synaptic transmission mediated via fast ionotropic receptors; i.e., excitatory (glutamatergic) AMPA and inhibitory GABA(A) receptors. Cross- and auto-spectra from the two regions were used to optimise three DCMs based on biologically plausible neural mass models and specific network architectures. Consistent with known extrinsic connectivity patterns in sensory hierarchies, we found that a model comprising forward connections from A1 to PAF and backward connections from PAF to A1 outperformed a model with forward connections from PAF to A1 and backward connections from A1 to PAF and a model with reciprocal lateral connections. The parameter estimates from the most plausible model indicated that the amplitude of fast glutamatergic excitatory postsynaptic potentials (EPSPs) and inhibitory postsynaptic potentials (IPSPs) behaved as predicted by previous neurophysiological studies. Specifically, with increasing levels of anaesthesia, glutamatergic EPSPs decreased linearly, whereas fast GABAergic IPSPs displayed a nonlinear (saturating) increase. The consistency of our model-based in vivo results with experimental in vitro results lends further validity to the capacity of DCM to infer on synaptic processes using macroscopic neurophysiological data.

An in vivo assay of synaptic function mediating human cognition

The contribution of dopamine to working memory has been studied extensively [1-3]. Here, we exploited its well characterized effects [1-3] to validate a novel human in vivo assay of ongoing synaptic [4, 5] processing. We obtained magnetoencephalographic (MEG) measurements from subjects performing a working memory (WM) task during a within-subject, placebo-controlled, pharmacological (dopaminergic) challenge. By applying dynamic causal modeling (DCM), a Bayesian technique for neuronal system identification [6], to MEG signals from prefrontal cortex, we demonstrate that it is possible to infer synaptic signaling by specific ion channels in behaving humans. Dopamine-induced enhancement of WM performance was accompanied by significant changes in MEG signal power, and a DCM assay disclosed related changes in synaptic signaling. By estimating the contribution of ionotropic receptors (AMPA, NMDA, and GABA(A)) to the observed spectral response, we demonstrate changes in their function commensurate with the synaptic effects of dopamine. The validity of our model is reinforced by a striking quantitative effect on NMDA and AMPA receptor signaling that predicted behavioral improvement over subjects. Our results provide a proof-of-principle demonstration of a novel framework for inferring, noninvasively, neuromodulatory influences on ion channel signaling via specific ionotropic receptors, providing a window on the hidden synaptic events mediating discrete psychological processes in humans.