Impact Factor 3.508
2017 JCR, Clarivate Analytics 2018

Frontiers journals are at the top of citation and impact metrics

Frontiers in Neurology


Original Research ARTICLE

Front. Neurol., 17 October 2016 |

Magnetoencephalography for the Detection of Intervention Effects of a Specific Nutrient Combination in Patients with Mild Alzheimer’s Disease: Results from an Exploratory Double-Blind, Randomized, Controlled Study

imageElisabeth C. W. van Straaten1,2*, imageHanneke de Waal3, imageMarieke M. Lansbergen2, imagePhilip Scheltens3, imageFernando Maestu4, imageRafal Nowak5, imageArjan Hillebrand1 and imageCornelis J. Stam1
  • 1Department of Clinical Neurophysiology, MEG Center, VU Medical Center, Amsterdam, Netherlands
  • 2Nutricia Advanced Medical Nutrition, Nutricia Research, Utrecht, Netherlands
  • 3Department of Neurology, Alzheimer Center, VU Medical Center, Amsterdam, Netherlands
  • 4Laboratory of Cognitive and Computational Neuroscience (UCM-UPM), Center for Biomedical Technology, Madrid, Spain
  • 5Magnetoencephalography Unit, Centro Medico Teknon, Barcelona, Spain

Synaptic loss is an early pathological finding in Alzheimer’s disease (AD) and correlates with memory impairment. Changes in macroscopic brain activity measured with electro- and magnetoencephalography (EEG and MEG) in AD indicate synaptic changes and may therefore serve as markers of intervention effects in clinical trials. EEG peak frequency and functional networks have shown, in addition to improved memory performance, to be sensitive to detect an intervention effect in mild AD patients of the medical food Souvenaid containing the specific nutrient combination Fortasyn® Connect, which is designed to enhance synapse formation and function. Here, we explore the value of MEG, with higher spatial resolution than EEG, in identifying intervention effects of the nutrient combination by comparing MEG spectral measures, functional connectivity, and networks between an intervention and a control group. Quantitative markers describing spectral properties, functional connectivity, and graph theoretical aspects of MEG from the exploratory 24-week, double-blind, randomized, controlled Souvenir II MEG sub-study (NTR1975, in drug naïve patients with mild AD were compared between a test group (n = 27), receiving Souvenaid, and a control group (n = 28), receiving an isocaloric control product. The groups were unbalanced at screening with respect to Mini-Mental State Examination. Peak frequencies of MEG were compared with EEG peak frequencies, recorded in the same patients at similar time points, were compared with respect to sensitivity to intervention effects. No consistent statistically significant intervention effects were detected. In addition, we found no difference in sensitivity between MEG and EEG peak frequency. This exploratory study could not unequivocally establish the value of MEG in detecting interventional effects on brain activity, possibly due to small sample size and unbalanced study groups. We found no indication that the difference could be attributed to a lack of sensitivity of MEG compared with EEG. MEG in randomized controlled trials is feasible but its value to disclose intervention effects of Souvenaid in mild AD patients needs to be studied further.


Alzheimer’s disease (AD) is the leading cause of dementia, and although it has an incompletely understood etiology, synaptic connectivity seems to be reduced already in early stages (1, 2). Synapses are the main structures for functional connectivity between neurons, and loss of this connectivity is related to impaired cognition in AD (3). Therefore, optimizing synapse formation and function may counteract some of the cognitive effects of the AD process.

Synapses consist of a large part of neuronal membrane, and formation and function of neuronal membrane can be enhanced by increasing the availability of specific nutrients that are required for the synthesis of phospholipids [for an overview, see Ref. (4)]. In many in vitro and in vivo studies, the combined administration of these nutrients has been shown to increase brain levels of phospholipids and synaptic proteins, to increase synaptic density, to enhance cholinergic neurotransmission and receptor functioning, and to improve functional brain connectivity and cognitive performance [e.g., Ref. (511)].

The medical food Souvenaid contains the specific nutrient combination Fortasyn Connect which was designed to counteract synapse loss and dysfunction in AD by supplying precursors and cofactors that are essential for neuronal membrane formation and maintenance and that are believed to be insufficiently available in AD, i.e., docosahexaenoic acid (DHA), eicosapentaenoic acid (EPA), uridine (in the form of uridine monophosphate), choline, phospholipids, folic acid, vitamins B6, B12, C, E, and selenium.

Souvenaid is intended as a medical food for oral consumption under medical supervision with the purpose of addressing disease-specific nutrient requirements. It has been found to positively affect memory function in mild AD (1214).

Neurophysiological measures aid in the interpretation of the mode of action of an intervention (15). A previous clinical study showed preserved EEG spectral and graph theory-based functional network measures in patients receiving Souvenaid compared with a control group in a double-blind, randomized, controlled trial, indicating preserved connectivity as a macroscopically measurable objective intervention effect (14). Recent research suggests that the AD effect is not homogeneously distributed over different brain regions and that, in terms of brain networks, especially highly connected hub regions are affected (1618). More spatially detailed information in these hub areas could therefore potentially disclose additional information on the mode of action of the intervention. Magnetoencephalography (MEG) combines a high temporal resolution with a relatively high spatial resolution. MEG signals are hardly disturbed by the skull, and absolute values can be used without the need for a reference. This allows the study of brain activity in great spatial detail and a reliable transformation to anatomical space.

Magnetoencephalography recordings in AD patients have mainly shown loss of functional connectivity of fast (>8 Hz) oscillatory brain activity and a less optimal network structure (1822). However, no intervention studies have been performed in AD patients to improve the abnormalities seen with MEG recordings. We therefore set out to examine in an exploratory setting the feasibility and value of MEG to detect intervention effects of a medical food on brain activity in mild AD patients.

Materials and Methods

Study Design and Subjects

Fifty-five patients aged ≥50 years, meeting the criteria for AD according to the National Institute of Neurological and Communicative Disorders and Stroke and the Alzheimer’s Disease and Related Disorders Association (NINCDS-ADRDA) criteria and with a Mini-Mental State Examination (MMSE) ≥20 were enrolled in the Souvenir II MEG sub-study, which had an exploratory, randomized, controlled, double-blind, parallel-group design. Participants were 1:1 randomly allocated to either the test (Souvenaid, containing Fortasyn Connect) or the isocaloric control product (without Fortasyn Connect) as a 125-ml daily drink for 24 weeks based on a computer-generated randomization list. Patient inclusion and exclusion criteria were identical to those previously described in the Souvenir II study (14). Assessments were done at baseline, week 12, and week 24, including MEG, EEG, blood sampling, safety, and compliance based on diaries. Participants were recruited from one center in Amsterdam (n = 47), one center in Madrid (n = 2), and one center in Barcelona (n = 6). The study was conducted in accordance with the Declaration of Helsinki, the International Conference on Harmonization (ICH) guidelines for Good Clinical Practice (GCP) as appropriate for nutritional products, and local legislation of the country in which the research was conducted. The Dutch Trial Registration number for this study is NTR1975. The ethical review boards of the local study sites approved the study. Written informed consent was obtained from patients and caregivers.

MEG Recording and Post-Processing

Figure 1 shows an overview of the MEG recording and analysis steps. Twenty-minute eyes-closed resting-state MEG data were recorded using whole head MEG systems inside a magnetically shielded room in three centers (Table 1). The head position within the scanner was determined by using four attached head position indicator coils. The positions of the indicator coils and the head shape were determined before the recordings using a digitizer (3Space Fast-Track, Polhemus, Colchester, VT, USA). The recorded MEG data were transferred to the Amsterdam center for central post-processing. Channels with poor signal quality (maximum of 10) were manually de-selected after visual inspection for artifacts by a trained technician, and additionally a spatio-temporal Signal-Space Separation (tSSS) filter (Elekta data) was applied for artifact removal (10 s window, activity–noise correlation set at 0.9) using MaxFilter software (Elekta-Neuromag Oy version 2.2.10) (23). Head movement and number of discarded channels did not differ between recording centers.


Figure 1. MEG recording and analysis pipeline. Arrows indicate steps in time. Source-space MEG of the patient was reconstructed using the signal-space MEG time series, the patient’s MRI, and an anatomical atlas (input data). MEG was used for frequency analysis (with outcome parameters: peak frequency and relative power) and the construction of a functional connectivity-based (PLI) adjacency matrix. From this matrix, mean PLI as well as weighted network measures gamma and lambda were computed and the minimum spanning tree (MST) matrix was derived. From the MST matrix, network metrics betweenness centrality, eccentricity, leaf fraction, diameter, and tree hierarchy were computed. All outcome measures were compared between groups.


Table 1. Descriptive data of the recording MEG laboratories.

A beamformer method was used to transfer the signal-space data of the Elekta MEGs (n = 49 subjects) to source-space for further analysis: the MEG signals, recorded at the sensors, were reconstructed to MEG signals originating from 78 sources [according to the 78 automated anatomical labeling (AAL) areas] (2426). Therefore, magnetic resonance imaging (MRI) T1 series of the patients were coregistered on the MEG data using the current head position and digitized scalp, as well as the outline of the scalp as obtained from the MRI outline data. For nine subjects, no suitable MRI was available. For these subjects, the best fitting MRI was identified out of the MRIs available from this study and used for matching with satisfying matching results. The MEGs recorded using 4D Neuroimaging (Barcelona, n = 6) were analyzed separately in the original signal-space, due to incompatible software requirements for source-space computation. Five artifact-free epochs of 4096 samples (down-sampled four times, 13.1–16.4 s) were selected from each recording by a trained MEG technician and analyzed with in-house developed open access software (BrainWave, version 0.9.72, C. J. Stam,

For the source-space MEGs, outcome measures included measures for analysis A (characterized by one average value per subject per time point, see Statistical Analysis) and measures for analysis B (with one value for each AAL region per MEG study, see Statistical Analysis). The measures for analysis A included mean occipital peak frequency (for comparison with the previous Souvenir II EEG study and defined as the median frequency between 4 and 13 Hz, averaged across the signals of the occipital regions), relative power in the different frequency bands (averaged over all brain regions), functional connectivity [assessed with the Phase Lag Index (PLI) and averaged over all brain regions], and several functional network measures [mean normalized clustering coefficient (gamma), the averaged normalized shortest path length (lambda), and minimum spanning tree (MST) network-derived measures (leaf fraction, diameter, tree hierarchy, maximum betweenness centrality, and maximum eccentricity)]. The measures for analysis B included non-averaged region-specific peak frequency, relative power, functional connectivity, and MST betweenness centrality and eccentricity (2730). Analyses were performed in six frequency bands (delta 0.5–4 Hz, theta 4–8 Hz, alpha1 8–10 Hz, alpha2 10–13 Hz, beta 13–30 Hz, and gamma 30–48 Hz). Functional connectivity quantifies the functional interactions between brain regions with high values indicating a strong functional connection and low values indicating a weak or absent functional connection. The MSTs were constructed in such a way that all network nodes were included in the graph using the connections with the highest functional connectivity and without the formation of loops (29). Only descriptive statistics with respect to relative power and peak frequency were analyzed for the signal-space MEG subsample (Barcelona).

Post Hoc Analysis

In addition to MEG, the participants underwent 20-min eyes-closed task-free EEG according to the protocol used in the previous Souvenaid II study just before or after the MEG measurement (14). In short, digital EEG was recorded according to the 10–20 system with 21 electrodes and sample frequency 500 Hz (BrainRT, OSG Belgium). To compare source-space MEG’s sensitivity in detecting an intervention effect to that of EEG, we additionally analyzed peak frequency based on sensor-level EEG data from the same study sample (n = 49) using a common or average reference, online filter settings high pass 0.16 Hz and low pass 70 Hz. As for MEG, EEG peak frequency was computed for the parieto-occipital electrodes (P3, P4, O1, and O2) as the median frequency between 4 and 13 Hz. Only peak frequency was assessed since this parameter had the most clearly interpretable significant results in the previous Souvenir II EEG study and is (contrary to functional connectivity and network measures) not dependent on number of electrodes/sensors.

Safety and Compliance

Safety assessments included the examination of patient medical history, recording of (serious) adverse events [(S)AE], and the monitoring of vital signs and additional laboratory parameters (liver panel and renal function). Study product compliance was assessed using a daily diary, which was completed by the subject or caregiver at home.

Magnetoencephalography analyses were done on the intent-to-treat (ITT) MEG population (n = 49), defined as all randomized patients (ITT) with source-space reconstructed MEGs. Safety analyses were performed on the all-subjects treated (AST) population (n = 55), defined as all randomized patients who received at least one dose of study product.

Statistical Analysis

No previous study results were available on intervention in AD patients to improve abnormalities seen with MEG recordings on which the sample size calculation could be based. With the proposed sample size of 20 patients per group, a treatment difference of about 1 SD can be detected with a power of 90% using t-test and a two-tailed alpha of 0.05.

Source-space MEG outcome parameters that yield one value per subject per time point (Analysis A: mean occipital peak frequency, mean PLI, gamma, lambda, leaf fraction, diameter, tree hierarchy, maximum betweenness centrality, and maximum eccentricity) were analyzed using mixed models for repeated measures (MMRM) with change-from-baseline as an outcome, using SAS® software (SAS Enterprise Guide 4.3 for Windows, SAS Institute Inc., Cary, NC, USA) and in line with the previous Souvenir II EEG sub-study. More specifically, a mixed model by specifying the within subject unstructured variance–covariance structure and including group (treatment arm), time (as continuous variable; treatment duration), and the group × time interaction as fixed effects, adjusted for baseline values (baseline was included as covariate) was used. “Site” was not included as random effect, because MEG data for the primary MEG analysis (based on “source-space”) were only available from two sites and the majority of the data were collected at one single site. Normality of the distribution was checked, and no transformations to correct for lack of normality were needed. The two degrees of freedom contrast (for group and group × time) describing the difference in trajectories over time between Souvenaid and control groups was taken as the primary indication of treatment effect during the intervention period.

Source-space MEG outcome parameters that yield one value per AAL region per subject per time point (analysis B: peak frequency, relative power, PLI, MST betweenness centrality, and MST eccentricity) were analyzed with two different approaches: (1) to adhere to the analysis method of the previous Souvenir II EEG sub-study and to take into account covariates and evolution over time, a statistical analysis with average values for 10 predefined brain areas [consisting of 2 (left and right) times 5 (frontal, central, temporal, parietal, and occipital) regions] to examine whether an interaction of treatment and time was constant over the 10 brain areas with brain area as additional fixed factor and (2) to use the higher spatial resolution of the MEG, a permutation analysis with each of the 78 AAL regions, correcting for multiple comparisons to indicate the brain regions responsible for group differences. With this method, we derived a null distribution for between-group differences by permuting group assignment and calculating a t-statistic after each permutation with an independent t-test (31). The maximum t-value across ROIs of each permutation was used to construct a distribution of maximum t-values of 1000 permutations. To correct for multiple comparisons, we used the single threshold test with p-value set at 0.05.


Subject Flow

Figure 2 and Table 2 show the flow and demographic data of subjects throughout the study for both study groups, separately for the total ITT population and the ITT MEG analysis population. Five subjects (9%) withdrew from the study early (n = 2 subjects from the test group and n = 3 subjects from the control group). Reasons for early withdrawal were withdrawal of consent (n = 2, 1 in each group), occurrence of SAEs (n = 2, 1 in each group), and lost to follow-up (n = 1 in the control group). Subjects were well balanced over the study groups with respect to baseline characteristics, except for the total MMSE score. Despite the randomization procedure, mean total MMSE score at screening was significantly higher for the control group compared with the test group in both ITT populations (total ITT population: t-test, p = 0.021; ITT MEG analysis population: t-test, p = 0.032), indicating worse cognition at screening in the test group as compared with the control group.


Figure 2. Subject flow chart. *MEG analysis population, excluding six subjects from Barcelona.


Table 2. Baseline demographics and characteristics of the intent-to-treat (ITT) populations.

MEG Results

Peak frequency and relative power of the source-space data were not different between the groups over the 24-week intervention period with any of the statistical methods (Appendices A and B in Supplementary Material). Mean functional connectivity as measured with PLI was not different between groups (MMRM, Appendix C in Supplementary Material), but with regional analysis using permutation testing, the test group showed significant lower theta band values in the left gyrus rectus (permutation test; p < 0.05). The network measures were generally not different between the groups except that the MMRM showed a group difference over time for normalized clustering coefficient gamma in the alpha-2 band [MMRM; F(2,43) = 3.91, p = 0.028] and tree hierarchy in the theta band [MMRM; F(2,43) = 5.38, p = 0.008], but at endpoint no differences were observed between study groups (Figure 3; Appendices D–J in Supplementary Material). In addition, the permutation analysis showed a significantly lower theta band MST betweenness centrality at baseline in the left gyrus rectus in the test group compared with the control group (permutation test; p < 0.05). Descriptive statistics for mean peak frequency and relative power of the MEGs in signal-space (from Barcelona, n = 6) are presented in Appendices K and L in Supplementary Material.


Figure 3. Mean normalized clustering coefficient (gamma) in the alpha 2 band (A) and minimum spanning tree hierarchy in the theta band (B) at baseline, week 12, and week 24. Data are raw means ± SEM.

Post Hoc Results

In agreement with the MEG results, EEG peak frequency analysis yielded no group differences, suggesting that MEG is not necessarily less sensitive than EEG to detect an intervention effect.

Post hoc analyses were conducted to test the influence of gender on the intervention effect on the MEG outcome parameters. We could not demonstrate an intervention effect in males or females when tested separately.

As significant differences were observed in screening MMSE between the study groups, additional statistical analyses were conducted for all MEG parameters adjusting for screening MMSE. The results from these confounder analyses were similar to the results from the primary MEG analyses, indicating no straightforward confounding effects of screening MMSE.

Safety and Compliance

Regarding safety, no clinically relevant differences between groups were found for the blood safety parameters, (S)AE occurrence, and vital signs assessments. The calculated subject reported compliance during 24 weeks was high with 91% in the test group and 98% in the control group (Student’s t-test, p = 0.19).



This Souvenir II MEG sub-study was a 24-week study in 55 patients with mild AD to explore whether MEG has the potential to detect effects resulting from 24 weeks’ intake of Souvenaid compared with a control product. Based on results from a previous EEG study, we expected group differences in the peak frequency, delta band functional connectivity, and beta band network parameters gamma and lambda (14, 32). However, we observed no endpoint differences in the MEG measures between the test group and the control group. Although some differences in trajectory over time were found (gamma in the alpha2 band and MST tree hierarchy in the theta band), these findings were not confirmed by differences at endpoint and did not reveal a consistent pattern that was anticipated based on the previous EEG results, existing knowledge on AD pathology, and the hypothesized mode of action of Souvenaid. Due to the number of tests, we anticipate that these results may have resulted from type 1 errors. The study product was safe and well tolerated, as has been shown previously in mild and in mild-to-moderate AD (13, 14).

Relationship between Souvenaid, Synapses, and Electrophysiological Changes in AD

The main mechanism for neuronal communication is synaptic transmission of electrical current. Synapses as specialized parts of the cell membrane are built to facilitate transmission of current by continuous excitatory and inhibitory input from sending axons to receiving dendrites. As a result, postsynaptic potentials are constantly changing, and the net fluctuation of the electrical field can be measured as oscillations of the magnetic field with MEG and of the electrical field with EEG. The most consistent electrophysiological finding in AD is slowing of the oscillatory activity (33). The exact mechanism of this phenomenon is unknown, but model studies indicate that loss of synaptic connectivity plays a role (18, 34). In addition, breakdown of functional network integrity has also been found in EEG (35, 36) and could be modeled by decreasing coupling between neuronal assemblies in a computer model (37). The effects of the Souvenaid on brain function in early AD have been assessed before with EEG and preserved peak frequency and network organization were related to the nutritional intervention (14, 32). Given the effects of the Fortasyn Connect nutrients on synaptic structure and function and on brain connectivity in animal studies (5, 6, 10, 11) and the known changes on EEG in AD patients, it was hypothesized that improved synaptic function due to the intervention may be picked up by EEG as preserved oscillatory frequency and preserved functional network organization. In essence, MEG and EEG both measure the macroscopic effects of postsynaptic potential changes. It was, therefore, expected that the same neurophysiological principles and intervention effects apply to MEG.

Control Group Results

Our control group did not behave as expected. From previous neurophysiological studies, it is clear that AD is related to noticeable macroscopic neurophysiological changes: the peak frequency and functional connectivity decrease over time, and network parameters indicate a more random organization of functional connections (15). In our control group, we could not observe most of these changes during the course of the study, indicating that our control group was not typical during the course of the study. On the other hand, 24 weeks may be too short to detect neurophysiological decline in a relative small group of AD subjects and therefore longer study duration may reveal the typical AD-related MEG changes.

MEG in Relation to EEG

Although the choice of the imaging technique has a known influence on the results (38), we did not find indications of lower sensitivity of MEG compared with EEG. EEG was part of the study protocol, and in a post hoc analysis, we analyzed EEG peak frequency from the same study population at the same time points as the MEGs using the same software and comparable (i.e., occipital) cortical areas. Peak frequency is one of the most robust and most investigated neurophysiological measures in AD and behaves in a similar way in EEG and MEG with lower values in AD compared with controls (15, 20, 3941). In addition, EEG peak frequency has previously been identified as a sensitive measure for group differences after 24 weeks of intervention with Souvenaid (14). However, in the current study, EEG peak frequency did not differ between the test and control groups over the course of the study and at endpoint. Therefore, we could not attribute the lack of group difference to a lower sensitivity of MEG compared with EEG. On the other hand, MEG, which is more costly, more technical demanding, and less widely available than EEG, did not perform superior to EEG in this study. MEG may theoretically outperform EEG given several technical differences: magnetic fields detected with MEG do not suffer from filtering effects of surrounding tissue as much as the electric fields detected with EEG. Therefore, the estimation of sources, especially when located in deeper brain structures is more reliable and the maximum spatial and temporal resolution is higher. Another potential advantage of MEG over EEG is the fact that no referential recording point is needed. The use of a reference has been shown to influence the results, especially in EEG functional connectivity and network studies. Despite these advantages, we could not establish a higher sensitivity for an intervention effect over 24 weeks.


Several methodological limitations can be identified. The data from different centers were not compatible when using the latest analysis techniques. The field of MEG analysis is rapidly developing, and the demands on the recordings and the analysis strategies change accordingly. At this moment, no uniform MEG data format is available, and recordings from different scanner types cannot be pre-processed similarly. Not all state-of-the-art techniques are available and validated for merging data from different MEG systems. In the current study, we decided to implement the latest analysis technique for Elekta MEGs at the cost of excluding several subjects from the analyses. This has led to a reduction of power, but we feel that the accuracy of the source-space analysis method compensated for this loss (26). However, next steps toward a uniform MEG data format will likely increase the applicability of MEG in clinical trials.

In the current study, we did not control for time of the recording. Brain oscillations as measured with EEG and MEG are known to vary over the course of 24 h, with the lowest power of the signal early in the morning (4244). However, this diurnal rhythm applies to absolute oscillatory power and to a large extend is similar for all frequency bands. In our study, we only used relative frequency measures for the spectral analyses. It can, therefore, be expected that circadian effects will largely cancel out in these types of analyses, and we feel that the present results are not influenced substantially by differences in time of recording.

This study was exploratory in nature with small subject groups, which has some known drawbacks when interpreting group differences, some of them being the larger risk of non-uniformity of the groups and the lack of statistical power. Unfortunately, the patient groups were found to be unbalanced with respect to screening MMSE, despite a careful computerized randomization process. Covariate analyses correcting for the baseline differences did not reveal confounding effects. However, these adjustments assumed a linear relationship between the baseline characteristic and the outcome, whereas in reality, the relationship between MMSE at screening and the outcome parameters might have been more complex. We controlled for some patient characteristics that have a well-known effect on brain oscillations, such as use of benzodiazepines and structural brain lesions, including infarctions and tumors. However, other, unknown, factors, such as brain atrophy or small vessel disease, might have biased the results. It, therefore, remains unclear to what extent the intervention effect was masked by group differences already present at screening.


Possibly due to methodological issues (unbalanced and small study groups), we were not able to assess the potential of MEG-based measures as markers of an intervention effect of a medical food in mild AD patients. Furthermore, we could not replicate the EEG results from the well-powered Souvenir II study (n = 259) in which a significant effect of Souvenaid compared with control product was shown on peak frequency, functional connectivity, and brain network organization in patients with mild AD (14, 32). We did not demonstrate consistent statistically significant differences between groups for additional regional-specific MEG outcome measures nor did we observe the typical AD deterioration over time in the control group. Substantially increasing the sample size for better-balanced study groups might give indications in this direction.

Magnetoencephalography did not seem less sensitive than EEG but with a greater ability of in-depth regional analysis. The potential of MEG analysis might be demonstrated in future, larger, clinical studies when the field has progressed toward a more unified data format. Moreover, newly developed MEG analyses with hippocampal and centroid voxel MEG virtual electrodes as well as other (directed) connectivity and network analysis can be applied to further characterize brain magnetic activity in AD.

Author Contributions

ES was involved in the design of the study, writing of the protocol, acquisition, analysis, and interpretation of the data, and wrote the protocol. HW contributed to the acquisition, analysis, and interpretation of the data. ML participated in design of the study, writing of the protocol, analysis, and interpretation of the data. PS and CS made contributions to the design of the work and interpretation of the data. FM, RN, and AH were involved in the design, acquisition, and interpretation of the data. All authors contributed to and have approved the final manuscript.

Conflict of Interest Statement

ML is an employee of the funding body. PS has received grant support (for the institution) from GE Healthcare, Danone Research, Piramal, and MERCK. In the past 2 years, he has received consultancy/speaker fees (paid to the institution) from Lilly, GE Healthcare, Novartis, Forum, Sanofi, Nutricia, Probiodrug, and EIP Pharma. The remaining coauthors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


The MEG sub-study of the Souvenir II study (Dutch Trial Register NTR1975) was funded by Nutricia Research. The sponsors were involved in the collection, analysis, and interpretation of the data and in writing the manuscript as indicated by the affiliations.

Supplementary Material

The Supplementary Material for this article can be found online at


1. Scheff SW, Price DA, Schmitt FA, Mufson EJ. Hippocampal synaptic loss in early Alzheimer’s disease and mild cognitive impairment. Neurobiol Aging (2006) 27(10):1372–84. doi:10.1016/j.neurobiolaging.2005.09.012

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Selkoe DJ. Alzheimer’s disease is a synaptic failure. Science (2002) 298(5594):789–91. doi:10.1126/science.1074069

CrossRef Full Text | Google Scholar

3. Terry RD, Masliah E, Salmon DP, Butters N, DeTeresa R, Hill R, et al. Physical basis of cognitive alterations in Alzheimer’s disease: synapse loss is the major correlate of cognitive impairment. Ann Neurol (1991) 30(4):572–80. doi:10.1002/ana.410300410

PubMed Abstract | CrossRef Full Text | Google Scholar

4. van Wijk N, Broersen LM, de Wilde MC, Hageman RJ, Groenendijk M, Sijben JW, et al. Targeting synaptic dysfunction in Alzheimer’s disease by administering a specific nutrient combination. J Alzheimers Dis (2014) 38(3):459–79. doi:10.3233/JAD-130998

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Cansev M, van Wijk N, Turkyilmaz M, Orhan F, Sijben JW, Broersen LM. Specific multi-nutrient enriched diet enhances hippocampal cholinergic transmission in aged rats. Neurobiol Aging (2015) 36(1):344–51. doi:10.1016/j.neurobiolaging.2014.07.021

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Janickova H, Rudajev V, Dolejsi E, Koivisto H, Jakubik J, Tanila H, et al. Lipid-based diets improve muscarinic neurotransmission in the hippocampus of transgenic APPswe/PS1dE9 mice. Curr Alzheimer Res (2015) 12(10):923–31. doi:10.2174/1567205012666151027130350

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Koivisto H, Grimm MO, Rothhaar TL, Berkecz R, Lütjohann DD, Giniatullina R, et al. Special lipid-based diets alleviate cognitive deficits in the APPswe/PS1dE9 transgenic mouse model of Alzheimer’s disease independent of brain amyloid deposition. J Nutr Biochem (2014) 25(2):157–69. doi:10.1016/j.jnutbio.2013.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Savelkoul PJ, Janickova H, Kuipers AA, Hageman RJ, Kamphuis PJ, Dolezal V, et al. A specific multi-nutrient formulation enhances M1 muscarinic acetylcholine receptor responses in vitro. J Neurochem (2012) 120(4):631–40. doi:10.1111/j.1471-4159.2011.07616.x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Wiesmann M, Jansen D, Zerbi V, Broersen LM, Garthe A, Kiliaan AJ. Improved spatial learning strategy and memory in aged Alzheimer AβPPswe/PS1dE9 mice on a multi-nutrient diet. J Alzheimers Dis (2013) 37(1):233–45. doi:10.3233/JAD-130179

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Wiesmann M, Zerbi V, Jansen D, Haast R, Lütjohann D, Broersen LM, et al. A dietary treatment improves cerebral blood flow and brain connectivity in aging apoE4 mice. Neural Plast (2016) 2016:6846721. doi:10.1155/2016/6846721

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Zerbi V, Jansen D, Wiesmann M, Fang X, Broersen LM, Veltien A, et al. Multinutrient diets improve cerebral perfusion and neuroprotection in a murine model of Alzheimer’s disease. Neurobiol Aging (2014) 35(3):600–13. doi:10.1016/j.neurobiolaging.2013.09.038

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Olde Rikkert MG, Verhey FR, Blesa R, von Arnim CA, Bongers A, Harrison J, et al. Tolerability and safety of Souvenaid in patients with mild Alzheimer’s disease: results of multi-center, 24-week, open-label extension study. J Alzheimers Dis (2015) 44(2):471–80. doi:10.3233/JAD-141305

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Scheltens P, Kamphuis PJ, Verhey FR, Olde Rikkert MG, Wurtman RJ, Wilkinson D, et al. Efficacy of a medical food in mild Alzheimer’s disease: a randomized, controlled trial. Alzheimers Dement (2010) 6(1):1.e–10.e. doi:10.1016/j.jalz.2009.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Scheltens P, Twisk JW, Blesa R, Scarpini E, von Arnim CA, Bongers A, et al. Efficacy of Souvenaid in mild Alzheimer’s disease: results from a randomized, controlled trial. J Alzheimers Dis (2012) 31(1):225–36. doi:10.3233/JAD-2012-121189

PubMed Abstract | CrossRef Full Text | Google Scholar

15. van Straaten EC, Scheltens P, Gouw AA, Stam CJ. Eyes-closed task-free electroencephalography in clinical trials for Alzheimer’s disease: an emerging method based upon brain dynamics. Alzheimers Res Ther (2014) 6(9):86. doi:10.1186/s13195-014-0086-x

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Buckner RL, Sepulcre J, Talukdar T, Krienen FM, Liu H, Hedden T, et al. Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer’s disease. J Neurosci (2009) 29(6):1860–73. doi:10.1523/JNEUROSCI.5062-08.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Dai Z, Yan C, Li K, Wang Z, Wang J, Cao M, et al. Identifying and mapping connectivity patterns of brain network hubs in Alzheimer’s disease. Cereb Cortex (2015) 25(10):3723–42. doi:10.1093/cercor/bhu246

CrossRef Full Text | Google Scholar

18. Stam CJ, de Haan W, Daffertshofer A, Jones BF, Manshanden I, van Cappellen van Walsum AM, et al. Graph theoretical analysis of magnetoencephalographic functional connectivity in Alzheimer’s disease. Brain (2009) 132(Pt 1):213–24. doi:10.1093/brain/awn262

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Alonso JF, Poza J, Mananas MA, Romero S, Fernandez A, Hornero R. MEG connectivity analysis in patients with Alzheimer’s disease using cross mutual information and spectral coherence. Ann Biomed Eng (2011) 39(1):524–36. doi:10.1007/s10439-010-0155-7

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Franciotti R, Iacono D, Della Penna S, Pizzella V, Torquati K, Onofrj M, et al. Cortical rhythms reactivity in AD, LBD and normal subjects: a quantitative MEG study. Neurobiol Aging (2006) 27(8):1100–9. doi:10.1016/j.neurobiolaging.2005.05.027

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Ranasinghe KG, Hinkley LB, Beagle AJ, Mizuiri D, Dowling AF, Honma SM, et al. Regional functional connectivity predicts distinct cognitive impairments in Alzheimer’s disease spectrum. Neuroimage Clin (2014) 5:385–95. doi:10.1016/j.nicl.2014.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Stam CJ, van Cappellen van Walsum AM, Pijnenburg YA, Berendse HW, de Munck JC, Scheltens P, et al. Generalized synchronization of MEG recordings in Alzheimer’s disease: evidence for involvement of the gamma band. J Clin Neurophysiol (2002) 19(6):562–74. doi:10.1097/00004691-200212000-00010

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Medvedovsky M, Taulu S, Bikmullina R, Ahonen A, Paetau R. Fine tuning the correlation limit of spatio-temporal signal space separation for magnetoencephalography. J Neurosci Methods (2009) 177(1):203–11. doi:10.1016/j.jneumeth.2008.09.035

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Gong G, He Y, Concha L, Lebel C, Gross DW, Evans AC, et al. Mapping anatomical connectivity patterns of human cerebral cortex using in vivo diffusion tensor imaging tractography. Cereb Cortex (2009) 19(3):524–36. doi:10.1093/cercor/bhn102

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage (2002) 15(1):273–89. doi:10.1006/nimg.2001.0978

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Hillebrand A, Barnes GR, Bosboom JL, Berendse HW, Stam CJ. Frequency-dependent functional connectivity within resting-state networks: an atlas-based MEG beamformer solution. Neuroimage (2012) 59(4):3909–21. doi:10.1016/j.neuroimage.2011.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Boersma M, Smit DJ, Boomsma DI, De Geus EJ, Delemarre-van de Waal HA, Stam CJ. Growing trees in child brains: graph theoretical analysis of electroencephalography-derived minimum spanning tree in 5- and 7-year-old children reflects brain maturation. Brain Connect (2013) 3(1):50–60. doi:10.1089/brain.2012.0106

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Dijkstra EW. A note on two problems in connexion with graphs. Numerische Math (1959) 1(1):269–71. doi:10.1007/BF01386390

CrossRef Full Text | Google Scholar

29. Nešetřil J, Milková E, Nešetřilová H. Otakar Borůvka on minimum spanning tree problem: translation of both the 1926 papers, comments, history. Discrete Math (2001) 233(1–3):3–36. doi:10.1016/S0012-365X(00)00224-7

CrossRef Full Text | Google Scholar

30. Stam CJ, Nolte G, Daffertshofer A. Phase lag index: assessment of functional connectivity from multi channel EEG and MEG with diminished bias from common sources. Hum Brain Mapp (2007) 28(11):1178–93. doi:10.1002/hbm.20346

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Nichols TE, Holmes AP. Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum Brain Mapp (2002) 15(1):1–25. doi:10.1002/hbm.1058

PubMed Abstract | CrossRef Full Text | Google Scholar

32. de Waal H, Stam CJ, Lansbergen MM, Wieggers RL, Kamphuis PJ, Scheltens P, et al. The effect of souvenaid on functional brain network organisation in patients with mild Alzheimer’s disease: a randomised controlled study. PLoS One (2014) 9(1):e86558. doi:10.1371/journal.pone.0086558

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Jeong J. EEG dynamics in patients with Alzheimer’s disease. Clin Neurophysiol (2004) 115(7):1490–505. doi:10.1016/j.clinph.2004.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Abuhassan K, Coyle D, Maguire LP. Investigating the neural correlates of pathological cortical networks in Alzheimer’s disease using heterogeneous neuronal models. IEEE Trans Biomed Eng (2012) 59(3):890–6. doi:10.1109/TBME.2011.2181843

PubMed Abstract | CrossRef Full Text | Google Scholar

35. de Haan W, Pijnenburg YA, Strijers RL, van der Made Y, van der Flier WM, Scheltens P, et al. Functional neural network analysis in frontotemporal dementia and Alzheimer’s disease using EEG and graph theory. BMC Neurosci (2009) 10:101. doi:10.1186/1471-2202-10-101

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Vecchio F, Miraglia F, Marra C, Quaranta D, Vita MG, Bramanti P, et al. Human brain networks in cognitive decline: a graph theoretical analysis of cortical connectivity from EEG data. J Alzheimers Dis (2014) 41(1):113–27. doi:10.3233/JAD-132087

PubMed Abstract | CrossRef Full Text | Google Scholar

37. de Haan W, Mott K, van Straaten EC, Scheltens P, Stam CJ. Activity dependent degeneration explains hub vulnerability in Alzheimer’s disease. PLoS Comput Biol (2012) 8(8):e1002582. doi:10.1371/journal.pcbi.1002582

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Tijms BM, Wink AM, de Haan W, van der Flier WM, Stam CJ, Scheltens P, et al. Alzheimer’s disease: connecting findings from graph theoretical studies of brain networks. Neurobiol Aging (2013) 34(8):2023–36. doi:10.1016/j.neurobiolaging.2013.02.020

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Berendse HW, Verbunt JP, Scheltens P, van Dijk BW, Jonkman EJ. Magnetoencephalographic analysis of cortical activity in Alzheimer’s disease: a pilot study. Clin Neurophysiol (2000) 111(4):604–12. doi:10.1016/S1388-2457(99)00309-0

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Osipova D, Ahveninen J, Jensen O, Ylikoski A, Pekkonen E. Altered generation of spontaneous oscillations in Alzheimer’s disease. Neuroimage (2005) 27(4):835–41. doi:10.1016/j.neuroimage.2005.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Verdoorn TA, McCarten JR, Arciniegas DB, Golden R, Moldauer L, Georgopoulos A, et al. Evaluation and tracking of Alzheimer’s disease severity using resting-state magnetoencephalography. J Alzheimers Dis (2011) 26(Suppl 3):239–55. doi:10.3233/JAD-2011-0056

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Cacot P, Tesolin B, Sebban C. Diurnal variations of EEG power in healthy adults. Electroencephalogr Clin Neurophysiol (1995) 94(5):305–12. doi:10.1016/0013-4694(94)00298-Y

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Cummings L, Dane A, Rhodes J, Lynch P, Hughes AM. Diurnal variation in the quantitative EEG in healthy adult volunteers. Br J Clin Pharmacol (2000) 50(1):21–6. doi:10.1046/j.1365-2125.2000.00205.x

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Higuchi S, Liu Y, Yuasa T, Maeda A, Motohashi Y. Diurnal variations in alpha power density and subjective sleepiness while performing repeated vigilance tasks. Clin Neurophysiol (2001) 112(6):997–1000. doi:10.1016/S1388-2457(01)00527-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: magnetoencephalography, brain networks, clinical trial, Alzheimer’s disease, medical nutrition

Citation: van Straaten ECW, de Waal H, Lansbergen MM, Scheltens P, Maestu F, Nowak R, Hillebrand A and Stam CJ (2016) Magnetoencephalography for the Detection of Intervention Effects of a Specific Nutrient Combination in Patients with Mild Alzheimer’s Disease: Results from an Exploratory Double-Blind, Randomized, Controlled Study. Front. Neurol. 7:161. doi: 10.3389/fneur.2016.00161

Received: 03 May 2016; Accepted: 13 September 2016;
Published: 17 October 2016

Edited by:

Jesus Avila, Autonomous University of Madrid, Spain

Reviewed by:

Rik Henson, University of Cambridge, UK
Jose L. Cantero, Pablo de Olavide University, Spain

Copyright: © 2016 van Straaten, de Waal, Lansbergen, Scheltens, Maestu, Nowak, Hillebrand and Stam. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Elisabeth C. W. van Straaten,