MATERIALS & METHODS
- Patients
We retrospectively reviewed the medical records of 32 de novo PD patients who underwent EEG between November 2015 and October 2021. PD was diagnosed according to the UK Brain Bank criteria. We collected the following information from patients: Hoehn and Yahr stage (H&Y stage), Unified Parkinson’s Disease Rating Scale (UPDRS) Part I–III, the Korean Version of the Mini-Mental State Examination (K-MMSE), the Korean Version of the Montreal Cognitive Assessment (K-MoCA), education level, disease duration, Beck Depression Inventory (BDI), Fatigue Severity Scale (FSS), and the Korean version of the Nonmotor Symptoms Scale (NMSS) [
6]. Depression was noted if the BDI score was 20 or higher [
1], and fatigue was based on an FSS score of four or more [
1]. The exclusion criteria were cognitive dysfunction (K-MMSE score less than the 2.5th percentile for age and educational-appropriate norm), and a history of stroke, epilepsy, pertinent head injury, or significant debilitating medical conditions [
7]. The Institutional Review Board of Hallym University Dongtan Sacred Heart Hospital approved this study (IRB No. 2021-07-016). The need for informed consent was waived because of the retrospective nature of the study. All methods were performed in accordance with the Declaration of Helsinki and followed relevant guidelines and regulations.
- EEG recording
Participants were asked to be relaxed and awake during the recording. Eyes-closed, resting-state EEG data were recorded for at least five minutes (Comet plus; Grass, West Warwick, RI, USA). The International 10–20 system was used for electrode placement, and 19 channels with referential montage were used: FP1, FP2, F7, F3, Fz, F4, F8, T3, C3, Cz, C4, T4, T5, P3, Pz, P4, T6, O1, and O2. One additional electrode applied to the left mastoid was used as a ground, and another electrode applied to the right mastoid was used as a reference. The EEG data were sampled at 200 Hz, and the bandpass filter was set between 1.0 Hz and 70 Hz. The impedance of every electrode was optimized to below 10 kΩ. Resting-state EEG was performed for 10 minutes. We selected three minutes of eyes-closed and artifact-free data based on visual inspection for further analysis. The mean time interval between clinical evaluation and EEG recording was 1.0 ± 2.4 days.
- Data preprocessing
Quantitative EEG (QEEG) processing and group analyses were conducted using iSyncBrain (iMediSync Inc., Seoul, Korea;
https://isyncbrain.com), a cloud-based, artificial intelligence EEG analysis platform [
8]. EEG preprocessing was performed to denoise all the data and minimize the effects of artifacts. During the first stage of EEG preprocessing, the signals were sampled at 250 Hz and filtered with a bandpass filter in the 1–45.5 Hz range. The EEG signals were then passed through a notch filter in preparation for downstream processing, including re-referencing (using common average reference), wrong epoch rejection (using artifact subspace reconstruction), and adaptive mixture independent component analysis (AMICA). AMICA helps provide a generic, unsupervised approach to identifying and modeling changes in EEG dynamics [
9]. Finally, artifacts identified via electromyography, cardiac signals, body movement, and electrooculography were removed to yield the clean QEEG normative data. To enhance the performance of iSyncBrain, we added an extra denoising process that included using a bandpass filter from 1 Hz to 45.5 Hz, re-referencing using a common average reference, and determining artifacts by AMICA [
10].
- Data processing and analyses
All EEG preprocessing processes, sensor-level data, source-level data calculations, and extractions were performed using a cloud-based EEG analysis platform (iSyncBrain). The power spectral density of the EEG rhythms was computed using Fast Fourier transform (FFT) analysis with 0.25 Hz of frequency resolution using iSyncBrain. Then, the signal was decomposed into the following frequency bands: delta (1–4 Hz), theta (4–8 Hz), alpha 1 (8–10 Hz), alpha 2 (10–12 Hz), beta 1 (12–15 Hz), beta 2 (15–20 Hz), beta 3 (20–30 Hz), and gamma (30–45 Hz). The alpha bands and beta bands were divided for more granular frequency analysis [
11-
13].
In the power spectrum analysis, the absolute and relative power of the seven frequency bands were measured using a FFT. Absolute band power is a spectral band power based on FFT, here provided by iSyncBrain. The relative power is the ratio of the power in a given frequency band to the power in the total frequency range. The topographical grouping of the electrodes was defined as frontal (Fp1, Fp2, F7, F3, Fz, F4, and F8), central (C3, Cz, and C4), temporal (T3, T4, T5, and T6), parietal (P3, Pz, and P4), or occipital (O1 and O2).
Source reconstructions were performed with sLORETA [
14] using the Colin 27 Head model [
15] with 68 ROI segmentations based on the Desikan-Killiany atlas [
16].
In the brain network analysis, we used a popular graph theory approach that provides quantitative undirected binary metrics for measuring brain connectivity networks [
17,
18]. A network is a mathematical representation of a real-world complex system and is defined by a collection of nodes (vertices) and links (edges) between pairs of nodes [
19]. Each ROI from the sLORETA analysis was considered a node. For the network analysis, we divided the brain into six regions and performed a network analysis for each region. We adopted the most commonly used network measures classified into global and local measurements. Global measurements included the global efficiency (GE), characteristic path length (CPL), clustering coefficient (CCO), and small-worldness (SW). The local measurements included the local efficiency (LE), degree centrality (DC), closeness centrality (CCE), and betweenness centrality (BC) [
19,
20]. Detailed explanations of these measurements are provided elsewhere [
19-
22]. Briefly, GE is the efficiency of remote information transmission in a network [
23]. It indicates the average effect of the relevant brain network and its overall information transmission capacity [
22]. The CPL measures the connectivity level of nodes in the network; lower values indicate more efficient information transfer between two nodes and better connectivity of the network [
22]. The CCO is the ratio of the number of current edges between adjacent nodes to all possible connected edges; a higher value indicates a more highly segregated network. SW is the ratio of the normalized CCO to the normalized CPL and shows networks that are significantly more clustered than random networks, i.e., the density and efficiency of information exchange in the networks. The LE is the average efficiency of information transfer between a node and its nearest neighbors [
24]. A higher LE means that the information transmission of adjacent nodes is more efficient and that the network is less differentiated [
22,
25]. DC is the number of adjacent nodes; a high DC means a more central node. CCE is the closeness of a node in the network to all other nodes (i.e., a measure of the relative importance of a node) [
26]. Higher values of a specific node mean that it is close to all other nodes. BC is the tendency of a single node to be more central than all other nodes and the quantification that acts as a bridge between two other connected nodes along its shortest path length [
27]. A higher BC value indicates a node closer to the center of the network. The brain network nodes represent the brain region, and the edges indicate connections between neurons or brain regions [
20].
- Statistical analysis
The statistical analysis was designed to examine the clinical characteristics and EEG data measured from 19 channels in two groups of PD patients: those without fatigue and those with fatigue. We used independent t tests (or Mann–Whitney U tests) and Fisher’s exact tests to compare the two groups, as appropriate. This analysis included a normality test using the Kolmogorov–Smirnov test for each group (
Supplementary Tables 1 and
2 in the online-only Data Supplement). We also used one-way analysis of covariance to compare the EEG data between the two groups after adjusting for the BDI score. We did not find significant differences between the two groups in terms of the BDI score, but it in a prior study, it was argued that the BDI score was strongly correlated with fatigue in PD patients [
3], and it is also worth considering that the
p value approached significance (
p = 0.059) (
Table 1). The EEG data, obtained by averaging the power spectrum within each five lobe regions, and by averaging network activity across six lobe regions, were compared between PD patients without fatigue and those with fatigue. Significance was determined at a
p value <0.05. Significance levels were Bonferroni-corrected for multiple comparisons:
p value <0.05/5 (0.01) for the power spectrum analysis and <0.05/6 (0.008333) for the network analysis. Demographic and clinical data were analyzed using IBM SPSS 25 statistics (IBM Corp., Armonk, NY, USA).
DISCUSSION
Our study demonstrated several differing features between PD patients with and without fatigue via QEEG analysis, especially when using graph theory. The main findings were as follows: 1) We did not find significantly different features in the power spectrum analysis between PD patients without fatigue and PD patients with fatigue. 2) In global measurements, the efficient network according to the frequency bands seemed to vary by brain region. In PD patients without fatigue, the theta band showed a better overall information transmission capacity in the frontal region (GE values), and the delta band showed a more efficient network, whereas the alpha1 band showed less effective connectivity in the temporal area (GE and CPL values). 3) Regarding local measurements, the results were similar to those of the global network analysis. In PD patients without fatigue, more centrality was observed in the frontal areas, and less centrality was observed in the left temporal areas. For the theta band, DC was greater in frontal regions (i.e., more central areas). CCE was greater for the theta band (i.e., relatively more important) in the parieto-occipital area and for the delta band in the temporal area in PD patients without fatigue. The alpha1 band was less prominent in the temporal areas (DC and CCE values).
There have been no prior power spectrum analysis studies in PD patients with fatigue. Although our results revealed no differences between PD patients without fatigue and those with fatigue, a meta-analysis on the influence of fatigue on brain activity revealed increases in overall EEG activity (particularly in the theta and alpha bands) in the fatigued state [
28]. The delta and alpha1 absolute powers in PD patients without fatigue seemed to be lower in all brain regions (
Figure 1). It is not known why the EEG spectrum differs between PD patients with and without fatigue. Because the decreased power is known to be associated with increased cortical activity or an alert state [
28,
29], the lower EEG power in PD patients without fatigue may be related to the more active cerebral function in the awake state. The higher EEG power (i.e., synchronization) in PD patients with fatigue indicates an idle state (i.e., a less active state).
Our network analysis using graph theory showed that the theta band in the frontal area was related to fatigue in PD patients because the network of the frequency band in the frontal area was less efficient in PD patients with fatigue than in PD patients without fatigue. The theta band in the frontal area is related to the recruitment of cognitive control during monitoring, the management of memory function, and error processing [
30]. These cognitive functions are reportedly related to fatigue in healthy people and people with chronic fatigue syndrome [
31,
32] or PD [
33]. Previous brain network studies support our results. Theta band functional coupling was weakened in PD patients with fatigue [
20,
34]. However, because there were also contradictory results in other studies, possibly due to different study designs [
20], more research may be needed.
It seemed that the delta band was also associated with fatigue because of the lower connectivity in the temporal area and the lower centrality in the frontal and temporal areas in PD patients with fatigue. We think that the delta band findings should be interpreted with caution because of possible EEG artifact contamination and previous contradictory reports. Some studies have shown that increased EEG power is correlated with fatigue, but others have not [
35-
37].
It seemed that the network activity for the alpha band in the temporal regions was a compensatory mechanism for fatigue in our study. A previous brain source analysis study may support our results. They showed greater source-current activity in the frontal-temporal-parietal regions in chronic fatigue syndrome patients [
38]. Since alpha bands are involved in attention and memory processes [
39], there may be an association between the efficiency of these brain functions and the fatigue process. Generally, alpha- and beta-band functional coupling are weakened in individuals with fatigue [
39,
40], but our results were not significant after multiple comparison correction.
We do not know the exact reason why there was no difference between patients without fatigue and those with fatigue in the power spectrum analysis but was a difference in the network analysis. A previous study also revealed significant differences in network measures between PD patients and healthy controls but this difference was not detected via power spectral analysis. They assumed that changes in the power spectrum were probably not evident in the early stages of PD because they studied de novo, early-stage PD patients, similar to our study [
41]. Generally, power spectrum analysis shows cerebral activity in localized brain areas, and network analysis shows long-range communication among several brain areas [
42]. Traditionally, localized brain areas have been investigated for neurologic symptoms, but network analysis of connected regions has gained much interest because it better explains neurological symptoms [
43,
44]. Brain connectivity using graph theory describes the characteristics of brain architecture, such as efficient communication between two brain areas. Another possibility is that fatigue in PD patients is less related to local brain oscillations and is a problem in the brain network. Brain oscillations play an essential role in both local activity and long-range communication. Many oscillatory bands in cerebral networks range from 0.05 Hz to 500 Hz, and their physiological roles are poorly understood. The function of local oscillations seems to differ from that of long-range communication, even in the same frequency bands [
45].
Our study has some limitations. First, the number of patients studied was relatively small. Nevertheless, we think that our EEG results, especially the network analysis (our examined parameters), showed consistent findings. Second, we did not perform detailed sleep evaluations to exclude patients with severe sleep problems. Sleep problems can be a confounding factor because they can affect fatigue. However, there was no significant difference in sleep scores on the NMSS between PD patients without and with fatigue (numbers 3, 5, and 6;
p > 0.100 in the
Supplementary Table 3 in the online-only Data Supplement). Third, it can be argued that the sample size was biased toward women, especially in the PD population. However, there was no significant difference in the sex ratio between the PD patients with and without fatigue. In addition, the female-to-male ratio of the incidence and prevalence of PD in Korea is higher (1.4:1 and 1.6:1, respectively) [
33]. This finding is different from that of Western population studies showing that PD is more predominant in men but is similar to that of many Asian population studies [
33].
In conclusion, our study showed that EEG variables can be biomarkers of fatigue in PD, and that PD patients with fatigue have less efficient brain networks in the frontal area than PD patients without fatigue. These findings help us understand the brain mechanisms associated with fatigue in PD. In addition, brain network analysis using graph theory is more valuable than power spectrum analysis in revealing the brain mechanisms related to fatigue.