Abstract
Objective The balance of cross-frequency coupling (CFC) over within-frequency coupling (WFC) can build a nonlinearity index (NI) that encapsulates the over-excitation of information flow between brain areas and across experimental time. The present study investigated for the very first time how the Greek High Phenolic Early Harvest Extra Virgin Olive Oil (HP-EH-EVOO) versus Moderate Phenolic (MP-EVOO) and Mediterranean Diet (MeDi) intervention in people with Mild Cognitive Impairment (MCI) could affect their spontaneous EEG dynamic connectivity.
Methods Fourty three subjects (14 in MeDi, 16 in MP-EVOO and 13 in HP-EH-EVOO) followed an EEG resting-state recording session (eyes-open and closed) before and after the treatment. Following our dominant coupling mode model (DoCM), we built a dynamic integrated dynamic functional connectivity graph (iDFCG) that tabulates both the functional strength and the DoCM of every pair of brain areas.
Results Signal spectrum within 1-13 Hz and theta/beta ratio have been decreased in the HP-EH-EVOO group in both conditions. FIDoCM has been improved after the intervention across groups and conditions but was more prominent in HP-EH-EVOO group (p < 0.001). Finally,we revealed a significant higher post-intervention reduction of NI (ΔNITotal and α) for the HP-EH-EVOO compared to the MP-EVOO and MeDi groups (p < 0.0001).
Conclusions Long-term intervention with HP-EH-EVOO reduced the over-excitation of information flow in spontaneous brain activity.
Significance Our study confirms the alteration of signal spectrum of EEG rhythms and dominant coupling mode due to the intervention with HP-EH-EVOO nutrition protocol.
Highlights
Non-pharmaceutical intervention based on HP-EH-EVOO in MCI reduces the over-excitation of information flow
Non-pharmaceutical intervention based on HP-EH-EVOO in MCI increases the human brain flexibility
Reconfiguration of dominant coupling modes in EEG resting-state due to the intervention is modulated by alpha frequency
1. Introduction
Mild cognitive impairment (MCI) is characterized as a multi-factorial syndrome that depends on psychological, neurobiological and social factors that all contributed to a high risk for developing dementia (Petersen et al., 1999). Clinical judgement of MCI demands multiple sessions and multiple domain tests to detect a cognitive decline and impairment across single or various cognitive functions (Winblad et al., 2004). Distinctions over multiple (amnestic or non-amnestic) and single domain MCI are more than important bearing the fact that amnestic-MCI single-domain, amnestic-MCI multiple-domain and mild Alzheimer’s disease (AD) could represent a continuum process that leads to deviations from normal aging (Brambati et al., 2009).
Besides the heterogeneity of pathophysiological mechanisms underlying MCI subjects and their diagnostic misidentification, there are no pharmaceutical therapeutic options available and approved by FDA or EMA. The pharmaceutical interventions include acetylcholinesterase-inhibitors and memantine which are approved only for AD and not for MCI and other medications which are not yet approved such as, ginkgo biloba, anti-inflammatory drugs, statins, platelet aggregation inhibitors (trifusal), piracetam and others (Karakaya et al., 2013). However, clinical trials based on the aforementioned substances targeting to neurofibrillary tangle formation and amyloid-beta accumulation failed to show a high efficient performance measures for AD patients with both primary and secondary cognitive parameters.
Recent evidences from non-pharmaceutical interventions in MCI and dementia subjects open new avenues over alternative options of treatments (Vlachos and Scarmeas, 2019). Recent studies reported promising evidences of alternative non-pharmaceutical treatments for MCI based on natural substances (e.g. Extra Virgin Olive Oil (EVOO) (Tsolaki et al 2020), saffron (Tsolaki et al., 2016)) and dietary interventions (e.g. the Mediterranean diet (MeDi)) (Andrich et al., 2015);. The MeDi is not a prescribed diet by a clinician but a general eating habituation from people located geographically around the Mediterranean sea. The traditional Mediterranean diet is characterized by high consumption of olive oil, vegetables, legumes, fruits and unprocessed cereals, moderate to high consumption of fish, low consumption of meat and meat products, and low to moderate consumption of dairy products. It is also characterized by moderate wine consumption. Scarmeas et al., reported that the adaptation of MeDi diet for a long period of time improved various cognitive functions and decreased the risk to develop MCI and also to progress from MCI to AD (Scarmeas et al., 2009). Moreover, healthy dietary habituations are strongly correlated with positive cognitive outcomes (Scarmeas et al., 2018).
Positive evidences of EVOO-enriched diet against AD pathology have been evaluated in both animals and humans. Researchers explored if the consumption of EVOO-enriched diet from TgSwDI mice could reduce the accumulation of amyloid- and tau-related brain levels and also the related alterations of cognitive functions. They reported that a long-term (6 months) consumption of EVOO containing diet provided a protective mechanism against AD pathology (Qosa et al., 2015). A clinical study investigated the potential positive role of EVOO-rich MeDi in 285 subjects with high vascular risk. The authors reported an improved cognitive functioning compared to a control diet. Additionally, participants followed EVOO-rich MeDi had a less number MCI subjects at follow up that the controls while participants followed a MeDi + Nuts didn’t differ from the controls in terms of the number of MCI subjects (Martinez-Lapiscina et al., 2013) after 7 years follow up.
The main ingredients of EVOO are glycerides of fatty acids (98%), mainly of mono-unsaturated fatty acids (MUFA) and especially of oleic acid. The remaining 2% include a variety of “minor compounds” like phenolic alcohols, phenolic acids, lignans, flavones and importantly secoiridoids and their derivatives (Servili et al., 2002, 2013). Natural phenols can modulate cell redox state (Singh et al., 2008) via a direct activity over enzymes, proteins, receptors and various signalling pathways (Kim et al., 2014; Goszcz et al.,2017). Another positive impact of natural phenols is their interference with biochemical homeostasis (Abuznait et al., 2013; Grossi et al., 2013) which further has impact to epigenetic alterations of the chromatin (Ayissi et al., 2014; Declerck et al., 2016).
There is strong evidence that MeDi prevents neurodegeneration due to the consumption of vegetables, legumes, and olive oil riched in polyphenols (Feart et al., 2010, 2013; Scarmeas et al., 2006, 2009, 2011). Phytochemicals which are chemical compounds produced by plants like flavonols, polyphenols support different protective biological activity of the cell. The repertoire of this activity could be anti-oxidant, anti-allergic, anti-inflammatory, anti-viral, anti-carcinogenic and anti-proliferative (Middleton, 1998; Hollman et al., 1999; Eastwood,1999). MeDi is characterized complementary to fruits and vegetables consumption by a consumption of EVOO in amounts ranging from 25 to 50 g/day (López-Miranda et al., 2008; Pitozzi et al., 2012; Qosa et al., 2015). Studies reported that polyphenols act as an anti-aggregation agent against the desposition of amyloid-β1-42 peptide plaques supporting an anti-toxic environment (Rigacci et al., 2011), while oleocanthal is capable of altering the oligomerization state of oligomers of amyloid-β1-42 peptide, which are toxic (Pitt et al., 2009). Both substances could have a great impact over the modifying treatment of AD.
A link between cognitive performance and dietary habits has been evaluated (Otaegui-Arrazola et al., 2014) while neurodegenerative disorders and cognitive decline have been associated to oxidative stress (von Bernhardi et al., 2012). A simple hypothesis from the aforementioned evidence of beneficial properties of antioxidant-rich foods is that their consumption could reduce the oxidative stress which could further protect from neurodegenerative diseases. MeDi is full of anti-oxidant rich foods with main origin from plants and has been associated with a health status (Sofi et al., 2010).
Normally, cells can counteract the oxidative insults by regulating their homeostasis. In neurodegenerative conditions related to aging, cell’s capacity to maintain a redox balance decreases which further has a consequence to mitochondrial dysfunction, neuronal injury and accumulation of radicals. We know that oxidative stress increases with aging (Finkel and Holbrook, 2000) and it is a significant age factor that makes neuronal systems more vulnerable to various neurodegenerative diseases (McCord, 2000). Oxidative stress causes mitochondrial injury which further causes loss of oxidative metabolism and an overproduction of free radicals. Oxidative stress makes a significant contribution in many human diseases such as AD (Lyras et al.,1999) and Parkinson disease (Cohen et al.,1984).
Our MICOIL Pilot Study was the first-ever reported longitudinal double-blind study which administered Greek HP-EH-EVOO to people with MCI for one year. The purpose of MICOIL was to evaluate the effect of a Greek HP-EH-EVOO+MeDi (Group 1) compared to a Greek MP-EVOO+MeDi (Group 2) and MeDi only (Group 3) on the multi cognitive status of patients with MCI, who are community-dwelling, Greek-speaking elderly, with age ranging from 60 to 80 (Tsolaki et al., 2020). The Greek HP-EH-EVOO received (50mL/day) together with MeDi instructions, the GreekMP-EVOO received (50mL/ day) together with MeDi instructions and MeDI group received only the MeDi instructions. The evaluated hypothesis that was that long-term intervention with either HP-EH-EVOO or MP-EVOO could be correlated with a significant improvement of cognitive functions compared to MeDi, independently of the presence of APOE ε4 carrier.
In the present study, we analysed electroencephalographic (EEG) resting-state recordings (both eyes-open and eyes-closed conditions) from a subset of participants from the three groups compared to the previous study focusing only on neuropsychological estimates (Tsolaki et al., 2020). EEG activity was recorded before and after the intervention with the three groups. We adopted a dynamic functional connectivity analytic pathway with the incorporation of our dominant coupling mode model (DoCM). Based on our DoCM, we incorporated both within frequency (intra-frequency) and between frequencies (cross-frequency) functional interactions under the same model. The main outcome of DoCM is the construction of an integrated dynamic functional connectivity graph (iDFCG) that tabulates both the functional strength and the DoCM of every pair of brain areas across experimental time (Dimitriadis, 2018). When two brain areas communicate under the same frequency content (intra-frequency) then this exchange of information is called linear. However, cross-frequency interactions are nonlinear information pathways and are highly active when two functionally distinct brain areas want to interact and exchange information (Hyafil et al., 2015). The ratio of cross-frequency over intra-frequency interactions can build a nonlinearity index (NI) of information change in our brain. In a resting-state EEG study, we showed in a group of adolescents with schizophrenia spectrum disorders (SSDs) that NI is significant higher compared to an age-matched healthy group and this reconfiguration of dominant coupling modes was driven by the alpha frequency (Dimitriadis, 2020). Here, we will adopt the same methodology assuming that the three groups will show a reduction of their individual NI after the treatment and also group1 and 2 administered with either Greek HP-EH-EVOO or Greek MP-EVOO will show a stronger reduction of their NI compared to the MeDi group 3. We hypothesized that the change of NI which is interpreted as a reduction of cross-frequency coupling (CFC) non-linear interactions will be driven by alpha frequency band which is a frequency of arousal at resting-state conditions (Barry et al., 2007). Eyes-open and eyes-closed are not two equivalent baseline ones and should be treated differently. Eyes-closed condition can be used as an arousal baseline while eyes-open condition as an activation baseline.
2. Material and Methods
2.1 MICOIL Study Design
From 1 December 2016 to 30 August 2018, participants were recruited from the Memory & Dementia clinic of the 3rd and 1st Neurology Departments of Aristotle University of Thessaloniki, Greece, and from the two Day Centers of the Greek Association of Alzheimer’s Disease and Related Disorders (GADRDA). The MICOIL study was realised in accordance with the Declaration of Helsinki and was approved by the GADRDA scientific & ethics committee (25/ 240 2016) (Clinical Trials Registration Number: NCT03 241 362996)1. All participants fulfilled the Petersen criteria for MCI (Petersen et al., 2009). Information about inclusion and exclusion criteria, study design, APOE genotype, randomization and allocation procedure of participants’ selection, withdrawals and description of every protocol followed for each group, we forward an interested reader to the first report of the MICOIL study (Tsolaki et al., 2020).
2.2 Participants
Demographics of the subjects that participated on EEG data collections is given in Table 1. Following a Kruskal-Wallis non-parametric test, we didn’t detect any effect over the groups in terms of the age, gender ratio and education level.
2.3 Neuropsychological Assessment
The neuropsychological battery included the Greek version of Mini-Mental State Examination (MMSE) [27, 28] to assess the general cognitive function, Rivermead Behavioral Memory Test-Story Immediate and Delayed recall [29, 30] for episodic memory, Rey Osterrieth Complex Figure Test copy and delayed recall [31, 32] which measures visuospatial long-term memory and executive functions, Trail Making Test parts A & B [33], to examine visuospatial ability, attention and executive functions, Alzheimer Disease Assessment Scale-Cognition (ADAS-Cog) [34, 35] to assess the severity of cognitive dysfunction, Wechsler Memory Scales Digit Span Forward and Backward [36, 37] to assess attention and working memory, Letter and Category Fluency Test [38] for assessing phonemic and semantic fluency and Clock-drawing Test [39] which measures visuo-spatial orientation, understanding of verbal instructions, abstract thinking, planning, concentration, executive and visuo-spatial skills. Depressive symptoms were assessed by the Geriatric Depression Scale [40, 41] using a cut-off score of < 6 at baseline. We also used the Neuropsychiatric Inventory [42, 43] for the assessment of other neuropsychiatric symptoms, since it is a critical component for the evaluation of the MCI subjects because their distress can cause or exacerbate cognitive problems. For further details regarding the neuropsychological assessments, see the first paper focusing on neuropsychological estimates (Tsolaki et al., 2020).
2.4 EEG Recordings
EEG experiments are performed in cognitive intact elderly subjects and subjects with suspect of neurodegenerative disorder and cognitive deficits. To reduce the time of the EEG montage, we should use a minimum of 19 standard exploring electrodes placed according to basic international 10-20 system (i.e. Fp1, Fp2, F7, F3, Fz, F4, F8, T3, C3, Cz, C4, T4, T5, P3, Pz, P4, T6, O1, and O2). A higher number of EEG electrodes extending the basic 10-20 system is appreciated. For all EEG recordings, there was inclusion of two ear (A1 and A2). Resistance of the EEG electrodes should be lower than 5 KOhm.
EEG experiment included with the resting state EEG conditions:
Resting state eyes closed EEG recording for 5 minutes (ideally, most of these 5 minutes of EEG recordings should be characterized by a subject relaxed not showing voluntary or involuntary movements).
Resting state eyes open EEG recording for 5 minutes.
The sampling frequency was set to fs = 500 Hz.
EEG time series were re-referenced to the average reference electrode (Nunez and Srinivasan, 2006) before pre-processing steps. The program of experiments has been approved by the local ethical committee of Alzheimer Hellas (no. 25/21-6-2016).
2.5 Artifact reduction with independent component analysis (ICA) and Wavelet Decomposition
An important pre-processing step before the estimation of dynamic functional connectivity graph (dFCG) is the denoising of EEG recordings. We adopted an already established algorithmic artefact reduction method that combines independent component analysis (ICA) and Wavelet Decomposition. For further details, see our recent study which describes in detail the adopted denoising approach (Dimitriadis, 2020). Briefly, the majority of studies adopted ICA that produces N independent components (ICs) where N denotes the number of EEG sensors. Afterward, we rejected ICs as artifactual components based on the time course of ICs, the topology of the related weights per sensor and also estimating measures like entropy, kurtosis and skewness to further support our decision. The main drawback of this approach especially for resting-state is that we have to zeroing one IC as a whole even if any kind of artifact is detected in specific epochs covering a small percentage of the total experimental time. For that reason, we decided in our previous study to decompose the time course of artifactual IC in epochs of 1 sec using wavelet decomposition. This approach gives us an advantage over traditional approach avoiding rejecting an IC that encapsulates also critical true brain activity mixed with a number of artifactual epochs (Dimitriadis, 2020).
2.6 Signal Power Analysis
We estimated power spectral density (PSD) using pwelch MATLAB function independently for every EEG sensor for each subject before and after the intervention and in both conditions. Then, we averaged PSD estimations across EEG sensors to characterize every subject’s PSD profile for every condition and in both pre and post intervention period. The PSD analysis was realized up to 45 Hz. We also estimated the frontal theta/beta ratio as a marker of attentional control (Putman et al., 2014; Clarke et al.,2019).
2.7 The integrated dynamic functional connectivity graph (IDFCG) based on iPLV
In this section, we will describe briefly our dominant intrinsic coupling model (DoCM) demonstrated in the whole repertoire of functional neuroimaging modalities (Antonakakis et al., 2017; Dimitriadis et al., 2017a-c, 2018a-d, 2019; Dimitriadis, 2020; Marimpis et al., 2020). The majority of studies so far were investigating functional connectivity independently for every frequency band while they focused mainly on within-frequency (intra-frequency) interactions. Last years, an increased amount of research studies exploring also functional connectivity between brain areas oscillating on a different frequency the so-called cross-frequency coupling (CFC) or inter-frequency interactions. DoCM model incorporates all possible functional connectivity modes under a same framework. The main aim of DoCM model is to detect the dominant coupling mode between every pair of EEG sensors and across temporal segments. Here, we adopted intra-frequency phase-to-phase and inter-frequency phase-to-amplitude potential coupling modes (CFC). We estimated dynamic functional connectivity graph (dFCG) within and between the seven studying frequency bands {δ, θ, α1, α2, β1, β2, γ} defined, respectively, within the ranges {0.5–4 Hz; 4–8 Hz; 8–10 Hz; 10–13 Hz; 13–20 Hz; 20–30 Hz; 30–48 Hz}. For this computation, we employed the EEG activity from the 19 EEG sensors in both resting-state conditions. EEG recordings were bandpass filtered using a 3rd order zero-phase Butterworth filter employing the filtfilt MATLAB function. The width of the temporal window was set equal to 500ms (or 250 samples) and moved forward across experimental time with a step equals to 100 ms (50 samples) which can encapsulate both slow and fast oscillations (Dimitriadis et al., 2013a,b, 2015a,b, 2016a,b,c, 2017a,b, 2018b,c; Marimpis et al., 2020). We finally analyzed 75 secs across subjects and conditions as a common experimental time across the participants.
For every pair of EEG sensors and for every temporal segment, we estimated the seven intra-frequency phase-to-phase possible interactions, one for every frequency band and twenty-one inter-frequency phase-to-amplitude interactions between every possible pair of the seven frequency bands. Here, we adopted iPLV as the proper connectivity estimator for both intra and inter-frequency coupling modes as in our previous studies. In total, we estimated twenty-eight possible coupling modes leading to twenty eight dFCG per subject and condition. A dFCG is a 3D matrix of dimensions [temporal segments x sensors x sensors]. Then, we followed a surrogate analysis with main scope to detect the dominant coupling mode per pair of EEG sensors and for each temporal segment (for further details see Dimitriads, 2020). This approach leads to the integrated dFCG (DIFCG) that keeps both the strength and also the dominant coupling mode per pair of EEG sensors and across experimental time. The IDFCG is a pair of 3D matrices of size [temporal segments x sensors x sensors] where one keeps the strength and the other the dominant coupling mode. The strength is a value of range [0,1] based on the adopted connectivity estimator iPLV while the dominant coupling mode is encoded with an integer from 1 up to 28 e.g 1 for d, 2 for θ, 3 for α1, …, 26 for β1-β2, 27 for β1-γ and 28 for β2 - γ.
Fig. 1 is created in analogy to previous studies to further exemplify the concept of DoCM model in not familiar readers (Dimitriadis, 2020; Marimpis et al., 2020). Fig.1A illustrates how DoCM is defined for the first two temporal segment of eyes-closed condition from the first subject of group A between Fp1 and Fp2 EEG sensors. Following the important step of surrogate analysis, we revealed α1-β2 cross-frequency phase-to-amplitude coupling (PAC) as the dominant coupling mode (DoCM) among 28 possible coupling modes for the first two temporal segments for this particular EEG pair of sensors. On the right side of frequency-dependent pairs of time series, we demonstrated as a matrix the functional strength measured with iPLV OF every possible coupling mode. The main diagonal of this matrix tabulates the intra-frequency phase-to-phase coupling modes while the off-diagonal stores the cross-frequency phase-to-amplitude coupling modes. Fig.1B illustrates the temporal evolution of DoCM for the Fp1-Fp2 EEG pair. The color encodes the functional strength of the coupling while the y-axis refer to the detected DoCM. On the right side of this semantic time series, we showed the probability distribution (PD) of DoCM across experimental time. From this representation, α1-α1 intra-frequency phase-to-phase coupling was the most representative across DoCM. Both the flexibility Index (FI) and PD are estimated from semantic time series as the one demonstrated in Fig. 1B.
2.8 Semantic features derived from the evolution of DoCM
This section describes in detail the semantic features that can be extracted from the 2nd 3D tensor that preserves the DoCM across EEG sensor space and experimental time.
2.8.1. Flexibility index (FI)
We estimated Flexibility index (FI), a measure that quantifies the transition rate of DoCM between every pair of EEG sensors (Dimitriadis and Salis, 2017; Dimitriadis et al., 2018d; Dimitriadis, 2020). We employed the 2nd 3D tensor of the DIFCG that tabulates the semantic information of DoCM across the brain and experimental time to estimate FI. An example of what information this tensor tabulates is shown in Fig.1B.
This metric will be called hereafter FIDoCM and it is defined as: We estimated the ΔFIDoCM between pre and post intervention FIDoCM We counted a transition only when a coupling exist between two consecutive temporal segments. The term (temporal segments – 1) equals the pairs of temporal segments where a functional coupling mode exists. FI measure counts only the change of the DoCM across time and not the preferred transition between specific pairs of frequency coupling modes. FIDoCM gets higher values for higher “transitions” of DoCM between neighboring in time temporal segments. Fig. 1B illustrates how FIDoCM is estimated for the Fp1-Fp2 EEG pair. The outcome of FI estimation for every pair of EEG sensors is a matrix of size 19 × 19. Then, we estimated the nodal FI as the mean of every row of this matrix. Finally, we estimated the global FI by averaging the 19 nodal FI.
2.8.2. Spatiotemporal distribution of DoCM—Comodulograms
Based on the 2nd 3D DIFCG that keeps the semantic information of the preferred dominant coupling mode, we can tabulate in a frequencies × frequencies matrix the probability distribution (PD) of observing each of the DoCM frequencies across 7 (intra-frequency) + 21 (crossfrequency coupling) = 28 possible coupling modes. We estimated PD across sensor space between the EEG sensor pairs and also across temporal segments.
The spatiotemporal PD tabulated in a matrix is called hereafter comodulogram and an example is demonstrated in Fig. 1B (Antonakakis et al., 2016, 2017a, 2017b; Dimitriadis and Salis, 2017; Dimitriadis et al., 2018a, 2018b, 2018c).
2.8.3. Nonlinearity index (NI) based on DoCM
When two brain areas communicates within the same frequency (intra-frequency) then this communication is linear. However, when two brain areas communicate via cross-frequency coupling then this pathway is called nonlinear and plays a pivotal important role in inter-areal communication (Chen et al., 2009, 2010; Dimitriadis, 2020). The distinction between within-frequency coupling (linear coupling) and cross-frequency coupling (non-linear coupling) has been validated via biophysical modeling (Chen et al., 2009, 2010). Nonlinear information pathways are highly active when two brain areas want to exchange information (Hyafil et al., 2015).
As it was mentioned in section 2.7.2, we estimated the PD of DoCM across EEG sensor space and at every temporal segment. The PD is divided in two sections, the intra-frequency interactions (7 in total) which are tabulated in the main diagonal in the comodulogram (Fig.1B) and the cross-frequency interactions in the off diagonal (21 in total).
The ratio of the sum of this 21 cross-frequency PD values versus the 7 intra-frequency values define our nonlinearity index (NI) (Dimitriadis, 2020) which is described in eq. 2. In our study as in a previous one (Dimitriadis, 2020), we assumed that any improvement (reduction) of NI after the intervention should be leaded by the α frequency due to the resting-state condition. To quantify the driving role of α frequency, we estimated NI as the ratio of the sum of PDs between α1 and {α2,β1,β2,γ} (4 PDs related to 4 cross-frequency coupling pairs) and between α2 and {β1,β2,γ} (3 PDs related to 3 cross-frequency coupling pairs) with the sum of PDs related to α1 and α2 within frequencies interactions (2 PDs related to 2 intra-frequency couplings) (eq. 3). The higher the NI the higher is the contribution of CFC to the DoCM and so the higher is the nonlinear communication between brain areas. The outcome of this process is a time series of size equal to the number of temporal segments which is hereafter will be called dynamic NI (dNI). Two dNI were estimated per subject and condition, one incorporating the whole repertoire of coupling modes and one targeting to the α frequency as a modulating frequency. The eq.3 and 4 report the estimation of NITotal and NIα for one temporal segment. To quantify the level of potential difference of dNI between the pre and post-condition, we estimated the delta difference Δ of the median values of dNI.
2.9 Statistical Analysis
We estimated Cohen’s d effect size for PSD as a mean across frequency bin between pre and post-intervention period per group and in both conditions. A Wilcoxon Signed Ranks Test has been applied to the theta/beta ratio between pre and post intervention period per group. We also performed statistical analysis between pairs of groups for detecting group differences in terms of ΔNI. We also applied statistical analysis over PD of DoCM independently per group between pre and post intervention condition to uncover reorganization of DoCM due to intervention. For both statistical comparisons over ΔNI and PD of DoCM, we adopted the one-way ANOVA with p < 0.05 corrected for multiple comparison.
We estimated the delta (Δ) difference of each neuropsychological assessment as for FI and NI with the following formula (post - pre)/pre. ΔNITotal or α from the eyes-open condition were used as response variable and neuropsychological assessments as predictors following a multi-linear regression analysis. Clock Drawing, clock copy, digit span forward and digit span backward were excluded from the analysis due to many stable findings in pre and post intervention period (Δ = 0). Our analysis has been realized solely for the eyes-open condition due to the fact that neuropsychological assessments have been done with eyes-open. For that reason, eyes-open condition is more appropriate than eyes-closed.
2.10 Power Analysis
We estimated effect size and actual power analysis a posteriori on the estimated measurements of global ΔFITotal and ΔNITotal or α. As a statistical test, we employed one-way ANOVA omnibus with fixed effects. Power (1 - b err prob) was set to 0.95 and α err prob to 0.05. We repeated the analysis independently for the four conditions ({eyes-open – eyes-closed}x{ΔNITotal,ΔNIα}).
2.11 Software
The analysis has been realized in the MATLAB environment (v2019b) using the signal processing toolbox. Fast ICA has been adopted from the fieldtrip toolbox. Dynamic functional connectivity has been based on in-house software provided on our github website: https://github.com/stdimitr/time_varying_PAC and https://github.com/stdimitr/docm_model
3. Results
3.1 Alterations of PSD after the intervention period across groups and conditions
Fig.2A-C illustrates the group averaged PSD for eyes-closed condition across the three groups and similarly Fig.2D-F demonstrates the group averaged PSD for eyes-open condition. We didn’t detect any significant difference between pre and post intervention periods across groups and conditions. However, we observed important trends where PSD was reduced across the studying spectrum (up to 45 Hz) mostly in HP-EH-EVOO and MP-EVOO groups.
In both conditions, we can clearly see large effect size for the HP-EH-EVOO, medium to large for MP-EVOO and small to medium for MeDi within 1 - 13 Hz (Table 2). A large effect size in beta band has been observed for the MP-EVOO group, and a small to medium effect for the HP-EH-EVOO and MeDi group. Table 2 summarizes the mean and standard deviation of the Cohen’s d effect size across frequency bins within specific frequency bands, the delta, theta, alpha and beta frequencies.
Interestingly, the frontal theta/beta ratio was significantly different only in HP-EH-EVOO in the eyes-open condition (p = 0.0032, zval = -2.7255; eyes – closed: p = 0.06, zval = -1.53). The theta/beta ratio was reduced in both conditions in HP-EH-EVOO group.
3.3 Improved ΔFIDoCM for all the Groups
ΔFIDoCM has been improved in the three groups in both conditions. However, the improvements are more prominent in HP-EH-EVOO and MP-EVOO groups (Fig.4).We didn’t observe any difference of the level of ΔFIDoCM between the two conditions. However, there is a trend of higher ΔFIDoCM values for HP-EH-EVOO and MP-EVOO in eyes-open compared to eyes-closed condition. Mean and standard deviations (st.d.) of the ΔFIDoCM are tabulated in Table 3.
3.4 ΔNI Group Differences Across Resting-State Conditions
Group-averaged dNITotal and group-averaged dNIα for the three groups and the two conditions are illustrated in Fig.5 and 6, correspondingly. Even though spontaneous activity is completely different across subjects, we reported group-averaged time-series in order to inform interested readers about the fluctuation of NI over time. Finally, we reported group-averaged ΔNITotal or α in Fig.7 where there are significant differences between groups in both conditions and in both ΔNITotal or α. HP-EH-EVOO group showed significant higher post-intervention reduction of the NI level in total and in α in both conditions. Interestingly, group-averaged ΔNIα were higher for eyes-open condition compared to eyes-closed in the three groups. Complementary, we observed interesting patterns for the rest of potential modulating frequencies. ΔNIδ was increased after the intervention in both conditions and across groups (Fig.8A,D) while ΔNIθ,β were reduced after the intervention in both conditions and across groups (Fig.8B,C,E,F). We detected the following significant differences between groups : i) the increment of ΔNIδ was higher for MeDi compared to HP-EH-EVOO group in both conditions (Fig.8A,D) and ii) the reduction of ΔNIβ was higher for HP-EH-EVOO group compared to MeDi group in eyes-closed condition (Fig.8C).
3.3 Power Analysis estimated over ΔFI and ΔNI findings
Table 4 reports the power analysis and the effect size of our ΔFI and ΔNI findings. It is clearly that our observations and statistical test are well powered by the sample size.
3.5 Multi-linear Regression Analysis
Multi-linear regression analysis for ΔNITotal or α responses with the Δ differences of eight neuropsychological assessments as potential predictors revealed interesting findings. Below,we reported the R2,F-stattisti and the relevant pvalue for the two multilinear regression analyses:
R2=0.1758 F=2.4620 pvalue= 0.0365 for ΔNITotal
R2=0.1217 F=1.8574 pvalue = 0.1070 for ΔNIα
β coefficients of both multi-linear regression analyses are illustrated in Fig.9. Our analysis untangled an interesting finding linking the reduction of ΔNITotal with Δ differences of eight neuropsychological variables between the baseline and the follow up periods. Similar analysis with ΔNIα didn’t reach the statistical level. However, the trend of β coefficients was similar with ΔNITotal (Fig.9).
4. Discussion
The present study is the first investigation of the effect of HP-EH-EVOO versus MP-EVOO and MeDi on the EEG brain dynamics at resting-state of people with MCI. In our previous study, we reported for the very first time that a long-term intervention with HP-EH-EVOO or MP-EVOO was associated with significant improvement in cognitive function compared to the MeDi group (Tsolaki et al., 2020). We adopted here our DoCM model which incorporates into a dFCG both the dominant coupling mode and also its functional strength between every pair of EEG sensros and acros experimental time (Antonakakis et al., 2017; Dimitriadis et al., 2017a-c, 2018a-d, 2019; Dimitriadis, 2020; Marimpis et al., 2020). We detected a significant higher post-intervention reduction of NI (ΔNITotal and α) for the HP-EH-EVOO compared to the MP-EVOO and MeDi groups. Our results support the reduction of the over-excitation of information flow across groups with a higher impact for the HP-EH-EVOO group.
Signal spectrum analysis didn’t reveal any significant difference between pre and post intervention periods across groups and conditions. However, we untangled an important trend By estimating the effect size within basic brain rhyhms (delta, theta, alpha and beta). In both conditions, we have observed a large effect size for the HP-EH-EVOO, medium to large for MP-EVOO and small to medium for MeDi within 1-13 Hz. A large effect size in beta band has been observed for the MP-EVOO group, and a small to medium effect for the HP-EH-EVOO and MeDi group.
Delta waves originate in medial frontal cortical areas and also in the insula an area associated with signals from our body (Mesulam and Mufson, 1982), the nucleus accumbens, and also the tegmental brainstem area” (Knyazev, 2012). Delta frequency was linked to an internal inhibition when it is important to inhibit all the irrelevant information from other sources. Delta frequencies up to 4 Hz originated in the frontal and anterior cingulate cortex have been linked to the inhibition of distracting information supporting the attentional demands of a task (Harmony, 2013). Global delta synchrony is suggested to be the result of GABAergic project neurons originating from the thalamic reticular nucleus (Herrera et al., 2019). The lateral geniculate nucleus, a second thalamic nucleus with a baseline attentional demands like in the eyes-open condition has been linked to the reduction of the global delta synchrony (Weynard et al., 2001). In our study, we observed a reduction of delta signal spectrum in both conditions mostly for the HP-EH-EVOO and secondly for the MP-EVOO underlined by large and medium effect size, correspondingly. Complemetary, ΔNIδ was increased after the intervention in both conditions and across groups while the increment of ΔNIδ was higher for MeDi compared to HP-EH-EVOO group in both conditions.
It is well-known that theta synchrony is involved in working memory processes (Sauseng et al., 2010). A recent study compared EEG signal power between mind wandering and cognitive tasks. They found an increased theta and a reduced beta power during mind wandering compared to cognitive tasks in a healthy group (Van Son et al., 2019). Especially, the defined frontal theta/beta ratio was increased and varied during mind wandering reflecting a reduced top-down attentional control favoring thoughts. Another study, reported a reduced frontal theta/beta ratio in healthy controls compared to a group with cognitive performance anxiety (Putman et al., 2014). The increased frontal theta/beta ratio is an index of top-down failure of attentional control. Here, we observed a significant reduction of frontal theta/beta ratio only for HP-EH-EVOO group in eyes-open condition. This reflects an increment of attentional control in HP-EH-EVOO group which was reflected to improved neuropsychological estimates. Moreover, we reported a reduction of theta signal spectrum in both conditions mostly for the HP-EH-EVOO and secondly for the MP-EVOO underlined by large and medium effect size, correspondingly. Group-averaged ΔNIθ was decreased after the intervention in both conditions and across groups.
Global alpha synchrony is an index of a healthy resting wakefulness where a subject can process any information receiving from the real world (Klimesch,1999). Alpha desynchronization can be realized a) via the activation of the visual system (eyes open compared to eyes closed condition) mediated by the reticular system (Volavka et al.,1967) and also b) due to changes of the cortical network and thalamo-cortical communication (Schurmann and Basar, 2001; Fisch, 2006). We observed a reduction of alpha signal spectrum in both conditions mostly for the HP-EH-EVOO and secondly for the MP-EVOO underlined by large and medium effect size, correspondingly. Group-averaged ΔNIα were higher for eyes-open condition compared to eyes-closed in the three groups. ΔNIα was decreased after the intervention in both conditions and mostly for HP-EH-EVOO group who showed significantly higher post-intervention reduction compared to the rest two groups. A recent study reported that EEG alpha power was significantly enhanced during mind wandering compared to a cognitive task underlying the its importance to explore temporal fluctuations of mind wandering (Compton et al., 2019).
Beta frequency power is involved at EEG resting-state networks in a healthy population (Mantini et al., 2007). Beta frequency has been linked to lower MMSE scores (Stam et al., 2003) while is part of an integrative biomarker that can predict the progression to AD (Poil et al., 2013). Complementary, a reduced beta power has been observed during mind wandering compared to cognitive tasks in a healthy group (Van Son et al., 2019). Here,we observed a reduction of beta signal spectrum in both conditions mostly for the HP-EH-EVOO and secondly for the MP-EVOO underlined by large and medium effect size, correspondingly. The reduction of ΔNIβ was higher for HP-EH-EVOO group compared to MeDi group in eyes-closed condition.
Fundamental changes of signa power and connectivity were observed in al the studied brain disorders/diseases with the direction of changes to be associated with the task (Basar et al., 2013). For that reason, brain oscillations can be used as possible biomarkers in clinical populations studies and also for the evaluation of treatment and intervention protocols. It is highly important to compare similar conditions and especially for resting-state conditions, eyes-closed should be trated as an arousal baseline while eyes-open condition as an activation baseline.
The fluctuation of dominant coupling modes quantified with FI and the related pre-post difference estimated with ΔFIDoCM showed an improvement in the three groups and in both conditions. We observed higher ΔFIDoCM values for HP-EH-EVOO and MP-EVOO in eyes-open compared to eyes-closed condition even though the findings didn’t reach significant level.
HP-EH-EVOO group showed significant higher post-intervention reduction of the NI level in total and in α in both conditions. Interestingly, group-averaged ΔNIα were higher for eyes-open condition compared to eyes-closed in the three groups. Complementary, we observed important to mention patterns for the rest of potential modulating frequencies. ΔNIδ was increased after the intervention in both conditions and across groups while ΔNIθ,β were both reduced after the intervention in both conditions and across groups. We detected the following significant differences between groups : i) the increment of ΔNIδ was higher for MeDi compared to HP-EH-EVOO group in both conditionsand ii) the reduction of ΔNIβ was higher for HP-EH-EVOO group compared to MeDi group in eyes-closed condition.
Following a multi-linear regression analysis with ΔNITotal or α as responses and the Δ differences of eight neuropsychological assessments as potential predictors untangled a link between reduction of ΔNITotal with Δ differences of eight neuropsychological variables. This is the very first study that reported a link between aberrant changes of DoCM quantified via the NI with neuropsychological estimates due to an intervention protocol.
The combination of the main key findings which are a) the increment of ΔFIDoCM and the reduction of post-intervention NI level in total and in α showed with negative values of ΔNItotal,α can support the following statement: Intervention protocol reduced the excitation of information flow expressed with the ratio of CFC over WFC which is driven mostly by α modulating frequency. A high ΔFIDoCM with a less involvement of CFC (reduction of NI) mostly modulating by α frequency means that dominant coupling modes fluctuate faster and are represented in a higher proportion over WFC compared to CFC.
This study complements the positive outcome of the first study that used HP-EH-EVOO. We further validated the improved cognitive measures published in the original MICOIL study (Tsolaki et al., 2020) by analysing the EEG resting-state recordings of a specific subgroup of subjects. Signal spectrum analysis and dynamic functional connectivity analysis was followed. Adopting our DoCM model, we explored how an intervention protocol tailored to MCI subjects can alter the multiplexity of dynamic functional connectivity in a positive direction (Buzsáki and Watson, 2012). The whole approach proved more sensitive to detect alterations of brain activity and connectivity compared to the trivial neuropsychological testing even on our small sample. The combination of neuropsychological estimates and EEG recordings is a promising low-cost and non-invasive monitoring assessment of inviduals at risk (Babiloni et al., 2010; Lizio et al., 2011)
Our study is novel and unique in both the intervention protocol and the adopted analytic pathway. However, we can report two basic limitations. The first one refers to our analysis using scalp EEG sensors. For that reason and in order to avoid and misleading interpretations of our findings in a local level, we reported our estimates across the whole EEG space. It would be very interesting to follow a similar protocol employing a high-density EEG system and work on virtual cortical space instead of scalp surface level. The second limitation is linked to the acceptable but small number of participants. Power analysis supported our findings but we could report that a higher number of participants will definitely be on the right direction.
5. Conclussions
In the present study, we evaluated for the very first time the effect of Greek High Phenolic Early Harvest Extra Virgin Olive Oil (HP-EH-EVOO) versus Moderate Phenolic (MP-EVOO) and Mediterranean Diet (MeDi) in people with MCI via the EEG resting-state analysis. Our analysis was unique in terms of combining both intra and cross-frequency interactions simultaneously under the DoCM model. We reported for the very first time a reconfiguration of phase driven DoCM in MCI subjects that followed a dietary protocol while this reconfiguration was more prominent for the HP-EH-EVOO group. Signal power reduced across the spectrum with most important findings up to 15 Hz for HP-EH-EVOO and MP-EVOO groups. Our analytic pathway can assist researchers of how they have to evaluate their intervention protocols tailored to MCI subjects but also in other target groups. Further analysis is needed to link our DoCM findings with the intervention protocol in a larger sample.
Data Availability
The dataset is not available to the public.
Author’s Contribution
SID: Conceptualization; Data curation; Formal analysis; Methodology; Writing - original draft
CL : EEG Data curation; Formal analysis; EEG Data collection. A.Ch,T.
EL. All Data collection (Demographics, Genetics, MRI, Blood examination
MC. Neuropsychological Examination
MT. Organization the methodology of the study, neurological examination of the patients, decision for the inclusion-exclusion criteria and final editing of the study
Acknowledgement
This work has been supported by the Greek Association of Alzheimer’s Disease and Related Disorders (GAADRD), theWorld Olive Centre for Health (WOCH) and Yanni’s Olive Grove company for the donation of HP-EH-EVOO: Yanni’s Fresh and MPEVOO: Yanni’s Selected. In addition, we would like to thank Panagiotis Diamantakos and Aimilia Riga kou for technical support in the analysis of EVOO, Hadar Halivni for editing the English language, and all clinicians who contributed to this study, as well as patients who took part in it. SD was supported by a MRC grant MR/K004360/1 (Behavioural and Neurophysiological Effects of Schizophrenia Risk Genes: A Multi-locus, Pathway Based Approach) and a MARIE-CURIE COFUND EU-UK Research Fellowship.
Footnotes
I have added the last author.