ORIGINAL RESEARCH article

Front. Aging Neurosci., 07 October 2021

Sec. Neurocognitive Aging and Behavior

Volume 13 - 2021 | https://doi.org/10.3389/fnagi.2021.680200

EEG Alpha and Beta Band Functional Connectivity and Network Structure Mark Hub Overload in Mild Cognitive Impairment During Memory Maintenance

  • 1. Department of Psychiatry and Psychotherapy, Semmelweis University, Budapest, Hungary

  • 2. Department of Neurology, National Institute of Clinical Neurosciences, Budapest, Hungary

  • 3. Department of Clinical Neurophysiology, Amsterdam Neuroscience, Vrije Universiteit Amsterdam, Amsterdam University Medical Centers, Amsterdam, Netherlands

  • 4. Department of Neurology, Alzheimer Center Amsterdam, Amsterdam Neuroscience, Vrije Universiteit Amsterdam, Amsterdam University Medical Centers, Amsterdam, Netherlands

Abstract

Background: While decreased alpha and beta-band functional connectivity (FC) and changes in network topology have been reported in Alzheimer’s disease, it is not yet entirely known whether these differences can mark cognitive decline in the early stages of the disease. Our study aimed to analyze electroencephalography (EEG) FC and network differences in the alpha and beta frequency band during visuospatial memory maintenance between Mild Cognitive Impairment (MCI) patients and healthy elderly with subjective memory complaints.

Methods: Functional connectivity and network structure of 17 MCI patients and 20 control participants were studied with 128-channel EEG during a visuospatial memory task with varying memory load. FC between EEG channels was measured by amplitude envelope correlation with leakage correction (AEC-c), while network analysis was performed by applying the Minimum Spanning Tree (MST) approach, which reconstructs the critical backbone of the original network.

Results: Memory load (increasing number of to-be-learned items) enhanced the mean AEC-c in the control group in both frequency bands. In contrast to that, after an initial increase, the MCI group showed significantly (p < 0.05) diminished FC in the alpha band in the highest memory load condition, while in the beta band this modulation was absent. Moreover, mean alpha and beta AEC-c correlated significantly with the size of medial temporal lobe structures in the entire sample. The network analysis revealed increased maximum degree, betweenness centrality, and degree divergence, and decreased diameter and eccentricity in the MCI group compared to the control group in both frequency bands independently of the memory load. This suggests a rerouted network in the MCI group with a more centralized topology and a more unequal traffic load distribution.

Conclusion: Alpha- and beta-band FC measured by AEC-c correlates with cognitive load-related modulation, with subtle medial temporal lobe atrophy, and with the disruption of hippocampal fiber integrity in the earliest stages of cognitive decline. The more integrated network topology of the MCI group is in line with the “hub overload and failure” framework and might be part of a compensatory mechanism or a consequence of neural disinhibition.

Introduction

Deteriorated working memory maintenance and the impairment of visuospatial memory are early symptoms of Mild Cognitive Impairment (MCI) and Alzheimer’s disease (AD) (; Parra et al., 2010; ; Moodley et al., 2015) and can serve as a sensitive marker of early cognitive decline (Tierney et al., 1996; Sano et al., 2011). Visuospatial memory tests, such as the Paired Associates Learning (PAL) test are considered especially effective in the early diagnosis of MCI (Sirály et al., 2013) and in the prediction of a higher risk of developing dementia in later life ().

Cognitive functions arise from the interactions between functionally connected regions of the brain (Rubinov and Sporns, 2010; Park and Friston, 2013; Stam, 2014). However, besides sufficient connections, proper cognitive functioning relies on an optimal organization of brain network () and the coordinated interaction of local information processing (“segregation”) and the long-range integration of this information (Sporns, 2013; Stam, 2014). A growing body of evidence suggests that healthy brain networks are cost-efficient small-world networks combining strong local connectivity with efficient long-distance connections (). Furthermore, it has been shown that brain network efficiency is related to cognitive performance (van den Heuvel et al., 2009) and network measures derived from electrophysiological data can discriminate cortical network features in healthy brain and neurodegenerative brain aging (Miraglia et al., 2017; Vecchio et al., 2017).

The pathological process of AD initially affects synaptic transmission with an overall disconnection (), which could be assessed using a network approach as the structural components of the brain form a complex network at different spatial scale (from neurons to anatomical regions) from which functional dynamics arise (Vecchio et al., 2017). The abnormal functional brain network in AD has been characterized by a loss of efficiency, disturbed community structure, and selective hub vulnerability in both structural and functional network studies (Tijms et al., 2013; Stam, 2014; Miraglia et al., 2017). Furthermore, the extent of network changes correlates with the extent of the underlying structural pathology, with the severity of the clinical symptoms, and with disease duration (Stam, 2014).

There is an increasing demand for functional markers of early cognitive decline to identify patient populations that have an increased risk of developing dementia as these individuals are the best applicants for therapeutic intervention. Previous EEG studies revealed potential spectral and functional connectivity (FC) biomarkers that are able to predict the future progression of cognitive decline (Moretti et al., 2011; Toth et al., 2014; Mazaheri et al., 2018; Sharma et al., 2019).

The assessment of functional connectivity and network topology can provide an integrative approach that can reflect progressive brain dysfunction in MCI and AD (Pievani et al., 2011; Stam, 2014; Hallett et al., 2020). Moreover, graph theory approach could provide a general language that could help us to understand how cortical atrophy and functional disruptions are linked together in the pathological processes of AD (; Stam, 2014; Miraglia et al., 2017; ) and to discover novel early diagnostic and predictive neurophysiological markers (Rossini et al., 2016; Horvath et al., 2018).

There is a considerable amount of literature reporting decreased resting-state functional connectivity in MCI and AD in the alpha- and beta frequency range (Stam et al., 2003; Stam, 2014; ; Koelewijn et al., 2017; Horvath et al., 2018; Núñez et al., 2019; ). Changes in memory task-related functional connectivity are much less investigated and former studies reported mixed results (Hogan et al., 2003; Pijnenburg et al., 2004; Jiang and Zheng, 2006; Hou et al., 2018). The conflicting results might be partly explained by differences in the diagnostic criteria of the study groups (clinical or biomarker-based, MCI or AD patients), sample size, and the choice of functional connectivity measure, some of which are not corrected for the effect of volume conduction, which might influence previous results (; Herreras, 2016).

Regarding the overall network structure, previous studies observed a progressive derangement of brain organization during the disease course causing a deviation from the optimal small-world architecture to a more random type configuration leading to a less efficient information transfer during resting state (; Stam et al., 2009; Stam, 2014; Wei et al., 2015; Miraglia et al., 2017), and cognitive tasks (Wei et al., 2015; ), firstly affecting alpha-band networks in MCI (Miraglia et al., 2017).

Former studies highlighted the role of hubs in network disturbances in MCI and AD (Stam, 2014), which are nodes with high values of relative importance—such as node degree or betweenness centrality—and take a central role in network organization by facilitating the optimal flow within healthy brain networks (van den Heuvel and Sporns, 2013; Stam, 2014). Hub regions have been found especially vulnerable in AD (Stam et al., 2009; ; ; Tijms et al., 2013; ; Stam, 2014; Miraglia et al., 2017; Yu et al., 2017) and disruption of the global network structure in AD has been explained by the overload and failure of hub nodes (; Stam, 2014). Throughout the disease progression neural activity, functional connectivity, and hub activity follow an inverted U shape: increasing in early MCI, followed by a decrease in late MCI and AD ().

From a network perspective, visuospatial memory in MCI is an area of particular interest, as neuronal networks associated with this cognitive function are particularly affected by the neuropathological process of AD (Pievani et al., 2011), especially frontoparietal and frontotemporal connections (). Moreover, using memory tasks enhances EEG abnormalities related to MCI and improves the classification accuracy of healthy subjects and patients (van der Hiele et al., 2007a; San-Martin et al., 2021). Therefore we applied a computerized implementation of a visuospatial memory task in the current study.

Our study aimed to analyze EEG functional connectivity and network differences in the alpha and beta frequency band during memory maintenance between MCI patients and healthy elderly with subjective memory complaints.

Former studies reported decreased alpha and beta-band AEC-c in AD (Koelewijn et al., 2017; Núñez et al., 2019; ), therefore we hypothesized a decreased alpha- and beta-band functional connectivity in MCI patients and we expected that the memory load-related modulation of global functional connectivity will be less prominent in the MCI patients than the control subjects, since their reduced available cognitive capacity.

In accordance with the early increase of network integration suggested by the “hub overload and failure” framework (Stam, 2014) and based on previous MST network studies (; Lopez et al., 2017; Wang et al., 2018) we hypothesized a more centralized network topology in MCI patients. As hub nodes are exposed to an increased traffic load in a more centralized network, this transition might lead to the overload and subsequent failure of these hub nodes and the disturbance of the modular system of the network (Stam, 2014). Therefore, the shift to a more integrated network configuration might reflect the increased vulnerability of brain networks in MCI.

Materials and Methods

Participants and Clinical Measures

The study was carried out in the Department of Psychiatry and Psychotherapy, Semmelweis University, Budapest, Hungary. EEG was recorded from 17 MCI patients and 20 healthy control participants during a visuospatial memory task. Among them, structural MRI data of 13 MCI patient and 13 control participant and diffusion-weighted MRI (DW-MRI) data of 10 MCI patient and 17 control participant was available (10 MCI patient and 13 healthy control subject had both structural and functional MRI data). Participants had subjective memory complaints and applied to take part in a cognitive training program announced among general practitioners and in a Retirement Home (The study is registered at ClinicalTrials.gov, the identifier is “NCT02310620”). Every participant underwent a regular psychiatric assessment to evaluate possible excluding comorbidity. After that, cognitive functions were assessed with neuropsychological tests to specify the diagnosis [Addenbrooke’s Cognitive Examination (ACE), Rey Auditory Verbal Learning Test (RAVLT), Trail Making Test (TMT)]. Participants were not financially compensated for their participation but received a detailed written feedback on their performance on the neuropsychological tests.

The diagnostic procedure of MCI was based on the Petersen criteria (Petersen, 2004), including subjective memory complaints corroborated by an informant, preserved everyday activities, memory impairment based on a standard neuropsychological test, preserved global cognitive functions, and the exclusion of dementia. For the detailed assessment of memory impairment, we applied the Rey Auditory Verbal Learning Test (RAVLT) (Strauss et al., 2006). Attention, executive functions, and cognitive flexibility were examined with the Trail Making Test (TMT) Part A and Part B (Tombaugh, 2004; Strauss et al., 2006), global cognitive performance was estimated with the Addenbrooke’s Cognitive Examination (ACE) (Mathuranath et al., 2000). For the differentiation between MCI and healthy controls, we applied a cut-off score of 1 SD under population mean standardized for age and gender/education in these neuropsychological tests. Participants, who scored under the cut-off value in the delayed recall subscore or the total score of RAVLT or the TMT Part B or the ACE, were put into the MCI group. Subjects with dementia were excluded based on cognitive impairment according to the Mini-Mental State Examination (MMSE) scores standardized for age and education (Strauss et al., 2006) and on the loss of ability to perform activities of daily living. The Geriatric Depression Scale (GDS) was used to assess depressive symptoms (Yesavage, 1988), while anxiety symptoms were measured by the Spielberger State-Trait Anxiety Inventory (STAI) (Spielberger et al., 1970). Exclusion criteria were history of head trauma with loss of consciousness, prior CNS infection, epileptic seizure, clinically significant brain lesions (stroke, severe periventricular white matter disease, clinically significant white matter infarcts), multiple sclerosis or other demyelinating disorders, hydrocephalus, untreated vitamin B12 deficiency, untreated hypothyroidism, syphilis or HIV infection, mental retardation, major depression, schizophrenia, other acute psychiatric disorder, electroconvulsive therapy, renal insufficiency, liver disease, significant systemic medical illness, alcohol, or substance use dependency. Demographic and neuropsychological data are summarized in Table 1.

TABLE 1

control (n = 20)MCI (n = 17)p-value
Age [Mean (SD)]65.2 (6.9)69.9 (6.5)p = 0.04
Educationa15%/15%/70%18%/18%/65%n.s.*
Gender (female)70%41.2%n.s.*
Rey Auditory Verbal Learning Test 1–5 sumb54.3 (7.8)40.0 (11.3)p < 0.0001
Rey Auditory Verbal Learning Test delayed recallc (MCI: n = 13)11.4 (2.6)7.2 (4.4)p = 0.007
ACE total scored94.9 (2.9)86.2 (8.3)p = 0.0006
ACE VL/OM-ratioe2.6 (0.3)2.8 (0.6)n.s.*
Mini mental state examination total scoref29 (1.2)27.9 (1.4)p = 0.02
Trail Making Test Part Ag34.9 (10.8)70.6 (52.9)p = 0.006
Trail Making Test Part Bg (MCI: n = 16)69.0 (22.7)143.5 (69.2)p < 0.0001
Geriatric Depression Scale scoreh (Control: n = 19)3.6 (2.9)4.3 (3.5)n.s.*
STAI scorei39.4 (11.0)35.8 (9.4)n.s.*

Demographic data and results of basic neuropsychological tests.

MCI, Mild cognitive impairment; ACE, Addenbrooke’s cognitive examination; STAI, State-trait anxiety inventory.

aParticipants were categorized into three education groups: 1 = less than 12 years; 2 = high school graduation (12 years education); 3 = more than 12 years of education.

bSum of all words in the first five trials. The maximum score is 75.

cThe maximum score is 15.

dThe maximum score is 100.

eVL/OM: verbal fluency and language points/orientation and delayed recall ratio can be defined based on ACE. A result below 2.2 indicate frontotemporal dementia and a result over 3.2 indicate Alzheimer’s disease.

fThe maximum score is 30.

gTime needed for completing the task in seconds.

hThe maximum score is 15.

iThe maximum score is 80.

*n.s. (not significant) = p > 0.05.

EEG Paradigm and Procedures

EEG examinations were carried out on weekdays between 10 a.m. and 4 p.m. Participants were seated in a dimly lit, sound-attenuated room. All participants had normal or corrected-to-normal vision.

To measure visuospatial memory, during the EEG recording participants performed an implementation of the PAL test used in several neuropsychological test batteries (Sirály et al., 2013). White windows and colored shapes sized 2.65 cm × 2.65 cm were presented as stimuli on a computer screen at approximately 50 cm distance with Presentation 13.0 software (Neurobehavioral Systems, Inc.; Albany, CA). At the onset of each trial, eight blank windows appeared on the screen for 1,500 ms. After that, two, three, or four random windows opened up sequentially for 1,500 ms with abstract shapes shown in them, separated by a fixation cross for 450–500 ms. Meanwhile, other windows remain blank depending on the difficulty level. For the retention period, a fixation cross appeared for 3,800–4,000 ms. During the retrieval period, the previously shown shapes reappeared in the windows, and participants were instructed to indicate by clicking with the mouse (yes-right/no-left) whether the shapes popped up in the same positions they saw them before (Figure 1). The test consisted of 72 trials in total (32 two-item, 24 three-item, 16 four-item). The response assignment was counterbalanced across trials. Efficiency was measured by response accuracy.

FIGURE 1

It was carefully monitored that the participants understood the instructions and stayed alert during the session to bypass the possible distorting effect of extended eye closure on the EEG activity, especially in the alpha frequency range (). For the same purpose, participants completed the task in three parts separated by a 3-min rest period.

EEG Recording and Processing

EEG was recorded from DC with a low-pass filter at 100 Hz using a high-density 128-channel BioSemi ActiveTwo amplifier (Metting van Rijn et al., 1990). Electrode caps had an equidistant layout and covered the whole head according to the Biosemi equiradial montage. Eye movements were monitored with EOG electrodes placed below the left and above the right external canthi. Data were digitized at a sampling rate of 1,024 Hz. Built-in and self-developed functions as well as the freeware EEGLAB toolbox () in the Matlab (MathWorks, Natick, MA) development environment was used for subsequent off-line data analyses. EEG was re-referenced to the common average reference and filtered off-line between 0.5 and 45 Hz using zero-phase shift forward, and reverse IIR Butterworth filter. As four channels (P2, FT7h, P7, P9) were exceptionally noisy across multiple subjects, they were removed from the recordings for all subjects prior to the analysis. Epochs from 500 ms pre-stimulus to 4,500 ms post-stimulus for the retention period were extracted from the continuous EEG. Removal of muscle, blinking, and eye movement artifacts (detected by EOG) were performed by the Multiple Artifact Rejection Algorithm (MARA), a machine-learning algorithm that evaluates the ICA- (Independent Component Analysis) derived components (Winkler et al., 2011, 2014), Furthermore, epochs with a voltage exceeding ±100 μV on any channel were rejected from the analysis. After artifact rejection, the average number of trials in the control group and in the MCI group were 71.7 (SD = 0.7) and 71.1 (SD = 2.3) for the retention condition, respectively.

EEG Data Analysis

After artifact rejection, epochs of 4,000 ms duration of the retention period (from 200 ms pre-stimulus to 3,800 ms post-stimulus, sampling rate 1,024 Hz, 4,096 time points) were extracted from the EEG recording, which, based on previous studies, we assumed to be sufficient to measure oscillatory activity in the alpha and beta frequency band (). EEG connectivity analyses were performed with open-access software BrainWave (version 0.9.152.12.26; available at http://home.kpn.nl/stam7883/brainwave.html). Functional connectivity between EEG channels was analyzed by measuring the amplitude envelope correlation with leakage correction (AEC-c) calculated for all EEG epochs of each subject, after having band-pass filtered the EEG time-series in the alpha (8–13 Hz) and beta (13–30 Hz) frequency band. The amplitude envelope correlation (AEC) measures the linear correlations of the envelopes of the band-pass filtered and Hilbert-transformed signals (). The leakage-corrected version of the AEC (Hipp et al., 2012) uses a pair-wise symmetric orthogonalization procedure before the calculations of the AEC to remove zero-lag correlation correlations that could be attributed to spurious connectivity caused by volume conduction. Therefore, it is considered a reliable measure of genuine functional connectivity (; Hipp et al., 2012; ; ). Connectivity metrics were averaged over epochs creating values for each electrode at the patient level. Global functional connectivity values were calculated by averaging the AEC-c of all electrodes.

We carried out a spectral analysis to assess whether the detected effects were solely driven by differences in spectral power or peak frequency. Relative power in alpha and beta frequency band and peak frequency (Hz; dominant frequency between 4 and 13 Hz) were calculated with the BrainWave software using Fast Fourier Transformation.

Graph-Theoretical Analysis

The graph-theoretical representation of the functional connectivity matrix was constructed by the Minimum Spanning Tree (MST), which is a simplified representation of the core network containing the strongest and most relevant “backbone” connections (Stam et al., 2014; Tewarie et al., 2015) that can reflect topological changes (Tewarie et al., 2015). Former studies pointed out that graph theoretical measures are dependent on network size and density, which can make the comparison across different groups and conditions by using conventional network analytical methods challenging (van Wijk et al., 2010; ; Stam et al., 2014). The MST calculation overcomes the bias of network density and degree without any additional normalization step by forming an acyclic subnetwork using the strongest available connections without forming loops and connecting all nodes with a fixed number [(number of nodes) - 1) of edges. MST graphs were generated for each participant, epoch for alpha and beta frequency band, based on the full connectivity matrix constructed from the AEC-c values obtained for each pair of electrodes. MST metrics were averaged over epochs for each subject.

Two extreme topologies of MST can be distinguished: a path-like and a star-like shape. In a path, all nodes are linked to exactly two other nodes, except the two nodes at the extremities of the tree. These nodes are connected to only one other node and are referred to as the leaves of the tree. In the case of a star shape, all but one node are linked to a central node (Stam et al., 2014). Between these two shapes, MST-s can have various configurations (Figure 2).

FIGURE 2

The diameter of the tree is the maximum number of edges between any two nodes of the network. Leaf fraction is the number of nodes with exactly one connection divided by the total number of nodes of the tree. Degree refers to the number of edges connected to a node. Betweenness centrality (BC) of a node refers to the normalized fraction of all paths connecting two nodes that pass through the selected node, and it characterizes the “hubness” of the node within the network. The eccentricity of a node denotes the longest shortest path to any other node in the MST. Degree divergence (kappa—κ) measures the broadness of the degree distribution, which shows high value in networks with high-degree hubs, and it is related to the resilience of the network against attacks. In an MST the most efficient communication can be achieved in a star-like configuration, as it has the shortest possible average path length between two arbitrary nodes, however, in this case, the central node might easily be overloaded. This trade-off between large-scale integration and the overload of central nodes is quantified by the tree hierarchy. The optimal MST topology balances efficiency and node load.

Global and node-specific parameters were computed with the Brainwave software, based on the measures described by previous studies (Stam et al., 2014; Tewarie et al., 2015), summarized in Table 2. Degree, betweenness centrality, and eccentricity were calculated for each node separately, and the maximum degree, maximum BC, and mean eccentricity were included in the statistical analysis as global characteristics of the MST. Global MST network parameters were averaged across epochs.

TABLE 2

MeasureExplanationFormula
Nodes (N)Number of nodes
Links (M)Number of links/maximum leaf number
Degree (k)Number of links for a given node. Nodes with a high degree may be considered hubs. We used the maximum degree to characterize the strength of the most important node of the network.
Leaf fraction (Lf)Fraction of leaf nodes (L) in the MST where a leaf node is defined as a node with only one connection. It describes to what extent the network has a central organization. A high leaf fraction indicates, that communication is largely dependent on hub nodes.Lf=L/M
DiameterLongest distance between any two nodes in an MST, normalized by the total number of connections. In a network with a low diameter, information is efficiently processed between remote brain regions. The diameter is also related to the leaf number: the value of the diameter decreases when the leaf number increases.D=d/M
EccentricityLongest shortest path from a reference node to any other node in the MST. Eccentricity is low if the node is located in the center of the tree. Eccentricity of the network describes how efficient information is communicated from the least central node.
Betweenness centrality (BC)Fraction of all shortest paths that pass through a particular node. BC ranges between 0 (leaf node) and 1 (central node in a star-like network). Nodes with a high BC are considered hub nodes based on their importance for global communication in the network. The BC of the tree was characterized by the maximum BC value, which describes the importance of the most central node and it is a measure of central network organization. ρih is the number of shortest paths between h and j, and ρih(i) is the number of shortest paths between h and j that pass through i
Degree divergence (κ)Measure of the broadness of the degree distribution. Related to resilience against attacks, epidemic spreading and the synchronizability of complex networks
Tree hierarchy (TH)Quantifies the trade-off between large scale integration in the MST and the overload of central nodes It characterizes the hypothesized optimal topology of brain network organization, where information is transferred between brain regions in the fewest possible steps, while preventing information overload of central brain regions.

Explanation of concepts and terminology based on Tewarie et al. (2015) and Van Dellen et al. (2015).

MR Image Acquisition and Processing and Diffusion Tensor Fitting

The obtained structural gray matter volumetric (cortical thickness and subcortical brain structure volumes) and the diffusion-weighted data were previously published by our study group (; ).

Participants underwent a routine brain MR examination, producing high-resolution anatomical images used for analysis. Image acquisitions were made at the MR Research Center, Semmelweis University, Budapest on a 3 Tesla Philips Achieva clinical MRI scanner equipped with an 8-channel SENSE head coil. High resolution, whole-brain anatomical images were obtained using a T1 weighted 3-dimensional spoiled gradient echo (T1W 3D TFE) sequence. 180 contiguous slices were acquired from each subject with the following imaging parameters: TR = 9.7 ms; TE = 4.6 ms; flip angle = 8°; FOV of 240 mm × 240 mm; voxel size of 1.0 × 1.0 × 1.0 mm. Brain DW-MRI images were collected with a single shot SE-EPI sequence, with b = 800 s/mm2 diffusion weighting in 32 directions and one b = 0 image. In-plane resolution was 1.67 × 1.67 mm; whole-brain coverage was achieved with 70 consecutive, 2 mm thick axial slices; TR = 9,660 ms repetition time, TE = 75.6 ms echo time, and 90° flip angle was used; the total acquisition time was 8:32 min.

Cortical reconstruction, volumetric segmentation and parcellation of the MRI data into standardized region of interest (ROIs) were performed automatically by Freesurfer 5.3 image analysis suite1 (see details in ), however, segmentation and cortical models were checked and corrected manually on each subject. Volumetric measurements were normalized by dividing by the intracranial volume (ICV) also computed during the Freesurfer pipeline, while cortical thickness measurements were included in the analysis without further normalization based on previous results (Westman et al., 2013).

DWI data were preprocessed using the Matlab-based ExploreDTI software package (Leemans et al., 2009). Processing steps included coordinate system transformation, rigid body transformations for correcting subject motion, non-rigid transformations for correcting susceptibility-related and EPI-induced distortions, with the local rotation of the b-matrix (the diffusion weighting directions) to avoid angular inaccuracies (Leemans and Jones, 2009). The high-resolution T1-weighted images were used as templates for registration to correct the distortions inherent to the EPI-acquisition method (Jezzard et al., 1998); thereby DW-images were spatially aligned to the T1W images. After tensor fitting, using the RESTORE (Robust Estimation of Tensors by Outlier Rejection) () algorithm, two voxel-wise DTI-measures, fractional anisotropy (FA) and mean diffusivity (MD) (Pierpaoli and Basser, 1996; ; ) were calculated from the tensor eigenvalues, following their well-established definitions, to be used in voxel-level and ROI-based analyses (see for details on tensor fitting and DTI scalar calculations).

Statistical Analysis

Demographic characteristics, results of the neuropsychological tests, and response accuracy of the study groups were compared with independent samples t-tests, Mann-Whitney U tests, or χ2 tests where appropriate. Normal distribution of variables was tested using the Kolmogorov–Smirnov test.

Group comparisons of global functional connectivity and MST metrics were performed on the EEG from three levels of memory load conditions (two-item, three-item, four-item), while we used the average of these conditions for the correlational analysis with the size of medial temporal lobe structures and hippocampal fiber integrity.

Functional connectivity and network parameters of the two study groups were tested by two-way analysis of covariance (ANCOVA) of the study group (HC vs. MCI) × memory load (two- vs. three vs. four-item sequence). All the main effects including age as a covariate and two-way interactions were included in the ANCOVA model. Statistical significance was determined at p < 0.05.

Post-hoc pairwise contrasts were conducted to investigate the interactions. Since between-group comparisons were evaluated over three levels of memory load, Hochberg correction for multiple comparisons was applied to the post-hoc contrasts (Hochberg, 1988; Hochberg and Benjamini, 1990). To characterize the magnitude of the reported effects we reported the values of effect size (Cohen’s d) ().

Structural and DW-MRI results were derived from previously published parts of our study (Siraly et al., 2015; ; ). We followed a ROI-based approach and assessed the association between functional connectivity and early-stage medial temporal lobe atrophy and hippocampal fiber integrity as these are important early markers of MCI (Márquez and Yassa, 2019). As the MRI results of some participants were outlier values, we applied the Spearman correlation in the analyses which is robust against the effect of outliers.

Results

Demographic and Neuropsychological Characteristics

In total, 17 MCI patients (mean age 69.9 ± 6.5 years; 7 females) and 20 healthy control participants (mean age 65.2 ± 6.9 years; 14 females) were included in the study. Groups did not differ with regard to gender, level of education, depressive symptoms (GDS score), and anxiety symptoms (STAI-score). However, MCI patients were older than the control participants, therefore statistical tests were corrected for age as a covariate. Furthermore, patients with MCI had a significantly lower score on the neuropsychological tests (ACE, MMSE, RAVLT, MMSE) than the control participants (Table 1).

Behavioral Results

In the PAL task response accuracy of the MCI patients showed a trend level decrease compared to the control group (MCI: mean = 77.2% SD = 21.2, HC: mean 88.4 = % SD = 7.2, U = 106.5, Z = 1.9, p = 0.05, Cohen’s d = 0.8). The control group had a significantly lower score in the high memory load (four-item) condition compared to the low memory load (two-item) (Z = –3.4, p = 0.0006) and to the medium memory load (three-item) condition (Z = 2.9, p = 0.0041), while in the MCI group no significant memory load-related differences were observed in response accuracy.

Functional Connectivity in the Alpha Band

During the retention period of the PAL test memory load had a significant modulatory effect on alpha AEC-c [F(2, 34) = 5.92, p = 0.006] (Figure 3). Furthermore, a trend-level interaction of group and memory load was observed [F(2, 34) = 3.03, p = 0.06]. The mean alpha AEC-c and the topography of average connectedness (i.e., mean AEC-c of each electrode) are shown in Figure 3. Post-hoc analysis of this interaction revealed, that the memory load-related modulation of AEC-c followed different dynamics in the two study groups: in the control group compared to the low memory load condition (two-item), a significantly increased mean AEC-c was observable in the medium memory load condition (three-item; t = 2.59, df = 34, p = 0.01, Cohen’s d = 0.4) and in the high memory load condition (four-item; t = 2.88, df = 34, p = 0.007, Cohen’s d = 0.4) and these memory load-related differences remained significant after correction for multiple comparisons.

FIGURE 3

In contrast to that, the MCI group showed a significantly increased mean AEC-c in the medium memory load condition compared to low memory load (t = 2.28, df = 34, p = 0.03, Cohen’s d = 0.3), while in the high memory load condition a significantly diminished mean functional connectivity was observable compared to medium memory load (t = 2.5, df = 34, p = 0.02, Cohen’s d = 0.3), however, these differences became trend level after correction for multiple comparisons (corrected p = 0.06 and 0.05, respectively). Study group and age did not have a significant effect on alpha functional connectivity (p > 0.05).

Functional Connectivity in the Beta Band

During the retention period of the PAL test study group, memory load and age did not have a significant effect on beta AEC-c (p > 0.05). The mean beta AEC-c and the topography of average connectedness (i.e., mean AEC-c of each electrode) are shown in Figure 4. Interaction of study group and memory load was not significant, however, the post-hoc analysis revealed that in the control group mean beta functional connectivity in the high memory load condition was significantly increased compared to the low memory load condition (t = 2.82, df = 34, p = 0.008, Cohen’s d = 0.4), which remained significant after correction for multiple comparisons, while in the MCI group no memory load-related differences were observable (Figure 4).

FIGURE 4

Correlational Analysis of Alpha and Beta Functional Connectivity and the Size and Fiber Integrity of the Medial Temporal Lobe Structures

Correlational analysis of mean functional connectivity averaged over all conditions and structural and DW-MRI results of medial temporal lobe structures (relative hippocampal volume, cortical thickness of the parahippocampal and the entorhinal gyrus, mean diffusivity (MD), and fractional anisotropy (FA) of the right and left cingulum—hippocampal subdivision) was performed on the entire sample. Mean alpha and beta AEC-c showed a significant positive correlation with the total relative hippocampal volume (alpha AEC-c: Spearman r = 0.47, p = 0.02, beta AEC-c: Spearman r = 0.54, p = 0.004), and with the cortical thickness of the parahippocampal gyrus (alpha AEC-c: Spearman r = 0.40, p = 0.04, beta AEC-c: Spearman r = 0.48, p = 0.01) and a significant negative correlation with the mean diffusivity of the right cingulum—hippocampal subdivision (alpha AEC-c: Spearman r = –0.41, p = 0.03, beta AEC-c: Spearman r = –0.50, p = 0.008). Furthermore, mean beta AEC-c correlated significantly with the cortical thickness of the entorhinal gyrus (beta AEC-c: Spearman r = 0.44, p = 0.02) (Figure 5).

FIGURE 5

Correlations of the mean beta AEC-c and structural MRI results were driven by the MCI group (relative hippocampal volume: Spearman r = 0.72, p = 0.008, parahippocampal gyrus Spearman r = 0.60, p = 0.04, entorhinal gyrus Spearman r = 0.70, p = 0.01). Moreover, correlations between the mean alpha and beta AEC-c and the mean diffusivity of the right hippocampal cingulum were driven by the MCI group (alpha AEC-c: Spearman r = –0.83, p = 0.003, beta AEC-c: Spearman r = –0.67, p = 0.03). Detailed results of the correlational analysis with stratified by diagnosis can be found in Supplementary Table 1.

Spectral Analysis

Our results showed that while study group [F(1, 34) = 0.02, p = 0.88] and age [F(1, 34) = 1.07, p = 0.30] did not have a significant effect on relative alpha power, memory load had a modulatory effect on relative alpha power [F(2, 34) = 4.04, p = 0.03]. Interaction of study group and memory load showed a trend level effect [F(2, 34) = 3.13, p = 0.06]. The post-hoc analysis revealed that in the control group the relative alpha power in the high memory load condition was significantly increased compared to the low memory load condition (t = 3.69, df = 34, p = 0.0006, Cohen’s d = 0.3), which remained significant after correction for multiple comparisons.

In the beta band neither study group [F(1, 34) = 1.26, p = 0.27] nor memory load [F(2, 34) = 0.41, p = 0.67] or age [F(1, 34) = 0.90, p = 0.35] had a significant effect on relative beta power. Interaction of study group and memory load showed a trend level effect [F(2, 34) = 0.06, p = 0.94]. The post-hoc analysis revealed no significant effects.

Furthermore, study groups did not have a significantly different peak frequency [F(1, 34) = 0.21, p = 0.65], however memory load had a significant modulatory effect on the peak frequency values [F(2, 34) = 6.44, p = 0.043]. Interaction of group and memory load showed a trend level effect [F(2, 34) = 2.97, p = 0.07], however, the post-hoc analysis revealed that in the control group the peak frequency in the high memory load condition was significantly increased compared to the low memory load condition (t = 3.87, df = 34, p = 0.0005, Cohen’s d = 0.3), which remained significant after correction for multiple comparisons. Furthermore, there was no significant difference between the two study groups regarding the mean peak frequency values (averaged over memory loads) (MCI: mean = 8.5 Hz, SD = 1.4, HC: mean = 8.2 Hz, SD = 1.3, t = 0.75, df = 35, p = 0.46). Distribution of mean relative power in the alpha and beta frequency band and peak frequency by study groups can be found in Supplementary Figure 1.

Minimum Spanning Tree Parameters in the Alpha Band

The network analysis (calculated over all memory load conditions) indicated that the MCI group had a significantly decreased MST diameter compared to the control group [F(1, 34) = 5.36, p = 0.03]. Furthermore, a decreased eccentricity was observed in the MCI group [F(1, 34) = 4.85, p = 0.03]. However, age also had a significant mean effect on these parameters [F(1, 34) = 4.64, p = 0.04 and F(1, 34) = 4.14, p = 0.05, respectively].

The MCI group had a significantly increased maximum MST degree [F(1, 34) = 5.69, p = 0.02], degree divergence [F(1, 34) = 6.12, p = 0.02], and maximum betweenness centrality [F(1, 34) = 7.37, p = 0.01] (Figure 6) compared to the control group. Furthermore, memory load had a significant modulatory effect on betweenness centrality [F(2, 34) = 3.53, p = 0.04], indicating a significantly increased BC in the medium memory load condition compared to the low memory load condition (t = 2.6, df = 34, p = 0.01). Leaf fraction and tree hierarchy did not differ significantly in the two groups. Group-average MSTs of the two study groups (average of all levels of memory load) are shown in Figure 6, represented in sensor space. The central hub (the node with the most connections) was the right temporal electrode T8 in both study groups. Detailed results of the MST analysis are summarized in Supplementary Table 2.

FIGURE 6

Minimum Spanning Tree Parameters in the Beta Band

The network analysis (calculated over all memory load conditions) indicated that the MCI group had a significantly decreased MST diameter compared to the control group [F(1, 34) = 4.58, p = 0.04]. Meanwhile, a decreased eccentricity was observed in the MCI group [F(1, 34) = 5.62, p = 0.02].

The MCI group had a significantly increased maximum MST degree [F(1, 34) = 7.55, p = 0.01], degree divergence [F(1, 34) = 7.15, p = 0.01], and maximum betweenness centrality [F(1, 34) = 6.95, p = 0.01] (Figure 6) compared to the control group. However, age also had a significant mean effect on maximum MST degree and degree divergence [F(1, 34) = 5.2, p = 0.03 and F(1, 34) = 5.44, p = 0.03, respectively]. There was no significant difference in leaf fraction and tree hierarchy. Group-average MSTs of the two study groups (average of all levels of memory load) are shown in Figure 6, represented in sensor space. The central hub (the node with the most connections) in the control group was the right temporal electrode T8, while in the MCI group it was the left frontal-temporal electrode FT7. Detailed results of the MST analysis are summarized in Supplementary Table 2.

Discussion

In this study, we aimed to examine functional connectivity and network structure during memory maintenance in MCI patients and healthy controls.

We used the orthogonalized Amplitude Envelope Correction (AEC-c) for the measurement of functional connectivity, which corrects to the effect of volume conduction, is independent of relative power (), gives reliable estimates of the underlying network topology (Lai et al., 2018), and has been found the most sensitive measure of functional connectivity in the alpha and beta frequency bands (Hipp et al., 2012). Moreover, the AEC-c produced the most reproducible and valid results in AD compared with other measures of functional connectivity, such as mean global coherence (Coh), imaginary coherence (iCoh), phase locking value (PLV), phase lag index (PLI), weighted PLI (wPLI), and the AEC without leakage-correction ().

Alpha and beta-band oscillatory synchrony play an important role in cognitive tasks by mediating top-down directed influences on task-relevant cortical areas (). Furthermore, alpha synchronization is instrumental in controlling the flow of information especially in the thalamo-cortical and cortico-cortical networks and in the timing of working memory-related processing by the modulation of neural excitability (Klimesch, 1999, 2012; Pfurtscheller and Lopes da Silva, 1999; Palva and Palva, 2011; Wang et al., 2014; Miraglia et al., 2016; Wianda and Ross, 2019) while beta oscillations have been linked to the active maintenance of newly acquired information for further task requirements (Onton et al., 2005; ; Missonnier et al., 2007; ; ) and to the facilitation of long-range connections in cortical networks (Kopell et al., 2000; Varela et al., 2001; ; ; ; Kilavik et al., 2013), especially during attentional and memory processes () and endogenous content reactivation (Spitzer and Haegens, 2017). However, many of these previous studies measured functional connectivity by amplitude correlation, whose exact relation to the correlation of amplitude envelopes (which has been used as the measure of connectivity in the present study) and the exact mechanism underlying communication by amplitude envelopes are still not fully understood.

According to our results, memory load modulated the mean functional connectivity in the alpha band, but it had a different modulatory effect in the two study groups. In the control group, increasing task difficulty enhanced the mean functional connectivity. In contrast to that, after an initial increase of the mean AEC-c in the medium memory load condition the MCI group showed significantly diminished functional connectivity in the high memory load condition. The control group showed a similar memory load-related increase in the mean AEC-c in the beta band, while in the MCI group this modulation was absent.

Cognitive load-related increase of alpha and beta band functional connectivity is in line with previous studies which observed the same phenomenon in healthy subjects as well as in MCI and AD patients in alpha and beta band (Pijnenburg et al., 2004; Palva et al., 2010; Wianda and Ross, 2019) and in broadband (Jiang and Zheng, 2006). However, other studies observed the diminishment of task-related increase of coherence of alpha oscillations in AD patients in the dementia phase (based on clinical diagnosis) (Hidasi et al., 2007).

Our results suggest that the AEC-c is a sensitive measure that can follow the modulation of functional connectivity by cognitive demand, especially in the alpha frequency band, which has been associated with attentional functions (Palva et al., 2010; Sato et al., 2018; Marzetti et al., 2019).

We observed that the initial increase of alpha-band functional connectivity in the medium memory load condition was followed by a decrease in the high memory load condition in the MCI group. As former studies linked the increase of alpha-band functional connectivity to enhanced cognitive demand (Pijnenburg et al., 2004; Palva et al., 2010; Wianda and Ross, 2019), and we found a similar modulation in the control group, we hypothesized that the initial increase of alpha connectivity in the MCI group in the medium memory load condition indicates the increased utilization of working memory. However, due to the limited cognitive reserve, MCI patients are unable to act likewise in the high memory load condition as task difficulty exceeds their cognitive capacity. Therefore, the reduction of alpha functional connectivity in the high memory load condition might indicate the reduced cognitive reserve and the impairment of working memory maintenance in MCI, although this was not reflected by a decrease in task performance.

In the beta band, while increasing memory load enhanced functional connectivity in the control group, the MCI group did not show memory load-related modulation. This might indicate a more extensive failure of working memory maintenance in MCI in the beta frequency band. Another possible explanation might be the general “slowing” of EEG in MCI (), namely that task-related dynamics of higher frequency bands got shifted toward lower frequency bands in MCI. The study group did not have a significant modulatory effect on the mean (whole head) AEC-c, which is in line with a previous study on alpha coherence during a working memory task which did not find significant differences between aMCI patients and healthy older adults (van der Hiele et al., 2007b). However, in MCI patients some studies found increased alpha and beta synchronization (Pijnenburg et al., 2004; Jiang, 2005; Jiang and Zheng, 2006), which has been attributed to compensatory mechanisms (; ).

The relatively small sample size of our study may be a factor contributing to the observed lack of between-group differences in mean functional connectivity. A different potential explanation might be that MCI is characterized by increased and decreased functional connectivity in different cortical regions simultaneously (Lopez et al., 2017), and therefore during averaging these changes might be smoothed out. Moreover, most of the previous studies assessed eyes-closed resting-state recordings, while we analyzed eyes-open task-related EEG, which might be another influencing factor.

Alpha and beta AEC-c showed a significant positive correlation with the size of medial temporal lobe structures and a significant negative correlation with the mean diffusivity (MD, a scalar measure of overall water diffusion) of the right hippocampal cingulum in the entire sample. MD increases in the presence of tissue damage and is typically used to assess the microstructural integrity of gray matter (Stebbins and Murphy, 2009). Elevation of hippocampal, parahippocampal, and temporal lobe MD are considered as important early markers of neuronal loss and disruption of myelin sheaths in MCI and AD (Kantarci et al., 2001, 2005; ; Ray et al., 2006; Stebbins and Murphy, 2009; Zhang et al., 2014).

Therefore, our results suggest that functional connectivity and specifically the AEC-c in the alpha and beta frequency band can reflect the subtle medial temporal lobe atrophy and the disruption of hippocampal fiber integrity in the earliest stages of cognitive decline. This is also corroborated by the fact, that these correlations were driven by the MCI subjects, who had a more pronounced hippocampal degeneration compared to the control group.

There is some evidence that changes in MD are more typical in MCI whereas as changes in MD and fractional anisotropy (FA, a measure of the directionality of diffusion) are more typical in AD (Rogalski et al., 2009; Stebbins and Murphy, 2009), which might be the reason why hippocampal FA did not show a significant correlation with mean functional connectivity. Altogether, our results are in line with previous DTI studies reporting a correlation between alpha-band functional connectivity and fiber tract integrity reduction in MCI and mild to moderate AD patients (Teipel et al., 2009; Vecchio et al., 2015).

The two study groups did not differ significantly regarding relative alpha and beta power and peak frequency, in contrast to former studies reporting a widespread decrease of alpha activity in the prefrontal, temporal, parietal, and occipital cortices during the n-back task (San-Martin et al., 2021). In the control group, we found a significant increase of alpha power in the high memory load compared to the low memory load condition in line with former studies (Jensen et al., 2002; Tuladhar et al., 2007; Palva et al., 2011), while this modulatory effect was absent in the MCI group. In the beta band, we did not detect a memory load-related modulation. However, in contrast to the functional connectivity results, the outcome of the power spectrum analysis did not show a load-related modulation of the MCI group in the alpha band (an initial increase followed by a decrease parallel with the enhancement of cognitive load). Therefore, we conclude that the detected differences in functional connectivity are not entirely the consequence of differences of spectral properties, although we cannot rule out the possibility that it might have an influence on the results, especially in the alpha frequency band, which might be a potential limitation.

We performed the network analysis by applying the MST approach, which provides an unbiased reconstruction of the critical backbone of the original network (Stam, 2014; Stam et al., 2014; Tewarie et al., 2015; Wang et al., 2018; Musaeus et al., 2019), and can capture the subtle changes of network topology in MCI more sensitively than traditional graph theoretical measures (Lopez et al., 2017).

We found a decreased MST diameter and eccentricity and increased maximum degree, degree divergence, and maximum betweenness centrality in the MCI group, suggesting a more centralized and integrated network topology compared to the control subjects both in the alpha as well as in the beta frequency band. Our results are in line with former studies, which reported increased BC values and node degree in MCI and AD patients (; Lopez et al., 2017).

The central hub (the node with the most connections) of the group-averaged MST network was the temporal electrode T8 in the alpha band in both study groups and in the beta band in the control group, while it was the left frontal-temporal electrode FT7 in the beta-band MST of the MCI group. The right superior temporal gyrus has been previously identified as an important hub region during working memory maintenance (Park et al., 2011) based on cross-frequency power correlations, however, as our analysis was performed on sensor-space data, we are not able to make precise assumptions about the exact spatial locations of the nodes of the networks. We found a left and slight frontal shift of hub location in the MCI group in the beta frequency band. Interestingly, a frontal shift of hub location (center of mass of BC) was observed with increasing disease severity in AD patients and has been attributed to the earlier impact of the disease pathology on the posterior regions ().

Previous studies interpreted the global network disturbances in MCI and AD by the “hub overload and failure” framework, which states that the initial disturbance of nodes leads to the abnormal rerouting of the information flow in the network to hub nodes with higher centrality leading to an increase of traffic load, and eventually to an overload and subsequent failure of these hub nodes. This stage might also coincide with the initial ascending phase in early MCI of the inverted U shape course of hub activity (). This initial increase of hub activity and the transition to a more integrated network topology might be part of a compensatory mechanism, but it might as well be a part of the degeneration process itself due to the early impairment of inhibitory neurons (disinhibition) (). Subsequently, in the chronic “hub failure” phase, these overloaded hubs break down and the rerouting is constrained locally to nodes with a lower level in the hierarchy in the remaining part of the network. This stage also corresponds to the descending phase of the trajectory of hub activity in late MCI and AD (). This will eventually lead to the disturbance of the modular system of the network (Stam, 2014).

The global network topology reflects this by an initial increase of centralization and a shift from local to global processing followed by a decrease of centrality (Stam, 2014). This transition has been confirmed by fMRI as well, where the MST of MCI patients showed a more star-like topology, while the MST of AD patients deviated toward a more line-like topology compared to healthy controls (Wang et al., 2018).

Our results suggest that brain networks of MCI patients show a transient shift to a more centralized, star-like topology to compensate for the initial impairments in accordance with the “hub overload” stage, and complement former EEG studies, which reported the deviation of the network topology from the optimal small-world architecture to a more random type configuration (Wei et al., 2015) and the shifting of the MST toward a more decentralized, line-like structure of AD patients in the “hub failure” stage during resting state (Yu et al., 2016; Peraza et al., 2018; ) and cognitive tasks ().

Interestingly, while functional connectivity sensitively reflected changes in cognitive demand, MST network measures did not show significant memory load-related modulation except for maximum BC in the alpha band. This suggests that the AEC-c might be a more state-like attribute, which reflects cognitive demand, while MST network parameters are more trait-like characteristics of MCI and are less dependent on the actual cognitive state.

Limitations

The present study was limited by the small sample size and a slight age difference between groups, therefore statistical tests were corrected for age as a covariate. However, age had a significant effect on some of the network parameters, which limits the generalizability of our results. Moreover, the PAL test did not have an equal number of trials in the different difficulty levels, which might have influenced the signal-to-noise ratio of the EEG analysis.

Furthermore, we analyzed global functional connectivity to assess robust differences that could be considered as potential biomarkers of cognitive decline. However, a regional analysis focusing especially on the connectivity of the working memory network [prefrontal cortex, the parietal and temporal lobe, and task-specific posterior areas (Ranganath, 2006; )] could have provided a more detailed picture of the exact topological distribution of MCI-related differences.

Moreover, in this study, we did not analyze functional connectivity in the theta-band, as the AEC-c produces less reliable and reproducible results in the theta band in contrast to the alpha and beta-band and therefore, it has been suggested that for the assessment of theta-band functional connectivity phase-based measures (PLI) should be used instead of amplitude-based measures (). However, this might be a potential limitation since frontal midline theta activity is an important marker of working memory processing (Jensen and Tesche, 2002; ; Sauseng et al., 2010; Kardos et al., 2014).

Furthermore, we performed a scalp-level EEG analysis, which does not allow inferences in terms of underlying neuroanatomy as the location of EEG channels do not relate trivially to the location of the underlying sources, which is a further limitation. It has been suggested, that results derived from scalp-level EEG network should be interpreted cautiously, however, the AEC-c may allow for more reliable estimates of the underlying global network organization compared to metrics that do not correct for the effect of volume conduction (Lai et al., 2018).

Furthermore, while the diagnosis of MCI patients was based on a detailed clinical examination, cerebrospinal fluid biomarkers were not available during the diagnostic procedure. Therefore, AD as the underlying cause of the cognitive disturbance could not be fully proven. Finally, follow-up data is not yet available to examine the predictive value of functional connectivity and network structure in the conversion rate to dementia. Our study provides a cross-sectional view of the changes in functional connectivity and network topology during working memory maintenance in MCI, although further studies in AD biomarker-proven subjects and applying similar paradigms are required to verify our results.

Conclusion

Our results suggest that the AEC-c sensitively reflects cognitive load-related modulation and impairment of memory retention in MCI. Moreover, alpha and beta-band AEC-c showed a significant correlation with the size of medial temporal lobe structures and with the mean diffusivity of the right hippocampal cingulum, therefore, the AEC-c can reflect subtle medial temporal lobe atrophy and the disruption of hippocampal fiber integrity in the earliest stages of cognitive decline.

Furthermore, the MST network topology of the MCI group showed a more centralized and integrated configuration compared to the healthy control subject, which is in line with the “hub overload and failure” framework, and might be part of a compensatory mechanism or a consequence of neural disinhibition.

Therefore, the assessment of EEG functional connectivity and network structure in the alpha and beta frequency range may provide a useful complementary diagnostic tool for the early detection of cognitive impairment and might be a step toward establishing functional biomarkers (Sharma et al., 2019). However, future research applying similar paradigms is required to further develop and confirm these initial findings by using follow-up data to determine the predictive value of functional connectivity measures and network parameters for future conversion to dementia.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by the National Scientific and Ethical Committee, Budapest, Hungary. The patients/participants provided their written informed consent to participate in this study.

Author contributions

GC designed the study, wrote the protocol, and contributed to the writing of all sections. ZF participated in the execution of measurements, managed the literature searches, undertook the statistical analysis, prepared the figures, and wrote the first draft of the manuscript. AH participated in the execution of measurements and contributed to the writing of all sections. ZH, AG, and CS contributed to the conceptualization of the study and the writing of all sections. All authors reviewed the manuscript.

Funding

This study was supported by the “Ambient Assisted Living Joint Programme (AAL) – Call 2” grant (project identifier AAL_08-1-2011-0005 M3W) (http://www.aal-europe.eu/), by the ÚNKP - New National Excellence Program of the Ministry of Human Capacities, by the Janos Bolyai Research Scholarship of the Hungarian Academy of Sciences (BO/00249/17, bo_78_20_2020), by the Predoctoral Fellowship Program of Semmelweis University (EFOP-3.6.3-VEKOP-16-2017-00009), by the Higher Education Institutional Excellence Program of the Ministry of Human Capacities in Hungary, within the framework of the Neurology thematic program of Semmelweis University, by the EU Joint Programme - Neurodegenerative Disease Research (JPND, www.jpnd.eu) project (National Research, Development and Innovation, Hungary, 2019-2.1.7-ERA-NET-2020-00006), and by the Hungarian Scientific Research Fund 2019 of the National Research, Development and Innovation Office (PD- 132652).

Conflict of interest

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2021.680200/full#supplementary-material

Supplementary Figure 1

Distribution of mean relative power in the alpha and beta frequency band and peak frequency in the study groups.

Supplementary Table 1

Correlations between the temporal lobe structures and mean functional connectivity stratified by diagnosis.

Supplementary Table 2

MST metrics in the alpha and beta frequency band.

References

  • 1

    AlexanderA. L.HurleyS. A.SamsonovA. A.AdluruN.HosseinborA. P.MossahebiP.et al (2011). Characterization of cerebral white matter properties using quantitative magnetic resonance imaging stains.Brain Connect.1423446. 10.1089/brain.2011.0071

  • 2

    BabiloniC.LizioR.MarzanoN.CapotostoP.SoricelliA.TriggianiA. I.et al (2016). Brain neural synchronization and functional coupling in Alzheimer’s disease as revealed by resting state EEG rhythms.Int. J. Psychophysiol.10388102. 10.1016/j.ijpsycho.2015.02.008

  • 3

    BajoR.MaestúF.NevadoA.SanchoM.GutiérrezR.CampoP.et al (2010). Functional connectivity in mild cognitive impairment during a memory task: implications for the disconnection hypothesis.J. Alzheimer’s Dis.22183193. 10.3233/jad-2010-100177

  • 4

    BarryR. J.ClarkeA. R.JohnstoneS. J.MageeC. A.RushbyJ. A. (2007). EEG differences between eyes-closed and eyes-open resting conditions.Clin. Neurophysiol.11827652773. 10.1016/j.clinph.2007.07.028

  • 5

    BasserP. J.PierpaoliC. (2011). Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. 1996.J. Magn. Reson.213560570. 10.1016/j.jmr.2011.09.022

  • 6

    BenchenaneK.TiesingaP.BattagliaF. (2011). Oscillations in the prefrontal cortex: a gateway to memory and attention.Curr. Opin. Neurobiol.21475485. 10.1016/j.conb.2011.01.004

  • 7

    BirdC.ChanD.HartleyT.PijnenburgY.RossorM.BurgessN. (2010). Topographical short-term memory differentiates alzheimer’s disease from frontotemporal lobar degeneration.Hippocampus2011541169. 10.1002/hipo.20715

  • 8

    BlackwellA.SahakianB.VeseyR.SempleJ.RobbinsT.HodgesJ. (2004). Detecting dementia: novel neuropsychological markers of preclinical Alzheimer’s disease.Dementia Geriatric Cogn. Dis.174248.

  • 9

    BrielsC.SchoonhovenD.StamC.De WaalH.ScheltensP.GouwA. (2020). Reproducibility of EEG functional connectivity in Alzheimer’s disease.Alzheimer’s Res. Ther.12:68.

  • 10

    BrookesM. J.HaleJ. R.ZumerJ. M.StevensonC. M.FrancisS. T.BarnesG. R.et al (2011). Measuring functional connectivity using MEG: methodology and comparison with fcMRI.Neuroimage5610821104. 10.1016/j.neuroimage.2011.02.054

  • 11

    BrunsA.EckhornR.JokeitH.EbnerA. (2000). Amplitude Envelope correlation detects coupling among incoherent brain signals.Neuroreport1115091514. 10.1097/00001756-200005150-00029

  • 12

    BullmoreE.SpornsO. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems.Nat. Rev. Neurosci.10186198. 10.1038/nrn2575

  • 13

    BullmoreE.SpornsO. (2012). The economy of brain network organization.Nat. Rev. Neurosci.13336349. 10.1038/nrn3214

  • 14

    CampoP.PochC. (2012). Neocortical-hippocampal dynamics of working memory in healthy and diseased brain states based on functional connectivity.Front. Hum. Neurosci.6:36.

  • 15

    ChangL. C.JonesD. K.PierpaoliC. (2005). RESTORE: robust estimation of tensors by outlier rejection.Magn. Reson. Med.5310881095. 10.1002/mrm.20426

  • 16

    ChenY.HuangX. (2015). Modulation of alpha and beta oscillations during an n-back task with varying temporal memory load.Front. Psychol.6:2031.

  • 17

    ColcloughG.WoolrichM.TewarieP.BrookesM.QuinnA.SmithS. (2016). How reliable are MEG resting-state connectivity metrics?NeuroImage138284293. 10.1016/j.neuroimage.2016.05.070

  • 18

    CrossleyN.MechelliA.ScottJ.CarlettiF.FoxP.McguireP.et al (2014). The hubs of the human connectome are generally implicated in the anatomy of brain disorders.Brain13723822395. 10.1093/brain/awu132

  • 19

    CsuklyG.SirályE.FodorZ.HorváthA.SalaczP.HidasiZ.et al (2016). The differentiation of amnestic type MCI from the non-amnestic types by structural MRI.Front. Aging Neurosci.8:52.

  • 20

    D’AmelioM.RossiniP. (2012). Brain excitability and connectivity of neuronal assemblies in Alzheimer’s disease: from animal models to human findings.Progr. Neurobiol.994260. 10.1016/j.pneurobio.2012.07.001

  • 21

    DasS.PuthankattilS. D. (2020). Complex network analysis of MCI-AD EEG signals under cognitive and resting state.Brain Res.1735:146743. 10.1016/j.brainres.2020.146743

  • 22

    DauwelsJ.VialatteF.CichockiA. (2010a). Diagnosis of Alzheimer’s disease from EEG signals: where are we standing?Curr. Alzheimer Res.7487505. 10.2174/156720510792231720

  • 23

    DauwelsJ.VialatteF.MushaT.CichockiA. (2010b). A comparative study of synchrony measures for the early diagnosis of Alzheimer’s disease based on EEG.NeuroImage49668693. 10.1016/j.neuroimage.2009.06.056

  • 24

    de HaanW.MottK.Van StraatenE. C.ScheltensP.StamC. J. (2012). Activity dependent degeneration explains hub vulnerability in Alzheimer’s disease.PLoS Comput. Biol.8:e1002582. 10.1371/journal.pcbi.1002582

  • 25

    de HaanW.PijnenburgY. A.StrijersR. L.Van Der MadeY.Van Der FlierW. M.ScheltensP.et al (2009). Functional neural network analysis in frontotemporal dementia and Alzheimer’s disease using EEG and graph theory.BMC Neurosci.10:101. 10.1186/1471-2202-10-101

  • 26

    de WaalH.StamC. J.LansbergenM. M.WieggersR. L.KamphuisP. J.ScheltensP.et al (2014). The effect of souvenaid on functional brain network organisation in patients with mild Alzheimer’s disease: a randomised controlled study.PLoS One9:e86558. 10.1371/journal.pone.0086558

  • 27

    DeiberM. P.MissonnierP.BertrandO.GoldG.Fazio-CostaL.IbanezV.et al (2007). Distinction between perceptual and attentional processing in working memory tasks: a study of phase-locked and induced oscillatory brain dynamics.J. Cogn. Neurosci.19158172. 10.1162/jocn.2007.19.1.158

  • 28

    DelbeuckX.Van Der LindenM.ColletteF. (2003). Alzheimer’s disease as a disconnection syndrome?Neuropsychol. Rev.137992.

  • 29

    DelormeA.MakeigS. (2004). EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis.J. Neurosci. Methods134921. 10.1016/j.jneumeth.2003.10.009

  • 30

    DonnerT. H.SiegelM. (2011). A framework for local cortical oscillation patterns.Trends Cogn. Sci.15191199. 10.1016/j.tics.2011.03.007

  • 31

    DouwL.Van DellenE.GouwA.GriffaA.De HaanW.Van Den HeuvelM.et al (2019). The road ahead in clinical network neuroscience.Network Neurosci.3969993. 10.1162/netn_a_00103

  • 32

    EngelA. K.FriesP. (2010). Beta-band oscillations–signalling the status quo?Curr. Opin. Neurobiol.20156165. 10.1016/j.conb.2010.02.015

  • 33

    EngelsM. M.StamC. J.Van Der FlierW. M.ScheltensP.De WaalH.Van StraatenE. C. (2015). Declining functional connectivity and changing hub locations in Alzheimer’s disease: an EEG study.BMC Neurol15:145. 10.1186/s12883-015-0400-7

  • 34

    FellgiebelA.WilleP.MüllerM. J.WintererG.ScheurichA.VucurevicG.et al (2004). Ultrastructural hippocampal and white matter alterations in mild cognitive impairment: a diffusion tensor imaging study.Dement Geriatr Cogn. Dis.18101108. 10.1159/000077817

  • 35

    FergusonC. J. (2009). An effect size primer: a guide for clinicians and researchers.Professional. Psychol. Res. Practice40532538. 10.1037/a0015808

  • 36

    FodorZ.SiralyE.HorvathA.SalaczP.HidasiZ.CsibriE.et al (2018). Decreased event-related beta synchronization during memory maintenance marks early cognitive decline in mild cognitive impairment.J. Alzheimers Dis.63489502. 10.3233/jad-171079

  • 37

    FornitoA.ZaleskyA.BreakspearM. (2013). Graph analysis of the human connectome: promise, progress, and pitfalls.NeuroImage80426444. 10.1016/j.neuroimage.2013.04.087

  • 38

    FraschiniM.DemuruM.CrobeA.MarrosuF.StamC.HillebrandA. (2016). The effect of epoch length on estimated EEG functional connectivity and brain network organisation.J. Neural. Eng.13:036015. 10.1088/1741-2560/13/3/036015

  • 39

    FriesP. (2015). Rhythms for cognition: communication through coherence.Neuron88220235. 10.1016/j.neuron.2015.09.034

  • 40

    GillisM. M.QuinnK. M.PhillipsP. A.HampsteadB. M. (2013). Impaired retention is responsible for temporal order memory deficits in mild cognitive impairment.Acta Psychol.1438895. 10.1016/j.actpsy.2013.03.001

  • 41

    GriesmayrB.GruberW. R.KlimeschW.SausengP. (2010). Human frontal midline theta and its synchronization to gamma during a verbal delayed match to sample task.Neurobiol. Learn. Mem.93208215. 10.1016/j.nlm.2009.09.013

  • 42

    GyebnárG.SzabóÁSirályE.FodorZ.SákovicsA.SalaczP.et al (2018). What can DTI tell about early cognitive impairment?Differentiation between MCI subtypes and healthy controls by diffusion tensor imaging.Psychiatry Res. Neuroimaging2724657. 10.1016/j.pscychresns.2017.10.007

  • 43

    HallettM.De HaanW.DecoG.DenglerR.Di IorioR.GalleaC.et al (2020). Human brain connectivity: clinical applications for clinical neurophysiology.Clin. Neurophysiol.13116211651. 10.1016/j.clinph.2020.03.031

  • 44

    HerrerasO. (2016). Local field potentials: myths and misunderstandings.Front. Neural. Circ.10:101.

  • 45

    HidasiZ.CziglerB.SalaczP.CsibriE.MolnárM. (2007). Changes of EEG spectra and coherence following performance in a cognitive task in Alzheimer’s disease.Int. J. Psychophysiol.6252260. 10.1016/j.ijpsycho.2007.05.002

  • 46

    HippJ.HawellekD.CorbettaM.SiegelM.EngelA. (2012). Large-scale cortical correlation structure of spontaneous oscillatory activity.Nat. Neurosci.15884890. 10.1038/nn.3101

  • 47

    HochbergY. (1988). A sharper Bonferroni procedure for multiple tests of significance.Biometrika75800802. 10.1093/biomet/75.4.800

  • 48

    HochbergY.BenjaminiY. (1990). More powerful procedures for multiple significance testing.Stat. Med.9811818. 10.1002/sim.4780090710

  • 49

    HoganM.SwanwickG.KaiserJ.RowanM.LawlorB. (2003). Memory-related EEG power and coherence reductions in mild Alzheimer’s disease.Int. J. Psychophysiol.49147163. 10.1016/s0167-8760(03)00118-1

  • 50

    HorvathA.SzucsA.CsuklyG.SakovicsA.StefanicsG.KamondiA. (2018). EEG and ERP biomarkers of Alzheimer’s disease: a critical review.Front. Biosci.23:183220. 10.2741/4587

  • 51

    HouF.LiuC.YuZ.XuX.ZhangJ.PengC.et al (2018). Age-related alterations in electroencephalography connectivity and network topology during n-back working memory task.Front. Hum. Neurosci.12:484.

  • 52

    JensenO.GelfandJ.KouniosJ.LismanJ. E. (2002). Oscillations in the alpha band (9-12 Hz) increase with memory load during retention in a short-term memory task.Cerebral. Cortex12877882. 10.1093/cercor/12.8.877

  • 53

    JensenO.TescheC. D. (2002). Frontal theta activity in humans increases with memory load in a working memory task.Eur. J. Neurosci.1513951399. 10.1046/j.1460-9568.2002.01975.x

  • 54

    JezzardP.BarnettA. S.PierpaoliC. (1998). Characterization of and correction for eddy current artifacts in echo planar diffusion imaging.Magn. Reson. Med.39801812. 10.1002/mrm.1910390518

  • 55

    JiangZ. (2005). Abnormal cortical functional connections in Alzheimer’s disease: analysis of inter- and intra-hemispheric EEG coherence.J. Zhejiang Univ. Sci. B6259264. 10.1631/jzus.2005.b0259

  • 56

    JiangZ.ZhengL. (2006). Inter- and intra-hemispheric EEG coherence in patients with mild cognitive impairment at rest and during working memory task.J. Zhejiang Univ. Sci. B7357364. 10.1631/jzus.2006.b0357

  • 57

    KantarciK.JackC. R.Jr.XuY. C.CampeauN. G.O’brienP. C.SmithG. E.et al (2001). Mild cognitive impairment and Alzheimer disease: regional diffusivity of water.Radiology219101107. 10.1148/radiology.219.1.r01ap14101

  • 58

    KantarciK.PetersenR. C.BoeveB. F.KnopmanD. S.WeigandS. D.O’brienP. C.et al (2005). DWI predicts future progression to Alzheimer disease in amnestic mild cognitive impairment.Neurology64902904. 10.1212/01.wnl.0000153076.46126.e9

  • 59

    KardosZ.TóthB.BohaR.FileB.MolnárM. (2014). Age-related changes of frontal-midline theta is predictive of efficient memory maintenance.Neuroscience273152162. 10.1016/j.neuroscience.2014.04.071

  • 60

    KilavikB.ZaepffelM.BrovelliA.MackayW.RiehleA. (2013). The Ups and Downs of β oscillations in sensorimotor cortex.Exp. Neurol.2451526. 10.1016/j.expneurol.2012.09.014

  • 61

    KlimeschW. (1999). EEG alpha and theta oscillations reflect cognitive and memory performance: a review and analysis.Brain Res. Brain Res. Rev.29169195.

  • 62

    KlimeschW. (2012). Alpha-band oscillations, attention, and controlled access to stored information.Trends Cogn. Sci.16606617. 10.1016/j.tics.2012.10.007

  • 63

    KoelewijnL.BompasA.TalesA.BrookesM. J.MuthukumaraswamyS. D.BayerA.et al (2017). Alzheimer’s disease disrupts alpha and beta-band resting-state oscillatory network connectivity.Clin. Neurophysiol.12823472357. 10.1016/j.clinph.2017.04.018

  • 64

    KopellN.ErmentroutG.WhittingtonM.TraubR. (2000). Gamma rhythms and beta rhythms have different synchronization properties.Proc. Natl. Acad. Sci. U.S.A.9718671872. 10.1073/pnas.97.4.1867

  • 65

    LaiM.DemuruM.HillebrandA.FraschiniM. (2018). A comparison between scalp- and source-reconstructed EEG networks.Sci. Rep.8:12269.

  • 66

    LeemansA.JeurissenB.SijbersJ.JonesD. K. (2009). “ExploreDTI: A graphical toolbox for processing, analyzing, and visualizing diffusion MR data,” in Proceeding of the 17th Annual Meeting of Intl Soc Mag Reson Med17.

  • 67

    LeemansA.JonesD. K. (2009). The B-matrix must be rotated when correcting for subject motion in DTI data.Magn. Reson. Med.6113361349. 10.1002/mrm.21890

  • 68

    LopezM. E.EngelsM. M. A.Van StraatenE. C. W.BajoR.DelgadoM. L.ScheltensP.et al (2017). MEG beamformer-based reconstructions of functional networks in mild cognitive impairment.Front. Aging Neurosci.9:107.

  • 69

    MárquezF.YassaM. A. (2019). Neuroimaging biomarkers for Alzheimer’s disease.Mol. Neurodegener14:21.

  • 70

    MarzettiL.BastiA.ChellaF.D’andreaA.SyrjäläJ.PizzellaV. (2019). Brain functional connectivity through phase coupling of neuronal oscillations: a perspective from magnetoencephalography.Front. Neurosci.13:964.

  • 71

    MathuranathP. S.NestorP. J.BerriosG. E.RakowiczW.HodgesJ. R. (2000). A brief cognitive test battery to differentiate Alzheimer’s disease and frontotemporal dementia.Neurology5516131620. 10.1212/01.wnl.0000434309.85312.19

  • 72

    MazaheriA.SegaertK.OlichneyJ.YangJ.-C.NiuY.-Q.ShapiroK.et al (2018). EEG oscillations during word processing predict MCI conversion to Alzheimer’s disease.NeuroImage Clin.17188197. 10.1016/j.nicl.2017.10.009

  • 73

    Metting van RijnA. C.PeperA.GrimbergenC. A. (1990). High-quality recording of bioelectric events. part 1. interference reduction, theory and practice.Med. Biol. Eng. Comput.28389397. 10.1007/bf02441961

  • 74

    MiragliaF.VecchioF.BramantiP.RossiniP. M. (2016). EEG characteristics in ”eyes-open” versus ”eyes-closed” conditions: Small-world network architecture in healthy aging and age-related brain degeneration.Clin. Neurophysiol.12712611268. 10.1016/j.clinph.2015.07.040

  • 75

    MiragliaF.VecchioF.RossiniP. (2017). Searching for signs of aging and dementia in EEG through network analysis.Behav. Brain Res.317292300. 10.1016/j.bbr.2016.09.057

  • 76

    MissonnierP.DeiberM. P.GoldG.HerrmannF. R.MilletP.MichonA.et al (2007). Working memory load-related electroencephalographic parameters can differentiate progressive from stable mild cognitive impairment.Neuroscience150346356. 10.1016/j.neuroscience.2007.09.009

  • 77

    MoodleyK.MinatiL.ContarinoV.PrioniS.WoodR.CooperR.et al (2015). Diagnostic differentiation of mild cognitive impairment due to alzheimer’s disease using a hippocampus-dependent test of spatial memory.Hippocampus25939951. 10.1002/hipo.22417

  • 78

    MorettiD. V.FrisoniG. B.FracassiC.PievaniM.GeroldiC.BinettiG.et al (2011). MCI patients’ EEGs show group differences between those who progress and those who do not progress to AD.Neurobiol. Aging32563571. 10.1016/j.neurobiolaging.2009.04.003

  • 79

    MusaeusC. S.EngedalK.HøghP.JelicV.MørupM.NaikM.et al (2019). Oscillatory connectivity as a diagnostic marker of dementia due to Alzheimer’s disease.Clin. Neurophysiol.13018891899. 10.1016/j.clinph.2019.07.016

  • 80

    NúñezP.PozaJ.GómezC.Rodríguez-GonzálezV.HillebrandA.Tola-ArribasM.et al (2019). Characterizing the fluctuations of dynamic resting-state electrophysiological functional connectivity: reduced neuronal coupling variability in mild cognitive impairment and dementia due to Alzheimer’s disease.J. Neural. Eng.16:056030. 10.1088/1741-2552/ab234b

  • 81

    OntonJ.DelormeA.MakeigS. (2005). Frontal midline EEG dynamics during working memory.Neuroimage27341356. 10.1016/j.neuroimage.2005.04.014

  • 82

    PalvaJ. M.MontoS.KulashekharS.PalvaS. (2010). Neuronal synchrony reveals working memory networks and predicts individual memory capacity.Proc. Natl. Acad. Sci. U.S.A.10775807585. 10.1073/pnas.0913113107

  • 83

    PalvaS.KulashekharS.HamalainenM.PalvaJ. M. (2011). Localization of cortical phase and amplitude dynamics during visual working memory encoding and retention.J. Neurosci.3150135025. 10.1523/jneurosci.5592-10.2011

  • 84

    PalvaS.PalvaJ. (2011). Functional roles of alpha-band phase synchronization in local and large-scale cortical networks.Front. Psychol.2:204.

  • 85

    ParkH.FristonK. (2013). Structural and functional brain networks: from connections to cognition.Science342:1238411. 10.1126/science.1238411

  • 86

    ParkH.KangE.KangH.KimJ. S.JensenO.ChungC. K.et al (2011). Cross-frequency power correlations reveal the right superior temporal gyrus as a hub region during working memory maintenance.Brain Connect1460472. 10.1089/brain.2011.0046

  • 87

    ParraM.AbrahamsS.LogieR.MéndezL.LoperaF.Della SalaS. (2010). Visual short-term memory binding deficits in familial Alzheimer’s disease.Brain J. Neurol.13327022713. 10.1093/brain/awq148

  • 88

    PerazaL. R.CromartyR.KobelevaX.FirbankM. J.KillenA.GraziadioS.et al (2018). Electroencephalographic derived network differences in Lewy body dementia compared to Alzheimer’s disease patients.Sci. Rep.8:4637.

  • 89

    PetersenR. C. (2004). Mild cognitive impairment as a diagnostic entity.J. Intern. Med.256183194. 10.1111/j.1365-2796.2004.01388.x

  • 90

    PfurtschellerG.Lopes da SilvaF. H. (1999). Event-related EEG/MEG synchronization and desynchronization: basic principles.Clin. Neurophysiol.11018421857. 10.1016/s1388-2457(99)00141-8

  • 91

    PierpaoliC.BasserP. J. (1996). Toward a quantitative assessment of diffusion anisotropy.Magn. Reson. Med.36893906. 10.1002/mrm.1910360612

  • 92

    PievaniM.De HaanW.WuT.SeeleyW. W.FrisoniG. B. (2011). Functional network disruption in the degenerative dementias.Lancet Neurol.10829843. 10.1016/s1474-4422(11)70158-2

  • 93

    PijnenburgY. A.v d MadeY.van Cappellen van WalsumA. M.KnolD. L.ScheltensP.StamC. J. (2004). EEG synchronization likelihood in mild cognitive impairment and Alzheimer’s disease during a working memory task.Clin. Neurophysiol.11513321339. 10.1016/j.clinph.2003.12.029

  • 94

    RanganathC. (2006). Working memory for visual objects: complementary roles of inferior temporal, medial temporal, and prefrontal cortex.Neuroscience139277289. 10.1016/j.neuroscience.2005.06.092

  • 95

    RayK. M.WangH.ChuY.ChenY. F.BertA.HassoA. N.et al (2006). Mild cognitive impairment: apparent diffusion coefficient in regional gray matter and white matter structures.Radiology241197205. 10.1148/radiol.2411051051

  • 96

    RogalskiE. J.MurphyC. M.Detoledo-MorrellL.ShahR. C.MoseleyM. E.BammerR.et al (2009). Changes in parahippocampal white matter integrity in amnestic mild cognitive impairment: a diffusion tensor imaging study.Behav. Neurol.215161. 10.1155/2009/408037

  • 97

    RossiniP.Di IorioR.GranataG.MiragliaF.VecchioF. (2016). From mild cognitive impairment to Alzheimer’s disease: a new perspective in the ”land” of human brain reactivity and connectivity.J. Alzheimer’s Dis.5313891393. 10.3233/jad-160482

  • 98

    RubinovM.SpornsO. (2010). Complex network measures of brain connectivity: uses and interpretations.NeuroImage5210591069. 10.1016/j.neuroimage.2009.10.003

  • 99

    San-MartinR.JohnsE.Quispe MamaniG.TavaresG.PhillipsN. A.FragaF. J. (2021). A method for diagnosis support of mild cognitive impairment through EEG rhythms source location during working memory tasks.Biomed. Signal Proc. Control66:102499. 10.1016/j.bspc.2021.102499

  • 100

    SanoM.RamanR.EmondJ.ThomasR. G.PetersenR.SchneiderL. S.et al (2011). Adding delayed recall to the Alzheimer disease assessment scale is useful in studies of mild cognitive impairment but not Alzheimer disease.Alzheimer Dis. Assoc. Disord25122127. 10.1097/wad.0b013e3181f883b7

  • 101

    SatoJ.MossadS.WongS.HuntB.DunkleyB.SmithM.et al (2018). Alpha keeps it together: alpha oscillatory synchrony underlies working memory maintenance in young children.Dev. Cogn. Neurosci.34114123. 10.1016/j.dcn.2018.09.001

  • 102

    SausengP.GriesmayrB.FreunbergerR.KlimeschW. (2010). Control mechanisms in working memory: a possible function of EEG theta oscillations.Neurosci. Biobehav. Rev.3410151022. 10.1016/j.neubiorev.2009.12.006

  • 103

    SharmaN.KolekarM. H.JhaK.KumarY. (2019). EEG and cognitive biomarkers based mild cognitive impairment diagnosis.IRBM40113121. 10.1016/j.irbm.2018.11.007

  • 104

    SiralyE.SzaboA.SzitaB.KovacsV.FodorZ.MarosiC.et al (2015). Monitoring the early signs of cognitive decline in elderly by computer games: an MRI study.PLoS One10:e0117918. 10.1371/journal.pone.0117918

  • 105

    SirályE.SzitaB.KovácsV.CsibriÉHidasiZ.SalaczP.et al (2013). Differentiation between mild cognitive impairment and healthy elderly population using neuropsychological tests.Neuropsychopharmacol. Hung.15139146.

  • 106

    SpielbergerC.GorsuchR.LusheneR. (1970). Manual for the State-Trait Anxiety Inventory.Santa Clara, CA: Palo Alto.

  • 107

    SpitzerB.HaegensS. (2017). Beyond the status quo: a role for beta oscillations in endogenous content (Re)activation.eNeuro4:ENEURO.170ENEURO.117.

  • 108

    SpornsO. (2013). Network attributes for segregation and integration in the human brain.Curr. Opin. Neurobiol.23162171. 10.1016/j.conb.2012.11.015

  • 109

    StamC.Van Der MadeY.PijnenburgY.ScheltensP. (2003). EEG Synchronization in mild cognitive impairment and Alzheimer’s disease.Acta Neurol. Scand.1089096. 10.1034/j.1600-0404.2003.02067.x

  • 110

    StamC. J. (2014). Modern network science of neurological disorders.Nat. Rev. Neurosci.15683695. 10.1038/nrn3801

  • 111

    StamC. J.De HaanW.DaffertshoferA.JonesB. F.ManshandenI.Van Cappellen Van WalsumA. M.et al (2009). Graph theoretical analysis of magnetoencephalographic functional connectivity in Alzheimer’s disease.Brain132213224. 10.1093/brain/awn262

  • 112

    StamC. J.TewarieP.Van DellenE.Van StraatenE. C.HillebrandA.Van MieghemP. (2014). The trees and the forest: characterization of complex brain networks with minimum spanning trees.Int. J. Psychophysiol.92129138. 10.1016/j.ijpsycho.2014.04.001

  • 113

    StamC. J.van StraatenE. C. (2012). The organization of physiological brain networks.Clin. Neurophysiol.12310671087. 10.1016/j.clinph.2012.01.011

  • 114

    StebbinsG. T.MurphyC. M. (2009). Diffusion tensor imaging in Alzheimer’s disease and mild cognitive impairment.Behav. Neurol.21:915041.

  • 115

    StraussE.ShermanE.SpreenO. (2006). A Compendium of Neuropsychological Tests, 3rd Edn. Oxford: Oxford University Press, 168188.

  • 116

    TeipelS.PogarellO.MeindlT.DietrichO.SydykovaD.HunklingerU.et al (2009). Regional networks underlying interhemispheric connectivity: an EEG and DTI study in healthy ageing and amnestic mild cognitive impairment.Hum. Brain Mapp.3020982119. 10.1002/hbm.20652

  • 117

    TewarieP.Van DellenE.HillebrandA.StamC. J. (2015). The minimum spanning tree: an unbiased method for brain network analysis.Neuroimage104177188. 10.1016/j.neuroimage.2014.10.015

  • 118

    TierneyM. C.SzalaiJ. P.SnowW. G.FisherR. H.NoresA.NadonG.et al (1996). Prediction of probable Alzheimer’s disease in memory-impaired patients: a prospective longitudinal study.Neurology46661665. 10.1212/wnl.46.3.661

  • 119

    TijmsB. M.WinkA. M.De HaanW.Van Der FlierW. M.StamC. J.ScheltensP.et al (2013). Alzheimer’s disease: connecting findings from graph theoretical studies of brain networks.Neurobiol. Aging3420232036.

  • 120

    TombaughT. N. (2004). Trail making Test A and B: normative data stratified by age and education.Arch. Clin. Neuropsychol.19203214. 10.1016/s0887-6177(03)00039-8

  • 121

    TothB.FileB.BohaR.KardosZ.HidasiZ.GaalZ. A.et al (2014). EEG network connectivity changes in mild cognitive impairmentPreliminary results.Int J Psychophysiol9217. 10.1016/j.ijpsycho.2014.02.001

  • 122

    TuladharA. M.Ter HuurneN.SchoffelenJ. M.MarisE.OostenveldR.JensenO. (2007). Parieto-occipital sources account for the increase in alpha activity with working memory load.Hum. Brain Mapp.28785792. 10.1002/hbm.20306

  • 123

    Van DellenE.De WaalH.Van Der FlierW. M.LemstraA. W.SlooterA. J.SmitsL. L.et al (2015). Loss of EEG network efficiency is related to cognitive impairment in dementia with lewy bodies.Mov. Disord.3017851793. 10.1002/mds.26309

  • 124

    van DellenE.DouwL.HillebrandA.De Witt HamerP. C.BaayenJ. C.HeimansJ. J.et al (2014). Epilepsy surgery outcome and functional network alterations in longitudinal MEG: a minimum spanning tree analysis.Neuroimage86354363. 10.1016/j.neuroimage.2013.10.010

  • 125

    van den HeuvelM. P.SpornsO. (2013). Network hubs in the human brain.Trends Cogn. Sci.17683696. 10.1016/j.tics.2013.09.012

  • 126

    van den HeuvelM. P.StamC. J.KahnR. S.Hulshoff PolH. E. (2009). Efficiency of functional brain networks and intellectual performance.J. Neurosci.2976197624. 10.1523/jneurosci.1443-09.2009

  • 127

    van der HieleK.VeinA. A.KramerC. G.ReijntjesR. H.Van BuchemM. A.WestendorpR. G.et al (2007a). Memory activation enhances EEG abnormality in mild cognitive impairment.Neurobiol. Aging288590. 10.1016/j.neurobiolaging.2005.11.006

  • 128

    van der HieleK.VeinA. A.ReijntjesR. H.WestendorpR. G.BollenE. L.Van BuchemM. A.et al (2007b). EEG correlates in the spectrum of cognitive decline.Clin. Neurophysiol.11819311939. 10.1016/j.clinph.2007.05.070

  • 129

    van LutterveldR.Van DellenE.PalP.YangH.StamC. J.BrewerJ. (2017). Meditation is associated with increased brain network integration.Neuroimage1581825. 10.1016/j.neuroimage.2017.06.071

  • 130

    van WijkB. C.StamC. J.DaffertshoferA. (2010). Comparing brain networks of different size and connectivity density using graph theory.PLoS One5:e13701. 10.1371/journal.pone.0013701

  • 131

    VarelaF.LachauxJ.RodriguezE.MartinerieJ. (2001). The brainweb: phase synchronization and large-scale integration.Nat. Rev. Neurosci.2229239. 10.1038/35067550

  • 132

    VecchioF.MiragliaF.CurcioG.AltavillaR.ScrasciaF.GiambattistelliF.et al (2015). Cortical brain connectivity evaluated by graph theory in dementia: a correlation study between functional and structural data.J. Alzheimer’s Dis.45745756. 10.3233/jad-142484

  • 133

    VecchioF.MiragliaF.Maria RossiniP. (2017). Connectome: graph theory application in functional brain network architecture.Clin. Neurophysiol. Pract.2206213. 10.1016/j.cnp.2017.09.003

  • 134

    WangB.MiaoL.NiuY.CaoR.LiD.YanP.et al (2018). Abnormal functional brain networks in mild cognitive impairment and Alzheimer’s disease: a minimum spanning tree analysis.J. Alzheimer’s Dis.6510931107. 10.3233/jad-180603

  • 135

    WangR.WangJ.YuH.WeiX.YangC.DengB. (2014). Decreased coherence and functional connectivity of electroencephalograph in Alzheimer’s disease.Chaos24:033136. 10.1063/1.4896095

  • 136

    WeiL.LiY.YangX.XueQ.WangY. (2015). Altered characteristic of brain networks in mild cognitive impairment during a selective attention task: an EEG study.Int. J. Psychophysiol.98816. 10.1016/j.ijpsycho.2015.05.015

  • 137

    WestmanE.AguilarC.MuehlboeckJ. S.SimmonsA. (2013). Regional magnetic resonance imaging measures for multivariate analysis in Alzheimer’s disease and mild cognitive impairment.Brain Topogr26923. 10.1007/s10548-012-0246-x

  • 138

    WiandaE.RossB. (2019). The roles of alpha oscillation in working memory retention.Brain Behav.9:1263.

  • 139

    WinklerI.BrandlS.HornF.WaldburgerE.AllefeldC.TangermannM. (2014). Robust artifactual independent component classification for bci practitioners.J. Neural. Eng.11:035013. 10.1088/1741-2560/11/3/035013

  • 140

    WinklerI.HaufeS.TangermannM. (2011). Automatic classification of artifactual ICA-components for artifact removal in EEG signals.Behav. Brain Functions7:30. 10.1186/1744-9081-7-30

  • 141

    YesavageJ. A. (1988). Geriatric depression scale.Psychopharmacol. Bull24709711.

  • 142

    YuM.EngelsM.HillebrandA.Van StraatenE.GouwA.TeunissenC.et al (2017). Selective impairment of hippocampus and posterior hub areas in alzheimer’s disease: an MEG-based multiplex network study.Brain14014661485. 10.1093/brain/awx050

  • 143

    YuM.GouwA. A.HillebrandA.TijmsB. M.StamC. J.Van StraatenE. C.et al (2016). Different functional connectivity and network topology in behavioral variant of frontotemporal dementia and Alzheimer’s disease: an EEG study.Neurobiol. Aging42150162. 10.1016/j.neurobiolaging.2016.03.018

  • 144

    ZhangB.XuY.ZhuB.KantarciK. (2014). The role of diffusion tensor imaging in detecting microstructural changes in prodromal Alzheimer’s disease.CNS Neurosci. Ther.2039. 10.1111/cns.12166

Summary

Keywords

mild cognitive impairment (MCI), electroencephalography (EEG), working memory (WM), functional connectivity, functional networks, minimum spanning tree (MST)

Citation

Fodor Z, Horváth A, Hidasi Z, Gouw AA, Stam CJ and Csukly G (2021) EEG Alpha and Beta Band Functional Connectivity and Network Structure Mark Hub Overload in Mild Cognitive Impairment During Memory Maintenance. Front. Aging Neurosci. 13:680200. doi: 10.3389/fnagi.2021.680200

Received

13 March 2021

Accepted

20 September 2021

Published

07 October 2021

Volume

13 - 2021

Edited by

Panagiotis D. Bamidis, Aristotle University of Thessaloniki, Greece

Reviewed by

Carmen Jiménez-Mesa, University of Granada, Spain; Márk Molnár, Hungarian Academy of Sciences (MTA), Hungary

Updates

Copyright

*Correspondence: Gábor Csukly, ;

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics