Functional Connectivity Estimates Research Articles

Online functional connectivity analysis of large all-to-all networks in MNE Scan

Abstract The analysis of electroencephalography (EEG)/magnetoencephalography (MEG) functional connectivity has become an important tool in neuroscience. Especially the high time resolution of EEG/MEG enables important insight into the functioning of the human brain. To date, functional connectivity is commonly estimated offline, that is, after the conclusion of the experiment. However, online computation of functional connectivity has the potential to enable unique experimental paradigms. For example, changes of functional connectivity due to learning processes could be tracked in real time and the experiment be adjusted based on these observations. Furthermore, the connectivity estimates can be used for neurofeedback applications or the instantaneous inspection of measurement results. In this study, we present the implementation and evaluation of online sensor and source space functional connectivity estimation in the open-source software MNE Scan. Online capable implementations of several functional connectivity metrics were established in the Connectivity library within MNE-CPP and made available as a plugin in MNE Scan. Online capability was achieved by enforcing multithreading and high efficiency for all computations, so that repeated computations were avoided wherever possible, which allows for a major speed-up in the case of overlapping intervals. We present comprehensive performance evaluations of these implementations proving the online capability for the computation of large all-to-all functional connectivity networks. As a proof of principle, we demonstrate the feasibility of online functional connectivity estimation in the evaluation of somatosensory evoked brain activity

Investigating Functional Connectivity in Adolescent Depression and Suicide Attempt during Neurofeedback Sessions: A Multivariate Random Covariance Model Approach

Background A recent neurofeedback functional magnetic resonance imaging (NF, fMRI) study on depressed vs. healthy adolescents elicited differential functional connectivity (FC) amongst brain regions of interest (ROIs). Previous results employed univariate methods and included only two seed areas of FC (amygdala and hippocampus). In this study, we propose a new multivariate analysis for whole-network FC estimation. Methods Primary analyses concerned a pre-identified network of 17 salient ROIs reflecting key regions in self-processing and emotion regulation. A random covariance model (RCM) was applied to jointly estimate participant- and group-specific connectivity, where FC was measured by partial correlation conditioned on or adjusted for rest-of-network connectivity patterns. Secondary analyses concerned participant-specific network association with mental functioning changes and the AAL3 whole-brain atlas. Results New findings suggested that depressed adolescents with a suicide attempt expressed significantly higher positive FC between the left temporal gyrus and the left amygdala during NF, compared to negative FC in non-attempting depressed youth, while healthy controls displayed negative FC between the insula, inferior frontal gyrus to inferior parietal lobe connection, compared to mild negative connectivity in depressed adolescents. Previous cross-hemispheric findings in depressed vs. healthy adolescents were corroborated. Conclusion A multivariate RCM uncovered key ROI-pairwise connections differentiating FC patterns between depressed youth vs. healthy controls and among depressed youth, with and without a suicide attempt. Findings were strengthened by enhanced inference vs. univariate methods, and corroboration of previous NF secondary analyses demonstrated future utility for participant-specific study in association with clinical outcomes and/or whole-brain analyses with larger sample sizes.

Dynamicity of brain network organization & their community architecture as characterizing features for classification of common mental disorders from whole-brain connectome

The urgency of addressing common mental disorders (bipolar disorder, attention-deficit hyperactivity disorder (ADHD), and schizophrenia) arises from their significant societal impact. Developing strategies to support psychiatrists is crucial. Previous studies focused on the relationship between these disorders and changes in the resting-state functional connectome’s modularity, often using static functional connectivity (sFC) estimation. However, understanding the dynamic reconfiguration of resting-state brain networks with rich temporal structure is essential for comprehending neural activity and addressing mental health disorders. This study proposes an unsupervised approach combining spatial and temporal characterization of brain networks to classify common mental disorders using fMRI timeseries data from two cohorts (N = 408 participants). We employ the weighted stochastic block model to uncover mesoscale community architecture differences, providing insights into network organization. Our approach overcomes sFC limitations and biases in community detection algorithms by modelling the functional connectome’s temporal dynamics as a landscape, quantifying temporal stability at whole-brain and network levels. Findings reveal individuals with schizophrenia exhibit less assortative community structure and participate in multiple motif classes, indicating less specialized network organization. Patients with schizophrenia and ADHD demonstrate significantly reduced temporal stability compared to healthy controls. This study offers insights into functional connectivity (FC) patterns’ spatiotemporal organization and their alterations in common mental disorders, highlighting the potential of temporal stability as a biomarker.

Estimation of static and dynamic functional connectivity in resting-state fMRI using zero-frequency resonator.

Resting-state functional magnetic resonance imaging (rs-fMRI) is increasingly being used to infer the functional organization of the brain. Blood oxygen level-dependent (BOLD) features related to spontaneous neuronal activity, are yet to be clearly understood. Prior studies have hypothesized that rs-fMRI is spontaneous event-related and these events convey crucial information about the neuronal activity in estimating resting state functional connectivity (FC). Attempts have been made to extract these temporal events using a predetermined threshold. However, the thresholding methods in addition to being very sensitive to noise, may consider redundant events or exclude the low-valued inflection points. Here, we extract the event-related temporal onsets from the rs-fMRI time courses using a zero-frequency resonator (ZFR). The ZFR reflects the transient behavior of the BOLD events at its output. The conditional rate (CR) of the BOLD events occurring in a time course with respect to a seed time course is used to derive static FC. The temporal activity around the estimated events called high signal-to-noise ratio (SNR) segments are also obtained in the rs-fMRI time course and are then used to compute static and dynamic FCs during rest. Coactivation pattern (CAP) is the dynamic FC obtained using the high SNR segments driven by the ZFR. The static FC demonstrates that the ZFR-based CR distinguishes the coactivation and non-coactivation scores well in the distribution. CAP analysis demonstrated the stable and longer dwell time dominant resting state functional networks with high SNR segments driven by the ZFR. Static and dynamic FC analysis underpins that the ZFR-driven temporal onsets of BOLD events derive reliable and consistent FCs in the resting brain using a subset of the time points.

In-Scanner Thoughts shape Resting-state Functional Connectivity: how participants "rest" matters.

Resting-state fMRI (rs-fMRI) scans-namely those lacking experimentally-controlled stimuli or cognitive demands-are often used to identify aberrant patterns of functional connectivity (FC) in clinical populations. To minimize interpretational uncertainty, researchers control for across-cohort disparities in age, gender, co-morbidities, and head motion. Yet, studies rarely, if ever, consider the possibility that systematic differences in inner experience (i.e., what subjects think and feel during the scan) may directly affect FC measures. Here we demonstrate that is the case using a rs-fMRI dataset comprising 471 scans annotated with experiential data. Wide-spread significant differences in FC are observed between scans that systematically differ in terms of reported in-scanner experience. Additionally, we show that FC can successfully predict specific aspects of in-scanner experience in a manner similar to how it predicts demographics, cognitive abilities, clinical outcomes and labels. Together, these results highlight the key role of in-scanner experience in shaping rs-fMRI estimates of FC.

The effects of spatial leakage correction on the reliability of EEG-based functional connectivity networks.

Electroencephalography (EEG) functional connectivity (FC) estimates are confounded by the volume conduction problem. This effect can be greatly reduced by applying FC measures insensitive to instantaneous, zero-lag dependencies (corrected measures). However, numerous studies showed that FC measures sensitive to volume conduction (uncorrected measures) exhibit higher reliability and higher subject-level identifiability. We tested how source reconstruction contributed to the reliability difference of EEG FC measures on a large (n = 201) resting-state data set testing eight FC measures (including corrected and uncorrected measures). We showed that the high reliability of uncorrected FC measures in resting state partly stems from source reconstruction: idiosyncratic noise patterns define a baseline resting-state functional network that explains a significant portion of the reliability of uncorrected FC measures. This effect remained valid for template head model-based, as well as individual head model-based source reconstruction. Based on our findings we made suggestions how to best use spatial leakage corrected and uncorrected FC measures depending on the main goals of the study.

Multiscale Modes of Functional Brain Connectivity.

Information processing in the brain spans from localised sensorimotor processes to higher-level cognition that integrates across multiple regions. Interactions between and within these subsystems enable multiscale information processing. Despite this multiscale characteristic, functional brain connectivity is often either estimated based on 10-30 distributed modes or parcellations with 100-1000 localised parcels, both missing across-scale functional interactions. We present Multiscale Probabilistic Functional Modes (mPFMs), a new mapping which comprises modes over various scales of granularity, thus enabling direct estimation of functional connectivity within- and across-scales. Crucially, mPFMs emerged from data-driven multilevel Bayesian modelling of large functional MRI (fMRI) populations. We demonstrate that mPFMs capture both distributed brain modes and their co-existing subcomponents. In addition to validating mPFMs using simulations and real data, we show that mPFMs can predict ~900 personalised traits from UK Biobank more accurately than current standard techniques. Therefore, mPFMs can offer a paradigm shift in functional connectivity modelling and yield enhanced fMRI biomarkers for traits and diseases.

EEGminer: discovering interpretable features of brain activity with learnable filters

Objective. The patterns of brain activity associated with different brain processes can be used to identify different brain states and make behavioural predictions. However, the relevant features are not readily apparent and accessible. Our aim is to design a system for learning informative latent representations from multichannel recordings of ongoing EEG activity. Approach: We propose a novel differentiable decoding pipeline consisting of learnable filters and a pre-determined feature extraction module. Specifically, we introduce filters parameterized by generalized Gaussian functions that offer a smooth derivative for stable end-to-end model training and allow for learning interpretable features. For the feature module, we use signal magnitude and functional connectivity estimates. Main results. We demonstrate the utility of our model on a new EEG dataset of unprecedented size (i.e. 721 subjects), where we identify consistent trends of music perception and related individual differences. Furthermore, we train and apply our model in two additional datasets, specifically for emotion recognition on SEED and workload classification on simultaneous task EEG workload. The discovered features align well with previous neuroscience studies and offer new insights, such as marked differences in the functional connectivity profile between left and right temporal areas during music listening. This agrees with the specialisation of the temporal lobes regarding music perception proposed in the literature. Significance. The proposed method offers strong interpretability of learned features while reaching similar levels of accuracy achieved by black box deep learning models. This improved trustworthiness may promote the use of deep learning models in real world applications. The model code is available at https://github.com/SMLudwig/EEGminer/.

Neuroimaging alterations and relapse in early-stage psychosis.

Recent reports have indicated that symptom exacerbation after a period of improvement, referred to as relapse, in early-stage psychosis could result in brain changes and poor disease outcomes. We hypothesized that substantial neuroimaging alterations may exist among patients who experience relapse in early-stage psychosis. We studied patients with psychosis within 2 years after the first psychotic event and healthy controls. We divided patients into 2 groups, namely those who did not experience relapse between disease onset and the magnetic resonance imaging (MRI) scan (no-relapse group) and those who did experience relapse between these 2 timings (relapse group). We analyzed 3003 functional connectivity estimates between 78 regions of interest (ROIs) derived from resting-state functional MRI data by adjusting for demographic and clinical confounding factors. We studied 85 patients, incuding 54 in the relapse group and 31 in the no-relapse group, along with 94 healthy controls. We observed significant differences in 47 functional connectivity estimates between the relapse and control groups after multiple comparison corrections, whereas no differences were found between the no-relapse and control groups. Most of these pathological signatures (64%) involved the thalamus. The Jonckheere-Terpstra test indicated that all 47 functional connectivity changes had a significant cross-group progression from controls to patients in the no-relapse group to patients in the relapse group. Longitudinal studies are needed to further validate the involvement and pathological importance of the thalamus in relapse. We observed pathological differences in neuronal connectivity associated with relapse in early-stage psychosis, which are more specifically associated with the thalamus. Our study implies the importance of considering neurobiological mechanisms associated with relapse in the trajectory of psychotic disorders.

Static and dynamic fMRI-derived functional connectomes represent largely similar information.

Functional connectivity (FC) of blood oxygen level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate), and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

Toward stability of dynamic FC estimates in neuroimaging and electrophysiology: Solutions and limits.

Time-varying functional connectivity (FC) methods are used to map the spatiotemporal organization of brain activity. However, their estimation can be unstable, in the sense that different runs of the inference may yield different solutions. But to draw meaningful relations to behavior, estimates must be robust and reproducible. Here, we propose two solutions using the hidden Markov model (HMM) as a descriptive model of time-varying FC. The first, best ranked HMM, involves running the inference multiple times and selecting the best model based on a quantitative measure combining fitness and model complexity. The second, hierarchical-clustered HMM, generates stable cluster state time series by applying hierarchical clustering to the state time series obtained from multiple runs. Experimental results on fMRI and magnetoencephalography data demonstrate that these approaches substantially improve the stability of time-varying FC estimations. Overall, hierarchical-clustered HMM is preferred when the inference variability is high, while the best ranked HMM performs better otherwise.

MEG Dynamic functional connectivity as a longitudinal biomarker for Alzheimer’s Disease relatives

AbstractBackgroundAlzheimer’s disease (AD) is the most common cause of dementia. Although it presents a long preclinical period, treatment has yet to be found, and individuals at risk, like first‐degree relatives, are being studied to identify early detection markers to maximise the effectiveness of treatments. Functional connectivity changes have been observed even at this early stage (Ramirez‐Toraño et al., 2021), but most disregard its temporal dimension. This study investigates the value of a simple estimation of dynamic functional connectivity (dFC) in a longitudinal follow‐up of individuals with varying risk of the disease, as this parameter is expected to decrease across the disease continuum (Grieder et al, 2018).Method95 cognitively unimpaired participants with a family history of AD (FH+) and 39 without it (FH‐) underwent two longitudinal 5‐minute resting‐state magnetoencephalography recordings, approximately 3 years apart. FC in the low alpha band (8‐10 Hz) was estimated in 4‐second segments using phase locking value in the source space, and the standard deviation (STD) across epochs was calculated. Then, a STD ratio was derived dividing the STD values of the follow‐up session by those of the baseline session, and then averaged over nodes belonging to the same regions of the automatic anatomical labelling atlas. Left‐tailed t‐tests between the STD ratio of both groups in each ROI link (n_comparisons = 3160) were calculated and corrected for multiple comparisons with a false discovery rate using q = 0.20.Result13 links show a significantly lower dFC in the FH+ group than in the FH‐ group (Figure 1). Interestingly, these links involve regions where amyloid‐beta deposits early in the AD continuum (Mattsson et al., 2019), and are strongly related to the default mode network (Palmqvist et al. 2017).ConclusionThis study suggests that a simple estimation of dFC like STD may be useful in early detection pathological patterns characteristic of AD, even in cognitively unimpaired individuals at risk. The fact that dFC decreases were detected in areas where amyloid deposits appear at an early stage might indicate that these changes are related to changes in the excitation/inhibition ratio, which will be further studied in upcoming investigations.

Multi-channel EEG-based sleep staging using brain functional connectivity and domain adaptation

Objective. Sleep stage recognition has essential clinical value for evaluating human physical/mental condition and diagnosing sleep-related diseases. To conduct a five-class (wake, N1, N2, N3 and rapid eye movement) sleep staging task, twenty subjects with recorded six-channel electroencephalography (EEG) signals from the ISRUC-SLEEP dataset is used. Approach. Unlike the exist methods ignoring the channel coupling relationship and non-stationarity characteristics, we developed a brain functional connectivity method to provide a new insight for multi-channel analysis. Furthermore, we investigated three frequency-domain features: two functional connectivity estimations, i.e. synchronization likelihood (SL) and wavelet-based correlation (WC) among four frequency bands, and energy ratio (ER) related to six frequency bands, respectively. Then, the Gaussian support vector machine (SVM) method was used to predict the five sleep stages. The performance of the applied features is evaluated in both subject dependence experiment by ten-fold cross validation and subject independence experiment by leave-one-subject-out cross-validation, respectively. Main results. In subject dependence experiment, the results showed that the fused feature (fusion of SL, WC and ER features) contributes significant gain the performance of SVM classifier, where the mean of classification accuracy can achieve 83.97% ± 1.04%. However, in subject-independence experiment, the individual differences EEG patterns across subjects leads to inferior accuracy. Five typical domain adaptation (DA) methods were applied to reduce the discrepancy of feature distributions by selecting the optimal subspace dimension. Results showed that four DA methods can significantly improve the mean accuracy by 1.89%–5.22% compared to the baseline accuracy 57.44% in leave-one-subject-out cross-validation. Significance. Compared with traditional time-frequency and nonlinear features, brain functional connectivity features can capture the correlation between different brain regions. For the individual EEG response differences, domain adaptation methods can transform features to improve the performance of sleep staging algorithms.

Similarity between structural and proxy estimates of brain connectivity

Functional magnetic resonance and diffusion weighted imaging have so far made a major contribution to delineation of the brain connectome at the macroscale. While functional connectivity (FC) was shown to be related to structural connectivity (SC) to a certain degree, their spatial overlap is unknown. Even less clear are relations of SC with estimates of connectivity from inter-subject covariance of regional F18-fluorodeoxyglucose uptake (FDGcov) and grey matter volume (GMVcov). Here, we asked to what extent SC underlies three proxy estimates of brain connectivity: FC, FDGcov and GMVcov. Simultaneous PET/MR acquisitions were performed in 56 healthy middle-aged individuals. Similarity between four networks was assessed using Spearman correlation and convergence ratio (CR), a measure of spatial overlap. Spearman correlation coefficient was 0.27 for SC-FC, 0.40 for SC-FDGcov, and 0.15 for SC-GMVcov. Mean CRs were 51% for SC-FC, 48% for SC-FDGcov, and 37% for SC-GMVcov. These results proved to be reproducible and robust against image processing steps. In sum, we found a relevant similarity of SC with FC and FDGcov, while GMVcov consistently showed the weakest similarity. These findings indicate that white matter tracts underlie FDGcov to a similar degree as FC, supporting FDGcov as estimate of functional brain connectivity.

Impact of sample size and regression of tissue-specific signals on effective connectivity within the core default mode network.

Interactions within brain networks are inherently directional, which are inaccessible to classical functional connectivity estimates from resting-state functional magnetic resonance imaging (fMRI) but can be detected using spectral dynamic causal modeling (DCM). The sample size and unavoidable presence of nuisance signals during fMRI measurement are the two important factors influencing the stability of group estimates of connectivity parameters. However, most recent studies exploring effective connectivity (EC) have been conducted with small sample sizes and minimally pre-processed datasets. We explore the impact of these two factors by analyzing clean resting-state fMRI data from 330 unrelated subjects from the Human Connectome Project database. We demonstrate that both the stability of the model selection procedures and the inference of connectivity parameters are highly dependent on the sample size. The minimum sample size required for stable DCM is approximately 50, which may explain the variability of the DCM results reported so far. We reveal a stable pattern of EC within the core default mode network computed for large sample sizes and demonstrate that the use of subject-specific thresholded whole-brain masks for tissue-specific signals regression enhances the detection of weak connections.

PSIICOS projection optimality for EEG and MEG based functional coupling detection

Functional connectivity is crucial for cognitive processes in the healthy brain and serves as a marker for a range of neuropathological conditions. Non-invasive exploration of functional coupling using temporally resolved techniques such as MEG allows for a unique opportunity of exploring this fundamental brain mechanism.The indirect nature of MEG measurements complicates the estimation of functional coupling due to the volume conduction and spatial leakage effects. In the previous work (Ossadtchi et al., 2018), we introduced PSIICOS, a method that for the first time allowed us to suppress the volume conduction effect and yet retain information about functional networks whose nodes are coupled with close to zero or zero mutual phase lag.In this paper, we demonstrate analytically that the PSIICOS projection is optimal in achieving a controllable trade-off between suppressing mutual spatial leakage and retaining information about zero- or close to zero-phase coupled networks. We also derive an alternative solution using the regularization-based inverse of the mutual spatial leakage matrix and show its equivalence to the original PSIICOS.We then discuss how PSIICOS solution to the functional connectivity estimation problem can be incorporated into the conventional source estimation framework. Instead of sources, the unknowns are the elementary dyadic networks and their activation time series are formalized by the corresponding source-space cross-spectral coefficients. This view on connectivity estimation as a regression problem opens up new opportunities for formulating a set of principled estimators based on the rich intuition accumulated in the neuroimaging community.

Kernel-Based Regularized EEGNet Using Centered Alignment and Gaussian Connectivity for Motor Imagery Discrimination

Brain–computer interfaces (BCIs) from electroencephalography (EEG) provide a practical approach to support human–technology interaction. In particular, motor imagery (MI) is a widely used BCI paradigm that guides the mental trial of motor tasks without physical movement. Here, we present a deep learning methodology, named kernel-based regularized EEGNet (KREEGNet), leveled on centered kernel alignment and Gaussian functional connectivity, explicitly designed for EEG-based MI classification. The approach proactively tackles the challenge of intrasubject variability brought on by noisy EEG records and the lack of spatial interpretability within end-to-end frameworks applied for MI classification. KREEGNet is a refinement of the widely accepted EEGNet architecture, featuring an additional kernel-based layer for regularized Gaussian functional connectivity estimation based on CKA. The superiority of KREEGNet is evidenced by our experimental results from binary and multiclass MI classification databases, outperforming the baseline EEGNet and other state-of-the-art methods. Further exploration of our model’s interpretability is conducted at individual and group levels, utilizing classification performance measures and pruned functional connectivities. Our approach is a suitable alternative for interpretable end-to-end EEG-BCI based on deep learning.

The confound of hemodynamic response function variability in human resting-state functional MRI studies.

Functional magnetic resonance imaging (fMRI) is an indirect measure of neural activity with the hemodynamic response function (HRF) coupling it with unmeasured neural activity. The HRF, modulated by several non-neural factors, is variable across brain regions, individuals and populations. Yet, a majority of human resting-state fMRI connectivity studies continue to assume a non-variable HRF. In this article, with supportive prior evidence, we argue that HRF variability cannot be ignored as it substantially confounds within-subject connectivity estimates and between-subjects connectivity group differences. We also discuss its clinical relevance with connectivity impairments confounded by HRF aberrations in several disorders. We present limited data on HRF differences between women and men, which resulted in a 15.4% median error in functional connectivity estimates in a group-level comparison. We also discuss the implications of HRF variability for fMRI studies in the spinal cord. There is a need for more dialogue within the community on the HRF confound, and we hope that our article is a catalyst in the process.

Cardiovascular and metabolic health is associated with functional brain connectivity in middle-aged and older adults: Results from the Human Connectome Project-Aging study

Several cardiovascular and metabolic indicators, such as cholesterol and blood pressure have been associated with altered neural and cognitive health as well as increased risk of dementia and Alzheimer's disease in later life. In this cross-sectional study, we examined how an aggregate index of cardiovascular and metabolic risk factor measures was associated with correlation-based estimates of resting-state functional connectivity (FC) across a broad adult age-span (36–90+ years) from 930 volunteers in the Human Connectome Project Aging (HCP-A). Increased (i.e., worse) aggregate cardiometabolic scores were associated with reduced FC globally, with especially strong effects in insular, medial frontal, medial parietal, and superior temporal regions. Additionally, at the network-level, FC between core brain networks, such as default-mode and cingulo-opercular, as well as dorsal attention networks, showed strong effects of cardiometabolic risk. These findings highlight the lifespan impact of cardiovascular and metabolic health on whole-brain functional integrity and how these conditions may disrupt higher-order network integrity.

Reliable evaluation of functional connectivity and graph theory measures in source-level EEG: How many electrodes are enough?

ObjectiveUsing EEG to characterise functional brain networks through graph theory has gained significant interest in clinical and basic research. However, the minimal requirements for reliable measures remain largely unaddressed. Here, we examined functional connectivity estimates and graph theory metrics obtained from EEG with varying electrode densities. MethodsEEG was recorded with 128 electrodes in 33 participants. The high-density EEG data were subsequently subsampled into three sparser montages (64, 32, and 19 electrodes). Four inverse solutions, four measures of functional connectivity, and five graph theory metrics were tested. ResultsThe correlation between the results obtained with 128-electrode and the subsampled montages decreased as a function of the number of electrodes. As a result of decreased electrode density, the network metrics became skewed: mean network strength and clustering coefficient were overestimated, while characteristic path length was underestimated. ConclusionsSeveral graph theory metrics were altered when electrode density was reduced. Our results suggest that, for optimal balance between resource demand and result precision, a minimum of 64 electrodes should be utilised when graph theory metrics are used to characterise functional brain networks in source-reconstructed EEG data. SignificanceCharacterisation of functional brain networks derived from low-density EEG warrants careful consideration.