To the best of our knowledge, this is the first study to investigate the brain functional connectivity in patients with PSE using source-level EEG analysis and graph theory. The main finding of this study was that there were significant alterations of functional connectivity in stroke patients with PSE compared to those without PSE. In addition, the differences of functional connectivity measures between patients with and without PSE were dependent on EEG frequency bands. There were no differences in the functional connectivity between the groups, especially in the delta, theta, and alpha bands. In the beta band, the radius and diameter were increased, and in the low gamma band, the radius was increased in patients with PSE compared to those without PSE. In the high gamma band, the radius, diameter, eccentricity, and characteristic path length were increased, whereas the average strength, global efficiency, local efficiency, mean clustering coefficient, and transitivity were decreased in the patients with PSE than in those without PSE.
Several cases have demonstrated functional connectivity alterations in patients with stroke. (Caliandro et al. 2017; Vecchio et al. 2019) Graph characteristics based on the EEG in 30 patients with acute stroke were investigated in a previous study. It found that brain network rearrangements, mainly decreased small-worldness, were detected in patients with stroke compared with healthy subjects. (Caliandro et al. 2017) Another study examined whether acute cerebellar stroke could cause alterations in brain network architecture as defined by cortical EEG rhythm. The authors demonstrated that the network remodeling characteristics were independent of the size of the ischemic lesion, and that cerebellar acute stroke induced significant changes in the balance of local segregation and global integration. (Vecchio et al. 2019) However, there was no research investigating the alterations of function connectivity in patients with PSE.
We demonstrated the alterations of functional connectivity in patients with PSE compared to stroke patients without PSE. The radius, diameter, eccentricity, and characteristic path length were increased, whereas the average strength, global efficiency, local efficiency, mean clustering coefficient, and transitivity were decreased in the patients with PSE than those without PSE, especially in the beta and gamma bands. The radius of a graph represents the minimum eccentricity of any graph vertex, while the diameter of a graph represents the maximum eccentricity of any graph vertex. (Collantoni et al. 2022; Falsaperla et al. 2021; Farahani et al. 2019; Mijalkov et al. 2017; Thomas et al. 2016; Zenil et al. 2018) The eccentricity is defined as the greatest distance between two vertices. The maximum distance between a vertex and every other vertex is regarded as the vertex's eccentricity. The average distance between pairs of vertices is known as the characteristic path length of the graph. The maximum number of edges that connect any two adjacent vertices in a multigraph. The global efficiency is a measure of the effectiveness of distant information transfer in a network and is defined as the inverse of the typical path length between all nodes. The local efficiency is a measure of the average efficiency of information transfer within local subgraphs or neighborhoods. It is defined as the inverse of the shortest average path length between all neighbors of a given node. The clustering coefficient quantifies the degree to which graph nodes tend to cluster together. The transitivity is the probability that adjacent network nodes are interconnected, revealing the existence of densely interconnected communities. (Collantoni et al. 2022; Falsaperla et al. 2021; Farahani et al. 2019; Mijalkov et al. 2017; Thomas et al. 2016; Zenil et al. 2018) Therefore, all these results suggest that functional segregation and integration are decreased in patients with PSE compared to those without PSE. Our finding in the present study suggests that post-stroke epileptogenesis is associated with alterations of the functional brain connectivity. This is consistent with a previous study with a rat model demonstrating that post-stroke epileptogenesis was related with a rearrangement of the intrinsic excitability across neurons. (Vera and Lippmann 2021)
Interestingly, the difference in functional connectivity between the two groups varied according to the EEG frequency bands. The differences of functional connectivity between the groups were revealed, especially in the beta and gamma bands. In the delta, theta, and alpha bands, all of the functional connectivity measures were not different between the groups. Brain oscillation has distinct function at different frequencies in various neural processing such as perceptual, sensorimotor and cognitive function. (Wang et al. 2017) Low frequency oscillations (delta, theta and alpha frequency bands) are associated with the integration across wide spatially large-scale networks, while high frequency oscillations (beta and gamma frequency bands) are linked to a more limited topographic area. (Wang et al. 2017) Beta and gamma frequency oscillations are known to occur in the inhibitory gamma-aminobutyric acid (GABA) interneuron network, and they are verified to relate to cognitive and executive functions. (Missonnier et al. 2010; Rossiter et al. 2014; Wang et al. 2017) The GABAergic interneurons receive excitatory inputs through N-methyl-D-aspartate (NMDA)- receptors and dysfunction of NMDA neurotransmission may result in abnormalities of gamma oscillation. (González-Burgos and Lewis 2012) The alterations of brain functional connectivity in patients with PSE in beta and gamma frequency bands might reflect the imbalance of GABAergic interneuron network induced by NMDA-receptor over-activation after stroke. (Nicolo et al. 2019)
Since the process of information in brain is encoded at a very short time scale (milliseconds to seconds), EEG is known as a good tool for evaluating the resting function with good time resolution. (Van Diessen et al. 2015) In addition, EEG is more advantageous than functional MRI in demonstrating mechanism of signal process among neuronal assemblies which do not modify energy consumption such as synchronization, coherence or phase locking-unlocking without any change of firing frequency, which reflect resting function of brain. (Rossini et al. 2019) However, there are some problems when analyzing sensor-level in performing connectivity analysis with EEG: field spread and volume conduction. (Hatlestad-Hall et al. 2021; Van Diessen et al. 2015) In addition, it is difficult to ensure resolution in terms of neuroanatomy in sensor level EEG analysis, as access to the source-level is required. (Lai et al. 2018b) Several solutions including Minimum Norm(MN), Dynamical Statistical Parametric Mapping(dSPM), Standard Maximum Entropy on the Mean(cMEM), Local AutoRegressive Average(LAURA) and LORETA, are introduced to solve the inverse problem that occurs during the reconstruction of EEG from the sensor-level to the source-level analysis. (Michel and Brunet 2019) These source reconstruction solutions assume the concept of a lead field that each signal of the scalp electrodes is produced by a source of unitary strength. (Van Diessen et al. 2015)
In the present study, we chose the combination of sLORETA for inverse solution and phase locking value for connectivity measure to analyze source-level EEG connectivity. In a recent study comparing source-level functional connectivity analysis and classical source localization in patients with epilepsy, source-level functional connectivity analysis showed high similarity to real interictal epileptic networks. (Hassan et al. 2017) In addition, it is important to select the optimal combination of inverse solution and connectivity measure. (Rizkallah et al. 2019) As a feature of the epileptic network, functionally connected brain regions have hyper-synchronization phenomena, and phase locking value is one of the viable methods because it has the advantage of nonlinearity network analysis. (Hatlestad-Hall et al. 2021; Rizkallah et al. 2019; Van Diessen et al. 2015)
This study has some limitations. First, the sample size was relatively small, and the study was conducted at a single center. Second, this study was a cross-sectional study. It is necessary to confirm whether there is a longitudinal difference in functional connectivity in patients with or without PSE, respectively. Third, time interval between stroke onset and EEG taken was different between the two groups. EEG were performed after more than one late-onset seizure was identified in patients with PSE, whereas in most of the patients without PSE, EEG were performed within a week after the stroke diagnosis. However, the effects of anti-seizure medications on EEG could be excluded because EEG were performed before anti-seizure medications were administered. Lastly, we only used 23 electrodes were placed using the international 10–20 system. In general, high density EEG is advantageous considering the forward problem from source to sensor. (Van Diessen et al. 2015)