Music in the brain.

Music is ubiquitous across human cultures - as a source of affective and pleasurable experience, moving us both physically and emotionally - and learning to play music shapes both brain structure and brain function. Music processing in the brain - namely, the perception of melody, harmony and rhythm - has traditionally been studied as an auditory phenomenon using passive listening paradigms. However, when listening to music, we actively generate predictions about what is likely to happen next. This enactive aspect has led to a more comprehensive understanding of music processing involving brain structures implicated in action, emotion and learning. Here we review the cognitive neuroscience literature of music perception. We show that music perception, action, emotion and learning all rest on the human brain's fundamental capacity for prediction - as formulated by the predictive coding of music model. This Review elucidates how this formulation of music perception and expertise in individuals can be extended to account for the dynamics and underlying brain mechanisms of collective music making. This in turn has important implications for human creativity as evinced by music improvisation. These recent advances shed new light on what makes music meaningful from a neuroscientific perspective.

Two distinct networks for encoding goals and forms of action: An effective connectivity study.

Goal-directed actions are characterized by two main features: the content (i.e., the action goal) and the form, called vitality forms (VF) (i.e., how actions are executed). It is well established that both the action content and the capacity to understand the content of another's action are mediated by a network formed by a set of parietal and frontal brain areas. In contrast, the neural bases of action forms (e.g., gentle or rude actions) have not been characterized. However, there are now studies showing that the observation and execution of actions endowed with VF activate, in addition to the parieto-frontal network, the dorso-central insula (DCI). In the present study, we established-using dynamic causal modeling (DCM)-the direction of information flow during observation and execution of actions endowed with gentle and rude VF in the human brain. Based on previous fMRI studies, the selected nodes for the DCM comprised the posterior superior temporal sulcus (pSTS), the inferior parietal lobule (IPL), the premotor cortex (PM), and the DCI. Bayesian model comparison showed that, during action observation, two streams arose from pSTS: one toward IPL, concerning the action goal, and one toward DCI, concerning the action vitality forms. During action execution, two streams arose from PM: one toward IPL, concerning the action goal and one toward DCI concerning action vitality forms. This last finding opens an interesting question concerning the possibility to elicit VF in two distinct ways: cognitively (from PM to DCI) and affectively (from DCI to PM).

Collective behavior from surprise minimization.

Collective motion is ubiquitous in nature; groups of animals, such as fish, birds, and ungulates appear to move as a whole, exhibiting a rich behavioral repertoire that ranges from directed movement to milling to disordered swarming. Typically, such macroscopic patterns arise from decentralized, local interactions among constituent components (e.g., individual fish in a school). Preeminent models of this process describe individuals as self-propelled particles, subject to self-generated motion and "social forces" such as short-range repulsion and long-range attraction or alignment. However, organisms are not particles; they are probabilistic decision-makers. Here, we introduce an approach to modeling collective behavior based on active inference. This cognitive framework casts behavior as the consequence of a single imperative: to minimize surprise. We demonstrate that many empirically observed collective phenomena, including cohesion, milling, and directed motion, emerge naturally when considering behavior as driven by active Bayesian inference-without explicitly building behavioral rules or goals into individual agents. Furthermore, we show that active inference can recover and generalize the classical notion of social forces as agents attempt to suppress prediction errors that conflict with their expectations. By exploring the parameter space of the belief-based model, we reveal nontrivial relationships between the individual beliefs and group properties like polarization and the tendency to visit different collective states. We also explore how individual beliefs about uncertainty determine collective decision-making accuracy. Finally, we show how agents can update their generative model over time, resulting in groups that are collectively more sensitive to external fluctuations and encode information more robustly.

From active affordance to active inference: vertical integration of cognition in the cerebral cortex through dual subcortical control systems.

In previous papers, we proposed that the dorsal attention system's top-down control is regulated by the dorsal division of the limbic system, providing a feedforward or impulsive form of control generating expectancies during active inference. In contrast, we proposed that the ventral attention system is regulated by the ventral limbic division, regulating feedback constraints and error-correction for active inference within the neocortical hierarchy. Here, we propose that these forms of cognitive control reflect vertical integration of subcortical arousal control systems that evolved for specific forms of behavior control. The feedforward impetus to action is regulated by phasic arousal, mediated by lemnothalamic projections from the reticular activating system of the lower brainstem, and then elaborated by the hippocampus and dorsal limbic division. In contrast, feedback constraint-based on environmental requirements-is regulated by the tonic activation furnished by collothalamic projections from the midbrain arousal control centers, and then sustained and elaborated by the amygdala, basal ganglia, and ventral limbic division. In an evolutionary-developmental analysis, understanding these differing forms of active affordance-for arousal and motor control within the subcortical vertebrate neuraxis-may help explain the evolution of active inference regulating the cognition of expectancy and error-correction within the mammalian 6-layered neocortex.

Bistable perception, precision and neuromodulation.

Bistable perception follows from observing a static, ambiguous, (visual) stimulus with two possible interpretations. Here, we present an active (Bayesian) inference account of bistable perception and posit that perceptual transitions between different interpretations (i.e. inferences) of the same stimulus ensue from specific eye movements that shift the focus to a different visual feature. Formally, these inferences are a consequence of precision control that determines how confident beliefs are and change the frequency with which one can perceive-and alternate between-two distinct percepts. We hypothesized that there are multiple, but distinct, ways in which precision modulation can interact to give rise to a similar frequency of bistable perception. We validated this using numerical simulations of the Necker cube paradigm and demonstrate the multiple routes that underwrite the frequency of perceptual alternation. Our results provide an (enactive) computational account of the intricate precision balance underwriting bistable perception. Importantly, these precision parameters can be considered the computational homologs of particular neurotransmitters-i.e. acetylcholine, noradrenaline, dopamine-that have been previously implicated in controlling bistable perception, providing a computational link between the neurochemistry and perception.

Neural Mechanisms and Psychology of Psychedelic Ego Dissolution.

Neuroimaging studies of psychedelics have advanced our understanding of hierarchical brain organization and the mechanisms underlying their subjective and therapeutic effects. The primary mechanism of action of classic psychedelics is binding to serotonergic 5-HT2A receptors. Agonist activity at these receptors leads to neuromodulatory changes in synaptic efficacy that can have a profound effect on hierarchical message-passing in the brain. Here, we review the cognitive and neuroimaging evidence for the effects of psychedelics: in particular, their influence on selfhood and subject-object boundaries-known as ego dissolution-surmised to underwrite their subjective and therapeutic effects. Agonism of 5-HT2A receptors, located at the apex of the cortical hierarchy, may have a particularly powerful effect on sentience and consciousness. These effects can endure well after the pharmacological half-life, suggesting that psychedelics may have effects on neural plasticity that may play a role in their therapeutic efficacy. Psychologically, this may be accompanied by a disarming of ego resistance that increases the repertoire of perceptual hypotheses and affords alternate pathways for thought and behavior, including those that undergird selfhood. We consider the interaction between serotonergic neuromodulation and sentience through the lens of hierarchical predictive coding, which speaks to the value of psychedelics in understanding how we make sense of the world and specific predictions about effective connectivity in cortical hierarchies that can be tested using functional neuroimaging. SIGNIFICANCE STATEMENT: Classic psychedelics bind to serotonergic 5-HT2A receptors. Their agonist activity at these receptors leads to neuromodulatory changes in synaptic efficacy, resulting in a profound effect on information processing in the brain. Here, we synthesize an abundance of brain imaging research with pharmacological and psychological interpretations informed by the framework of predictive coding. Moreover, predictive coding is suggested to offer more sophisticated interpretations of neuroimaging findings by bridging the role between the 5-HT2A receptors and large-scale brain networks.

Computational psychiatry: from synapses to sentience.

This review considers computational psychiatry from a particular viewpoint: namely, a commitment to explaining psychopathology in terms of pathophysiology. It rests on the notion of a generative model as underwriting (i) sentient processing in the brain, and (ii) the scientific process in psychiatry. The story starts with a view of the brain-from cognitive and computational neuroscience-as an organ of inference and prediction. This offers a formal description of neuronal message passing, distributed processing and belief propagation in neuronal networks; and how certain kinds of dysconnection lead to aberrant belief updating and false inference. The dysconnections in question can be read as a pernicious synaptopathy that fits comfortably with formal notions of how we-or our brains-encode uncertainty or its complement, precision. It then considers how the ensuing process theories are tested empirically, with an emphasis on the computational modelling of neuronal circuits and synaptic gain control that mediates attentional set, active inference, learning and planning. The opportunities afforded by this sort of modelling are considered in light of in silico experiments; namely, computational neuropsychology, computational phenotyping and the promises of a computational nosology for psychiatry. The resulting survey of computational approaches is not scholarly or exhaustive. Rather, its aim is to review a theoretical narrative that is emerging across subdisciplines within psychiatry and empirical scales of investigation. These range from epilepsy research to neurodegenerative disorders; from post-traumatic stress disorder to the management of chronic pain, from schizophrenia to functional medical symptoms.

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.

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.

Longitudinal motor system changes from acute to chronic spinal cord injury.

BACKGROUND AND PURPOSE: In acute spinal cord injury (SCI), magnetic resonance imaging (MRI) reveals tissue bridges and neurodegeneration for 2 years. This 5-year study aims to track initial lesion changes, subsequent neurodegeneration, and their impact on recovery. METHODS: This prospective longitudinal study enrolled acute SCI patients and healthy controls who were assessed clinically-and by MRI-regularly from 3 days postinjury up to 60 months. We employed histologically cross-validated quantitative MRI sequences sensitive to volume, myelin, and iron changes, thereby reflecting indirectly processes of neurodegeneration and neuroinflammation. General linear models tracked lesion and remote changes in volume, myelin- and iron-sensitive magnetic resonance indices over 5 years. Associations between lesion, degeneration, and recovery (using the Spinal Cord Independence Measure [SCIM] questionnaire and the International Standards for Neurological Classification of Spinal Cord Injury total motor score) were assessed. RESULTS: Patients' motor scores improved by an average of 12.86 (95% confidence interval [CI] = 6.70-19.00) points, and SCIM by 26.08 (95% CI = 17.00-35.20) points. Within 3-28 days post-SCI, lesion size decreased by more than two-thirds (3 days: 302.52 ± 185.80 mm2 , 28 days: 76.77 ± 88.62 mm2 ), revealing tissue bridges. Cervical cord and corticospinal tract volumes transiently increased in SCI patients by 5% and 3%, respectively, accompanied by cervical myelin decreases and iron increases. Over time, progressive atrophy was observed in both regions, which was linked to early lesion dynamics. Tissue bridges, reduced swelling, and myelin content decreases were predictive of long-term motor score recovery and improved SCIM score. CONCLUSIONS: Studying acute changes and their impact on longer follow-up provides insights into SCI trajectory, highlighting the importance of acute intervention while indicating the potential to influence outcomes in the later stages.

Predictive coding and stochastic resonance as fundamental principles of auditory phantom perception.

Mechanistic insight is achieved only when experiments are employed to test formal or computational models. Furthermore, in analogy to lesion studies, phantom perception may serve as a vehicle to understand the fundamental processing principles underlying healthy auditory perception. With a special focus on tinnitus-as the prime example of auditory phantom perception-we review recent work at the intersection of artificial intelligence, psychology and neuroscience. In particular, we discuss why everyone with tinnitus suffers from (at least hidden) hearing loss, but not everyone with hearing loss suffers from tinnitus. We argue that intrinsic neural noise is generated and amplified along the auditory pathway as a compensatory mechanism to restore normal hearing based on adaptive stochastic resonance. The neural noise increase can then be misinterpreted as auditory input and perceived as tinnitus. This mechanism can be formalized in the Bayesian brain framework, where the percept (posterior) assimilates a prior prediction (brain's expectations) and likelihood (bottom-up neural signal). A higher mean and lower variance (i.e. enhanced precision) of the likelihood shifts the posterior, evincing a misinterpretation of sensory evidence, which may be further confounded by plastic changes in the brain that underwrite prior predictions. Hence, two fundamental processing principles provide the most explanatory power for the emergence of auditory phantom perceptions: predictive coding as a top-down and adaptive stochastic resonance as a complementary bottom-up mechanism. We conclude that both principles also play a crucial role in healthy auditory perception. Finally, in the context of neuroscience-inspired artificial intelligence, both processing principles may serve to improve contemporary machine learning techniques.

Neurophysiological consequences of synapse loss in progressive supranuclear palsy.

Synaptic loss occurs early in many neurodegenerative diseases and contributes to cognitive impairment even in the absence of gross atrophy. Currently, for human disease there are few formal models to explain how cortical networks underlying cognition are affected by synaptic loss. We advocate that biophysical models of neurophysiology offer both a bridge from preclinical to clinical models of pathology and quantitative assays for experimental medicine. Such biophysical models can also disclose hidden neuronal dynamics generating neurophysiological observations such as EEG and magnetoencephalography. Here, we augment a biophysically informed mesoscale model of human cortical function by inclusion of synaptic density estimates as captured by 11C-UCB-J PET, and provide insights into how regional synapse loss affects neurophysiology. We use the primary tauopathy of progressive supranuclear palsy (Richardson's syndrome) as an exemplar condition, with high clinicopathological correlations. Progressive supranuclear palsy causes a marked change in cortical neurophysiology in the presence of mild cortical atrophy and is associated with a decline in cognitive functions associated with the frontal lobe. Using parametric empirical Bayesian inversion of a conductance-based canonical microcircuit model of magnetoencephalography data, we show that the inclusion of regional synaptic density-as a subject-specific prior on laminar-specific neuronal populations-markedly increases model evidence. Specifically, model comparison suggests that a reduction in synaptic density in inferior frontal cortex affects superficial and granular layer glutamatergic excitation. This predicted individual differences in behaviour, demonstrating the link between synaptic loss, neurophysiology and cognitive deficits. The method we demonstrate is not restricted to progressive supranuclear palsy or the effects of synaptic loss: such pathology-enriched dynamic causal models can be used to assess the mechanisms of other neurological disorders, with diverse non-invasive measures of pathology, and is suitable to test the effects of experimental pharmacology.

Altered effective connectivity in sensorimotor cortices is a signature of severity and clinical course in depression.

Functional neuroimaging research on depression has traditionally targeted neural networks associated with the psychological aspects of depression. In this study, instead, we focus on alterations of sensorimotor function in depression. We used resting-state functional MRI data and dynamic causal modeling (DCM) to assess the hypothesis that depression is associated with aberrant effective connectivity within and between key regions in the sensorimotor hierarchy. Using hierarchical modeling of between-subject effects in DCM with parametric empirical Bayes we first established the architecture of effective connectivity in sensorimotor cortices. We found that in (interoceptive and exteroceptive) sensory cortices across participants, the backward connections are predominantly inhibitory, whereas the forward connections are mainly excitatory in nature. In motor cortices these parities were reversed. With increasing depression severity, these patterns are depreciated in exteroceptive and motor cortices and augmented in the interoceptive cortex, an observation that speaks to depressive symptomatology. We established the robustness of these results in a leave-one-out cross-validation analysis and by reproducing the main results in a follow-up dataset. Interestingly, with (nonpharmacological) treatment, depression-associated changes in backward and forward effective connectivity partially reverted to group mean levels. Overall, altered effective connectivity in sensorimotor cortices emerges as a promising and quantifiable candidate marker of depression severity and treatment response.

Shared Protentions in Multi-Agent Active Inference.

In this paper, we unite concepts from Husserlian phenomenology, the active inference framework in theoretical biology, and category theory in mathematics to develop a comprehensive framework for understanding social action premised on shared goals. We begin with an overview of Husserlian phenomenology, focusing on aspects of inner time-consciousness, namely, retention, primal impression, and protention. We then review active inference as a formal approach to modeling agent behavior based on variational (approximate Bayesian) inference. Expanding upon Husserl's model of time consciousness, we consider collective goal-directed behavior, emphasizing shared protentions among agents and their connection to the shared generative models of active inference. This integrated framework aims to formalize shared goals in terms of shared protentions, and thereby shed light on the emergence of group intentionality. Building on this foundation, we incorporate mathematical tools from category theory, in particular, sheaf and topos theory, to furnish a mathematical image of individual and group interactions within a stochastic environment. Specifically, we employ morphisms between polynomial representations of individual agent models, allowing predictions not only of their own behaviors but also those of other agents and environmental responses. Sheaf and topos theory facilitates the construction of coherent agent worldviews and provides a way of representing consensus or shared understanding. We explore the emergence of shared protentions, bridging the phenomenology of temporal structure, multi-agent active inference systems, and category theory. Shared protentions are highlighted as pivotal for coordination and achieving common objectives. We conclude by acknowledging the intricacies stemming from stochastic systems and uncertainties in realizing shared goals.

Natural language syntax complies with the free-energy principle.

Natural language syntax yields an unbounded array of hierarchically structured expressions. We claim that these are used in the service of active inference in accord with the free-energy principle (FEP). While conceptual advances alongside modelling and simulation work have attempted to connect speech segmentation and linguistic communication with the FEP, we extend this program to the underlying computations responsible for generating syntactic objects. We argue that recently proposed principles of economy in language design-such as "minimal search" criteria from theoretical syntax-adhere to the FEP. This affords a greater degree of explanatory power to the FEP-with respect to higher language functions-and offers linguistics a grounding in first principles with respect to computability. While we mostly focus on building new principled conceptual relations between syntax and the FEP, we also show through a sample of preliminary examples how both tree-geometric depth and a Kolmogorov complexity estimate (recruiting a Lempel-Ziv compression algorithm) can be used to accurately predict legal operations on syntactic workspaces, directly in line with formulations of variational free energy minimization. This is used to motivate a general principle of language design that we term Turing-Chomsky Compression (TCC). We use TCC to align concerns of linguists with the normative account of self-organization furnished by the FEP, by marshalling evidence from theoretical linguistics and psycholinguistics to ground core principles of efficient syntactic computation within active inference.

A primer on Variational Laplace (VL).

This article details a scheme for approximate Bayesian inference, which has underpinned thousands of neuroimaging studies since its introduction 15 years ago. Variational Laplace (VL) provides a generic approach to fitting linear or non-linear models, which may be static or dynamic, returning a posterior probability density over the model parameters and an approximation of log model evidence, which enables Bayesian model comparison. VL applies variational Bayesian inference in conjunction with quadratic or Laplace approximations of the evidence lower bound (free energy). Importantly, update equations do not need to be derived for each model under consideration, providing a general method for fitting a broad class of models. This primer is intended for experimenters and modellers who may wish to fit models to data using variational Bayesian methods, without assuming previous experience of variational Bayes or machine learning. Accompanying code demonstrates how to fit different kinds of model using the reference implementation of the VL scheme in the open-source Statistical Parametric Mapping (SPM) software package. In addition, we provide a standalone software function that does not require SPM, in order to ease translation to other fields, together with detailed pseudocode. Finally, the supplementary materials provide worked derivations of the key equations.

Attentional effects on local V1 microcircuits explain selective V1-V4 communication.

Selective attention implements preferential routing of attended stimuli, likely through increasing the influence of the respective synaptic inputs on higher-area neurons. As the inputs of competing stimuli converge onto postsynaptic neurons, presynaptic circuits might offer the best target for attentional top-down influences. If those influences enabled presynaptic circuits to selectively entrain postsynaptic neurons, this might explain selective routing. Indeed, when two visual stimuli induce two gamma rhythms in V1, only the gamma induced by the attended stimulus entrains gamma in V4. Here, we modelled induced responses with a Dynamic Causal Model for Cross-Spectral Densities and found that selective entrainment can be explained by attentional modulation of intrinsic V1 connections. Specifically, local inhibition was decreased in the granular input layer and increased in the supragranular output layer of the V1 circuit that processed the attended stimulus. Thus, presynaptic attentional influences and ensuing entrainment were sufficient to mediate selective routing.

Coherent, time-shifted patterns of microstructural plasticity during motor-skill learning.

Motor skill learning relies on neural plasticity in the motor and limbic systems. However, the spatial and temporal characteristics of these changes-and their microstructural underpinnings-remain unclear. Eighteen healthy males received 1 h of training in a computer-based motion game, 4 times a week, for 4 consecutive weeks, while 14 untrained participants underwent scanning only. Performance improvements were observed in all trained participants. Serial myelin- and iron-sensitive multiparametric mapping at 3T during this period of intensive motor skill acquisition revealed temporally and spatially distributed, performance-related microstructural changes in the grey and white matter across a corticospinal-cerebellar-hippocampal circuit. Analysis of the trajectory of these transient changes suggested time-shifted cascades of plasticity from the dominant sensorimotor system to the contralateral hippocampus. In the cranial corticospinal tracts, changes in myelin-sensitive metrics during training in the posterior limb of the internal capsule were of greater magnitude in those who trained their upper limbs vs. lower limb trainees. Motor skill learning is associated with waves of grey and white matter plasticity, across a broad sensorimotor network.

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

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

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.