Testing non-linearity and directedness of interactions between neural groups in the macaque inferotemporal cortex

Connectivity Analysis in EEG Data: A Tutorial Review of the State of the Art and Emerging Trends

Understanding how different areas of the human brain communicate with each other is a crucial issue in neuroscience. The concepts of structural, functional and effective connectivity have been widely exploited to describe the human connectome, consisting of brain networks, their structural connections and functional interactions. Despite high-spatial-resolution imaging techniques such as functional magnetic resonance imaging (fMRI) being widely used to map this complex network of multiple interactions, electroencephalographic (EEG) recordings claim high temporal resolution and are thus perfectly suitable to describe either spatially distributed and temporally dynamic patterns of neural activation and connectivity. In this work, we provide a technical account and a categorization of the most-used data-driven approaches to assess brain-functional connectivity, intended as the study of the statistical dependencies between the recorded EEG signals. Different pairwise and multivariate, as well as directed and non-directed connectivity metrics are discussed with a pros-cons approach, in the time, frequency, and information-theoretic domains. The establishment of conceptual and mathematical relationships between metrics from these three frameworks, and the discussion of novel methodological approaches, will allow the reader to go deep into the problem of inferring functional connectivity in complex networks. Furthermore, emerging trends for the description of extended forms of connectivity (e.g., high-order interactions) are also discussed, along with graph-theory tools exploring the topological properties of the network of connections provided by the proposed metrics. Applications to EEG data are reviewed. In addition, the importance of source localization, and the impacts of signal acquisition and pre-processing techniques (e.g., filtering, source localization, and artifact rejection) on the connectivity estimates are recognized and discussed. By going through this review, the reader could delve deeply into the entire process of EEG pre-processing and analysis for the study of brain functional connectivity and learning, thereby exploiting novel methodologies and approaches to the problem of inferring connectivity within complex networks.

Parallel factor analysis for multidimensional decomposition of functional near-infrared spectroscopy data

SIGNIFICANCE

Current techniques for data analysis in functional near-infrared spectroscopy (fNIRS), such as artifact correction, do not allow to integrate the information originating from both wavelengths, considering only temporal and spatial dimensions of the signal's structure. Parallel factor analysis (PARAFAC) has previously been validated as a multidimensional decomposition technique in other neuroimaging fields.

AIM

We aimed to introduce and validate the use of PARAFAC for the analysis of fNIRS data, which is inherently multidimensional (time, space, and wavelength).

APPROACH

We used data acquired in 17 healthy adults during a verbal fluency task to compare the efficacy of PARAFAC for motion artifact correction to traditional two-dimensional decomposition techniques, i.e., target principal (tPCA) and independent component analysis (ICA). Correction performance was further evaluated under controlled conditions with simulated artifacts and hemodynamic response functions.

RESULTS

PARAFAC achieved significantly higher improvement in data quality as compared to tPCA and ICA. Correction in several simulated signals further validated its use and promoted it as a robust method independent of the artifact's characteristics.

CONCLUSIONS

This study describes the first implementation of PARAFAC in fNIRS and provides validation for its use to correct artifacts. PARAFAC is a promising data-driven alternative for multidimensional data analyses in fNIRS and this study paves the way for further applications.

Phase-based causality analysis with partial mutual information from mixed embedding

Instantaneous phases extracted from multivariate time series can retain information about the relationships between the underlying mechanisms that generate the series. Although phases have been widely used in the study of nondirectional coupling and connectivity, they have not found similar appeal in the study of causality. Herein, we present a new method for phase-based causality analysis, which combines ideas from the mixed embedding technique and the information-theoretic approach to causality in coupled oscillatory systems. We then use the introduced method to investigate causality in simulated datasets of bivariate, unidirectionally paired systems from combinations of Rössler, Lorenz, van der Pol, and Mackey-Glass equations. We observe that causality analysis using the phases can capture the true causal relation for coupling strength smaller than the analysis based on the amplitudes can capture. On the other hand, the causality estimation based on the phases tends to have larger variability, which is attributed more to the phase extraction process than the actual phase-based causality method. In addition, an application on real electroencephalographic data from an experiment on elicited human emotional states reinforces the usefulness of phases in causality identification.

Multivariate autoregressive model estimation for high dimensional intracranial electrophysiological data

Fundamental to elucidating the functional organization of the brain is the assessment of causal interactions between different brain regions. Multivariate autoregressive (MVAR) modeling techniques applied to multisite electrophysiological recordings are a promising avenue for identifying such causal links. They estimate the degree to which past activity in one or more brain regions is predictive of another region's present activity, while simultaneously accounting for the mediating effects of other regions. Including as many mediating variables as possible in the model has the benefit of drastically reducing the odds of detecting spurious causal connectivity. However, effective bounds on the number of MVAR model coefficients that can be estimated reliably from limited data make exploiting the potential of MVAR models challenging for even modest numbers of recording sites. Here, we utilize well-established dimensionality-reduction techniques to fit MVAR models to human intracranial data from ∼100 - 200 recording sites spanning dozens of anatomically and functionally distinct cortical regions. First, we show that high dimensional MVAR models can be successfully estimated from long segments of data and yield plausible connectivity profiles. Next, we use these models to generate synthetic data with known ground-truth connectivity to explore the utility of applying principal component analysis and group least absolute shrinkage and selection operator (gLASSO) to reduce the number of parameters (connections) during model fitting to shorter data segments. We show that gLASSO is highly effective for recovering ground-truth connectivity in the limited data regime, capturing important features of connectivity for high-dimensional models with as little as 10 seconds of data. The methods presented here have broad applicability to the analysis of high-dimensional time series data in neuroscience, facilitating the elucidation of the neural basis of sensation, cognition, and arousal.

Effects of Acute Stress on the Oscillatory Activity of the Hippocampus-Amygdala-Prefrontal Cortex Network

Displaying a stress response to threatening stimuli is essential for survival. These reactions must be adjusted to be adaptive. Otherwise, even mental illnesses may develop. Describing the physiological stress response may contribute to distinguishing the abnormal responses that accompany the pathology, which may help to improve the development of both diagnoses and treatments. Recent advances have elucidated many of the processes and structures involved in stress response management; however, there is still much to unravel regarding this phenomenon. The main aim of the present research is to characterize the response of three brain areas deeply involved in the stress response (i.e., to an acute stressful experience). Specifically, the electrophysiological activity of the infralimbic division of the medial prefrontal cortex (IL), the basolateral nucleus of the amygdala (BLA), and the dorsal hippocampus (dHPC) was recorded after the infusion of 0.5 µl of corticosterone-releasing factor into the dorsal raphe nucleus (DRN), a procedure which has been validated as a paradigm to cause acute stress. This procedure induced a delayed reduction in slow waves in the three structures, and an increase in faster oscillations, such as those in theta, beta, and gamma bands. The mutual information at low theta frequencies between the BLA and the IL increased, and the delta and slow wave mutual information decreased. The low theta-mid gamma phase-amplitude coupling increased within BLA, as well as between BLA and IL. This electrical pattern may facilitate the activation of these structures, in response to the stressor, and memory consolidation.

Stimulus evoked causality estimation in stereo-EEG

Objective.Stereo-electroencephalography (SEEG) has recently gained importance in analyzing brain functions. Its high temporal resolution and spatial specificity make it a powerful tool to investigate the strength, direction, and spectral content of brain networks interactions, especially when these connections are stimulus-evoked. However, choosing the best approach to evaluate the flow of information is not trivial, due to the lack of validated methods explicitly designed for SEEG.Approach.We propose a novel non-parametric statistical test for event-related causality (ERC) assessment on SEEG recordings. Here, we refer to the ERC as the causality evoked by a particular part of the stimulus (a response window (RW)). We also present a data surrogation method to evaluate the performance of a causality estimation algorithm. We finally validated our pipeline using surrogate SEEG data derived from an experimentally collected dataset, and compared the most used and successful measures to estimate effective connectivity, belonging to the Geweke-Granger causality framework.Main results.Here we show that our workflow correctly identified all the directed connections in the RW of the surrogate data and prove the robustness of the procedure against synthetic noise with amplitude exceeding physiological-plausible values. Among the causality measures tested, partial directed coherence performed best.Significance.This is the first non-parametric statistical test for ERC estimation explicitly designed for SEEG datasets. The pipeline, in principle, can also be applied to the analysis of any type of time-varying estimator, if there exists a clearly defined RW.

Directed Connectivity Analysis of the Brain Network in Mathematically Gifted Adolescents

The neurocognitive characteristics of mathematically gifted adolescents are characterized by highly developed functional interactions between the right hemisphere and excellent cognitive control of the prefrontal cortex, enhanced frontoparietal cortex, and posterior parietal cortex. However, it is still unclear when and how these cortical interactions occur. In this paper, we used directional coherence analysis based on Granger causality to study the interactions between the frontal brain area and the posterior brain area in the mathematical frontoparietal network system during deductive reasoning tasks. Specifically, the scalp electroencephalography (EEG) signal was first converted into a cortical dipole source signal to construct a Granger causality network over the θ-band and γ-band ranges. We constructed the binary Granger causality network at the 40 pairs of cortical nodes in the frontal lobe and parietal lobe across the θ-band and the γ-band, which were selected as regions of interest (ROI). We then used graph theory to analyze the network differences. It was found that, in the process of reasoning tasks, the frontoparietal regions of the mathematically gifted show stronger working memory information processing at the θ-band. Additionally, in the middle and late stages of the conclusion period, the mathematically talented individuals have less information flow in the anterior and posterior parietal regions of the brain than the normal subjects. We draw the conclusion that the mathematically gifted brain frontoparietal network appears to have more “automated” information processing during reasoning tasks.

Electrophysiological Brain Connectivity: Theory and Implementation

We review the theory and algorithms of electrophysiological brain connectivity analysis. This tutorial is aimed at providing an introduction to brain functional connectivity from electrophysiological signals, including electroencephalography (EEG), magnetoencephalography (MEG), electrocorticography (ECoG), stereoelectroencephalography (SEEG). Various connectivity estimators are discussed, and algorithms introduced. Important issues for estimating and mapping brain functional connectivity with electrophysiology are discussed.

A Time-Varying Connectivity Analysis from Distributed EEG Sources: A Simulation Study

Time-varying connectivity analysis based on sources reconstructed using inverse modeling of electroencephalographic (EEG) data is important to understand the dynamic behaviour of the brain. We simulated cortical data from a visual spatial attention network with a time-varying connectivity structure, and then simulated the propagation to the scalp to obtain EEG data. Distributed EEG source modeling using sLORETA was applied. We compared different dipole (representing a source) selection strategies based on their time series in a region of interest. Next, we estimated multivariate autoregressive (MVAR) parameters using classical Kalman filter and general linear Kalman filter approaches followed by the calculation of partial directed coherence (PDC). MVAR parameters and PDC values for the selected sources were compared with the ground-truth. We found that the best strategy to extract the time series of a region of interest was to select a dipole with time series showing the highest correlation with the average time series in the region of interest. Dipole selection based on power or based on the largest singular value offer comparable alternatives. Among the different Kalman filter approaches, the use of a general linear Kalman filter was preferred to estimate PDC based connectivity except when only a small number of trials are available. In the latter case, the classical Kalman filter can be an alternative.

Causal relationships between neurons of the nucleus incertus and the hippocampal theta activity in the rat

KEY POINTS

The nucleus incertus is a key node of the brainstem circuitry involved in hippocampal theta rhythmicity. Synchronisation exists between the nucleus incertus and hippocampal activities during theta periods. By the Granger causality analysis, we demonstrated a directional information flow between theta rhythmical neurons in the nucleus incertus and the hippocampus in theta-on states. The electrical stimulation of the nucleus incertus is also able to evoke a phase reset of the hippocampal theta wave. Our data suggest that the nucleus incertus is a key node of theta generation and the modulation network.

ABSTRACT

In recent years, a body of evidence has shown that the nucleus incertus (NI), in the dorsal tegmental pons, is a key node of the brainstem circuitry involved in hippocampal theta rhythmicity. Ascending reticular brainstem system activation evokes hippocampal theta rhythm with coupled neuronal activity in the NI. In a recent paper, we showed three populations of neurons in the NI with differential firing during hippocampal theta activation. The objective of this work was to better evaluate the causal relationship between the activity of NI neurons and the hippocampus during theta activation in order to further understand the role of the NI in the theta network. A Granger causality analysis was run to determine whether hippocampal theta activity with sensory-evoked theta depends on the neuronal activity of the NI, or vice versa. The analysis showed causal interdependence between the NI and the hippocampus during theta activity, whose directional flow depended on the different neuronal assemblies of the NI. Whereas type I and II NI neurons mainly acted as receptors of hippocampal information, type III neuronal activity was the predominant source of flow between the NI and the hippocampus in theta states. We further determined that the electrical activation of the NI was able to reset hippocampal waves with enhanced theta-band power, depending on the septal area. Collectively, these data suggest that hippocampal theta oscillations after sensory activation show dependence on NI neuron activity, which could play a key role in establishing optimal conditions for memory encoding.

Multiscale modeling of brain dynamics: from single neurons and networks to mathematical tools

The extreme complexity of the brain naturally requires mathematical modeling approaches on a large variety of scales; the spectrum ranges from single neuron dynamics over the behavior of groups of neurons to neuronal network activity. Thus, the connection between the microscopic scale (single neuron activity) to macroscopic behavior (emergent behavior of the collective dynamics) and vice versa is a key to understand the brain in its complexity. In this work, we attempt a review of a wide range of approaches, ranging from the modeling of single neuron dynamics to machine learning. The models include biophysical as well as data-driven phenomenological models. The discussed models include Hodgkin-Huxley, FitzHugh-Nagumo, coupled oscillators (Kuramoto oscillators, Rössler oscillators, and the Hindmarsh-Rose neuron), Integrate and Fire, networks of neurons, and neural field equations. In addition to the mathematical models, important mathematical methods in multiscale modeling and reconstruction of the causal connectivity are sketched. The methods include linear and nonlinear tools from statistics, data analysis, and time series analysis up to differential equations, dynamical systems, and bifurcation theory, including Granger causal connectivity analysis, phase synchronization connectivity analysis, principal component analysis (PCA), independent component analysis (ICA), and manifold learning algorithms such as ISOMAP, and diffusion maps and equation-free techniques. WIREs Syst Biol Med 2016, 8:438-458. doi: 10.1002/wsbm.1348 For further resources related to this article, please visit the WIREs website.

A Nonlinear Causality Estimator Based on Non-Parametric Multiplicative Regression

Causal prediction has become a popular tool for neuroscience applications, as it allows the study of relationships between different brain areas during rest, cognitive tasks or brain disorders. We propose a nonparametric approach for the estimation of nonlinear causal prediction for multivariate time series. In the proposed estimator, C NPMR , Autoregressive modeling is replaced by Nonparametric Multiplicative Regression (NPMR). NPMR quantifies interactions between a response variable (effect) and a set of predictor variables (cause); here, we modified NPMR for model prediction. We also demonstrate how a particular measure, the sensitivity Q, could be used to reveal the structure of the underlying causal relationships. We apply C NPMR on artificial data with known ground truth (5 datasets), as well as physiological data (2 datasets). C NPMR correctly identifies both linear and nonlinear causal connections that are present in the artificial data, as well as physiologically relevant connectivity in the real data, and does not seem to be affected by filtering. The Sensitivity measure also provides useful information about the latent connectivity.The proposed estimator addresses many of the limitations of linear Granger causality and other nonlinear causality estimators. C NPMR is compared with pairwise and conditional Granger causality (linear) and Kernel-Granger causality (nonlinear). The proposed estimator can be applied to pairwise or multivariate estimations without any modifications to the main method. Its nonpametric nature, its ability to capture nonlinear relationships and its robustness to filtering make it appealing for a number of applications.

Identification of Dominant Excitation Patterns and Sources of Atrial Fibrillation by Causality Analysis

Burden of atrial fibrillation (AF) can be reduced by ablation of sources of electrical impulses driving AF but driver identification is still challenging. This study presents a new methodology based on causality analysis that allows identifying the hierarchically dominant areas driving AF. Identification of dominant propagation patterns was achieved by computing causal relations between intracardiac multi-electrode catheter recordings of four paroxysmal AF patients during sinus rhythm, pacing and AF. In addition, realistic mathematical models of the atria during AF were used to validate the methodology both in the presence and absence of dominant frequency (DF) gradients. During electrical pacing, sources of propagation patterns detected by causality analysis were consistent with the location of the stimulating catheter. During AF, propagation patterns presented temporal variability, but a dominant direction accounted for significantly more propagations than other directions (49 ± 15% vs. 14 ± 13% or less, p < 0.01). Both in patients with a DF gradient and in mathematical models, causal maps allowed the identification of sites responsible for maintenance of AF. Causal maps allowed the identification of atrial dominant sites. In particular, causality analysis resulted in stable dominant cause-effect propagation directions during AF and could serve as a guide for performing ablation procedures in AF patients.

Foundational perspectives on causality in large-scale brain networks

A profusion of recent work in cognitive neuroscience has been concerned with the endeavor to uncover causal influences in large-scale brain networks. However, despite the fact that many papers give a nod to the important theoretical challenges posed by the concept of causality, this explosion of research has generally not been accompanied by a rigorous conceptual analysis of the nature of causality in the brain. This review provides both a descriptive and prescriptive account of the nature of causality as found within and between large-scale brain networks. In short, it seeks to clarify the concept of causality in large-scale brain networks both philosophically and scientifically. This is accomplished by briefly reviewing the rich philosophical history of work on causality, especially focusing on contributions by David Hume, Immanuel Kant, Bertrand Russell, and Christopher Hitchcock. We go on to discuss the impact that various interpretations of modern physics have had on our understanding of causality. Throughout all this, a central focus is the distinction between theories of deterministic causality (DC), whereby causes uniquely determine their effects, and probabilistic causality (PC), whereby causes change the probability of occurrence of their effects. We argue that, given the topological complexity of its large-scale connectivity, the brain should be considered as a complex system and its causal influences treated as probabilistic in nature. We conclude that PC is well suited for explaining causality in the brain for three reasons: (1) brain causality is often mutual; (2) connectional convergence dictates that only rarely is the activity of one neuronal population uniquely determined by another one; and (3) the causal influences exerted between neuronal populations may not have observable effects. A number of different techniques are currently available to characterize causal influence in the brain. Typically, these techniques quantify the statistical likelihood that a change in the activity of one neuronal population affects the activity in another. We argue that these measures access the inherently probabilistic nature of causal influences in the brain, and are thus better suited for large-scale brain network analysis than are DC-based measures. Our work is consistent with recent advances in the philosophical study of probabilistic causality, which originated from inherent conceptual problems with deterministic regularity theories. It also resonates with concepts of stochasticity that were involved in establishing modern physics. In summary, we argue that probabilistic causality is a conceptually appropriate foundation for describing neural causality in the brain.

Functional brain connectivity from EEG in epilepsy: seizure prediction and epileptogenic focus localization

Today, neuroimaging techniques are frequently used to investigate the integration of functionally specialized brain regions in a network. Functional connectivity, which quantifies the statistical dependencies among the dynamics of simultaneously recorded signals, allows to infer the dynamical interactions of segregated brain regions. In this review we discuss how the functional connectivity patterns obtained from intracranial and scalp electroencephalographic (EEG) recordings reveal information about the dynamics of the epileptic brain and can be used to predict upcoming seizures and to localize the seizure onset zone. The added value of extracting information that is not visibly identifiable in the EEG data using functional connectivity analysis is stressed. Despite the fact that many studies have showed promising results, we must conclude that functional connectivity analysis has not made its way into clinical practice yet.

Reliability of Dynamic Causal Modeling using the Statistical Parametric Mapping Toolbox

Dynamic causal modeling (DCM) is a recently developed approach for effective connectivity measurement in the brain. It has attracted considerable attention in recent years and quite widespread used to investigate brain connectivity in response to different tasks as well as auditory, visual, and somatosensory stimulation. This method uses complex algorithms, and currently the only implementation available is the Statistical Parametric Mapping (SPM8) toolbox with functionality for use on EEG and fMRI. The objective of the current work is to test the robustness of the toolbox when applied to EEG, by comparing results obtained from various versions of the software and operating systems when using identical datasets. Contrary to expectations, it was found that estimated connectivities were not consistent between different operating systems, the version of SPM8, or the version of MATLAB being used. The exact cause of this problem is not clear, but may relate to the high number of parameters in the model. Caution is thus recommended when interpreting the results of DCM estimated with the SPM8 software.

Successful reconstruction of a physiological circuit with known connectivity from spiking activity alone

Identifying the structure and dynamics of synaptic interactions between neurons is the first step to understanding neural network dynamics. The presence of synaptic connections is traditionally inferred through the use of targeted stimulation and paired recordings or by post-hoc histology. More recently, causal network inference algorithms have been proposed to deduce connectivity directly from electrophysiological signals, such as extracellularly recorded spiking activity. Usually, these algorithms have not been validated on a neurophysiological data set for which the actual circuitry is known. Recent work has shown that traditional network inference algorithms based on linear models typically fail to identify the correct coupling of a small central pattern generating circuit in the stomatogastric ganglion of the crab Cancer borealis. In this work, we show that point process models of observed spike trains can guide inference of relative connectivity estimates that match the known physiological connectivity of the central pattern generator up to a choice of threshold. We elucidate the necessary steps to derive faithful connectivity estimates from a model that incorporates the spike train nature of the data. We then apply the model to measure changes in the effective connectivity pattern in response to two pharmacological interventions, which affect both intrinsic neural dynamics and synaptic transmission. Our results provide the first successful application of a network inference algorithm to a circuit for which the actual physiological synapses between neurons are known. The point process methodology presented here generalizes well to larger networks and can describe the statistics of neural populations. In general we show that advanced statistical models allow for the characterization of effective network structure, deciphering underlying network dynamics and estimating information-processing capabilities.

To appreciate how neural circuits control behaviors, we must understand two things. First, how the neurons comprising the circuit are connected, and second, how neurons and their connections change after learning or in response to neuromodulators. Neuronal connectivity is difficult to determine experimentally, whereas neuronal activity can often be readily measured. We describe a statistical model to estimate circuit connectivity directly from measured activity patterns. We use the timing relationships between observed spikes to predict synaptic interactions between simultaneously observed neurons. The model estimate provides each predicted connection with a curve that represents how strongly, and at which temporal delays, one circuit element effectively influences another. These curves are analogous to synaptic interactions of the level of the membrane potential of biological neurons and share some of their features such as being inhibitory or excitatory. We test our method on recordings from the pyloric circuit in the crab stomatogastric ganglion, a small circuit whose connectivity is completely known beforehand, and find that the predicted circuit matches the biological one — a result other techniques failed to achieve. In addition, we show that drug manipulations impacting the circuit are revealed by this technique. These results illustrate the utility of our analysis approach for inferring connections from neural spiking activity.

Measuring information-transfer delays

In complex networks such as gene networks, traffic systems or brain circuits it is important to understand how long it takes for the different parts of the network to effectively influence one another. In the brain, for example, axonal delays between brain areas can amount to several tens of milliseconds, adding an intrinsic component to any timing-based processing of information. Inferring neural interaction delays is thus needed to interpret the information transfer revealed by any analysis of directed interactions across brain structures. However, a robust estimation of interaction delays from neural activity faces several challenges if modeling assumptions on interaction mechanisms are wrong or cannot be made. Here, we propose a robust estimator for neuronal interaction delays rooted in an information-theoretic framework, which allows a model-free exploration of interactions. In particular, we extend transfer entropy to account for delayed source-target interactions, while crucially retaining the conditioning on the embedded target state at the immediately previous time step. We prove that this particular extension is indeed guaranteed to identify interaction delays between two coupled systems and is the only relevant option in keeping with Wiener’s principle of causality. We demonstrate the performance of our approach in detecting interaction delays on finite data by numerical simulations of stochastic and deterministic processes, as well as on local field potential recordings. We also show the ability of the extended transfer entropy to detect the presence of multiple delays, as well as feedback loops. While evaluated on neuroscience data, we expect the estimator to be useful in other fields dealing with network dynamics.

Measuring connectivity in linear multivariate processes: definitions, interpretation, and practical analysis

This tutorial paper introduces a common framework for the evaluation of widely used frequency-domain measures of coupling (coherence, partial coherence) and causality (directed coherence, partial directed coherence) from the parametric representation of linear multivariate (MV) processes. After providing a comprehensive time-domain definition of the various forms of connectivity observed in MV processes, we particularize them to MV autoregressive (MVAR) processes and derive the corresponding frequency-domain measures. Then, we discuss the theoretical interpretation of these MVAR-based connectivity measures, showing that each of them reflects a specific time-domain connectivity definition and how this results in the description of peculiar aspects of the information transfer in MV processes. Furthermore, issues related to the practical utilization of these measures on real-time series are pointed out, including MVAR model estimation and significance assessment. Finally, limitations and pitfalls arising from model mis-specification are discussed, indicating possible solutions and providing practical recommendations for a safe computation of the connectivity measures. An example of estimation of the presented measures from multiple EEG signals recorded during a combined visuomotor task is also reported, showing how evaluation of coupling and causality in the frequency domain may help describing specific neurophysiological mechanisms.

Cross validation for selection of cortical interaction models from scalp EEG or MEG

A cross-validation (CV) method based on state-space framework is introduced for comparing the fidelity of different cortical interaction models to the measured scalp electroencephalogram (EEG) or magnetoencephalography (MEG) data being modeled. A state equation models the cortical interaction dynamics and an observation equation represents the scalp measurement of cortical activity and noise. The measured data are partitioned into training and test sets. The training set is used to estimate model parameters and the model quality is evaluated by computing test data innovations for the estimated model. Two CV metrics normalized mean square error and log-likelihood are estimated by averaging over different training/test partitions of the data. The effectiveness of this method of model selection is illustrated by comparing two linear modeling methods and two nonlinear modeling methods on simulated EEG data derived using both known dynamic systems and measured electrocorticography data from an epilepsy patient.

A short history of causal modeling of fMRI data

Twenty years ago, the discovery of the blood oxygen level dependent (BOLD) contrast and invention of functional magnetic resonance imaging (MRI) not only allowed for enhanced analyses of regional brain activity, but also laid the foundation for novel approaches to studying effective connectivity, which is essential for mechanistically interpretable accounts of neuronal systems. Dynamic causal modeling (DCM) and Granger causality (G-causality) modeling have since become the most frequently used techniques for inferring effective connectivity from fMRI data. In this paper, we provide a short historical overview of these approaches, describing milestones of their development from our subjective perspectives.

TRENTOOL: a Matlab open source toolbox to analyse information flow in time series data with transfer entropy

Background

Transfer entropy (TE) is a measure for the detection of directed interactions. Transfer entropy is an information theoretic implementation of Wiener's principle of observational causality. It offers an approach to the detection of neuronal interactions that is free of an explicit model of the interactions. Hence, it offers the power to analyze linear and nonlinear interactions alike. This allows for example the comprehensive analysis of directed interactions in neural networks at various levels of description. Here we present the open-source MATLAB toolbox TRENTOOL that allows the user to handle the considerable complexity of this measure and to validate the obtained results using non-parametrical statistical testing. We demonstrate the use of the toolbox and the performance of the algorithm on simulated data with nonlinear (quadratic) coupling and on local field potentials (LFP) recorded from the retina and the optic tectum of the turtle (Pseudemys scripta elegans) where a neuronal one-way connection is likely present.

Results

In simulated data TE detected information flow in the simulated direction reliably with false positives not exceeding the rates expected under the null hypothesis. In the LFP data we found directed interactions from the retina to the tectum, despite the complicated signal transformations between these stages. No false positive interactions in the reverse directions were detected.

Conclusions

TRENTOOL is an implementation of transfer entropy and mutual information analysis that aims to support the user in the application of this information theoretic measure. TRENTOOL is implemented as a MATLAB toolbox and available under an open source license (GPL v3). For the use with neural data TRENTOOL seamlessly integrates with the popular FieldTrip toolbox.

More discussions for granger causality and new causality measures

Granger causality (GC) has been widely applied in economics and neuroscience to reveal causality influence of time series. In our previous paper (Hu et al., in IEEE Trans on Neural Netw, 22(6), pp. 829-844, 2011), we proposed new causalities in time and frequency domains and particularly focused on new causality in frequency domain by pointing out the shortcomings/limitations of GC or Granger-alike causality metrics and the advantages of new causality. In this paper we continue our previous discussions and focus on new causality and GC or Granger-alike causality metrics in time domain. Although one strong motivation was introduced in our previous paper (Hu et al., in IEEE Trans on Neural Netw, 22(6), pp. 829-844, 2011) we here present additional motivation for the proposed new causality metric and restate the previous motivation for completeness. We point out one property of conditional GC in time domain and the shortcomings/limitations of conditional GC which cannot reveal the real strength of the directional causality among three time series. We also show the shortcomings/limitations of directed causality (DC) or normalize DC for multivariate time series and demonstrate it cannot reveal real causality at all. By calculating GC and new causality values for an example we demonstrate the influence of one of the time series on the other is linearly increased as the coupling strength is linearly increased. This fact further supports reasonability of new causality metric. We point out that larger instantaneous correlation does not necessarily mean larger true causality (e.g., GC and new causality), or vice versa. Finally we conduct analysis of statistical test for significance and asymptotic distribution property of new causality metric by illustrative examples.

Effective connectivity: influence, causality and biophysical modeling

This is the final paper in a Comments and Controversies series dedicated to "The identification of interacting networks in the brain using fMRI: Model selection, causality and deconvolution". We argue that discovering effective connectivity depends critically on state-space models with biophysically informed observation and state equations. These models have to be endowed with priors on unknown parameters and afford checks for model Identifiability. We consider the similarities and differences among Dynamic Causal Modeling, Granger Causal Modeling and other approaches. We establish links between past and current statistical causal modeling, in terms of Bayesian dependency graphs and Wiener-Akaike-Granger-Schweder influence measures. We show that some of the challenges faced in this field have promising solutions and speculate on future developments.

Complex propagation patterns characterize human cortical activity during slow-wave sleep

Cortical electrical activity during nonrapid eye movement (non-REM) sleep is dominated by slow-wave activity (SWA). At larger spatial scales (∼2-30 cm), investigated by scalp EEG recordings, SWA has been shown to propagate globally over wide cortical regions as traveling waves, which has been proposed to serve as a temporal framework for neural plasticity. However, whether SWA dynamics at finer spatial scales also reflects the orderly propagation has not previously been investigated in humans. To reveal the local, finer spatial scale (∼1-6 cm) patterns of SWA propagation during non-REM sleep, electrocorticographic (ECoG) recordings were conducted from subdurally implanted electrode grids and a nonlinear correlation technique [mutual information (MI)] was implemented. MI analysis revealed spatial maps of correlations between cortical areas demonstrating SWA propagation directions, speed, and association strength. Highest correlations, indicating significant coupling, were detected during the initial positive-going deflection of slow waves. SWA propagated predominantly between adjacent cortical areas, albeit spatial noncontinuities were also frequently observed. MI analysis further uncovered significant convergence and divergence patterns. Areas receiving the most convergent activity were similar to those with high divergence rate, while reciprocal and circular propagation of SWA was also frequent. We hypothesize that SWA is characterized by distinct attributes depending on the spatial scale observed. At larger spatial scales, the orderly SWA propagation dominates; at the finer scale of the ECoG recordings, non-REM sleep is characterized by complex SWA propagation patterns.

Dynamics of large-scale cortical interactions at high gamma frequencies during word production: event related causality (ERC) analysis of human electrocorticography (ECoG)

Intracranial EEG studies in humans have shown that functional brain activation in a variety of functional-anatomic domains of human cortex is associated with an increase in power at a broad range of high gamma (>60Hz) frequencies. Although these electrophysiological responses are highly specific for the location and timing of cortical processing and in animal recordings are highly correlated with increased population firing rates, there has been little direct empirical evidence for causal interactions between different recording sites at high gamma frequencies. Such causal interactions are hypothesized to occur during cognitive tasks that activate multiple brain regions. To determine whether such causal interactions occur at high gamma frequencies and to investigate their functional significance, we used event-related causality (ERC) analysis to estimate the dynamics, directionality, and magnitude of event-related causal interactions using subdural electrocorticography (ECoG) recorded during two word production tasks: picture naming and auditory word repetition. A clinical subject who had normal hearing but was skilled in American Signed Language (ASL) provided a unique opportunity to test our hypothesis with reference to a predictable pattern of causal interactions, i.e. that language cortex interacts with different areas of sensorimotor cortex during spoken vs. signed responses. Our ERC analyses confirmed this prediction. During word production with spoken responses, perisylvian language sites had prominent causal interactions with mouth/tongue areas of motor cortex, and when responses were gestured in sign language, the most prominent interactions involved hand and arm areas of motor cortex. Furthermore, we found that the sites from which the most numerous and prominent causal interactions originated, i.e. sites with a pattern of ERC "divergence", were also sites where high gamma power increases were most prominent and where electrocortical stimulation mapping interfered with word production. These findings suggest that the number, strength and directionality of event-related causal interactions may help identify network nodes that are not only activated by a task but are critical to its performance.

Model-based assessment of baroreflex and cardiopulmonary couplings during graded head-up tilt

We propose a multivariate dynamical adjustment (MDA) modeling approach to assess the strength of baroreflex and cardiopulmonary couplings from spontaneous cardiovascular variabilities. Open loop MDA (OLMDA) and closed loop MDA (CLMDA) models were compared. The coupling strength was assessed during progressive sympathetic activation induced by graded head-up tilt. Both OLMDA and CLMDA models suggested that baroreflex coupling progressively increased with tilt table inclination. Only CLMDA model indicated that cardiopulmonary coupling due to the direct link from respiration to heart period gradually decreased with tilt table angles, while that due to the indirect link mediated by systolic arterial pressure progressively increased.

Causality analysis of neural connectivity: critical examination of existing methods and advances of new methods

Granger causality (GC) is one of the most popular measures to reveal causality influence of time series and has been widely applied in economics and neuroscience. Especially, its counterpart in frequency domain, spectral GC, as well as other Granger-like causality measures have recently been applied to study causal interactions between brain areas in different frequency ranges during cognitive and perceptual tasks. In this paper, we show that: 1) GC in time domain cannot correctly determine how strongly one time series influences the other when there is directional causality between two time series, and 2) spectral GC and other Granger-like causality measures have inherent shortcomings and/or limitations because of the use of the transfer function (or its inverse matrix) and partial information of the linear regression model. On the other hand, we propose two novel causality measures (in time and frequency domains) for the linear regression model, called new causality and new spectral causality, respectively, which are more reasonable and understandable than GC or Granger-like measures. Especially, from one simple example, we point out that, in time domain, both new causality and GC adopt the concept of proportion, but they are defined on two different equations where one equation (for GC) is only part of the other (for new causality), thus the new causality is a natural extension of GC and has a sound conceptual/theoretical basis, and GC is not the desired causal influence at all. By several examples, we confirm that new causality measures have distinct advantages over GC or Granger-like measures. Finally, we conduct event-related potential causality analysis for a subject with intracranial depth electrodes undergoing evaluation for epilepsy surgery, and show that, in the frequency domain, all measures reveal significant directional event-related causality, but the result from new spectral causality is consistent with event-related time-frequency power spectrum activity. The spectral GC as well as other Granger-like measures are shown to generate misleading results. The proposed new causality measures may have wide potential applications in economics and neuroscience.

Transfer entropy in magnetoencephalographic data: quantifying information flow in cortical and cerebellar networks

The analysis of cortical and subcortical networks requires the identification of their nodes, and of the topology and dynamics of their interactions. Exploratory tools for the identification of nodes are available, e.g. magnetoencephalography (MEG) in combination with beamformer source analysis. Competing network topologies and interaction models can be investigated using dynamic causal modelling. However, we lack a method for the exploratory investigation of network topologies to choose from the very large number of possible network graphs. Ideally, this method should not require a pre-specified model of the interaction. Transfer entropy--an information theoretic implementation of Wiener-type causality--is a method for the investigation of causal interactions (or information flow) that is independent of a pre-specified interaction model. We analysed MEG data from an auditory short-term memory experiment to assess whether the reconfiguration of networks implied in this task can be detected using transfer entropy. Transfer entropy analysis of MEG source-level signals detected changes in the network between the different task types. These changes prominently involved the left temporal pole and cerebellum--structures that have previously been implied in auditory short-term or working memory. Thus, the analysis of information flow with transfer entropy at the source-level may be used to derive hypotheses for further model-based testing.

A novel approach to propagation pattern analysis in intracardiac atrial fibrillation signals

The purpose of this study is to investigate propagation patterns in intracardiac signals recorded during atrial fibrillation (AF) using an approach based on partial directed coherence (PDC), which evaluates directional coupling between multiple signals in the frequency domain. The PDC is evaluated at the dominant frequency of AF signals and tested for significance using a surrogate data procedure specifically designed to assess causality. For significantly coupled sites, the approach allows also to estimate the delay in propagation. The methods potential is illustrated with two simulation scenarios based on a detailed ionic model of the human atrial myocyte as well as with real data recordings, selected to present typical propagation mechanisms and recording situations in atrial tachyarrhythmias. In both simulation scenarios the significant PDCs correctly reflect the direction of coupling and thus the propagation between all recording sites. In the real data recordings, clear propagation patterns are identified which agree with previous clinical observations. Thus, the results illustrate the ability of the novel approach to identify propagation patterns from intracardiac signals during AF, which can provide important information about the underlying AF mechanisms, potentially improving the planning and outcome of arrhythmia ablation.

Transfer entropy--a model-free measure of effective connectivity for the neurosciences

Understanding causal relationships, or effective connectivity, between parts of the brain is of utmost importance because a large part of the brain's activity is thought to be internally generated and, hence, quantifying stimulus response relationships alone does not fully describe brain dynamics. Past efforts to determine effective connectivity mostly relied on model based approaches such as Granger causality or dynamic causal modeling. Transfer entropy (TE) is an alternative measure of effective connectivity based on information theory. TE does not require a model of the interaction and is inherently non-linear. We investigated the applicability of TE as a metric in a test for effective connectivity to electrophysiological data based on simulations and magnetoencephalography (MEG) recordings in a simple motor task. In particular, we demonstrate that TE improved the detectability of effective connectivity for non-linear interactions, and for sensor level MEG signals where linear methods are hampered by signal-cross-talk due to volume conduction.

Transfer entropy--a model-free measure of effective connectivity for the neurosciences

Understanding causal relationships, or effective connectivity, between parts of the brain is of utmost importance because a large part of the brain’s activity is thought to be internally generated and, hence, quantifying stimulus response relationships alone does not fully describe brain dynamics. Past efforts to determine effective connectivity mostly relied on model based approaches such as Granger causality or dynamic causal modeling. Transfer entropy (TE) is an alternative measure of effective connectivity based on information theory. TE does not require a model of the interaction and is inherently non-linear. We investigated the applicability of TE as a metric in a test for effective connectivity to electrophysiological data based on simulations and magnetoencephalography (MEG) recordings in a simple motor task. In particular, we demonstrate that TE improved the detectability of effective connectivity for non-linear interactions, and for sensor level MEG signals where linear methods are hampered by signal-cross-talk due to volume conduction.

Dynamic Granger causality based on Kalman filter for evaluation of functional network connectivity in fMRI data

Increasing interest in understanding dynamic interactions of brain neural networks leads to formulation of sophisticated connectivity analysis methods. Recent studies have applied Granger causality based on standard multivariate autoregressive (MAR) modeling to assess the brain connectivity. Nevertheless, one important flaw of this commonly proposed method is that it requires the analyzed time series to be stationary, whereas such assumption is mostly violated due to the weakly nonstationary nature of functional magnetic resonance imaging (fMRI) time series. Therefore, we propose an approach to dynamic Granger causality in the frequency domain for evaluating functional network connectivity in fMRI data. The effectiveness and robustness of the dynamic approach was significantly improved by combining a forward and backward Kalman filter that improved estimates compared to the standard time-invariant MAR modeling. In our method, the functional networks were first detected by independent component analysis (ICA), a computational method for separating a multivariate signal into maximally independent components. Then the measure of Granger causality was evaluated using generalized partial directed coherence that is suitable for bivariate as well as multivariate data. Moreover, this metric provides identification of causal relation in frequency domain, which allows one to distinguish the frequency components related to the experimental paradigm. The procedure of evaluating Granger causality via dynamic MAR was demonstrated on simulated time series as well as on two sets of group fMRI data collected during an auditory sensorimotor (SM) or auditory oddball discrimination (AOD) tasks. Finally, a comparison with the results obtained from a standard time-invariant MAR model was provided.

Analyzing the connectivity between regions of interest: an approach based on cluster Granger causality for fMRI data analysis

The identification, modeling, and analysis of interactions between nodes of neural systems in the human brain have become the aim of interest of many studies in neuroscience. The complex neural network structure and its correlations with brain functions have played a role in all areas of neuroscience, including the comprehension of cognitive and emotional processing. Indeed, understanding how information is stored, retrieved, processed, and transmitted is one of the ultimate challenges in brain research. In this context, in functional neuroimaging, connectivity analysis is a major tool for the exploration and characterization of the information flow between specialized brain regions. In most functional magnetic resonance imaging (fMRI) studies, connectivity analysis is carried out by first selecting regions of interest (ROI) and then calculating an average BOLD time series (across the voxels in each cluster). Some studies have shown that the average may not be a good choice and have suggested, as an alternative, the use of principal component analysis (PCA) to extract the principal eigen-time series from the ROI(s). In this paper, we introduce a novel approach called cluster Granger analysis (CGA) to study connectivity between ROIs. The main aim of this method was to employ multiple eigen-time series in each ROI to avoid temporal information loss during identification of Granger causality. Such information loss is inherent in averaging (e.g., to yield a single "representative" time series per ROI). This, in turn, may lead to a lack of power in detecting connections. The proposed approach is based on multivariate statistical analysis and integrates PCA and partial canonical correlation in a framework of Granger causality for clusters (sets) of time series. We also describe an algorithm for statistical significance testing based on bootstrapping. By using Monte Carlo simulations, we show that the proposed approach outperforms conventional Granger causality analysis (i.e., using representative time series extracted by signal averaging or first principal components estimation from ROIs). The usefulness of the CGA approach in real fMRI data is illustrated in an experiment using human faces expressing emotions. With this data set, the proposed approach suggested the presence of significantly more connections between the ROIs than were detected using a single representative time series in each ROI.

Quantifying auditory event-related responses in multichannel human intracranial recordings

Multichannel intracranial recordings are used increasingly to study the functional organization of human cortex. Intracranial recordings of event-related activity, or electrocorticography (ECoG), are based on high density electrode arrays implanted directly over cortex, combining good temporal and spatial resolution. Developing appropriate statistical methods for analyzing event-related responses in these high dimensional ECoG datasets remains a major challenge for clinical and systems neuroscience. We present a novel methodological framework that combines complementary, existing methods adapted for statistical analysis of auditory event-related responses in multichannel ECoG recordings. This analytic framework integrates single-channel (time-domain, time–frequency) and multichannel analyses of event-related ECoG activity to determine statistically significant evoked responses, induced spectral responses, and effective (causal) connectivity. Implementation of this quantitative approach is illustrated using multichannel ECoG data from recent studies of auditory processing in patients with epilepsy. Methods described include a time–frequency matching pursuit algorithm adapted for modeling brief, transient cortical spectral responses to sound, and a recently developed method for estimating effective connectivity using multivariate autoregressive modeling to measure brief event-related changes in multichannel functional interactions. A semi-automated spatial normalization method for comparing intracranial electrode locations across patients is also described. The individual methods presented are published and readily accessible. We discuss the benefits of integrating multiple complementary methods in a unified and comprehensive quantitative approach. Methodological considerations in the analysis of multichannel ECoG data, including corrections for multiple comparisons are discussed, as well as remaining challenges in the development of new statistical approaches.

Nonlinear connectivity by Granger causality

The communication among neuronal populations, reflected by transient synchronous activity, is the mechanism underlying the information processing in the brain. Although it is widely assumed that the interactions among those populations (i.e. functional connectivity) are highly nonlinear, the amount of nonlinear information transmission and its functional roles are not clear. The state of the art to understand the communication between brain systems are dynamic causal modeling (DCM) and Granger causality. While DCM models nonlinear couplings, Granger causality, which constitutes a major tool to reveal effective connectivity, and is widely used to analyze EEG/MEG data as well as fMRI signals, is usually applied in its linear version. In order to capture nonlinear interactions between even short and noisy time series, a few approaches have been proposed. We review them and focus on a recently proposed flexible approach has been recently proposed, consisting in the kernel version of Granger causality. We show the application of the proposed approach on EEG signals and fMRI data.

The effect of filtering on Granger causality based multivariate causality measures

In the past, causality measures based on Granger causality have been suggested for assessing directionality in neural signals. In frequency domain analyses (power or coherence) of neural data, it is common to preprocess the time series by filtering or decimating. However, in other fields, it has been shown theoretically that filtering in combination with Granger causality may lead to spurious or missed causalities. We investigated whether this result translates to multivariate causality methods derived from Granger causality with (a) a simulation study and (b) an application to magnetoencephalographic data. To this end, we performed extensive simulations of the effect of applying different filtering techniques and evaluated the performance of five different multivariate causality measures in combination with two numerical significance measures (random permutation and leave one out method). The analysis included three of the most widely used filters (high-pass, low-pass, notch filter), four different filter types (Butterworth, Chebyshev I and II, elliptic filter), variation of filter order, decimating and interpolation. The simulation results suggest that preprocessing without a strong prior about the artifact to be removed disturbs the information content and time ordering of the data and leads to spurious and missed causalities. Only if apparent artifacts like a current or movement artifact are present, filtering out the respective disturbance seems advisable. While oversampling poses no problem, decimation by a factor greater than the minimum time shift between the time series may lead to wrong inferences. In general, the multivariate causality measures are very sensitive to data preprocessing.

Directed information measure for quantifying the information flow in the brain

In neurophysiology, it is important to determine the causal relationships between neuronal sites. The major problem with existing methods for quantifying the causality in the brain, e.g. Granger causality, is that they assume an multi-variate autoregressive signal model for the multi-channel EEG signals and do not take the nonlinear dependencies between neuronal oscillations into account. In this paper, we propose to quantify the causality of the interactions based on the directed information (DI) criterion, which measures the information flow between two signals over time. The proposed measure is applied to real EEG data from control and schizophrenic subject groups. The significance of the computed DI values are verified by Fourier bootstrapping technique.

A MATLAB toolbox for Granger causal connectivity analysis

Assessing directed functional connectivity from time series data is a key challenge in neuroscience. One approach to this problem leverages a combination of Granger causality analysis and network theory. This article describes a freely available MATLAB toolbox--'Granger causal connectivity analysis' (GCCA)--which provides a core set of methods for performing this analysis on a variety of neuroscience data types including neuroelectric, neuromagnetic, functional MRI, and other neural signals. The toolbox includes core functions for Granger causality analysis of multivariate steady-state and event-related data, functions to preprocess data, assess statistical significance and validate results, and to compute and display network-level indices of causal connectivity including 'causal density' and 'causal flow'. The toolbox is deliberately small, enabling its easy assimilation into the repertoire of researchers. It is however readily extensible given proficiency with the MATLAB language.

GABAergic neurons of the medial septum lead the hippocampal network during theta activity

Information processing in the hippocampus critically relies on its reciprocal interaction with the medial septum (MS). Synchronization of the septo-hippocampal system was demonstrated during both major hippocampal activity states, the regular theta rhythm and the large amplitude irregular activity. Previous experimental and modeling data suggest that the MS provides rhythmic drive to the hippocampus, and hippocampo-septal feedback synchronizes septal pacemaker units. However, this view has recently been questioned based on the possibility of intrahippocampal theta genesis. Previously, we identified putative pacemaker neurons expressing parvalbumin (PV) and/or the pacemaker hyperpolarization-activated and cyclic nucleotide-gated nonselective cation channel (HCN) in the MS. In this study, by analyzing the temporal relationship of activity between the PV/HCN-containing medial septal neurons and hippocampal local field potential, we aimed to uncover whether the sequence of events during theta formation supports the classic view of septal drive or the challenging theory of hippocampal pacing of theta. Importantly, by implementing a circular statistical method, a temporal lead of these septal neurons over the hippocampus was observed on the course of theta synchronization. Moreover, the activity of putative hippocampal interneurons also preceded hippocampal local field theta, but by a shorter time period compared with PV/HCN-containing septal neurons. Using the concept of mutual information, the action potential series of PV/HCN-containing neurons shared higher amount of information with hippocampal field oscillation than PV/HCN-immunonegative cells. Thus, a pacemaker neuron population of the MS leads hippocampal activity, presumably via the synchronization of hippocampal interneurons.

Confounding effects of indirect connections on causality estimation

Addressing the issue of effective connectivity, this study focuses on effects of indirect connections on inferring stable causal relations: partial transfer entropy. We introduce a Granger causality measure based on a multivariate version of transfer entropy. The statistic takes into account the influence of the rest of the network (environment) on observed coupling between two given nodes. This formalism allows us to quantify, for a specific pathway, the total amount of indirect coupling mediated by the environment. We show that partial transfer entropy is a more sensitive technique to identify robust causal relations than its bivariate equivalent. In addition, we demonstrate the confounding effects of the variation in indirect coupling on the detectability of robust causal links. Finally, we consider the problem of model misspecification and its effect on the robustness of the observed connectivity patterns, showing that misspecifying the model may be an issue even for model-free information-theoretic approach.

Parallel versus serial processing dependencies in the perisylvian speech network: a Granger analysis of intracranial EEG data

In this work, we apply Granger causality analysis to high spatiotemporal resolution intracranial EEG (iEEG) data to examine how different components of the left perisylvian language network interact during spoken language perception. The specific focus is on the characterization of serial versus parallel processing dependencies in the dominant hemisphere dorsal and ventral speech processing streams. Analysis of iEEG data from a large, 64-electrode grid implanted over the left perisylvian region in a single right-handed patient showed a consistent pattern of direct posterior superior temporal gyrus influence over sites distributed over the entire ventral pathway for words, non-words, and phonetically ambiguous items that could be interpreted either as words or non-words. For the phonetically ambiguous items, this pattern was overlayed by additional dependencies involving the inferior frontal gyrus, which influenced activation measured at electrodes located in both ventral and dorsal stream speech structures. Implications of these results for understanding the functional architecture of spoken language processing and interpreting the role of the posterior superior temporal gyrus in speech perception are discussed.

Causal modelling and brain connectivity in functional magnetic resonance imaging

Recent advances in data analysis and modeling allow the use of fMRI data to ask not just which brain regions are involved in various cognitive and perceptual tasks, but also how they communicate with each other. Karl Friston examines two different state-of-the-art approaches to modeling brain connectivity using neuroimaging.

Identifying neural drivers with functional MRI: an electrophysiological validation

Whether functional magnetic resonance imaging (fMRI) allows the identification of neural drivers remains an open question of particular importance to refine physiological and neuropsychological models of the brain, and/or to understand neurophysiopathology. Here, in a rat model of absence epilepsy showing spontaneous spike-and-wave discharges originating from the first somatosensory cortex (S1BF), we performed simultaneous electroencephalographic (EEG) and fMRI measurements, and subsequent intracerebral EEG (iEEG) recordings in regions strongly activated in fMRI (S1BF, thalamus, and striatum). fMRI connectivity was determined from fMRI time series directly and from hidden state variables using a measure of Granger causality and Dynamic Causal Modelling that relates synaptic activity to fMRI. fMRI connectivity was compared to directed functional coupling estimated from iEEG using asymmetry in generalised synchronisation metrics. The neural driver of spike-and-wave discharges was estimated in S1BF from iEEG, and from fMRI only when hemodynamic effects were explicitly removed. Functional connectivity analysis applied directly on fMRI signals failed because hemodynamics varied between regions, rendering temporal precedence irrelevant. This paper provides the first experimental substantiation of the theoretical possibility to improve interregional coupling estimation from hidden neural states of fMRI. As such, it has important implications for future studies on brain connectivity using functional neuroimaging.

Our understanding of how the brain works relies on the development of neuropsychological models, which describe how brain activity is coordinated among different regions during the execution of a given task. Knowing the directionality of information transfer between connected regions, and in particular distinguishing neural drivers, or the source of forward connections in the brain, from other brain regions, is critical to refine models of the brain. However, whether functional magnetic resonance imaging (fMRI), the most common technique for imaging brain function, allows one to identify neural drivers remains an open question. Here, we used a rat model of absence epilepsy, a form of nonconvulsive epilepsy that occurs during childhood in humans, showing spontaneous spike-and-wave discharges (nonconvulsive seizures) originating from the first somatosensory cortex, to validate several functional connectivity measures derived from fMRI. Standard techniques estimating interactions directly from fMRI data failed because blood flow dynamics varied between regions. However, we were able to identify the neural driver of spike-and-wave discharges when hemodynamic effects were explicitly removed using appropriate modelling. This study thus provides the first experimental substantiation of the theoretical possibility to improve interregional coupling estimation from hidden neural states of fMRI. As such, it has important implications for future studies on connectivity in the functional neuroimaging literature.

Neural long-range interactions can be distinguished from hemodynamic confounds in functional magnetic resonance imaging using new data analysis techniques that will allow experimental validation of models of brain function.