Assessing cortical functional connectivity by partial directed coherence: simulations and application to real data

Causal pattern recovery from neural spike train data using the Snap Shot Score

We present a new approach to learning directed information flow networks from multi-channel spike train data. A novel scoring function, the Snap Shot Score, is used to assess potential networks with respect to their quality of causal explanation for the data. Additionally, we suggest a generic concept of plausibility in order to assess network learning techniques under partial observability conditions. Examples demonstrate the assessment of networks with the Snap Shot Score, and neural network simulations show its performance in complex situations with partial observability. We discuss the application of the new score to real data and indicate how it can be modified to suit other neural data types.

Choice of multivariate autoregressive model order affecting real network functional connectivity estimate

OBJECTIVE

A realistic simulation exploiting real cortical sources identified from non-invasive extra-cranial recordings in healthy subjects has been considered in order to select the most robust procedure for choosing the correct order of multivariate autoregressive (MVAR) models. Different signal-to-noise ratios filter settings and sampling rates were also tested on the estimate of functional connectivity among the network nodes, in simulated and real cases.

METHODS

Starting from magnetoencephalographic recordings, cortical sources in primary sensorimotor areas of the hand were obtained by functional source separation (FSS). Different criteria for the choice of the model order were compared in the simulated network constructed through one of the FSS-extracted sources and its noise-added delayed copies. In two real cases, a validation of the model order (not known a priori) choice was obtained by comparing the time-frequency properties as depicted by classical non-parametric and MVAR methods at rest, during isometric contraction (stationary states) and while dynamically responding to a sensory stimulation (transient state). For completeness, the whole set of MVAR functional connectivity measures was taken into account, to assess the most suitable for our network description.

RESULTS

That the use of an incorrect model order distorts network functional connectivity estimate was documented both in the realistic simulation and in the two real cases. The Minimal Description Length and Schwartz Bayesian Criterion were selected as the most robust for MVAR model order choice. Partial directed coherence (PDC) was the most suitable method for time-frequency connectivity estimate in the simulated as well as in the real cases, both in stationary and transient states. Moreover, the results of MVAR-based connectivity estimate depend on filter setting in the real case.

CONCLUSIONS

The most robust procedure for choosing the correct MVAR model order was provided. The adjunctive comparison of MVAR to classical methods is recommended to validate the choice in the real case.

SIGNIFICANCE

Correct MVAR model order choice and band filtering play an important role for the correct network connectivity estimate.

Directed cortical information flow during human object recognition: analyzing induced EEG gamma-band responses in brain's source space

The increase of induced gamma-band responses (iGBRs; oscillations >30 Hz) elicited by familiar (meaningful) objects is well established in electroencephalogram (EEG) research. This frequency-specific change at distinct locations is thought to indicate the dynamic formation of local neuronal assemblies during the activation of cortical object representations. As analytically power increase is just a property of a single location, phase-synchrony was introduced to investigate the formation of large-scale networks between spatially distant brain sites. However, classical phase-synchrony reveals symmetric, pair-wise correlations and is not suited to uncover the directionality of interactions. Here, we investigated the neural mechanism of visual object processing by means of directional coupling analysis going beyond recording sites, but rather assessing the directionality of oscillatory interactions between brain areas directly. This study is the first to identify the directionality of oscillatory brain interactions in source space during human object recognition and suggests that familiar, but not unfamiliar, objects engage widespread reciprocal information flow. Directionality of cortical information-flow was calculated based upon an established Granger-Causality coupling-measure (partial-directed coherence; PDC) using autoregressive modeling. To enable comparison with previous coupling studies lacking directional information, phase-locking analysis was applied, using wavelet-based signal decompositions. Both, autoregressive modeling and wavelet analysis, revealed an augmentation of iGBRs during the presentation of familiar objects relative to unfamiliar controls, which was localized to inferior-temporal, superior-parietal and frontal brain areas by means of distributed source reconstruction. The multivariate analysis of PDC evaluated each possible direction of brain interaction and revealed widespread reciprocal information-transfer during familiar object processing. In contrast, unfamiliar objects entailed a sparse number of only unidirectional connections converging to parietal areas. Considering the directionality of brain interactions, the current results might indicate that successful activation of object representations is realized through reciprocal (feed-forward and feed-backward) information-transfer of oscillatory connections between distant, functionally specific brain areas.

Estimate of causality between independent cortical spatial patterns during movement volition in spinal cord injured patients

Static hemodynamic or neuroelectric images of brain regions activated during particular tasks do not convey the information of how these regions communicate to each other. Cortical connectivity estimation aims at describing these interactions as connectivity patterns which hold the direction and strength of the information flow between cortical areas. In this study, we attempted to estimate the causality between distributed cortical systems during a movement volition task in preparation for execution of simple movements by a group of normal healthy subjects and by a group of Spinal Cord Injured (SCI) patients. To estimate the causality between the spatial distributed patterns of cortical activity in the frequency domain, we applied a series of processing steps on the recorded EEG data. From the high-resolution EEG recordings we estimated the cortical waveforms for the regions of interest (ROIs), each representing a selected sensor group population. The solutions of the linear inverse problem returned a series of cortical waveforms for each ROI considered and for each trial analyzed. For each subject, the cortical waveforms were then subjected to Independent Component Analysis (ICA) pre-processing. The independent components obtained by the application of the ThinICA algorithm were further processed by a Partial Directed Coherence algorithm, in order to extract the causality between spatial cortical patterns of the estimated data. The source-target cortical dependencies found in the group of normal subjects were relatively similar in all frequency bands analyzed. For the normal subjects we observed a common source pattern in an ensemble of cortical areas including the right parietal and right lip primary motor areas and bilaterally the primary foot and posterior SMA areas. The target of this cortical network, in the Granger-sense of causality, was shown to be a smaller network composed mostly by the primary foot motor areas and the posterior SMA bilaterally. In the case of the SCI population, both the source and the target cortical patterns had larger sizes than in the normal population. The source cortical areas included always the primary foot and lip motor areas, often bilaterally. In addition, the right parietal area and the bilateral premotor area 6 were also involved. Again, the patterns remained substantially stable across the different frequency bands analyzed. The target cortical patterns observed in the SCI population had larger extensions when compared to the normal ones, since in most cases they involved the bilateral activation of the primary foot movement areas as well as the SMA, the primary lip areas and the parietal cortical areas.