Next Article in Journal
Real-World fNIRS Brain Activity Measurements during Ashtanga Vinyasa Yoga
Next Article in Special Issue
Functional Connectivity-Derived Optimal Gestational-Age Cut Points for Fetal Brain Network Maturity
Previous Article in Journal
Repetitive Transcranial Magnetic Stimulation as an Add-On Treatment for Cognitive Impairment in Alzheimer’s Disease and Its Impact on Self-Rated Quality of Life and Caregiver’s Burden
Previous Article in Special Issue
Time-Frequency Characterization of Resting Brain in Bipolar Disorder during Euthymia—A Preliminary Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Brain Sciences
Volume 11
Issue 6
10.3390/brainsci11060741
brainsci-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detection of Resting-State Functional Connectivity from High-Density Electroencephalography Data: Impact of Head Modeling Strategies

by
Gaia Amaranta Taberna
1,
Jessica Samogin
1,
Marco Marino
1,2 and
Dante Mantini
1,2,*
1
Research Center for Motor Control and Neuroplasticity, KU Leuven, 3001 Leuven, Belgium
2
Brain Imaging and Neural Dynamics Research Group, IRCCS San Camillo Hospital, 30126 Venice, Italy
*
Author to whom correspondence should be addressed.
Brain Sci. 2021, 11(6), 741; https://doi.org/10.3390/brainsci11060741
Submission received: 31 March 2021 / Revised: 23 May 2021 / Accepted: 31 May 2021 / Published: 3 June 2021
(This article belongs to the Special Issue Neural Networks and Connectivity among Brain Regions)
Browse Figures
Versions Notes

Abstract

:
Recent technological advances have been permitted to use high-density electroencephalography (hdEEG) for the estimation of functional connectivity and the mapping of resting-state networks (RSNs). The reliable estimate of activity and connectivity from hdEEG data relies on the creation of an accurate head model, defining how neural currents propagate from the cortex to the sensors placed over the scalp. To the best of our knowledge, no study has been conducted yet to systematically test to what extent head modeling accuracy impacts on EEG-RSN reconstruction. To address this question, we used 256-channel hdEEG data collected in a group of young healthy participants at rest. We first estimated functional connectivity in EEG-RSNs by means of band-limited power envelope correlations, using neural activity estimated with an optimized analysis workflow. Then, we defined a series of head models with different levels of complexity, specifically testing the effect of different electrode positioning techniques and head tissue segmentation methods. We observed that robust EEG-RSNs can be obtained using a realistic head model, and that inaccuracies due to head tissue segmentation impact on RSN reconstruction more than those due to electrode positioning. Additionally, we found that EEG-RSN robustness to head model variations had space and frequency specificity. Overall, our results may contribute to defining a benchmark for assessing the reliability of hdEEG functional connectivity measures.

1. Introduction

Electroencephalography (EEG) is a non-invasive neurophysiological technique that measures the variation in electrical potentials by means of electrodes positioned over the scalp [1]. This variation in potentials is associated with neural activity that is generated inside the brain, referred to as source space, and propagates to the surface, named sensor space [1]. The analysis of EEG recordings can be used as a proxy to better understand neural mechanisms in health and disease. So far, most EEG studies have relied on analyses conducted in the sensor space, to make inferences on neural responses to specific stimuli (event-related potential analysis) [2,3] or modulations of neural oscillations during task performance (event-related synchronization/desynchronization analysis) [4,5,6,7].
Over the last years, there has been a growing interest in using electrophysiological techniques to study functional interactions between distant brain regions, namely functional connectivity, especially in the resting-state condition. Different resting-state networks (RSNs), among which are the default mode (DMN), the dorsal (DAN) and the ventral (VAN) attention, the language (LN), the somatomotor (SMN) and the visual (VN) networks, have been detected first using magnetoencephalography [8,9] and then EEG [10,11,12,13]. The use of high-density EEG (hdEEG) systems capable of recording more than 100 signals at different scalp locations have been crucial to reliably measure functional connectivity and map RSNs [14,15,16,17]. Specifically, dedicated analysis workflows for hdEEG data were developed to improve the reconstruction of neural activity in the brain [10,11,12,18,19,20,21], including tools for signal pre-processing, head modeling and source localization. The use of an appropriate head modeling strategy is deemed fundamental to take advantage of the superior scalp coverage provided by hdEEG [22,23,24,25]. The localization of neural generators requires indeed the modeling of neural activity propagation from the sources, located inside the brain, to the potentials measured by EEG sensors over the scalp. The linear relationship between source activity and EEG measures is mathematically modeled by means of the leadfield matrix [17,26]. Notably, an accurate realistic head model has been found to have a direct positive impact on the performance of source localization algorithms [14,16,17,27,28,29]. Two factors primarily contribute to the accuracy of head modeling: head image segmentation and EEG electrode positioning.
The segmentation of the head image, which can be obtained using the magnetic resonance (MR) imaging technique, permits the modeling of the topological and physical properties of each head tissue compartment. To this end, different MR segmentation methods were proposed, whose accuracy depends on the space where such segmentation is performed—individual [22,25,30] or template [31,32,33,34]—and on the number of considered tissues, which typically ranges from 3 to 12 [10,22,25,30,35,36,37,38]. In a recent study, we introduced MR-TIM, a computational method for automated head image segmentation in 12 tissue classes [25]. Importantly, we demonstrated the superior performance of MR-TIM compared to existing methods, which can identify 12 tissue classes by warping a segmented template image to individual space [10], or can segment the real MR image into 3 to 5 tissue classes directly in the individual space [35,38].
The accuracy of electrode positions is also very important for the generation of the leadfield matrix. Electrode positions can be obtained before or after the EEG measurements using digitization [39,40], photogrammetry [41,42] and 3D scanning [43,44,45,46] techniques. These are influenced to a variable extent by single-electrode localization errors [47] as well as by errors related to the co-registration of the whole set of electrodes to the head shape [48]. To the best of our knowledge, no study has yet assessed the impact of different electrode positioning, as well as of head tissue segmentation techniques, on functional connectivity measured on source-reconstructed hdEEG data.
In this study, we used real hdEEG data to examine the RSN reconstruction performance that can be achieved thanks to an optimized version of our hdEEG analysis workflow [49], including newly-developed methods for realistic MR segmentation in 12 head tissue compartments [25] and for leadfield matrix generation by finite difference modeling [50]. The connectivity patterns that we obtained were used as a reference to assess to what extent EEG-RSNs are sensitive to head model inaccuracies. We first tested the impact of random and systematic errors in electrode positioning, and—from this analysis—we derived the potential reduction in RSN reconstruction performance associated with the use of 3D scanning and digitization techniques for electrode localization. We then extended our performance analysis to different head tissue segmentations based on MR data (with twelve or with three tissues, either based on the individual MR image or on a template image). This allowed us to evaluate whether electrode positioning or head tissue segmentation most strongly influences RSN reconstructions from hdEEG data. Finally, we assessed if all RSNs are influenced in a similar manner by the methods used for electrode positioning and head tissue segmentation, or if—alternatively—some RSNs are less robust than others.

2. Materials and Methods

2.1. EEG Experiment

Eyes-open resting-state data were collected in a cohort of 19 healthy young adults (age 28 ± 6 years, 5 males) and were previously used in other studies from our group [11,12,49]. Ethical approval was granted by the Ethics Committee of ETH Zurich; informed consent was obtained from all participants and the experiment was performed in accordance with the relevant guidelines and regulations. For each participant, we recorded hdEEG signals for 5 min at 1000 Hz sampling rate using a 256-channel HydroCel Geodesic Sensor Net by Electrical Geodesics (Eugene, OR, USA). Electromyograms and vertical and horizontal electrooculograms (EOG) were collected in addition to the hdEEG signals. Positions of the EEG sensors and of three landmarks (nasion, left and right preauricular points) were localized using a Geodesic Photogrammetry System (GPS) [42]. In a separate session, we acquired a T1-weighted whole-head anatomical image using a Philips Ingenia 3T MR scanner (Best, The Netherlands) with a turbo field echo sequence. Scanning parameters were: TR = 8.25 ms, TE = 3.8 ms, flip angle = 8°, voxel size = 1 mm3 isotropic.

2.2. Standard EEG Data Analysis

Functional connectivity was estimated from hdEEG data using an automated analysis workflow developed in the MATLAB® (MathWorks Inc., Natick, MA, USA) environment. The workflow consisted of four main steps: pre-processing of EEG signals; creation of individual head model; reconstruction of cortical activity; and seed-based connectivity analysis. Each step is briefly described in the following sections; for a more detailed explanation refer to [10,12,49].

2.2.1. EEG Signal Pre-Processing

In order to clean EEG recordings from noise and biological artefacts, we first identified channels with poor signal quality and corrected them by spatially interpolating the signals from their neighbors, as defined using the FieldTrip toolbox (http://www.fieldtriptoolbox.org/, accessed on 2 July 2019). We band-pass filtered the data in 1–80 Hz and then removed EOG, muscular, movement and biological artefacts using Independent Component Analysis (ICA). Independent components that were not classified as artefactual were linearly mixed to reconstruct artefact-free channel data [51]. Finally, the pre-processed EEG data were re-referenced using the Reference Electrode Standardization Technique [51,52,53,54].

2.2.2. Realistic Head Model Creation

An individual head model was reconstructed using the MR image of the participant and the EEG sensor positions. The individual MR image was segmented in 12 head compartments (skin; eyes; muscle; fat; spongy bone; compact bone; cortical grey matter; cerebellar grey matter; cortical white matter; cerebellar white matter; cerebrospinal fluid; brain stem) using the MR-TIM software [25]. The conductivity value of each tissue compartment was assigned based on previous literature [55]. Electrode positions were extracted from the GPS data, and projected over the head surface to minimize localization errors. Afterwards, a 3D regular 6 mm grid overlapping with the cortical/subcortical and cerebellar compartments was generated, to define all possible dipolar sources. Finally, a whole-head conductivity model was created using the optimized finite difference method described in [50]. Based on this model, the leadfield matrix, expressing the linear relationship between scalp EEG data and source-space neural activity, was calculated.

2.2.3. Source Activity Reconstruction

The pre-processed EEG signals were then combined with the realistic head model to compute the cortical three-dimensional distribution of current density using the exact low-resolution brain electromagnetic tomography (eLORETA) [56]. The eLORETA method is a weighted minimum norm inverse solution, where the weights are unique and the inverse solution provides exact localization for any point source in the brain. In this way, cortical neural activity was estimated with high-temporal resolution in a 6 mm homogeneous grid constrained to the grey matter.

2.2.4. Functional Connectivity

We reconstructed functional interactions within and between six large-scale RSNs: DMN, DAN, VAN, LN, SMN and VN. In order to estimate the frequency-specific interactions, we selected a total of 21 nodes within the main regions of the six RSNs (Figure 1 and Table S1), based on previous studies [49,57,58,59,60].
Node coordinates, defined in MNI space, were projected onto the participant’s cortical space. We defined 6 mm radius spherical regions of interest (ROIs) centered over the seed coordinates. Time courses of these seed ROIs were decomposed in the time–frequency domain using the short-time Fourier transform, with Hamming windows of 2 s and a 50% overlap between consecutive windows. Connectivity values were estimated between each pair of seeds using band-limited power envelope correlations. Prior to that, we performed a frequency-by-frequency orthogonalization [61] to remove spurious interactions at zero-lag due to signal leakage. Logarithmic-transformed signal-orthogonalized power time-courses were correlated and resulting correlations were converted to z-values using the Fisher’s transform [61,62]. The average connectivity between the seeds of the same network was defined as the intra-network connectivity. Similarly, the average connectivity between the seeds associated with two different networks was defined as the inter-network connectivity [63,64]. We examined the strength of network interactions corresponding to delta (1–4 Hz), theta (4–8 Hz), alpha (8–13 Hz), beta (13–30 Hz) and gamma (30–80 Hz) bands, by averaging the z-values associated with the frequencies within the relevant range. In this manner, we derived band-limited RSN correlation matrices. Finally, for each network and frequency band separately, we used a one-tailed Wilcoxon signed rank test to contrast intra- and inter-network connectivity values.

2.3. Impact of Head Modeling Strategies

Several test models were created, by varying the reference analysis workflow described above with respect to a single parameter at the time. These tests were focused on assessing electrode localization errors and head tissue segmentation accuracy, respectively. For each test model, the corresponding analysis workflow was re-run on the hdEEG data, to produce new connectivity results for: (1) band-limited RSN correlation matrices and (2) the contrast between intra- and inter-network connectivity for each RSN and frequency band. To test the reliability of each model, we calculated Spearman’s correlation values between the connectivity matrices reconstructed from the reference analysis workflow and those obtained from the specific test model. The correlation values were transformed to z-values using the Fisher transformation to improve data normality.

Generation of Test Models

In order to evaluate the impact of electrode localization errors on the connectivity results, we created eight sets of synthetic electrode positions by adding displacements to the reference electrode positions (Figure 2). These displacements differed in magnitude and type. In particular, we defined six possible rotations around the anteroposterior, the mediolateral and the longitudinal axes of the head. The amplitude of the rotation determined the displacement of each electrode, which was defined according to previous studies [45,48,65], and was set to approximately: 0.25, 0.5, 0.75 and 1 cm. Within each dataset, the same direction of rotation was applied to all the electrodes to simulate the co-registration error (“Systematic”), whereas a random rotation was applied to each electrode for the localization error (“Random”). Rotations were randomized across datasets.
As a second analysis step, we simulated the impact of electrode position errors that are typical of 3D scanning [43,44,45] and digitization techniques [39,40]. The synthetic set of electrode positions for 3D scanning (3D Scan) was created by considering a systematic error equal to 0.25 cm and no random error; for the digitization one (Digitizer) we included a systematic error of 0.5 cm and a random error of 0.25 cm [44,45,65] (Table S2). The amplitude of the rotation was constant for all the datasets, while the direction of the rotation was randomized across participants.
We also examined the impact of the head tissue segmentation on the connectivity results, by using three different segmentation approaches with decreasing levels of complexity (Figure 3). First, we applied the warping of template segmentation (WTS) technique [10,25] on the individual MR image to obtain a “12-layer WTS” segmentation, with the same compartments defined by MR-TIM [25]. Then, the “3-layer template” segmentation was generated by clustering different compartments of the 12-layer template (brain: cortical and cerebellar white and grey matter, brainstem, cerebrospinal fluid; skull: compacta and spongiosa; scalp: muscle, fat, eyes, skin) [66]. Finally, we used the WTS technique to generate a “3-layer WTS” segmentation in individual space, based on the 3-layer template segmentation (Figure S1). For the 12-layer and 3-layer WTS we used the reference electrode positions in individual space; conversely, we used template electrode positions of the 256-channel HydroCel Geodesic Sensor Net for the segmented image in template space.
In addition, we created matrices contrasting intra- and inter-network connectivity. To this end, we used a one-tailed Wilcoxon signed rank test for the Spearman’s correlation values derived from 12-layer WTS, 3-layer WTS, 3-layer template, 3D Scan and Digitizer test models. These correlation values were obtained using the corresponding connectivity matrix derived from the standard analysis workflow as reference. With such comparisons, we aimed at understanding how robust the connectivity pattern of each network in each frequency band was when varying either the electrode positioning or the head tissue segmentation technique.

2.4. Statistical Analysis

For the analysis of electrode localization errors, we used a four-way analysis of variance (ANOVA) on the z-values to estimate to what extent the following factors influenced the connectivity values: network (DMN; DAN; VAN; LN; SMN; VN), frequency band (delta; theta; alpha; beta; gamma), error magnitude (0.25; 0.5; 0.75; 1 cm) and error type (systematic; random). A multiple comparison test was also performed to investigate the differences among estimated marginal group means, for each significant factor [67]. Furthermore, we compared correlation values corresponding to systematic and random localization errors to investigate in more detail under which conditions the two errors yield significantly different connectivity results. To this purpose, we used a two-tailed Wilcoxon signed rank paired test with significance level set to p < 0.05, corrected for multiple comparisons using the Bonferroni method. In order to test the impact of the head tissue segmentation method on connectivity, we performed a three-way ANOVA with factors: network (DMN; DAN; VAN; LN; SMN; VN), frequency band (delta; theta; alpha; beta; gamma) and segmentation method (12-layer WTS; 3-layer WTS; 3-layer template). Finally, we used a two-tailed Wilcoxon signed-rank paired test to check for differences in connectivity values between the head tissue segmentation and electrode positioning techniques. The significance level was set to p < 0.05, corrected for multiple comparisons using the Bonferroni method.

3. Results

In this study, we aimed at investigating the impact of the head model accuracy on EEG-RSN connectivity estimation, particularly focusing on the role of different electrode localization techniques and head tissue segmentation methods.
First, we reconstructed the functional connectivity values using our reference data. In particular, we quantified the connectivity between pairs of RSNs for each frequency band (Figure 4), and the differences between intra- and inter-network connectivity for each network and frequency band (Figure 5 and Figure S2). We observed remarkable similarities with connectivity results obtained in our recent study [49]. This was a valuable starting point for a detailed analysis of the impact of different head models on RSN connectivity.
The ANOVA performed on the z-values extracted for electrode positioning errors revealed significant modulations of all four factors under investigation (p < 0.001 for RSN, frequency band, error magnitude and type) (Table 1). Post-hoc tests showed that DMN, DAN and SMN were the most robust networks, whereas VAN and LN were the most sensitive to electrode positioning accuracy (Figure S3A). Similarly, alpha and delta bands were the most and least robust frequency bands, respectively (Figure S3B). As expected, our analyses confirmed that the difference between reference and test connectivity data increased as the magnitude of the positioning error was larger (Figure S3C). The Wilcoxon signed rank test yielded significantly lower correlations with random errors compared to systematic ones, for a magnitude level larger than 0.25 mm (Figure S3D and Figure 6 and Table S3).
The ANOVA run on the z-values obtained with different head tissue segmentation methods showed significant modulations (p < 0.001) of all three factors under investigation: RSN, frequency band and segmentation method (Table 2). Once again, post-hoc tests revealed that DMN, DAN and SMN were more robust in regard to head tissue segmentation variations than VAN and LN (Figure S4A), and that alpha and delta bands were the least and most sensitive frequency bands, respectively (Figure S4B). The 3-layer template method was found to be the least accurate. Notably, no significant differences were observed between the 12- and 3-layer WTS methods (Figure S4C).
When comparing the two factors contributing to head model inaccuracy, i.e., the electrode positioning technique and the head tissue segmentation method, we observed that the former one had significantly less impact on the estimation of the connectivity matrices (Figure 7). Indeed, the z-values for 3D scanning and digitizer techniques were always significantly higher (p < 0.001) than those obtained for the different head tissue segmentation methods (Table S4).
Lastly, we contrasted intra- and inter-network connectivity values for each resting-state network and frequency band, and examined the impact of head model inaccuracy (Figure S5). A correlation analysis conducted using the standard analysis workflow as a reference confirmed that the head tissue segmentation had a stronger impact on connectivity estimation than electrode positioning (Figure 8 and Figure S6). In line with other analyses conducted in this study (Figures S3 and S4), VAN and LN were found to be the most sensitive networks to variations in the head modeling strategy, whereas DAN was the most robust one.

4. Discussion

In recent years, an increasing number of studies have performed source-level functional connectivity analyses using hdEEG [1,10,18,19,49,68,69]. The reliability of the source activity reconstruction from hdEEG data, which is needed for estimating functional interactions between brain regions, does not depend solely on the density and coverage of the montage, but also on the workflow used for data analysis [14,16,17,27,70]. The primary goal of this study was to specifically examine the robustness of EEG-RSN reconstruction using different electrode localization and head tissue segmentation approaches.

4.1. Impact of Electrode Localization Error

The precision of electrode localization can be generally ascribed to two main factors: the single-electrode localization and the overall co-registration of the sensor positions over the head shape [48]. The results of our analyses show that random errors in electrode positioning (i.e., single-electrode localization errors) had generally a larger impact than the systematic errors (i.e., sensor co-registration errors) on the EEG connectivity results (Figure 6). The effect of random and systematic errors was significantly different for average displacements equal to 0.5 cm, with the former one clearly hindering the reconstruction of EEG connectivity patterns. For greater error magnitudes, instead, the effects of the two error types were more comparable, although systematic errors always had a more limited impact on EEG connectivity measures than random errors. Notably, for displacements as low as 0.25 cm there was no significant difference between the random and systematic errors, and relatively high correlations with our benchmark results. This finding is in line with previous literature, suggesting that electrode positioning errors smaller than 0.5 cm are negligible for an adequate brain source reconstruction [71,72].
When collecting electrode locations to be included in the individual head model, commonly-used electrode positioning techniques can be affected by both random and systematic errors, but their respective magnitude can vary consistently [41,42,45,46,48,65,73]. In this study, we modeled electrode positioning errors of the 3D scanning, assuming that the related measurements were affected by small systematic errors (0.25 mm) and that random errors were negligible. We also modeled electrode positioning errors of the digitizing technique as prone to both error types. The results of our study clearly show that the different measurement accuracy of these two techniques [44,45,48] resulted in evident differences in the EEG-RSN connectivity results (Figure 7 and Table S4).

4.2. Impact of Head Tissue Segmentation

The accuracy of the tissue segmentation can be related to the space in which the segmentation is performed (individual, as in MR-TIM, or template, as in WTS) and to the number of considered tissues. Our results evidence that warping a template segmentation on the individual MR image (12- or 3-layer WTS) [10] was a substantially less accurate approach than performing head tissue segmentation directly on the individual MR image itself [25] (Figure 7 and Table S4). We expected an improvement in performance when discriminating twelve tissue compartments (12-layer WTS) instead of only three (3-layer WTS). However, our tests did not reveal any significant difference when fewer tissues were included in the head model. This result may simply suggest that the WTS technique may not be sufficiently accurate under several conditions, and that it would be preferable to use MR-TIM for head tissue segmentation [25].

4.3. Differential Impact of Electrode Localization and Head Tissue Segmentation

As already mentioned, the head model used for reconstructing cortical activity combines information about electrode locations and head tissue distribution. In our study, we aimed to test which of these two factors has a stronger impact on EEG connectivity results. Our simulations yielded significantly higher correlations with the reference connectivity values when using different electrode positioning techniques (3D Scan and Digitizer), rather than when considering different head tissue segmentation methods (12- or 3-layer WTS) (Figure 7 and Table S4). In particular, suboptimal performance was achieved when using the template head model (3-layer template), created from the electrode positions and the template MR image defined in MNI space. This finding highlights the importance of including subject-specific information in the head model, namely real electrode positions and individual tissue segmentation [27,28,29].

4.4. Analysis of Robustness for Different Frequency Bands and Networks

The use of hdEEG has recently enabled the investigation of frequency-dependent functional interactions of several cortical regions and networks [10,18,49]. The results of our study clearly reveal that variations in EEG connectivity due to head model errors depend on the network and the frequency band considered (Figure 5 and Figures S2–S5). Specifically, the alpha band was the most robust in regard to head model variations, whereas the delta band was the most sensitive one (Figures S3 and S4). It should be noted that we analyzed the electrophysiological activity of participants in the resting state, a condition in which the alpha rhythm dominates the interactions among all brain regions [11,74]. Therefore, the robustness of the functional connectivity patterns in each band could vary according to the specific experimental conditions. In particular, the strongest brain oscillations are expected to be the least influenced by changes in the head model. We also found that DMN, DAN and SMN were relatively robust, whereas VAN and LN had large connectivity variations in the different test conditions (Figure 8 and Figures S3, S4 and S6). It should be considered that the number and the spatial distribution of network nodes over the cortex may influence the robustness of RSN reconstruction: VAN and LN are strongly lateralized networks with main nodes in the same hemisphere, whereas the other networks have core regions on both hemispheres. The robustness of the EEG-RSN reconstruction may also be influenced by other factors related to the RSN spatial distribution: DMN and DAN encompass deeper brain regions surrounded by uniformly distributed tissue interfaces, whereas key regions of VAN and LN are located closer to air–tissue interfaces, which are more likely to cause distortions in the electrical fields and therefore to affect the accuracy of the source reconstruction.

4.5. Study Limitations

A number of limitations in this study must be considered. First, we acknowledge the fact that we did not have the possibility to locate the EEG electrodes directly from the MR image, which would have assured error-free localization [73,75]. To minimize the effect of any electrode positioning error associated with the GPS technique, we projected the obtained electrode coordinates over the head shape [42,73] and manually corrected any apparent imprecision in electrode positioning. Additionally, we did not collect real electrode positions with several devices to compare our reference with; we therefore simulated positioning errors for 3D scanning and digitizing techniques, respectively, using error type and magnitude reported in previous studies. Similarly, a manually segmented MR image was not available to be used in our study. The analysis workflow that we used as a benchmark included MR-TIM, which was proven to outperform other commonly-used head tissue segmentation techniques [25]. We would also like to mention that resting-state EEG recordings used in this study were collected for 5 min in each participant [10,16], and it would be important to test if EEG-RSN reconstruction is stable for different EEG recording lengths. Furthermore, EEG connectivity results were based on time–frequency analyses conducted using short-time Fourier transform with uniform window length for all frequency bands, as in our previous studies [10,12,49]. The Wavelet transform could be alternatively used, such that the window length is adaptively defined based on the frequency; this may increase the sensitivity of the connectivity analyses for lower-frequency bands. Lastly, we assessed the impact of different head models on EEG connectivity by calculating the similarity of EEG connectivity results with respect to a standard configuration. This solution permitted us to conduct only relative comparisons among our test datasets; to the best of our knowledge, an absolute threshold level against which assessing the correspondence of two connectivity profiles can hardly be defined.

5. Conclusions

Network connectivity estimation from hdEEG data can be performed using an electrophysiological source imaging workflow that includes signal pre-processing, head model creation and brain activity reconstruction at the cortical level. In this study, we characterized to what extent head model variations associated with electrode positioning and head tissue segmentation influence EEG connectivity patterns. Our results suggest that inaccuracies in tissue segmentation have a stronger impact than those related to the electrode localization. Additionally, we observed that strongly lateralized networks, such as VAN and LN, were substantially less robust than bilateral networks. We hope that the present study may provide important information for researchers planning new EEG studies, supporting their choice of a suitable electrode positioning and head tissue segmentation approach, respectively. We suggest that future studies could perform test–retest analyses, which offer an alternative way to estimate the reliability of EEG connectivity results with respect to different head modeling approaches. Finally, it may be interesting to extend our findings by testing the impact of head modeling on EEG connectivity using other datasets, collected, for instance, in healthy older adults or in neurological patients.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/brainsci11060741/s1, Table S1: Definition of the 21 seeds used for analysis; Table S2: Definition of the simulated electrode localization errors due to different positioning techniques; Table S3: Correspondence between the reference connectivity matrices and those obtained when introducing electrode localization errors; Table S4: Comparative analysis of the z-values expressing the correspondence between the reference connectivity matrices and the matrices obtained when introducing errors in the head model; Figure S1: Representation of the head tissues in template space; Figure S2: Comparison of intra- and inter-network connectivity values for each frequency band and network; Figure S3: Estimated marginal means plots for each of the four factors included in the ANOVA test for electrode localization errors; Figure S4: Estimated marginal means plots for each of the three factors included in the ANOVA test for the data processed with different head tissue segmentation methods; Figure S5: Contrasts between intra- and inter-network connectivity, for different head models; Figure S6: Scatterplots of intra- and inter-network connectivity differences for different head models, compared to the corresponding values obtained using the reference analysis workflow.

Author Contributions

Conceptualization, D.M. and G.A.T.; Methodology, D.M., G.A.T. and J.S.; Software, G.A.T. and J.S.; Validation, D.M., G.A.T., J.S. and M.M.; Formal Analysis, D.M., G.A.T. and J.S.; Investigation, D.M., G.A.T., J.S. and M.M.; Resources, D.M.; Data Curation, G.A.T. and J.S.; Writing—Original Draft Preparation, G.A.T., J.S. and M.M.; Writing—Review and Editing, D.M.; Visualization, G.A.T. and J.S.; Supervision, D.M.; Project Administration, D.M.; Funding Acquisition, D.M. and M.M. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the KU Leuven Special Research Fund (grant C16/15/070), Research Foundation—Flanders (FWO) (grants G0F76.16 N, G0936.16 N, EOS.30446199 and postdoctoral fellowship 1211820N to M.M.) and Italian Ministry of Health (grant RF-2018-12366899). The computing resources and services used in the study were provided by the VSC (Flemish Supercomputer Center), funded by the FWO and the Flemish Government.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of ETH Zurich (EC 2013-N-06).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data are available from the corresponding author (D.M.) upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Michel, C.M.; Murray, M.M. Towards the utilization of EEG as a brain imaging tool. Neuroimage 2012, 61, 371–385. [Google Scholar] [CrossRef] [PubMed]
  2. Pfurtscheller, G.; Aranibar, A. Event-related cortical desynchronization detected by power measurements of scalp EEG. Electroencephalogr. Clin. Neurophysiol. 1977, 42, 817–826. [Google Scholar] [CrossRef]
  3. Pfurtscheller, G. Event-related synchronization (ERS): An electrophysiological correlate of cortical areas at rest. Electroencephalogr. Clin. Neurophysiol. 1992, 83, 62–69. [Google Scholar] [CrossRef]
  4. Frund, I.; Schadow, J.; Busch, N.A.; Korner, U.; Herrmann, C.S. Evoked gamma oscillations in human scalp EEG are test-retest reliable. Clin. Neurophysiol. 2007, 118, 221–227. [Google Scholar] [CrossRef] [PubMed]
  5. Morash, V.; Bai, O.; Furlani, S.; Lin, P.; Hallett, M. Classifying EEG signals preceding right hand, left hand, tongue, and right foot movements and motor imageries. Clin. Neurophysiol. 2008, 119, 2570–2578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Pfurtscheller, G.; Lopes da Silva, F.H. Event-related EEG/MEG synchronization and desynchronization: Basic principles. Clin. Neurophysiol. 1999, 110, 1842–1857. [Google Scholar] [CrossRef]
  7. Ruiz, Y.; Pockett, S.; Freeman, W.J.; Gonzalez, E.; Li, G. A method to study global spatial patterns related to sensory perception in scalp EEG. J. Neurosci. Methods 2010, 191, 110–118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Brookes, M.J.; Woolrich, M.; Luckhoo, H.; Price, D.; Hale, J.R.; Stephenson, M.C.; Barnes, G.R.; Smith, S.M.; Morris, P.G. Investigating the electrophysiological basis of resting state networks using magnetoencephalography. Proc. Natl. Acad. Sci. USA 2011, 108, 16783–16788. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. de Pasquale, F.; Della Penna, S.; Snyder, A.Z.; Lewis, C.; Mantini, D.; Marzetti, L.; Belardinelli, P.; Ciancetta, L.; Pizzella, V.; Romani, G.L.; et al. Temporal dynamics of spontaneous MEG activity in brain networks. Proc. Natl. Acad. Sci. USA 2010, 107, 6040–6045. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, Q.; Farahibozorg, S.; Porcaro, C.; Wenderoth, N.; Mantini, D. Detecting large-scale networks in the human brain using high-density electroencephalography. Hum. Brain Mapp. 2017, 38, 4631–4643. [Google Scholar] [CrossRef] [Green Version]
  11. Marino, M.; Arcara, G.; Porcaro, C.; Mantini, D. Hemodynamic Correlates of Electrophysiological Activity in the Default Mode Network. Front. Neurosci. 2019, 13, 1060. [Google Scholar] [CrossRef]
  12. Samogin, J.; Liu, Q.; Marino, M.; Wenderoth, N.; Mantini, D. Shared and connection-specific intrinsic interactions in the default mode network. Neuroimage 2019, 200, 474–481. [Google Scholar] [CrossRef]
  13. Sockeel, S.; Schwartz, D.; Pelegrini-Issac, M.; Benali, H. Large-Scale Functional Networks Identified from Resting-State EEG Using Spatial ICA. PLoS ONE 2016, 11, e0146845. [Google Scholar] [CrossRef] [PubMed]
  14. Hassan, M.; Dufor, O.; Merlet, I.; Berrou, C.; Wendling, F. EEG source connectivity analysis: From dense array recordings to brain networks. PLoS ONE 2014, 9, e105041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Klamer, S.; Elshahabi, A.; Lerche, H.; Braun, C.; Erb, M.; Scheffler, K.; Focke, N.K. Differences between MEG and high-density EEG source localizations using a distributed source model in comparison to fMRI. Brain Topogr. 2015, 28, 87–94. [Google Scholar] [CrossRef]
  16. Liu, Q.; Ganzetti, M.; Wenderoth, N.; Mantini, D. Detecting Large-Scale Brain Networks Using EEG: Impact of Electrode Density, Head Modeling and Source Localization. Front. Neuroinform. 2018, 12, 4. [Google Scholar] [CrossRef]
  17. Michel, C.M.; Murray, M.M.; Lantz, G.; Gonzalez, S.; Spinelli, L.; Grave de Peralta, R. EEG source imaging. Clin. Neurophysiol. 2004, 115, 2195–2222. [Google Scholar] [CrossRef] [PubMed]
  18. Bola, M.; Sabel, B.A. Dynamic reorganization of brain functional networks during cognition. Neuroimage 2015, 114, 398–413. [Google Scholar] [CrossRef] [PubMed]
  19. Brodbeck, V.; Spinelli, L.; Lascano, A.M.; Wissmeier, M.; Vargas, M.I.; Vulliemoz, S.; Pollo, C.; Schaller, K.; Michel, C.M.; Seeck, M. Electroencephalographic source imaging: A prospective study of 152 operated epileptic patients. Brain 2011, 134, 2887–2897. [Google Scholar] [CrossRef]
  20. Canuet, L.; Ishii, R.; Pascual-Marqui, R.D.; Iwase, M.; Kurimoto, R.; Aoki, Y.; Ikeda, S.; Takahashi, H.; Nakahachi, T.; Takeda, M. Resting-state EEG source localization and functional connectivity in schizophrenia-like psychosis of epilepsy. PLoS ONE 2011, 6, e27863. [Google Scholar] [CrossRef] [Green Version]
  21. Zhao, M.; Marino, M.; Samogin, J.; Swinnen, S.P.; Mantini, D. Hand, foot and lip representations in primary sensorimotor cortex: A high-density electroencephalography study. Sci. Rep. 2019, 9, 19464. [Google Scholar] [CrossRef]
  22. Li, K.; Papademetris, X.; Tucker, D.M. BrainK for Structural Image Processing: Creating Electrical Models of the Human Head. Comput. Intell. Neurosci. 2016, 2016, 1349851. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Montes-Restrepo, V.; van Mierlo, P.; Strobbe, G.; Staelens, S.; Vandenberghe, S.; Hallez, H. Influence of skull modeling approaches on EEG source localization. Brain Topogr. 2014, 27, 95–111. [Google Scholar] [CrossRef] [Green Version]
  24. Ramon, C.; Haueisen, J.; Schimpf, P.H. Influence of head models on neuromagnetic fields and inverse source localizations. Biomed. Eng. Online 2006, 5, 55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Taberna, G.A.; Samogin, J.; Mantini, D. Automated Head Tissue Modelling Based on Structural Magnetic Resonance Images for Electroencephalographic Source Reconstruction. Neuroinformatics 2021. [Google Scholar] [CrossRef]
  26. Hallez, H.; Vanrumste, B.; Grech, R.; Muscat, J.; De Clercq, W.; Vergult, A.; D’Asseler, Y.; Camilleri, K.P.; Fabri, S.G.; Van Huffel, S.; et al. Review on solving the forward problem in EEG source analysis. J. Neuroeng. Rehabil. 2007, 4, 46. [Google Scholar] [CrossRef]
  27. Akalin Acar, Z.; Makeig, S. Effects of forward model errors on EEG source localization. Brain Topogr. 2013, 26, 378–396. [Google Scholar] [CrossRef] [Green Version]
  28. Cespedes-Villar, Y.; Martinez-Vargas, J.D.; Castellanos-Dominguez, G. Influence of Patient-Specific Head Modeling on EEG Source Imaging. Comput. Math. Methods Med. 2020, 2020, 5076865. [Google Scholar] [CrossRef]
  29. Von Ellenrieder, N.; Muravchik, C.H.; Wagner, M.; Nehorai, A. Effect of head shape variations among individuals on the EEG/MEG forward and inverse problems. IEEE Trans. Biomed. Eng. 2009, 56, 587–597. [Google Scholar] [CrossRef] [PubMed]
  30. Wagner, S.; Rampersad, S.M.; Aydin, U.; Vorwerk, J.; Oostendorp, T.F.; Neuling, T.; Herrmann, C.S.; Stegeman, D.F.; Wolters, C.H. Investigation of tDCS volume conduction effects in a highly realistic head model. J. Neural. Eng. 2014, 11, 016002. [Google Scholar] [CrossRef] [PubMed]
  31. Collins, D.L.; Zijdenbos, A.P.; Kollokian, V.; Sled, J.G.; Kabani, N.J.; Holmes, C.J.; Evans, A.C. Design and construction of a realistic digital brain phantom. IEEE Trans. Med. Imaging 1998, 17, 463–468. [Google Scholar] [CrossRef]
  32. Huang, Y.; Parra, L.C.; Haufe, S. The New York Head-A precise standardized volume conductor model for EEG source localization and tES targeting. Neuroimage 2016, 140, 150–162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Mazziotta, J.; Toga, A.; Evans, A.; Fox, P.; Lancaster, J.; Zilles, K.; Woods, R.; Paus, T.; Simpson, G.; Pike, B.; et al. A probabilistic atlas and reference system for the human brain: International Consortium for Brain Mapping (ICBM). Philos. Trans. R. Soc. Lond. B Biol. Sci. 2001, 356, 1293–1322. [Google Scholar] [CrossRef] [PubMed]
  34. Valdes-Hernandez, P.A.; von Ellenrieder, N.; Ojeda-Gonzalez, A.; Kochen, S.; Aleman-Gomez, Y.; Muravchik, C.; Valdes-Sosa, P.A. Approximate average head models for EEG source imaging. J. Neurosci. Methods 2009, 185, 125–132. [Google Scholar] [CrossRef]
  35. Fischl, B. FreeSurfer. Neuroimage 2012, 62, 774–781. [Google Scholar] [CrossRef] [Green Version]
  36. Holdefer, R.N.; Sadleir, R.; Russell, M.J. Predicted current densities in the brain during transcranial electrical stimulation. Clin. Neurophysiol. 2006, 117, 1388–1397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Homma, S.; Musha, T.; Nakajima, Y.; Okamoto, Y.; Blom, S.; Flink, R.; Hagbarth, K.E.; Mostrom, U. Location of electric current sources in the human brain estimated by the dipole tracing method of the scalp-skull-brain (SSB) head model. Electroencephalogr. Clin. Neurophysiol. 1994, 91, 374–382. [Google Scholar] [CrossRef]
  38. Oostenveld, R.; Fries, P.; Maris, E.; Schoffelen, J.M. FieldTrip: Open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput. Intell. Neurosci. 2011, 2011, 156869. [Google Scholar] [CrossRef]
  39. Le, J.; Lu, M.; Pellouchoud, E.; Gevins, A. A rapid method for determining standard 10/10 electrode positions for high resolution EEG studies. Electroencephalogr. Clin. Neurophysiol. 1998, 106, 554–558. [Google Scholar] [CrossRef]
  40. Whalen, C.; Maclin, E.L.; Fabiani, M.; Gratton, G. Validation of a method for coregistering scalp recording locations with 3D structural MR images. Hum. Brain Mapp. 2008, 29, 1288–1301. [Google Scholar] [CrossRef]
  41. Clausner, T.; Dalal, S.S.; Crespo-Garcia, M. Photogrammetry-Based Head Digitization for Rapid and Accurate Localization of EEG Electrodes and MEG Fiducial Markers Using a Single Digital SLR Camera. Front. Neurosci. 2017, 11, 264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Russell, G.S.; Jeffrey Eriksen, K.; Poolman, P.; Luu, P.; Tucker, D.M. Geodesic photogrammetry for localizing sensor positions in dense-array EEG. Clin. Neurophysiol. 2005, 116, 1130–1140. [Google Scholar] [CrossRef]
  43. Homolle, S.; Oostenveld, R. Using a structured-light 3D scanner to improve EEG source modeling with more accurate electrode positions. J. Neurosci. Methods 2019, 326, 108378. [Google Scholar] [CrossRef] [PubMed]
  44. Taberna, G.A.; Guarnieri, R.; Mantini, D. SPOT3D: Spatial positioning toolbox for head markers using 3D scans. Sci. Rep. 2019, 9, 12813. [Google Scholar] [CrossRef] [Green Version]
  45. Taberna, G.A.; Marino, M.; Ganzetti, M.; Mantini, D. Spatial localization of EEG electrodes using 3D scanning. J. Neural. Eng. 2019, 16, 026020. [Google Scholar] [CrossRef] [PubMed]
  46. Koessler, L.; Cecchin, T.; Caspary, O.; Benhadid, A.; Vespignani, H.; Maillard, L. EEG-MRI co-registration and sensor labeling using a 3D laser scanner. Ann. Biomed. Eng. 2011, 39, 983–995. [Google Scholar] [CrossRef]
  47. Koessler, L.; Maillard, L.; Benhadid, A.; Vignal, J.P.; Braun, M.; Vespignani, H. Spatial localization of EEG electrodes. Neurophysiol. Clin. 2007, 37, 97–102. [Google Scholar] [CrossRef]
  48. Shirazi, S.Y.; Huang, H.J. More Reliable EEG Electrode Digitizing Methods Can Reduce Source Estimation Uncertainty, but Current Methods Already Accurately Identify Brodmann Areas. Front. Neurosci. 2019, 13, 1159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Samogin, J.; Marino, M.; Porcaro, C.; Wenderoth, N.; Dupont, P.; Swinnen, S.P.; Mantini, D. Frequency-dependent functional connectivity in resting state networks. Hum. Brain Mapp. 2020, 41, 5187–5198. [Google Scholar] [CrossRef]
  50. Cuartas Morales, E.; Acosta-Medina, C.D.; Castellanos-Dominguez, G.; Mantini, D. A Finite-Difference Solution for the EEG Forward Problem in Inhomogeneous Anisotropic Media. Brain Topogr. 2019, 32, 229–239. [Google Scholar] [CrossRef]
  51. Mantini, D.; Franciotti, R.; Romani, G.L.; Pizzella, V. Improving MEG source localizations: An automated method for complete artifact removal based on independent component analysis. Neuroimage 2008, 40, 160–173. [Google Scholar] [CrossRef]
  52. Liu, Q.; Balsters, J.H.; Baechinger, M.; van der Groen, O.; Wenderoth, N.; Mantini, D. Estimating a neutral reference for electroencephalographic recordings: The importance of using a high-density montage and a realistic head model. J. Neural. Eng. 2015, 12, 056012. [Google Scholar] [CrossRef] [PubMed]
  53. Yao, D. A method to standardize a reference of scalp EEG recordings to a point at infinity. Physiol. Meas. 2001, 22, 693–711. [Google Scholar] [CrossRef]
  54. Yao, D.; Wang, L.; Oostenveld, R.; Nielsen, K.D.; Arendt-Nielsen, L.; Chen, A.C. A comparative study of different references for EEG spectral mapping: The issue of the neutral reference and the use of the infinity reference. Physiol. Meas. 2005, 26, 173–184. [Google Scholar] [CrossRef]
  55. Haueisen, J.; Ramon, C.; Eiselt, M.; Brauer, H.; Nowak, H. Influence of tissue resistivities on neuromagnetic fields and electric potentials studied with a finite element model of the head. IEEE Trans. Biomed. Eng. 1997, 44, 727–735. [Google Scholar] [CrossRef]
  56. Pascual-Marqui, R.D.; Lehmann, D.; Koukkou, M.; Kochi, K.; Anderer, P.; Saletu, B.; Tanaka, H.; Hirata, K.; John, E.R.; Prichep, L.; et al. Assessing interactions in the brain with exact low-resolution electromagnetic tomography. Philos. Trans. A Math. Phys. Eng. Sci. 2011, 369, 3768–3784. [Google Scholar] [CrossRef] [PubMed]
  57. Deng, Y.; Shi, L.; Lei, Y.; Liang, P.; Li, K.; Chu, W.C.; Wang, D.; Alzheimer’s Disease Neuroimaging, I. Mapping the “What” and “Where” Visual Cortices and Their Atrophy in Alzheimer’s Disease: Combined Activation Likelihood Estimation with Voxel-Based Morphometry. Front. Hum. Neurosci. 2016, 10, 333. [Google Scholar] [CrossRef]
  58. Fishbein, D.H.; Eldreth, D.L.; Hyde, C.; Matochik, J.A.; London, E.D.; Contoreggi, C.; Kurian, V.; Kimes, A.S.; Breeden, A.; Grant, S. Risky decision making and the anterior cingulate cortex in abstinent drug abusers and nonusers. Brain Res. Cogn. Brain Res. 2005, 23, 119–136. [Google Scholar] [CrossRef]
  59. Kaufmann, C.; Wehrle, R.; Wetter, T.C.; Holsboer, F.; Auer, D.P.; Pollmacher, T.; Czisch, M. Brain activation and hypothalamic functional connectivity during human non-rapid eye movement sleep: An EEG/fMRI study. Brain 2006, 129, 655–667. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. van Buuren, M.; Vink, M.; Kahn, R.S. Default-mode network dysfunction and self-referential processing in healthy siblings of schizophrenia patients. Schizophr. Res. 2012, 142, 237–243. [Google Scholar] [CrossRef] [Green Version]
  61. Hipp, J.F.; Hawellek, D.J.; Corbetta, M.; Siegel, M.; Engel, A.K. Large-scale cortical correlation structure of spontaneous oscillatory activity. Nat. Neurosci. 2012, 15, 884–890. [Google Scholar] [CrossRef] [Green Version]
  62. de Pasquale, F.; Della Penna, S.; Snyder, A.Z.; Marzetti, L.; Pizzella, V.; Romani, G.L.; Corbetta, M. A cortical core for dynamic integration of functional networks in the resting human brain. Neuron 2012, 74, 753–764. [Google Scholar] [CrossRef] [Green Version]
  63. King, B.R.; van Ruitenbeek, P.; Leunissen, I.; Cuypers, K.; Heise, K.F.; Santos Monteiro, T.; Hermans, L.; Levin, O.; Albouy, G.; Mantini, D.; et al. Age-Related Declines in Motor Performance are Associated With Decreased Segregation of Large-Scale Resting State Brain Networks. Cereb. Cortex 2018, 28, 4390–4402. [Google Scholar] [CrossRef]
  64. Newton, A.T.; Morgan, V.L.; Rogers, B.P.; Gore, J.C. Modulation of steady state functional connectivity in the default mode and working memory networks by cognitive load. Hum. Brain Mapp. 2011, 32, 1649–1659. [Google Scholar] [CrossRef] [Green Version]
  65. Vema Krishna Murthy, S.; MacLellan, M.; Beyea, S.; Bardouille, T. Faster and improved 3-D head digitization in MEG using Kinect. Front. Neurosci. 2014, 8, 326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Iacono, M.I.; Neufeld, E.; Akinnagbe, E.; Bower, K.; Wolf, J.; Vogiatzis Oikonomidis, I.; Sharma, D.; Lloyd, B.; Wilm, B.J.; Wyss, M.; et al. MIDA: A Multimodal Imaging-Based Detailed Anatomical Model of the Human Head and Neck. PLoS ONE 2015, 10, e0124126. [Google Scholar] [CrossRef]
  67. Searle, S.R.; Speed, F.M.; Milliken, G.A. Population Marginal Means in the Linear Model: An Alternative to Least Squares Means. Am. Stat. 1980, 34, 216–221. [Google Scholar] [CrossRef]
  68. Coito, A.; Michel, C.M.; van Mierlo, P.; Vulliemoz, S.; Plomp, G. Directed Functional Brain Connectivity Based on EEG Source Imaging: Methodology and Application to Temporal Lobe Epilepsy. IEEE Trans. Biomed. Eng. 2016, 63, 2619–2628. [Google Scholar] [CrossRef]
  69. He, B.; Dai, Y.; Astolfi, L.; Babiloni, F.; Yuan, H.; Yang, L. eConnectome: A MATLAB toolbox for mapping and imaging of brain functional connectivity. J. Neurosci. Methods 2011, 195, 261–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Michel, C.M.; Brunet, D. EEG Source Imaging: A Practical Review of the Analysis Steps. Front. Neurol. 2019, 10, 325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Beltrachini, L.; von Ellenrieder, N.; Muravchik, C.H. General bounds for electrode mislocation on the EEG inverse problem. Comput. Methods Programs Biomed. 2011, 103, 1–9. [Google Scholar] [CrossRef] [PubMed]
  72. Wang, Y.; Gotman, J. The influence of electrode location errors on EEG dipole source localization with a realistic head model. Clin. Neurophysiol. 2001, 112, 1777–1780. [Google Scholar] [CrossRef]
  73. Marino, M.; Liu, Q.; Brem, S.; Wenderoth, N.; Mantini, D. Automated detection and labeling of high-density EEG electrodes from structural MR images. J. Neural. Eng. 2016, 13, 056003. [Google Scholar] [CrossRef]
  74. De Munck, J.C.; Goncalves, S.I.; Huijboom, L.; Kuijer, J.P.; Pouwels, P.J.; Heethaar, R.M.; Lopes da Silva, F.H. The hemodynamic response of the alpha rhythm: An EEG/fMRI study. Neuroimage 2007, 35, 1142–1151. [Google Scholar] [CrossRef]
  75. Scarff, C.J.; Reynolds, A.; Goodyear, B.G.; Ponton, C.W.; Dort, J.C.; Eggermont, J.J. Simultaneous 3-T fMRI and high-density recording of human auditory evoked potentials. Neuroimage 2004, 23, 1129–1142. [Google Scholar] [CrossRef] [PubMed]
Brainsci 11 00741 g001 550
Figure 1. Anatomical positions of the 21 seeds used for analysis, subdivided into the corresponding networks. Default mode network (DMN): posterior cingulate cortex (PCC), medial prefrontal cortex (MPFC), left/right angular gyrus (lANG/rANG); Dorsal attention network (DAN): left/right frontal eye field (lFEF/rFEF), left/right inferior parietal sulcus (lIPS/rIPS); Ventral attention network (VAN): right temporo-parietal junction (rTPJ), right inferior frontal gyrus (rIFG); Language network (LN, green): left temporo-parietal junction (lTPJ), left inferior frontal gyrus (lIFG); Somatomotor network (SMN): supplementary motor area (SMA), left/right primary somatosensory cortex (lS1/rS1), left/right secondary somatosensory cortex (lS2/rS2); Visual network (VN): left/right human ventral visual area 4 (lV4v/rV4v), left/right visual area 5 (lV5/rV5).
Figure 1. Anatomical positions of the 21 seeds used for analysis, subdivided into the corresponding networks. Default mode network (DMN): posterior cingulate cortex (PCC), medial prefrontal cortex (MPFC), left/right angular gyrus (lANG/rANG); Dorsal attention network (DAN): left/right frontal eye field (lFEF/rFEF), left/right inferior parietal sulcus (lIPS/rIPS); Ventral attention network (VAN): right temporo-parietal junction (rTPJ), right inferior frontal gyrus (rIFG); Language network (LN, green): left temporo-parietal junction (lTPJ), left inferior frontal gyrus (lIFG); Somatomotor network (SMN): supplementary motor area (SMA), left/right primary somatosensory cortex (lS1/rS1), left/right secondary somatosensory cortex (lS2/rS2); Visual network (VN): left/right human ventral visual area 4 (lV4v/rV4v), left/right visual area 5 (lV5/rV5).
Brainsci 11 00741 g001
Brainsci 11 00741 g002 550
Figure 2. Example of EEG electrode positions and their simulated misplacement. From top to bottom: original electrode positions for a representative participant; positions with a systematic error of 1 cm; positions with a random error of 1 cm.
Figure 2. Example of EEG electrode positions and their simulated misplacement. From top to bottom: original electrode positions for a representative participant; positions with a systematic error of 1 cm; positions with a random error of 1 cm.
Brainsci 11 00741 g002
Brainsci 11 00741 g003 550
Figure 3. Example of head tissue segmentations used in the study. From left to right: T1-weighted individual MR image; 12-layer segmentation using the MR-TIM toolbox; 12-layer segmentation using WTS; 3-layer segmentation using WTS.
Figure 3. Example of head tissue segmentations used in the study. From left to right: T1-weighted individual MR image; 12-layer segmentation using the MR-TIM toolbox; 12-layer segmentation using WTS; 3-layer segmentation using WTS.
Brainsci 11 00741 g003
Brainsci 11 00741 g004 550
Figure 4. Functional connectivity matrices for the following frequency bands: delta (1–4 Hz), theta (4–8 Hz), alpha (8–13 Hz), beta (13–30 Hz) and gamma (30–80 Hz). In each matrix, the diagonal values represent the group average intra-network connectivity, whereas the off-diagonal values represent the group average inter-network connectivity.
Figure 4. Functional connectivity matrices for the following frequency bands: delta (1–4 Hz), theta (4–8 Hz), alpha (8–13 Hz), beta (13–30 Hz) and gamma (30–80 Hz). In each matrix, the diagonal values represent the group average intra-network connectivity, whereas the off-diagonal values represent the group average inter-network connectivity.
Brainsci 11 00741 g004
Brainsci 11 00741 g005 550
Figure 5. Intra- and inter-network connectivity differences per network and frequency band. A comparison was performed by means of a Wilcoxon signed-rank paired test; the color bar represents the z-value obtained through the test, where the red colors stand for intra-network connectivity values higher than inter-network ones, and the blue colors vice-versa.
Figure 5. Intra- and inter-network connectivity differences per network and frequency band. A comparison was performed by means of a Wilcoxon signed-rank paired test; the color bar represents the z-value obtained through the test, where the red colors stand for intra-network connectivity values higher than inter-network ones, and the blue colors vice-versa.
Brainsci 11 00741 g005
Brainsci 11 00741 g006 550
Figure 6. Correlation between the reference connectivity matrices and those obtained when introducing electrode localization errors. For each error, magnitude is tested in the range 0.25 to 1 cm, and distinction is made between systematic (dark grey) and random (light grey) errors. Black dots represent the outliers. The asterisks define the significance level of the two-tailed Wilcoxon signed-rank paired test between error types, for each error magnitude: * for p < 0.05, ** for p < 0.01, **** for p < 0.0001.
Figure 6. Correlation between the reference connectivity matrices and those obtained when introducing electrode localization errors. For each error, magnitude is tested in the range 0.25 to 1 cm, and distinction is made between systematic (dark grey) and random (light grey) errors. Black dots represent the outliers. The asterisks define the significance level of the two-tailed Wilcoxon signed-rank paired test between error types, for each error magnitude: * for p < 0.05, ** for p < 0.01, **** for p < 0.0001.
Brainsci 11 00741 g006
Brainsci 11 00741 g007 550
Figure 7. Correlation between the reference connectivity matrices and those obtained when introducing errors in the head model. These errors were related to the electrode localization technique (dark grey) or to the head tissue segmentation method (light grey). Black dots represent the outliers.
Figure 7. Correlation between the reference connectivity matrices and those obtained when introducing errors in the head model. These errors were related to the electrode localization technique (dark grey) or to the head tissue segmentation method (light grey). Black dots represent the outliers.
Brainsci 11 00741 g007
Brainsci 11 00741 g008 550
Figure 8. Impact of head modeling on the intra- and inter-network connectivity differences, for each RSN. This was assessed by calculating the correlation with values obtained using the reference analysis workflow. Two electrode localization techniques (3D Scan and Digitizer—on the left) and three head tissue segmentation methods (12-layer WTS, 3-layer WTS and 3-layer template—on the right) were tested.
Figure 8. Impact of head modeling on the intra- and inter-network connectivity differences, for each RSN. This was assessed by calculating the correlation with values obtained using the reference analysis workflow. Two electrode localization techniques (3D Scan and Digitizer—on the left) and three head tissue segmentation methods (12-layer WTS, 3-layer WTS and 3-layer template—on the right) were tested.
Brainsci 11 00741 g008
Table
Table 1. ANOVA testing for electrode localization errors. Four factors were assessed: resting-state network (DMN, DAN, VAN, LN, SMN and VN), frequency band (delta; theta; alpha; beta; gamma), error magnitude (0.25; 0.5; 0.75; 1 cm) and error type (systematic; random).
Table 1. ANOVA testing for electrode localization errors. Four factors were assessed: resting-state network (DMN, DAN, VAN, LN, SMN and VN), frequency band (delta; theta; alpha; beta; gamma), error magnitude (0.25; 0.5; 0.75; 1 cm) and error type (systematic; random).
Electrode Localization ErrorsdfSum SquaresMean SquareFp-Value
Network56.581.3279.39<0.001
Band418.104.53273.19<0.001
Error magnitude37.282.43146.54<0.001
Error type10.460.4627.72<0.001
Residuals2263.740.02
Total23936.16
Table
Table 2. ANOVA testing for head tissue segmentation methods. Three factors were assessed: resting-state network (DMN; DAN; VAN; LN; SMN; VN), frequency band (delta; theta; alpha; beta; gamma) and segmentation method (12-layer WTS; 3-layer WTS; 3-layer template).
Table 2. ANOVA testing for head tissue segmentation methods. Three factors were assessed: resting-state network (DMN; DAN; VAN; LN; SMN; VN), frequency band (delta; theta; alpha; beta; gamma) and segmentation method (12-layer WTS; 3-layer WTS; 3-layer template).
Head Tissue Segmentation MethodsdfSum SquaresMean SquareFp-Value
Network51.180.2419.52<0.001
Band45.831.46120.42<0.001
Segmentation method20.290.1512.17<0.001
Residuals780.940.01
Total898.24
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Taberna, G.A.; Samogin, J.; Marino, M.; Mantini, D. Detection of Resting-State Functional Connectivity from High-Density Electroencephalography Data: Impact of Head Modeling Strategies. Brain Sci. 2021, 11, 741. https://doi.org/10.3390/brainsci11060741

AMA Style

Taberna GA, Samogin J, Marino M, Mantini D. Detection of Resting-State Functional Connectivity from High-Density Electroencephalography Data: Impact of Head Modeling Strategies. Brain Sciences. 2021; 11(6):741. https://doi.org/10.3390/brainsci11060741

Chicago/Turabian Style

Taberna, Gaia Amaranta, Jessica Samogin, Marco Marino, and Dante Mantini. 2021. "Detection of Resting-State Functional Connectivity from High-Density Electroencephalography Data: Impact of Head Modeling Strategies" Brain Sciences 11, no. 6: 741. https://doi.org/10.3390/brainsci11060741

APA Style

Taberna, G. A., Samogin, J., Marino, M., & Mantini, D. (2021). Detection of Resting-State Functional Connectivity from High-Density Electroencephalography Data: Impact of Head Modeling Strategies. Brain Sciences, 11(6), 741. https://doi.org/10.3390/brainsci11060741

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

Article Metrics

Back to TopTop