Dynamic causal modelling of lateral interactions in the visual cortex

Individual/Peak Gamma Frequency: What Do We Know?

In recent years, the concept of individualized measures of electroencephalographic (EEG) activity has emerged. Gamma-band activity plays an important role in many sensory and cognitive processes. Thus, peak frequency in the gamma range has received considerable attention. However, peak or individual gamma frequency (IGF) is rarely used as a primary measure of interest; consequently, little is known about its nature and functional significance. With this review, we attempt to comprehensively overview available information on the functional properties of peak gamma frequency, addressing its relationship with certain processes and/or modulation by various factors. Here, we show that IGFs seem to be related to various endogenous and exogenous factors. Broad functional aspects that are related to IGF might point to the differences in underlying mechanisms. Therefore, research utilizing different types of stimulation for IGF estimation and covering several functional aspects in the same population is required. Moreover, IGFs span a wide range of frequencies (30-100 Hz). This could be partly due to the variability of methods used to extract the measures of IGF. In order to overcome this issue, further studies aiming at the optimization of IGF extraction would be greatly beneficial.

Dynamic causal modelling shows a prominent role of local inhibition in alpha power modulation in higher visual cortex

During resting-state EEG recordings, alpha activity is more prominent over the posterior cortex in eyes-closed (EC) conditions compared to eyes-open (EO). In this study, we characterized the difference in spectra between EO and EC conditions using dynamic causal modelling. Specifically, we investigated the role of intrinsic and extrinsic connectivity-within the visual cortex-in generating EC-EO alpha power differences over posterior electrodes. The primary visual cortex (V1) and the bilateral middle temporal visual areas (V5) were equipped with bidirectional extrinsic connections using a canonical microcircuit. The states of four intrinsically coupled subpopulations-within each occipital source-were also modelled. Using Bayesian model selection, we tested whether modulations of the intrinsic connections in V1, V5 or extrinsic connections (or a combination thereof) provided the best evidence for the data. In addition, using parametric empirical Bayes (PEB), we estimated group averages under the winning model. Bayesian model selection showed that the winning model contained both extrinsic connectivity modulations, as well as intrinsic connectivity modulations in all sources. The PEB analysis revealed increased extrinsic connectivity during EC. Overall, we found a reduction in the inhibitory intrinsic connections during EC. The results suggest that the intrinsic modulations in V5 played the most important role in producing EC-EO alpha differences, suggesting an intrinsic disinhibition in higher order visual cortex, during EC resting state.

New Therapeutics in Alzheimer's Disease Longitudinal Cohort study (NTAD): study protocol

INTRODUCTION

With the pressing need to develop treatments that slow or stop the progression of Alzheimer's disease, new tools are needed to reduce clinical trial duration and validate new targets for human therapeutics. Such tools could be derived from neurophysiological measurements of disease.

METHODS AND ANALYSIS

The New Therapeutics in Alzheimer's Disease study (NTAD) aims to identify a biomarker set from magneto/electroencephalography that is sensitive to disease and progression over 1 year. The study will recruit 100 people with amyloid-positive mild cognitive impairment or early-stage Alzheimer's disease and 30 healthy controls aged between 50 and 85 years. Measurements of the clinical, cognitive and imaging data (magnetoencephalography, electroencephalography and MRI) of all participants will be taken at baseline. These measurements will be repeated after approximately 1 year on participants with Alzheimer's disease or mild cognitive impairment, and clinical and cognitive assessment of these participants will be repeated again after approximately 2 years. To assess reliability of magneto/electroencephalographic changes, a subset of 30 participants with mild cognitive impairment or early-stage Alzheimer's disease will also undergo repeat magneto/electroencephalography 2 weeks after baseline. Baseline and longitudinal changes in neurophysiology are the primary analyses of interest. Additional outputs will include atrophy and cognitive change and estimated numbers needed to treat each arm of simulated clinical trials of a future disease-modifying therapy.

ETHICS AND DATA STATEMENT

The study has received a favourable opinion from the East of England Cambridge Central Research Ethics Committee (REC reference 18/EE/0042). Results will be disseminated through internal reports, peer-reviewed scientific journals, conference presentations, website publication, submission to regulatory authorities and other publications. Data will be made available via the Dementias Platform UK Data Portal on completion of initial analyses by the NTAD study group.

Can the causal role of brain oscillations be studied through rhythmic brain stimulation?

Many studies have investigated the causal relevance of brain oscillations using rhythmic stimulation, either through direct-brain or sensory stimulation. Yet, how intrinsic rhythms interact with the externally generated rhythm is largely unknown. We presented a flickered (60 Hz) visual grating or its correspondent unflickered stimulus in a psychophysical change detection task during simultaneous magnetoencephalography recordings to humans to test the effect of visual entrainment on induced gamma oscillations. Notably, we generally observed the coexistence of the broadband induced gamma rhythm with the entrained flicker rhythm (reliably measured in each participant), with the peak frequency of the induced response remaining unaltered in approximately half of participants—relatively independently of their native frequency. However, flicker increased broadband induced gamma power, and this was stronger in participants with a native frequency closer to the flicker frequency (resonance) and led to strong phase entrainment. Presence of flicker did not change behavior itself but profoundly altered brain behavior correlates across the sample: While broadband induced gamma oscillations correlated with reaction times for unflickered stimuli (as known previously), for the flicker, the amplitude of the entrained flicker rhythm (but no more the induced oscillation) correlated with reaction times. This, however, strongly depended on whether a participant's peak frequency shifted to the entrained rhythm. Our results suggests that rhythmic brain stimulation leads to a coexistence of two partially independent oscillations with heterogeneous effects across participants on the downstream relevance of these rhythms for behavior. This may explain the inconsistency of findings related to external entrainment of brain oscillations and poses further questions toward causal manipulations of brain oscillations in general.

Conductance-based dynamic causal modeling: A mathematical review of its application to cross-power spectral densities

Dynamic Causal Modeling (DCM) is a Bayesian framework for inferring on hidden (latent) neuronal states, based on measurements of brain activity. Since its introduction in 2003 for functional magnetic resonance imaging data, DCM has been extended to electrophysiological data, and several variants have been developed. Their biophysically motivated formulations make these models promising candidates for providing a mechanistic understanding of human brain dynamics, both in health and disease. However, due to their complexity and reliance on concepts from several fields, fully understanding the mathematical and conceptual basis behind certain variants of DCM can be challenging. At the same time, a solid theoretical knowledge of the models is crucial to avoid pitfalls in the application of these models and interpretation of their results. In this paper, we focus on one of the most advanced formulations of DCM, i.e. conductance-based DCM for cross-spectral densities, whose components are described across multiple technical papers. The aim of the present article is to provide an accessible exposition of the mathematical background, together with an illustration of the model's behavior. To this end, we include step-by-step derivations of the model equations, point to important aspects in the software implementation of those models, and use simulations to provide an intuitive understanding of the type of responses that can be generated and the role that specific parameters play in the model. Furthermore, all code utilized for our simulations is made publicly available alongside the manuscript to allow readers an easy hands-on experience with conductance-based DCM.

Tiagabine induced modulation of oscillatory connectivity and activity match PET-derived, canonical GABA-A receptor distributions

As the most abundant inhibitory neurotransmitter in the mammalian brain, γ-aminobutyric acid (GABA) plays a crucial role in shaping the frequency and amplitude of oscillations, which suggests a role for GABA in shaping the topography of functional connectivity and activity. This study explored the effects of pharmacologically blocking the reuptake of GABA (increasing local concentrations) using the GABA transporter 1 (GAT1) blocker, tiagabine (15 mg). In a placebo-controlled crossover design, we collected resting magnetoencephalography (MEG) recordings from 15 healthy individuals prior to, and at 1-, 3- and 5- hours post, administration of tiagabine and placebo. We quantified whole brain activity and functional connectivity in discrete frequency bands. Drug-by-session (2 × 4) analysis of variance in connectivity revealed interaction and main effects. Post-hoc permutation testing of each post-drug recording vs. respective pre-drug baseline revealed consistent reductions of a bilateral occipital network spanning theta, alpha and beta frequencies, across 1- 3- and 5- hour recordings following tiagabine only. The same analysis applied to activity revealed significant increases across frontal regions, coupled with reductions in posterior regions, across delta, theta, alpha and beta frequencies. Crucially, the spatial distribution of tiagabine-induced changes overlap with group-averaged maps of the distribution of GABA_A receptors, from flumazenil (FMZ-V_T) PET, demonstrating a link between GABA availability, GABA_A receptor distribution, and low-frequency network oscillations. Our results indicate that the relationship between PET receptor distributions and MEG effects warrants further exploration, since elucidating the nature of this relationship may uncover electrophysiologically-derived maps of oscillatory activity as sensitive, time-resolved, and targeted receptor-mapping tools for pharmacological imaging.

Heritable functional architecture in human visual cortex

•

We analyzed retinotopic maps from monozygotic and dizygotic twin pairs.

•

Visual field maps in V1-V3 are more similar in monozygotic twins.

•

Heritability is greater in V1 and V3 for polar angle and population receptive field sizes.

•

Eccentricity maps show lesser degree of heritability.

•

Further evidence for link between cortical morphology and topology of retinotopic maps.

How much of the functional organization of our visual system is inherited? Here we tested the heritability of retinotopic maps in human visual cortex using functional magnetic resonance imaging. We demonstrate that retinotopic organization shows a closer correspondence in monozygotic (MZ) compared to dizygotic (DZ) twin pairs, suggesting a partial genetic determination. Using population receptive field (pRF) analysis to examine the preferred spatial location and selectivity of these neuronal populations, we estimate a heritability around 10–20% for polar angle preferences and spatial selectivity, as quantified by pRF size, in extrastriate areas V2 and V3. Our findings are consistent with heritability in both the macroscopic arrangement of visual regions and stimulus tuning properties of visual cortex. This could constitute a neural substrate for variations in a range of perceptual effects, which themselves have been found to be at least partially genetically determined. These findings also add convergent evidence for the hypothesis that functional map topology is linked with cortical morphology.

Differences in visually induced MEG oscillations reflect differences in deep cortical layer activity

Neural activity is organized at multiple scales, ranging from the cellular to the whole brain level. Connecting neural dynamics at different scales is important for understanding brain pathology. Neurological diseases and disorders arise from interactions between factors that are expressed in multiple scales. Here, we suggest a new way to link microscopic and macroscopic dynamics through combinations of computational models. This exploits results from statistical decision theory and Bayesian inference. To validate our approach, we used two independent MEG datasets. In both, we found that variability in visually induced oscillations recorded from different people in simple visual perception tasks resulted from differences in the level of inhibition specific to deep cortical layers. This suggests differences in feedback to sensory areas and each subject's hypotheses about sensations due to differences in their prior experience. Our approach provides a new link between non-invasive brain imaging data, laminar dynamics and top-down control.

Statistical decision theory and multiscale analyses of human brain data

BACKGROUND

In the era of Big Data, large scale electrophysiological data from animal and human studies are abundant. These data contain information at multiple spatiotemporal scales. However, current approaches for the analysis of electrophysiological data often focus on a single spatiotemporal scale only.

NEW METHOD

We discuss a multiscale approach for the analysis of electrophysiological data. This is based on combining neural models that describe brain data at different scales. It allows us to make laminar-specific inferences about neurobiological properties of cortical sources using non invasive human electrophysiology data.

RESULTS

We provide a mathematical proof of this approach using statistical decision theory. We also consider its extensions to brain imaging studies including data from the same subjects performing different tasks. As an illustration, we show that changes in gamma oscillations between different people might originate from differences in recurrent connection strengths of inhibitory interneurons in layers 5/6.

COMPARISON WITH EXISTING METHODS

This is a new approach that follows up on our recent work. It is different from other approaches where the scale of spatiotemporal dynamics is fixed.

CONCLUSIONS

We discuss a multiscale approach for the analysis of human MEG data. This uses a neural mass model that includes constraints informed by a compartmental model. This has two advantages. First, it allows us to find differences in cortical laminar dynamics and understand neurobiological properties like neuromodulation, excitation to inhibition balance etc. using non invasive data. Second, it allows us to validate macroscale models by exploiting animal data.

Working Memory Load Modulates Neuronal Coupling

There is a severe limitation in the number of items that can be held in working memory. However, the neurophysiological limits remain unknown. We asked whether the capacity limit might be explained by differences in neuronal coupling. We developed a theoretical model based on Predictive Coding and used it to analyze Cross Spectral Density data from the prefrontal cortex (PFC), frontal eye fields (FEF), and lateral intraparietal area (LIP). Monkeys performed a change detection task. The number of objects that had to be remembered (memory load) was varied (1-3 objects in the same visual hemifield). Changes in memory load changed the connectivity in the PFC-FEF-LIP network. Feedback (top-down) coupling broke down when the number of objects exceeded cognitive capacity. Thus, impaired behavioral performance coincided with a break-down of Prediction signals. This provides new insights into the neuronal underpinnings of cognitive capacity and how coupling in a distributed working memory network is affected by memory load.

Dissociable neural effects of temporal expectations due to passage of time and contextual probability

The human brain is equipped with complex mechanisms to track the changing probability of events in time. While the passage of time itself usually leads to a mounting expectation, context can provide additional information about when events are likely to happen. In this study we dissociate these two sources of temporal expectation in terms of their neural correlates and underlying brain connectivity patterns. We analysed magnetoencephalographic (MEG) data acquired from N = 24 healthy participants listening to auditory stimuli. These stimuli could be presented at different temporal intervals but occurred most often at intermediate intervals, forming a contextual probability distribution. Evoked MEG response amplitude was sensitive to both passage of time (time elapsed since the cue) and contextual probability, albeit at different latencies: the effects of passage of time were observed earlier than the effects of context. The underlying sources of MEG activity were also different across the two types of temporal prediction: the effects of passage of time were localised to early auditory regions and superior temporal gyri, while context was additionally linked to activity in inferior parietal cortices. Finally, these differences were modelled using biophysical (dynamic causal) modelling: passage of time was explained in terms of widespread gain modulation and decreased prediction error signalling at lower levels of the hierarchy, while contextual expectation led to more localised gain modulation and decreased prediction error signalling at higher levels of the hierarchy. These results present a comprehensive account of how independent sources of temporal prediction may be differentially expressed in cortical circuits.

Bayesian Modelling of Induced Responses and Neuronal Rhythms

Neural rhythms or oscillations are ubiquitous in neuroimaging data. These spectral responses have been linked to several cognitive processes; including working memory, attention, perceptual binding and neuronal coordination. In this paper, we show how Bayesian methods can be used to finesse the ill-posed problem of reconstructing-and explaining-oscillatory responses. We offer an overview of recent developments in this field, focusing on (i) the use of MEG data and Empirical Bayes to build hierarchical models for group analyses-and the identification of important sources of inter-subject variability and (ii) the construction of novel dynamic causal models of intralaminar recordings to explain layer-specific activity. We hope to show that electrophysiological measurements contain much more spatial information than is often thought: on the one hand, the dynamic causal modelling of non-invasive (low spatial resolution) electrophysiology can afford sub-millimetre (hyper-acute) resolution that is limited only by the (spatial) complexity of the underlying (dynamic causal) forward model. On the other hand, invasive microelectrode recordings (that penetrate different cortical layers) can reveal laminar-specific responses and elucidate hierarchical message passing and information processing within and between cortical regions at a macroscopic scale. In short, the careful and biophysically grounded modelling of sparse data enables one to characterise the neuronal architectures generating oscillations in a remarkable detail.

Not All Predictions Are Equal: "What" and "When" Predictions Modulate Activity in Auditory Cortex through Different Mechanisms

Using predictions based on environmental regularities is fundamental for adaptive behavior. While it is widely accepted that predictions across different stimulus attributes (e.g., time and content) facilitate sensory processing, it is unknown whether predictions across these attributes rely on the same neural mechanism. Here, to elucidate the neural mechanisms of predictions, we combine invasive electrophysiological recordings (human electrocorticography in 4 females and 2 males) with computational modeling while manipulating predictions about content ("what") and time ("when"). We found that "when" predictions increased evoked activity over motor and prefrontal regions both at early (∼180 ms) and late (430-450 ms) latencies. "What" predictability, however, increased evoked activity only over prefrontal areas late in time (420-460 ms). Beyond these dissociable influences, we found that "what" and "when" predictability interactively modulated the amplitude of early (165 ms) evoked responses in the superior temporal gyrus. We modeled the observed neural responses using biophysically realistic neural mass models, to better understand whether "what" and "when" predictions tap into similar or different neurophysiological mechanisms. Our modeling results suggest that "what" and "when" predictability rely on complementary neural processes: "what" predictions increased short-term plasticity in auditory areas, whereas "when" predictability increased synaptic gain in motor areas. Thus, content and temporal predictions engage complementary neural mechanisms in different regions, suggesting domain-specific prediction signaling along the cortical hierarchy. Encoding predictions through different mechanisms may endow the brain with the flexibility to efficiently signal different sources of predictions, weight them by their reliability, and allow for their encoding without mutual interference.SIGNIFICANCE STATEMENT Predictions of different stimulus features facilitate sensory processing. However, it is unclear whether predictions of different attributes rely on similar or different neural mechanisms. By combining invasive electrophysiological recordings of cortical activity with experimental manipulations of participants' predictions about content and time of acoustic events, we found that the two types of predictions had dissociable influences on cortical activity, both in terms of the regions involved and the timing of the observed effects. Further, our biophysical modeling analysis suggests that predictability of content and time rely on complementary neural processes: short-term plasticity in auditory areas and synaptic gain in motor areas, respectively. This suggests that predictions of different features are encoded with complementary neural mechanisms in different brain regions.

Generative models for clinical applications in computational psychiatry

Despite the success of modern neuroimaging techniques in furthering our understanding of cognitive and pathophysiological processes, translation of these advances into clinically relevant tools has been virtually absent until now. Neuromodeling represents a powerful framework for overcoming this translational deadlock, and the development of computational models to solve clinical problems has become a major scientific goal over the last decade, as reflected by the emergence of clinically oriented neuromodeling fields like Computational Psychiatry, Computational Neurology, and Computational Psychosomatics. Generative models of brain physiology and connectivity in the human brain play a key role in this endeavor, striving for computational assays that can be applied to neuroimaging data from individual patients for differential diagnosis and treatment prediction. In this review, we focus on dynamic causal modeling (DCM) and its use for Computational Psychiatry. DCM is a widely used generative modeling framework for functional magnetic resonance imaging (fMRI) and magneto-/electroencephalography (M/EEG) data. This article reviews the basic concepts of DCM, revisits examples where it has proven valuable for addressing clinically relevant questions, and critically discusses methodological challenges and recent methodological advances. We conclude this review with a more general discussion of the promises and pitfalls of generative models in Computational Psychiatry and highlight the path that lies ahead of us. This article is categorized under: Neuroscience > Computation Neuroscience > Clinical Neuroscience.

A model of individualized canonical microcircuits supporting cognitive operations

Major cognitive functions such as language, memory, and decision-making are thought to rely on distributed networks of a large number of basic elements, called canonical microcircuits. In this theoretical study we propose a novel canonical microcircuit model and find that it supports two basic computational operations: a gating mechanism and working memory. By means of bifurcation analysis we systematically investigate the dynamical behavior of the canonical microcircuit with respect to parameters that govern the local network balance, that is, the relationship between excitation and inhibition, and key intrinsic feedback architectures of canonical microcircuits. We relate the local behavior of the canonical microcircuit to cognitive processing and demonstrate how a network of interacting canonical microcircuits enables the establishment of spatiotemporal sequences in the context of syntax parsing during sentence comprehension. This study provides a framework for using individualized canonical microcircuits for the construction of biologically realistic networks supporting cognitive operations.

Neurophysiologically-informed markers of individual variability and pharmacological manipulation of human cortical gamma

The ability to quantify synaptic function at the level of cortical microcircuits from non-invasive data would be enormously useful in the study of neuronal processing in humans and the pathophysiology that attends many neuropsychiatric disorders. Here, we provide proof of principle that one can estimate inter-and intra-laminar interactions among specific neuronal populations using induced gamma responses in the visual cortex of human subjects - using dynamic causal modelling based upon the canonical microcircuit (CMC; a simplistic model of a cortical column). Using variability in induced (spectral) responses over a large cohort of normal subjects, we find that the predominant determinants of gamma responses rest on recurrent and intrinsic connections between superficial pyramidal cells and inhibitory interneurons. Furthermore, variations in beta responses were mediated by inter-subject differences in the intrinsic connections between deep pyramidal cells and inhibitory interneurons. Interestingly, we also show that increasing the self-inhibition of superficial pyramidal cells suppresses the amplitude of gamma activity, while increasing its peak frequency. This systematic and nonlinear relationship was only disclosed by modelling the causes of induced responses. Crucially, we were able to validate this form of neurophysiological phenotyping by showing a selective effect of the GABA re-uptake inhibitor tiagabine on the rate constants of inhibitory interneurons. Remarkably, we were able to recover the pharmacodynamics of this effect over the course of several hours on a per subject basis. These findings speak to the possibility of measuring population specific synaptic function - and its response to pharmacological intervention - to provide subject-specific biomarkers of mesoscopic neuronal processes using non-invasive data. Finally, our results demonstrate that, using the CMC as a proxy, the synaptic mechanisms that underlie the gain control of neuronal message passing within and between different levels of cortical hierarchies may now be amenable to quantitative study using non-invasive (MEG) procedures.

The Cumulative Effects of Predictability on Synaptic Gain in the Auditory Processing Stream

Stimulus predictability can lead to substantial modulations of brain activity, such as shifts in sustained magnetic field amplitude, measured with magnetoencephalography (MEG). Here, we provide a mechanistic explanation of these effects using MEG data acquired from healthy human volunteers (N = 13, 7 female). In a source-level analysis of induced responses, we established the effects of orthogonal predictability manipulations of rapid tone-pip sequences (namely, sequence regularity and alphabet size) along the auditory processing stream. In auditory cortex, regular sequences with smaller alphabets induced greater gamma activity. Furthermore, sequence regularity shifted induced activity in frontal regions toward higher frequencies. To model these effects in terms of the underlying neurophysiology, we used dynamic causal modeling for cross-spectral density and estimated slow fluctuations in neural (postsynaptic) gain. Using the model-based parameters, we accurately explain the sensor-level sustained field amplitude, demonstrating that slow changes in synaptic efficacy, combined with sustained sensory input, can result in profound and sustained effects on neural responses to predictable sensory streams.

SIGNIFICANCE STATEMENT Brain activity can be strongly modulated by the predictability of stimuli it is currently processing. An example of such a modulation is a shift in sustained magnetic field amplitude, measured with magnetoencephalography. Here, we provide a mechanistic explanation of these effects. First, we establish the oscillatory neural correlates of independent predictability manipulations in hierarchically distinct areas of the auditory processing stream. Next, we use a biophysically realistic computational model to explain these effects in terms of the underlying neurophysiology. Finally, using the model-based parameters describing neural gain modulation, we can explain the previously unexplained effects observed at the sensor level. This demonstrates that slow modulations of synaptic gain can result in profound and sustained effects on neural activity.

Abnormal frontoparietal synaptic gain mediating the P300 in patients with psychotic disorder and their unaffected relatives

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. Hum Brain Mapp 38:3262-3276, 2017. © 2017 Wiley Periodicals, Inc.

Expectation violation and attention to pain jointly modulate neural gain in somatosensory cortex

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.

Linking canonical microcircuits and neuronal activity: Dynamic causal modelling of laminar recordings

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.

Cortical Dynamic Causality Network for Auditory-Motor Tasks

Motor preparation and execution require the interactions of a large-scale brain network, while the study of the dynamic changes of their interactions could uncover the underlying neural mechanism of the corresponding information processing. This dynamic analysis requires high temporal resolution of the recorded signals. Electroencephalogram (EEG) with high temporal resolution has been widely used in related studies. However, studies based on scalp EEG always lead to distorted results, due to scalp volume conduction, compared with that of cortically recorded signals. In the current study, the dynamic networks of motor preparation and execution are investigated using Go/No-go tasks performed with the left/right hand. In the analysis, the EEG source localization and dynamic causal model are combined together to investigate the neural processes of motor preparation and execution. The results show that similar network patterns with nodes distributed in the bilateral occipital lobe, bilateral temporal lobe, bilateral dorsolateral prefrontal cortex, and contralateral supplementary motor area could be revealed for both the Go and No-go tasks. Statistical testing further indicates that stronger couplings with the supplementary motor area could be found in Go and right-hand response tasks compared with No-go and left-hand response tasks, respectively. The findings in the current study demonstrate that the information exchange within the motor related brain networks plays an important role for motor related functions, i.e., the different motor functions may have the different information exchange and processing network patterns.

Gamma Frequency and the Spatial Tuning of Primary Visual Cortex

Visual stimulation produces oscillatory gamma responses in human primary visual cortex (V1) that also relate to visual perception. We have shown previously that peak gamma frequency positively correlates with central V1 cortical surface area. We hypothesized that people with larger V1 would have smaller receptive fields and that receptive field size, not V1 area, might explain this relationship. Here we set out to test this hypothesis directly by investigating the relationship between fMRI estimated population receptive field (pRF) size and gamma frequency in V1. We stimulated both the near-center and periphery of the visual field using both large and small stimuli in each location and replicated our previous finding of a positive correlation between V1 surface area and peak gamma frequency. Counter to our expectation, we found that between participants V1 size (and not PRF size) accounted for most of the variability in gamma frequency. Within-participants we found that gamma frequency increased, rather than decreased, with stimulus eccentricity directly contradicting our initial hypothesis.

Neuronal interactions between mentalising and action systems during indirect request processing

Human communication relies on the ability to process linguistic structure and to map words and utterances onto our environment. Furthermore, as what we communicate is often not directly encoded in our language (e.g. in the case of irony, jokes or indirect requests), we need to extract additional cues to infer the beliefs and desires of our conversational partners. Although the functional interplay between language and the ability to mentalise has been discussed in theoretical accounts in the past, the neurobiological underpinnings of these dynamics are currently not well understood. Here, we address this issue using functional imaging (fMRI). Participants listened to question-reply dialogues. In these dialogues, a reply is interpreted as a direct reply, an indirect reply or a request for action, depending on the question. We show that inferring meaning from indirect replies engages parts of the mentalising network (mPFC) while requests for action also activate the cortical motor system (IPL). Subsequent connectivity analysis using Dynamic Causal Modelling (DCM) revealed that this pattern of activation is best explained by an increase in effective connectivity from the mentalising network (mPFC) to the action system (IPL). These results are an important step towards a more integrative understanding of the neurobiological basis of indirect speech processing.

Probabilistic delay differential equation modeling of event-related potentials

"Dynamic causal models" (DCMs) are a promising approach in the analysis of functional neuroimaging data due to their biophysical interpretability and their consolidation of functional-segregative and functional-integrative propositions. In this theoretical note we are concerned with the DCM framework for electroencephalographically recorded event-related potentials (ERP-DCM). Intuitively, ERP-DCM combines deterministic dynamical neural mass models with dipole-based EEG forward models to describe the event-related scalp potential time-series over the entire electrode space. Since its inception, ERP-DCM has been successfully employed to capture the neural underpinnings of a wide range of neurocognitive phenomena. However, in spite of its empirical popularity, the technical literature on ERP-DCM remains somewhat patchy. A number of previous communications have detailed certain aspects of the approach, but no unified and coherent documentation exists. With this technical note, we aim to close this gap and to increase the technical accessibility of ERP-DCM. Specifically, this note makes the following novel contributions: firstly, we provide a unified and coherent review of the mathematical machinery of the latent and forward models constituting ERP-DCM by formulating the approach as a probabilistic latent delay differential equation model. Secondly, we emphasize the probabilistic nature of the model and its variational Bayesian inversion scheme by explicitly deriving the variational free energy function in terms of both the likelihood expectation and variance parameters. Thirdly, we detail and validate the estimation of the model with a special focus on the explicit form of the variational free energy function and introduce a conventional nonlinear optimization scheme for its maximization. Finally, we identify and discuss a number of computational issues which may be addressed in the future development of the approach.

Impaired prefrontal synaptic gain in people with psychosis and their relatives during the mismatch negativity

The mismatch negativity (MMN) evoked potential, a preattentive brain response to a discriminable change in auditory stimulation, is significantly reduced in psychosis. Glutamatergic theories of psychosis propose that hypofunction of NMDA receptors (on pyramidal cells and inhibitory interneurons) causes a loss of synaptic gain control. We measured changes in neuronal effective connectivity underlying the MMN using dynamic causal modeling (DCM), where the gain (excitability) of superficial pyramidal cells is explicitly parameterised. EEG data were obtained during a MMN task—for 24 patients with psychosis, 25 of their first‐degree unaffected relatives, and 35 controls—and DCM was used to estimate the excitability (modeled as self‐inhibition) of (source‐specific) superficial pyramidal populations. The MMN sources, based on previous research, included primary and secondary auditory cortices, and the right inferior frontal gyrus. Both patients with psychosis and unaffected relatives (to a lesser degree) showed increased excitability in right inferior frontal gyrus across task conditions, compared to controls. Furthermore, in the same region, both patients and their relatives showed a reversal of the normal response to deviant stimuli; that is, a decrease in excitability in comparison to standard conditions. Our results suggest that psychosis and genetic risk for the illness are associated with both context‐dependent (condition‐specific) and context‐independent abnormalities of the excitability of superficial pyramidal cell populations in the MMN paradigm. These abnormalities could relate to NMDA receptor hypofunction on both pyramidal cells and inhibitory interneurons, and appear to be linked to the genetic aetiology of the illness, thereby constituting potential endophenotypes for psychosis. Hum Brain Mapp 37:351–365, 2016. © 2015 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

Neural field simulator: two-dimensional spatio-temporal dynamics involving finite transmission speed

Neural Field models (NFM) play an important role in the understanding of neural population dynamics on a mesoscopic spatial and temporal scale. Their numerical simulation is an essential element in the analysis of their spatio-temporal dynamics. The simulation tool described in this work considers scalar spatially homogeneous neural fields taking into account a finite axonal transmission speed and synaptic temporal derivatives of first and second order. A text-based interface offers complete control of field parameters and several approaches are used to accelerate simulations. A graphical output utilizes video hardware acceleration to display running output with reduced computational hindrance compared to simulators that are exclusively software-based. Diverse applications of the tool demonstrate breather oscillations, static and dynamic Turing patterns and activity spreading with finite propagation speed. The simulator is open source to allow tailoring of code and this is presented with an extension use case.

The enhanced information flow from visual cortex to frontal area facilitates SSVEP response: evidence from model-driven and data-driven causality analysis

The neural mechanism of steady-state visual evoked potentials (SSVEP) is still not clearly understood. Especially, only certain frequency stimuli can evoke SSVEP. Our previous network study reveals that 8 Hz stimulus that can evoke strong SSVEP response shows the enhanced linkage strength between frontal and visual cortex. To further probe the directed information flow between the two cortex areas for various frequency stimuli, this paper develops a causality analysis based on the inversion of double columns model using particle swarm optimization (PSO) to characterize the directed information flow between visual and frontal cortices with the intracranial rat electroencephalograph (EEG). The estimated model parameters demonstrate that the 8 Hz stimulus shows the enhanced directional information flow from visual cortex to frontal lobe facilitates SSVEP response, which may account for the strong SSVEP response for 8 Hz stimulus. Furthermore, the similar finding is replicated by data-driven causality analysis. The inversion of neural mass model proposed in this study may be helpful to provide the new causality analysis to link the physiological model and the observed datasets in neuroscience and clinical researches.

Model-based estimation of intra-cortical connectivity using electrophysiological data

This paper provides a new method for model-based estimation of intra-cortical connectivity from electrophysiological measurements. A novel closed-form solution for the connectivity function of the Amari neural field equations is derived as a function of electrophysiological observations. The resultant intra-cortical connectivity estimate is driven from experimental data, but constrained by the mesoscopic neurodynamics that are encoded in the computational model. A demonstration is provided to show how the method can be used to image physiological mechanisms that govern cortical dynamics, which are normally hidden in clinical data from epilepsy patients. Accurate estimation performance is demonstrated using synthetic data. Following the computational testing, results from patient data are obtained that indicate a dominant increase in surround inhibition prior to seizure onset that subsides in the cases when the seizures spread.

Attentional Enhancement of Auditory Mismatch Responses: a DCM/MEG Study

Despite similar behavioral effects, attention and expectation influence evoked responses differently: Attention typically enhances event-related responses, whereas expectation reduces them. This dissociation has been reconciled under predictive coding, where prediction errors are weighted by precision associated with attentional modulation. Here, we tested the predictive coding account of attention and expectation using magnetoencephalography and modeling. Temporal attention and sensory expectation were orthogonally manipulated in an auditory mismatch paradigm, revealing opposing effects on evoked response amplitude. Mismatch negativity (MMN) was enhanced by attention, speaking against its supposedly pre-attentive nature. This interaction effect was modeled in a canonical microcircuit using dynamic causal modeling, comparing models with modulation of extrinsic and intrinsic connectivity at different levels of the auditory hierarchy. While MMN was explained by recursive interplay of sensory predictions and prediction errors, attention was linked to the gain of inhibitory interneurons, consistent with its modulation of sensory precision.

LFP and oscillations-what do they tell us?

•

A brief treatment of dynamic coordination in terms of predictive coding.

•

Understanding synchronous message passing in terms of hierarchical predictive coding.

•

Characterising cortical gain control with the dynamic causal modelling of neural fields.

•

Characterising pathophysiological oscillations with dynamic causal modelling of neural masses.

This review surveys recent trends in the use of local field potentials—and their non-invasive counterparts—to address the principles of functional brain architectures. In particular, we treat oscillations as the (observable) signature of context-sensitive changes in synaptic efficacy that underlie coordinated dynamics and message-passing in the brain. This rich source of information is now being exploited by various procedures—like dynamic causal modelling—to test hypotheses about neuronal circuits in health and disease. Furthermore, the roles played by neuromodulatory mechanisms can be addressed directly through their effects on oscillatory phenomena. These neuromodulatory or gain control processes are central to many theories of normal brain function (e.g. attention) and the pathophysiology of several neuropsychiatric conditions (e.g. Parkinson's disease).

Reversal of cortical information flow during visual imagery as compared to visual perception

The role of bottom-up and top-down connections during visual perception and the formation of mental images was examined by analyzing high-density EEG recordings of brain activity using two state-of-the-art methods for assessing the directionality of cortical signal flow: state-space Granger causality and dynamic causal modeling. We quantified the directionality of signal flow in an occipito-parieto-frontal cortical network during perception of movie clips versus mental replay of the movies and free visual imagery. Both Granger causality and dynamic causal modeling analyses revealed an increased top-down signal flow in parieto-occipital cortices during mental imagery as compared to visual perception. These results are the first direct demonstration of a reversal of the predominant direction of cortical signal flow during mental imagery as compared to perception.

Competitive interactions in sensorimotor cortex: oscillations express separation between alternative movement targets

Choice behavior is influenced by factors such as reward and number of alternatives but also by physical context, for instance, the relative position of alternative movement targets. At small separation, speeded eye or hand movements are more likely to land between targets (spatial averaging) than at larger separation. Neurocomputational models explain such behavior in terms of cortical activity being preshaped by the movement environment. Here, we manipulate target separation, as a determinant of motor cortical activity in choice behavior, to address neural mechanisms of response selection. Specifically, we investigate whether context-induced changes in the balance of cooperative and competitive interactions between competing groups of neurons are expressed in the power spectrum of sensorimotor rhythms. We recorded magnetoencephalography while participants were precued to two possible movement target locations at different angles of separation (30, 60, or 90°). After a delay, one of the locations was cued as the target for a joystick pointing movement. We found that late delay-period movement-preparatory activity increased more strongly for alternative targets at 30 than at 60 or 90° of separation. This nonlinear pattern was evident in slow event-related fields as well as in beta- and low-gamma-band suppression. A comparable pattern was found within an earlier window for theta-band synchronization. We interpret the late delay effects in terms of increased movement-preparatory activity when there is greater overlap and hence less competition between groups of neurons encoding two response alternatives. Early delay-period theta-band synchronization may reflect covert response activation relevant to behavioral spatial averaging effects.

Contrast gain control and horizontal interactions in V1: a DCM study

Using high-density electrocorticographic recordings – from awake-behaving monkeys – and dynamic causal modelling, we characterised contrast dependent gain control in visual cortex, in terms of synaptic rate constants and intrinsic connectivity. Specifically, we used neural field models to quantify the balance of excitatory and inhibitory influences; both in terms of the strength and spatial dispersion of horizontal intrinsic connections. Our results allow us to infer that increasing contrast increases the sensitivity or gain of superficial pyramidal cells to inputs from spiny stellate populations. Furthermore, changes in the effective spatial extent of horizontal coupling nuance the spatiotemporal filtering properties of cortical laminae in V1 — effectively preserving higher spatial frequencies. These results are consistent with recent non-invasive human studies of contrast dependent changes in the gain of pyramidal cells elaborating forward connections — studies designed to test specific hypotheses about precision and gain control based on predictive coding. Furthermore, they are consistent with established results showing that the receptive fields of V1 units shrink with increasing visual contrast.

•

A new observation model suitable for ECoG recordings.

•

An canonical microcircuit field model.

•

A DCM treatment of multiple experimental conditions and trial-specific effects.

•

Excitation-inhibition balance in terms of strength and dispersion.

Neural masses and fields in dynamic causal modeling

Dynamic causal modeling (DCM) provides a framework for the analysis of effective connectivity among neuronal subpopulations that subtend invasive (electrocorticograms and local field potentials) and non-invasive (electroencephalography and magnetoencephalography) electrophysiological responses. This paper reviews the suite of neuronal population models including neural masses, fields and conductance-based models that are used in DCM. These models are expressed in terms of sets of differential equations that allow one to model the synaptic underpinnings of connectivity. We describe early developments using neural mass models, where convolution-based dynamics are used to generate responses in laminar-specific populations of excitatory and inhibitory cells. We show that these models, though resting on only two simple transforms, can recapitulate the characteristics of both evoked and spectral responses observed empirically. Using an identical neuronal architecture, we show that a set of conductance based models-that consider the dynamics of specific ion-channels-present a richer space of responses; owing to non-linear interactions between conductances and membrane potentials. We propose that conductance-based models may be more appropriate when spectra present with multiple resonances. Finally, we outline a third class of models, where each neuronal subpopulation is treated as a field; in other words, as a manifold on the cortical surface. By explicitly accounting for the spatial propagation of cortical activity through partial differential equations (PDEs), we show that the topology of connectivity-through local lateral interactions among cortical layers-may be inferred, even in the absence of spatially resolved data. We also show that these models allow for a detailed analysis of structure-function relationships in the cortex. Our review highlights the relationship among these models and how the hypothesis asked of empirical data suggests an appropriate model class.

Top-down influences on local networks: basic theory with experimental implications

The response of a population of cortical neurons to an external stimulus depends not only on the receptive field properties of the neurons, but also the level of arousal and attention or goal-oriented cognitive biases that guide information processing. These top-down effects on cortical neurons bias the output of the neurons and affect behavioral outcomes such as stimulus detection, discrimination, and response time. In any physiological study, neural dynamics are observed in a specific brain state; the background state partly determines neuronal excitability. Experimental studies in humans and animal models have also demonstrated that slow oscillations (typically in the alpha or theta bands) modulate the fast oscillations (gamma band) associated with local networks of neurons. Cross-frequency interaction is of interest as a mechanism for top-down or bottom up interactions between systems at different spatial scales. We develop a generic model of top-down influences on local networks appropriate for comparison with EEG. EEG provides excellent temporal resolution to investigate neuronal oscillations but is space-averaged on the cm scale. Thus, appropriate EEG models are developed in terms of population synaptic activity. We used the Wilson–Cowan population model to investigate fast (gamma band) oscillations generated by a local network of excitatory and inhibitory neurons. We modified the Wilson–Cowan equations to make them more physiologically realistic by explicitly incorporating background state variables into the model. We found that the population response is strongly influenced by the background state. We apply the model to reproduce the modulation of gamma rhythms by theta rhythms as has been observed in animal models and human ECoG and EEG studies. The concept of a dynamic background state presented here using the Wilson–Cowan model can be readily applied to incorporate top-down modulation in more detailed models of specific cortical systems.

High-frequency brain activity and muscle artifacts in MEG/EEG: a review and recommendations

In recent years high-frequency brain activity in the gamma-frequency band (30-80 Hz) and above has become the focus of a growing body of work in MEG/EEG research. Unfortunately, high-frequency neural activity overlaps entirely with the spectral bandwidth of muscle activity (~20-300 Hz). It is becoming appreciated that artifacts of muscle activity may contaminate a number of non-invasive reports of high-frequency activity. In this review, the spectral, spatial, and temporal characteristics of muscle artifacts are compared with those described (so far) for high-frequency neural activity. In addition, several of the techniques that are being developed to help suppress muscle artifacts in MEG/EEG are reviewed. Suggestions are made for the collection, analysis, and presentation of experimental data with the aim of reducing the number of publications in the future that may contain muscle artifacts.

Enhanced stimulus-induced gamma activity in humans during propofol-induced sedation

Stimulus-induced gamma oscillations in the 30-80 Hz range have been implicated in a wide number of functions including visual processing, memory and attention. While occipital gamma-band oscillations can be pharmacologically modified in animal preparations, pharmacological modulation of stimulus-induced visual gamma oscillations has yet to be demonstrated in non-invasive human recordings. Here, in fifteen healthy humans volunteers, we probed the effects of the GABAA agonist and sedative propofol on stimulus-related gamma activity recorded with magnetoencephalography, using a simple visual grating stimulus designed to elicit gamma oscillations in the primary visual cortex. During propofol sedation as compared to the normal awake state, a significant 60% increase in stimulus-induced gamma amplitude was seen together with a 94% enhancement of stimulus-induced alpha suppression and a simultaneous reduction in the amplitude of the pattern-onset evoked response. These data demonstrate, that propofol-induced sedation is accompanied by increased stimulus-induced gamma activity providing a potential window into mechanisms of gamma-oscillation generation in humans.

From oscillatory transcranial current stimulation to scalp EEG changes: a biophysical and physiological modeling study

Both biophysical and neurophysiological aspects need to be considered to assess the impact of electric fields induced by transcranial current stimulation (tCS) on the cerebral cortex and the subsequent effects occurring on scalp EEG. The objective of this work was to elaborate a global model allowing for the simulation of scalp EEG signals under tCS. In our integrated modeling approach, realistic meshes of the head tissues and of the stimulation electrodes were first built to map the generated electric field distribution on the cortical surface. Secondly, source activities at various cortical macro-regions were generated by means of a computational model of neuronal populations. The model parameters were adjusted so that populations generated an oscillating activity around 10 Hz resembling typical EEG alpha activity. In order to account for tCS effects and following current biophysical models, the calculated component of the electric field normal to the cortex was used to locally influence the activity of neuronal populations. Lastly, EEG under both spontaneous and tACS-stimulated (transcranial sinunoidal tCS from 4 to 16 Hz) brain activity was simulated at the level of scalp electrodes by solving the forward problem in the aforementioned realistic head model. Under the 10 Hz-tACS condition, a significant increase in alpha power occurred in simulated scalp EEG signals as compared to the no-stimulation condition. This increase involved most channels bilaterally, was more pronounced on posterior electrodes and was only significant for tACS frequencies from 8 to 12 Hz. The immediate effects of tACS in the model agreed with the post-tACS results previously reported in real subjects. Moreover, additional information was also brought by the model at other electrode positions or stimulation frequency. This suggests that our modeling approach can be used to compare, interpret and predict changes occurring on EEG with respect to parameters used in specific stimulation configurations.