In Bayesian brain theories, hierarchically related prediction errors (PEs) play a central role for predicting sensory inputs and inferring their underlying causes, e.g., the probabilistic structure of the environment and its volatility. Notably, PEs at different hierarchical levels may be encoded by different neuromodulatory transmitters. Here, we tested this possibility in computational fMRI studies of audio-visual learning. Using a hierarchical Bayesian model, we found that low-level PEs about visual stimulus outcome were reflected by widespread activity in visual and supramodal areas but also in the midbrain. In contrast, high-level PEs about stimulus probabilities were encoded by the basal forebrain. These findings were replicated in two groups of healthy volunteers. While our fMRI measures do not reveal the exact neuron types activated in midbrain and basal forebrain, they suggest a dichotomy between neuromodulatory systems, linking dopamine to low-level PEs about stimulus outcome and acetylcholine to more abstract PEs about stimulus probabilities.

Models of effective connectivity characterize the influence that neuronal populations exert over each other. Additionally, some approaches, for example Dynamic Causal Modelling (DCM) and variants of Structural Equation Modelling, describe how effective connectivity is modulated by experimental manipulations. Mathematically, both are based on bilinear equations, where the bilinear term models the effect of experimental manipulations on neuronal interactions. The bilinear framework, however, precludes an important aspect of neuronal interactions that has been established with invasive electrophysiological recording studies; i.e., how the connection between two neuronal units is enabled or gated by activity in other units. These gating processes are critical for controlling the gain of neuronal populations and are mediated through interactions between synaptic inputs (e.g. by means of voltage-sensitive ion channels). They represent a key mechanism for various neurobiological processes, including top-down (e.g. attentional) modulation, learning and neuromodulation. This paper presents a nonlinear extension of DCM that models such processes (to second order) at the neuronal population level. In this way, the modulation of network interactions can be assigned to an explicit neuronal population. We present simulations and empirical results that demonstrate the validity and usefulness of this model. Analyses of synthetic data showed that nonlinear and bilinear mechanisms can be distinguished by our extended DCM. When applying the model to empirical fMRI data from a blocked attention to motion paradigm, we found that attentioninduced increases in V5 responses could be best explained as a gating of the V1→V5 connection by activity in posterior parietal cortex. Furthermore, we analysed fMRI data from an event-related binocular rivalry paradigm and found that interactions amongst percept-selective visual areas were modulated by activity in the middle frontal gyrus. In both practical examples, Bayesian model selection favoured the nonlinear models over corresponding bilinear ones.

This work demonstrates a fast, sensitive method of characterizing the dynamic performance of MR gradient systems. The accuracy of gradient time-courses is often compromised by field imperfections of various causes, including eddy currents and mechanical oscillations. Characterizing these perturbations is instrumental for corrections by pre-emphasis or post hoc signal processing. Herein, a gradient chain is treated as a linear time-invariant system, whose impulse response function is determined by measuring field responses to known gradient inputs. Triangular inputs are used to probe the system and response measurements are performed with a dynamic field camera consisting of NMR probes. In experiments on a whole-body MR system, it is shown that the proposed method yields impulse response functions of high temporal and spectral resolution. Besides basic properties such as bandwidth and delay, it also captures subtle features such as mechanically induced field oscillations. For validation, measured response functions were used to predict gradient field evolutions, which was achieved with an error below 0.2%. The field camera used records responses of various spatial orders simultaneously, rendering the method suitable also for studying cross-responses and dynamic shim systems. It thus holds promise for a range of applications, including pre-emphasis optimization, quality assurance, and image reconstruction.