- 1Rehabilitation Medical Department, Tianjin Union Medical Centre, Tianjin, China
- 2Department of Health and Exercise Science, Tianjin University of Sport, Tianjin, China
- 3Rehabilitation Medical Department, Tianjin University of Traditional Chinese Medicine, Tianjin, China
- 4Department of Biomedical Engineering, College of Precision Instrument and Opto-Electronics Engineering, Tianjin University, Tianjin, China
Post-stroke depression (PSD) is the most common stroke-related emotional disorder, and it severely affects the recovery process. However, more than half cases are not correctly diagnosed. This study was designed to develop a new method to assess PSD using EEG signal to analyze the specificity of PSD patients' brain network. We have 107 subjects attended in this study (72 stabilized stroke survivors and 35 non-depressed healthy subjects). A Hamilton Depression Rating Scale (HDRS) score was determined for all subjects before EEG data collection. According to HDRS score, the 72 patients were divided into 3 groups: post-stroke non-depression (PSND), post-stroke mild depression (PSMD) and post-stroke depression (PSD). Mutual information (MI)-based graph theory was used to analyze brain network connectivity. Statistical analysis of brain network characteristics was made with a threshold of 10–30% of the strongest MIs. The results showed significant weakened interhemispheric connections and lower clustering coefficient in post-stroke depressed patients compared to those in healthy controls. Stroke patients showed a decreasing trend in the connection between the parietal-occipital and the frontal area as the severity of the depression increased. PSD subjects showed abnormal brain network connectivity and network features based on EEG, suggesting that MI-based brain network may have the potential to assess the severity of depression post stroke.
Introduction
Post-stroke depression (PSD) is among the most frequent neuropsychiatric consequences of cerebral ischemia (Cojocaru et al., 2013). PSD is an abnormal negative emotional response caused by loss, disappointment or failure. PSD has a significant negative impact on the rehabilitation of stroke (Ghose et al., 2005), thus seriously affecting the patient's future quality of life (Bays, 2001; Ayerbe et al., 2013; Chen et al., 2014) and delaying or even hindering the process for rehabilitation and return to society. Approximately one-third of stroke patients have aphasia (Berthier, 2005; Engelter et al., 2006), and approximately 70% will have cognitive impairment (Nys et al., 2007). Aphasia and cognitive impairment make it difficult to obtain the changes in patients' emotions and interests, which pose a great challenge for the diagnosis of PSD. There are few guidelines for the assessment, treatment and prevention of PSD (Babkair, 2017), and more than half cases are not correctly diagnosed.
Depression was thought to be the result of a dysregulation in the ability of brain cells to communicate with each other (Cai et al., 2013). Researchers have found abnormalities in the transmission of excitatory signals between cells in depression. Restoring normal brain communication is one mechanism underlying the successful function of antidepressant drugs such as serotonin, which is a key factor in depression remission (Cai et al., 2013). Disrupted network connectivity has been found in some core major depressive disorder (MDD) networks (Brakowski et al., 2017). Previous findings in geriatric depression have also strongly suggested “brain network dysfunction” as the best explanatory model for understanding the biological mechanism of depression (Drevets et al., 2008). All of the possible etiologies of late-life depression result in different depressive symptoms by disturbing the dynamics and functions of different brain networks (Tadayonnejad and Ajilore, 2014). Impairment of the affective regulatory pathway has been suggested as a possible pathogenic factor related to vascular disease according to previous studies (Alexopoulos et al., 1997a,b). We suggest that PSD patients' abnormal connectivity among brain areas could be driving this pathogenesis, which may appear as “disconnection” symptoms.
Functional connectivity in the human brain can be represented as a network using electroencephalography (EEG) signals (Rathee et al., 2017). One of the functional connectivity measures for analyzing EEG is Mutual information (MI) which is a non-directional connectivity measure. It enables the estimation of both linear and non-linear statistical dependencies between time series and can be used to detect functional coupling (Wang et al., 2009). Because neural dynamics almost certainly includes many highly nonlinear processes, MI analysis may be helpful in understanding and quantifying the nonlinear transmission of information within the brain (Jeong et al., 2001). Abnormal cortical connections using MI have been found in nervous system diseases, such as Alzheimer's, schizophrenia and Parkinson's (Coronel et al., 2017; Yin et al., 2017).
Graph theory has played an integral role in recent efforts to understand the function of complex systems including brain networks. Importantly, graph-based representations of brain networks can quantitatively describe the connectivity of different brain regions. It has been applied to understand brain networks and emerged as a powerful analytic tool for brain connectivity. Using this method, many researchers have studied the structural and functional networks of the brain and the network anomalies caused by neuropsychiatric disorders (Schreiber, 2000; Bernhardt et al., 2013; Rathee et al., 2017). In brain networks, different connections represent different paths of information transfer. This study aimed to analyze the features of MI-based undirected and weighted brain network to explore the abnormal brain connectivity of the stroke patients with different degrees of depression.
Materials and Methods
Participants
This study was performed in the Department of Rehabilitation, Tianjin Union Medical Center, Tianjin, China. All participants were right-handed and native speakers of Mandarin Chinese. The hospital ethics committee approved the study. All participants were informed of the aims and protocols of the experiments.
This study involved 35 healthy controls (HC) and 72 stroke patients. The HC group had no history of neurological or psychiatric disease. All patients were divided into three groups based on their Hamilton Depression Rating Scale (HDRS) score. The patients in the post-stroke non-depression group (PSND), post-stroke mild-depression group (PSMD), and post-stroke depression group (PSD) have HDRS scores of ≤ 5, 6–20, and >20, respectively. Other demographic and general subject characteristics are listed in Table 1.
EEG Recording and Preprocessing
The subjects were seated in a resting state with their eyes closed for 5 min in a quiet environment. The EEG was recorded at 16 scalp loci (Fp1, Fp2, F3, F4, F7, F8, C3, C4, T3, T4, P3, P4, O1, O2, T5, and T6) in compliance with the international 10–20 system using a NicoletOne digital video electroencephalograph made by US. The skin resistance at each site was <10 kΩ. EEG data were collected for 300 s at a rate of 250 Hz. Data containing artifacts were removed in an off line analysis. We also used independent component analysis (ICA) to identify and remove residual ocular activity (Fanciullacci et al., 2017). The EEG signals were re-referenced to the bilateral mastoid electrodes (A1 and A2), and removed each channels baseline from continuous EEG data by using the routine pop_rmbase (EEGLAB). Then a Hamming windowed sinc FIR filter was used to filter the data with a bandwidth of 0.1–100 Hz by using the routine pop_eegfiltnew (EEGLAB).
As previous studies have proved that the infinity reference was proper for EEG network analysis (Qin et al., 2010), we changed linked earlobes to infinity reference using a reference electrode standardization technique (REST) (Dong et al., 2017; Yao, 2017). REST is used for the approximate standardization of the reference of scalp EEG recordings to a point at infinity that, being far from all possible neural sources, acts like a neutral virtual reference(Marzetti et al., 2007). Numerous studies have shown that REST is the most accurate reference method for brain network analysis (Yao, 2001; Qin et al., 2010). A REST toolbox which developed by Dong et al. (2017) were used in this study.
Multivariate Causal Analysis of Data
In information theory, MI is a measure of the statistical dependence between two random variables (Ince et al., 2017). The average amount of information obtained from any observation of X = {xi} is the entropy H of a system:
where PX (xi) is the probability that an isolated measurement will find the system in the ith element of the bin. We evaluated these probabilities PX(xi) by constructing a histogram (from 1,250 data points) of the variations of the measurement xi.
Before any measurement of X, this information is called uncertainty. Under the condition Y = yj, H(X) has to be replaced by the conditional uncertainty on X
Where PXY(xi, yj) is the joint probability density for the measurements of X and Y that produce the values X and Y. H(X|Y = yj) indicates the amount of uncertainty in a measurement of x, given that y has been measured and found to be yj. From this, we get the mean conditional uncertainty on X over yj, under the condition that Y is known
where
So we define the MI as the amount by which a measurement of Y reduces the uncertainty of X. The MI is as follows: MIXY = H(X)−H(X|Y) = H(X)+H(Y)−H(X, Y) = MIYXwhich can be rewritten as:
MI has the maximum value when the two time series are completely the same. If one system is completely independent of the other, the MI is zero (Na et al., 2002). The principal difficulty in calculating the MI from experimental data is estimating PXY(x, y) from histograms, selecting different sampling bins has a great influence on the accuracy of MI (Jeong et al., 2001). In this study, we took logarithm with base e, and 11 bins were adopted for 1,250 samples, which can provide a stable estimate.
In this study, EEG data were segmented into 5-s epochs (1,250 data points), and a total of 60 epochs for each channel were analyzed using the MI. 60 MI values can be obtained between any two channels, and the mean index value of the 60 MIs were performed as the final MI indices. All routines above were implemented in MATLAB (MathWorks, Inc.).
Graphical Description of the Network
Graph theory has proven very useful in statistics as a way to describe the dependent relations between random variables (Salvador et al., 2005). In graph theory, a network is reduced to an abstract description as a set of nodes connected by edges (or lines) (Bassett and Bullmore, 2009). The edges can be directed or undirected and weighted or unweighted.
The nodes and edges of a brain graph can be empirically defined in many ways. In this study, we used 16 leads as nodes and constructed the cortical undirected network graph by using the calculated MI as the edge of the network. By using the topological properties of networks, we analyzed the characteristics of brain networks in different subjects and then explored the abnormal connectivity of the brain in patients with depression after stroke. The setup process for the brain network is shown in Figure 1.
Figure 1. Schematic diagram of a brain network based on mutual information (MI) technology. (A) nodes location; (B) acquisition of 16-lead EEG data; (C) calculation of MI; (D) formation of a MI matrix; (E) adjacency matrix with a certain threshold; (F) the brain network diagram was formed according to the relationship among the 16 leads.
Topological properties of a brain network can be described using some graph measures based on Graph theory, such as clustering coefficient and betweenness centrality. Clustering coefficient is one key topological metric which quantifies degree of collectivization of one network. The clustering coefficient of one node measures the connecting size of its adjacent edges. The calculation formula for clustering coefficient Ci of node i is shown
Where, ki is the number of all adjacent nodes of node i, ei is the number of connected edges between all neighboring nodes of node i. One node had value 0, while which only has a neighbor or none. The mean clustering coefficients of all the nodes represent the network's coefficient. The betweenness centrality is used to describe the role and status of one node to the network. Higher betweenness centrality indicates more important status and the corresponding node is a core node for the network. The caculation formula for betweeness centrality is shown
σjk(i) is the number of shortest path from node j to node k, which passing node i. In this study, clustering coefficients and betweenness centrality were calculated by binary MI matrices (elements above the threshold were defined as 1, otherwise defined as 0) at each threshold.
Simulation of MI-Based Brain Network
Using MI to assess statistical dependence between two EEG signals, there can be contamination of spurious connectivity caused by volume conduction. In order to solve this problem, we used a surrogate data method to conduct a simulation study. We generated a dataset which has the same structure with our EEG data using Matlab code provided by Stefan Haufe et al. (Fonov et al., 2009, 2011; Haufe et al., 2013). In this dataset a linear time-lagged information flow from the left hemisphere (brain area below C3) to the right hemisphere (area below C4) is simulated by means of a bivariate AR model. This flow is to be detected as the only true time-lagged interaction happening in the data. We have established the MI brain network of this dataset, the result is shown in Figure 2. It can be seen from Figure 2 that MI can better reflect the true connection between the corresponding brain regions and suppress spurious connectivity. This is basically consistent with the connectivity between simulated EEG sensor measurements estimated by phase-slope index (PSI) in Stefan Haufe et al. (2013). The difference between the two methods is that MI has no directionality, and PSI can reflect the direction of information flow.
Figure 2. The brain network of the simulation dataset. Information flow from the left (below C3) to the right (below C4) source is modeled by means of a bivariate AR model.
Results
As shown in Table 1, four groups showed no significant difference in other demographic and clinical features except for HDRS.
There are 16 channels' time series of 300 s duration for each subject. These time series were analyzed in sequential windows of 5 s duration, yielding 16 time series with a length of 60 epochs. For each subject, this approach yielded 120 unique MIs (from the 16 × 16 MI matrix removing diagonal and symmetric data).
Figure 3 shows the rank-ordered average MIs for unthresholded MI matrices of four groups. We can find that more than 70% of MIs are between 0.05 and 0.2, they capture only a small amount of the common variance (the square of the MI < 4%) in the underlying dynamics, and also the difference between the four curves in this range is not significant in the Figure 3. According to the study of Rubinov et al. (2009), the 10–30% of the strongest MIs are more likely to reflect the underlying network architecture. Selected a certain range is also more convenient for us to find patterns in complex brain networks, so the following analysis mainly focused on 10–30% of the strongest MIs. There is no significant difference between four groups in the Figure 3. MIs in healthy people are on average slightly greater than in the other three groups for a range of rank-ordered means.
Figure 3. Edges from mean unthresholded MI matrices of four subject groups rank-ordered by MI values. There are a total of 120 unique correlations in each unthresholded MI matrix. The strongest 10–30% of these MIs are considered for subsequent graph analysis. MIs in healthy people are on average slightly greater than in the other three groups for a range of rank-ordered means.
Figure 4 shows the average MI matrices for four groups, which were thresholded such that 20% of the strongest edges are presented. The MI matrices had the same number of elements after thresholding. The white matrix elements represent functional connectivity. The connection between the parietal-occipital area and the frontal area shows a decreasing trend as the severity of the disease increases (white squares in the Figure 4).
Figure 4. The average MI matrices of four subject groups thresholded such that only 20% of the strongest weights are preserved. The white matrix elements represent functional connectivity. The key difference areas are marked with boxes.
Figure 5 shows the brain networks based on the average MI matrices of four groups thresholded such that only 20% of the strongest edges are preserved. Different colors represent the size of the betweenness centrality of the nodes, that is, the importance of each node in the network. The connection between the left and right brain is weakened as the degree of depression increases. And the internal connections of each hemisphere have been enhanced correspondingly.
Figure 5. The brain networks based on average MI matrices of four subject groups at threshold 20% of the strongest edges. Different colors indicate the betweenness centrality of the each node.
The nodes color shows that the core nodes of PSD were more scattered than the other three groups. This may affect the degree of clustering between nodes, and the subsequent analysis proves this by the statistical results of clustering coefficients. The subsequent statistical analysis results showed that there was no significant difference in the betweenness centrality of each node among the four groups at almost all thresholds.
According to the brain network characteristics of Figure 5, we performed statistical analysis of the relevant topological properties. Post-stroke depressed subjects showed weaker connections between the left and right hemisphere. Table 2 shows the statistical significance of the edge numbers between left and right cerebral hemispheres across 10–30% thresholds, assessed at each threshold by One-Way ANOVA. Post-hoc group comparisons were performed using Least Significant Difference (LSD) or Tamhane's T2 (IBM SPSS Statistics 19), according to whether the variance meets the condition of homogeneity. The significant results (p < 0.05) are in bold. If the p values are all above 0.05 at all thresholds between two groups will not be shown here. As can be seen from the table, there are some significant differences between the HC group and the PSMD and the PSD group at about half number of the thresholds. In particular, there is a clear difference between HC and PSD at a large continuous threshold range (19–29%). We noticed no significant difference between the HC group and the PSND group, but there are differences between PSND and PSD at several thresholds.
Table 3 shows clustering coefficients in four groups with thresholds of 10–30% of strongest edges. The difference in clustering is significant at 11 of the 21 thresholds between HC and PSD, most assemble at higher thresholds, and there is no significant difference between any other two groups. The clustering of HC is higher than that of PSD in Table 3. It indicates that healthy people's EEG signals have a higher degree of clustering.
Discussion
In this study, we examined the brain network performance in post-stroke depressed patients using the EEG-MI. We found that stroke patients with different degrees of depression showed different connection features. These features may be helpful in the diagnosis of PSD. Our results showed significant weakened connections between the left and right cerebral hemispheres in stroke patients compared to those in healthy controls, and this feature is more obvious with the deepening of the degree of depression. This suggests that depression affects the information communication between the left and right hemispheres in stroke patients. Among the stroke patients, the core nodes of PSD were more scattered than the other three groups. The connections between the parietal-occipital area and the frontal area showed a decreasing trend as the severity of the depression increases. Post-stroke depressed patients have a lower clustering coefficient than healthy subjects, with a significant difference at one-half thresholds.
The basal ganglia proved to play key roles in cortical and subcortical connected circuits, including the frontal, premotor and motor networks (Draganski et al., 2008; Thomas, 2009; Lao et al., 2016). This area may receive multiple cortical inputs in the presence of oscillatory activity and produce a high frequency drive back to the cerebral cortex, especially the supplementary motor area (Williams et al., 2002). Dysfunction of the frontal-parietal-occipital network in stroke patients may result from an organic lesion of the basal ganglia.
For depressed patients following stroke, the interhemispheric interaction was found to be highly disturbed in this study. Yamada et al. (1995) found that depressed patients showed lower frontal interhemispheric coherences than normal controls in each EEG band, and EEG power and coherence in presenile and senile depression. Wei et al. (2010) get similar findings with the above research. Furthermore, they found the inter-hemispheric coherence was correlated with some emotional processing. A decreased interhemispheric modulation was found in patients with major depression (Bajwa et al., 2008; Wu, 2014), which is consistent with our findings. Slow interhemispheric switching mechanisms in mood disorders may explain the weakened hemispheric information flow in PSD patients.
The frontal lobe plays a regulatory role in emotional cognition, the connection between the parietal-occipital and the frontal was decreased in depression in this study. Previous studies have reported aberrant EEG performance, such as increased slow activity in the frontal areas (Grin-Yatsenko et al., 2009, 2010), in depressed patients. Depressed older adults were found to have decreased frontal and parietal activation during some working memory tasks (Dumas and Newhouse, 2015). Weakened prefrontal and frontal connections may suggest decreased activation of the cortico-limbic circuit, which is related to symptoms such anhedonia or blunted affect (Fingelkurts and Fingelkurts, 2006). Some studies found that local information flow in the frontal-parietal-occipital network was related to the level of sedation (Rathee et al., 2016). For most stroke patients, the main symptoms of depression are decreased interest and retardation, which may cause the performance in the frontal-parietal-occipital network to become similar to that with sedation.
Post-stroke depressed patients exhibit lower clustering coefficients and more diffuse distribution of core nodes. The hypothesis of nerve loop connectivity injury has been used to explain the incidence of depression in some studies. Specifically, the pathogenesis of depression has certain neuroanatomical mechanisms. Damage to certain brain-related areas disrupts the neural pathway of emotional regulation, resulting in depressive episodes (Greicius et al., 2007; Alexopoulos et al., 2008). Previous studies found abnormal connectivity of neural circuits in depressed subjects. Studies also found that antidepressant drugs can restore this connection, which identified the relationship between the incidence of depression and nerve connection disorders (Cai et al., 2013; Tadayonnejad et al., 2014; Gudayol-Ferre et al., 2015). We suggest that abnormal communication in emotion-related brain areas results in disconnection in PSD subjects, and this phenomenon is also related to the damaged “core node.”
Previous studies found dopamine-dependent changes in the functional connectivity between the basal ganglia and cerebral cortex (Williams et al., 2002). As depressive disorders were considered a syndrome of cortical-subcortical dysrhythmia (Fingelkurts and Fingelkurts, 2015), a basal ganglia lesion should disrupt the normal cortical-subcortical neural pathway, which regulates emotions. Our results support the conclusion that post-stroke depressed subjects demonstrated abnormal brain network connectivity and that network features determined based on EEG may be utilized as reliable biomarkers for the effective assessment of PSD in the future.
There are some limitations of the present study. Only 16 EEG channels were used in the present study, which limited the network nodes. We plan to collect 64-channel EEGs in the future to obtain a more precise network. Another limitation of the study was the different locations of hemispheric lesions in participants. This study contained both left and right hemispheric lesioned patients, which may confound the current results. As the left and right hemispheres have different roles in emotional processing, a depressed mood following different hemispheric lesions may result from different brain disconnections. In future studies, we plan to investigate differences in the brain network in post-stroke depressed subjects with left and right hemispheric lesions. There can be contamination of spurious connectivity caused by volume conduction using MI algorithm. We will try more methods such as phase lag index (Stam et al., 2010) and imaginary part of coherency (Nolte et al., 2004) to cope with this limitation in later studies.
Ethics Statement
This study was carried out in accordance with the recommendations of Tianjin Union Medical Center committee with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the Tianjin Union Medical Center committee.
Author Contributions
CS did the EEG analysis work and write the article. FY made all the figures and tables. CW worked for the depression diagnosis and severity assessment of poststroke patients. ZW did the work of statistical analyses of the data. YZ recruited all the subjects and made the selection for the study. DM verified the results of the article. JD designed the research plan and offered the electroencephalograph acquisition equipment.
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.
Acknowledgments
This study was supported by National Natural Science Foundation of China (81630051), Scientific and Technological Projects of Tianjin Municipal Commission of Health and Family Planning (14KG107), Research projects in key areas of Tianjin traditional Chinese medicine (2017009), and Tianjin Union Medical Centre project (2017YJ011), Tianjin, China.
References
Alexopoulos, G. S., Meyers, B. S., Young, R. C., Campbell, S., Silbersweig, D., and Charlson, M. (1997a). “Vascular depression” hypothesis. Arch. Gen. Psychiatry 54, 915–922.
Alexopoulos, G. S., Meyers, B. S., Young, R. C., Kakuma, T., Silbersweig, D., and Charlson, M. (1997b). Clinically defined vascular depression. Am. J. Psychiatry 154, 562–565. doi: 10.1176/ajp.154.4.562
Alexopoulos, G. S., Murphy, C. F., Gunning-Dixon, F. M., Latoussakis, V., Kanellopoulos, D., Klimstra, S., et al. (2008). Microstructural white matter abnormalities and remission of geriatric depression. Am. J. Psychiatry 165, 238–244. doi: 10.1176/appi.ajp.2007.07050744
Ayerbe, L., Ayis, S., Wolfe, C. D. A., and Rudd, A. G. (2013). Natural history, predictors and outcomes of depression after stroke: systematic review and meta-analysis. Br. J. Psychiatry 202, 14–21. doi: 10.1192/bjp.bp.111.107664
Babkair, L. A. (2017). Risk factors for poststroke depression: an integrative review. J. Neurosci. Nurs. 49, 73–84. doi: 10.1097/JNN.0000000000000271
Bajwa, S., Bermpohl, F., Rigonatti, S. P., Pascual-Leone, A., Boggio, P. S., and Fregni, F. (2008). Impaired interhemispheric interactions in patients with major depression. J. Nerv. Mental Dis. 196, 671–677. doi: 10.1097/NMD.0b013e318183f86f
Bassett, D. S., and Bullmore, E. T. (2009). Human brain networks in health and disease. Curr. Opin. Neurol. 22, 340–347. doi: 10.1097/WCO.0b013e32832d93dd
Bays, C. L. (2001). Quality of life of stroke survivors: a research synthesis. J. Neurosci. Nurs. 33, 310–316. doi: 10.1097/01376517-200112000-00005
Bernhardt, B. C., Hong, S., Bernasconi, A., and Bernasconi, N. (2013). Imaging structural and functional brain networks in temporal lobe epilepsy. Front. Hum. Neurosci. 7:624. doi: 10.3389/fnhum.2013.00624
Berthier, M. L. (2005). Poststroke aphasia: epidemiology, pathophysiology and treatment. Drugs Aging 22, 163–182. doi: 10.2165/00002512-200522020-00006
Brakowski, J., Spinelli, S., Dörig, N., Bosch, O. G., Manoliu, A., Holtforth, M. G., et al. (2017). Resting state brain network function in major depression-depression symptomatology, antidepressant treatment effects, future research. J. Psychiatr. Res. 92, 147–159. doi: 10.1016/j.jpsychires.2017.04.007
Cai, X., Kallarackal, A. J., Kvarta, M. D., Goluskin, S., Gaylor, K., Bailey, A. M., et al. (2013). Local potentiation of excitatory synapses by serotonin and its alteration in rodent models of depression. Nat. Neurosci. 16, 464–472. doi: 10.1038/nn.3355
Chen, J., Venkat, P., Zacharek, A., and Chopp, M. (2014). Neurorestorative therapy for stroke. Front. Hum. Neurosci. 8:382. doi: 10.3389/fnhum.2014.00382
Cojocaru, G. R., Popawagner, A., Stanciulescu, E. C., Babadan, L., and Buga, A. M. (2013). Post-stroke depression and the aging brain. J. Mol. Psychiatry 1:14. doi: 10.1186/2049-9256-1-14
Coronel, C., Garn, H., Waser, M., Deistler, M., Benke, T., Dal-Bianco, P., et al. (2017). Quantitative EEG markers of entropy and auto mutual information in relation to MMSE scores of probable Alzheimer's disease patients. Entropy 19:130. doi: 10.3390/e19030130
Dong, L., Li, F., Liu, Q., Wen, X., Lai, Y., Xu, P., et al. (2017). MATLAB toolboxes for reference electrode standardization technique (REST) of scalp EEG. Front. Neurosci. 11:601. doi: 10.3389/fnins.2017.00601
Draganski, B., Kherif, F., Kloppel, S., Cook, P. A., Alexander, D. C., Parker, G. J., et al. (2008). Evidence for segregated and integrative connectivity patterns in the human Basal ganglia. J. Neurosci. 28, 7143–7152. doi: 10.1523/JNEUROSCI.1486-08.2008
Drevets, W. C., Price, J. L., and Furey, M. L. (2008). Brain structural and functional abnormalities in mood disorders: implications for neurocircuitry models of depression. Brain Struct. Funct. 213, 93–118. doi: 10.1007/s00429-008-0189-x
Dumas, J. A., and Newhouse, P. A. (2015). Impaired working memory in geriatric depression: an fMRI study. Am. J. Geriatr. Psychiatry 23, 433–436. doi: 10.1016/j.jagp.2014.09.011
Engelter, S. T., Gostynski, M., Papa, S., Frei, M., Born, C., Ajdacic-Gross, V., et al. (2006). Epidemiology of aphasia attributable to first ischemic stroke: incidence, severity, fluency, etiology, and thrombolysis. Stroke 37, 1379–1384. doi: 10.1161/01.STR.0000221815.64093.8c
Fanciullacci, C., Bertolucci, F., Lamola, G., Panarese, A., Artoni, F., Micera, S., et al. (2017). Delta power is higher and more symmetrical in ischemic stroke patients with cortical involvement. Front. Hum. Neurosci. 11:385. doi: 10.3389/fnhum.2017.00385
Fingelkurts, A. A., and Fingelkurts, A. A. (2006). Timing in cognition and EEG brain dynamics: discreteness versus continuity. Cogn. Process. 7, 135–162. doi: 10.1007/s10339-006-0035-0
Fingelkurts, A. A., and Fingelkurts, A. A. (2015). Altered structure of dynamic electroencephalogram oscillatory pattern in major depression. Biol. Psychiatry 77, 1050–1060. doi: 10.1016/j.biopsych.2014.12.011
Fonov, V., Evans, A. C., Botteron, K., Almli, C. R., McKinstry, R. C., and Collins, D. L. (2011). Unbiased average age-appropriate atlases for pediatric studies. Neuroimage 54, 313–327. doi: 10.1016/j.neuroimage.2010.07.033
Fonov, V., Evans, A., McKinstry, R., Almli, C., and Collins, D. (2009). Unbiased nonlinear average age-appropriate brain templates from birth to adulthood. Neuroimage 47:S102. doi: 10.1016/S1053-8119(09)70884-5
Ghose, S. S., Williams, L. S., and Swindle, R. W. (2005). Depression and other mental health diagnoses after stroke increase inpatient and outpatient medical utilization three years poststroke. Med. Care 43, 1259–1264. doi: 10.1097/01.mlr.0000185711.50480.13
Greicius, M. D., Flores, B. H., Menon, V., Glover, G. H., Solvason, H. B., Kenna, H., et al. (2007). Resting-state functional connectivity in major depression: abnormally increased contributions from subgenual cingulate cortex and thalamus. Biol. Psychiatry 62, 429–437. doi: 10.1016/j.biopsych.2006.09.020
Grin-Yatsenko, V. A., Baas, I., Ponomarev, V. A., and Kropotov, J. D. (2009). EEG power spectra at early stages of depressive disorders. J. Clin. Neurophysiol. 26, 401–406. doi: 10.1097/WNP.0b013e3181c298fe
Grin-Yatsenko, V. A., Baas, I., Ponomarev, V. A., and Kropotov, J. D. (2010). Independent component approach to the analysis of EEG recordings at early stages of depressive disorders. Clin. Neurophysiol. 121, 281–289. doi: 10.1016/j.clinph.2009.11.015
Gudayol-Ferre, E., Pero-Cebollero, M., Gonzalez-Garrido, A. A., and Guardia-Olmos, J. (2015). Changes in brain connectivity related to the treatment of depression measured through fMRI: a systematic review. Front. Hum. Neurosci. 9:582. doi: 10.3389/fnhum.2015.00582
Haufe, S., Nikulin, V. V., Müller, K. R., and Nolte, G. (2013). A critical assessment of connectivity measures for EEG data: a simulation study. Neuroimage 64, 120–133. doi: 10.1016/j.neuroimage.2012.09.036
Ince, R. A. A., Giordano, B. L., Kayser, C., Rousselet, G. A., Gross, J., and Schyns, P. G. (2017). A statistical framework for neuroimaging data analysis based on mutual information estimated via a gaussian copula. Hum. Brain Mapp. 38, 1541–1573. doi: 10.1002/hbm.23471
Jeong, J., Gore, J. C., and Peterson, B. S. (2001). Mutual information analysis of the eeg in patients with alzheimer's disease. Clin. Neurophysiol. 112, 827–835. doi: 10.1016/S1388-2457(01)00513-2
Lao, Y., Wang, Y., Shi, J., Ceschin, R., Nelson, M. D., Panigrahy, A., et al. (2016). Thalamic alterations in preterm neonates and their relation to ventral striatum disturbances revealed by a combined shape and pose analysis. Brain Struct. Funct. 221, 487–506. doi: 10.1007/s00429-014-0921-7
Marzetti, L., Nolte, G., Perrucci, M. G., Romani, G. L., and Del Gratta, C. (2007). The use of standardized infinity reference in EEG coherency studies. Neuroimage 36, 48–63. doi: 10.1016/j.neuroimage.2007.02.034
Na, S. H., Jin, S.-H., Kim, S. Y., and Ham, B.-J. (2002). EEG in schizophrenic patients: mutual information analysis. Clin. Neurophysiol. 113, 1954–1960. doi: 10.1016/S1388-2457(02)00197-9
Nolte, G., Bai, O., Wheaton, L., Mari, Z., Vorbach, S., and Hallett, M. (2004). Identifying true brain interaction from eeg data using the imaginary part of coherency. Clin. Neurophysiol. 115, 2292–2307. doi: 10.1016/j.clinph.2004.04.029
Nys, G. M., van Zandvoort, M. J., de Kort, P. L., Jansen, B. P., de Haan, E. H., and Kappelle, L. J. (2007). Cognitive disorders in acute stroke: prevalence and clinical determinants. Cerebrovasc. Dis. 23, 408–416. doi: 10.1159/000101464
Qin, Y., Xu, P., and Yao, D. (2010). A comparative study of different references for EEG default mode network: the use of the infinity reference. Clin. Neurophysiol. 121, 1981–1991. doi: 10.1016/j.clinph.2010.03.056
Rathee, D., Cecotti, H., and Prasad, G. (2016). “Estimation of effective fronto-parietal connectivity during motor imagery using partial granger causality analysis,” in Proceedings of International Joint Conference on Neural Networks 2016–Octob (Vancouver, BC).
Rathee, D., Cecotti, H., and Prasad, G. (2017). “Propofol-induced sedation diminishes the strength of frontal-parietal-occipital EEG network,” in Proceedings of Annual International Conference of the IEEE Engineering Medical and Biology Society EMBS (Jeju Island), 4463–4466.
Rubinov, M., Knock, S. A., Stam, C. J., Micheloyannis, S., Harris, A. W. F., Williams, L. M., et al. (2009). Small-world properties of nonlinear brain activity in schizophrenia. Hum. Brain Mapp. 30, 403–416. doi: 10.1002/hbm.20517
Salvador, R., Suckling, J., Schwarzbauer, C., and Bullmore, E. (2005). Undirected graphs of frequency-dependent functional connectivity in whole brain networks. Phil. Trans. R. Soc. B, 937–946. doi: 10.1098/rstb.2005.1645
Schreiber, T. (2000). Measuring information transfer. Phys. Rev. Lett. 85, 461–464. doi: 10.1103/PhysRevLett.85.461
Stam, C. J., Nolte, G., and Daffertshofer, A. (2010). Phase lag index: assessment of functional connectivity from multi channel eeg and meg with diminished bias from common sources. Hum. Brain Mapp. 28, 1178–1193. doi: 10.1002/hbm.20346
Tadayonnejad, R., and Ajilore, O. (2014). Brain network dysfunction in late-life depression: a literature review. J. Geriatr. Psychiatry Neurol. 27, 5–12. doi: 10.1177/0891988713516539
Tadayonnejad, R., Yang, S., Kumar, A., and Ajilore, O. (2014). Multimodal brain connectivity analysis in unmedicated late-life depression. PLoS ONE 9:e96033. doi: 10.1371/journal.pone.0096033
Thomas, E. J. (2009). Brain imaging correlates of cognitive impairment in depression. Front. Hum. Neurosci. 3:30. doi: 10.3389/neuro.09.030.2009
Wang, Z. J., Lee, P. W.-H., and McKeown, M. J. (2009). A novel segmentation, mutual information network framework for EEG analysis of motor tasks. Biomed. Eng. Online 8:9. doi: 10.1186/1475-925X-8-9
Wei, L., Li, Y., Tang, Y., and Wang, J. (2010). [Brain functional connectivity in depression during emotion processing]. J. Biomed. Eng. 27:961.
Williams, D., Tijssen, M., Van Bruggen, G., Bosch, A., Insola, A., Di Lazzaro, V., et al. (2002). Dopamine-dependent changes in the functional connectivity between basal ganglia and cerebral cortex in humans. Brain, 125, 1558–1569. doi: 10.1093/brain/awf156
Wu, Y. T. (2014). Decreased inter-hemispheric connectivity in anterior sub-network of default mode network and cerebellum: significant findings in major depressive disorder. Int. J. Neuropsychopharmacol. 17, 1935–1942. doi: 10.1017/S1461145714000947
Yamada, M., Kimura, M., Mori, T., and Endo, S. (1995). [EEG power and coherence in presenile and senile depression. Characteristic findings related to differences between anxiety type and retardation type. Nihon Ika Daigaku Zasshi 62:176. doi: 10.1272/jnms1923.62.176
Yao, D. (2001). A method to standardize a reference of scalp EEG recordings to a point at infinity. Physiol. Meas. 22, 693–711. doi: 10.1088/0967-3334/22/4/305
Yao, D. (2017). Is the surface potential integral of a dipole in a volume conductor always zero? A cloud over the average reference of EEG and ERP. Brain Topogr. 30, 161–171. doi: 10.1007/s10548-016-0543-x
Keywords: post-stroke depression (PSD), electroencephalography (EEG), mutual information (MI), graph theory, brain network
Citation: Sun C, Yang F, Wang C, Wang Z, Zhang Y, Ming D and Du J (2018) Mutual Information-Based Brain Network Analysis in Post-stroke Patients With Different Levels of Depression. Front. Hum. Neurosci. 12:285. doi: 10.3389/fnhum.2018.00285
Received: 18 January 2018; Accepted: 25 June 2018;
Published: 17 July 2018.
Edited by:
Raffaella Franciotti, Università degli Studi G. d'Annunzio Chieti e Pescara, ItalyReviewed by:
Dezhong Yao, University of Electronic Science and Technology of China, ChinaLiYu Huang, Xidian University, China
Copyright © 2018 Sun, Yang, Wang, Wang, Zhang, Ming and Du. 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) and the copyright owner(s) 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: Jingang Du, cm15eTEyNmR1QHNpbmEuY29t