Original Research ARTICLE
Combining Partial Directed Coherence and Graph Theory to Analyse Effective Brain Networks of Different Mental Tasks
- 1School of Life Science and Technology, Xidian University, Xi'an, China
- 2Department of Radiology, Medical Physics, University Medical Center Freiburg, Freiburg, Germany
- 3School of Electronic Engineering, Xidian University, Xi'an, China
Purpose: The aim of this study is to qualify the network properties of the brain networks between two different mental tasks (play task or rest task) in a healthy population.
Methods and Materials: EEG signals were recorded from 19 healthy subjects when performing different mental tasks. Partial directed coherence (PDC) analysis, based on Granger causality (GC), was used to assess the effective brain networks during the different mental tasks. Moreover, the network measures, including degree, degree distribution, local and global efficiency in delta, theta, alpha, and beta rhythms were calculated and analyzed.
Results: The local efficiency is higher in the beta frequency and lower in the theta frequency during play task whereas the global efficiency is higher in the theta frequency and lower in the beta frequency in the rest task.
Significance: This study reveals the network measures during different mental states and efficiency measures may be used as characteristic quantities for improvement in attentional performance.
Executive functioning (EF) is an umbrella term that refers to high level cognitive functions, which is of great importance to maintain effective goal directed behavior (Welsh and Pennington, 1988; Hughes and Graham, 2002; Alvarez and Emory, 2006). During the last decade it has become a very active research field with over 12000 scientific published articles. EF incorporates (a) the capability to inhibit or delay a particular response, (b) control the mental representation of the task with working memory, and (c) plans for action sequences (Pennington and Ozonoff, 1996; Piek et al., 2004).
One challenge in this field is the question whether executive functions are modulated by the frontal lobes (Welsh, 2002; Alvarez and Emory, 2006), producing ambiguity of definition. A mainstream view is that individuals who perform badly on executive function have “frontal lobe deficit” (Duke and Kaszniak, 2000; Stuss, 2000). However, some reserchers found that persons with frontal lesions can still perform like healthy people on executive function tests (Shallice and Burgess, 1991; Ahola et al., 1996). Moreover, persons with diffuse lesions or non-frontal lesions were even reported to perform similarly to those who had frontal lesions on these tests (Anderson et al., 1991; Axelrod et al., 1996). Furthermore, based on a meta-analytic review of previous woks, Alvarez and Emory (2006) claimed that frontal brain, as well as non-frontal brain regions, are both essential to intact executive functions which should be regarded as a more integrative process, and not merely limited to frontal brain regions.
Brain networks can be investigated with diffusion magnetic resonance imaging (DTI), functional MRI (fMRI), magnetoencephalography (MEG), and electroencephalography (EEG) (Bullmore and Sporns, 2009; Rubinov and Sporns, 2010; Lin et al., 2014). However, in contrast to the other neuroimaging techniques, EEG still remains the cheapest and most widespread technique so far. In particular, its high resolution in time domain allows the observation of signal oscillation in a scale of miliseconds.
In previous attempts to investigate the brain network organization with EEG signals, a large number of conventional methods, such as the techniques of time-domain, frequency-domain, correlation, nonlinear dynamics, and synchronism were proposed (Hjorth, 1970; Aertsen and Gerstein, 1985; Sandkühler and Bhattacharya, 2008; Knyazeva et al., 2010). These methods clarified the network features in the analysis of directionless functional connectivity. However, they were not able to provide direct insights into the directional flows of information among network nodes. In fact, the human brain is a complex dynamic system and constantly responds to the external stimuli (Fox et al., 2005). Thus, it is crucial to explore the directional information flows of the instantaneous interactions of the brain among several nodes.
Previous studies have provided several kinds of sophisticated EEG analysis techniques to study the directional interactions between any pair of network nodes or cortical areas (Baccalá and Sameshima, 2001; Kamiński et al., 2001; Schelter et al., 2006; Supp et al., 2007; De Vico Fallani et al., 2010, 2007). For example, partial directed coherence (PDC) on the basis of the Granger causality (GC) (Granger, 1969) and the multivariate autoregressive (MVAR) model (Blinowska et al., 2004) has attracted more attention since it can be adopted to characterize the effective connections, representing direct or indirect causal influences of one region to another (Rubinov and Sporns, 2010). Using graph theoretical approach (Sporns et al., 2004; Bullmore and Sporns, 2009; Van Dellen et al., 2009), topological properties of complex networks were characterized (Ahmadlou et al., 2012; Zhang et al., 2013, 2015; Gast et al., 2014; Gupta and Falk, 2015; Liu et al., 2015; Youssofzadeh et al., 2016).
The experiment in this study was designed to investigate the attention aspect of the effective function. The aim of the current work is to investigate the brain networks using electrodes related to different mental tasks in a healthy population. Moreover, the question is whether the effective network, evaluated by graph analysis, may be different based on the above chosen electrodes. In particular, the main experimental questions we would like to know from the study are the following:
1. Is the efficiency index significantly different in the obtained brain networks in different mental states?
2. If it does exist, is such a difference related to a specific frequency?
To the best of our knowledge, this is the first study that combined PDC and graph theory to evaluate the effective networks properties toward different mental states in a healthy population. In particular, the efficiency indexes with regard to the four EEG rhythms (delta, theta, alpha, and beta) were closely investigated. Moreover, it shed light on the brain network toward executive functions in a healthy population.
Methods and Materials
The nineteen university students participating in the study were recruited from Xidian University (mean age ±SD = 24.25 ± 1.57 years; right-handed; male; visual acuity or corrected visual acuity no less than 1.0). They had a full-scale Raven's Standard Progressive Matrices score of 100 or above. Written, informed consent was obtained from each subject after the explanation of the study, which was approved by the Ethical committee of the institute of life science and technology of Xidian University. The subjects had no personal history of neurological disease, clinically significant head trauma, a substance use disorder, or psychiatric illness. The participants were not allowed to take any medication, coffee, and alcohol for a minimum of 4 h before the experiment.
Two mental tasks have been designed for participants as follows (Yan et al., 2007):
A) Play task: participants were comfortably seated on a sofa in front of a computer, about 72 cm from the screen. Each subject was asked to focus the attention on a game called “feeding frenzy.” The participants had to use the mouse to control “their fish.” The controlled fish swam around the sea, adding points by eating smaller fish but avoiding bigger fish simultaneously. One session lasted for 3 min. To ensure the quality of the data, another session was recorded.
B) Rest task: all the surroundings and requirements were the same as in the play task. The only difference was that the subjects were asked to sit in silence while watching the still image of the game on the screen and should not think about anything special. One session lasted for 3 min. Two sessions were recorded.
EEG data were recorded at 16 scalp loci (FP1, FP2, F3, F4, C3, C4, P3, P4, O1, O2, F7, F8, T3, T4, T5, T6, and A1+A2 as reference) in line with the international 10–20 system to linked earlobes. The experiments were performed in a normal chamber without noise interference. The electrodes and the positions connecting with the electrodes were defatted with alcohol to reduce the impedance and interference. Impedance was kept below 10 kΩ. EEG signals of each subject were recorded with a 32-channel EEG system (NT9200, Beijing, China) when performing different mental tasks. The signals were recorded with a bandpass of 0.05–30 Hz, 50 Hz notch, and digitized at 1000 Hz.
The whole data processing work flow is summarized in Figure 1. The pre-processing included: (1) The EEG data were down-sampled to 250 Hz to reduce the volume of the data. (2) The frequent swallowing and winking fragments of the EEG data were diminished by means of independent component analysis (ICA) employing the function runica from EEGLAB (Delorme and Makeig, 2004).
Electrodes Selection by LMD and SVM
Some investigators suggested that (Pardo et al., 1991; Buschman and Miller, 2007) the activation of frontal and parietal cortical areas are involved in the process of attention. However, it is not sure yet how many of the electrodes should be accurately included. To reduce the number of electrodes for future clinical use and eliminate a source of subjectivity, in this work, the technique combining LMD and SVM was used to classify the two mental tasks for each subject.
LMD (Smith, 2005; Park et al., 2011) is an iterative method to demodulate amplitude and frequency modulated signals. It decomposes the EEG signal into a set of components or functions. Each of the components is the product of an envelope signal and a frequency modulated signal. In our study, the energy of the first four components was summed as the input vector for the SVM (Burges, 1998; Noble, 2006; Hsu et al., 2008). Our SVM model was developed in LIBSVM (version 3.1) with Radial basis function (RBF) kernel parameters. To choose the optimal parameters for gamma-g and cost-c, which affected the performance of SVM model greatly, the genetic algorithm (GA) was used for every data set in a cross validation manner. In particular, leave-one-out cross validation was used in our study to ensure the results reliable and repeatable.
The results of SVM classification for all the EEG electrodes (Table 1) show that different EEG electrodes result in different classification precisions, and 10 electrodes (F3, F4, C3, C4, P3, P4, T3, T4, T5, and T6) with classification precision higher than 60% were chosen as the network nodes.
Partial Directed Coherence (PDC)
The concept of PDC, proposed by Baccala and Sameshima in 2001 (Baccalá and Sameshima, 2001), put forward a new frequency-domain method for the description of GC.
Generally, a MVAR model with m channels of EEG signals and order p is defined by,
Where X(t) is the vector of m channels of EEG signals at time t and the matrix A(r) contains the rth order AR parameters. E(t) denotes the estimated error, which is supposed to be an uncorrelated Gaussian process with zero mean. The fitting results of MVAR model will be influenced by the model order p and Akaike information criterion (AIC) can be used to determine the proper model order p for our study.
Once the coefficients of the MVAR model were adequately estimated, A(f) can be obtained by . Further, the transfer function of the m channels of EEG signals can be quantified by , where the elements of A(f) are defined as:
Finally, the PDC from channel j to i can be given by,
where āi(f)(i = 1, 2, …M) is the ith column of the matrix Ā(f) and PDCij denotes the direction and intensity of the information flow from channel j to i at the frequency of f.
Given that the sliding time window techniques might disclose higher resolution dynamics (Ding et al., 2000), a sliding time window of 2 s of EEG data with 50% overlapped was used in our study. The EEGLAB toolbox, which provides functions to calculate AIC (Cui et al., 2008; Janacek, 2010), was used to determine the proper order of the MVAR model for each subject. Furthermore, a bootstrap method (Kamiński et al., 2001) was used to determine the statistical significance regarding to the above interactive networks.
Creating Adjacency Matrices
So far, the connectivity matrices have been created with values of PDC between all pairwise associations for each EEG sub-band (delta, theta, alpha, and beta) in different mental states. Further, a threshold was applied to convert the connectivity matrices to adjacency matrices.
The choice of threshold is crucial in the process of generating an adjacency matrix from the association matrix. Conceivably, different thresholds will result in graphs of different connection densities or sparsities. Although there are currently several techniques able to select the threshold, e.g., surrogate-data technique (Pritchard et al., 2014), significance level based technique (De Vico Fallani et al., 2007), and so on, there is currently no consent on which one is optimal.
Here, a method based on significance level was employed to determine the threshold T in this study. The effective network was converted into a directed unweighted graph and the connection matrix containing PDC values for each directed pair of nodes into an adjacency matrix A by means of a threshold T. The threshold T can be viewed as the ratio of the number of real effective connections to the number of all possible connections in this network, and T indicates the number of the most powerful interactions to be considered. In order to explore the features of the effective networks at different connection densities, the threshold T was slid within the interval 0.1–0.9 with a step of 0.05 and the average values of local efficiency in beta rhythm were obtained under the related threshold between the two mental states.
In this paper, degree, degree distribution and efficiency indexes were calculated with the Brain Connectivity Toolbox (Rubinov and Sporns, 2010) under Matlab R2010a. In particular, efficiency is a quantity that measures how efficiently the information is exchanged over the network (Latora and Marchiori, 2001; Fan et al., 2002; Goulas et al., 2015). The efficiency eij in the communication between vertices i and j is defined as eij = 1∕dij ∀i, j, where dij is the shortest distance between the nodes i and j. If there is no path in the graph between the nodes i and j then dij = ∞, and therefore eij = 0.
The global efficiency of A is defined as: , where is the number of vertices in the graph. The local properties of A can be characterized by evaluating for each node i the efficiency of Ai (the sub-graph of the neighbors of i). The local efficiency is the average efficiency of the local sub-graphs, i.e., . This property tells us the “fault-tolerant” of the network system, to be more specific, how efficient is the communication between the first neighbors of i when i is removed (De Vico Fallani et al., 2007; Klados et al., 2012, 2013).
In the study, the optimal model order for each subject obtained from AIC was in the range between 10 and 16. The connectivity matrixes were created with the values of PDC between all pairwise electrodes for each EEG sub-band following a subject specific MVAR order in different mental states.
Figure 2 denotes the average values of the local efficiency in the two analyzed mental states in the beta rhythm at different thresholds. Red line represents values for the local efficiency in the subjects during the play task whereas black line refers to the local efficiency during the rest task. The short vertical lines denote the 95% confidence intervals for each threshold in the interval 0.1–0.9 with a step of 0.05. Within the sub-interval of threshold in which efficiency remains significantly different (p < 0.05) between the two mental states, 0.35 was used to illustrate the following results.
Figure 2. Average values of local efficiency at various values of the threshold during the mental tasks for the beta rhythm. The red line refers to the values in the play task and the black line represents the rest task. The short vertical lines denote the 95% confidence intervals for each threshold in the interval 0.1 to 0.9 with a step of 0.05.
Degrees and Degree Distributions
Figure 3 shows the degrees of the subjects in the frequency rhythms analyzed during different mental states. The delta, theta, alpha, and beta rhythms during the play task and rest task are shown in Figures 3A–D, respectively. A high degree-in denotes that a region is more likely to be influenced by a number of other brain regions whereas a high degree-out suggests a large quantity of functional targets. The results show that the degree-in follows an overall downward trend along the electrodes F3, F4, C3, C4, P3, P4, T3, T4, T5, and T6 in both mental states. However, notable flows outgoing (efferent) from the left temporal region (T3) are present in the play task, and flows outgoing from the right temporal region (T4) exist in the rest task. On the other hand, notable flows incoming (afferent) to left frontal region (F3) are present in the play task, and flows incoming to the right frontal region (F4) exist in the rest task.
Figure 3. Average degrees of the play and rest task in the four rhythms delta (A), theta (B), alpha (C), and beta (D). The red lines represent the play task and the black lines refer to the rest task. The stars represent the degree in and the squares refer to the degree out. The x-axes represent the network nodes and the y-axes refer to the average degree of all the subjects for each network node.
Figure 4 depicts the average degree distributions for the two mental tasks in various rhythms. The delta, theta, alpha, and beta rhythms during the play task and rest task are shown in Figures 4A–D, respectively. In order to compare the results, the normalized histogram values were taken into consideration. Right-skew tail of in-degree distributions in the play task in beta frequency indicates the presence of few nodes with a very high level of out-going connections.
Figure 4. Degree distributions of the play task and rest task in four rhythms delta (A), theta (B), alpha (C), and beta (D). Green bars represent the in-degree distribution and blue bars refer to the out-degree distribution in the play task. The purple bars represent the in-degree distribution in the rest task and the pink bars represent the out-degree distribution in the rest task. Histogram values were normalized to the number of the elements within the network (10 nodes).
Local Efficiency and Global Efficiency
To investigate the inter-individual variance of the results, the local and global efficiency were calculated for the networks estimated from each volunteer during the two mental states (play and rest) in the four rhythms (delta, theta, alpha, and beta). Furthermore, analysis of variance (ANOVA) was performed on the efficiency indexes to investigate the relationships between efficiency indexes, frequency bands, and different mental states.
Figure 5 shows the average local efficiency for various rhythms between two mental tasks. A statistically significant difference (p < 0.05) was noted between the two mental states. Red line represents the average values of the local efficiency during the play task, black line refers to the rest task, while vertical bars denote 95% confidence intervals.
Figure 5. Average values of local efficiency for various rhythms between two mental tasks. A statistically significant difference (p < 0.05) was noted between them. Red line represents the average values of the local efficiency during the play task, black line refers to the rest task, while vertical bars denote 95% confidence intervals.
ANOVA was further performed on the local efficiency. It showed no significant differences for the factor frequency band (F = 0.46, p = 0.7105) but significant differences for the factor group (F = 13.4, p = 0.0003). The interaction between group and band were found to be not significant (F = 0.77, p = 0.5131).
Figure 6 describes the average global efficiency for various rhythms between two mental tasks. There is a statistically significant difference (p < 0.05) between the two mental states.
Figure 6. Average values of global efficiency for various rhythms during the two mental tasks. A statistically significant difference (p < 0.05) was noted between them. Red line represents the average values of the local efficiency during the play task, black line refers to the rest task, while vertical bars denote 95% confidence intervals.
ANOVA was also performed on the global efficiency. The results showed no differences for the factor frequency band (F = 0.32, p = 0.8097) but significant differences for the factor group (F = 8.93, p = 0.0033). Moreover, the interaction between group and band were found to be not significant (F = 0.06, p = 0.9789).
Assessment of Electrodes Selection
Using the technique combining the LMD and SVM, we quantitatively selected the electrodes most related to different tasks. Results demonstrated that 10 electrodes (F3, F4, C3, C4, P3, P4, T3, T4, T5, and T6) had a predicted accuracy higher than 60%, therefore, we consider that the information from these electrodes can reflect the different mental states better than the other electrodes.
Generally, the first step in the construction of a structural or functional brain network is the definition of network nodes (Bullmore and Sporns, 2009). These can be defined either from EEG electrodes or cortical areas with a head model or anatomically selected regions from histological data, MRI or DTI (Bullmore and Sporns, 2009; Rubinov and Sporns, 2010). For our costumed EEG experiment, the data driven method combining LMD and SVM was used to quantitatively select the proper EEG electrodes. LMD proposed by Smith (2005) has been applied in a set of scalp EEG visual perception data. The idea that it can decompose signals into a series of components or functions characterized by the frequency seems interesting for EEGs, since it is widely accepted that different frequencies of EEGs (gamma, alpha, theta, and delta) are related to certain aspects of cognitive function.
SVM has been extensively used recently in classification in EEG (Hosseini and Khalilzadeh, 2010; Subasi and Gursoy, 2010; Temko et al., 2011; Li et al., 2013; Wang et al., 2014; Hosseini et al., 2015) and fMRI (LaConte et al., 2005; Mourão-Miranda et al., 2005; De Martino et al., 2007; Altmann et al., 2015). Given that the performance of SVM depends on the choice of gamma-g and cost-c, the GA was conducted on every dataset in a cross validation manner in order to select the optimal parameters for the SVM model. In this way, one can get the highest predicted accuracy as well as smallest c, since a higher c may contribute to “overfitting.” More precisely, the leave-one-out cross validation was used to ensure that the results from the experiments can be reproducible and reliable.
The technique combining LMD and SVM resulted in 10 chosen electrodes (F3, F4, C3, C4, P3, P4, T3, T4, T5, and T6). According to the layout of international 10–20 system, F3 and F4 electrodes are located in the frontal lobe. As is well know that, frontal lobe associates with attention, short-term memory, motivation, and planning (Stuss and Knight, 2002). Thus, the described electrodes make sense from the physiological point of view. C3 and C4 electrodes were located on the edge of the frontal lobe, so that these two electrodes should be included too. Furthermore, in our play task, the subject had to control the move of the “fish” in the game by a mouse, therefore, it was not unexpected that P3 and P4 electrodes showed a good classification: These electrodes were located in the parietal lobe, which is known for sensory information among various modalities, including the sense of touch in the somatosensory cortex (Fogassi and Luppino, 2005). Moreover, T3, T4, T5, and T6, which were located in the temporal or posterior temporal lobe, showed a high distinction among the two different mental states probably because in the play task, one had to open the eyes and concentrate on the game; on the other hand the temporal lobe is related to the process visual memories into interpretations.
The technique of combining LMD and SVM in this study was applied in the data from each subject and 10 electrodes were defined as the nodes of the brain network. In this way, a major source of subjectivity was eliminated and no prior knowledge was required to define the ROI.
Effective Connectivity in Different Mental States
In this work, PDC and graph theory were applied to study the effective brain networks in normal subjects during different mental tasks. Combining SVM with LMD, we identified the network nodes. The distinctive properties of the effective networks were adequately analyzed from the perspective of the degree, degree distribution, as well as local and global efficiency. As mentioned above, the combination of these techniques improves the capability of detecting relevant features of the effective brain networks.
From the analysis performed on the effective network from a healthy population in different mental states, we found that nodes in the frontal lobe (F3 and F4) in both mental states showed overall the highest in-degree. It suggests that there is a great influence from other regions of the brain. This means, the brain regions around these electrodes may act as “hub” for the in-flow of information in different mental states. Removal of this region from the estimated network is more likely to cause a collapsing of the whole system. We also demonstrated notable flows outgoing from the left temporal region (T3) in the play task, and flows ingoing from the right temporal region (T4) in the rest task. Moreover, notable flows ingoing to left the frontal region (F3) exists in the play task, and flows ingoing to the frontal region (F4) exist in the rest task. The phenomenon above may be explained by the division of the left and right hemispheres. More specifically, the left hemisphere takes more responsibility for logical thinking, such as planning a best way to escape from dangerous situation in the play task.
Although the average trend of degree distributions for the two mental tasks in various rhythms (Figure 4) seems a “scale free” property, it is difficult to formally assert this conclusion, because the size of the networks estimated in the study is not big enough to achieve a reliable degree distribution.
The higher values of local efficiency during the play task (Figure 5) result from the involvement related to the execution of “feeding frenzy,” which may imply that the networks have a larger level of the internal organization and fault tolerance (Pfurtscheller and Lopes da Silva, 1999; Sivan et al., 1999; Crucitti et al., 2003). The downward trend of local efficiency from the delta rhythm to the beta rhythm in the rest task indicates a more intense mental activity in the other rhythms than the beta rhythms. The properties above may suggest a more concentrated attentiveness as a local response in the play task. Furthermore, the higher global efficiency indexes in the play task suggest stronger information transfer ability compared with the rest mental states.
Another method toward evaluating causal relations in neural system is directed transfer function (DTF) (Kamiński et al., 2001). It also estimates the parameters for the MVAR models which derive from the multi-channel EEG signals and then transformed into the frequency domain. It has been successfully applied in several actual neurobiological recordings (De Vico Fallani et al., 2007; Liu et al., 2010; Nan et al., 2010; Dissanayaka et al., 2015). In particular, the combination of DTF with short window AMVAR method (Ding et al., 2000) is able to get the dynamic causal relations among cortical regions or EEG electrodes. However, in the simulation (not shown here), when the value of DTF between two interactions is not zero, it does not mean that there is a causal relationship between the two interactions (Baccalá and Sameshima, 2001). In other words, DTF is prone to be affected by possible alternative interactions or unpredictable factors. On the other hand, PDC which has been widely used recently (Schelter et al., 2006; Supp et al., 2007; Michels et al., 2013; Silchenko et al., 2013; Pester et al., 2015; Youssofzadeh et al., 2016), illustrates the relative strength of direct interactions and avoids the inverse operation, therefore reducing the computation time. With this consideration, PDC is more reliable and faster to quantify the causal interactions among multi-channel EEG signals in our study.
Limitation and Future Implications
We have to mention that the predicted accuracies from electrodes (Table 1) were overall not high. It may due to two factors. First, the experiments were carried out in a small chamber without electrical interference. We have mentioned that the experiment is part of the “personalized neurofeedback treatment and brain mechanisms project.” We therefore tried to set all the conditions similar to the home environment for the further comparison between a healthy population and patients. The second reason is that we placed the position of each electrode one by one manually according to layout of international 10–20 system; therefore, individual difference may exist between different experiment operators.
There are several implications of the present study. Firstly, the participants in our study were university students. Previous studies showed that the executive function is related to the neurobiological development of the brain (De Luca and Leventer, 2008). Furthermore, Liu et al. (Liu et al., 2015) has reported the properties of attention-related functional networks between ADHD and normal control children with EEG. Therefore, it would be interesting to find how the efficiency indexes are for ADHD Children from the information flow perspective. However, this is not the scope of the present study.
Secondly, the threshold we used in our study was based on the paired t-test according to the local efficiency during different mental states in the beta frequency, which was considered to be the most prominent band in our experiment. This threshold could destroy the cohesion of the network. Since there is currently no consent on how to choose an optimal threshold, we would like to skip this step by analysing weighted graph directly in the future study.
What is more, the fact that the local efficiency is higher in the beta frequency and lower in the theta frequency during play task while the global efficiency is higher in the theta frequency and lower in the beta frequency in the rest task suggests that these network properties could be indexes to classify different mental states in a healthy population. Similarly, Ahmadlou et al. (2012) showed that the short path lengths and clustering coefficient can distinguish the ADHD brain from non-ADHD from the delta EEG sub-band. Petti et al. (Petti et al., 2013) used the small-worldness to classify the populations between young and mid-aged adults with EEG.
This study aimed to investigate the effective networks in a healthy population in different mental states without the prior knowledge for the region of interest for the experiment design. We have demonstrated that the graph theory combined with the PDC can describe the local and the global properties in the effective connections estimated from the EEG signals. The network measures, such as degrees together with their distributions, local and global efficiency, were applied to investigate the properties of the effective network between the different mental tasks.
The present study can be also regarded as the step of investigation of executive function of attention, to which eventually aim at investigation from ADHD population in the future. It is expected that our study will have useful and valuable implications for preclinical studies.
Conceived and designed the experiments: DH, AR, and QL. Performed the experiments: YZ and ZY. Analyzed the data: DH and SJ. Contributed reagents/materials/analysis tools: KV. Wrote the paper: DH and LH.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by the National Natural Science Foundation of China under grant NOs.81071221 and 31271063; the Fundamental Science Research Funds for the Central Universities under grant NO. NSIY131409. The authors thank Dr. E. Farrher for proofreading the manuscript.
Altmann, A., Schröter, M. S., Spoormaker, V. I., Kiem, S. A., Jordan, D., Ilg, R., et al. (2015). Validation of non-REM sleep stage decoding from resting state fMRI using linear support vector machines. Neuroimage 125, 544–555. doi: 10.1016/j.neuroimage.2015.09.072
Anderson, S. W., Damasio, H., Jones, R. D., and Tranel, D. (1991). Wisconsin card sorting test performance as a measure of frontal lobe damage. J. Clin. Exp. Neuropsychol. 13, 909–922. doi: 10.1080/01688639108405107
Axelrod, B. N., Goldman, R. S., Heaton, R. K., Curtiss, G., Thompson, L. L., Chelune, G. J., et al. (1996). Discriminability of the Wisconsin Card sorting test using the standardization sample. J. Clin. Exp. Neuropsychol. 18, 338–342. doi: 10.1080/01688639608408991
Blinowska, K. J., Kuś, R., and Kamiński, M. (2004). Granger causality and information flow in multivariate processes. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70:050902. doi: 10.1103/PhysRevE.70.050902
Crucitti, P., Latora, V., Marchiori, M., and Rapisarda, A. (2003). Efficiency of scale-free networks: error and attack tolerance. Phys. A Stat. Mech. Appl. 320, 622–642. doi: 10.1016/S0378-4371(02)01545-5
Cui, J., Xu, L., Bressler, S. L., Ding, M., and Liang, H. (2008). BSMART: a Matlab/C toolbox for analysis of multichannel neural time series. Neural Netw. 21, 1094–1104. doi: 10.1016/j.neunet.2008.05.007
Delorme, A., and Makeig, S. (2004). EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 134, 9–21. doi: 10.1016/j.jneumeth.2003.10.009
De Luca, C. R., and Leventer, R. J. (2008). “Developmental trajectories of executive functions across the lifespan,” in Executive Functions and the Frontal Lobes: A Lifespan Perspective, eds V. Anderson, R. Jacobs, and P. J. Anderson (New York, NY: Taylor & Francis), 3–21.
De Martino, F., Gentile, F., Esposito, F., Balsi, M., Di Salle, F., Goebel, R., et al. (2007). Classification of fMRI independent components using IC-fingerprints and support vector machine classifiers. Neuroimage 34, 177–194. doi: 10.1016/j.neuroimage.2006.08.041
De Vico Fallani, F., Astolfi, L., Cincotti, F., Mattia, D., Marciani, M. G., Salinari, S., et al. (2007). Cortical functional connectivity networks in normal and spinal cord injured patients: evaluation by graph analysis. Hum. Brain Mapp. 28, 1334–1346. doi: 10.1002/hbm.20353
De Vico Fallani, F., Nicosia, V., Sinatra, R., Astolfi, L., Cincotti, F., Mattia, D., et al. (2010). Defecting or not defecting: how to “read” human behavior during cooperative games by EEG measurements. PLoS ONE 5:e14187. doi: 10.1371/journal.pone.0014187
Ding, M., Bressler, S. L., Yang, W., and Liang, H. (2000). Short-window spectral analysis of cortical event-related potentials by adaptive multivariate autoregressive modeling: data preprocessing, model validation, and variability assessment. Biol. Cybern. 83, 35–45. doi: 10.1007/s004229900137
Dissanayaka, C., Ben-Simon, E., Gruberger, M., Maron-Katz, A., Sharon, H., Hendler, T., et al. (2015). Comparison between human awake, meditation and drowsiness EEG activities based on directed transfer function and MVDR coherence methods. Med. Biol. Eng. Comput. 53, 599–607. doi: 10.1007/s11517-015-1272-0
Fan, J., McCandliss, B. D., Sommer, T., Raz, A., and Posner, M. I. (2002). Testing the efficiency and independence of attentional networks. J. Cogn. Neurosci. 14, 340–347. doi: 10.1162/089892902317361886
Fox, M. D., Snyder, A. Z., Vincent, J. L., Corbetta, M., Van Essen, D. C., and Raichle, M. E. (2005). The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc. Natl. Acad. Sci. U.S.A. 102, 9673–9678. doi: 10.1073/pnas.0504136102
Gast, H., Müller, M., Rummel, C., Roth, C., Mathis, J., Schindler, K., et al. (2014). Epileptic seizures as condensed sleep: an analysis of network dynamics from electroencephalogram signals. J. Sleep Res. 23, 268–273. doi: 10.1111/jsr.12113
Gupta, R., and Falk, T. H. (2015). “Affective state characterization based on electroencephalography graph-theoretic features,” in 2015 7th International IEEE/EMBS Conference on Neural Engineering (NER) (Montpellier: IEEE), 577–580. doi: 10.1109/ner.2015.7146688
Hosseini, S. A., and Khalilzadeh, M. A. (2010). “Emotional stress recognition system using EEG and psychophysiological signals: using new labelling process of EEG signals in emotional stress state,” in 2010 International Conference on Biomedical Engineering and Computer Science (ICBECS) (Yantai), 1–6. doi: 10.1109/ICBECS.2010.5462520
Hosseini, S. A., Khalilzadeh, M. A., Naghibi-Sistani, M. B., and Homam, S. M. (2015). Emotional stress recognition using a new fusion link between electroencephalogram and peripheral signals. Iran. J. Neurol. 14, 142–151.
Kamiński, M., Ding, M., Truccolo, W. A., and Bressler, S. L. (2001). Evaluating causal relations in neural systems: granger causality, directed transfer function and statistical assessment of significance. Biol. Cybern. 85, 145–157. doi: 10.1007/s004220000235
Klados, M. A., Kanatsouli, K., Antoniou, I., Babiloni, F., Tsirka, V., Bamidis, P. D., et al. (2013). A Graph theoretical approach to study the organization of the cortical networks during different mathematical tasks. PLoS ONE 8:e71800. doi: 10.1371/journal.pone.0071800
Klados, M. A., Lithari, C., Antoniou, I., Semertzidou, A., Bratsas, C., Micheloyannis, S., et al. (2012). “Towards a graph theoretical approach to study gender lateralization effect in mathematical thinking,” in 2012 IEEE 12th International Conference on Bioinformatics & Bioengineering (BIBE) (Larnaca: IEEE), 666–670. doi: 10.1109/bibe.2012.6399746
Knyazeva, M. G., Jalili, M., Brioschi, A., Bourquin, I., Fornari, E., Hasler, M., et al. (2010). Topography of EEG multivariate phase synchronization in early Alzheimer's disease. Neurobiol. Aging 31, 1132–1144. doi: 10.1016/j.neurobiolaging.2008.07.019
LaConte, S., Strother, S., Cherkassky, V., Anderson, J., and Hu, X. (2005). Support vector machines for temporal classification of block design fMRI data. Neuroimage 26, 317–329. doi: 10.1016/j.neuroimage.2005.01.048
Lin, P., Sun, J., Yu, G., Wu, Y., Yang, Y., Liang, M., et al. (2014). Global and local brain network reorganization in attention-deficit/hyperactivity disorder. Brain Imaging Behav. 8, 558–569. doi: 10.1007/s11682-013-9279-3
Liu, J. P., Zhang, C., and Zheng, C. X. (2010). Estimation of the cortical functional connectivity by directed transfer function during mental fatigue. Appl. Ergon. 42, 114–121. doi: 10.1016/j.apergo.2010.05.008
Liu, T., Chen, Y., Lin, P., and Wang, J. (2015). Small-world brain functional networks in children with attention-deficit/hyperactivity disorder revealed by EEG synchrony. Clin. EEG Neurosci. 46, 183–191. doi: 10.1177/1550059414523959
Michels, L., Muthuraman, M., Lüchinger, R., Martin, E., Anwar, A. R., Raethjen, J., et al. (2013). Developmental changes of functional and directed resting-state connectivities associated with neuronal oscillations in EEG. Neuroimage 81, 231–242. doi: 10.1016/j.neuroimage.2013.04.030
Mourão-Miranda, J., Bokde, A. L. W., Born, C., Hampel, H., and Stetter, M. (2005). Classifying brain states and determining the discriminating activation patterns: support vector machine on functional MRI data. Neuroimage 28, 980–995. doi: 10.1016/j.neuroimage.2005.06.070
Nan, Y., Wang, J., Xue, S. A., Sheng, H., Jiao, Y., and Wang, J. (2010). Analysis of propagation of multi-channel EEG in the test of sustained attention. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2010, 1666–1669. doi: 10.1109/iembs.2010.5626692
Pester, B., Ligges, C., Leistritz, L., Witte, H., and Schiecke, K. (2015). Advanced insights into functional brain connectivity by combining tensor decomposition and partial directed coherence. PLoS One 10:e0129293. doi: 10.1371/journal.pone.0129293
Petti, M., Toppi, J., Pichiorri, F., Cincotti, F., Salinari, S., Babiloni, F., et al. (2013). Aged-related changes in brain activity classification with respect to age by means of graph indexes. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2013, 4350–4353. doi: 10.1109/EMBC.2013.6610509
Pfurtscheller, G., and Lopes da Silva, F. H. (1999). Event-related EEG/MEG synchronization and desynchronization: basic principles. Clin. Neurophysiol. 110, 1842–1857. doi: 10.1016/S1388-2457(99)00141-8
Piek, J. P., Dyck, M. J., Nieman, A., Anderson, M., Hay, D., Smith, L. M., et al. (2004). The relationship between motor coordination, executive functioning and attention in school aged children. Arch. Clin. Neuropsychol. 19, 1063–1076. doi: 10.1016/j.acn.2003.12.007
Pritchard, W. S., Laurienti, P. J., Burdette, J. H., and Hayasaka, S. (2014). Functional brain networks formed using cross-sample entropy are scale free. Brain Connect. 4, 454–464. doi: 10.1089/brain.2013.0217
Schelter, B., Winterhalder, M., Eichler, M., Peifer, M., Hellwig, B., Guschlbauer, B., et al. (2006). Testing for directed influences among neural signals using partial directed coherence. J. Neurosci. Methods 152, 210–219. doi: 10.1016/j.jneumeth.2005.09.001
Silchenko, A. N., Adamchic, I., Hauptmann, C., and Tass, P. A. (2013). Impact of acoustic coordinated reset neuromodulation on effective connectivity in a neural network of phantom sound. Neuroimage 77, 133–147. doi: 10.1016/j.neuroimage.2013.03.013
Stuss, D. (2000). Wisconsin Card Sorting Test performance in patients with focal frontal and posterior brain damage: effects of lesion location and test structure on separable cognitive processes. Neuropsychologia 38, 388–402. doi: 10.1016/S0028-3932(99)00093-7
Supp, G. G., Schlögl, A., Trujillo-Barreto, N., Müller, M. M., and Gruber, T. (2007). Directed cortical information flow during human object recognition: analyzing induced EEG gamma-band responses in brain's source space. PLoS One 2:e684. doi: 10.1371/journal.pone.0000684
Temko, A., Thomas, E., Marnane, W., Lightbody, G., and Boylan, G. (2011). EEG-based neonatal seizure detection with Support Vector Machines. Clin. Neurophysiol. 122, 464–473. doi: 10.1016/j.clinph.2010.06.034
Van Dellen, E., Douw, L., Baayen, J. C., Heimans, J. J., Ponten, S. C., Vandertop, W. P., et al. (2009). Long-term effects of temporal lobe epilepsy on local neural networks: a graph theoretical analysis of corticography recordings. PLoS One 4:e8081. doi: 10.1371/journal.pone.0008081
Welsh, M. C. (2002). “Developmental and clinical variations in executive functions,” in Developmental Variations in Learning: Applications to Social, Executive Function, Language, and Reading Skills, eds D. L. Molfese and V. J. Molfese (Hove: Psychology Press), 139–185.
Yan, N., Wang, J., Wei, N., and Zong, L. (2007). Feature Exaction and Classification of Attention Related Electroencephalographic Signals Based on Sample Entropy. Xi'an: Hsi-An Chiao Tung Ta Hsueh/Journal Xi'an Jiaotong University.
Youssofzadeh, V., Prasad, G., Naeem, M., and Wong-Lin, K. (2016). Temporal information of directed causal connectivity in multi-trial ERP data using partial granger causality. Neuroinformatics 14, 99–120. doi: 10.1007/s12021-015-9281-6
Zhang, R., Yao, D., Valdés-Sosa, P. A., Li, F., Li, P., Zhang, T., et al. (2015). Efficient resting-state EEG network facilitates motor imagery performance. J. Neural Eng. 12:066024. doi: 10.1088/1741-2560/12/6/066024
Keywords: electroencephalogram (EEG), partial directed coherence, graph theory, effective networks, mental tasks
Citation: Huang D, Ren A, Shang J, Lei Q, Zhang Y, Yin Z, Li J, von Deneen KM and Huang L (2016) Combining Partial Directed Coherence and Graph Theory to Analyse Effective Brain Networks of Different Mental Tasks. Front. Hum. Neurosci. 10:235. doi: 10.3389/fnhum.2016.00235
Received: 04 October 2015; Accepted: 04 May 2016;
Published: 20 May 2016.
Edited by:Yong He, Beijing Normal University, China
Reviewed by:Lester Melie-Garcia, Lausanne University Hospital, Switzerland
Manousos A. Klados, Max Planck Institute for Human Cognitive & Brain Sciences, Germany
Wei Liao, University of Electronic Science and Technology of China, China
Copyright © 2016 Huang, Ren, Shang, Lei, Zhang, Yin, Li, von Deneen and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Liyu Huang, firstname.lastname@example.org
†These authors have contributed equally to this work.