Search

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 (A rest ) and task states (A task ), (ii) cluster phenotypes of task-related changes in effective connectivity (B task ) across participants, (iii) identify edges (B task ) showing high inter-individual effective connectivity differences and (iv) associate reaction times with the similarity between B task and A rest in these edges. We found a strong correlation between A rest and A task over subjects but a marked difference between B task and A rest . We further observed a strong clustering of individuals in terms of B task , which was not apparent in A rest . The task-related effective connectivity B task 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 B task and A rest 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. [ABSTRACT FROM AUTHOR]

Many analysis methods exist to extract graphs of functional connectivity from neuronal networks. Confidence in the results is limited because, (i) different methods give different results, (ii) parameter setting directly influences the final result, and (iii) systematic evaluation of the results is not always performed. Here, we introduce MULAN (MULtiple method ANalysis), which assumes an ensemble based approach combining multiple analysis methods and fuzzy logic to extract graphs with the most probable structure. In order to reduce the dependency on parameter settings, we determine the best set of parameters using a genetic algorithm on simulated datasets, whose temporal structure is similar to the experimental one. After a validation step, the selected set of parameters is used to analyze experimental data. The final step cross-validates experimental subsets of data and provides a direct estimate of the most likely graph and our confidence in the proposed connectivity. A systematic evaluation validates our strategy against empirical stereotactic electroencephalography (SEEG) and functional magnetic resonance imaging (fMRI) data. [ABSTRACT FROM AUTHOR]

The neural substrate of bodily ownership can be disclosed by the rubber hand illusion (RHI); namely, the illusory self-attribution of an artificial hand that is induced by synchronous tactile stimulation of the subject's hand that is hidden from view. Previous studies have pointed to the premotor cortex (PMC) as a pivotal area in such illusions. To investigate the effective connectivity between – and within – sensory and premotor areas involved in bodily perceptions, we used dynamic causal modeling of touch-evoked responses in 13 healthy subjects. Each subject's right hand was stroked while viewing their own hand (“REAL”), or an artificial hand presented in an anatomically plausible (“CONGRUENT”) or implausible (“INCONGRUENT”) position. Bayesian model comparison revealed strong evidence for a differential involvement of the PMC in the generation of touch-evoked responses under the three conditions, confirming a crucial role of PMC in bodily self-attribution. In brief, the extrinsic (forward) connection from left occipital cortex to left PMC was stronger for CONGRUENT and INCONGRUENT as compared to REAL, reflecting the augmentation of bottom-up visual input when multisensory integration is challenged. Crucially, intrinsic connectivity in the primary somatosensory cortex (S1) was attenuated in the CONGRUENT condition, during the illusory percept. These findings support predictive coding models of the functional architecture of multisensory integration (and attenuation) in bodily perceptual experience. [ABSTRACT FROM AUTHOR]

This technical note describes some Bayesian procedures for the analysis of group studies that use nonlinear models at the first (within-subject) level – e.g., dynamic causal models – and linear models at subsequent (between-subject) levels. Its focus is on using Bayesian model reduction to finesse the inversion of multiple models of a single dataset or a single (hierarchical or empirical Bayes) model of multiple datasets. These applications of Bayesian model reduction allow one to consider parametric random effects and make inferences about group effects very efficiently (in a few seconds). We provide the relatively straightforward theoretical background to these procedures and illustrate their application using a worked example. This example uses a simulated mismatch negativity study of schizophrenia. We illustrate the robustness of Bayesian model reduction to violations of the (commonly used) Laplace assumption in dynamic causal modelling and show how its recursive application can facilitate both classical and Bayesian inference about group differences. Finally, we consider the application of these empirical Bayesian procedures to classification and prediction. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

In this technical note we compare the performance of four gradient-free MCMC samplers (random walk Metropolis sampling, slice-sampling, adaptive MCMC sampling and population-based MCMC sampling with tempering) in terms of the number of independent samples they can produce per unit computational time. For the Bayesian inversion of a single-node neural mass model, both adaptive and population-based samplers are more efficient compared with random walk Metropolis sampler or slice-sampling; yet adaptive MCMC sampling is more promising in terms of compute time. Slice-sampling yields the highest number of independent samples from the target density — albeit at almost 1000% increase in computational time, in comparison to the most efficient algorithm (i.e., the adaptive MCMC sampler). [ABSTRACT FROM AUTHOR]

Human faces may signal relevant information and are therefore analysed rapidly and effectively by the brain. However, the precise mechanisms and pathways involved in rapid face processing are unclear. One view posits a role for a subcortical connection between early visual sensory regions and the amygdala, while an alternative account emphasises cortical mediation. To adjudicate between these functional architectures, we recorded magnetoencephalographic (MEG) evoked fields in human subjects to presentation of faces with varying emotional valence. Early brain activity was better explained by dynamic causal models containing a direct subcortical connection to the amygdala irrespective of emotional modulation. At longer latencies, models without a subcortical connection had comparable evidence. Hence, our results support the hypothesis that a subcortical pathway to the amygdala plays a role in rapid sensory processing of faces, in particular during early stimulus processing. This finding contributes to an understanding of the amygdala as a behavioural relevance detector. [ABSTRACT FROM AUTHOR]

This technical paper offers a critical re-evaluation of (spectral) Granger causality measures in the analysis of biological timeseries. Using realistic (neural mass) models of coupled neuronal dynamics, we evaluate the robustness of parametric and nonparametric Granger causality. Starting from a broad class of generative (state-space) models of neuronal dynamics, we show how their Volterra kernels prescribe the second-order statistics of their response to random fluctuations; characterised in terms of cross-spectral density, cross-covariance, autoregressive coefficients and directed transfer functions. These quantities in turn specify Granger causality — providing a direct (analytic) link between the parameters of a generative model and the expected Granger causality. We use this link to show that Granger causality measures based upon autoregressive models can become unreliable when the underlying dynamics is dominated by slow (unstable) modes — as quantified by the principal Lyapunov exponent. However, nonparametric measures based on causal spectral factors are robust to dynamical instability. We then demonstrate how both parametric and nonparametric spectral causality measures can become unreliable in the presence of measurement noise. Finally, we show that this problem can be finessed by deriving spectral causality measures from Volterra kernels, estimated using dynamic causal modelling. [ABSTRACT FROM AUTHOR]

This paper examines intrinsic brain networks in light of recent developments in the characterisation of resting state fMRI timeseries -- and simulations of neuronal fluctuations based upon the connectome. Its particular focus is on patterns or modes of distributed activity that underlie functional connectivity. We first demonstrate that the eigenmodes of functional connectivity - or covariance among regions or nodes - are the same as the eigenmodes of the underlying effective connectivity, provided we limit ourselves to symmetrical connections. This symmetry constraint is motivated by appealing to proximity graphs based upon multidimensional scaling. Crucially, the principal modes of functional connectivity correspond to the dynamically unstable modes of effective connectivity that decay slowly and show long term memory. Technically, these modes have small negative Lyapunov exponents that approach zero from below. Interestingly, the superposition of modes - whose exponents are sampled from a power law distribution - produces classical 1/f (scale free) spectra. We conjecture that the emergence of dynamical instability - that underlies intrinsic brain networks - is inevitable in any system that is separated from external states by a Markov blanket. This conjecture appeals to a free energy formulation of nonequilibrium steady-state dynamics. The common theme that emerges from these theoretical considerations is that endogenous fluctuations are dominated by a small number of dynamically unstable modes. We use this as the basis of a dynamic causal model (DCM) of resting state fluctuations -- as measured in terms of their complex cross spectra. In this model, effective connectivity is parameterised in terms of eigenmodes and their Lyapunov exponents -- that can also be interpreted as locations in a multidimensional scaling space. Model inversion provides not only estimates of edges or connectivity but also the topography and dimensionality of the underlying scaling space. Here, we focus on conceptual issues with simulated fMRI data and provide an illustrative application using an empirical multi-region timeseries. [ABSTRACT FROM AUTHOR]

Abstract: This technical note introduces a dynamic causal model (DCM) for resting state fMRI time series based upon observed functional connectivity—as measured by the cross spectra among different brain regions. This DCM is based upon a deterministic model that generates predicted crossed spectra from a biophysically plausible model of coupled neuronal fluctuations in a distributed neuronal network or graph. Effectively, the resulting scheme finds the best effective connectivity among hidden neuronal states that explains the observed functional connectivity among haemodynamic responses. This is because the cross spectra contain all the information about (second order) statistical dependencies among regional dynamics. In this note, we focus on describing the model, its relationship to existing measures of directed and undirected functional connectivity and establishing its face validity using simulations. In subsequent papers, we will evaluate its construct validity in relation to stochastic DCM and its predictive validity in Parkinson's and Huntington's disease. [Copyright &y& Elsevier]

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. [ABSTRACT FROM AUTHOR]

This technical note describes a variational or Bayesian implementation of representational similarity analysis (RSA) and pattern component modelling (PCM). It considers RSA and PCM as Bayesian model comparison procedures that assess the evidence for stimulus or condition-specific patterns of responses distributed over voxels or channels. On this view, one can use standard variational inference procedures to quantify the contributions of particular patterns to the data, by evaluating second-order parameters or hyperparameters. Crucially, this allows one to use parametric empirical Bayes (PEB) to infer which patterns are consistent among subjects. At the between-subject level, one can then assess the evidence for different (combinations of) hypotheses about condition-specific effects using Bayesian model comparison. Alternatively, one can select a single hypothesis that best explains the pattern of responses using Bayesian model selection. This note rehearses the technical aspects of within and between-subject RSA using a worked example, as implemented in the Statistical Parametric Mapping (SPM) software. En route , we highlight the connection between univariate and multivariate analyses of neuroimaging data and the sorts of analyses that are possible using component modelling and representational similarity analysis. • We introduce variational RSA, a method for multivariate analysis in neuroimaging. • This treats RSA a standard covariance component estimation problem. • An efficient estimation scheme, variational Laplace, is used to estimate parameters. • Bayesian model comparison is used to optimally test for mixtures of effects. • We illustrate the approach using simulated and empirical fMRI data. [ABSTRACT FROM AUTHOR]

Abstract Context-sensitive and activity-dependent fluctuations in connectivity underlie functional integration in the brain and have been studied widely in terms of synaptic plasticity, learning and condition-specific (e.g., attentional) modulations of synaptic efficacy. This dynamic aspect of brain connectivity has recently attracted a lot of attention in the resting state fMRI community. To explain dynamic functional connectivity in terms of directed effective connectivity among brain regions, we introduce a novel method to identify dynamic effective connectivity using spectral dynamic causal modelling (spDCM). We used parametric empirical Bayes (PEB) to model fluctuations in directed coupling over consecutive windows of resting state fMRI time series. Hierarchical PEB can model random effects on connectivity parameters at the second (between-window) level given connectivity estimates from the first (within-window) level. In this work, we used a discrete cosine transform basis set or eigenvariates (i.e., expression of principal components) to model fluctuations in effective connectivity over windows. We evaluated the ensuing dynamic effective connectivity in terms of the consistency of baseline connectivity within default mode network (DMN), using the resting state fMRI from Human Connectome Project (HCP). To model group-level baseline and dynamic effective connectivity for DMN, we extended the PEB approach by conducting a multilevel PEB analysis of between-session and between-subject group effects. Model comparison clearly spoke to dynamic fluctuations in effective connectivity – and the dynamic functional connectivity these changes explain. Furthermore, baseline effective connectivity was consistent across independent sessions – and notably more consistent than estimates based upon conventional models. This work illustrates the advantage of hierarchical modelling with spDCM, in characterizing the dynamics of effective connectivity. Highlights • We describe efficient estimation of dynamics in resting state effective connectivity. • Spectral DCM and PEB are used to model fluctuations in neuronal coupling over time. • Dynamics in responses are explained in terms of its causes (effective connectivity). • Baseline and dynamic components of the default mode connectivity are identified. [ABSTRACT FROM AUTHOR]

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

• We model selective visual attention in macaques using Dynamic Causal Modeling. • Intrinsic V1 modulation can explain attention effects in V1-V4 communication. • Modulation of superficial and granular inhibition is key to induce the effects. • Those modulations increase V1-V4 communication in a feedforward manner. 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. [Display omitted] [ABSTRACT FROM AUTHOR]

Dynamic Causal Modelling (DCM) is the predominant method for inferring effective connectivity from neuroimaging data. In the 15 years since its introduction, the neural models and statistical routines in DCM have developed in parallel, driven by the needs of researchers in cognitive and clinical neuroscience. In this guide, we step through an exemplar fMRI analysis in detail, reviewing the current implementation of DCM and demonstrating recent developments in group-level connectivity analysis. In the appendices, we detail the theory underlying DCM and the assumptions (i.e., priors) in the models. In the first part of the guide (current paper), we focus on issues specific to DCM for fMRI. This is accompanied by all the necessary data and instructions to reproduce the analyses using the SPM software. In the second part (in a companion paper), we move from subject-level to group-level modelling using the Parametric Empirical Bayes framework, and illustrate how to test for commonalities and differences in effective connectivity across subjects, based on imaging data from any modality. • This guide walks through a group effective connectivity study using DCM and PEB. • Part 1, presented here, covers first level analysis using DCM for fMRI. • It clarifies the specific neural and haemodynamic models in DCM and their priors. • An accompanying dataset is provided with step-by-step analysis instructions. [ABSTRACT FROM AUTHOR]

This paper provides a worked example of using Dynamic Causal Modelling (DCM) and Parametric Empirical Bayes (PEB) to characterise inter-subject variability in neural circuitry (effective connectivity). It steps through an analysis in detail and provides a tutorial style explanation of the underlying theory and assumptions (i.e, priors). The analysis procedure involves specifying a hierarchical model with two or more levels. At the first level, state space models (DCMs) are used to infer the effective connectivity that best explains a subject's neuroimaging timeseries (e.g. fMRI, MEG, EEG). Subject-specific connectivity parameters are then taken to the group level, where they are modelled using a General Linear Model (GLM) that partitions between-subject variability into designed effects and additive random effects. The ensuing (Bayesian) hierarchical model conveys both the estimated connection strengths and their uncertainty (i.e., posterior covariance) from the subject to the group level; enabling hypotheses to be tested about the commonalities and differences across subjects. This approach can also finesse parameter estimation at the subject level, by using the group-level parameters as empirical priors. The preliminary first level (subject specific) DCM for fMRI analysis is covered in a companion paper. Here, we detail group-level analysis procedures that are suitable for use with data from any neuroimaging modality. This paper is accompanied by an example dataset, together with step-by-step instructions demonstrating how to reproduce the analyses. • This guide walks through a group effective connectivity study using DCM and PEB. • It explains recently developed tools for hierarchical Bayesian modelling. • The appendices clarify the technical detail of the PEB framework and its priors. • An accompanying dataset is provided with step-by-step analysis instructions. [ABSTRACT FROM AUTHOR]

• We developed neurochemistry-enriched dynamic causal model of resting-state MEG data using MRS measurements of neurotransmitter concentrations. • Neurotransmitter levels, place prior constraints on synaptic connectivity, in the biologically informed model of MEG data. • Comparing evidence of group DCMs using Bayesian model selection procedure reveal synaptic connections that are likely to be informed by MRS data. 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. [ABSTRACT FROM AUTHOR]

Abstract Parkinson's disease (PD) is a neurodegenerative condition in which aberrant oscillatory synchronization of neuronal activity at beta frequencies (15–35 Hz) across the cortico-basal ganglia-thalamocortical circuit is associated with debilitating motor symptoms, such as bradykinesia and rigidity. Mounting evidence suggests that the magnitude of beta synchrony in the parkinsonian state fluctuates over time, but the mechanisms by which thalamocortical circuitry regulates the dynamic properties of cortical beta in PD are poorly understood. Using the recently developed generic Dynamic Causal Modelling (DCM) framework, we recursively optimized a set of plausible models of the thalamocortical circuit (n = 144) to infer the neural mechanisms that best explain the transitions between low and high beta power states observed in recordings of field potentials made in the motor cortex of anesthetized Parkinsonian rats. Bayesian model comparison suggests that upregulation of cortical rhythmic activity in the beta-frequency band results from changes in the coupling strength both between and within the thalamus and motor cortex. Specifically, our model indicates that high levels of cortical beta synchrony are mainly achieved by a delayed (extrinsic) input from thalamic relay cells to deep pyramidal cells and a fast (intrinsic) input from middle pyramidal cells to superficial pyramidal cells. From a clinical perspective, our study provides insights into potential therapeutic strategies that could be utilized to modulate the network mechanisms responsible for the enhancement of cortical beta in PD. Specifically, we speculate that cortical stimulation aimed to reduce the enhanced excitatory inputs to either the superficial or deep pyramidal cells could be a potential non-invasive therapeutic strategy for PD. Highlights • Coupling changes within and between circuit nodes lead to cortical beta enhancement. • Input propagation delays play a crucial role in the up-regulation of cortical beta. • Beta power could be modulated by altering lamina specific inputs. [ABSTRACT FROM AUTHOR]

Abstract 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. Highlights • We present a generic framework that affords a flexible implementation of DCM. • Different cortical and subcortical nodes can now be combined within a network. • User specified modulatory effects on intrinsic connections are accommodated. • We used a cortical-basal ganglia DCM to study beta oscillations in Parkinsonism. • Dopamine reduces synaptic efficacy within hyper and indirect pathways. [ABSTRACT FROM AUTHOR]

Abstract 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 regional activity for different cognitive conditions are explained by modifications of the inputs to the brain or its connectivity? We observe that transitions of fMRI activity between areas convey information about the tasks performed by 19 subjects, watching a movie versus a black screen (rest). We use a model-based framework that explains this spatiotemporal functional connectivity pattern by the local variability for 66 cortical regions and the network effective connectivity between them. We find that, among the estimated model parameters, movie viewing affects to a larger extent the local activity, which we interpret as extrinsic changes related to the increased stimulus load. However, detailed changes in the effective connectivity preserve a balance in the propagating activity and select specific pathways such that high-level brain regions integrate visual and auditory information, in particular boosting the communication between the two brain hemispheres. These findings speak to a dynamic coordination underlying the functional integration in the brain. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The neural processing and experience of pain are influenced by both expectations and attention. For example, the amplitude of event-related pain responses is enhanced by both novel and unexpected pain, and by moving the focus of attention towards a painful stimulus. Under predictive coding, this congruence can be explained by appeal to a precision-weighting mechanism, which mediates bottom-up and top-down attentional processes by modulating the influence of feedforward and feedback signals throughout the cortical hierarchy. The influence of expectation and attention on pain processing can be mapped onto changes in effective connectivity between or within specific neuronal populations, using a canonical microcircuit (CMC) model of hierarchical processing. We thus implemented a CMC within dynamic causal modelling for magnetoencephalography in human subjects, to investigate how expectation violation and attention to pain modulate intrinsic (within-source) and extrinsic (between-source) connectivity in the somatosensory hierarchy. This enabled us to establish whether both expectancy and attentional processes are mediated by a similar precision-encoding mechanism within a network of somatosensory, frontal and parietal sources. We found that both unexpected and attended pain modulated the gain of superficial pyramidal cells in primary and secondary somatosensory cortex. This modulation occurred in the context of increased lateralized recurrent connectivity between somatosensory and fronto-parietal sources, driven by unexpected painful occurrences. Finally, the strength of effective connectivity parameters in S1, S2 and IFG predicted individual differences in subjective pain modulation ratings. Our findings suggest that neuromodulatory gain control in the somatosensory hierarchy underlies the influence of both expectation violation and attention on cortical processing and pain perception. [ABSTRACT FROM AUTHOR]

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

This paper shows that it is possible to estimate the subjective precision (inverse variance) of Bayesian beliefs during oculomotor pursuit. Subjects viewed a sinusoidal target, with or without random fluctuations in its motion. Eye trajectories and magnetoencephalographic (MEG) data were recorded concurrently. The target was periodically occluded, such that its reappearance caused a visual evoked response field (ERF). Dynamic causal modelling (DCM) was used to fit models of eye trajectories and the ERFs. The DCM for pursuit was based on predictive coding and active inference, and predicts subjects' eye movements based on their (subjective) Bayesian beliefs about target (and eye) motion. The precisions of these hierarchical beliefs can be inferred from behavioural (pursuit) data. The DCM for MEG data used an established biophysical model of neuronal activity that includes parameters for the gain of superficial pyramidal cells, which is thought to encode precision at the neuronal level. Previous studies (using DCM of pursuit data) suggest that noisy target motion increases subjective precision at the sensory level: i.e., subjects attend more to the target's sensory attributes. We compared (noisy motion-induced) changes in the synaptic gain based on the modelling of MEG data to changes in subjective precision estimated using the pursuit data. We demonstrate that imprecise target motion increases the gain of superficial pyramidal cells in V1 (across subjects). Furthermore, increases in sensory precision – inferred by our behavioural DCM – correlate with the increase in gain in V1, across subjects. This is a step towards a fully integrated model of brain computations, cortical responses and behaviour that may provide a useful clinical tool in conditions like schizophrenia. [ABSTRACT FROM AUTHOR]

The body underlies our sense of self, emotion, and agency. Signals arising from the skin convey warmth, social touch, and the physical characteristics of external stimuli. Surprising or unexpected tactile sensations can herald events of motivational salience, including imminent threats (e.g., an insect bite) and hedonic rewards (e.g., a caressing touch). Awareness of such events is thought to depend upon the hierarchical integration of body-related mismatch responses by the anterior insula. To investigate this possibility, we measured brain activity using functional magnetic resonance imaging, while healthy participants performed a roving tactile oddball task. Mass-univariate analysis demonstrated robust activations in limbic, somatosensory, and prefrontal cortical areas previously implicated in tactile deviancy, body awareness, and cognitive control. Dynamic Causal Modelling revealed that unexpected stimuli increased the strength of forward connections along a caudal to rostral hierarchy—projecting from thalamic and somatosensory regions towards insula, cingulate and prefrontal cortices. Within this ascending flow of sensory information, the AIC was the only region to show increased backwards connectivity to the somatosensory cortex, augmenting a reciprocal exchange of neuronal signals. Further, participants who rated stimulus changes as easier to detect showed stronger modulation of descending PFC to AIC connections by deviance. These results suggest that the AIC coordinates hierarchical processing of tactile prediction error. They are interpreted in support of an embodied predictive coding model where AIC mediated body awareness is involved in anchoring a global neuronal workspace. [ABSTRACT FROM AUTHOR]

Clinical assessments of brain function rely upon visual inspection of electroencephalographic waveform abnormalities in tandem with functional magnetic resonance imaging. However, no current technology proffers in vivo assessments of activity at synapses, receptors and ion-channels, the basis of neuronal communication. Using dynamic causal modeling we compared electrophysiological responses from two patients with distinct monogenic ion channelopathies and a large cohort of healthy controls to demonstrate the feasibility of assaying synaptic-level channel communication non-invasively. Synaptic channel abnormality was identified in both patients (100% sensitivity) with assay specificity above 89%, furnishing estimates of neurotransmitter and voltage-gated ion throughput of sodium, calcium, chloride and potassium. This performance indicates a potential novel application as an adjunct for clinical assessments in neurological and psychiatric settings. More broadly, these findings indicate that biophysical models of synaptic channels can be estimated non-invasively, having important implications for advancing human neuroimaging to the level of non-invasive ion channel assays. [ABSTRACT FROM AUTHOR]

Recently, there has been a lot of interest in characterising the connectivity of resting state brain networks. Most of the literature uses functional connectivity to examine these intrinsic brain networks. Functional connectivity has well documented limitations because of its inherent inability to identify causal interactions. Dynamic causal modelling (DCM) is a framework that allows for the identification of the causal (directed) connections among neuronal systems — known as effective connectivity. This technical note addresses the validity of a recently proposed DCM for resting state fMRI – as measured in terms of their complex cross spectral density – referred to as spectral DCM. Spectral DCM differs from (the alternative) stochastic DCM by parameterising neuronal fluctuations using scale free (i.e., power law) forms, rendering the stochastic model of neuronal activity deterministic. Spectral DCM not only furnishes an efficient estimation of model parameters but also enables the detection of group differences in effective connectivity, the form and amplitude of the neuronal fluctuations or both. We compare and contrast spectral and stochastic DCM models with endogenous fluctuations or state noise on hidden states. We used simulated data to first establish the face validity of both schemes and show that they can recover the model (and its parameters) that generated the data. We then used Monte Carlo simulations to assess the accuracy of both schemes in terms of their root mean square error. We also simulated group differences and compared the ability of spectral and stochastic DCMs to identify these differences. We show that spectral DCM was not only more accurate but also more sensitive to group differences. Finally, we performed a comparative evaluation using real resting state fMRI data (from an open access resource) to study the functional integration within default mode network using spectral and stochastic DCMs. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. 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. [ABSTRACT FROM AUTHOR]

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. • 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. [ABSTRACT FROM AUTHOR]

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. • 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. [ABSTRACT FROM AUTHOR]

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. • Multimodal DCM shows how the same neuronal activity causes multiple measurements. • Bayesian fusion of EEG/fMRI resolves conditional dependencies between parameters. • Information gain quantifies the added benefits of multimodal Bayesian fusion. [ABSTRACT FROM AUTHOR]

Models of coupled phase oscillators are used to describe a wide variety of phenomena in neuroimaging. These models typically rest on the premise that oscillator dynamics do not evolve beyond their respective limit cycles, and hence that interactions can be described purely in terms of phase differences. Whilst mathematically convenient, the restrictive nature of phase-only models can limit their explanatory power. We therefore propose a generalisation of dynamic causal modelling that incorporates both phase and amplitude. This allows for the separate quantifications of phase and amplitude contributions to the connectivity between neural regions. We show, using model-generated data and simulations of coupled pendula, that phase-amplitude models can describe strongly coupled systems more effectively than their phase-only counterparts. We relate our findings to four metrics commonly used in neuroimaging: the Kuramoto order parameter, cross-correlation, phase-lag index, and spectral entropy. We find that, with the exception of spectral entropy, the phase-amplitude model is able to capture all metrics more effectively than the phase-only model. We then demonstrate, using local field potential recordings in rodents and functional magnetic resonance imaging in macaque monkeys, that amplitudes in oscillator models play an important role in describing neural dynamics in anaesthetised brain states. • We present a dynamic causal model that incorporates both phase and amplitude. • Phase-amplitude models can describe strongly coupled systems more effectively. • Phase-amplitude models capture common neuroimaging metrics more accurately. • Oscillator amplitudes are important in quantifying anaesthetised brain states. [ABSTRACT FROM AUTHOR]