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

Kernel-nonlinear-PDC extends Partial Directed Coherence to detecting nonlinear causal coupling

Here we investigate a new concept, kernel-nonlinear-Partial Directed Coherence, whereby a kernel feature space representation of the data allows detecting nonlinear causal links that are otherwise undetectable through linear modeling. We show that adequate connectivity detection is achievable by applying asympotic decision criteria similar to the ones developed for linear models.

The effect of normalization of Partial Directed Coherence on the statistical assessment of connectivity patterns: a simulation study

Partial Directed Coherence (PDC) is a spectral multivariate estimator for effective connectivity, relying on the concept of Granger causality. Even if its original definition derived directly from information theory, two modifies were introduced in order to provide better physiological interpretations of the estimated networks: i) normalization of the estimator according to rows, ii) squared transformation. In the present paper we investigated the effect of PDC normalization on the performances achieved by applying the statistical validation process on investigated connectivity patterns under different conditions of Signal to Noise ratio (SNR) and amount of data available for the analysis. Results of the statistical analysis revealed an effect of PDC normalization only on the percentages of type I and type II errors occurred by using Shuffling procedure for the assessment of connectivity patterns. No effects of the PDC formulation resulted on the performances achieved during the validation process executed instead by means of Asymptotic Statistic approach. Moreover, the percentages of both false positives and false negatives committed by Asymptotic Statistic are always lower than those achieved by Shuffling procedure for each type of normalization.

Brain-computer interface boosts motor imagery practice during stroke recovery

OBJECTIVE

Motor imagery (MI) is assumed to enhance poststroke motor recovery, yet its benefits are debatable. Brain-computer interfaces (BCIs) can provide instantaneous and quantitative measure of cerebral functions modulated by MI. The efficacy of BCI-monitored MI practice as add-on intervention to usual rehabilitation care was evaluated in a randomized controlled pilot study in subacute stroke patients.

METHODS

Twenty-eight hospitalized subacute stroke patients with severe motor deficits were randomized into 2 intervention groups: 1-month BCI-supported MI training (BCI group, n = 14) and 1-month MI training without BCI support (control group; n = 14). Functional and neurophysiological assessments were performed before and after the interventions, including evaluation of the upper limbs by Fugl-Meyer Assessment (FMA; primary outcome measure) and analysis of oscillatory activity and connectivity at rest, based on high-density electroencephalographic (EEG) recordings.

RESULTS

Better functional outcome was observed in the BCI group, including a significantly higher probability of achieving a clinically relevant increase in the FMA score (p < 0.03). Post-BCI training changes in EEG sensorimotor power spectra (ie, stronger desynchronization in the alpha and beta bands) occurred with greater involvement of the ipsilesional hemisphere in response to MI of the paralyzed trained hand. Also, FMA improvements (effectiveness of FMA) correlated with the changes (ie, post-training increase) at rest in ipsilesional intrahemispheric connectivity in the same bands (p < 0.05).

INTERPRETATION

The introduction of BCI technology in assisting MI practice demonstrates the rehabilitative potential of MI, contributing to significantly better motor functional outcomes in subacute stroke patients with severe motor impairments.

Quantitative Assessment of the Training Improvement in a Motor-Cognitive Task by Using EEG, ECG and EOG Signals

Generally, the training evaluation methods consist in experts supervision and qualitative check of the operator's skills improvement by asking them to perform specific tasks and by verifying the final performance. The aim of this work is to find out if it is possible to obtain quantitative information about the degree of the learning process throughout the training period by analyzing neuro-physiological signals, such as the electroencephalogram, the electrocardiogram and the electrooculogram. In fact, it is well known that such signals correlate with a variety of cognitive processes, e.g. attention, information processing, and working memory. A group of 10 subjects have been asked to train daily with the NASA multi-attribute-task-battery. During such training period the neuro-physiological, behavioral and subjective data have been collected. In particular, the neuro-physiological signals have been recorded on the first (T1), on the third (T3) and on the last training day (T5), while the behavioral and subjective data have been collected every day. Finally, all these data have been compared for a complete overview of the learning process and its relations with the neuro-physiological parameters. It has been shown how the integration of brain activity, in the theta and alpha frequency bands, with the autonomic parameters of heart rate and eyeblink rate could be used as metric for the evaluation of the learning progress, as well as the final training level reached by the subjects, in terms of request of cognitive resources.

Dynamic directed interictal connectivity in left and right temporal lobe epilepsy

OBJECTIVE

There is increasing evidence that epileptic activity involves widespread brain networks rather than single sources and that these networks contribute to interictal brain dysfunction. We investigated the fast-varying behavior of epileptic networks during interictal spikes in right and left temporal lobe epilepsy (RTLE and LTLE) at a whole-brain scale using directed connectivity.

METHODS

In 16 patients, 8 with LTLE and 8 with RTLE, we estimated the electrical source activity in 82 cortical regions of interest (ROIs) using high-density electroencephalography (EEG), individual head models, and a distributed linear inverse solution. A multivariate, time-varying, and frequency-resolved Granger-causal modeling (weighted Partial Directed Coherence) was applied to the source signal of all ROIs. A nonparametric statistical test assessed differences between spike and baseline epochs. Connectivity results between RTLE and LTLE were compared between RTLE and LTLE and with neuropsychological impairments.

RESULTS

Ipsilateral anterior temporal structures were identified as key drivers for both groups, concordant with the epileptogenic zone estimated invasively. We observed an increase in outflow from the key driver already before the spike. There were also important temporal and extratemporal ipsilateral drivers in both conditions, and contralateral only in RTLE. A different network pattern between LTLE and RTLE was found: in RTLE there was a much more prominent ipsilateral to contralateral pattern than in LTLE. Half of the RTLE patients but none of the LTLE patients had neuropsychological deficits consistent with contralateral temporal lobe dysfunction, suggesting a relationship between connectivity changes and cognitive deficits.

SIGNIFICANCE

The different patterns of time-varying connectivity in LTLE and RTLE suggest that they are not symmetrical entities, in line with our neuropsychological results. The highest outflow region was concordant with invasive validation of the epileptogenic zone. This enhanced characterization of dynamic connectivity patterns could better explain cognitive deficits and help the management of epilepsy surgery candidates.

Neuroelectrical correlates of trustworthiness and dominance judgments related to the observation of political candidates

The present research investigates the neurophysiological activity elicited by fast observations of faces of real candidates during simulated political elections. We used simultaneous recording of electroencephalographic (EEG) signals as well as galvanic skin response (GSR) and heart rate (HR) as measurements of central and autonomic nervous systems. Twenty healthy subjects were asked to give judgments on dominance, trustworthiness, and a preference of vote related to the politicians' faces. We used high-resolution EEG techniques to map statistical differences of power spectral density (PSD) cortical activity onto a realistic head model as well as partial directed coherence (PDC) and graph theory metrics to estimate the functional connectivity networks and investigate the role of cortical regions of interest (ROIs). Behavioral results revealed that judgment of dominance trait is the most predictive of the outcome of the simulated elections. Statistical comparisons related to PSD and PDC values highlighted an asymmetry in the activation of frontal cortical areas associated with the valence of the judged trait as well as to the probability to cast the vote. Overall, our results highlight the existence of cortical EEG features which are correlated with the prediction of vote and with the judgment of trustworthy and dominant faces.

Dynamic connectivity among cortical layers in local and large-scale sensory processing

Cortical processing of sensory stimuli typically recruits multiple areas, but how each area dynamically incorporates activity from other areas is not well understood. We investigated interactions between cortical columns of bilateral primary sensory regions (S1s) in rats by recording local field potentials and multi-unit activity simultaneously in both S1s with electrodes positioned at each cortical layer. Using dynamic connectivity analysis based on Granger-causal modeling, we found that, shortly after whisker stimulation (< 10 ms), contralateral S1 (cS1) already relays activity to granular and infragranular layers of S1 in the other hemisphere, after which cS1 shows a pattern of within-column interactions that directs activity upwards toward superficial layers. This pattern of predominant upward driving was also observed in S1 ipsilateral to stimulation, but at longer latencies. In addition, we found that interactions between the two S1s most strongly target granular and infragranular layers. Taken together, the results suggest a possible mechanism for how cortical columns integrate local and large-scale neocortical computation by relaying information from deeper layers to local processing in superficial layers.

Epileptic focus localization based on resting state interictal MEG recordings is feasible irrespective of the presence or absence of spikes

OBJECTIVE

To investigate whether epileptogenic focus localization is possible based on resting state connectivity analysis of magnetoencephalographic (MEG) data.

METHODS

A multivariate autoregressive (MVAR) model was constructed using the sensor space data and was projected to the source space using lead field and inverse matrix. The generalized partial directed coherence was estimated from the MVAR model in the source space. The dipole with the maximum information inflow was hypothesized to be within the epileptogenic focus.

RESULTS

Applying the focus localization algorithm (FLA) to the interictal MEG recordings from five patients with neocortical epilepsy, who underwent presurgical evaluation for the identification of epileptogenic focus, we were able to correctly localize the focus, on the basis of maximum interictal information inflow in the presence or absence of interictal epileptic spikes in the data, with three out of five patients undergoing resective surgery and being seizure free since.

CONCLUSION

Our preliminary results suggest that accurate localization of the epileptogenic focus may be accomplished using noninvasive spontaneous "resting-state" recordings of relatively brief duration and without the need to capture definite interictal and/or ictal abnormalities.

SIGNIFICANCE

Epileptogenic focus localization is possible through connectivity analysis of resting state MEG data irrespective of the presence/absence of spikes.

Assessing dynamic spectral causality by lagged adaptive directed transfer function and instantaneous effect factor

It is of significance to assess the dynamic spectral causality among physiological signals. Several practical estimators adapted from spectral Granger causality have been exploited to track dynamic causality based on the framework of time-varying multivariate autoregressive (tvMVAR) models. The nonzero covariance of the model's residuals has been used to describe the instantaneous effect phenomenon in some causality estimators. However, for the situations with Gaussian residuals in some autoregressive models, it is challenging to distinguish the directed instantaneous causality if the sufficient prior information about the "causal ordering" is missing. Here, we propose a new algorithm to assess the time-varying causal ordering of tvMVAR model under the assumption that the signals follow the same acyclic causal ordering for all time lags and to estimate the instantaneous effect factor (IEF) value in order to track the dynamic directed instantaneous connectivity. The time-lagged adaptive directed transfer function (ADTF) is also estimated to assess the lagged causality after removing the instantaneous effect. In this study, we first investigated the performance of the causal-ordering estimation algorithm and the accuracy of IEF value. Then, we presented the results of IEF and time-lagged ADTF method by comparing with the conventional ADTF method through simulations of various propagation models. Statistical analysis results suggest that the new algorithm could accurately estimate the causal ordering and give a good estimation of the IEF values in the Gaussian residual conditions. Meanwhile, the time-lagged ADTF approach is also more accurate in estimating the time-lagged dynamic interactions in a complex nervous system after extracting the instantaneous effect. In addition to the simulation studies, we applied the proposed method to estimate the dynamic spectral causality on real visual evoked potential (VEP) data in a human subject. Its usefulness in time-variant spectral causality assessment was demonstrated through the mutual causality investigation of brain activity during the VEP experiments.

Frontal EEG theta changes assess the training improvements of novices in flight simulation tasks

The aim of the study is to analyze the variation of the EEG power spectra in theta band when a novice starts to learn a new task. In particular, the goal is to find out the differences from the beginning of the training to the session in which the performance level is good enough for considering him/her able to complete the task without any problems. While the novices were engaged in the flight simulation tasks we recorded the brain activity by using high resolution EEG techniques as well as neurophysiologic variables such as heart rate (HR) and eye blinks rate (EBR). Results show clear changes in the EEG power spectra in theta band over the frontal brain areas, either over the left, the midline and the right side, during the learning process of the task. These results are also supported by the autonomic signals of HR and EBR, by the performances' trends and by the questionnaires for the evaluation of the perceived workload level.

Enhanced effective connectivity in mild occipital stroke patients with hemianopia

Plasticity-based spontaneous recovery and rehabilitation intervention of stroke-induced hemianopia have drawn great attention in recent years. However, the underlying neural mechanism remains unknown. This study aims to investigate brain network disruption and reorganization in hemianopia patients due to mild occipital stroke. Resting-state networks were constructed from 12 hemianopia patients with right occipital infarct by partial directed coherence analysis of multi-channel electroencephalograms. Compared with control subjects, the patients presented enhanced connectivity owing to newly formed connections. Compensational connections mostly originated from the peri-infarct area and targeted contralesional frontal, central, and parietal cortices. These new ipsilesional-to-contralesional inter-hemispheric connections coordinately presented significant correlation with the extent of vision loss. The enhancement of connectivity might be the neural substrate for brain plasticity in stroke-induced hemianopia and may shed light on plasticity-based recovery or rehabilitation.

The physiological plausibility of time-varying Granger-causal modeling: normalization and weighting by spectral power

Time-varying connectivity methods are increasingly used to study directed interactions between brain regions from electrophysiological signals. These methods often show good results in simulated data but it is unclear to what extent connectivity results obtained from real data are physiologically plausible. Here we introduce a benchmark approach using multichannel somatosensory evoked potentials (SEPs) measured across rat cortex, where the structural and functional connectivity is relatively simple and well-understood. Rat SEPs to whisker stimulation are exclusively initiated by contralateral primary sensory cortex (S1), at known latencies, and with activity spread from S1 to specific cortical regions. This allows for a comparison of time-varying connectivity measures according to fixed criteria. We thus evaluated the performance of time-varying Partial Directed Coherence (PDC) and the Directed Transfer Function (DTF), comparing row- and column-wise normalization and the effect of weighting by the power spectral density (PSD). The benchmark approach revealed clear differences between methods in terms of physiological plausibility, effect size and temporal resolution. The results provide a validation of time-varying directed connectivity methods in an animal model and suggest a driving role for ipsilateral S1 in the later part of the SEP. The benchmark SEP dataset is made freely available.

Performance comparison between gPDC and PCMI for measuring directionality of neural information flow

BACKGROUND

General partial directed coherence (gPDC) and permutation conditional mutual information (PCMI) have been widely used to analyze neural activities. These two algorithms are representative of linear and nonlinear methods, respectively. However, there is little known about the difference between their performances in measurements of neural information flow (NIF).

NEW METHOD

Comparison of these two approaches was effectively performed based on the neural mass model (NMM) and real local field potentials.

RESULTS

The results showed that the sensitivity of PCMI was more robust than that of gPDC. The coupling strengths calculated by PCMI were closer to theoretical values in the bidirectional mode of NMM. Furthermore, there was a small Coefficient of Variance (C.V.) for the PCMI results. The gPDC was more sensitive to alterations in the directionality index or the coupling strength of NMM; the gPDC method was more likely to detect a difference between two distinct types of coupling strengths compared to that of PCMI, and gPDC performed well in the identification of the coupling strength in the unidirectional mode.

COMPARISON TO EXISTING METHOD(S)

A comparison between gPDC and PCMI was performed and the advantages of the approaches are discussed.

CONCLUSIONS

The performance of the PCMI is better than that of gPDC in measuring the characteristics of connectivity between neural populations. However, gPDC is recommended to distinguish the differences in connectivity between two states in the same pathway or to detect the coupling strength of the unidirectional mode, such as the hippocampal CA3-CA1 pathway.

Individual cortical connectivity changes after stroke: a resampling approach to enable statistical assessment at single-subject level

One of the main limitations commonly encountered when dealing with the estimation of brain connectivity is the difficulty to perform a statistical assessment of significant changes in brain networks at a single-subject level. This is mainly due to the lack of information about the distribution of the connectivity estimators at different conditions. While group analysis is commonly adopted to perform a statistical comparison between conditions, it may impose major limitations when dealing with the heterogeneity expressed by a given clinical condition in patients. This holds true particularly for stroke when seeking for quantitative measurements of the efficacy of any rehabilitative intervention promoting recovery of function. The need is then evident of an assessment which may account for individual pathological network configuration associated with different level of patients' response to treatment; such network configuration is highly related to the effect that a given brain lesion has on neural networks. In this study we propose a resampling-based approach to the assessment of statistically significant changes in cortical connectivity networks at a single subject level. First, we provide the results of a simulation study testing the performances of the proposed approach under different conditions. Then, to show the sensitivity of the method, we describe its application to electroencephalographic (EEG) data recorded from two post-stroke patients who showed different clinical recovery after a rehabilitative intervention.

Determination of ECoG information flow activity based on Granger causality and Hilbert transformation

Analysis of directional information flow patterns among different regions of the brain is important for investigating the relation between ECoG (electrocorticographic) and mental activity. The objective is to study and evaluate the information flow activity at different frequencies in the primary motor cortex. We employed Granger causality for capturing the future state of the propagation path and direction between recording electrode sites on the cerebral cortex. A grid covered the right motor cortex completely due to its size (approx. 8 cm×8 cm) but grid area extends to the surrounding cortex areas. During the experiment, a subject was asked to imagine performing two activities: movement of the left small finger and/or movement of the tongue. The time series of the electrical brain activity was recorded during these trials using an 8×8 (0.016-300 Hz band with) ECoG platinum electrode grid, which was placed on the contralateral (right) motor cortex. For detection of information flow activity and communication frequencies among the electrodes, we have proposed a method based on following steps: (i) calculation of analytical time series such as amplitude and phase difference acquired from Hilbert transformation, (ii) selection of frequency having highest interdependence for the electrode pairs for the concerned time series over a sliding window in which we assumed time series were stationary, (iii) calculation of Granger causality values for each pair with selected frequency. The information flow (causal influence) activity and communication frequencies between the electrodes in grid were determined and shown successfully. It is supposed that information flow activity and communication frequencies between the electrodes in the grid are approximately the same for the same pattern. The successful employment of Granger causality and Hilbert transformation for the detection of the propagation path and direction of each component of ECoG among different sub-cortex areas were capable of determining the information flow (causal influence) activity and communication frequencies between the populations of neurons successfully.

Time-variant partial directed coherence for analysing connectivity: a methodological study

For the past decade, the detection and quantification of interactions within and between physiological networks has become a priority-in-common between the fields of biomedicine and computer science. Prominent examples are the interaction analysis of brain networks and of the cardiovascular-respiratory system. The aim of the study is to show how and to what extent results from time-variant partial directed coherence analysis are influenced by some basic estimator and data parameters. The impacts of the Kalman filter settings, the order of the autoregressive (AR) model, signal-to-noise ratios, filter procedures and volume conduction were investigated. These systematic investigations are based on data derived from simulated connectivity networks and were performed using a Kalman filter approach for the estimation of the time-variant multivariate AR model. Additionally, the influence of electrooculogram artefact rejection on the significance and dynamics of interactions in 29 channel electroencephalography recordings, derived from a photic driving experiment, is demonstrated. For artefact rejection, independent component analysis was used. The study provides rules to correctly apply particular methods that will aid users to achieve more reliable interpretations of the results.

Describing relevant indices from the resting state electrophysiological networks

The "Default Mode Network" concept was defined, in fMRI field, as a consistent pattern, involving some regions of the brain, which is active during resting state activity and deactivates during attention demanding or goal-directed tasks. Several fMRI studies described its features also correlating the deactivations with the attentive load required for the task execution. Despite the efforts in EEG field, aiming at correlating the spectral features of EEG signals with DMN, an electrophysiological correlate of the DMN hasn't yet been found. In this study we used advanced techniques for functional connectivity estimation for describing the neuroelectrical properties of DMN. We analyzed the connectivity patterns elicited during the rest condition by 55 healthy subjects by means of Partial Directed Coherence. We extracted some graph indexes in order to describe the properties of the resting network in terms of local and global efficiencies, symmetries and influences between different regions of the scalp. Results highlighted the presence of a consistent network, elicited by more than 70% of analyzed population, involving mainly frontal and parietal regions. The properties of the resting network are uniform among the population and could be used for the construction of a normative database for the identification of pathological conditions.

Towards the time varying estimation of complex brain connectivity networks by means of a General Linear Kalman Filter approach

One of the main limitations of the brain functional connectivity estimation methods based on Autoregressive Modeling, like the Granger Causality family of estimators, is the hypothesis that only stationary signals can be included in the estimation process. This hypothesis precludes the analysis of transients which often contain important information about the neural processes of interest. On the other hand, previous techniques developed for overcoming this limitation are affected by problems linked to the dimension of the multivariate autoregressive model (MVAR), which prevents from analysing complex networks like those at the basis of most cognitive functions in the brain. The General Linear Kalman Filter (GLKF) approach to the estimation of adaptive MVARs was recently introduced to deal with a high number of time series (up to 60) in a full multivariate analysis. In this work we evaluated the performances of this new method in terms of estimation quality and adaptation speed, by means of a simulation study in which specific factors of interest were systematically varied in the signal generation to investigate their effect on the method performances. The method was then applied to high density EEG data related to an imaginative task. The results confirmed the possibility to use this approach to study complex connectivity networks in a full multivariate and adaptive fashion, thus opening the way to an effective estimation of complex brain connectivity networks.

EEG frontal asymmetry related to pleasantness of music perception in healthy children and cochlear implanted users

Interestingly, the international debate about the quality of music fruition for cochlear implanted users does not take into account the hypothesis that bilateral users could perceive music in a more pleasant way with respect to monolateral users. In this scenario, the aim of the present study was to investigate if cerebral signs of pleasantness during music perception in healthy child are similar to those observed in monolateral and in bilateral cochlear implanted users. In fact, previous observations in literature on healthy subjects have indicated that variations of the frontal EEG alpha activity are correlated with the perceived pleasantness of the sensory stimulation received (approach-withdrawal theory). In particular, here we described differences between cortical activities estimated in the alpha frequency band for a healthy child and in patients having a monolateral or a bilateral cochlear implant during the fruition of a musical cartoon. The results of the present analysis showed that the alpha EEG asymmetry patterns observed in a healthy child and that of a bilateral cochlear implanted patient are congruent with the approach-withdrawal theory. Conversely, the scalp topographic distribution of EEG power spectra in the alpha band resulting from the monolateral cochlear user presents a different EEG pattern from the normal and bilateral implanted patients. Such differences could be explained at the light of the approach-withdrawal theory. In fact, the present findings support the hypothesis that a monolateral cochlear implanted user could perceive the music in a less pleasant way when compared to a healthy subject or to a bilateral cochlear user.

Modulation of local field potential power of the subthalamic nucleus during isometric force generation in patients with Parkinson's disease

Investigations of local field potentials of the subthalamic nucleus of patients with Parkinson's disease have provided evidence for pathologically exaggerated oscillatory beta-band activity (13-30 Hz) which is amenable to physiological modulation by, e.g., voluntary movement. Previous functional magnetic resonance imaging studies in healthy controls have provided evidence for an increase of subthalamic nucleus blood-oxygenation-level-dependant signal in incremental force generation tasks. However, the modulation of neuronal activity by force generation and its relationship to peripheral feedback remain to be elucidated. We hypothesised that beta-band activity in the subthalamic nucleus is modulated by incremental force generation. Subthalamic nucleus local field potentials were recorded intraoperatively in 13 patients with Parkinson's disease (37 recording sites) during rest and five incremental isometric force generation conditions of the arm with applied loads of 0-400 g (in 100-g increments). Repeated measures analysis of variance (ANOVA) revealed a modulation of local field potential (LFP) power in the upper beta-band (in 24-30 Hz; F(₃.₀₄₂)=4.693, p=0.036) and the gamma-band (in 70-76 Hz; F(₄)=4.116, p=0.036). Granger-causality was computed with the squared partial directed coherence and showed no significant modulation during incremental isometric force generation. Our findings indicate that the upper beta- and gamma-band power of subthalamic nucleus local field potentials are modulated by the physiological task of force generation in patients with Parkinson's disease. This modulation seems to be not an effect of a modulation of peripheral feedback.

Measuring neurophysiological signals in aircraft pilots and car drivers for the assessment of mental workload, fatigue and drowsiness

This paper reviews published papers related to neurophysiological measurements (electroencephalography: EEG, electrooculography EOG; heart rate: HR) in pilots/drivers during their driving tasks. The aim is to summarise the main neurophysiological findings related to the measurements of pilot/driver's brain activity during drive performance and how particular aspects of this brain activity could be connected with the important concepts of "mental workload", "mental fatigue" or "situational awareness". Review of the literature suggests that exists a coherent sequence of changes for EEG, EOG and HR variables during the transition from normal drive, high mental workload and eventually mental fatigue and drowsiness. In particular, increased EEG power in theta band and a decrease in alpha band occurred in high mental workload. Successively, increased EEG power in theta as well as delta and alpha bands characterise the transition between mental workload and mental fatigue. Drowsiness is also characterised by increased blink rate and decreased HR values. The detection of such mental states is actually performed "offline" with accuracy around 90% but not online. A discussion on the possible future applications of findings provided by these neurophysiological measurements in order to improve the safety of the vehicles will be also presented.

The principle of coherence in multi-level brain information processing

Synchronisation has become one of the major scientific tools to explain biological order at many levels of organisation. In systems neuroscience, synchronised subthreshold and suprathreshold oscillatory neuronal activity within and between distributed neuronal assemblies is acknowledged as a fundamental mode of neuronal information processing. Coherent neuronal oscillations correlate with all basic cognitive functions, mediate local and long-range neuronal communication and affect synaptic plasticity. However, it remains unclear how the very fast and complex changes of functional neuronal connectivity necessary for cognition, as mediated by dynamic patterns of neuronal synchrony, could be explained exclusively based on the well-established synaptic mechanisms. A growing body of research indicates that the intraneuronal matrix, composed of cytoskeletal elements and their binding proteins, structurally and functionally connects the synapses within a neuron, modulates neurotransmission and memory consolidation, and is hypothesised to be involved in signal integration via electric signalling due to its charged surface. Theoretical modelling, as well as emerging experimental evidence indicate that neuronal cytoskeleton supports highly cooperative energy transport and information processing based on molecular coherence. We suggest that long-range coherent dynamics within the intra- and extracellular filamentous matrices could establish dynamic ordered states, capable of rapid modulations of functional neuronal connectivity via their interactions with neuronal membranes and synapses. Coherence may thus represent a common denominator of neurophysiological and biophysical approaches to brain information processing, operating at multiple levels of neuronal organisation, from which cognition may emerge as its cardinal manifestation.

How the statistical validation of functional connectivity patterns can prevent erroneous definition of small-world properties of a brain connectivity network

The application of Graph Theory to the brain connectivity patterns obtained from the analysis of neuroelectrical signals has provided an important step to the interpretation and statistical analysis of such functional networks. The properties of a network are derived from the adjacency matrix describing a connectivity pattern obtained by one of the available functional connectivity methods. However, no common procedure is currently applied for extracting the adjacency matrix from a connectivity pattern. To understand how the topographical properties of a network inferred by means of graph indices can be affected by this procedure, we compared one of the methods extensively used in Neuroscience applications (i.e. fixing the edge density) with an approach based on the statistical validation of achieved connectivity patterns. The comparison was performed on the basis of simulated data and of signals acquired on a polystyrene head used as a phantom. The results showed (i) the importance of the assessing process in discarding the occurrence of spurious links and in the definition of the real topographical properties of the network, and (ii) a dependence of the small world properties obtained for the phantom networks from the spatial correlation of the neighboring electrodes.

Testing the asymptotic statistic for the assessment of the significance of Partial Directed Coherence connectivity patterns

Partial Directed Coherence (PDC) is a powerful tool to estimate a frequency domain description of Granger causality between multivariate time series. One of the main limitation of this estimator, however, has been so far the criteria used to assess the statistical significance, which have been obtained through surrogate data approach or arbitrarily imposed thresholds. The aim of this work is to test the performances of a validation approach based on the rigorous asymptotic distributions of PDC, recently proposed in literature. The performances of this method, defined in terms of percentages of false positives and false negatives, were evaluated by means of a simulation study taking into account factors like the Signal to Noise Ratio (SNR) and the amount of data available for the estimation and the use of different methods for the statistical corrections for multiple comparisons. Results of the Analysis Of Variance (ANOVA) performed on false positives and false negatives revealed a strong dependency of the performances from all the factors investigated. In particular, results indicate an amount of Type I errors below 7% for all conditions, while Type II errors are below 10% when the SNR is at least 1, the data length of at least 50 seconds and the appropriate correction for multiple comparisons is applied.

Study of the functional hyperconnectivity between couples of pilots during flight simulation: an EEG hyperscanning study

Brain Hyperscanning, i.e. the simultaneous recording of the cerebral activity of different human subjects involved in interaction tasks, is a very recent field of Neuroscience aiming at understanding the cerebral processes generating and generated by social interactions. This approach allows the observation and modeling of the neural signature specifically dependent on the interaction between subjects, and, even more interestingly, of the functional links existing between the activities in the brains of the subjects interacting together. In this EEG hyperscanning study we explored the functional hyperconnectivity between the activity in different scalp sites of couples of Civil Aviation Pilots during different phases of a flight reproduced in a flight simulator. Results shown a dense network of connections between the two brains in the takeoff and landing phases, when the cooperation between them is maximal, in contrast with phases during which the activity of the two pilots was independent, when no or quite few links were shown. These results confirms that the study of the brain connectivity between the activity simultaneously acquired in human brains during interaction tasks can provide important information about the neural basis of the "spirit of the group".

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.

Cortical activity and functional hyperconnectivity by simultaneous EEG recordings from interacting couples of professional pilots

Controlling an aircraft during a flight is a compelling condition, which requires a strict and well coded interaction between the crew. The interaction level between the Captain and the First Officer changes during the flight, ranging from a maximum (during takeoff and landing, as well as in case of a failure of the instrumentation or other emergency situations) to a minimum during quiet mid-flight. In this study, our aim is to investigate the neural correlates of different kinds and levels of interaction between couples of professional crew members by means of the innovative technique called brain hyperscanning, i.e. the simultaneous recording of the hemodynamic or neuroelectrical activity of different human subjects involved in interaction tasks. This approach allows the observation and modeling of the neural signature specifically dependent on the interaction between subjects, and, even more interestingly, of the functional links existing between the brain activities of the subjects interacting together. In this EEG hyperscanning study, different phases of a flight were reproduced in a professional flight simulator, which allowed, on one side, to reproduce the ecological setting of a real flight, and, on the other, to keep under control the different levels of interaction induced in the crew by means of systematic and simulated failures of the aircraft instrumentation. Results of the procedure of linear inverse estimation, together with functional hyperconnectivity estimated by means of Partial Directed Coherence, showed a dense network of connections between the activity in the two brains in the takeoff and landing phases, when the cooperation between the crew is maximal, while conversely no significant links were shown during the phases in which the activity of the two pilots was independent.

Assessment of mental fatigue during car driving by using high resolution EEG activity and neurophysiologic indices

Driving tasks are vulnerable to the effects of sleep deprivation and mental fatigue, diminishing driver's ability to respond effectively to unusual or emergent situations. Physiological and brain activity analysis could help to understand how to provide useful feedback and alert signals to the drivers for avoiding car accidents. In this study we analyze the insurgence of mental fatigue or drowsiness during car driving in a simulated environment by using high resolution EEG techniques as well as neurophysiologic variables such as heart rate (HR) and eye blinks rate (EBR). Results suggest that it is possible to introduce a EEG-based cerebral workload index that it is sensitive to the mental efforts of the driver during drive tasks of different levels of difficulty. Workload index was based on the estimation of increase of EEG power spectra in the theta band over prefrontal areas and the simultaneous decrease of EEG power spectra over parietal areas in alpha band during difficult drive conditions. Such index could be used in a future to assess on-line the mental state of the driver during the drive task.

Extending transfer entropy improves identification of effective connectivity in a spiking cortical network model

Transfer entropy (TE) is an information-theoretic measure which has received recent attention in neuroscience for its potential to identify effective connectivity between neurons. Calculating TE for large ensembles of spiking neurons is computationally intensive, and has caused most investigators to probe neural interactions at only a single time delay and at a message length of only a single time bin. This is problematic, as synaptic delays between cortical neurons, for example, range from one to tens of milliseconds. In addition, neurons produce bursts of spikes spanning multiple time bins. To address these issues, here we introduce a free software package that allows TE to be measured at multiple delays and message lengths. To assess performance, we applied these extensions of TE to a spiking cortical network model (Izhikevich, 2006) with known connectivity and a range of synaptic delays. For comparison, we also investigated single-delay TE, at a message length of one bin (D1TE), and cross-correlation (CC) methods. We found that D1TE could identify 36% of true connections when evaluated at a false positive rate of 1%. For extended versions of TE, this dramatically improved to 73% of true connections. In addition, the connections correctly identified by extended versions of TE accounted for 85% of the total synaptic weight in the network. Cross correlation methods generally performed more poorly than extended TE, but were useful when data length was short. A computational performance analysis demonstrated that the algorithm for extended TE, when used on currently available desktop computers, could extract effective connectivity from 1 hr recordings containing 200 neurons in ∼5 min. We conclude that extending TE to multiple delays and message lengths improves its ability to assess effective connectivity between spiking neurons. These extensions to TE soon could become practical tools for experimentalists who record hundreds of spiking neurons.

On the use of EEG or MEG brain imaging tools in neuromarketing research

Here we present an overview of some published papers of interest for the marketing research employing electroencephalogram (EEG) and magnetoencephalogram (MEG) methods. The interest for these methodologies relies in their high-temporal resolution as opposed to the investigation of such problem with the functional Magnetic Resonance Imaging (fMRI) methodology, also largely used in the marketing research. In addition, EEG and MEG technologies have greatly improved their spatial resolution in the last decades with the introduction of advanced signal processing methodologies. By presenting data gathered through MEG and high resolution EEG we will show which kind of information it is possible to gather with these methodologies while the persons are watching marketing relevant stimuli. Such information will be related to the memorization and pleasantness related to such stimuli. We noted that temporal and frequency patterns of brain signals are able to provide possible descriptors conveying information about the cognitive and emotional processes in subjects observing commercial advertisements. These information could be unobtainable through common tools used in standard marketing research. We also show an example of how an EEG methodology could be used to analyze cultural differences between fruition of video commercials of carbonated beverages in Western and Eastern countries.

Reliability of multivariate causality measures for neural data

In the past decade several multivariate causality measures based on Granger causality have been suggested to assess directionality of neural signals. To date, however, a detailed evaluation of the reliability of these measures is largely missing. We systematically evaluated the performance of five different causality measures (squared partial directed coherence (sPDC), partial directed coherence (PDC), directed transfer function (DTF), direct directed transfer function (dDTF) and transfer function) depending upon data length, noise level, coupling strength, and model order and performed simulations based on four different neural data recording procedures (magnetoencephalography, electroencephalography, electromyography, intraoperative local field potentials). Moreover, we analyzed the effect of two common numerical methods to determine the significance of the particular causality measure (random permutation and the leave one out method (LOOM)). The simulations showed the sPDC combined with the LOOM to be the most reliable and robust choice for assessing directionality in neural data. While DTF and H by construction were unable to distinguish between direct and indirect connections, the dDTF also failed this test. Finally, we applied the causality measures to a real data set. This showed the usefulness of our simulation results for practical applications in order to draw correct inferences and distinguish between conflicting evidence obtained with different causality measures.

Large-scale cortical networks estimated from scalp EEG signals during performance of goal-directed motor tasks

The evaluation of the topological properties of brain networks is an emergent research topic, since the estimated cerebral connectivity patterns often have relatively large size and complex structure. Since a graph is a mathematical representation of a network, the use of a theoretical graph approach would describe concisely the topological features of the functional brain connectivity network estimated using neuroimaging techniques. In the present study, we analyze the changes in brain synchronization networks using high-resolution EEG signals obtained during performance of a complex goal-directed visuomotor task. Our results show that the cortical network is more stable when subjects reach the goal than when they fail by hitting an obstacle. These findings suggest the presence of a possible cerebral "marker" for motor actions that result in successful reaching of a target.

Spectral EEG frontal asymmetries correlate with the experienced pleasantness of TV commercial advertisements

The aim of this research is to analyze the changes in the EEG frontal activity during the observation of commercial videoclips. In particular, we aimed to investigate the existence of EEG frontal asymmetries in the distribution of the signals' power spectra related to experienced pleasantness of the video, as explicitly rated by the eleven experimental subjects investigated. In the analyzed population, maps of Power spectral density (PSD) showed an asymmetrical increase of theta and alpha activity related to the observation of pleasant (unpleasant) advertisements in the left (right) hemisphere. A correlation analysis revealed that the increase of PSD at left frontal sites is negatively correlated with the degree of pleasantness perceived. Conversely, the de-synchronization of left alpha frontal activity is positively correlated with judgments of high pleasantness. Moreover, our data presented an increase of PSD related to the observation of unpleasant commercials, which resulted higher with respect to the one elicited by pleasant advertisements.

Semiparametric detection of nonlinear causal coupling using partial directed coherence

Infering causal relationships from observed time series has attracted much recent attention. In cases of nonlinear coupling, adequate inference is often hindered by the need to specify coupling details that call for many parameters and global minimization of nonconvex functions. In this paper we use an example to investigate a new concept, termed here running entropy mapping, whereby time series are mapped onto other entropy related time sequences whose analysis via a linear parametric time series methods, such as partial directed coherence, is able to expose the presence of formerly linearly undetectable causal relationships.

Extended causal modeling to assess Partial Directed Coherence in multiple time series with significant instantaneous interactions

The Partial Directed Coherence (PDC) and its generalized formulation (gPDC) are popular tools for investigating, in the frequency domain, the concept of Granger causality among multivariate (MV) time series. PDC and gPDC are formalized in terms of the coefficients of an MV autoregressive (MVAR) model which describes only the lagged effects among the time series and forsakes instantaneous effects. However, instantaneous effects are known to affect linear parametric modeling, and are likely to occur in experimental time series. In this study, we investigate the impact on the assessment of frequency domain causality of excluding instantaneous effects from the model underlying PDC evaluation. Moreover, we propose the utilization of an extended MVAR model including both instantaneous and lagged effects. This model is used to assess PDC either in accordance with the definition of Granger causality when considering only lagged effects (iPDC), or with an extended form of causality, when we consider both instantaneous and lagged effects (ePDC). The approach is first evaluated on three theoretical examples of MVAR processes, which show that the presence of instantaneous correlations may produce misleading profiles of PDC and gPDC, while ePDC and iPDC derived from the extended model provide here a correct interpretation of extended and lagged causality. It is then applied to representative examples of cardiorespiratory and EEG MV time series. They suggest that ePDC and iPDC are better interpretable than PDC and gPDC in terms of the known cardiovascular and neural physiologies.

Influences of brain development and ageing on cortical interactive networks

OBJECTIVE

To study the effect of brain development and ageing on the pattern of cortical interactive networks.

METHODS

By causality analysis of multichannel electroencephalograph (EEG) with partial directed coherence (PDC), we investigated the different neural networks involved in the whole cortex as well as the anterior and posterior areas in three age groups, i.e., children (0-10 years), mid-aged adults (26-38 years) and the elderly (56-80 years).

RESULTS

By comparing the cortical interactive networks in different age groups, the following findings were concluded: (1) the cortical interactive network in the right hemisphere develops earlier than its left counterpart in the development stage; (2) the cortical interactive network of anterior cortex, especially at C3 and F3, is demonstrated to undergo far more extensive changes, compared with the posterior area during brain development and ageing; (3) the asymmetry of the cortical interactive networks declines during ageing with more loss of connectivity in the left frontal and central areas.

CONCLUSIONS

The age-related variation of cortical interactive networks from resting EEG provides new insights into brain development and ageing.

SIGNIFICANCE

Our findings demonstrated that the PDC analysis of EEG is a powerful approach for characterizing the cortical functional connectivity during brain development and ageing.

Cross-correlation of instantaneous amplitudes of field potential oscillations: a straightforward method to estimate the directionality and lag between brain areas

Researchers performing multi-site recordings are often interested in identifying the directionality of functional connectivity and estimating lags between sites. Current techniques for determining directionality require spike trains or involve multivariate autoregressive modeling. However, it is often difficult to sample large numbers of spikes from multiple areas simultaneously, and modeling can be sensitive to noise. A simple, model-independent method to estimate directionality and lag using local field potentials (LFPs) would be of general interest. Here we describe such a method using the cross-correlation of the instantaneous amplitudes of filtered LFPs. The method involves four steps. First, LFPs are band-pass filtered; second, the instantaneous amplitude of the filtered signals is calculated; third, these amplitudes are cross-correlated and the lag at which the cross-correlation peak occurs is determined; fourth, the distribution of lags obtained is tested to determine if it differs from zero. This method was applied to LFPs recorded from the ventral hippocampus and the medial prefrontal cortex in awake behaving mice. The results demonstrate that the hippocampus leads the mPFC, in good agreement with the time lag calculated from the phase locking of mPFC spikes to vHPC LFP oscillations in the same dataset. We also compare the amplitude cross-correlation method to partial directed coherence, a commonly used multivariate autoregressive model-dependent method, and find that the former is more robust to the effects of noise. These data suggest that the cross-correlation of instantaneous amplitude of filtered LFPs is a valid method to study the direction of flow of information across brain areas.

A graph-theoretical approach in brain functional networks. Possible implications in EEG studies

BACKGROUND

Recently, it was realized that the functional connectivity networks estimated from actual brain-imaging technologies (MEG, fMRI and EEG) can be analyzed by means of the graph theory, that is a mathematical representation of a network, which is essentially reduced to nodes and connections between them.

METHODS

We used high-resolution EEG technology to enhance the poor spatial information of the EEG activity on the scalp and it gives a measure of the electrical activity on the cortical surface. Afterwards, we used the Directed Transfer Function (DTF) that is a multivariate spectral measure for the estimation of the directional influences between any given pair of channels in a multivariate dataset. Finally, a graph theoretical approach was used to model the brain networks as graphs. These methods were used to analyze the structure of cortical connectivity during the attempt to move a paralyzed limb in a group (N=5) of spinal cord injured patients and during the movement execution in a group (N=5) of healthy subjects.

RESULTS

Analysis performed on the cortical networks estimated from the group of normal and SCI patients revealed that both groups present few nodes with a high out-degree value (i.e. outgoing links). This property is valid in the networks estimated for all the frequency bands investigated. In particular, cingulate motor areas (CMAs) ROIs act as ''hubs'' for the out fl ow of information in both groups, SCI and healthy. Results also suggest that spinal cord injuries affect the functional architecture of the cortical network sub-serving the volition of motor acts mainly in its local feature property.In particular, a higher local efficiency El can be observed in the SCI patients for three frequency bands, theta (3-6 Hz), alpha (7-12 Hz) and beta (13-29 Hz).By taking into account all the possible pathways between different ROI couples, we were able to separate clearly the network properties of the SCI group from the CTRL group. In particular, we report a sort of compensatory mechanism in the SCI patients for the Theta (3-6 Hz) frequency band, indicating a higher level of "activation" Omega within the cortical network during the motor task. The activation index is directly related to diffusion, a type of dynamics that underlies several biological systems including possible spreading of neuronal activation across several cortical regions.

CONCLUSIONS

The present study aims at demonstrating the possible applications of graph theoretical approaches in the analyses of brain functional connectivity from EEG signals. In particular, the methodological aspects of the i) cortical activity from scalp EEG signals, ii) functional connectivity estimations iii) graph theoretical indexes are emphasized in the present paper to show their impact in a real application.

On the use of dynamic Bayesian networks in reconstructing functional neuronal networks from spike train ensembles

Coordination among cortical neurons is believed to be a key element in mediating many high-level cortical processes such as perception, attention, learning, and memory formation. Inferring the structure of the neural circuitry underlying this coordination is important to characterize the highly nonlinear, time-varying interactions between cortical neurons in the presence of complex stimuli. In this work, we investigate the applicability of dynamic Bayesian networks (DBNs) in inferring the effective connectivity between spiking cortical neurons from their observed spike trains. We demonstrate that DBNs can infer the underlying nonlinear and time-varying causal interactions between these neurons and can discriminate between mono- and polysynaptic links between them under certain constraints governing their putative connectivity. We analyzed conditionally Poisson spike train data mimicking spiking activity of cortical networks of small and moderately large size. The performance was assessed and compared to other methods under systematic variations of the network structure to mimic a wide range of responses typically observed in the cortex. Results demonstrate the utility of DBN in inferring the effective connectivity in cortical networks.

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.

Causality between local field potentials of the subthalamic nucleus and electromyograms of forearm muscles in Parkinson's disease

Deep brain stimulation of the subthalamic nucleus is an effective treatment for Parkinson's disease, although its precise mechanisms remain poorly understood. To gain further insight into the mechanisms underlying deep brain stimulation, we analysed the causal relationship between forearm muscle activity and local field potentials derived from the subthalamic nucleus. In 19 patients suffering from Parkinson's disease of the akinetic-rigid subtype, we calculated the squared partial directed coherence between muscles of the contralateral forearm and the subthalamic nucleus or zona incerta during both a rest and a hold condition of the arm. For both recording regions, data analysis revealed that, during the rest condition, electromyographic activity was significantly more often 'Granger-causal' for the local field potentials than the opposite causation. In contrast, during the hold condition, no significant difference was found in the occurrence of causalities. Contrary to the existing basal ganglia model and the current concept of Parkinson's disease pathophysiology, we found the subthalamic nucleus to receive more 'afferences' than it emitted 'efferences', suggesting that its role is more complex than a simple driving nucleus in the basal ganglia loop. Therefore, the effect of deep brain stimulation in the subthalamic nucleus could, at least in part, result from a blockade of pathological afferent input.

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.