Next Article in Journal
Comparison Between Performance of Fluorite Flotation Under Different Depressants Reagents in Two Pieces of Laboratory Equipment
Previous Article in Journal
Self-Assembled Three-Dimensional Microporous rGO/PNT/Fe3O4 Hydrogel Sorbent for Magnetic Preconcentration of Multi-Residue Insecticides
Previous Article in Special Issue
Cortical Reorganization after Rehabilitation in a Patient with Conduction Aphasia Using High-Density EEG
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An eLORETA Longitudinal Analysis of Resting State EEG Rhythms in Alzheimer’s Disease

DICEAM Department, Mediterranea University of Reggio Calabria, Via Graziella Feo di Vito, 89060 Reggio Calabria, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(16), 5666; https://doi.org/10.3390/app10165666
Submission received: 7 July 2020 / Revised: 3 August 2020 / Accepted: 13 August 2020 / Published: 15 August 2020
(This article belongs to the Special Issue Advances in Biomedical Signal Processing)

Abstract

:
Alzheimer’s disease (AD) is a degenerative brain disorder which is the most common cause of dementia. As there is no cure for AD, an early diagnosis is essential to slow down the progression of the disease with a proper pharmacological treatment. Electroencephalography (EEG) represents a valid tool for studying AD. EEG signals of AD patients are characterized by a “slowing”, meaning the power increases in low frequencies (delta and theta) and decreases in higher frequency (alpha and beta), compared to normal elderly. The purpose of our study is the computation of the power current density in eight patients, who were diagnosed with MCI at time T0 and mild AD at time T1 (four months later), starting from the brain active source reconstruction. The novelty is that we employed the eLORETA algorithm, unlike the previous studies which used the old version of the algorithm named LORETA. It is also the first longitudinal study which considers such a short time period to explore the evolution of the disease. Five patients out of eight showed an increasing power in delta and theta bands. Seven patients exhibited a lower activation in alpha 1 and beta 2 bands. Finally, six patients revealed a decreased power in alpha 2 and beta 1 bands. These findings are consistent with those reported in literature. On the other hand, the discrepancy of some outcome could be related to a not yet severe stage of the disease. In our opinion, this study could represent a good starting point for more detailed future investigation.

1. Introduction

Alzheimer’s disease (AD) is a brain disorder causing a progressive decline in cognitive functions. The diagnosis of the disease is based, initially, on the evaluation of neuropsychological and behavioral symptoms related to disorders of memory, language, and visual-spatial perception [1].
There are four stages which characterize the evolution of the disease: In the first stage, called Mild Cognitive Impairment (MCI), the subject shows short-term memory loss and difficulty in remembering even words of daily use. These symptoms, however, do not substantially affect everyday life. This condition is not always irreversible because not all MCI subjects convert to AD—it happens for 50–60% of them. During the next two stages (Mild and Moderate AD) the patient shows difficulty in recognizing relatives and friends, and suffers from spatial and temporal disorientation. Moreover, the first behavioral changes appear, such as anxiety, insomnia, and personality changes. During the last stage—Severe AD—the patient no longer recognizes anyone; does not understand words; and has difficulty in dressing, eating, and swallowing. The subject is completely dependent on caregivers. Death occurs because of the further worsening of these symptoms [2,3,4]. To date, there is no cure for Alzheimer’s disease, but several medications are used for delaying the onset of the symptoms and improving the patient’s quality of life. Therefore, early diagnosis of AD is very important because the medications are effective only if taken during the first stage. However, the diagnosis is hard because the first symptoms are often mistaken for normal aging [5].
During the last decades, neuroimaging methods—such as magnetic resonance imaging (MRI) and positron emission tomography (PET)—have played an even more important role in diagnosing AD [6,7,8]. Nevertheless, most of them are expensive and rather invasive, thus the need of a tool which overcomes these drawbacks. Electroencephalography (EEG) represents a valid technique to detect any pathological alterations in the brain [9,10]. EEG is a low cost and noninvasive instrument, which records the brain electrical activity by means of electrodes distributed on the scalp. In particular, the signal is produced by the sum of postsynaptic potentials generated simultaneously by population of cortical pyramidal neurons. The EEG signal derives from the contribution of all active sources because of the conductivity of the scalp, the skull and the brain. Therefore, the choice of a realistic head model is essential for a correct EEG analysis [11]. EEG is characterized by an excellent temporal resolution, on the order of milliseconds; conversely, the spatial resolution is poor.
EEG has been successfully employed for investigating AD [12,13,14,15,16]. Several studies have proven that the EEG frequency bands are affected by AD in a specific way. In particular, the spectral analysis shows a “slowing” of EEG signals in MCI and AD patients: the power increases in low frequencies (delta and theta) and decreases in higher frequency (alpha and beta), compared to normal elderly (Nold) [17,18,19,20,21,22]. Moreover, the EEG signals of AD patients are characterized by less complexity and a reduction of synchronization between different areas of the brain [5,23,24].
The current flows generated by neurons can be modeled with dipoles. The localization of cortical sources which produce the EEG rhythms requires the solution of the so-called inverse problem. It is an ill-posed problem which has infinite possible solutions, because the number of the known parameters (electrodes) is lower than the number of the unknown parameters (current density of the sources). Therefore, the solutions are regularized using mathematical and neuroanatomical constraints, in order to choose the most likely solution. One of the most employed algorithm for solving the EEG inverse problem is Low-Resolution Electromagnetic Tomography (LORETA) [25]. It belongs to the high temporal resolution neuroimaging method and provides a 3D reconstruction of the cortical electrical activity. Over the years, LORETA has been widely employed for studying the spectral power density of EEG signals in AD. Babiloni carried out several studies about this topic. In [26], the results showed that the power in alpha 1 declined in AD in the central, parietal, temporal, and limbic lobes as compared to Nold. Alpha 1 power in the occipital lobe was lower when the disease was at a more severe stage. There was also an “anteriorization” of alpha 1 rhythm, as the power decreased much more in the parieto-occipital area than in the frontal one. As for the theta band, there was no significant difference between mild AD and Nold. This study was extended in [27], as it included also MCI. The findings revealed that, in comparison to Nold and mild AD groups, MCI group had intermediate power in delta (occipital area) and alpha 1 (parietal, occipital, temporal, and limbic areas) bands and greater power in alpha 2 (occipital area) band. In [28], the correlation between the disease and the hippocampal volume was explored. For the alpha 1 band, the activity was maximum in MCI with greater hippocampal volume, intermediate in MCI with smaller hippocampal volume, and minimum in AD. The study conducted in [29] showed that for MCI and AD subjects, posterior delta and alpha EEG rhythms do not decline when white matter vascular lesion increased. In [30], it was found a positive correlation between the gray matter volume and the alpha 1 power and a negative correlation with delta power in MCI and AD. Canuet and colleagues [31] reported a significant reduction of alpha 1 activity in the parieto-occipital regions for AD as compared to healthy subjects. In [32], the authors found a positive correlation between the severity of the disease and delta power and a negative correlation with alpha 1 power. The study in [33] showed a greater power in delta band in large areas, including the frontal lobe, in AD in comparison to controls. All the above-mentioned papers refer to cross-sectional studies. In the literature, there are only three longitudinal studies that applied LORETA. In [34], the EEG power of a group of MCI patients was evaluated during a follow-up of approximately 14 months. Twenty-four subjects converted to AD, whereas 45 subjects remained stable. The results showed that the EEG power for delta (temporal area), theta (temporal, parietal, and occipital areas) and alpha 1 (central, temporal, parietal, occipital, and limbic areas) bands was significantly higher in MCI who converted to AD, as compared to MCI stable. The remaining two works were carried out by Babiloni and colleagues. In [35], a group of MCI patients were examined about 1 year after the first EEG recording. The subjects were divided into three groups, according to their MMSE score variation, as “stable”, “decreased”, and “increased”. The only significant difference was detected for the alpha 1 band: the temporal, parietal, and occipital power density was greater in the stable group than in the decreased and increased groups, suggesting that preserved resting cortical activity at alpha band is related to a long-term stable cognitive function in MCI subjects. Resting state EEG of a mild AD group was investigated in [36]. The second recordings were carried out after 1 year. The alpha rhythm was characterized by a widespread power decrease as well as the temporal, parietal and occipital beta 1 band. It was also shown an increased power of delta sources in all brain areas except for the central region.
Motivated by the encouraging results of the above-mentioned works, at least at group level, we propose a preliminary study to evaluate whether resting state EEG sources can be markers for the prediction of the disease progression, from MCI to AD, in a single subject. In particular, the purpose of our work is to explore the power current density in eight patients, who were diagnosed with MCI at time T0 and mild AD at time T1 (four months later). Starting from the brain active source reconstruction, the eLORETA algorithm was used for the computation of the power spectral density for each EEG rhythm and the results were grouped according to the brain lobes. The use of the eLORETA algorithm represents a novelty as compared to the above-mentioned studies, which employ the old version named LORETA. It is also the first longitudinal study that considers such a short time period to examine the evolution of the disease. Most results are consistent with those reported in literature. On the other hand, in our opinion the discrepancy of some outcome could be related to a still mild stage of the disease at time T1.

2. Materials and Methods

2.1. Subjects

Eight amnestic MCI patients (four males and four females) were enrolled by a team of experts at the IRCCS Centro Neurolesi Bonino-Pulejo of Messina (Italy). The subjects underwent a series of examination at time T0 and at time T1 (four months later). The doctors agreed that patients converted from MCI (at T0) to AD (at T1) on the basis of the Montreal Cognitive Assessment (MoCA) scores. The study excluded the subjects suffering from epilepsy, neurological or psychiatric diseases related to cognitive impairment, traumatic brain injuries, systemic diseases. Subjects who took psychoactive drugs were also excluded. The study followed a clinical protocol approved by the the local Ethics Committee (Prot. E29/16). The patients or their caregivers were notified of the procedure and signed an informed consent form. The demographic data of the subjects are reported in Table 1.

2.2. EEG Data

The EEG signals were acquired in resting state with eyes closed but the subjects were awake. The data were recorded using a Brain-Quick System (Micromed) from nineteen channels, according to the international 10–20 system: Fp1, Fp2, F3, F4, C3, C4, P3, P4, O1, O2, F7, F8, T3, T4, T5, T6, Fz, Cz, and Pz, with linked earlobe (A1-A2) reference. The artifactual segments were detected and marked by the doctors during the EEG acquisition. Those segments were not considered for the analysis. The sampling rate was 256 Hz. It was applied a notch filter at 50 Hz. Data were transformed to a common average reference montage and divided into artifact-free non-overlapping epochs of 256 samples (1s).

2.3. LORETA

LORETA was the first linear distributed method for solving the EEG inverse problem which provided a valid solution from the neurophysiological point of view [25]. LORETA, in fact, solved the matter of the large source localization error of the minimum norm method, which misplaced deep source onto the outermost cortex [37]. LORETA employs a three-shell spherical head model (MNI305) registered to the Talairach human brain atlas [38,39]. The head model is made up of the scalp, skull, and brain compartment. The brain compartment, limited to the gray matter and hippocampus, is discretized in 2394 voxels at 7 mm spatial resolution. Each source is placed on a voxel and is characterized by a current density vector with unknown components. Therefore, the potential recorded on the scalp derives from the electrical activity at each point. LORETA seeks the smoothest of all possible infinite solutions by means of a discrete spatial Laplacian operator. Therefore, the LORETA images represent the electrical activity at each voxel as power of the computed current density. Because of the smoothing, the location of the maximal activity is maintained but the images are “blurred”. The average localization error is one voxel.
The expression of the general inverse problem is [40]
min J F W
with
F W = Φ KJ 2 + α J T WJ
where Φ R N E x 1 is the vector of the electric potential differences measured at N E electrodes with respect to a single common reference electrode; K R N E x ( 3 N V ) is the lead field matrix corresponding to N v voxels, J R ( 3 N V ) x 1 is the current density, and α > 0 is the Tikhonov regularization parameter [41].
The solution is
J ^ W = T W Φ
here T W is the pseudoinverse, given by
T W = W 1 K T ( KW 1 K T + α H ) +
The matrix H is the average reference operator. For LORETA, the matrix W represents the squared spatial discrete Laplacian operator.
After a few years, the sLORETA algorithm was introduced [42]. Unlike LORETA, the brain volume is made up of 6239 voxels at 5 mm spatial resolution. The electric potential lead field used for the inverse solution is computed with the boundary element method applied to the MNI152 template [43,44]. Moreover, sLORETA abandons the use of the laplacian operator and is based on the standardization of the current density. Differently from the Dale method [45], which takes into account only the variation due to measurement noise, in sLORETA the estimated current density considers the variation of the actual sources and, if any, the variation due to noisy measurements. The electric potential variance is given by
S Φ = KS J K T + S Φ N o i s e
where S Φ N o i s e is the variance of the noisy measurements and S J is the actual source variance (prior). The latter, according to the Bayesian formulation, is equal to the identity matrix. Therefore, the variance of the estimated current density can be expressed as
S J ^ = TS Φ T T = T ( K K T + α H ) T T = K T [ KK T + α H ] + K
Finally, sLORETA is given by the estimates of standardized current density power:
σ v T σ v = j ^ v T [ S J ^ ] v 1 j ^ v
where [ S J ^ ] v is the vth 3 × 3 diagonal matrix in [ S J ^ ] . As for the localization accuracy, sLORETA has zero localization error in ideal (no-noise) conditions and a lower localization error as compared with both the minimum norm and the Dale methods in the case of noisy measurements [37,45].
The latest algorithm of the LORETA family was eLORETA [46]. It introduces a weight matrix W which considers the deeper sources more adequately. Therefore, in the case of eLORETA, the matrix W in Equation (4), is defined as
W v = [ K v T ( KW 1 K T + α H ) + K v ] 1 / 2
The matrix W v R 3 × 3 is the vth diagonal subblock of W . The complete mathematical formulation of the three algorithms is reported in [25,40,42,46].
eLORETA outperforms sLORETA in terms of resolution and localization error, as reported in [47]. In particular, seventeen subjects underwent a visual stimulus and the corresponding ERP data were collected by 25 electrodes. The active source reconstruction was performed with both sLORETA and eLORETA methods. The comparison between the resulting images at four different time instants showed that eLORETA removes less significant sources and generates less blurred images as compared with sLORETA.
In our paper, the computation of the power current density was conducted by means of the eLORETA algorithm implemented in the LORETA-KEY software (v20200414). The power current density was evaluated by eLORETA for each patient, at time T0 and T1, for the following frequency sub-bands: delta (0.5–4 Hz), theta (4–8 Hz), alpha 1 (8–10.5 Hz), alpha 2 (10.5–13 Hz), beta 1 (13–18 Hz), and beta 2 (18–30 Hz). The comparison between the values at T0 and T1 was carried out by grouping the eLORETA solutions by frontal, occipital, parietal and temporal lobes.

3. Results

The methodology explained in Section 2.3 was applied to the EEG recordings and produced the following results. The EEG preprocessing and the grouping of power current density values according to the four regions of interest were executed with Matlab (R2018a). A statistical analysis was carried out for determining whether the differences between values at T0 and T1 were significant. For this purpose, the Wilcoxon rank sum test [48] was performed under the null hypothesis that, for each patient, the medians of the current density values at time T0 do not differ from those at time T1. The significance level was set at 5%, which means that there is a 5% risk of concluding that a difference exists when there is no actual difference. So, the difference between T0 and T1 is statistically significant when the p-value is less than 0.05. Table A1, Table A2 and Table A3 (see Appendix A) report the p-values resulting from the statistical analysis for each sub-band and each region of interest. F. L., O. L., P. L., and T. L. stand for Frontal Lobe, Occipital Lobe, Parietal Lobe, Temporal Lobe, respectively. The p-values statistically significant which are consistent with the expected results are highlighted in gray. Patient 03 shows a widespread increase of power density in delta and theta bands, and a decrease in alpha 1 (parietal region) and beta 2 (temporal lobe) bands. Patient 32 exhibits a greater activation of delta sources in all regions; an increasing power in theta band in the occipital, parietal, and temporal lobes; and a lower activity in alpha 1, alpha 2, and beta 1 bands in the frontal lobe and in beta 2 band in the frontal and occipital lobes. Patient 41 shows a greater power in delta (occipital and parietal areas) and theta (occipital and temporal regions) bands, and a lower activation in alpha 1 (frontal lobe), alpha 2 (frontal, parietal and temporal lobes), beta 1, and beta 2 bands in all lobes. Patients 51 and 71 reveal a widespread decrease of power density for alpha 1, alpha 2, beta 1, and beta 2 bands. As for patient 164, the results are the same except for the parietal lobe in alpha 1 band: in this case the power density increases. Patient 180 displays a widespread greater power in delta and theta bands. Patient 184 shows an increased activation in delta (temporal lobe) and theta (parietal and temporal lobes) bands; the power density decreases in the frontal and occipital regions for alpha 1 and alpha 2 bands, in the occipital lobe for beta 1 band, in the frontal, occipital, and parietal lobe for beta 2 band. The bar graph in Figure 1 summarizes the statistically significant results highlighted in gray in Table A1, Table A2 and Table A3. The graph reveals that most significant variations were detected in higher frequency bands. In particular, six patients show a decreasing power density in alpha 1 and alpha 2 (frontal lobe) bands, in addition to beta 2 (frontal and occipital lobes) band. As for delta and theta bands, the power density increases in four subjects at most, except for the temporal lobe in theta band: in this case, the increase is detected in five patients. Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9 show the boxplots of the power density values of each patient at time T0 and T1, for each frequency sub-band and each region of interest. From the analysis of the results, our hypothesis is that the outcome not consistent with previous findings are due to the passage from MCI to a mild AD condition. Therefore, the not yet severe stage of the disease, together with the inter-subject variability, could lead to the above-mentioned discrepancy of some outcome.

4. Discussion and Conclusions

Alzheimer’s disease is a neurological disorder which is the most common cause of dementia. As there is no cure for AD, early diagnosis is very important. In fact, a pharmacological treatment in the early stage (MCI) can slow down the evolution of the disease. The potential of EEG as a tool for studying AD has been widely explored and proved in the last years. In fact, the use of the EEG is advisable as compared to imaging techniques because it is easy-to-use, not expensive, comfortable, and safe for the patient.
Over the years, several studies have revealed that EEG signals of AD patients become slower, meaning that the spectral power is shifted towards lower frequencies [49]. In our paper, we conducted a longitudinal study involving eight patients who converted from MCI to AD after four months. The doctors made the diagnosis according to the MoCA scores. The novelty is that the spectral power analysis was carried out by means of eLORETA algorithm, which should provide a better source localization accuracy compared to LORETA. The previous works, both cross-sectional and longitudinal, were all carried out using LORETA [26,27,28,29,30,31,32,33,34,35,36]. The time period of follow-up was relatively short, so the condition at time T1 should be considered as slight worse in comparison with MCI condition at time T0. Finally, unlike the three longitudinal studies reported in [34,35,36], whose results are at group level, our work wants to investigate whether resting state EEG sources can be a marker to follow the progression of the disease in a single AD subject. In our study, five patients (03, 32, 41, 180, and 184) out of eight showed an increasing power in delta and theta bands. All patients, except for 180, exhibited a lower activation in alpha 1 and beta 2 bands. Finally, all patients, except for 03 and 180, revealed a decreased power in alpha 2 and beta 1 bands. Most of our findings are consistent with those of previous studies. On the other hand, the discrepancy of some outcomes could be related to the not yet severe stage of the disease at time T1. This topic needs to be further explored. As for the future development, it could be interesting to apply the eLORETA method to high-density EEG, which should improve the source localization accuracy [50]. In addition, it would be advisable to test the method on a larger number of patients for the validation of results. In order to further confirm the better performance of eLORETA as compared to both LORETA and sLORETA, as already proved in [46,47], a comparative methodological study could be carried out. Moreover, as in recent years several new algorithms for solving the EEG inverse problem have been developed, a comparison between them and eLORETA would be essential for testing the eLORETA effectiveness [51]. Finally, the combination of source reconstruction with complex network analysis could also provide a further instrument for a better understanding of the evolution of the disease. In conclusion, our work could represent a first step for more detailed investigations, in order to make the analysis of power current density a valid tool for evaluating the cognitive decline in AD.

Author Contributions

Conceptualization, S.D. and F.L.F.; Methodology, S.D. and F.L.F.; Formal Analysis, S.D.; Software, S.D.; Project administration, F.L.F.; Writing—original draft preparation, S.D.; Writing—review and editing, F.L.F.; Supervision, F.L.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors thank the doctors of IRCCS Centro Neurolesi Bonino-Pulejo of Messina (Italy) for their insightful comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. p-values resulting from the Wilcoxon rank sum test for delta and theta sub-bands.
Table A1. p-values resulting from the Wilcoxon rank sum test for delta and theta sub-bands.
Patient IDDeltaTheta
F. L.O. L.P. L.T. L.F. L.O. L.P. L.T. L.
Pt 038.75 × 10 222 1.12 × 10 235 2.82 × 10 180 7.55 × 10 264 1.54 × 10 23 1.49 × 10 161 5.13 × 10 13 4.16 × 10 48
Pt 322.38 × 10 301 6.10 × 10 207 02.12 × 10 240 2.73 × 10 8 2.90 × 10 244 1.82 × 10 234 2.54 × 10 160
Pt 413.22 × 10 14 2.78 × 10 55 2.37 × 10 18 3.382 × 10 29 06.27 × 10 101 8.97 × 10 6 5.55 × 10 5
Pt 5104.60 × 10 250 0004.60 × 10 250 00
Pt 718.15 × 10 184 2.31 × 10 232 4.56 × 10 303 1.37 × 10 148 7.94 × 10 319 5.59 × 10 222 6.23 × 10 309 1.13 × 10 134
Pt 1643.87× 10 261 9.26 × 10 214 4.17 × 10 210 1.55 × 10 147 01.39 × 10 168 4.16 × 10 189 2.62 × 10 171
Pt 18002.66 × 10 242 0004.60 × 10 250 00
Pt 1841.05 × 10 126 7.47 × 10 124 7.11 × 10 15 3.80 × 10 4 1.15 × 10 44 3.22 × 10 187 4.31 × 10 131 6.08 × 10 30
Table A2. p-values resulting from the Wilcoxon rank sum test for alpha 1 and alpha 2 sub-bands.
Table A2. p-values resulting from the Wilcoxon rank sum test for alpha 1 and alpha 2 sub-bands.
Patient IDAlpha 1Alpha 2
F. L.O. L.P. L.T. L.F. L.O. L.P. L.T. L.
Pt 037.81 × 10 70 0.05113.42 × 10 38 1.23 × 10 47 4.07 × 10 27 4.38 × 10 111 0.48368.77 × 10 8
Pt 322.30 × 10 63 7.24 × 10 36 0.65182.06 × 10 15 3.78 × 10 14 2.45 × 10 69 3.65 × 10 27 2.24 × 10 33
Pt 4104.60 × 10 250 1.37 × 10 97 5.93 × 10 88 02.07 × 10 118 8.614 × 10 90 5.95 × 10 5
Pt 5108.09 × 10 250 0005.44 × 10 242 4.74 × 10 242 0
Pt 711.69 × 10 302 4.72 × 10 206 5.86 × 10 90 1.25 × 10 150 02.97 × 10 225 1.63 × 10 164 1.55 × 10 248
Pt 1641.46 × 10 61 5.18 × 10 7 2.09 × 10 18 2.84 × 10 18 5.62 × 10 317 5.06 × 10 81 1.50 × 10 39 6.83 × 10 74
Pt 18001.81 × 10 9 0001.32 × 10 75 00
Pt 1841.45 × 10 68 4.03 × 10 92 2.72 × 10 127 1.16 × 10 31 4.80 × 10 42 1.28 × 10 217 2.05 × 10 90 0.7671
Table A3. p-values resulting from the Wilcoxon rank sum test for beta 1 and beta 2 sub-bands.
Table A3. p-values resulting from the Wilcoxon rank sum test for beta 1 and beta 2 sub-bands.
Patient IDBeta 1Beta 2
F. L.O. L.P. L.T. L.F. L.O. L.P. L.T. L.
Pt 030.00155.75 × 10 79 8.44 × 10 9 0.86070.47861.66 × 10 11 0.02794.77 × 10 6
Pt 321.02 × 10 21 8.53 × 10 13 7.93 × 10 40 5.12 × 10 11 1.17 × 10 34 0.00373.85 × 10 105 0.0128
Pt 4109.77 × 10 248 0004.60 × 10 250 00
Pt 517.31 × 10 268 3.02 × 10 228 6.65 × 10 241 01.83 × 10 191 6.56 × 10 229 9.64 × 10 273 1.41 × 10 321
Pt 716.83 × 10 260 6.19 × 10 232 6.19 × 10 232 1.02 × 10 176 1.58 × 10 196 9.14 × 10 229 03.09 × 10 141
Pt 1642.77 × 10 80 4.62 × 10 87 3.14 × 10 19 1.23 × 10 68 3.13 × 10 181 8.86 × 10 214 1.23 × 10 132 2.67 × 10 279
Pt 18000.43700001.84 × 10 149 00
Pt 1840.00165.54 × 10 200 2.60 × 10 55 6.59 × 10 45 1.45 × 10 33 1.05 × 10 26 2.67 × 10 4 1.47 × 10 78

References

  1. Weintraub, S.; Wicklund, A.H.; Salmon, D.P. The neuropsychological profile of Alzheimer disease. Cold Spring Harb. Perspect. Med. 2012, 2, a006171. [Google Scholar] [CrossRef]
  2. Sperling, R.A.; Aisen, P.S.; Beckett, L.A.; Bennett, D.A.; Craft, S.; Fagan, A.M.; Iwatsubo, T.; Jack, C.R., Jr.; Kaye, J.; Montine, T.J.; et al. Toward defining the preclinical stages of Alzheimer’s disease: Recommendations from the National Institute on Aging-Alzheimer’s Association workgroups on diagnostic guidelines for Alzheimer’s disease. Alzheimers Dement. 2011, 7, 280–292. [Google Scholar] [CrossRef] [Green Version]
  3. Albert, M.S.; DeKosky, S.T.; Dickson, D.; Dubois, B.; Feldman, H.H.; Fox, N.C.; Gamst, A.; Holtzman, D.M.; Jagust, W.J.; Petersen, R.C.; et al. The diagnosis of mild cognitive impairment due to Alzheimer’s disease: Recommendations from the National Institute on Aging-Alzheimer’s Association workgroups on diagnostic guidelines for Alzheimer’s disease. Alzheimers Dement. 2011, 7, 270–279. [Google Scholar] [CrossRef] [Green Version]
  4. McKhann, G.M.; Knopman, D.S.; Chertkow, H.; Hyman, B.T.; Jack, C.R., Jr.; Kawas, C.H.; Klunk, W.E.; Koroshetz, W.J.; Manly, J.J.; Mayeux, R.; et al. The diagnosis of dementia due to Alzheimer’s disease: Recommendations from the National Institute on Aging-Alzheimer’s Association workgroups on diagnostic guidelines for Alzheimer’s disease. Alzheimers Dement. 2011, 7, 263–269. [Google Scholar] [CrossRef] [Green Version]
  5. Rossini, P.; Di Iorio, R.; Vecchio, F.; Anfossi, M.; Babiloni, C.; Bozzali, M.; Bruni, A.; Cappa, S.; Escudero, J.; Fraga, F.; et al. Early diagnosis of alzheimer’s disease: The role of biomarkers including advanced EEG signals analysis. An IFCN-sponsored panel of Experts. Clin. Neurophysiol. 2020, 131, 1287–1310. [Google Scholar] [CrossRef]
  6. McGinnis, S.M. Neuroimaging in neurodegenerative dementias. In Seminars in Neurology; Thieme Medical Publishers: New York, NY, USA, 2012; Volume 32, pp. 347–360. [Google Scholar]
  7. Bozzali, M.; Serra, L.; Cercignani, M. Quantitative MRI to understand Alzheimer’s disease pathophysiology. Curr. Opin. Neurol. 2016, 29, 437–444. [Google Scholar] [CrossRef] [PubMed]
  8. Iaccarino, L.; Chiotis, K.; Alongi, P.; Almkvist, O.; Wall, A.; Cerami, C.; Bettinardi, V.; Gianolli, L.; Nordberg, A.; Perani, D. A cross-validation of FDG-and amyloid-PET biomarkers in mild cognitive impairment for the risk prediction to dementia due to Alzheimer’s disease in a clinical setting. J. Alzheimers Dis. 2017, 59, 603–614. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Pillai, J.; Sperling, M.R. Interictal EEG and the diagnosis of epilepsy. Epilepsia 2006, 47, 14–22. [Google Scholar] [CrossRef] [PubMed]
  10. Mammone, N.; De Salvo, S.; Ieracitano, C.; Marino, S.; Cartella, E.; Bramanti, A.; Giorgianni, R.; Morabito, F.C. Compressibility of High-Density EEG Signals in Stroke Patients. Sensors 2018, 18, 4107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Wolters, C.H.; Anwander, A.; Tricoche, X.; Weinstein, D.; Koch, M.A.; Macleod, R.S. Influence of tissue conductivity anisotropy on EEG/MEG field and return current computation in a realistic head model: A simulation and visualization study using high-resolution finite element modeling. NeuroImage 2006, 30, 813–826. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Morabito, F.C.; Campolo, M.; Ieracitano, C.; Ebadi, J.M.; Bonanno, L.; Bramanti, A.; Desalvo, S.; Mammone, N.; Bramanti, P. Deep convolutional neural networks for classification of mild cognitive impaired and Alzheimer’s disease patients from scalp EEG recordings. In Proceedings of the 2016 IEEE 2nd International Forum on Research and Technologies for Society and Industry Leveraging a better tomorrow (RTSI), Bologna, Italy, 7–9 September 2016; pp. 1–6. [Google Scholar]
  13. Mammone, N.; Ieracitano, C.; Adeli, H.; Bramanti, A.; Morabito, F.C. Permutation Jaccard distance-based hierarchical clustering to estimate EEG network density modifications in MCI subjects. IEEE Trans. Neural Networks Learn. Syst. 2018, 29, 5122–5135. [Google Scholar] [CrossRef] [PubMed]
  14. Mammone, N.; De Salvo, S.; Bonanno, L.; Ieracitano, C.; Marino, S.; Marra, A.; Bramanti, A.; Morabito, F.C. Brain network analysis of compressive sensed high-density EEG signals in AD and MCI subjects. IEEE Trans. Ind. Inform. 2018, 15, 527–536. [Google Scholar] [CrossRef]
  15. La Foresta, F.; Morabito, F.C.; Marino, S.; Dattola, S. High-density EEG signal processing based on active-source reconstruction for brain network analysis in alzheimer’s disease. Electronics 2019, 8, 1031. [Google Scholar] [CrossRef] [Green Version]
  16. Ieracitano, C.; Mammone, N.; Hussain, A.; Morabito, F.C. A novel multi-modal machine learning based approach for automatic classification of EEG recordings in dementia. Neural Netw. 2020, 123, 176–190. [Google Scholar] [CrossRef] [PubMed]
  17. Huang, C.; Wahlund, L.O.; Dierks, T.; Julin, P.; Winblad, B.; Jelic, V. Discrimination of Alzheimer’s disease and mild cognitive impairment by equivalent EEG sources: A cross-sectional and longitudinal study. Clin. Neurophysiol. 2000, 111, 1961–1967. [Google Scholar] [CrossRef]
  18. Moretti, D.V.; Babiloni, C.; Binetti, G.; Cassetta, E.; Dal Forno, G.; Ferreric, F.; Ferri, R.; Lanuzza, B.; Miniussi, C.; Nobili, F.; et al. Individual analysis of EEG frequency and band power in mild Alzheimer’s disease. Clin. Neurophysiol. 2004, 115, 299–308. [Google Scholar] [CrossRef]
  19. Van der Hiele, K.; Vein, A.; Reijntjes, R.; Westendorp, R.; Bollen, E.; Van Buchem, M.; Van Dijk, J.; Middelkoop, H. EEG correlates in the spectrum of cognitive decline. Clin. Neurophysiol. 2007, 118, 1931–1939. [Google Scholar] [CrossRef]
  20. Czigler, B.; Csikós, D.; Hidasi, Z.; Gaál, Z.A.; Csibri, É.; Kiss, É.; Salacz, P.; Molnár, M. Quantitative EEG in early Alzheimer’s disease patients—power spectrum and complexity features. Int. J. Psychophysiol. 2008, 68, 75–80. [Google Scholar] [CrossRef]
  21. Gasser, U.S.; Rousson, V.; Hentschel, F.; Sattel, H.; Gasser, T. Alzheimer disease versus mixed dementias: An EEG perspective. Clin. Neurophysiol. 2008, 119, 2255–2259. [Google Scholar] [CrossRef]
  22. Hsiao, F.J.; Wang, Y.J.; Yan, S.H.; Chen, W.T.; Lin, Y.Y. Altered oscillation and synchronization of default-mode network activity in mild Alzheimer’s disease compared to mild cognitive impairment: An electrophysiological study. PLoS ONE 2013, 8, e68792. [Google Scholar] [CrossRef]
  23. Jeong, J. EEG dynamics in patients with Alzheimer’s disease. Clin. Neurophysiol. 2004, 115, 1490–1505. [Google Scholar] [CrossRef] [PubMed]
  24. Babiloni, C.; Lizio, R.; Marzano, N.; Capotosto, P.; Soricelli, A.; Triggiani, A.I.; Cordone, S.; Gesualdo, L.; Del Percio, C. Brain neural synchronization and functional coupling in Alzheimer’s disease as revealed by resting state EEG rhythms. Int. J. Psychophysiol. 2016, 103, 88–102. [Google Scholar] [CrossRef] [PubMed]
  25. Marqui, R.P.; Michel, C.M.; Lehmann, D. Low-resolution electromagnetic tomography—A new method for localizing electrical activity in the brain. Int. J. Psychophysiol. 1994, 18, 49–65. [Google Scholar] [CrossRef]
  26. Babiloni, C.; Binetti, G.; Cassetta, E.; Cerboneschi, D.; Dal Forno, G.; Del Percio, C.; Ferreri, F.; Ferri, R.; Lanuzza, B.; Miniussi, C.; et al. Mapping distributed sources of cortical rhythms in mild Alzheimer’s disease. A multicentric EEG study. Neuroimage 2004, 22, 57–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Babiloni, C.; Binetti, G.; Cassetta, E.; Dal Forno, G.; Del Percio, C.; Ferreri, F.; Ferri, R.; Frisoni, G.; Hirata, K.; Lanuzza, B.; et al. Sources of cortical rhythms change as a function of cognitive impairment in pathological aging: A multicenter study. Clin. Neurophysiol. 2006, 117, 252–268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Babiloni, C.; Frisoni, G.B.; Pievani, M.; Vecchio, F.; Lizio, R.; Buttiglione, M.; Geroldi, C.; Fracassi, C.; Eusebi, F.; Ferri, R.; et al. Hippocampal volume and cortical sources of EEG alpha rhythms in mild cognitive impairment and Alzheimer disease. Neuroimage 2009, 44, 123–135. [Google Scholar] [CrossRef] [PubMed]
  29. Babiloni, C.; Lizio, R.; Carducci, F.; Vecchio, F.; Redolfi, A.; Marino, S.; Tedeschi, G.; Montella, P.; Guizzaro, A.; Esposito, F.; et al. Resting state cortical electroencephalographic rhythms and white matter vascular lesions in subjects with Alzheimer’s disease: An Italian multicenter study. J. Alzheimers Dis. 2011, 26, 331–346. [Google Scholar] [CrossRef]
  30. Babiloni, C.; Carducci, F.; Lizio, R.; Vecchio, F.; Baglieri, A.; Bernardini, S.; Cavedo, E.; Bozzao, A.; Buttinelli, C.; Esposito, F.; et al. Resting state cortical electroencephalographic rhythms are related to gray matter volume in subjects with mild cognitive impairment and Alzheimer’s disease. Hum. Brain Mapp. 2013, 34, 1427–1446. [Google Scholar] [CrossRef]
  31. Canuet, L.; Tellado, I.; Couceiro, V.; Fraile, C.; Fernandez-Novoa, L.; Ishii, R.; Takeda, M.; Cacabelos, R. Resting-state network disruption and APOE genotype in Alzheimer’s disease: A lagged functional connectivity study. PLoS ONE 2012, 7, e46289. [Google Scholar] [CrossRef]
  32. Gianotti, L.R.; Künig, G.; Lehmann, D.; Faber, P.L.; Pascual-Marqui, R.D.; Kochi, K.; Schreiter-Gasser, U. Correlation between disease severity and brain electric LORETA tomography in Alzheimer’s disease. Clin. Neurophysiol. 2007, 118, 186–196. [Google Scholar] [CrossRef]
  33. Nishida, K.; Yoshimura, M.; Isotani, T.; Yoshida, T.; Kitaura, Y.; Saito, A.; Mii, H.; Kato, M.; Takekita, Y.; Suwa, A.; et al. Differences in quantitative EEG between frontotemporal dementia and Alzheimer’s disease as revealed by LORETA. Clin. Neurophysiol. 2011, 122, 1718–1725. [Google Scholar] [CrossRef] [PubMed]
  34. Rossini, P.; Del Percio, C.; Pasqualetti, P.; Cassetta, E.; Binetti, G.; Dal Forno, G.; Ferreri, F.; Frisoni, G.; Chiovenda, P.; Miniussi, C.; et al. Conversion from mild cognitive impairment to Alzheimer’s disease is predicted by sources and coherence of brain electroencephalography rhythms. Neuroscience 2006, 143, 793–803. [Google Scholar] [CrossRef] [PubMed]
  35. Babiloni, C.; Frisoni, G.B.; Vecchio, F.; Lizio, R.; Pievani, M.; Cristina, G.; Fracassi, C.; Vernieri, F.; Rodriguez, G.; Nobili, F.; et al. Stability of clinical condition in mild cognitive impairment is related to cortical sources of alpha rhythms: An electroencephalographic study. Hum. Brain Mapp. 2011, 32, 1916–1931. [Google Scholar] [CrossRef] [PubMed]
  36. Babiloni, C.; Lizio, R.; Del Percio, C.; Marzano, N.; Soricelli, A.; Salvatore, E.; Ferri, R.; Cosentino, F.I.; Tedeschi, G.; Montella, P.; et al. Cortical sources of resting state EEG rhythms are sensitive to the progression of early stage Alzheimer’s disease. J. Alzheimers Dis. 2013, 34, 1015–1035. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Hämäläinen, M.S.; Ilmoniemi, R.J. Interpreting magnetic fields of the brain: Minimum norm estimates. Med. Biol. Eng. Comput. 1994, 32, 35–42. [Google Scholar] [CrossRef]
  38. Talairach, J.; Tournoux, P. Co-Planar Stereotaxic Atlas of the Human Brain; Theime: New York, NY, USA, 1988. [Google Scholar]
  39. Collins, D.L.; Neelin, P.; Peters, T.M.; Evans, A.C. Automatic 3D intersubject registration of MR volumetric data in standardized Talairach space. J. Comput. Assist. Tomogr. 1994, 18, 192–205. [Google Scholar] [CrossRef]
  40. Tong, S.; Thakor, N.V. Quantitative EEG Analysis Methods and Clinical Applications; Artech House: Norwood, MA, USA, 2009. [Google Scholar]
  41. Tikhonov, A.N.; Arsenin, V.Y. Solutions of Ill-Posed Problems; Wiley: New York, NY, USA, 1977; pp. 1–30. [Google Scholar]
  42. Pascual-Marqui, R.D. Standardized low-resolution brain electromagnetic tomography (sLORETA): Technical details. Methods Find. Exp. Clin. Pharmacol. 2002, 24, 5–12. [Google Scholar]
  43. Mazziotta, J.; Toga, A.; Evans, A.; Fox, P.; Lancaster, J.; Zilles, K.; Woods, R.; Paus, T.; Simpson, G.; Pike, B.; et al. A probabilistic atlas and reference system for the human brain: International Consortium for Brain Mapping (ICBM). Philos. Trans. R. Soc. Lond. Ser. Biol. Sci. 2001, 356, 1293–1322. [Google Scholar] [CrossRef]
  44. Fuchs, M.; Kastner, J.; Wagner, M.; Hawes, S.; Ebersole, J.S. A standardized boundary element method volume conductor model. Clin. Neurophysiol. 2002, 113, 702–712. [Google Scholar] [CrossRef] [Green Version]
  45. Dale, A.M.; Liu, A.K.; Fischl, B.R.; Buckner, R.L.; Belliveau, J.W.; Lewine, J.D.; Halgren, E. Dynamic statistical parametric mapping: Combining fMRI and MEG for high-resolution imaging of cortical activity. Neuron 2000, 26, 55–67. [Google Scholar] [CrossRef] [Green Version]
  46. Pascual-Marqui, R.D. Discrete, 3D distributed, linear imaging methods of electric neuronal activity. Part 1: Exact, zero error localization. arXiv 2007, arXiv:0710.3341. [Google Scholar]
  47. Jatoi, M.A.; Kamel, N.; Malik, A.S.; Faye, I. EEG based brain source localization comparison of sLORETA and eLORETA. Australas. Phys. Eng. Sci. Med. 2014, 37, 713–721. [Google Scholar] [CrossRef] [PubMed]
  48. Gibbons, J.; Chakraborti, S. Nonparametric Statistical Inference; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  49. Dauwels, J.; Vialatte, F.; Cichocki, A. Diagnosis of Alzheimer’s disease from EEG signals: Where are we standing? Curr. Alzheimer Res. 2010, 7, 487–505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Dattola, S.; Morabito, F.C.; Mammone, N.; La Foresta, F. Findings about LORETA Applied to High-Density EEG—A Review. Electronics 2020, 9, 660. [Google Scholar] [CrossRef] [Green Version]
  51. Jatoi, M.A.; Kamel, N.; Malik, A.S.; Faye, I.; Begum, T. A survey of methods used for source localization using EEG signals. Biomed. Signal Process. Control 2014, 11, 42–52. [Google Scholar] [CrossRef]
Figure 1. Bar graph of the statistically significant results highlighted in gray in Table A1, Table A2 and Table A3.
Figure 1. Bar graph of the statistically significant results highlighted in gray in Table A1, Table A2 and Table A3.
Applsci 10 05666 g001
Figure 2. Boxplot of the power current density of patient 03 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Figure 2. Boxplot of the power current density of patient 03 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Applsci 10 05666 g002
Figure 3. Boxplot of the power current density of patient 32 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Figure 3. Boxplot of the power current density of patient 32 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Applsci 10 05666 g003
Figure 4. Boxplot of the power current density of patient 41 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Figure 4. Boxplot of the power current density of patient 41 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Applsci 10 05666 g004
Figure 5. Boxplot of the power current density of patient 51 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Figure 5. Boxplot of the power current density of patient 51 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Applsci 10 05666 g005
Figure 6. Boxplot of the power current density of patient 71 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Figure 6. Boxplot of the power current density of patient 71 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Applsci 10 05666 g006
Figure 7. Boxplot of the power current density of patient 164 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Figure 7. Boxplot of the power current density of patient 164 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Applsci 10 05666 g007
Figure 8. Boxplot of the power current density of patient 180 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Figure 8. Boxplot of the power current density of patient 180 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Applsci 10 05666 g008
Figure 9. Boxplot of the power current density of patient 184 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Figure 9. Boxplot of the power current density of patient 184 at time T0 and T1 for each frequency sub-band and each region of interest (frontal, occipital, parietal, and temporal).
Applsci 10 05666 g009
Table 1. Demographic data of the recruited subjects.
Table 1. Demographic data of the recruited subjects.
Patient IDGenderAge
Pt 03M68
Pt 32M78
Pt 41M78
Pt 51F72
Pt 71F79
Pt 164M76
Pt 180F78
Pt 184F69

Share and Cite

MDPI and ACS Style

Dattola, S.; La Foresta, F. An eLORETA Longitudinal Analysis of Resting State EEG Rhythms in Alzheimer’s Disease. Appl. Sci. 2020, 10, 5666. https://doi.org/10.3390/app10165666

AMA Style

Dattola S, La Foresta F. An eLORETA Longitudinal Analysis of Resting State EEG Rhythms in Alzheimer’s Disease. Applied Sciences. 2020; 10(16):5666. https://doi.org/10.3390/app10165666

Chicago/Turabian Style

Dattola, Serena, and Fabio La Foresta. 2020. "An eLORETA Longitudinal Analysis of Resting State EEG Rhythms in Alzheimer’s Disease" Applied Sciences 10, no. 16: 5666. https://doi.org/10.3390/app10165666

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop