A systematic evaluation of source reconstruction of resting MEG of the human brain with a new high-resolution atlas: Performance, precision, and parcellation

Resting state electroencephalographic alpha rhythms are sensitive to Alzheimer's disease mild cognitive impairment progression at a 6-month follow-up

Are posterior resting-state electroencephalographic (rsEEG) alpha rhythms sensitive to the Alzheimer's disease mild cognitive impairment (ADMCI) progression at a 6-month follow-up? Clinical, cerebrospinal, neuroimaging, and rsEEG datasets in 52 ADMCI and 60 Healthy old seniors (equivalent groups for demographic features) were available from an international archive (www.pdwaves.eu). The ADMCI patients were arbitrarily divided into two groups: REACTIVE and UNREACTIVE, based on the reduction (reactivity) in the posterior rsEEG alpha eLORETA source activities from the eyes-closed to eyes-open condition at ≥ -10% and -10%, respectively. 75% of the ADMCI patients were REACTIVE. Compared to the UNREACTIVE group, the REACTIVE group showed (1) less abnormal posterior rsEEG source activity during the eyes-closed condition and (2) a decrease in that activity at the 6-month follow-up. These effects could not be explained by neuroimaging and neuropsychological biomarkers of AD. Such a biomarker might reflect abnormalities in cortical arousal in quiet wakefulness to be used for clinical studies in ADMCI patients using 6-month follow-ups.

Dorsal brain activity reflects the severity of menopausal symptoms

OBJECTIVE

The severity of menopausal symptoms, despite being triggered by hormonal imbalance, does not directly correspond to hormone levels in the blood; thus, the level of unpleasantness is assessed using subjective questionnaires in clinical practice. To provide better treatments, alternative objective assessments have been anticipated to support medical interviews and subjective assessments. This study aimed to develop a new objective measurement for assessing unpleasantness.

METHODS

Fourteen participants with menopausal symptoms and two age-matched participants who visited our outpatient section were enrolled. Resting-state brain activity was measured using magnetoencephalography. The level of unpleasantness of menopausal symptoms was measured using the Kupperman Kohnenki Shogai Index. The blood level of follicle-stimulating hormone and luteinizing hormone were also measured. Correlation analyses were performed between the oscillatory power of brain activity, index score, and hormone levels.

RESULTS

The level of unpleasantness of menopausal symptoms was positively correlated with high-frequency oscillatory powers in the parietal and bordering cortices (alpha; P = 0.016, beta; P = 0.015, low gamma; P = 0.010). The follicle-stimulating hormone blood level was correlated with high-frequency oscillatory powers in the dorsal part of the cortex (beta; P = 0.008, beta; P = 0.005, low gamma; P = 0.017), whereas luteinizing hormone blood level was not correlated.

CONCLUSION

Resting-state brain activity can serve as an objective measurement of unpleasantness associated with menopausal symptoms, which aids the selection of appropriate treatment and monitors its outcome.

From bench to bedside: Overview of magnetoencephalography in basic principle, signal processing, source localization and clinical applications

Magnetoencephalography (MEG) is a non-invasive technique that can precisely capture the dynamic spatiotemporal patterns of the brain by measuring the magnetic fields arising from neuronal activity along the order of milliseconds. Observations of brain dynamics have been used in cognitive neuroscience, the diagnosis of neurological diseases, and the brain-computer interface (BCI). In this study, we outline the basic principle, signal processing, and source localization of MEG, and describe its clinical applications for cognitive assessment, the diagnoses of neurological diseases and mental disorders, preoperative evaluation, and the BCI. This review not only provides an overall perspective of MEG, ranging from practical techniques to clinical applications, but also enhances the prevalent understanding of neural mechanisms. The use of MEG is expected to lead to significant breakthroughs in neuroscience.

Magnetoencephalography Atlas Viewer for Dipole Localization and Viewing

Magnetoencephalography (MEG) is a noninvasive neuroimaging technique widely recognized for epilepsy and tumor mapping. MEG clinical reporting requires a multidisciplinary team, including expert input regarding each dipole's anatomic localization. Here, we introduce a novel tool, the "Magnetoencephalography Atlas Viewer" (MAV), which streamlines this anatomical analysis. The MAV normalizes the patient's Magnetic Resonance Imaging (MRI) to the Montreal Neurological Institute (MNI) space, reverse-normalizes MNI atlases to the native MRI, identifies MEG dipole files, and matches dipoles' coordinates to their spatial location in atlas files. It offers a user-friendly and interactive graphical user interface (GUI) for displaying individual dipoles, groups, coordinates, anatomical labels, and a tri-planar MRI view of the patient with dipole overlays. It evaluated over 273 dipoles obtained in clinical epilepsy subjects. Consensus-based ground truth was established by three neuroradiologists, with a minimum agreement threshold of two. The concordance between the ground truth and MAV labeling ranged from 79% to 84%, depending on the normalization method. Higher concordance rates were observed in subjects with minimal or no structural abnormalities on the MRI, ranging from 80% to 90%. The MAV provides a straightforward MEG dipole anatomic localization method, allowing a nonspecialist to prepopulate a report, thereby facilitating and reducing the time of clinical reporting.

Mindfulness meditation styles differently modulate source-level MEG microstate dynamics and complexity

Background

The investigation of mindfulness meditation practice, classically divided into focused attention meditation (FAM), and open monitoring meditation (OMM) styles, has seen a long tradition of theoretical, affective, neurophysiological and clinical studies. In particular, the high temporal resolution of magnetoencephalography (MEG) or electroencephalography (EEG) has been exploited to fill the gap between the personal experience of meditation practice and its neural correlates. Mounting evidence, in fact, shows that human brain activity is highly dynamic, transiting between different brain states (microstates). In this study, we aimed at exploring MEG microstates at source-level during FAM, OMM and in the resting state, as well as the complexity and criticality of dynamic transitions between microstates.

Methods

Ten right-handed Theravada Buddhist monks with a meditative expertise of minimum 2,265 h participated in the experiment. MEG data were acquired during a randomized block design task (6 min FAM, 6 min OMM, with each meditative block preceded and followed by 3 min resting state). Source reconstruction was performed using eLORETA on individual cortical space, and then parcellated according to the Human Connect Project atlas. Microstate analysis was then applied to parcel level signals in order to derive microstate topographies and indices. In addition, from microstate sequences, the Hurst exponent and the Lempel-Ziv complexity (LZC) were computed.

Results

Our results show that the coverage and occurrence of specific microstates are modulated either by being in a meditative state or by performing a specific meditation style. Hurst exponent values in both meditation conditions are reduced with respect to the value observed during rest, LZC shows significant differences between OMM, FAM, and REST, with a progressive increase from REST to FAM to OMM.

Discussion

Importantly, we report changes in brain criticality indices during meditation and between meditation styles, in line with a state-like effect of meditation on cognitive performance. In line with previous reports, we suggest that the change in cognitive state experienced in meditation is paralleled by a shift with respect to critical points in brain dynamics.

Oscillatory characteristics of resting-state magnetoencephalography reflect pathological and symptomatic conditions of cognitive impairment

Background

Dementia and mild cognitive impairment are characterised by symptoms of cognitive decline, which are typically assessed using neuropsychological assessments (NPAs), such as the Mini-Mental State Examination (MMSE) and Frontal Assessment Battery (FAB). Magnetoencephalography (MEG) is a novel clinical assessment technique that measures brain activities (summarised as oscillatory parameters), which are associated with symptoms of cognitive impairment. However, the relevance of MEG and regional cerebral blood flow (rCBF) data obtained using single-photon emission computed tomography (SPECT) has not been examined using clinical datasets. Therefore, this study aimed to investigate the relationships among MEG oscillatory parameters, clinically validated biomarkers computed from rCBF, and NPAs using outpatient data retrieved from hospital records.

Methods

Clinical data from 64 individuals with mixed pathological backgrounds were retrieved and analysed. MEG oscillatory parameters, including relative power (RP) from delta to high gamma bands, mean frequency, individual alpha frequency, and Shannon's spectral entropy, were computed for each cortical region. For SPECT data, three pathological parameters-'severity', 'extent', and 'ratio'-were computed using an easy z-score imaging system (eZIS). As for NPAs, the MMSE and FAB scores were retrieved.

Results

MEG oscillatory parameters were correlated with eZIS parameters. The eZIS parameters associated with Alzheimer's disease pathology were reflected in theta power augmentation and slower shift of the alpha peak. Moreover, MEG oscillatory parameters were found to reflect NPAs. Global slowing and loss of diversity in neural oscillatory components correlated with MMSE and FAB scores, whereas the associations between eZIS parameters and NPAs were sparse.

Conclusion

MEG oscillatory parameters correlated with both SPECT (i.e. eZIS) parameters and NPAs, supporting the clinical validity of MEG oscillatory parameters as pathological and symptomatic indicators. The findings indicate that various components of MEG oscillatory characteristics can provide valuable pathological and symptomatic information, making MEG data a rich resource for clinical examinations of patients with cognitive impairments. SPECT (i.e. eZIS) parameters showed no correlations with NPAs. The results contributed to a better understanding of the characteristics of electrophysiological and pathological examinations for patients with cognitive impairments, which will help to facilitate their co-use in clinical application, thereby improving patient care.

Connectivity-based Meta-Bands: A new approach for automatic frequency band identification in connectivity analyses

The majority of electroencephalographic (EEG) and magnetoencephalographic (MEG) studies filter and analyse neural signals in specific frequency ranges, known as "canonical" frequency bands. However, this segmentation, is not exempt from limitations, mainly due to the lack of adaptation to the neural idiosyncrasies of each individual. In this study, we introduce a new data-driven method to automatically identify frequency ranges based on the topological similarity of the frequency-dependent functional neural network. The resting-state neural activity of 195 cognitively healthy subjects from three different databases (MEG: 123 subjects; EEG1: 27 subjects; EEG2: 45 subjects) was analysed. In a first step, MEG and EEG signals were filtered with a narrow-band filter bank (1 Hz bandwidth) from 1 to 70 Hz with a 0.5 Hz step. Next, the connectivity in each of these filtered signals was estimated using the orthogonalized version of the amplitude envelope correlation to obtain the frequency-dependent functional neural network. Finally, a community detection algorithm was used to identify communities in the frequency domain showing a similar network topology. We have called this approach the "Connectivity-based Meta-Bands" (CMB) algorithm. Additionally, two types of synthetic signals were used to configure the hyper-parameters of the CMB algorithm. We observed that the classical approaches to band segmentation are partially aligned with the underlying network topologies at group level for the MEG signals, but they are missing individual idiosyncrasies that may be biasing previous studies, as revealed by our methodology. On the other hand, the sensitivity of EEG signals to reflect this underlying frequency-dependent network structure is limited, revealing a simpler frequency parcellation, not aligned with that defined by the "canonical" frequency bands. To the best of our knowledge, this is the first study that proposes an unsupervised band segmentation method based on the topological similarity of functional neural network across frequencies. This methodology fully accounts for subject-specific patterns, providing more robust and personalized analyses, and paving the way for new studies focused on exploring the frequency-dependent structure of brain connectivity.

Effect of channel density, inverse solutions and connectivity measures on EEG resting-state networks reconstruction: A simulation study

Along with the study of brain activity evoked by external stimuli, the past two decades witnessed an increased interest in characterizing the spontaneous brain activity occurring during resting conditions. The identification of connectivity patterns in this so-called "resting-state" has been the subject of a great number of electrophysiology-based studies, using the Electro/Magneto-Encephalography (EEG/MEG) source connectivity method. However, no consensus has been reached yet regarding a unified (if possible) analysis pipeline, and several involved parameters and methods require cautious tuning. This is particularly challenging when different analytical choices induce significant discrepancies in results and drawn conclusions, thereby hindering the reproducibility of neuroimaging research. Hence, our objective in this study was to shed light on the effect of analytical variability on outcome consistency by evaluating the implications of parameters involved in the EEG source connectivity analysis on the accuracy of resting-state networks (RSNs) reconstruction. We simulated, using neural mass models, EEG data corresponding to two RSNs, namely the default mode network (DMN) and dorsal attentional network (DAN). We investigated the impact of five channel densities (19, 32, 64, 128, 256), three inverse solutions (weighted minimum norm estimate (wMNE), exact low-resolution brain electromagnetic tomography (eLORETA), and linearly constrained minimum variance (LCMV) beamforming) and four functional connectivity measures (phase-locking value (PLV), phase-lag index (PLI), and amplitude envelope correlation (AEC) with and without source leakage correction), on the correspondence between reconstructed and reference networks. We showed that, with different analytical choices related to the number of electrodes, source reconstruction algorithm, and functional connectivity measure, high variability is present in the results. More specifically, our results show that a higher number of EEG channels significantly increased the accuracy of the reconstructed networks. Additionally, our results showed significant variability in the performance of the tested inverse solutions and connectivity measures. Such methodological variability and absence of analysis standardization represent a critical issue for neuroimaging studies that should be prioritized. We believe that this work could be useful for the field of electrophysiology connectomics, by increasing awareness regarding the challenge of variability in methodological approaches and its implications on reported results.

The temporal dynamics of Large-Scale brain network changes in disorders of consciousness: A Microstate-Based study

BACKGROUND AND OBJECTIVES

The resting-state brain is composed of several discrete networks, which remain stable for 10-100 ms. These functional microstates are considered the building blocks of spontaneous consciousness. Electroencephalography (EEG) microstate analysis may provide insight into the altered brain dynamics underlying consciousness recovery in patients with disorders of consciousness (DOC). We aimed to analyze microstates in the resting-state EEG source space in patients with DOC, the relationship between state-specific features and consciousness levels, and the corresponding patterns of microstates and functional networks.

METHODS

We obtained resting-state EEG data from 84 patients with DOC (27 in a minimally conscious state [MCS] and 57 in a vegetative state [VS] or with unresponsive wakefulness syndrome). We conducted a microstate analysis of the resting-state (EEG) source space and developed a state-transition analysis protocol for patients with DOC.

RESULTS

We identified seven microstates with distinct spatial distributions of cortical activation. Multivariate pattern analyses revealed that different functional connectivity patterns were associated with source-level microstates. There were significant differences in the microstate properties, including spatial activation patterns, temporal dynamics, state shifts, and connectivity construction, between the MCS and VS groups.

DISCUSSION

Our findings suggest that consciousness depends on complex dynamics within the brain and may originate from the anterior cortex.

NLGC: Network localized Granger causality with application to MEG directional functional connectivity analysis

Identifying the directed connectivity that underlie networked activity between different cortical areas is critical for understanding the neural mechanisms behind sensory processing. Granger causality (GC) is widely used for this purpose in functional magnetic resonance imaging analysis, but there the temporal resolution is low, making it difficult to capture the millisecond-scale interactions underlying sensory processing. Magnetoencephalography (MEG) has millisecond resolution, but only provides low-dimensional sensor-level linear mixtures of neural sources, which makes GC inference challenging. Conventional methods proceed in two stages: First, cortical sources are estimated from MEG using a source localization technique, followed by GC inference among the estimated sources. However, the spatiotemporal biases in estimating sources propagate into the subsequent GC analysis stage, may result in both false alarms and missing true GC links. Here, we introduce the Network Localized Granger Causality (NLGC) inference paradigm, which models the source dynamics as latent sparse multivariate autoregressive processes and estimates their parameters directly from the MEG measurements, integrated with source localization, and employs the resulting parameter estimates to produce a precise statistical characterization of the detected GC links. We offer several theoretical and algorithmic innovations within NLGC and further examine its utility via comprehensive simulations and application to MEG data from an auditory task involving tone processing from both younger and older participants. Our simulation studies reveal that NLGC is markedly robust with respect to model mismatch, network size, and low signal-to-noise ratio, whereas the conventional two-stage methods result in high false alarms and mis-detections. We also demonstrate the advantages of NLGC in revealing the cortical network-level characterization of neural activity during tone processing and resting state by delineating task- and age-related connectivity changes.

The interindividual variability of multimodal brain connectivity maintains spatial heterogeneity and relates to tissue microstructure

Humans differ from each other in a wide range of biometrics, but to what extent brain connectivity varies between individuals remains largely unknown. By combining diffusion-weighted imaging (DWI) and magnetoencephalography (MEG), this study characterizes the inter-subject variability (ISV) of multimodal brain connectivity. Structural connectivity is characterized by higher ISV in association cortices including the core multiple-demand network and lower ISV in the sensorimotor cortex. MEG ISV exhibits frequency-dependent signatures, and the extent of MEG ISV is consistent with that of structural connectivity ISV in selective macroscopic cortical clusters. Across the cortex, the ISVs of structural connectivity and beta-band MEG functional connectivity are negatively associated with cortical myelin content indexed by the quantitative T1 relaxation rate measured by high-resolution 7 T MRI. Furthermore, MEG ISV from alpha to gamma bands relates to the hindrance and restriction of the white-matter tissue estimated by DWI microstructural models. Our findings depict the inter-relationship between the ISV of brain connectivity from multiple modalities, and highlight the role of tissue microstructure underpinning the ISV.

Generative Models of Brain Dynamics

This review article gives a high-level overview of the approaches across different scales of organization and levels of abstraction. The studies covered in this paper include fundamental models in computational neuroscience, nonlinear dynamics, data-driven methods, as well as emergent practices. While not all of these models span the intersection of neuroscience, AI, and system dynamics, all of them do or can work in tandem as generative models, which, as we argue, provide superior properties for the analysis of neuroscientific data. We discuss the limitations and unique dynamical traits of brain data and the complementary need for hypothesis- and data-driven modeling. By way of conclusion, we present several hybrid generative models from recent literature in scientific machine learning, which can be efficiently deployed to yield interpretable models of neural dynamics.

Whole-Brain Network Models: From Physics to Bedside

Computational neuroscience has come a long way from its humble origins in the pioneering work of Hodgkin and Huxley. Contemporary computational models of the brain span multiple spatiotemporal scales, from single neuronal compartments to models of social cognition. Each spatial scale comes with its own unique set of promises and challenges. Here, we review models of large-scale neural communication facilitated by white matter tracts, also known as whole-brain models (WBMs). Whole-brain approaches employ inputs from neuroimaging data and insights from graph theory and non-linear systems theory to model brain-wide dynamics. Over the years, WBM models have shown promise in providing predictive insights into various facets of neuropathologies such as Alzheimer's disease, Schizophrenia, Epilepsy, Traumatic brain injury, while also offering mechanistic insights into large-scale cortical communication. First, we briefly trace the history of WBMs, leading up to the state-of-the-art. We discuss various methodological considerations for implementing a whole-brain modeling pipeline, such as choice of node dynamics, model fitting and appropriate parcellations. We then demonstrate the applicability of WBMs toward understanding various neuropathologies. We conclude by discussing ways of augmenting the biological and clinical validity of whole-brain models.

Globally elevated excitation-inhibition ratio in children with autism spectrum disorder and below-average intelligence

Background

Altered neuronal excitation–inhibition (E–I) balance is strongly implicated in ASD. However, it is not known whether the direction and degree of changes in the E–I ratio in individuals with ASD correlates with intellectual disability often associated with this developmental disorder. The spectral slope of the aperiodic 1/f activity reflects the E–I balance at the scale of large neuronal populations and may uncover its putative alternations in individuals with ASD with and without intellectual disability.

Methods

Herein, we used magnetoencephalography (MEG) to test whether the 1/f slope would differentiate ASD children with average and below–average (< 85) IQ. MEG was recorded at rest with eyes open/closed in 49 boys with ASD aged 6–15 years with IQ ranging from 54 to 128, and in 49 age-matched typically developing (TD) boys. The cortical source activity was estimated using the beamformer approach and individual brain models. We then extracted the 1/f slope by fitting a linear function to the log–log-scale power spectra in the high-frequency range.

Results

The global 1/f slope averaged over all cortical sources demonstrated high rank-order stability between the two conditions. Consistent with previous research, it was steeper in the eyes-closed than in the eyes-open condition and flattened with age. Regardless of condition, children with ASD and below-average IQ had flatter slopes than either TD or ASD children with average or above-average IQ. These group differences could not be explained by differences in signal-to-noise ratio or periodic (alpha and beta) activity.

Limitations

Further research is needed to find out whether the observed changes in E–I ratios are characteristic of children with below-average IQ of other diagnostic groups.

Conclusions

The atypically flattened spectral slope of aperiodic activity in children with ASD and below-average IQ suggests a shift of the global E–I balance toward hyper-excitation. The spectral slope can provide an accessible noninvasive biomarker of the E–I ratio for making objective judgments about treatment effectiveness in people with ASD and comorbid intellectual disability.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13229-022-00498-2.

MEG cortical microstates: Spatiotemporal characteristics, dynamic functional connectivity and stimulus-evoked responses

EEG microstate analysis is an approach to study brain states and their fast transitions in healthy cognition and disease. A key limitation of conventional microstate analysis is that it must be performed at the sensor level, and therefore gives limited anatomical insight. Here, we generalise the microstate methodology to be applicable to source-reconstructed electrophysiological data. Using simulations of a neural-mass network model, we first established the validity and robustness of the proposed method. Using MEG resting-state data, we uncovered ten microstates with distinct spatial distributions of cortical activation. Multivariate pattern analysis demonstrated that source-level microstates were associated with distinct functional connectivity patterns. We further demonstrated that the occurrence probability of MEG microstates were altered by auditory stimuli, exhibiting a hyperactivity of the microstate including the auditory cortex. Our results support the use of source-level microstates as a method for investigating brain dynamic activity and connectivity at the millisecond scale.

Optimizing EEG Source Reconstruction with Concurrent fMRI-Derived Spatial Priors

Reconstructing EEG sources involves a complex pipeline, with the inverse problem being the most challenging. Multiple inversion algorithms are being continuously developed, aiming to tackle the non-uniqueness of this problem, which has been shown to be partially circumvented by including prior information in the inverse models. Despite a few efforts, there are still current and persistent controversies regarding the inversion algorithm of choice and the optimal set of spatial priors to be included in the inversion models. The use of simultaneous EEG-fMRI data is one approach to tackle this problem. The spatial resolution of fMRI makes fMRI derived spatial priors very convenient for EEG reconstruction, however, only task activation maps and resting-state networks (RSNs) have been explored so far, overlooking the recent, but already accepted, notion that brain networks exhibit dynamic functional connectivity fluctuations. The lack of a systematic comparison between different source reconstruction algorithms, considering potentially more brain-informative priors such as fMRI, motivates the search for better reconstruction models. Using simultaneous EEG-fMRI data, here we compared four different inversion algorithms (minimum norm, MN; low resolution electromagnetic tomography, LORETA; empirical Bayes beamformer, EBB; and multiple sparse priors, MSP) under a Bayesian framework (as implemented in SPM), each with three different sets of priors consisting of: (1) those specific to the algorithm; (2) those specific to the algorithm plus fMRI task activation maps and RSNs; and (3) those specific to the algorithm plus fMRI task activation maps and RSNs and network modules of task-related dFC states estimated from the dFC fluctuations. The quality of the reconstructed EEG sources was quantified in terms of model-based metrics, namely the expectation of the posterior probability P(model|data) and variance explained of the inversion models, and the overlap/proportion of brain regions known to be involved in the visual perception tasks that the participants were submitted to, and RSN templates, with/within EEG source components. Model-based metrics suggested that model parsimony is preferred, with the combination MSP and priors specific to this algorithm exhibiting the best performance. However, optimal overlap/proportion values were found using EBB and priors specific to this algorithm and fMRI task activation maps and RSNs or MSP and considering all the priors (algorithm priors, fMRI task activation maps and RSNs and dFC state modules), respectively, indicating that fMRI spatial priors, including dFC state modules, might contain useful information to recover EEG source components reflecting neuronal activity of interest. Our main results show that providing fMRI spatial derived priors that reflect the dynamics of the brain might be useful to map neuronal activity more accurately from EEG-fMRI. Furthermore, this work paves the way towards a more informative selection of the optimal EEG source reconstruction approach, which may be critical in future studies.

Dementia ConnEEGtome: Towards multicentric harmonization of EEG connectivity in neurodegeneration

The proposal to use brain connectivity as a biomarker for dementia phenotyping can be potentiated by conducting large-scale multicentric studies using high-density electroencephalography (hd- EEG). Nevertheless, several barriers preclude the development of a systematic "ConnEEGtome" in dementia research. Here we review critical sources of variability in EEG connectivity studies, and provide general guidelines for multicentric protocol harmonization. We describe how results can be impacted by the choice for data acquisition, and signal processing workflows. The implementation of a particular processing pipeline is conditional upon assumptions made by researchers about the nature of EEG. Due to these assumptions, EEG connectivity metrics are typically applicable to restricted scenarios, e.g., to a particular neurocognitive disorder. "Ground truths" for the choice of processing workflow and connectivity analysis are impractical. Consequently, efforts should be directed to harmonizing experimental procedures, data acquisition, and the first steps of the preprocessing pipeline. Conducting multiple analyses of the same data and a proper integration of the results need to be considered in additional processing steps. Furthermore, instead of using a single connectivity measure, using a composite metric combining different connectivity measures brings a powerful strategy to scale up the replicability of multicentric EEG connectivity studies. These composite metrics can boost the predictive strength of diagnostic tools for dementia. Moreover, the implementation of multi-feature machine learning classification systems that include EEG-based connectivity analyses may help to exploit the potential of multicentric studies combining clinical-cognitive, molecular, genetics, and neuroimaging data towards a multi-dimensional characterization of the dementia.

A systematic evaluation of source reconstruction of resting MEG of the human brain with a new high-resolution atlas: Performance, precision, and parcellation

Noninvasive functional neuroimaging of the human brain can give crucial insight into the mechanisms that underpin healthy cognition and neurological disorders. Magnetoencephalography (MEG) measures extracranial magnetic fields originating from neuronal activity with high temporal resolution, but requires source reconstruction to make neuroanatomical inferences from these signals. Many source reconstruction algorithms are available, and have been widely evaluated in the context of localizing task-evoked activities. However, no consensus yet exists on the optimum algorithm for resting-state data. Here, we evaluated the performance of six commonly-used source reconstruction algorithms based on minimum-norm and beamforming estimates. Using human resting-state MEG, we compared the algorithms using quantitative metrics, including resolution properties of inverse solutions and explained variance in sensor-level data. Next, we proposed a data-driven approach to reduce the atlas from the Human Connectome Project's multi-modal parcellation of the human cortex based on metrics such as MEG signal-to-noise-ratio and resting-state functional connectivity gradients. This procedure produced a reduced cortical atlas with 230 regions, optimized to match the spatial resolution and the rank of MEG data from the current generation of MEG scanners. Our results show that there is no "one size fits all" algorithm, and make recommendations on the appropriate algorithms depending on the data and aimed analyses. Our comprehensive comparisons and recommendations can serve as a guide for choosing appropriate methodologies in future studies of resting-state MEG.

Mean-Field Modeling of Brain-Scale Dynamics for the Evaluation of EEG Source-Space Networks

Understanding the dynamics of brain-scale functional networks at rest and during cognitive tasks is the subject of intense research efforts to unveil fundamental principles of brain functions. To estimate these large-scale brain networks, the emergent method called "electroencephalography (EEG) source connectivity" has generated increasing interest in the network neuroscience community, due to its ability to identify cortical brain networks with satisfactory spatio-temporal resolution, while reducing mixing and volume conduction effects. However, no consensus has been reached yet regarding a unified EEG source connectivity pipeline, and several methodological issues have to be carefully accounted to avoid pitfalls. Thus, a validation toolbox that provides flexible "ground truth" models is needed for an objective methods/parameters evaluation and, thereby an optimization of the EEG source connectivity pipeline. In this paper, we show how a recently developed large-scale model of brain-scale activity, named COALIA, can provide to some extent such ground truth by providing realistic simulations of source-level and scalp-level activity. Using a bottom-up approach, the model bridges cortical micro-circuitry and large-scale network dynamics. Here, we provide an example of the potential use of COALIA to analyze, in the context of epileptiform activity, the effect of three key factors involved in the "EEG source connectivity" pipeline: (i) EEG sensors density, (ii) algorithm used to solve the inverse problem, and (iii) functional connectivity measure. Results showed that a high electrode density (at least 64 channels) is required to accurately estimate cortical networks. Regarding the inverse solution/connectivity measure combination, the best performance at high electrode density was obtained using the weighted minimum norm estimate (wMNE) combined with the weighted phase lag index (wPLI). Although those results are specific to the considered aforementioned context (epileptiform activity), we believe that this model-based approach can be successfully applied to other experimental questions/contexts. We aim at presenting a proof-of-concept of the interest of COALIA in the network neuroscience field, and its potential use in optimizing the EEG source-space network estimation pipeline.