Search

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.

We propose a computational neurology of movement based on the convergence of theoretical neurobiology and clinical neurology. A significant development in the former is the idea that we can frame brain function as a process of (active) inference, in which the nervous system makes predictions about its sensory data. These predictions depend upon an implicit predictive (generative) model used by the brain. This means neural dynamics can be framed as generating actions to ensure sensations are consistent with these predictions—and adjusting predictions when they are not. We illustrate the significance of this formulation for clinical neurology by simulating a clinical examination of the motor system using an upper limb coordination task. Specifically, we show how tendon reflexes emerge naturally under the right kind of generative model. Through simulated perturbations, pertaining to prior probabilities of this model’s variables, we illustrate the emergence of hyperreflexia and pendular reflexes, reminiscent of neurological lesions in the corticospinal tract and cerebellum. We then turn to the computational lesions causing hypokinesia and deficits of coordination. This in silico lesion-deficit analysis provides an opportunity to revisit classic neurological dichotomies (e.g. pyramidal versus extrapyramidal systems) from the perspective of modern approaches to theoretical neurobiology—and our understanding of the neurocomputational architecture of movement control based on first principles., Parr et al. report a series of numerical simulations—based on Bayesian message passing—of healthy and pathological motor control. They perform an in silico lesion-deficit analysis and demonstrate how the sorts of syndromes observed during a clinical neurological examination emerge from pathological computation.

Research in clinical neuroscience is founded on the idea that a better understanding of brain (dys)function will improve our ability to diagnose and treat neurological and psychiatric disorders. In recent years, neuroscience has converged on the notion that the brain is a 'prediction machine,' in that it actively predicts the sensory input that it will receive if one or another course of action is chosen. These predictions are used to select actions that will (most often, and in the long run) maintain the body within the narrow range of physiological states consistent with survival. This insight has given rise to an area of clinical computational neuroscience research that focuses on characterizing neural circuit architectures that can accomplish these predictive functions, and on how the associated processes may break down or become aberrant within clinical conditions. Here, we provide a brief review of examples of recent work on the application of predictive processing models of brain function to study clinical (psychiatric) disorders, with the aim of highlighting current directions and their potential clinical utility. We offer examples of recent conceptual models, formal mathematical models, and applications of such models in empirical research in clinical populations, with a focus on making this material accessible to clinicians without expertise in computational neuroscience. In doing so, we aim to highlight the potential insights and opportunities that understanding the brain as a prediction machine may offer to clinical research and practice.

Changes in neurovascular coupling are associated with both Alzheimer's disease and vascular dementia in later life, but this may be confounded by cerebrovascular risk. We hypothesized that hemodynamic latency would be associated with reduced cognitive functioning across the lifespan, holding constant demographic and cerebrovascular risk. In 387 adults aged 18-85 (mean = 48.82), dynamic causal modeling was used to estimate the hemodynamic response function in the left and right V1 and V3-ventral regions of the visual cortex in response to a simple checkerboard block design stimulus with minimal cognitive demands. The hemodynamic latency (transit time) in the visual cortex was used to predict general cognitive ability (Full-Scale IQ), controlling for demographic variables (age, race, education, socioeconomic status) and cerebrovascular risk factors (hypertension, alcohol use, smoking, high cholesterol, BMI, type 2 diabetes, cardiac disorders). Increased hemodynamic latency in the visual cortex predicted reduced cognitive function (p

The moment‐to‐moment focus of our mind's eye results from a complex interplay of voluntary and involuntary influences on attention. Previous neuroimaging studies suggest that the brain networks of voluntary versus involuntary attention can be segregated into a frontal‐versus‐parietal or a dorsal‐versus‐ventral partition—although recent work suggests that the dorsal network may be involved in both bottom‐up and top‐down attention. Research with nonhuman primates has provided evidence that a key distinction between top‐down and bottom‐up attention may be the direction of connectivity between frontal and parietal areas. Whereas typical fMRI connectivity analyses cannot disambiguate the direction of connections, dynamic causal modeling (DCM) can model directionality. Using DCM, we provide new evidence that directed connections within the dorsal attention network are differentially modulated for voluntary versus involuntary attention. These results suggest that the intraparietal sulcus exerts a baseline inhibitory effect on the frontal eye fields that is strengthened during exogenous orienting and attenuated during endogenous orienting. Furthermore, the attenuation from endogenous attention occurs even with salient peripheral cues when those cues are known to be counter predictive. Thus, directed connectivity between frontal and parietal regions of the dorsal attention network is highly influenced by the type of attention that is engaged.

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., NeuroImage, 237, ISSN:1053-8119, ISSN:1095-9572

This technical note introduces adiabatic dynamic causal modelling, a method for inferring slow changes in biophysical parameters that control fluctuations of fast neuronal states. The application domain we have in mind is inferring slow changes in variables (e.g., extracellular ion concentrations or synaptic efficacy) that underlie phase transitions in brain activity (e.g., paroxysmal seizure activity). The scheme is efficient and yet retains a biophysical interpretation, in virtue of being based on established neural mass models that are equipped with a slow dynamic on the parameters (such as synaptic rate constants or effective connectivity). In brief, we use an adiabatic approximation to summarise fast fluctuations in hidden neuronal states (and their expression in sensors) in terms of their second order statistics; namely, their complex cross spectra. This allows one to specify and compare models of slowly changing parameters (using Bayesian model reduction) that generate a sequence of empirical cross spectra of electrophysiological recordings. Crucially, we use the slow fluctuations in the spectral power of neuronal activity as empirical priors on changes in synaptic parameters. This introduces a circular causality, in which synaptic parameters underwrite fast neuronal activity that, in turn, induces activity-dependent plasticity in synaptic parameters. In this foundational paper, we describe the underlying model, establish its face validity using simulations and provide an illustrative application to a chemoconvulsant animal model of seizure activity.

This paper addresses perceptual synthesis by comparing responses evoked by visual stimuli before and after they are recognized, depending on prior exposure. Using magnetoencephalography, we analyzed distributed patterns of neuronal activity – evoked by Mooney figures – before and after they were recognized as meaningful objects. Recognition induced changes were first seen at 100–120 ​ms, for both faces and tools. These early effects – in right inferior and middle occipital regions – were characterized by an increase in power in the absence of any changes in spatial patterns of activity. Within a later 210–230 ​ms window, a quite different type of recognition effect appeared. Regions of the brain’s value system (insula, entorhinal cortex and cingulate of the right hemisphere for faces and right orbitofrontal cortex for tools) evinced a reorganization of their neuronal activity without an overall power increase in the region. Finally, we found that during the perception of disambiguated face stimuli, a face-specific response in the right fusiform gyrus emerged at 240–290 ​ms, with a much greater latency than the well-known N170m component, and, crucially, followed the recognition effect in the value system regions. These results can clarify one of the most intriguing issues of perceptual synthesis, namely, how a limited set of high-level predictions, which is required to reduce the uncertainty when resolving the ill-posed inverse problem of perception, can be available before category-specific processing in visual cortex. We suggest that a subset of local spatial features serves as partial cues for a fast re-activation of object-specific appraisal by the value system. The ensuing top-down feedback from value system to visual cortex, in particular, the fusiform gyrus enables high levels of processing to form category-specific predictions. This descending influence of the value system was more prominent for faces than for tools, the fact that reflects different dependence of these categories on value-related information., Highlights • Recognition effects following Mooney figure disambiguation were decoded from MEG data. • Neural changes occur in 3 stages that challenge the hierarchical view of perception. • Category-specific activity in the brain’s value system precedes that in visual cortex. • Face-specific response in the right fusiform gyrus is delayed for disambiguated faces. • The top-down influence of the value system is more prominent for faces than for tools.

Meta-awareness refers to the capacity to explicitly notice the current content of consciousness and has been identified as a key component for the successful control of cognitive states, such as the deliberate direction of attention. This paper proposes a formal model of meta-awareness and attentional control using hierarchical active inference. To do so, we cast mental action as policy selection over higher-level cognitive states and add a further hierarchical level to model meta-awareness states that modulate the expected confidence (precision) in the mapping between observations and hidden cognitive states. We simulate the example of mind-wandering and its regulation during a task involving sustained selective attention on a perceptual object. This provides a computational case study for an inferential architecture that is apt to enable the emergence of these central components of human phenomenology, namely, the ability to access and control cognitive states. We propose that this approach can be generalized to other cognitive states, and hence, this paper provides the first steps towards the development of a computational phenomenology of mental action and more broadly of our ability to monitor and control our own cognitive states. Future steps of this work will focus on fitting the model with qualitative, behavioural, and neural data.

Previous studies suggest that the practice of long-term (months to years) mindfulness meditation induces structural plasticity in gray matter. However, it remains unknown whether short-term (

This technical note presents a dynamic causal modelling (DCM) procedure for evaluating different models of neurovascular coupling in the human brain – using combined electromagnetic (M/EEG) and functional magnetic resonance imaging (fMRI) data. This procedure compares the evidence for biologically informed models of neurovascular coupling using Bayesian model comparison. First, fMRI data are used to localise regionally specific neuronal responses. The coordinates of these responses are then used as the location priors in a DCM of electrophysiological responses elicited by the same paradigm. The ensuing estimates of model parameters are then used to generate neuronal drive functions, which model pre- or post-synaptic activity for each experimental condition. These functions form the input to a model of neurovascular coupling, whose parameters are estimated from the fMRI data. Crucially, this enables one to evaluate different models of neurovascular coupling, using Bayesian model comparison – asking, for example, whether instantaneous or delayed, pre- or post-synaptic signals mediate haemodynamic responses. We provide an illustrative application of the procedure using a single-subject auditory fMRI and MEG dataset. The code and exemplar data accompanying this technical note are available through the statistical parametric mapping (SPM) software., Highlights • A method is introduced for Bayesian fusion of M/EEG and BOLD fMRI data using dynamic causal modelling. • The key novel contribution is the introduction of neural drive functions, which link models of the two modalities. • Bayesian model comparison is used to select plausible hypotheses about the function of neurovascular coupling.

The influence of global BOLD fluctuations on resting state functional connectivity in fMRI data remains a topic of debate, with little consensus. In this study, we assessed the effects of global signal regression (GSR) on effective connectivity within and between resting state networks (RSNs) - as estimated with dynamic causal modelling (DCM) for resting state fMRI (rsfMRI). DCM incorporates a forward (generative) model that quantifies the contribution of different types of noise (including global measurement noise), effective connectivity, and (neuro)vascular processes to functional connectivity measurements. DCM analyses were applied to two different designs; namely, longitudinal and cross-sectional designs. In the modelling of longitudinal designs, we considered four extensive longitudinal resting state fMRI datasets with a total number of 20 subjects. In the analysis of cross-sectional designs, we used rsfMRI data from 361 subjects from the Human Connectome Project. We hypothesized that (1) GSR would have no discernible impact on effective connectivity estimated with DCM, and (2) GSR would be reflected in the parameters representing global measurement noise. Additionally, we performed comparative analyses of information gain with and without GSR. Our results showed negligible to small effects of GSR on effective connectivity within small (separately estimated) RSNs. However, although the effect sizes were small, there was substantial to conclusive evidence for an effect of GSR on connectivity parameters. For between-network connectivity, we found two important effects: the effect of GSR on between-network effective connectivity (averaged over all connections) was negligible to small, while the effect of GSR on individual connections was non-negligible. In the cross-sectional (but not in the longitudinal) data, some connections showed substantial to conclusive evidence for an effect of GSR. Contrary to our expectations, we found either no effect (in the longitudinal designs) or a non-specific (cross-sectional design) effect of GSR on parameters characterising (global) measurement noise. Data without GSR were found to be more informative than data with GSR; however, in small resting state networks the precision of posterior estimates was greater after GSR. In conclusion, GSR is a minor concern in DCM studies; however, quantitative interpretation of between-network connections (as opposed to average between-network connectivity) and noise parameters should be treated with some caution. The Kullback-Leibler divergence of the posterior from the prior (i.e., information gain) - together with the precision of posterior estimates - might offer a useful measure to assess the appropriateness of GSR in resting state fMRI., Highlights • We investigated the effect of GSR on parameters estimated with DCM. • We found negligible to small effects of GSR in small-scale resting state networks. • Small effects of GSR on the net influence between networks were observed. • Results showed non-negligible effects of GSR on individual between-network connections. • The effect of GSR on noise parameters was either absent or unspecific, depending on design.

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

Rapid imaging techniques are increasingly used in functional MRI studies because they allow a greater number of samples to be acquired per unit time, thereby increasing statistical power. However, temporal correlations limit the increase in functional sensitivity and must be accurately accounted for to control the false‐positive rate. A common approach to accounting for temporal correlations is to whiten the data prior to estimating fMRI model parameters. Models of white noise plus a first‐order autoregressive process have proven sufficient for conventional imaging studies, but more elaborate models are required for rapidly sampled data. Here we show that when the “FAST” model implemented in SPM is used with a well‐controlled number of parameters, it can successfully prewhiten 80% of grey matter voxels even with volume repetition times as short as 0.35 s. We further show that the temporal signal‐to‐noise ratio (tSNR), which has conventionally been used to assess the relative functional sensitivity of competing imaging approaches, can be augmented to account for the temporal correlations in the time series. This amounts to computing the t‐score testing for the mean signal. We show in a visual perception task that unlike the tSNR weighted by the number of samples, the t‐score measure is directly related to the t‐score testing for activation when the temporal correlations are correctly modeled. This score affords a more accurate means of evaluating the functional sensitivity of different data acquisition options.

Although the relationship between resting-state functional connectivity and task-related activity has been addressed, the relationship between task and resting-state directed or effective connectivity – and its behavioral concomitants – remains elusive. We evaluated effective connectivity under an N-back working memory task in 24 participants using stochastic dynamic causal modelling (DCM) of 7 T fMRI data. We repeated the analysis using resting-state data, from the same subjects, to model connectivity among the same brain regions engaged by the N-back task. This allowed us to: (i) examine the relationship between intrinsic (task-independent) effective connectivity during resting (Arest) and task states (Atask), (ii) cluster phenotypes of task-related changes in effective connectivity (Btask) across participants, (iii) identify edges (Btask) showing high inter-individual effective connectivity differences and (iv) associate reaction times with the similarity between Btask and Arest in these edges. We found a strong correlation between Arest and Atask over subjects but a marked difference between Btask and Arest. We further observed a strong clustering of individuals in terms of Btask, which was not apparent in Arest. The task-related effective connectivity Btask varied highly in the edges from the parietal to the frontal lobes across individuals, so the three groups were clustered mainly by the effective connectivity within these networks. The similarity between Btask and Arest at the edges from the parietal to the frontal lobes was positively correlated with 2-back reaction times. This result implies that a greater change in context-sensitive coupling – from resting-state connectivity – is associated with faster reaction times. In summary, task-dependent connectivity endows resting-state connectivity with a context sensitivity, which predicts the speed of information processing during the N-back task.

Schizophrenia is a disorder characterized by functional dysconnectivity among distributed brain regions. However, it is unclear how causal influences among large-scale brain networks are disrupted in schizophrenia. In this study, we used dynamic causal modeling (DCM) to assess the hypothesis that there is aberrant directed (effective) connectivity within and between three key large-scale brain networks (the dorsal attention network, the salience network and the default mode network) in schizophrenia during a working memory task. Functional MRI data during an n-back task from 40 patients with schizophrenia and 62 healthy controls were analyzed. Using hierarchical modeling of between-subject effects in DCM with Parametric Empirical Bayes, we found that intrinsic (within-region) and extrinsic (between-region) effective connectivity involving prefrontal regions were abnormal in schizophrenia. Specifically, in patients (i) inhibitory self-connections in prefrontal regions of the dorsal attention network were decreased across task conditions; (ii) extrinsic connectivity between regions of the default mode network was increased; specifically, from posterior cingulate cortex to the medial prefrontal cortex; (iii) between-network extrinsic connections involving the prefrontal cortex were altered; (iv) connections within networks and between networks were correlated with the severity of clinical symptoms and impaired cognition beyond working memory. In short, this study revealed the predominance of reduced synaptic efficacy of prefrontal efferents and afferents in the pathophysiology of schizophrenia., Highlights • A first use of hierarchical modeling of effective connectivity to characterize large-scale networks in schizophrenia. • Intrinsic and extrinsic effective connectivity involving prefrontal regions were abnormal in schizophrenia. • Diagnostic connections could predict the severity of clinical symptoms and cognition in schizophrenia.

Objective Posttraumatic stress disorder (PTSD) is associated with decreased top-down emotion modulation from medial prefrontal cortex (mPFC) regions, a pathophysiology accompanied by hyperarousal and hyperactivation of the amygdala. By contrast, PTSD patients with the dissociative subtype (PTSD + DS) often exhibit increased mPFC top-down modulation and decreased amygdala activation associated with emotional detachment and hypoarousal. Crucially, PTSD and PTSD + DS display distinct functional connectivity within the PFC, amygdala complexes, and the periaqueductal gray (PAG), a region related to defensive responses/emotional coping. However, differences in directed connectivity between these regions have not been established in PTSD, PTSD + DS, or controls. Methods To examine directed (effective) connectivity among these nodes, as well as group differences, we conducted resting-state stochastic dynamic causal modeling (sDCM) pairwise analyses of coupling between the ventromedial (vm)PFC, the bilateral basolateral and centromedial (CMA) amygdala complexes, and the PAG, in 155 participants (PTSD [n = 62]; PTSD + DS [n = 41]; age-matched healthy trauma-unexposed controls [n = 52]). Results PTSD was characterized by a pattern of predominant bottom-up connectivity from the amygdala to the vmPFC and from the PAG to the vmPFC and amygdala. Conversely, PTSD + DS exhibited predominant top-down connectivity between all node pairs (from the vmPFC to the amygdala and PAG, and from the amygdala to the PAG). Interestingly, the PTSD + DS group displayed the strongest intrinsic inhibitory connections within the vmPFC. Conclusions These results suggest the contrasting symptom profiles of PTSD and its dissociative subtype (hyper- vs. hypo-emotionality, respectively) may be driven by complementary changes in directed connectivity corresponding to bottom-up defensive fear processing versus enhanced top-down regulation. Hum Brain Mapp 38:5551-5561, 2017. © 2017 Wiley Periodicals, Inc.

The development of large-scale network models that infer the effective (directed) connectivity among neuronal populations from neuroimaging data represents a key challenge for computational neuroscience. Dynamic causal models (DCMs) of neuroimaging and electrophysiological data are frequently used for inferring effective connectivity but are presently restricted to small graphs (typically up to 10 regions) in order to keep model inversion computationally feasible. Here, we present a novel variant of DCM for functional magnetic resonance imaging (fMRI) data that is suited to assess effective connectivity in large (whole-brain) networks. The approach rests on translating a linear DCM into the frequency domain and reformulating it as a special case of Bayesian linear regression. This paper derives regression DCM (rDCM) in detail and presents a variational Bayesian inversion method that enables extremely fast inference and accelerates model inversion by several orders of magnitude compared to classical DCM. Using both simulated and empirical data, we demonstrate the face validity of rDCM under different settings of signal-to-noise ratio (SNR) and repetition time (TR) of fMRI data. In particular, we assess the potential utility of rDCM as a tool for whole-brain connectomics by challenging it to infer effective connection strengths in a simulated whole-brain network comprising 66 regions and 300 free parameters. Our results indicate that rDCM represents a computationally highly efficient approach with promising potential for inferring whole-brain connectivity from individual fMRI data., NeuroImage, 155, ISSN:1053-8119, ISSN:1095-9572

Effective connectivity is commonly assessed using blood oxygenation level-dependent (BOLD) signals. In (Havlicek et al., 2015), we presented a novel, physiologically informed dynamic causal model (P-DCM) that extends current generative models. We demonstrated the improvements afforded by P-DCM in terms of the ability to model commonly observed neuronal and vascular transients in single regions. Here, we assess the ability of the novel and previous DCM variants to estimate effective connectivity among a network of five ROIs driven by a visuo-motor task. We demonstrate that connectivity estimates depend sensitively on the DCM used, due to differences in the modeling of hemodynamic response transients; such as the post-stimulus undershoot or adaptation during stimulation. In addition, using a novel DCM for arterial spin labeling (ASL) fMRI that measures BOLD and CBF signals simultaneously, we confirmed our findings (by using the BOLD data alone and in conjunction with CBF). We show that P-DCM provides better estimates of effective connectivity, regardless of whether it is applied to BOLD data alone or to ASL time-series, and that all new aspects of P-DCM (i.e. neuronal, neurovascular, hemodynamic components) constitute an improvement compared to those in the previous DCM variants. In summary, (i) accurate modeling of fMRI response transients is crucial to obtain valid effective connectivity estimates and (ii) any additional hemodynamic data, such as provided by ASL, increases the ability to disambiguate neuronal and vascular effects present in the BOLD signal.

The neuronal mechanisms underlying arithmetic calculations are not well understood but the differences between mental addition and subtraction could be particularly revealing. Using fMRI and dynamic causal modeling (DCM), this study aimed to identify the distinct neuronal architectures engaged by the cognitive processes of simple addition and subtraction. Our results revealed significantly greater activation during subtraction in regions along the dorsal pathway, including the left inferior frontal gyrus (IFG), middle portion of dorsolateral prefrontal cortex (mDLPFC), and supplementary motor area (SMA), compared with addition. Subsequent analysis of the underlying changes in connectivity - with DCM - revealed a common circuit processing basic (numeric) attributes and the retrieval of arithmetic facts. However, DCM showed that addition was more likely to engage (numeric) retrieval-based circuits in the left hemisphere, while subtraction tended to draw on (magnitude) processing in bilateral parietal cortex, especially the right intraparietal sulcus (IPS). Our findings endorse previous hypotheses about the differences in strategic implementation, dominant hemisphere, and the neuronal circuits underlying addition and subtraction. Moreover, for simple arithmetic, our connectivity results suggest that subtraction calls on more complex processing than addition: auxiliary phonological, visual, and motor processes, for representing numbers, were engaged by subtraction, relative to addition. Hum Brain Mapp 38:3210-3225, 2017. © 2017 Wiley Periodicals, Inc.

The “dysconnection hypothesis” of psychosis suggests that a disruption of functional integration underlies cognitive deficits and clinical symptoms. Impairments in the P300 potential are well documented in psychosis. Intrinsic (self-)connectivity in a frontoparietal cortical hierarchy during a P300 experiment was investigated. Dynamic Causal Modeling was used to estimate how evoked activity results from the dynamics of coupled neural populations and how neural coupling changes with the experimental factors. Twenty-four patients with psychotic disorder, twenty-four unaffected relatives, and twenty-five controls underwent EEG recordings during an auditory oddball paradigm. Sixteen frontoparietal network models (including primary auditory, superior parietal, and superior frontal sources) were analyzed and an optimal model of neural coupling, explaining diagnosis and genetic risk effects, as well as their interactions with task condition were identified. The winning model included changes in connectivity at all three hierarchical levels. Patients showed decreased self-inhibition—that is, increased cortical excitability—in left superior frontal gyrus across task conditions, compared with unaffected participants. Relatives had similar increases in excitability in left superior frontal and right superior parietal sources, and a reversal of the normal synaptic gain changes in response to targets relative to standard tones. It was confirmed that both subjects with psychotic disorder and their relatives show a context-independent loss of synaptic gain control at the highest hierarchy levels. The relatives also showed abnormal gain modulation responses to task-relevant stimuli. These may be caused by NMDA-receptor and/or GABAergic pathologies that change the excitability of superficial pyramidal cells and may be a potential biological marker for psychosis.

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.

Highlights • We decouple seen and felt hand postures during action via virtual reality. • Vision of the hand is either task-relevant or a distractor. • Task-relevance of vision is reflected by in- or decreases of occipital beta power. • DCM suggests underlying changes in cortical (visual) excitability. • Occipital beta may indicate the contextual gating of visual action feedback., In sensorimotor integration, the brain needs to decide how its predictions should accommodate novel evidence by ‘gating’ sensory data depending on the current context. Here, we examined the oscillatory correlates of this process by recording magnetoencephalography (MEG) data during a new task requiring action under intersensory conflict. We used virtual reality to decouple visual (virtual) and proprioceptive (real) hand postures during a task in which the phase of grasping movements tracked a target (in either modality). Thus, we rendered visual information either task-relevant or a (to-be-ignored) distractor. Under visuo-proprioceptive incongruence, occipital beta power decreased (relative to congruence) when vision was task-relevant but increased when it had to be ignored. Dynamic causal modeling (DCM) revealed that this interaction was best explained by diametrical, task-dependent changes in visual gain. These results suggest a crucial role for beta oscillations in the contextual gating (i.e., gain or precision control) of visual vs proprioceptive action feedback, depending on current behavioral demands., Graphical abstract Image, graphical abstract

Models of coupled 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 establish, using model-generated data and simulations of coupled pendula, that phase-only models perform well only under weak coupling conditions. We also show that, despite their higher complexity, 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., arXiv admin note: text overlap with arXiv:1812.06315

Neural models describe brain activity at different scales, ranging from single cells to whole brain networks. Here, we attempt to reconcile models operating at the microscopic (compartmental) and mesoscopic (neural mass) scales to analyse data from microelectrode recordings of intralaminar neural activity. Although these two classes of models operate at different scales, it is relatively straightforward to create neural mass models of ensemble activity that are equipped with priors obtained after fitting data generated by detailed microscopic models. This provides generative (forward) models of measured neuronal responses that retain construct validity in relation to compartmental models. We illustrate our approach using cross spectral responses obtained from V1 during a visual perception paradigm that involved optogenetic manipulation of the basal forebrain. We find that the resulting neural mass model can distinguish between activity in distinct cortical layers – both with and without optogenetic activation – and that cholinergic input appears to enhance (disinhibit) superficial layer activity relative to deep layers. This is particularly interesting from the perspective of predictive coding, where neuromodulators are thought to boost prediction errors that ascend the cortical hierarchy., Highlights • Biophysical mean field model fitted to data from different cortical layers. • Neural mass model combined with a laminar-specific forward model. • Fitted V1 data with and without optogenetic activation of the basal forebrain. • Test of Predictive Coding based on the role of neuromodulation. • Increase in prediction error activity due to cholinergic effects.

In this opinion paper, we describe a combined view of functional and effective brain connectivity along with the free-energy principle for investigating persistent disruptions in brain networks of patients with focal epilepsy. These changes are likely reflected in effective connectivity along the cortical hierarchy and construct the basis of increased local functional connectivity in focal epilepsy. We propose a testable framework based on dynamic causal modelling and functional connectivity analysis with the capacity of explaining commonly observed connectivity changes during interictal periods. We then hypothesise their possible relation with disrupted free-energy minimisation in the Bayesian brain. This may offer a new approach for neuroimaging to specifically develop and address hypotheses regarding the network pathomechanisms underlying epileptic phenotypes., Highlights • A combined view of functional brain connectivity and dynamic causal modeling of persistent networks changes in focal epilepsy. • Increased functional connectivity in focal epilepsy is likely reflected in effective connectivity of cortical hierarchy. • A testable framework based on dynamic causal modeling and functional connectivity analysis during interictal periods.

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., Author summary Epilepsy is a disease that affects over 50 million people worldwide. Current treatments include dangerous surgical procedures in which brain connections are severed, or even in which entire problem brain regions are removed. Pharmaceutical options are available, but only about one third of patients are responsive. However, even in these cases the drugs can cause such severe side effects that the patients sometimes choose to suffer seizures. We are proposing an innovative treatment of epilepsy that could be achieved by using non-invasive electrical stimulation. Specifically, we show that stimulation should be applied in such a way as to mirror the activity in a problem brain region, by targeting its neighbouring areas. This counterintuitive approach is based on a mathematical model in which this mirroring strategy is applied pre-emptively, i.e. long before the seizure has a chance to set in. The hope is that future clinical trials will be able to use this model to lessen the effect of seizures, or even prevent them from occurring in the first place.

This paper asks whether integrating multimodal EEG and fMRI data offers a better characterisation of functional brain architectures than either modality alone. This evaluation rests upon a dynamic causal model that generates both EEG and fMRI data from the same neuronal dynamics. We introduce the use of Bayesian fusion to provide informative (empirical) neuronal priors – derived from dynamic causal modelling (DCM) of EEG data – for subsequent DCM of fMRI data. To illustrate this procedure, we generated synthetic EEG and fMRI timeseries for a mismatch negativity (or auditory oddball) paradigm, using biologically plausible model parameters (i.e., posterior expectations from a DCM of empirical, open access, EEG data). Using model inversion, we found that Bayesian fusion provided a substantial improvement in marginal likelihood or model evidence, indicating a more efficient estimation of model parameters, in relation to inverting fMRI data alone. We quantified the benefits of multimodal fusion with the information gain pertaining to neuronal and haemodynamic parameters – as measured by the Kullback-Leibler divergence between their prior and posterior densities. Remarkably, this analysis suggested that EEG data can improve estimates of haemodynamic parameters; thereby furnishing proof-of-principle that Bayesian fusion of EEG and fMRI is necessary to resolve conditional dependencies between neuronal and haemodynamic estimators. These results suggest that Bayesian fusion may offer a useful approach that exploits the complementary temporal (EEG) and spatial (fMRI) precision of different data modalities. We envisage the procedure could be applied to any multimodal dataset that can be explained by a DCM with a common neuronal parameterisation.

A fundamental question in systems neuroscience is how endogenous neuronal activity self-organizes during particular brain states. Recent neuroimaging studies have demonstrated systematic relationships between resting-state and task-induced functional connectivity (FC). In particular, continuous task studies, such as movie watching, speak to alterations in coupling among cortical regions and enhanced fluctuations in FC compared to the resting-state. This suggests that FC may reflect systematic and large-scale reorganization of functionally integrated responses while subjects are watching movies. In this study, we characterized fluctuations in FC during resting-state and movie-watching conditions. We found that the FC patterns induced systematically by movie-watching can be explained with a single principal component. These condition-specific FC fluctuations overlapped with inter-subject synchronization patterns in occipital and temporal brain regions. However, unlike inter-subject synchronization, condition-specific FC patterns were characterized by increased correlations within frontal brain regions and reduced correlations between frontal-parietal brain regions. We investigated these condition-specific functional variations as a shorter time scale, using time-resolved FC. The time-resolved FC showed condition-specificity over time; notably when subjects watched both the same and different movies. To explain self-organisation of global FC through the alterations in local dynamics, we used a large-scale computational model. We found that condition-specific reorganization of FC could be explained by local changes that engendered changes in FC among higher-order association regions, mainly in frontal and parietal cortices. ispartof: NEUROIMAGE vol:184 pages:335-348 ispartof: location:United States status: published

Seizure activity in EEG recordings can persist for hours with seizure dynamics changing rapidly over time and space. To characterise the spatiotemporal evolution of seizure activity, large data sets often need to be analysed. Dynamic causal modelling (DCM) can be used to estimate the synaptic drivers of cortical dynamics during a seizure; however, the requisite (Bayesian) inversion procedure is computationally expensive. In this note, we describe a straightforward procedure, within the DCM framework, that provides efficient inversion of seizure activity measured with non-invasive and invasive physiological recordings; namely, EEG/ECoG. We describe the theoretical background behind a Bayesian belief updating scheme for DCM. The scheme is tested on simulated and empirical seizure activity (recorded both invasively and non-invasively) and compared with standard Bayesian inversion. We show that the Bayesian belief updating scheme provides similar estimates of time-varying synaptic parameters, compared to standard schemes, indicating no significant qualitative change in accuracy. The difference in variance explained was small (less than 5%). The updating method was substantially more efficient, taking approximately 5–10 min compared to approximately 1–2 h. Moreover, the setup of the model under the updating scheme allows for a clear specification of how neuronal variables fluctuate over separable timescales. This method now allows us to investigate the effect of fast (neuronal) activity on slow fluctuations in (synaptic) parameters, paving a way forward to understand how seizure activity is generated., Highlights • We describe a DCM procedure that provides efficient inversion of seizure activity. • Similar accuracy but substantially more efficient compared to standard DCM methods. • Physiological fluctuations over different timescales can be specified. • This scheme should contribute to understanding seizure activity using DCM.

In this technical note, we derive two MCMC (Markov chain Monte Carlo) samplers for dynamic causal models (DCMs). Specifically, we use (a) Hamiltonian MCMC (HMC-E) where sampling is simulated using Hamilton’s equation of motion and (b) Langevin Monte Carlo algorithm (LMC-R and LMC-E) that simulates the Langevin diffusion of samples using gradients either on a Euclidean (E) or on a Riemannian (R) manifold. While LMC-R requires minimal tuning, the implementation of HMC-E is heavily dependent on its tuning parameters. These parameters are therefore optimised by learning a Gaussian process model of the time-normalised sample correlation matrix. This allows one to formulate an objective function that balances tuning parameter exploration and exploitation, furnishing an intervention-free inference scheme. Using neural mass models (NMMs)—a class of biophysically motivated DCMs—we find that HMC-E is statistically more efficient than LMC-R (with a Riemannian metric); yet both gradient-based samplers are far superior to the random walk Metropolis algorithm, which proves inadequate to steer away from dynamical instability., Highlights • We compare 2 gradient-based MCMC methods for parameter inference in neural mass models. • These methods are Hamiltonian and Langevin Monte Carlo (HMC-E and LMC-R). • LMC-R is computationally more efficient while HMC-E proves to be statistically most efficient. • We propose combining gradient-based and gradient-free samplers for efficient Bayesian inference.

The functional MRI (fMRI) signal is an indirect measure of neuronal activity. In order to deconvolve the neuronal activity from the experimental fMRI data, biophysical generative models have been proposed describing the link between neuronal activity and the cerebral blood flow (the neurovascular coupling), and further the hemodynamic response and the BOLD signal equation. These generative models have been employed both for single brain area deconvolution and to infer effective connectivity in networks of multiple brain areas. In the current paper, we introduce a new fMRI model inspired by experimental observations about the physiological underpinnings of the BOLD signal and compare it with the generative models currently used in dynamic causal modeling (DCM), a widely used framework to study effective connectivity in the brain. We consider three fundamental aspects of such generative models for fMRI: (i) an adaptive two-state neuronal model that accounts for a wide repertoire of neuronal responses during and after stimulation; (ii) feedforward neurovascular coupling that links neuronal activity to blood flow; and (iii) a balloon model that can account for vascular uncoupling between the blood flow and the blood volume. Finally, we adjust the parameterization of the BOLD signal equation for different magnetic field strengths. This paper focuses on the form, motivation and phenomenology of DCMs for fMRI and the characteristics of the various models are demonstrated using simulations. These simulations emphasize a more accurate modeling of the transient BOLD responses-such as adaptive decreases to sustained inputs during stimulation and the post-stimulus undershoot. In addition, we demonstrate using experimental data that it is necessary to take into account both neuronal and vascular transients to accurately model the signal dynamics of fMRI data. By refining the models of the transient responses, we provide a more informed perspective on the underlying neuronal process and offer new ways of inferring changes in local neuronal activity and effective connectivity from fMRI. (C) 2015 The Authors. Published by Elsevier Inc.

Dynamic causal modeling (DCM) is a Bayesian framework for inferring effective connectivity among brain regions from neuroimaging data. While the validity of DCM has been investigated in various previous studies, the reliability of DCM parameter estimates across sessions has been examined less systematically. Here, we report results of a software comparison with regard to test-retest reliability of DCM for fMRI, using a challenging scenario where complex models with many parameters were applied to relatively few data points. Specifically, we examined the reliability of different DCM implementations (in terms of the intra-class correlation coefficient, ICC) based on fMRI data from 35 human subjects performing a simple motor task in two separate sessions, one month apart. We constructed DCMs of motor regions with fair to excellent reliability of conventional activation measures. Using classical DCM (cDCM) in SPM5, we found that the test-retest reliability of DCM results was high, both concerning the model evidence (ICC=0.94) and the model parameter estimates (median ICC=0.47). However, when using a more recent DCM version (DCM10 in SPM8), test-retest reliability was reduced notably. Analyses indicated that, in our particular case, the prior distributions played a crucial role in this change in reliability across software versions. Specifically, when using cDCM priors for model inversion in DCM10, this not only restored reliability but yielded even better results than in cDCM. Analyzing each component of the objective function in DCM, we found a selective change in the reliability of posterior mean estimates. This suggests that tighter regularization afforded by cDCM priors reduces the possibility of local extrema in the objective function. We conclude this paper with an outlook to ongoing developments for overcoming the software-dependency of reliability observed in this study, including global optimization and empirical Bayesian procedures.

Individual differences in cognitive efficiency, particularly in relation to working memory (WM), have been associated both with personality dimensions that reflect enduring regularities in brain configuration, and with short-term neural plasticity, that reflects task-related changes in brain connectivity. To elucidate the relationship of these two divergent mechanisms, we tested the hypothesis that personality dimensions, which reflect enduring aspects of brain configuration, inform about the neurobiological framework within which short-term, task-related plasticity, as measured by effective connectivity, can be facilitated or constrained. As WM consistently engages the dorsolateral prefrontal (DLPFC), parietal (PAR), and anterior cingulate cortex (ACC), we specified a WM network model with bidirectional, ipsilateral, and contralateral connections between these regions from a functional magnetic resonance imaging dataset obtained from 40 healthy adults while performing the 3-back WM task. Task-related effective connectivity changes within this network were estimated using Dynamic Causal Modelling. Personality was evaluated along the major dimensions of Neuroticism, Extraversion, Openness to Experience, Agreeableness, and Conscientiousness. Only two dimensions were relevant to task-dependent effective connectivity. Neuroticism and Conscientiousness respectively constrained and facilitated neuroplastic responses within the WM network. These results suggest individual differences in cognitive efficiency arise from the interplay between enduring and short-term plasticity in brain configuration.

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.

Our behavior entails a flexible and context-sensitive interplay between brain areas to integrate information according to goal-directed requirements. However, the neural mechanisms governing the entrainment of functionally specialized brain areas remain poorly understood. In particular, the question arises whether observed changes in the activity for different cognitive conditions are explained by modifications of the (intrinsic) inputs or (extrinsic) recurrent connectivity? Here, we show that fMRI transitions over successive time points - a new version of dynamic functional connectivity compared to usual definitions - convey information about the task performed by 19 subjects, namely watching a movie versus a black screen (rest). We use a theoretical framework that extracts this information to characterize the mechanisms underlying cortical coordination at the whole brain level. Our approach relies on our recent network model that was introduced for resting-state fMRI and estimates both the amplitude of input fluctuating activity to each of the 66 cortical regions and the strengths of cortico-cortical interactions between them. In our model, the change of condition does modify the (extrinsic) recurrent connectivity as much as the (intrinsic) inputs, which are putatively related to the movie stimulus. However, detailed changes in connectivity preserve a balance in the propagation of fluctuating activity and select specific pathways that integrate sensory information from the visual and auditory systems to high-level regions of the brain. These findings speak to a dynamic functional integration that underlies the hierarchical processing in the brain.

Psychostimulants, such as methylphenidate, are thought to improve information processing in motivation-reward and attention-activation networks by enhancing the effects of more relevant signals and suppressing those of less relevant ones; however the nature of such reciprocal influences remains poorly understood. To explore this question, we tested the effect of methylphenidate on performance and associated brain activity in the Anticipation, Conflict, Reward (ACR) task. Sixteen healthy adult volunteers, ages 21-45, were scanned twice using functional magnetic resonance imaging (fMRI) as they performed the ACR task under placebo and methylphenidate conditions. A three-way repeated measures analysis of variance, with cue (reward vs. non-reward), target (congruent vs. incongruent) and medication condition (methylphenidate vs. placebo) as the factors, was used to analyze behaviors on the task. Blood oxygen level dependent (BOLD) signals, reflecting task-related neural activity, were evaluated using linear contrasts. Participants exhibited significantly greater accuracy in the methylphenidate condition than the placebo condition. Compared with placebo, the methylphenidate condition also was associated with lesser task-related activity in components of attention-activation systems irrespective of the reward cue, and less task-related activity in components of the reward-motivation system, particularly the insula, during reward trials irrespective of target difficulty. These results suggest that methylphenidate enhances task performance by improving efficiency of information processing in both reward-motivation and in attention-activation systems.

In this paper, we revisit the problem of Bayesian model selection (BMS) at the group level. We originally addressed this issue in Stephan et al. (2009) , where models are treated as random effects that could differ between subjects, with an unknown population distribution. Here, we extend this work, by (i) introducing the Bayesian omnibus risk (BOR) as a measure of the statistical risk incurred when performing group BMS, (ii) highlighting the difference between random effects BMS and classical random effects analyses of parameter estimates, and (iii) addressing the problem of between group or condition model comparisons. We address the first issue by quantifying the chance likelihood of apparent differences in model frequencies. This leads to the notion of protected exceedance probabilities. The second issue arises when people want to ask “whether a model parameter is zero or not” at the group level. Here, we provide guidance as to whether to use a classical second-level analysis of parameter estimates, or random effects BMS. The third issue rests on the evidence for a difference in model labels or frequencies across groups or conditions. Overall, we hope that the material presented in this paper finesses the problems of group-level BMS in the analysis of neuroimaging and behavioural data.

Reward outcome signalling in the sensory cortex is held as important for linking stimuli to their consequences and for modulating perceptual learning in response to incentives. Evidence for reward outcome signalling has been found in sensory regions including the visual, auditory and somatosensory cortices across a range of different paradigms, but it is unknown whether the population of neurons signalling rewarding outcomes are the same as those processing predictive stimuli. We addressed this question using a multivariate analysis of high-resolution functional magnetic resonance imaging (fMRI), in a task where subjects were engaged in instrumental learning with visual predictive cues and auditory signalled reward feedback. We found evidence that outcome signals in sensory regions localise to the same areas involved in stimulus processing. These outcome signals are non-specific and we show that the neuronal populations involved in stimulus representation are not their exclusive target, in keeping with theoretical models of value learning. Thus, our results reveal one likely mechanism through which rewarding outcomes are linked to predictive sensory stimuli, a link that may be key for both reward and perceptual learning., Highlights • Rewarding outcomes reactivate sensory regions involved in stimulus representation. • Within these regions reward is not directed specifically to stimulus-coding neurons. • This resembles a teaching signal predicted by reinforcement learning models.

Lohmann et al. (this issue) make three unremarkable observations about model selection and use them to critique dynamic causal modelling—a Bayesian model selection procedure based on causal models of dynamical systems (Marreiros et al., 2010). In this response, we unpack their misconceptions and try to answer their questions.

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., Highlights ► Effects of hippocampal sclerosis (HS) on effective connectivity were evaluated. ► HS induced a network adjustment in the ipsilesional and the contralesional hemisphere. ► Strength in ipsilesional and contralesional couplings was negatively related. ► Contralesional couplings were inversely related to performance and level of education. ► These findings support the notion that network alterations underlie cognitive decline.

Conventionally, set-level inference on statistical parametric maps (SPMs) is based on the topological features of an excursion set above some threshold—for example, the number of clusters or Euler characteristic. The expected Euler characteristic—under the null hypothesis—can be predicted from an intrinsic measure or volume of the SPM, such as the resel counts or the Lipschitz–Killing curvatures (LKC). We propose a new approach that performs a null hypothesis omnibus test on an SPM, by testing whether its intrinsic volume (described by LKC coefficients) is different from the volume of the underlying residual fields: intuitively, whether the number of peaks in the statistical field (testing for signal) and the residual fields (noise) are consistent or not. Crucially, this new test requires no arbitrary feature-defining threshold but is nevertheless sensitive to distributed or spatially extended patterns. We show the similarities between our approach and conventional topological inference—in terms of false positive rate control and sensitivity to treatment effects—in two and three dimensional simulations. The test consistently improves on classical approaches for moderate (> 20) degrees of freedom. We also demonstrate the application to real data and illustrate the comparison of the expected and observed Euler characteristics over the complete threshold range., Highlights ► We calculate SPM intrinsic volume through the Euler Characteristic over thresholds. ► We compare these volume estimates to those obtained using classical methods. ► A multivariate test on these volume estimates gives a threshold free test of the null.

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., NeuroImage, 145, ISSN:1053-8119, ISSN:1095-9572

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., Highlights ► Novel model of basal ganglia-cortical system based on data from Parkinson’s patients. ► Model exposes criticality of connections to and from subthalamic nucleus in PD. ► Model able to replicate the effects of known surgical interventions in PD. ► Allows in silico tests of candidate therapeutic interventions in circuit disorders.

We demonstrate the capacity of dynamic causal modeling to characterize the nonlinear coupling among cortical sources that underlie time-frequency modulations in MEG data. Our experimental task involved the mental rotation of hand drawings that ten subjects used to decide if it was a right or left hand. Reaction times were shorter when the stimuli were presented with a small rotation angle (fast responses) compared to a large rotation angle (slow responses). The grand-averaged data showed that in both cases performance was accompanied by a marked increase in gamma activity in occipital areas and a concomitant decrease in alpha and beta power in occipital and motor regions. Modeling directed (cross) frequency interactions between the two regions revealed that after the stimulus induced a gamma increase and beta decrease in occipital regions, interactions with the motor area served to attenuate these modulations. The difference between fast and slow behavioral responses was manifest as an altered coupling strength in both forward and backward connections, which led to a less pronounced attenuation for more difficult (slow reaction time) trials. This was mediated by a (backwards) beta to gamma coupling from motor till occipital sources, whereas other interactions were mainly within the same frequency. Results are consistent with the theory of predictive coding and suggest that during motor imagery, the influence of motor areas on activity in occipital cortex co-determines performance. Our study illustrates the benefit of modeling experimental responses in terms of a generative model that can disentangle the contributions of intra-areal vis-à-vis inter-areal connections to time-frequency modulations during task performance. © 2013 Elsevier Inc.