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

NORSE/FIRES: how can we advance our understanding of this devastating condition?

Introduction

New onset refractory status epilepticus (NORSE) is a rare and devastating condition characterised by the sudden onset of refractory status epilepticus (RSE) without an identifiable acute or active structural, toxic, or metabolic cause in an individual without a pre-existing diagnosis of epilepsy. Febrile infection-related epilepsy syndrome (FIRES) is considered a subcategory of NORSE and presents following a febrile illness prior to seizure onset. NORSE/FIRES is associated with high morbidity and mortality in children and adults.

Methods and results

In this review we first briefly summarise the reported clinical, paraclinical, treatment and outcome data in the literature. We then report on existing knowledge of the underlying pathophysiology in relation to in vitro and in vivo pre-clinical seizure and epilepsy models of potential relevance to NORSE/FIRES.

Discussion

We highlight how pre-clinical models can enhance our understanding of FIRES/NORSE and propose future directions for research.

Global dynamics of neural mass models

Neural mass models are used to simulate cortical dynamics and to explain the electrical and magnetic fields measured using electro- and magnetoencephalography. Simulations evince a complex phase-space structure for these kinds of models; including stationary points and limit cycles and the possibility for bifurcations and transitions among different modes of activity. This complexity allows neural mass models to describe the itinerant features of brain dynamics. However, expressive, nonlinear neural mass models are often difficult to fit to empirical data without additional simplifying assumptions: e.g., that the system can be modelled as linear perturbations around a fixed point. In this study we offer a mathematical analysis of neural mass models, specifically the canonical microcircuit model, providing analytical solutions describing slow changes in the type of cortical activity, i.e. dynamical itinerancy. We derive a perturbation analysis up to second order of the phase flow, together with adiabatic approximations. This allows us to describe amplitude modulations in a relatively simple mathematical format providing analytic proof-of-principle for the existence of semi-stable states of cortical dynamics at the scale of a cortical column. This work allows for model inversion of neural mass models, not only around fixed points, but over regions of phase space that encompass transitions among semi or multi-stable states of oscillatory activity. Crucially, these theoretical results speak to model inversion in the context of multiple semi-stable brain states, such as the transition between interictal, pre-ictal and ictal activity in epilepsy.

Fully Bayesian estimation of virtual brain parameters with self-tuning Hamiltonian Monte Carlo

Virtual brain models are data-driven patient-specific brain models integrating individual brain imaging data with neural mass modeling in a single computational framework, capable of autonomously generating brain activity and its associated brain imaging signals. Along the example of epilepsy, we develop an efficient and accurate Bayesian methodology estimating the parameters linked to the extent of the epileptogenic zone. State-of-the-art advances in Bayesian inference using Hamiltonian Monte Carlo (HMC) algorithms have remained elusive for large-scale differential-equations based models due to their slow convergence. We propose appropriate priors and a novel reparameterization to facilitate efficient exploration of the posterior distribution in terms of computational time and convergence diagnostics. The methodology is illustrated for in-silico dataset and then, applied to infer the personalized model parameters based on the empirical stereotactic electroencephalography (SEEG) recordings of retrospective patients. This improved methodology may pave the way to render HMC methods sufficiently easy and efficient to use, thus applicable in personalized medicine.

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.

Neural Decoding of EEG Signals with Machine Learning: A Systematic Review

Electroencephalography (EEG) is a non-invasive technique used to record the brain's evoked and induced electrical activity from the scalp. Artificial intelligence, particularly machine learning (ML) and deep learning (DL) algorithms, are increasingly being applied to EEG data for pattern analysis, group membership classification, and brain-computer interface purposes. This study aimed to systematically review recent advances in ML and DL supervised models for decoding and classifying EEG signals. Moreover, this article provides a comprehensive review of the state-of-the-art techniques used for EEG signal preprocessing and feature extraction. To this end, several academic databases were searched to explore relevant studies from the year 2000 to the present. Our results showed that the application of ML and DL in both mental workload and motor imagery tasks has received substantial attention in recent years. A total of 75% of DL studies applied convolutional neural networks with various learning algorithms, and 36% of ML studies achieved competitive accuracy by using a support vector machine algorithm. Wavelet transform was found to be the most common feature extraction method used for all types of tasks. We further examined the specific feature extraction methods and end classifier recommendations discovered in this systematic review.

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.

On the influence of prior information evaluated by fully Bayesian criteria in a personalized whole-brain model of epilepsy spread

Individualized anatomical information has been used as prior knowledge in Bayesian inference paradigms of whole-brain network models. However, the actual sensitivity to such personalized information in priors is still unknown. In this study, we introduce the use of fully Bayesian information criteria and leave-one-out cross-validation technique on the subject-specific information to assess different epileptogenicity hypotheses regarding the location of pathological brain areas based on a priori knowledge from dynamical system properties. The Bayesian Virtual Epileptic Patient (BVEP) model, which relies on the fusion of structural data of individuals, a generative model of epileptiform discharges, and a self-tuning Monte Carlo sampling algorithm, is used to infer the spatial map of epileptogenicity across different brain areas. Our results indicate that measuring the out-of-sample prediction accuracy of the BVEP model with informative priors enables reliable and efficient evaluation of potential hypotheses regarding the degree of epileptogenicity across different brain regions. In contrast, while using uninformative priors, the information criteria are unable to provide strong evidence about the epileptogenicity of brain areas. We also show that the fully Bayesian criteria correctly assess different hypotheses about both structural and functional components of whole-brain models that differ across individuals. The fully Bayesian information-theory based approach used in this study suggests a patient-specific strategy for epileptogenicity hypothesis testing in generative brain network models of epilepsy to improve surgical outcomes.

Altered EEG markers of synaptic plasticity in a human model of NMDA receptor deficiency: Anti-NMDA receptor encephalitis

Plasticity of synaptic strength and density is a vital mechanism enabling memory consolidation, learning, and neurodevelopment. It is strongly dependent on the intact function of N-Methyl-d-Aspartate Receptors (NMDAR). The importance of NMDAR is further evident as their dysfunction is involved in many diseases such as schizophrenia, Alzheimer's disease, neurodevelopmental disorders, and epilepsies. Synaptic plasticity is thought to be reflected by changes of sleep slow wave slopes across the night, namely higher slopes after wakefulness at the beginning of sleep than after a night of sleep. Hence, a functional NMDAR deficiency should theoretically lead to altered overnight changes of slow wave slopes. Here we investigated whether pediatric patients with anti-NMDAR encephalitis, being a very rare but unique human model of NMDAR deficiency due to autoantibodies against receptor subunits, indeed show alterations in this sleep EEG marker for synaptic plasticity. We retrospectively analyzed 12 whole-night EEGs of 9 patients (age 4.3-20.8 years, 7 females) and compared them to a control group of 45 healthy individuals with the same age distribution. Slow wave slopes were calculated for the first and last hour of Non-Rapid Eye Movement (NREM) sleep (factor 'hour') for patients and controls (factor 'group'). There was a significant interaction between 'hour' and 'group' (p = 0.013), with patients showing a smaller overnight decrease of slow wave slopes than controls. Moreover, we found smaller slopes during the first hour in patients (p = 0.022), whereas there was no group difference during the last hour of NREM sleep (p = 0.980). Importantly, the distribution of sleep stages was not different between the groups, and in our main analyses of patients without severe disturbance of sleep architecture, neither was the incidence of slow waves. These possible confounders could therefore not account for the differences in the slow wave slope values, which we also saw in the analysis of the whole sample of EEGs. These results suggest that quantitative EEG analysis of slow wave characteristics may reveal impaired synaptic plasticity in patients with anti-NMDAR encephalitis, a human model of functional NMDAR deficiency. Thus, in the future, the changes of sleep slow wave slopes may contribute to the development of electrophysiological biomarkers of functional NMDAR deficiency and synaptic plasticity in general.

TAPAS: An Open-Source Software Package for Translational Neuromodeling and Computational Psychiatry

Psychiatry faces fundamental challenges with regard to mechanistically guided differential diagnosis, as well as prediction of clinical trajectories and treatment response of individual patients. This has motivated the genesis of two closely intertwined fields: (i) Translational Neuromodeling (TN), which develops "computational assays" for inferring patient-specific disease processes from neuroimaging, electrophysiological, and behavioral data; and (ii) Computational Psychiatry (CP), with the goal of incorporating computational assays into clinical decision making in everyday practice. In order to serve as objective and reliable tools for clinical routine, computational assays require end-to-end pipelines from raw data (input) to clinically useful information (output). While these are yet to be established in clinical practice, individual components of this general end-to-end pipeline are being developed and made openly available for community use. In this paper, we present the Translational Algorithms for Psychiatry-Advancing Science (TAPAS) software package, an open-source collection of building blocks for computational assays in psychiatry. Collectively, the tools in TAPAS presently cover several important aspects of the desired end-to-end pipeline, including: (i) tailored experimental designs and optimization of measurement strategy prior to data acquisition, (ii) quality control during data acquisition, and (iii) artifact correction, statistical inference, and clinical application after data acquisition. Here, we review the different tools within TAPAS and illustrate how these may help provide a deeper understanding of neural and cognitive mechanisms of disease, with the ultimate goal of establishing automatized pipelines for predictions about individual patients. We hope that the openly available tools in TAPAS will contribute to the further development of TN/CP and facilitate the translation of advances in computational neuroscience into clinically relevant computational assays.

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.

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.

Automated diagnosis of encephalitis in pediatric patients using EEG rhythms and slow biphasic complexes

Slow biphasic complexes (SBC) have been identified in the EEG of patients suffering for inflammatory brain diseases. Their amplitude, location and frequency of appearance were found to correlate with the severity of encephalitis. Other characteristics of SBCs and of EEG traces of patients could reflect the grade of pathology. Here, EEG rhythms are investigated together with SBCs for a better characterization of encephalitis. EEGs have been acquired from pediatric patients: ten controls and ten encephalitic patients. They were split by neurologists into five classes of different severity of the pathology. The relative power of EEG rhythms was found to change significantly in EEGs labeled with different severity scores. Moreover, a significant variation was found in the last seconds before the appearance of an SBC. This information and quantitative indexes characterizing the SBCs were used to build a binary classification decision tree able to identify the classes of severity. True classification rate of the best model was 76.1% (73.5% with leave-one-out test). Moreover, the classification errors were among classes with similar severity scores (precision higher than 80% was achieved considering three instead of five classes). Our classification method may be a promising supporting tool for clinicians to diagnose, assess and make the follow-up of patients with encephalitis.

Dynamic Causal Modeling of the Relationship between Cognition and Theta-alpha Oscillations in Adults with Down Syndrome

Individuals with Down syndrome (DS) show high inter-subject variability in cognitive ability and have an ultra-high risk of developing dementia (90% lifetime prevalence). Elucidating factors underlying variability in cognitive function can inform us about intellectual disability (ID) and may improve our understanding of factors associated with later cognitive decline. Increased neuronal inhibition has been posited to contribute to ID in DS. Combining electroencephalography (EEG) with dynamic causal modeling (DCM) provides a non-invasive method for investigating excitatory/inhibitory mechanisms. Resting-state EEG recordings were obtained from 36 adults with DS with no evidence of cognitive decline. Theta–alpha activity (4–13 Hz) was characterized in relation to general cognitive ability (raw Kaufmann’s Brief Intelligence Test second Edition (KBIT-2) score). Higher KBIT-2 was associated with higher frontal alpha peak amplitude and higher theta–alpha band power across distributed regions. Modeling this association with DCM revealed intrinsic self-inhibition was the key network parameter underlying observed differences in 4–13 Hz power in relation to KBIT-2 and age. In particular, intrinsic self-inhibition in right V1 was negatively correlated with KBIT-2. Results suggest intrinsic self-inhibition within the alpha network is associated with individual differences in cognitive ability in adults with DS, and may provide a potential therapeutic target for cognitive enhancement.

Characterization and Outcomes of Epileptic Seizures in Mexican Pediatric Patients With Anti-N-Methyl-D-Aspartate Receptor Encephalitis

Introduction

Anti-N-methyl-D-aspartate receptor (NMDAR) encephalitis is one of the most common autoimmune encephalitides. The frequency of anti-NMDAR encephalitis is known to exceed the frequency of any individual viral encephalitis in young subjects. Epileptic seizures are a cardinal symptom in anti-NMDAR encephalitis; a significant amount of pediatric patients exhibit seizures as the first symptom of the disease, and most of them will develop them during the acute phase. The use of antiepileptic drugs (AEDs) is a cornerstone of the treatment of these patients, but the choice of agent and duration of treatment is currently unknown.

Materials and methods

This was a single-center retrospective review case series of all pediatric patients with a confirmed diagnosis of anti-NMDAR encephalitis and epileptic seizures admitted to the National Institute of Pediatrics in Mexico City from January 2012 to July 2019.

Results

We included a total of 31 patients (males 64.5%, median age: 10 years). No patient showed evidence of teratoma; only 38% of cases had a viral prodrome. Most patients initially exhibited psychiatric symptoms (51%), but the leading cause in soliciting medical assistance was the presence of epileptic seizures (71%). About 85% of patients presented epileptic seizures during the course of the illness, predominantly focal onset seizures (42% focal to bilateral tonic-clonic seizures, 32% focal seizures with impaired awareness). Electroencephalogram (EEG) was abnormal in 97% of patients; the characteristic extreme delta brush pattern was found in 9% of patients. Two AEDs on average were required to control seizures during the acute stage. In six (19%) patients, human herpesvirus (HHV) was detected in cerebrospinal fluid (CSF); all of them had epileptic seizures, which were more resistant to pharmacological treatment during the acute phase, requiring a higher number of AED (median 2.5 vs. 2). The development of epilepsy after acute encephalitis was uncommon; at 24 months, only one patient continued to have epileptic seizures. One of the factors most closely related to the persistence of epileptic seizures was the inadequate response to immunotherapy after four weeks. The functional prognosis was generally good; at a two-year follow-up, only two (10%) patients had a significant disability [modified Rankin Scale (mRS) score: 3-5]; both patients had seizures at a one-year follow-up.

Conclusions

Sustained use of AEDs after the acute phase of anti-NMDAR encephalitis is controversial. We found that the continuation of AEDs after the acute phase could be considered in the following scenarios: status epilepticus (SE), inadequate response to immunotherapy at four weeks, and a high mRS score at discharge and during follow-up. In other cases, discontinuation of AED may be warranted. More studies are needed in our country to replicate these results.

Experimental Models of Neuroimmunological Disorders: A Review

Immune-mediated inflammatory diseases of the central nervous system (CNS) are a group of neurological disorders in which inflammation and/or demyelination are induced by cellular and humoral immune responses specific to CNS antigens. They include diseases such as multiple sclerosis (MS), neuromyelitis optica spectrum disorders (NMOSD), acute disseminated encephalomyelitis (ADEM) and anti-NMDA receptor encephalitis (NMDAR encephalitis). Over the years, many in vivo and in vitro models were used to study clinical, pathological, physiological and immunological features of these neuroimmunological disorders. Nevertheless, there are important aspects of human diseases that are not fully reproduced in the experimental models due to their technical limitations. In this review, we describe the preclinical models of neuroimmune disorders, and how they contributed to the understanding of these disorders and explore potential treatments. We also describe the purpose and limitation of each one, as well as the recent advances in this field.

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.

Dynamic effective connectivity

Metastability is a key source of itinerant dynamics in the brain; namely, spontaneous spatiotemporal reorganization of neuronal activity. This itinerancy has been the focus of numerous dynamic functional connectivity (DFC) analyses - developed to characterize the formation and dissolution of distributed functional patterns over time, using resting state fMRI. However, aside from technical and practical controversies, these approaches cannot recover the neuronal mechanisms that underwrite itinerant (e.g., metastable) dynamics-due to their descriptive, model-free nature. We argue that effective connectivity (EC) analyses are more apt for investigating the neuronal basis of metastability. To this end, we appeal to biologically-grounded models (i.e., dynamic causal modelling, DCM) and dynamical systems theory (i.e., heteroclinic sequential dynamics) to create a probabilistic, generative model of haemodynamic fluctuations. This model generates trajectories in the parametric space of EC modes (i.e., states of connectivity) that characterize functional brain architectures. In brief, it extends an established spectral DCM, to generate functional connectivity data features that change over time. This foundational paper tries to establish the model's face validity by simulating non-stationary fMRI time series and recovering key model parameters (i.e., transition probabilities among connectivity states and the parametric nature of these states) using variational Bayes. These data are further characterized using Bayesian model comparison (within and between subjects). Finally, we consider practical issues that attend applications and extensions of this scheme. Importantly, the scheme operates within a generic Bayesian framework - that can be adapted to study metastability and itinerant dynamics in any non-stationary time series.

Mechanism of visual network dysfunction in relapsing-remitting multiple sclerosis and its relation to cognition

OBJECTIVE

To investigate if changes in brain network function and connectivity contribute to the abnormalities in visual event related potentials (ERP) in relapsing-remitting multiple sclerosis (RRMS), and explore their relation to a decrease in cognitive performance.

METHODS

We evaluated 72 patients with RRMS and 89 healthy control subjects in a cross-sectional study. Visual ERP were generated using illusory and non-illusory stimuli and recorded using 21 EEG scalp electrodes. The measured activity was modelled using Dynamic Causal Modelling. The model network consisted of 4 symmetric nodes including the primary visual cortex (V1/V2) and the Lateral Occipital Complex. Patients and controls were tested with a neuropsychological test battery consisting of 18 cognitive tests covering six cognitive domains.

RESULTS

We found reduced cortical connectivity in bottom-up and interhemispheric connections to the right lateral occipital complex in patients (p < 0.001). Furthermore, interhemispherical connections were related to cognitive dysfunction in several domains (attention, executive function, visual perception and organization, processing speed and global cognition) for patients (p < 0.05). No relation was seen between cortical network connectivity and cognitive function in the healthy control subjects.

CONCLUSION

Changes in the functional connectivity to higher cortical regions provide a neurobiological explanation for the changes of the visual ERP in RRMS.

SIGNIFICANCE

This study suggests that changes in connectivity to higher cortical regions partly explain visual network dysfunction in RRMS where a lower interhemispheric connectivity may contribute to impaired cognitive function.

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.

Computational Neuropsychology and Bayesian Inference

Computational theories of brain function have become very influential in neuroscience. They have facilitated the growth of formal approaches to disease, particularly in psychiatric research. In this paper, we provide a narrative review of the body of computational research addressing neuropsychological syndromes, and focus on those that employ Bayesian frameworks. Bayesian approaches to understanding brain function formulate perception and action as inferential processes. These inferences combine ‘prior’ beliefs with a generative (predictive) model to explain the causes of sensations. Under this view, neuropsychological deficits can be thought of as false inferences that arise due to aberrant prior beliefs (that are poor fits to the real world). This draws upon the notion of a Bayes optimal pathology – optimal inference with suboptimal priors – and provides a means for computational phenotyping. In principle, any given neuropsychological disorder could be characterized by the set of prior beliefs that would make a patient’s behavior appear Bayes optimal. We start with an overview of some key theoretical constructs and use these to motivate a form of computational neuropsychology that relates anatomical structures in the brain to the computations they perform. Throughout, we draw upon computational accounts of neuropsychological syndromes. These are selected to emphasize the key features of a Bayesian approach, and the possible types of pathological prior that may be present. They range from visual neglect through hallucinations to autism. Through these illustrative examples, we review the use of Bayesian approaches to understand the link between biology and computation that is at the heart of neuropsychology.

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.

Allostatic Self-efficacy: A Metacognitive Theory of Dyshomeostasis-Induced Fatigue and Depression

This paper outlines a hierarchical Bayesian framework for interoception, homeostatic/allostatic control, and meta-cognition that connects fatigue and depression to the experience of chronic dyshomeostasis. Specifically, viewing interoception as the inversion of a generative model of viscerosensory inputs allows for a formal definition of dyshomeostasis (as chronically enhanced surprise about bodily signals, or, equivalently, low evidence for the brain's model of bodily states) and allostasis (as a change in prior beliefs or predictions which define setpoints for homeostatic reflex arcs). Critically, we propose that the performance of interoceptive-allostatic circuitry is monitored by a metacognitive layer that updates beliefs about the brain's capacity to successfully regulate bodily states (allostatic self-efficacy). In this framework, fatigue and depression can be understood as sequential responses to the interoceptive experience of dyshomeostasis and the ensuing metacognitive diagnosis of low allostatic self-efficacy. While fatigue might represent an early response with adaptive value (cf. sickness behavior), the experience of chronic dyshomeostasis may trigger a generalized belief of low self-efficacy and lack of control (cf. learned helplessness), resulting in depression. This perspective implies alternative pathophysiological mechanisms that are reflected by differential abnormalities in the effective connectivity of circuits for interoception and allostasis. We discuss suitably extended models of effective connectivity that could distinguish these connectivity patterns in individual patients and may help inform differential diagnosis of fatigue and depression in the future.

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.

Computational neuroimaging strategies for single patient predictions

Neuroimaging increasingly exploits machine learning techniques in an attempt to achieve clinically relevant single-subject predictions. An alternative to machine learning, which tries to establish predictive links between features of the observed data and clinical variables, is the deployment of computational models for inferring on the (patho)physiological and cognitive mechanisms that generate behavioural and neuroimaging responses. This paper discusses the rationale behind a computational approach to neuroimaging-based single-subject inference, focusing on its potential for characterising disease mechanisms in individual subjects and mapping these characterisations to clinical predictions. Following an overview of two main approaches - Bayesian model selection and generative embedding - which can link computational models to individual predictions, we review how these methods accommodate heterogeneity in psychiatric and neurological spectrum disorders, help avoid erroneous interpretations of neuroimaging data, and establish a link between a mechanistic, model-based approach and the statistical perspectives afforded by machine learning.

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.