Establishing brain states in neuroimaging data.

The definition of a brain state remains elusive, with varying interpretations across different sub-fields of neuroscience-from the level of wakefulness in anaesthesia, to activity of individual neurons, voltage in EEG, and blood flow in fMRI. This lack of consensus presents a significant challenge to the development of accurate models of neural dynamics. However, at the foundation of dynamical systems theory lies a definition of what constitutes the 'state' of a system-i.e., a specification of the system's future. Here, we propose to adopt this definition to establish brain states in neuroimaging timeseries by applying Dynamic Causal Modelling (DCM) to low-dimensional embedding of resting and task condition fMRI data. We find that ~90% of subjects in resting conditions are better described by first-order models, whereas ~55% of subjects in task conditions are better described by second-order models. Our work calls into question the status quo of using first-order equations almost exclusively within computational neuroscience and provides a new way of establishing brain states, as well as their associated phase space representations, in neuroimaging datasets.

Microbiota links to neural dynamics supporting threat processing.

There is growing recognition that the composition of the gut microbiota influences behaviour, including responses to threat. The cognitive-interoceptive appraisal of threat-related stimuli relies on dynamic neural computations between the anterior insular (AIC) and the dorsal anterior cingulate (dACC) cortices. If, to what extent, and how microbial consortia influence the activity of this cortical threat processing circuitry is unclear. We addressed this question by combining a threat processing task, neuroimaging, 16S rRNA profiling and computational modelling in healthy participants. Results showed interactions between high-level ecological indices with threat-related AIC-dACC neural dynamics. At finer taxonomic resolutions, the abundance of Ruminococcus was differentially linked to connectivity between, and activity within the AIC and dACC during threat updating. Functional inference analysis provides a strong rationale to motivate future investigations of microbiota-derived metabolites in the observed relationship with threat-related brain processes.

Estimating anisotropy directly via neural timeseries.

An isotropic dynamical system is one that looks the same in every direction, i.e., if we imagine standing somewhere within an isotropic system, we would not be able to differentiate between different lines of sight. Conversely, anisotropy is a measure of the extent to which a system deviates from perfect isotropy, with larger values indicating greater discrepancies between the structure of the system along its axes. Here, we derive the form of a generalised scalable (mechanically similar) discretized field theoretic Lagrangian that allows for levels of anisotropy to be directly estimated via timeseries of arbitrary dimensionality. We generate synthetic data for both isotropic and anisotropic systems and, by using Bayesian model inversion and reduction, show that we can discriminate between the two datasets - thereby demonstrating proof of principle. We then apply this methodology to murine calcium imaging data collected in rest and task states, showing that anisotropy can be estimated directly from different brain states and cortical regions in an empirical in vivo biological setting. We hope that this theoretical foundation, together with the methodology and publicly available MATLAB code, will provide an accessible way for researchers to obtain new insight into the structural organization of neural systems in terms of how scalable neural regions grow - both ontogenetically during the development of an individual organism, as well as phylogenetically across species.

Model-based prediction of muscarinic receptor function from auditory mismatch negativity responses.

Drugs affecting neuromodulation, for example by dopamine or acetylcholine, take centre stage among therapeutic strategies in psychiatry. These neuromodulators can change both neuronal gain and synaptic plasticity and therefore affect electrophysiological measures. An important goal for clinical diagnostics is to exploit this effect in the reverse direction, i.e., to infer the status of specific neuromodulatory systems from electrophysiological measures. In this study, we provide proof-of-concept that the functional status of cholinergic (specifically muscarinic) receptors can be inferred from electrophysiological data using generative (dynamic causal) models. To this end, we used epidural EEG recordings over two auditory cortical regions during a mismatch negativity (MMN) paradigm in rats. All animals were treated, across sessions, with muscarinic receptor agonists and antagonists at different doses. Together with a placebo condition, this resulted in five levels of muscarinic receptor status. Using a dynamic causal model - embodying a small network of coupled cortical microcircuits - we estimated synaptic parameters and their change across pharmacological conditions. The ensuing parameter estimates associated with (the neuromodulation of) synaptic efficacy showed both graded muscarinic effects and predictive validity between agonistic and antagonistic pharmacological conditions. This finding illustrates the potential utility of generative models of electrophysiological data as computational assays of muscarinic function. In application to EEG data of patients from heterogeneous spectrum diseases, e.g. schizophrenia, such models might help identify subgroups of patients that respond differentially to cholinergic treatments. SIGNIFICANCE STATEMENT: In psychiatry, the vast majority of pharmacological treatments affect actions of neuromodulatory transmitters, e.g. dopamine or acetylcholine. As treatment is largely trial-and-error based, one of the goals for computational psychiatry is to construct mathematical models that can serve as "computational assays" and infer the status of specific neuromodulatory systems in individual patients. This translational neuromodeling strategy has great promise for electrophysiological data in particular but requires careful validation. The present study demonstrates that the functional status of cholinergic (muscarinic) receptors can be inferred from electrophysiological data using dynamic causal models of neural circuits. While accuracy needs to be enhanced and our results must be replicated in larger samples, our current results provide proof-of-concept for computational assays of muscarinic function using EEG.

Conservation laws by virtue of scale symmetries in neural systems.

In contrast to the symmetries of translation in space, rotation in space, and translation in time, the known laws of physics are not universally invariant under transformation of scale. However, a special case exists in which the action is scale invariant if it satisfies the following two constraints: 1) it must depend upon a scale-free Lagrangian, and 2) the Lagrangian must change under scale in the same way as the inverse time, [Formula: see text]. Our contribution lies in the derivation of a generalised Lagrangian, in the form of a power series expansion, that satisfies these constraints. This generalised Lagrangian furnishes a normal form for dynamic causal models-state space models based upon differential equations-that can be used to distinguish scale symmetry from scale freeness in empirical data. We establish face validity with an analysis of simulated data, in which we show how scale symmetry can be identified and how the associated conserved quantities can be estimated in neuronal time series.

Neural diffusivity and pre-emptive epileptic seizure intervention.

The propagation of epileptic seizure activity in the brain is a widespread pathophysiology that, in principle, should yield to intervention techniques guided by mathematical models of neuronal ensemble dynamics. During a seizure, neural activity will deviate from its current dynamical regime to one in which there are significant signal fluctuations. In silico treatments of neural activity are an important tool for the understanding of how the healthy brain can maintain stability, as well as of how pathology can lead to seizures. The hope is that, contained within the mathematical foundations of such treatments, there lie potential strategies for mitigating instabilities, e.g. via external stimulation. Here, we demonstrate that the dynamic causal modelling neuronal state equation generalises to a Fokker-Planck formalism if one extends the framework to model the ways in which activity propagates along the structural connections of neural systems. Using the Jacobian of this generalised state equation, we show that an initially unstable system can be rendered stable via a reduction in diffusivity-i.e., by lowering the rate at which neuronal fluctuations disperse to neighbouring regions. We show, for neural systems prone to epileptic seizures, that such a reduction in diffusivity can be achieved via external stimulation. Specifically, we show that this stimulation should be applied in such a way as to temporarily mirror the activity profile of a pathological region in its functionally connected areas. This counter-intuitive method is intended to be used pre-emptively-i.e., in order to mitigate the effects of the seizure, or ideally even prevent it from occurring in the first place. We offer proof of principle using simulations based on functional neuroimaging data collected from patients with idiopathic generalised epilepsy, in which we successfully suppress pathological activity in a distinct sub-network prior to seizure onset. Our hope is that this technique can form the basis for future real-time monitoring and intervention devices that are capable of treating epilepsy in a non-invasive manner.

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.

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.

Locus Coeruleus tracking of prediction errors optimises cognitive flexibility: An Active Inference model.

The locus coeruleus (LC) in the pons is the major source of noradrenaline (NA) in the brain. Two modes of LC firing have been associated with distinct cognitive states: changes in tonic rates of firing are correlated with global levels of arousal and behavioural flexibility, whilst phasic LC responses are evoked by salient stimuli. Here, we unify these two modes of firing by modelling the response of the LC as a correlate of a prediction error when inferring states for action planning under Active Inference (AI). We simulate a classic Go/No-go reward learning task and a three-arm 'explore/exploit' task and show that, if LC activity is considered to reflect the magnitude of high level 'state-action' prediction errors, then both tonic and phasic modes of firing are emergent features of belief updating. We also demonstrate that when contingencies change, AI agents can update their internal models more quickly by feeding back this state-action prediction error-reflected in LC firing and noradrenaline release-to optimise learning rate, enabling large adjustments over short timescales. We propose that such prediction errors are mediated by cortico-LC connections, whilst ascending input from LC to cortex modulates belief updating in anterior cingulate cortex (ACC). In short, we characterise the LC/ NA system within a general theory of brain function. In doing so, we show that contrasting, behaviour-dependent firing patterns are an emergent property of the LC that translates state-action prediction errors into an optimal balance between plasticity and stability.

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.

Altered Prefrontal Excitation/Inhibition Balance and Prefrontal Output: Markers of Aging in Human Memory Networks.

Memory impairments and heightened prefrontal cortical (PFC) activity are hallmarks of cognitive and neurobiological human aging. While structural integrity of PFC gray matter and interregional white matter tracts are thought to impact memory processing, the balance of neurotransmitters within the PFC itself is less well understood. We used fMRI to establish whole-brain networks involved in a memory encoding task and dynamic causal models (DCMs) for fMRI to determine the causal relationships between these areas. These data revealed enhanced connectivity from PFC to medial temporal cortex that negatively correlated with recall ability. To better understand the intrinsic activity within the PFC, DCM for EEG was employed after continuous theta burst transcranial magnetic stimulation (TMS) to the PFC to assess the effect on excitatory/inhibitory (E/I) synaptic ratios and behavior. These data revealed that the young cohort had a stable E/I ratio that was unaffected by the TMS intervention, while the aged cohort exhibited lower E/I ratios driven by a greater intrinsic inhibitory tone. TMS to the aged cohort resulted in decreased intrinsic inhibition and a decrement in memory performance. These results demonstrate increased top-down influence of PFC upon medial temporal lobe in healthy aging that is associated with decreased memory and may be due to unstable local inhibitory tone within the PFC.

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.

Computational modelling of movement-related beta-oscillatory dynamics in human motor cortex.

Oscillatory activity in the beta range, in human primary motor cortex (M1), shows interesting dynamics that are tied to behaviour and change systematically in disease. To investigate the pathophysiology underlying these changes, we must first understand how changes in beta activity are caused in healthy subjects. We therefore adapted a canonical (repeatable) microcircuit model used in dynamic causal modelling (DCM) previously used to model induced responses in visual cortex. We adapted this model to accommodate cytoarchitectural differences between visual and motor cortex. Using biologically plausible connections, we used Bayesian model selection to identify the best model of measured MEG data from 11 young healthy participants, performing a simple handgrip task. We found that the canonical M1 model had substantially more model evidence than the generic canonical microcircuit model when explaining measured MEG data. The canonical M1 model reproduced measured dynamics in humans at rest, in a manner consistent with equivalent studies performed in mice. Furthermore, the changes in excitability (self-inhibition) necessary to explain beta suppression during handgrip were consistent with the attenuation of sensory precision implied by predictive coding. These results establish the face validity of a model that can be used to explore the laminar interactions that underlie beta-oscillatory dynamics in humans in vivo. Our canonical M1 model may be useful for characterising the synaptic mechanisms that mediate pathophysiological beta dynamics associated with movement disorders, such as stroke or Parkinson's disease.

Precision and neuronal dynamics in the human posterior parietal cortex during evidence accumulation.

Primate studies show slow ramping activity in posterior parietal cortex (PPC) neurons during perceptual decision-making. These findings have inspired a rich theoretical literature to account for this activity. These accounts are largely unrelated to Bayesian theories of perception and predictive coding, a related formulation of perceptual inference in the cortical hierarchy. Here, we tested a key prediction of such hierarchical inference, namely that the estimated precision (reliability) of information ascending the cortical hierarchy plays a key role in determining both the speed of decision-making and the rate of increase of PPC activity. Using dynamic causal modelling of magnetoencephalographic (MEG) evoked responses, recorded during a simple perceptual decision-making task, we recover ramping-activity from an anatomically and functionally plausible network of regions, including early visual cortex, the middle temporal area (MT) and PPC. Precision, as reflected by the gain on pyramidal cell activity, was strongly correlated with both the speed of decision making and the slope of PPC ramping activity. Our findings indicate that the dynamics of neuronal activity in the human PPC during perceptual decision-making recapitulate those observed in the macaque, and in so doing we link observations from primate electrophysiology and human choice behaviour. Moreover, the synaptic gain control modulating these dynamics is consistent with predictive coding formulations of evidence accumulation.

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.

Free energy, precision and learning: the role of cholinergic neuromodulation.

Acetylcholine (ACh) is a neuromodulatory transmitter implicated in perception and learning under uncertainty. This study combined computational simulations and pharmaco-electroencephalography in humans, to test a formulation of perceptual inference based upon the free energy principle. This formulation suggests that ACh enhances the precision of bottom-up synaptic transmission in cortical hierarchies by optimizing the gain of supragranular pyramidal cells. Simulations of a mismatch negativity paradigm predicted a rapid trial-by-trial suppression of evoked sensory prediction error (PE) responses that is attenuated by cholinergic neuromodulation. We confirmed this prediction empirically with a placebo-controlled study of cholinesterase inhibition. Furthermore, using dynamic causal modeling, we found that drug-induced differences in PE responses could be explained by gain modulation in supragranular pyramidal cells in primary sensory cortex. This suggests that ACh adaptively enhances sensory precision by boosting bottom-up signaling when stimuli are predictable, enabling the brain to respond optimally under different levels of environmental uncertainty.

Broadband cortical desynchronization underlies the human psychedelic state.

Psychedelic drugs produce profound changes in consciousness, but the underlying neurobiological mechanisms for this remain unclear. Spontaneous and induced oscillatory activity was recorded in healthy human participants with magnetoencephalography after intravenous infusion of psilocybin--prodrug of the nonselective serotonin 2A receptor agonist and classic psychedelic psilocin. Psilocybin reduced spontaneous cortical oscillatory power from 1 to 50 Hz in posterior association cortices, and from 8 to 100 Hz in frontal association cortices. Large decreases in oscillatory power were seen in areas of the default-mode network. Independent component analysis was used to identify a number of resting-state networks, and activity in these was similarly decreased after psilocybin. Psilocybin had no effect on low-level visually induced and motor-induced gamma-band oscillations, suggesting that some basic elements of oscillatory brain activity are relatively preserved during the psychedelic experience. Dynamic causal modeling revealed that posterior cingulate cortex desynchronization can be explained by increased excitability of deep-layer pyramidal neurons, which are known to be rich in 5-HT2A receptors. These findings suggest that the subjective effects of psychedelics result from a desynchronization of ongoing oscillatory rhythms in the cortex, likely triggered by 5-HT2A receptor-mediated excitation of deep pyramidal cells.

Deterioration in cognitive control related mPFC function underlying development of treatment resistance in early psychosis.

One third of people with psychosis become antipsychotic treatment-resistant and the underlying mechanisms remain unclear. We investigated whether altered cognitive control function is a factor underlying development of treatment resistance. We studied 50 people with early psychosis at a baseline visit (mean < 2 years illness duration) and follow-up visit (1 year later), when 35 were categorized at treatment-responsive and 15 as treatment-resistant. Participants completed an emotion-yoked reward learning task that requires cognitive control whilst undergoing fMRI and MR spectroscopy to measure glutamate levels from Anterior Cingulate Cortex (ACC). Changes in cognitive control related activity (in prefrontal cortex and ACC) over time were compared between treatment-resistant and treatment-responsive groups and related to glutamate. Compared to treatment-responsive, treatment-resistant participants showed blunted activity in right amygdala (decision phase) and left pallidum (feedback phase) at baseline which increased over time and was accompanied by a decrease in medial Prefrontal Cortex (mPFC) activity (feedback phase) over time. Treatment-responsive participants showed a negative relationship between mPFC activity and glutamate levels at follow-up, no such relationship existed in treatment-resistant participants. Reduced activity in right amygdala and left pallidum at baseline was predictive of treatment resistance at follow-up (67% sensitivity, 94% specificity). The findings suggest that deterioration in mPFC function over time, a key cognitive control region needed to compensate for an initial dysfunction within a social-emotional network, is a factor underlying development of treatment resistance in early psychosis. An uncoupling between glutamate and cognitive control related mPFC function requires further investigation that may present a future target for interventions.

Basal ganglia-cortical interactions in Parkinsonian patients.

Parkinson's disease is a common and debilitating condition, caused by aberrant activity in a complex basal ganglia-thalamocortical circuit. Therapeutic advances rely on characterising interactions in this circuit. However, recording electrophysiological responses over the entire circuit is impractical. Dynamic causal modelling offers large-scale models of predictive value based on a limited or partial sampling of complex networks. Using dynamic causal modelling, we determined the network changes underlying the pathological excess of beta oscillations that characterise the Parkinsonian state. We modelled data from five patients undergoing surgery for deep brain stimulation of more than one target. We found that connections to and from the subthalamic nucleus were strengthened and promoted beta synchrony, in the untreated compared to the treated Parkinsonian state. Dynamic causal modelling was able to replicate the effects of lesioning this nucleus and may provide a new means of directing the search for therapeutic targets.

Network reconfiguration and working memory impairment in mesial temporal lobe epilepsy.

Mesial temporal lobe epilepsy (mTLE) is the most prevalent form of focal epilepsy, and hippocampal sclerosis (HS) is considered the most frequent associated pathological finding. Recent connectivity studies have shown that abnormalities, either structural or functional, are not confined to the affected hippocampus, but can be found in other connected structures within the same hemisphere, or even in the contralesional hemisphere. Despite the role of hippocampus in memory functions, most of these studies have explored network properties at resting state, and in some cases compared connectivity values with neuropsychological memory scores. Here, we measured magnetoencephalographic responses during verbal working memory (WM) encoding in left mTLE patients and controls, and compared their effective connectivity within a frontotemporal network using dynamic causal modelling. Bayesian model comparison indicated that the best model included bilateral, forward and backward connections, linking inferior temporal cortex (ITC), inferior frontal cortex (IFC), and the medial temporal lobe (MTL). Test for differences in effective connectivity revealed that patients exhibited decreased ipsilesional MTL-ITC backward connectivity, and increased bidirectional IFC-MTL connectivity in the contralesional hemisphere. Critically, a negative correlation was observed between these changes in patients, with decreases in ipsilesional coupling among temporal sources associated with increases contralesional frontotemporal interactions. Furthermore, contralesional frontotemporal interactions were inversely related to task performance and level of education. The results demonstrate that unilateral sclerosis induced local and remote changes in the dynamic organization of a distributed network supporting verbal WM. Crucially, pre-(peri) morbid factors (educational level) were reflected in both cognitive performance and (putative) compensatory changes in physiological coupling.