Comparing MEG and high-density EEG for intrinsic functional connectivity mapping

Simultaneous invasive and non-invasive recordings in humans: A novel Rosetta stone for deciphering brain activity

Simultaneous noninvasive and invasive electrophysiological recordings provide a unique opportunity to achieve a comprehensive understanding of human brain activity, much like a Rosetta stone for human neuroscience. In this review we focus on the increasingly-used powerful combination of intracranial electroencephalography (iEEG) with scalp electroencephalography (EEG) or magnetoencephalography (MEG). We first provide practical insight on how to achieve these technically challenging recordings. We then provide examples from clinical research on how simultaneous recordings are advancing our understanding of epilepsy. This is followed by the illustration of how human neuroscience and methodological advances could benefit from these simultaneous recordings. We conclude with a call for open data sharing and collaboration, while ensuring neuroethical approaches and argue that only with a true collaborative approach the promises of simultaneous recordings will be fulfilled.

Non-invasive measurement of rat auditory evoked fields using an optically pumped atomic magnetometer: Effects of task manipulation

Optically pumped magnetometers (OPMs) have become a favorable tool for magnetoencephalography (MEG) measurement, offering a non-invasive method of measurement. OPMs do not require cryogenic environments, sensors can be more closely aligned with the brain. We employed a passive single-stimulus paradigm in conjunction with OPMs with a sensitivity of 20 fT/Hz to investigate the auditory response of rats to inter-stimulus interval (ISI) and frequencies, recording the rat auditory event-related magnetic fields (ERMFs). Our findings include: (1) Auditory evoked fields can be detected non-invasively by OPMs; (2) The amplitude of the rat auditory ERMFs varies with changes in ISI, with more pronounced amplitude changes observed after 5 s; (3) When the sound stimulus frequency is altered at the same ISI, the amplitude of the rats ERMFs changes with frequency, indicating significant differences in attention. Our method offers a valuable tool for the clinical application of a single stimulus paradigm and opens up a new avenue for research on the brain magnetic field detections.

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.

Navigating Pubertal Goldilocks: The Optimal Pace for Hierarchical Brain Organization

Adolescence is a timed process with an onset, tempo, and duration. Nevertheless, the temporal dimension, especially the pace of maturation, remains an insufficiently studied aspect of developmental progression. The primary objective is to estimate the precise influence of pubertal maturational tempo on the configuration of associative brain regions. To this end, the connection between maturational stages and the level of hierarchical organization of large-scale brain networks in 12-13-year-old females is analyzed. Skeletal maturity is used as a proxy for pubertal progress. The degree of maturity is defined by the difference between bone age and chronological age. To assess the level of hierarchical organization in the brain, the temporal dynamic of closed eye resting state high-density electroencephalography (EEG) in the alpha frequency range is analyzed. Different levels of hierarchical order are captured by the measured asymmetry in the directionality of information flow between different regions. The calculated EEG-based entropy production of participant groups is then compared with accelerated, average, and decelerated maturity. Results indicate that an average maturational trajectory optimally aligns with cerebral hierarchical order, and both accelerated and decelerated timelines result in diminished cortical organization. This suggests that a "Goldilocks rule" of brain development is favoring a particular maturational tempo.

CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

We present CiftiStorm, an electrophysiological source imaging (ESI) pipeline incorporating recently developed methods to improve forward and inverse solutions. The CiftiStorm pipeline produces Human Connectome Project (HCP) and megconnectome-compliant outputs from dataset inputs with varying degrees of spatial resolution. The input data can range from low-sensor-density electroencephalogram (EEG) or magnetoencephalogram (MEG) recordings without structural magnetic resonance imaging (sMRI) to high-density EEG/MEG recordings with an HCP multimodal sMRI compliant protocol. CiftiStorm introduces a numerical quality control of the lead field and geometrical corrections to the head and source models for forward modeling. For the inverse modeling, we present a Bayesian estimation of the cross-spectrum of sources based on multiple priors. We facilitate ESI in the T1w/FSAverage32k high-resolution space obtained from individual sMRI. We validate this feature by comparing CiftiStorm outputs for EEG and MRI data from the Cuban Human Brain Mapping Project (CHBMP) acquired with technologies a decade before the HCP MEG and MRI standardized dataset.

Using optically pumped magnetometers to replicate task-related responses in next generation magnetoencephalography

Optically pumped magnetometers (OPMs) offer a new wearable means to measure magnetoencephalography (MEG) signals, with many advantages compared to conventional systems. However, OPMs are an emerging technology, thus characterizing and replicating MEG recordings is essential. Using OPM-MEG and SQUID-MEG, this study investigated evoked responses, oscillatory power, and functional connectivity during emotion processing in 20 adults, to establish replicability across the two technologies. Five participants with dental fixtures were included to assess the validity of OPM-MEG recordings in those with irremovable metal. Replicable task-related evoked responses were observed in both modalities. Similar patterns of oscillatory power to faces were observed in both systems. Increased connectivity was found in SQUID-versus OPM-MEG in an occipital and parietal anchored network. Notably, high quality OPM-MEG data were retained in participants with metallic fixtures, from whom no useable data were collected using conventional MEG.

Exploring the limits of MEG spatial resolution with multipolar expansions

The advent of scalp magnetoencephalography (MEG) based on optically pumped magnetometers (OPMs) may represent a step change in the field of human electrophysiology. Compared to cryogenic MEG based on superconducting quantum interference devices (SQUIDs, placed 2-4 cm above scalp), scalp MEG promises significantly higher spatial resolution imaging but it also comes with numerous challenges regarding how to optimally design OPM arrays. In this context, we sought to provide a systematic description of MEG spatial resolution as a function of the number of sensors (allowing comparison of low- vs. high-density MEG), sensor-to-brain distance (cryogenic SQUIDs vs. scalp OPM), sensor type (magnetometers vs. gradiometers; single- vs. multi-component sensors), and signal-to-noise ratio. To that aim, we present an analytical theory based on MEG multipolar expansions that enables, once supplemented with experimental input and simulations, quantitative assessment of the limits of MEG spatial resolution in terms of two qualitatively distinct regimes. In the regime of asymptotically high-density MEG, we provide a mathematically rigorous description of how magnetic field smoothness constraints spatial resolution to a slow, logarithmic divergence. In the opposite regime of low-density MEG, it is sensor density that constraints spatial resolution to a faster increase following a square-root law. The transition between these two regimes controls how MEG spatial resolution saturates as sensors approach sources of neural activity. This two-regime model of MEG spatial resolution integrates known observations (e.g., the difficulty of improving spatial resolution by increasing sensor density, the gain brought by moving sensors on scalp, or the usefulness of multi-component sensors) and gathers them under a unifying theoretical framework that highlights the underlying physics and reveals properties inaccessible to simulations. We propose that this framework may find useful applications to benchmark the design of future OPM-based scalp MEG systems.

Parietal alpha underlies slower cognitive responses during interference processing in adolescents

Adolescence is a critical period when cognitive control is rapidly maturing across several core dimensions. Here, we evaluated how healthy adolescents (13-17 years of age, n = 44) versus young adults (18-25 years of age, n = 49) differ across a series of cognitive assessments with simultaneous electroencephalography (EEG) recordings. Cognitive tasks included selective attention, inhibitory control, working memory, as well as both non-emotional and emotional interference processing. We found that adolescents displayed significantly slower responses than young adults specifically on the interference processing tasks. Evaluation of EEG event-related spectral perturbations (ERSPs) on the interference tasks showed that adolescents consistently had greater event-related desynchronization in alpha/beta frequencies in parietal regions. Midline frontal theta activity was also greater in the flanker interference task in adolescents, suggesting greater cognitive effort. Parietal alpha activity predicted age-related speed differences during non-emotional flanker interference processing, and frontoparietal connectivity, specifically midfrontal theta - parietal alpha functional connectivity predicted speed effects during emotional interference. Overall, our neuro-cognitive results illustrate developing cognitive control in adolescents particularly for interference processing, predicted by differential alpha band activity and connectivity in parietal brain regions.

Impact of secondary particles on the magnetic field generated by a proton pencil beam: a finite-element analysis based on Geant4-DNA simulations

PURPOSE

To investigate the static magnetic field generated by a proton pencil beam as a candidate for range verification by means of Monte Carlo simulations, thereby improving upon existing analytical calculations. We focus on the impact of statistical current fluctuations and secondary protons and electrons.

METHODS

We considered a pulsed beam (10 μ ${\umu}$ s pulse duration) during the duty cycle with a peak beam current of 0.2 μ $\umu$ A and an initial energy of 100 MeV. We ran Geant4-DNA Monte Carlo simulations of a proton pencil beam in water and extracted independent particle phase spaces. We calculated longitudinal and radial current density of protons and electrons, serving as an input for a magnetic field estimation based on a finite element analysis in a cylindrical geometry. We made sure to allow for non-solenoidal current densities as is the case of a stopping proton beam.

RESULTS

The rising proton charge density toward the range is not perturbed by energy straggling and only lowered through nuclear reactions by up to 15%, leading to an approximately constant longitudinal current. Their relative low density however (at most 0.37 protons/mm³ for the 0.2  μ ${\umu}$ A current and a beam cross-section of 2.5 mm), gives rise to considerable current density fluctuations. The radial proton current resulting from lateral scattering and being two orders of magnitude weaker than the longitudinal current is subject to even stronger fluctuations. Secondary electrons with energies above 10 eV, that far outnumber the primary protons, reduce the primary proton current by only 10% due to their largely isotropic flow. A small fraction of electrons (<1%), undergoing head-on collisions, constitutes the relevant electron current. In the far-field, both contributions to the magnetic field strength (longitudinal and lateral) are independent of the beam spot size. We also find that the nuclear reaction-related losses cause a shift of 1.3 mm to the magnetic field profile relative to the actual range, which is further enlarged to 2.4 mm by the electron current (at a distance of ρ = 50 $\rho =50$  mm away from the central beam axis). For ρ > 45 $\rho &gt;45$  mm, the shift increases linearly. While the current density variations cause significant magnetic field uncertainty close to the central beam axis with a relative standard deviation (RSD) close to 100%, they average out at a distance of 10 cm, where the RSD of the total magnetic field drops below 2%.

CONCLUSIONS

With the small influence of the secondary electrons together with the low RSD, our analysis encourages an experimental detection of the magnetic field through sensitive instrumentation, such as optical magnetometry or SQUIDs.

Graph Theory Measures and Their Application to Neurosurgical Eloquence

Improving patient safety and preserving eloquent brain are crucial in neurosurgery. Since there is significant clinical variability in post-operative lesions suffered by patients who undergo surgery in the same areas deemed compensable, there is an unknown degree of inter-individual variability in brain 'eloquence'. Advances in connectomic mapping efforts through diffusion tractography allow for utilization of non-invasive imaging and statistical modeling to graphically represent the brain. Extending the definition of brain eloquence to graph theory measures of hubness and centrality may help to improve our understanding of individual variability in brain eloquence and lesion responses. While functional deficits cannot be immediately determined intra-operatively, there has been potential shown by emerging technologies in mapping of hub nodes as an add-on to existing surgical navigation modalities to improve individual surgical outcomes. This review aims to outline and review current research surrounding novel graph theoretical concepts of hubness, centrality, and eloquence and specifically its relevance to brain mapping for pre-operative planning and intra-operative navigation in neurosurgery.

Distinct roles for the anterior temporal lobe and angular gyrus in the spatiotemporal cortical semantic network

Semantic knowledge is supported by numerous brain regions, but the spatiotemporal configuration of the network that links these areas remains an open question. The hub-and-spokes model posits that a central semantic hub coordinates this network. In this study, we explored distinct aspects that define a semantic hub, as reflected in the spatiotemporal modulation of neural activity and connectivity by semantic variables, from the earliest stages of semantic processing. We used source-reconstructed electro/magnetoencephalography, and investigated the concreteness contrast across three tasks. In a whole-cortex analysis, the left anterior temporal lobe (ATL) was the only area that showed modulation of evoked brain activity from 100 ms post-stimulus. Furthermore, using Dynamic Causal Modeling of the evoked responses, we investigated effective connectivity amongst the candidate semantic hub regions, that is, left ATL, supramarginal/angular gyrus (SMG/AG), middle temporal gyrus, and inferior frontal gyrus. We found that models with a single semantic hub showed the highest Bayesian evidence, and the hub region was found to change from ATL (within 250 ms) to SMG/AG (within 450 ms) over time. Our results support a single semantic hub view, with ATL showing sustained modulation of neural activity by semantics, and both ATL and AG underlying connectivity depending on the stage of semantic processing.

A systematic data-driven approach to analyze sensor-level EEG connectivity: Identifying robust phase-synchronized network components using surface Laplacian with spectral-spatial PCA

Although conventional averaging across predefined frequency bands reduces the complexity of EEG functional connectivity (FC), it obscures the identification of resting-state brain networks (RSN) and impedes accurate estimation of FC reliability. Extending prior work, we combined scalp current source density (CSD; spherical spline surface Laplacian) and spectral-spatial PCA to identify FC components. Phase-based FC was estimated via debiased-weighted phase-locking index from CSD-transformed resting EEGs (71 sensors, 8 min, eyes open/closed, 35 healthy adults, 1-week retest). Spectral PCA extracted six robust alpha and theta components (86.6% variance). Subsequent spatial PCA for each spectral component revealed seven robust regionally focused (posterior, central, and frontal) and long-range (posterior-anterior) alpha components (peaks at 8, 10, and 13 Hz) and a midfrontal theta (6 Hz) component, accounting for 37.0% of FC variance. These spatial FC components were consistent with well-known networks (e.g., default mode, visual, and sensorimotor), and four were sensitive to eyes open/closed conditions. Most FC components had good-to-excellent internal consistency (odd/even epochs, eyes open/closed) and test-retest reliability (ICCs ≥ .8). Moreover, the FC component structure was generally present in subsamples (session × odd/even epoch, or smaller subgroups [n = 7-10]), as indicated by high similarity of component loadings across PCA solutions. Apart from systematically reducing FC dimensionality, our approach avoids arbitrary thresholds and allows quantification of meaningful and reliable network components that may prove to be of high relevance for basic and clinical research applications.

Combined electrophysiological and morphological phenotypes in patients with genetic generalized epilepsy and their healthy siblings

OBJECTIVE

Genetic generalized epilepsy is characterized by aberrant neuronal dynamics and subtle structural alterations. We evaluated whether a combination of magnetic and electrical neuronal signals and cortical thickness would provide complementary information about network pathology in GGE. We also investigated if these imaging phenotypes were present in healthy siblings of the patients to test for genetic influence.

METHODS

In this cross-sectional study, we analyzed five minutes of resting-state data acquired using electroencephalography (EEG) and magnetoencephalography (MEG) in patients, their siblings, and controls, matched for age and sex. We computed source-reconstructed power and connectivity in six frequency bands (1-40 Hz) and cortical thickness (derived from magnetic resonance imaging (MRI)). Group differences were assessed using permutation analysis of linear models for each modality separately and jointly for all modalities using a non-parametric combination.

RESULTS

Patients with GGE (n = 23) had higher power than controls (n = 35) in all frequencies, with a more posterior focus in MEG than EEG. Connectivity was also increased, particularly in frontotemporal and central regions in theta (strongest in EEG) and low beta frequencies (strongest in MEG), which was eminent in the joint EEG/MEG analysis. EEG showed weaker connectivity differences in higher frequencies, possibly related to drug effects. The inclusion of cortical thickness reinforced group differences in connectivity and power. Siblings (n = 18) had functional and structural patterns intermediate between those of patients and controls.

SIGNIFICANCE

EEG detected increased connectivity and power in GGE similar to MEG, but with different spectral sensitivity, highlighting the importance of theta and beta oscillations. Cortical thickness reductions in GGE corresponded to functional imaging patterns. Our multimodal approach extends the understanding of the resting-state in GGE and points to genetic underpinnings of the imaging markers studied, providing new insights into the causes and consequences of epilepsy.

Atypical resting-state functional brain connectivity in children with developmental coordination disorder

•

Atypical connectivity in children with developmental coordination disorder.

•

Stronger connectivity mainly found within the dorsal extrastriate network.

•

May reflect a brain trait of children with developmental coordination disorder.

•

This atypical connectivity is not associated with motor/visual perceptual abilities.

•

Lower visuomotor performance associated with stronger sensorimotor connectivity.

Children with developmental coordination disorder (DCD) present lower abilities to acquire and execute coordinated motor skills. DCD is frequently associated with visual perceptual (with or without motor component) impairments. This magnetoencephalography (MEG) study compares the brain resting-state functional connectivity (rsFC) and spectral power of children with and without DCD.

29 children with DCD and 28 typically developing (TD) peers underwent 2 × 5 min of resting-state MEG. Band-limited power envelope correlation and spectral power were compared between groups using a functional connectome of 59 nodes from eight resting-state networks. Correlation coefficients were calculated between fine and gross motor activity, visual perceptual and visuomotor abilities measures on the one hand, and brain rsFC and spectral power on the other hand. Nonparametric statistics were used.

Significantly higher rsFC between nodes of the visual, attentional, frontoparietal, default-mode and cerebellar networks was observed in the alpha (maximum statistics, p = .0012) and the low beta (p = .0002) bands in children with DCD compared to TD peers. Lower visuomotor performance (copying figures) was associated with stronger interhemispheric rsFC within sensorimotor areas and power in the cerebellum (right lobule VIII).

Children with DCD showed increased rsFC mainly in the dorsal extrastriate visual brain system and the cerebellum. However, this increase was not associated with their coordinated motor/visual perceptual abilities. This enhanced functional brain connectivity could thus reflect a characteristic brain trait of children with DCD compared to their TD peers. Moreover, an interhemispheric compensatory process might be at play to perform visuomotor task within the normative range.

Assessing spino-cortical proprioceptive processing in childhood unilateral cerebral palsy with corticokinematic coherence

OBJECTIVE

To develop an electrophysiological marker of proprioceptive spino-cortical tracts integrity based on corticokinematic coherence (CKC) in young children with unilateral cerebral palsy (UCP), in whom behavioral measures are not applicable.

METHODS

Electroencephalography (EEG) signals from 12 children with UCP aged 19 to 57 months were recorded using 128-channel EEG caps while their fingers were moved at 2 Hz by an experimenter, in separate sessions for the affected and non-affected hands. The coherence between movement kinematics and EEG signals (i.e., CKC) was computed at the sensor and source (using a realistic head model) levels. Peaks of CKC obtained for the affected and non-affected hands were compared for location and strength. The relation between CKC strength on the lesion-side, the lesion-type (cortico-subcortical vs. subcortical) and the level of manual ability were studied with 2-way repeated-measures ANOVA.

RESULTS

At the individual level, a significant CKC peak at the central area contralateral to the moved hand was found in all young children with their non-affected hand and in 8 out of 12 children with their affected hand. At the group level, CKC to the affected hand movements was weaker than CKC to the non-affected hand movements. This difference was influenced by the type of lesion, the effect being predominant in the subgroup (n = 5) with cortico-subcortical lesions.

CONCLUSION

CKC is measurable with EEG in young children with UCP and provides electrophysiological evidence for altered proprioceptive spino-cortical tracts on the lesioned brain hemisphere, particularly in children with cortico-subcortical lesions.

Microstates and power envelope hidden Markov modeling probe bursting brain activity at different timescales

State modeling of whole-brain electroencephalography (EEG) or magnetoencephalography (MEG) allows to investigate transient, recurring neurodynamical events. Two widely-used techniques are the microstate analysis of EEG signals and hidden Markov modeling (HMM) of MEG power envelopes. Both reportedly lead to similar state lifetimes on the 100 ms timescale, suggesting a common neural basis. To investigate whether microstates and power envelope HMM states describe the same neural dynamics, we used simultaneous MEG/EEG recordings at rest and compared the spatial signature and temporal activation dynamics of microstates and power envelope HMM states obtained separately from EEG and MEG. Results showed that microstates and power envelope HMM states differ both spatially and temporally. Microstates reflect sharp events of neural synchronization, whereas power envelope HMM states disclose network-level activity with 100–200 ms lifetimes. Further, MEG microstates do not correspond to the canonical EEG microstates but are better interpreted as split HMM states. On the other hand, both MEG and EEG HMM states involve the (de)activation of similar functional networks. Microstate analysis and power envelope HMM thus appear sensitive to neural events occurring over different spatial and temporal scales. As such, they represent complementary approaches to explore the fast, sub-second scale bursting electrophysiological dynamics in spontaneous human brain activity.

MEG Node Degree Differences in Patients with Focal Epilepsy vs. Controls-Influence of Experimental Conditions

Drug-resistant epilepsy can be most limiting for patients, and surgery represents a viable therapy option. With the growing research on the human connectome and the evidence of epilepsy being a network disorder, connectivity analysis may be able to contribute to our understanding of epilepsy and may be potentially developed into clinical applications. In this magnetoencephalographic study, we determined the whole-brain node degree of connectivity levels in patients and controls. Resting-state activity was measured at five frequency bands in 15 healthy controls and 15 patients with focal epilepsy of different etiologies. The whole-brain all-to-all imaginary part of coherence in source space was then calculated. Node degree was determined and parcellated and was used for further statistical evaluation. In comparison to controls, we found a significantly higher overall node degree in patients with lesional and non-lesional epilepsy. Furthermore, we examined the conditions of high/reduced vigilance and open/closed eyes in controls, to analyze whether patient node degree levels can be achieved. We evaluated intraclass-correlation statistics (ICC) to evaluate the reproducibility. Connectivity and specifically node degree analysis could present new tools for one of the most common neurological diseases, with potential applications in epilepsy diagnostics.

Improved spatio-temporal measurements of visually evoked fields using optically-pumped magnetometers

Recent developments in performance and practicality of optically-pumped magnetometers (OPMs) have enabled new capabilities in non-invasive brain function mapping through magnetoencephalography. In particular, the lack of cryogenic operating conditions allows for more flexible placement of sensor heads closer to the brain, leading to improved spatial resolution and source localisation capabilities. Through recording visually evoked brain fields (VEFs), we demonstrate that the closer sensor proximity can be exploited to improve temporal resolution. We use OPMs, and superconducting quantum interference devices (SQUIDs) for reference, to measure brain responses to flash and pattern reversal stimuli. We find highly reproducible signals with consistency across multiple participants, stimulus paradigms and sensor modalities. The temporal resolution advantage of OPMs is manifest in a twofold improvement, compared to SQUIDs. The capability for improved spatio-temporal signal tracing is illustrated by simultaneous vector recordings of VEFs in the primary and associative visual cortex, where a time lag on the order of 10–20 ms is consistently found. This paves the way for further spatio-temporal studies of neurophysiological signal tracking in visual stimulus processing, and other brain responses, with potentially far-reaching consequences for time-critical mapping of functionality in healthy and pathological brains.

MEGnet: Automatic ICA-based artifact removal for MEG using spatiotemporal convolutional neural networks

Magnetoencephalography (MEG) is a functional neuroimaging tool that records the magnetic fields induced by neuronal activity; however, signal from non-neuronal sources can corrupt the data. Eye-blinks, saccades, and cardiac activity are three of the most common sources of non-neuronal artifacts. They can be measured by affixing eye proximal electrodes, as in electrooculography (EOG), and chest electrodes, as in electrocardiography (ECG), however this complicates imaging setup, decreases patient comfort, and can induce further artifacts from movement. This work proposes an EOG- and ECG-free approach to identify eye-blinks, saccades, and cardiac activity signals for automated artifact suppression. The contribution of this work is three-fold. First, using a data driven, multivariate decomposition approach based on Independent Component Analysis (ICA), a highly accurate artifact classifier is constructed as an amalgam of deep 1-D and 2-D Convolutional Neural Networks (CNNs) to automate the identification and removal of ubiquitous whole brain artifacts including eye-blink, saccade, and cardiac artifacts. The specific architecture of this network is optimized through an unbiased, computer-based hyperparameter random search. Second, visualization methods are applied to the learned abstraction to reveal what features the model uses and to bolster user confidence in the model's training and potential for generalization. Finally, the model is trained and tested on both resting-state and task MEG data from 217 subjects, and achieves a new state-of-the-art in artifact detection accuracy of 98.95% including 96.74% sensitivity and 99.34% specificity on the held out test-set. This work automates MEG processing for both clinical and research use, adapts to the acquired acquisition time, and can obviate the need for EOG or ECG electrodes for artifact detection.

Investigating the spectral features of the brain meso-scale structure at rest

Recent studies provide novel insights into the meso-scale organization of the brain, highlighting the co-occurrence of different structures: classic assortative (modular), disassortative, and core-periphery. However, the spectral properties of the brain meso-scale remain mostly unexplored. To fill this knowledge gap, we investigated how the meso-scale structure is organized across the frequency domain. We analyzed the resting state activity of healthy participants with source-localized high-density electroencephalography signals. Then, we inferred the community structure using weighted stochastic block-model (WSBM) to capture the landscape of meso-scale structures across the frequency domain. We found that different meso-scale modalities co-exist and are diversely organized over the frequency spectrum. Specifically, we found a core-periphery structure dominance, but we also highlighted a selective increase of disassortativity in the low frequency bands (<8 Hz), and of assortativity in the high frequency band (30-50 Hz). We further described other features of the meso-scale organization by identifying those brain regions which, at the same time, (a) exhibited the highest degree of assortativity, disassortativity, and core-peripheriness (i.e., participation) and (b) were consistently assigned to the same community, irrespective from the granularity imposed by WSBM (i.e., granularity-invariance). In conclusion, we observed that the brain spontaneous activity shows frequency-specific meso-scale organization, which may support spatially distributed and local information processing.

Age of onset modulates resting-state brain network dynamics in Friedreich Ataxia

This magnetoencephalography (MEG) study addresses (i) how Friedreich ataxia (FRDA) affects the sub‐second dynamics of resting‐state brain networks, (ii) the main determinants of their dynamic alterations, and (iii) how these alterations are linked with FRDA‐related changes in resting‐state functional brain connectivity (rsFC) over long timescales. For that purpose, 5 min of resting‐state MEG activity were recorded in 16 FRDA patients (mean age: 27 years, range: 12–51 years; 10 females) and matched healthy subjects. Transient brain network dynamics was assessed using hidden Markov modeling (HMM). Post hoc median‐split, nonparametric permutations and Spearman rank correlations were used for statistics. In FRDA patients, a positive correlation was found between the age of symptoms onset (ASO) and the temporal dynamics of two HMM states involving the posterior default mode network (DMN) and the temporo‐parietal junctions (TPJ). FRDA patients with an ASO <11 years presented altered temporal dynamics of those two HMM states compared with FRDA patients with an ASO > 11 years or healthy subjects. The temporal dynamics of the DMN state also correlated with minute‐long DMN rsFC. This study demonstrates that ASO is the main determinant of alterations in the sub‐second dynamics of posterior associative neocortices in FRDA patients and substantiates a direct link between sub‐second network activity and functional brain integration over long timescales.

A Transcriptome Community-and-Module Approach of the Human Mesoconnectome

Graph analysis allows exploring transcriptome compartments such as communities and modules for brain mesostructures. In this work, we proposed a bottom-up model of a gene regulatory network to brain-wise connectome workflow. We estimated the gene communities across all brain regions from the Allen Brain Atlas transcriptome database. We selected the communities method to yield the highest number of functional mesostructures in the network hierarchy organization, which allowed us to identify specific brain cell functions (e.g., neuroplasticity, axonogenesis and dendritogenesis communities). With these communities, we built brain-wise region modules that represent the connectome. Our findings match with previously described anatomical and functional brain circuits, such the default mode network and the default visual network, supporting the notion that the brain dynamics that carry out low- and higher-order functions originate from the modular composition of a GRN complex network.

Resting-state functional brain connectivity is related to subsequent procedural learning skills in school-aged children

This magnetoencephalography (MEG) study investigates how procedural sequence learning performance is related to prior brain resting-state functional connectivity (rsFC), and to what extent sequence learning induces rapid changes in brain rsFC in school-aged children. Procedural learning was assessed in 30 typically developing children (mean age ± SD: 9.99 years ± 1.35) using a serial reaction time task (SRTT). During SRTT, participants touched as quickly and accurately as possible a stimulus sequentially or randomly appearing in one of the quadrants of a touchscreen. Band-limited power envelope correlation (brain rsFC) was applied to MEG data acquired at rest pre- and post-learning. Correlation analyses were performed between brain rsFC and sequence-specific learning or response time indices. Stronger pre-learning interhemispheric rsFC between inferior parietal and primary somatosensory/motor areas correlated with better subsequent sequence learning performance and faster visuomotor response time. Faster response time was associated with post-learning decreased rsFC within the dorsal extra-striate visual stream and increased rsFC between temporo-cerebellar regions. In school-aged children, variations in functional brain architecture at rest within the sensorimotor network account for interindividual differences in sequence learning and visuomotor performance. After learning, rapid adjustments in functional brain architecture are associated with visuomotor performance but not sequence learning skills.

Rapidly Decoding Image Categories From MEG Data Using a Multivariate Short-Time FC Pattern Analysis Approach

Recent advances in the development of multivariate analysis methods have led to the application of multivariate pattern analysis (MVPA) to investigate the interactions between brain regions using graph theory (functional connectivity, FC) and decode visual categories from functional magnetic resonance imaging (fMRI) data from a continuous multicategory paradigm. To estimate stable FC patterns from fMRI data, previous studies required long periods in the order of several minutes, in comparison to the human brain that categories visual stimuli within hundreds of milliseconds. Constructing short-time dynamic FC patterns in the order of milliseconds and decoding visual categories is a relatively novel concept. In this study, we developed a multivariate decoding algorithm based on FC patterns and applied it to magnetoencephalography (MEG) data. MEG data were recorded from participants presented with image stimuli in four categories (faces, scenes, animals and tools). MEG data from 17 participants demonstrate that short-time dynamic FC patterns yield brain activity patterns that can be used to decode visual categories with high accuracy. Our results show that FC patterns change over the time window, and FC patterns extracted in the time window of 0∼200 ms after the stimulus onset were most stable. Further, the categorizing accuracy peaked (the mean binary accuracy is above 78.6% at individual level) in the FC patterns estimated within the 0∼200 ms interval. These findings elucidate the underlying connectivity information during visual category processing on a relatively smaller time scale and demonstrate that the contribution of FC patterns to categorization fluctuates over time.

Measuring the cortical tracking of speech with optically-pumped magnetometers

During continuous speech listening, brain activity tracks speech rhythmicity at frequencies matching with the repetition rate of phrases (0.2-1.5 Hz), words (2-4 Hz) and syllables (4-8 Hz). Here, we evaluated the applicability of wearable MEG based on optically-pumped magnetometers (OPMs) to measure such cortical tracking of speech (CTS). Measuring CTS with OPMs is a priori challenging given the complications associated with OPM measurements at frequencies below 4 Hz, due to increased intrinsic interference and head movement artifacts. Still, this represents an important development as OPM-MEG provides lifespan compliance and substantially improved spatial resolution compared with classical MEG. In this study, four healthy right-handed adults listened to continuous speech for 9 min. The radial component of the magnetic field was recorded simultaneously with 45-46 OPMs evenly covering the scalp surface and fixed to an additively manufactured helmet which fitted all 4 participants. We estimated CTS with reconstruction accuracy and coherence, and determined the number of dominant principal components (PCs) to remove from the data (as a preprocessing step) for optimal estimation. We also identified the dominant source of CTS using a minimum norm estimate. CTS estimated with reconstruction accuracy and coherence was significant in all 4 participants at phrasal and word rates, and in 3 participants (reconstruction accuracy) or 2 (coherence) at syllabic rate. Overall, close-to-optimal CTS estimation was obtained when the 3 (reconstruction accuracy) or 10 (coherence) first PCs were removed from the data. Importantly, values of reconstruction accuracy (~0.4 for 0.2-1.5-Hz CTS and ~0.1 for 2-8-Hz CTS) were remarkably close to those previously reported in classical MEG studies. Finally, source reconstruction localized the main sources of CTS to bilateral auditory cortices. In conclusion, t his study demonstrates that OPMs can be used for the purpose of CTS assessment. This finding opens new research avenues to unravel the neural network involved in CTS across the lifespan and potential alterations in, e.g., language developmental disorders. Data also suggest that OPMs are generally suitable for recording neural activity at frequencies below 4 Hz provided PCA is used as a preprocessing step; 0.2-1.5-Hz being the lowest frequency range successfully investigated here.

BCI learning induces core-periphery reorganization in M/EEG multiplex brain networks

Brain-computer interfaces (BCIs) constitute a promising tool for communication and control. However, mastering non-invasive closed-loop systems remains a learned skill that is difficult to develop for a non-negligible proportion of users. The involved learning process induces neural changes associated with a brain network reorganization that remains poorly understood. To address this inter-subject variability, we adopted a multilayer approach to integrate brain network properties from electroencephalographic (EEG) and magnetoencephalographic (MEG) data resulting from a four-session BCI training program followed by a group of healthy subjects. Our method gives access to the contribution of each layer to multilayer network that tends to be equal with time. We show that regardless the chosen modality, a progressive increase in the integration of somatosensory areas in the α band was paralleled by a decrease of the integration of visual processing and working memory areas in the β band. Notably, only brain network properties in multilayer network correlated with future BCI scores in the α2 band: positively in somatosensory and decision-making related areas and negatively in associative areas. Our findings cast new light on neural processes underlying BCI training. Integrating multimodal brain network properties provides new information that correlates with behavioral performance and could be considered as a potential marker of BCI learning.

The relationship between EEG and fMRI connectomes is reproducible across simultaneous EEG-fMRI studies from 1.5T to 7T

Both electroencephalography (EEG) and functional Magnetic Resonance Imaging (fMRI) are non-invasive methods that show complementary aspects of human brain activity. Despite measuring different proxies of brain activity, both the measured blood-oxygenation (fMRI) and neurophysiological recordings (EEG) are indirectly coupled. The electrophysiological and BOLD signal can map the underlying functional connectivity structure at the whole brain scale at different timescales. Previous work demonstrated a moderate but significant correlation between resting-state functional connectivity of both modalities, however there is a wide range of technical setups to measure simultaneous EEG-fMRI and the reliability of those measures between different setups remains unknown. This is true notably with respect to different magnetic field strengths (low and high field) and different spatial sampling of EEG (medium to high-density electrode coverage). Here, we investigated the reproducibility of the bimodal EEG-fMRI functional connectome in the most comprehensive resting-state simultaneous EEG-fMRI dataset compiled to date including a total of 72 subjects from four different imaging centers. Data was acquired from 1.5T, 3T and 7T scanners with simultaneously recorded EEG using 64 or 256 electrodes. We demonstrate that the whole-brain monomodal connectivity reproducibly correlates across different datasets and that a moderate crossmodal correlation between EEG and fMRI connectivity of r ≈ 0.3 can be reproducibly extracted in low- and high-field scanners. The crossmodal correlation was strongest in the EEG-β frequency band but exists across all frequency bands. Both homotopic and within intrinsic connectivity network (ICN) connections contributed the most to the crossmodal relationship. This study confirms, using a considerably diverse range of recording setups, that simultaneous EEG-fMRI offers a consistent estimate of multimodal functional connectomes in healthy subjects that are dominantly linked through a functional core of ICNs across spanning across the different timescales measured by EEG and fMRI. This opens new avenues for estimating the dynamics of brain function and provides a better understanding of interactions between EEG and fMRI measures. This observed level of reproducibility also defines a baseline for the study of alterations of this coupling in pathological conditions and their role as potential clinical markers.

Measuring functional connectivity with wearable MEG

Optically-pumped magnetometers (OPMs) offer the potential for a step change in magnetoencephalography (MEG) enabling wearable systems that provide improved data quality, accommodate any subject group, allow data capture during movement and potentially reduce cost. However, OPM-MEG is a nascent technology and, to realise its potential, it must be shown to facilitate key neuroscientific measurements, such as the characterisation of brain networks. Networks, and the connectivities that underlie them, have become a core area of neuroscientific investigation, and their importance is underscored by many demonstrations of their disruption in brain disorders. Consequently, a demonstration of network measurements using OPM-MEG would be a significant step forward. Here, we aimed to show that a wearable 50-channel OPM-MEG system enables characterisation of the electrophysiological connectome. To this end, we measured connectivity in the resting state and during a visuo-motor task, using both OPM-MEG and a state-of-the-art 275-channel cryogenic MEG device. Our results show that resting-state connectome matrices from OPM and cryogenic systems exhibit a high degree of similarity, with correlation values >70%. In addition, in task data, similar differences in connectivity between individuals (scanned multiple times) were observed in cryogenic and OPM-MEG data, again demonstrating the fidelity of the OPM-MEG device. This is the first demonstration of network connectivity measured using OPM-MEG, and results add weight to the argument that OPMs will ultimately supersede cryogenic sensors for MEG measurement.

The dynamic modular fingerprints of the human brain at rest

The human brain is a dynamic modular network that can be decomposed into a set of modules, and its activity changes continually over time. At rest, several brain networks, known as Resting-State Networks (RSNs), emerge and cross-communicate even at sub-second temporal scale. Here, we seek to decipher the fast reshaping in spontaneous brain modularity and its relationships with RSNs. We use Electro/Magneto-Encephalography (EEG/MEG) to track the dynamics of modular brain networks, in three independent datasets (N = 568) of healthy subjects at rest. We show the presence of strikingly consistent RSNs, and a splitting phenomenon of some of these networks, especially the default mode network, visual, temporal and dorsal attentional networks. We also demonstrate that between-subjects variability in mental imagery is associated with the temporal characteristics of specific modules, particularly the visual network. Taken together, our findings show that large-scale electrophysiological networks have modularity-dependent dynamic fingerprints at rest.

Alterations in resting-state network dynamics along the Alzheimer's disease continuum

Human brain activity is intrinsically organized into resting-state networks (RSNs) that transiently activate or deactivate at the sub-second timescale. Few neuroimaging studies have addressed how Alzheimer's disease (AD) affects these fast temporal brain dynamics, and how they relate to the cognitive, structural and metabolic abnormalities characterizing AD. We aimed at closing this gap by investigating both brain structure and function using magnetoencephalography (MEG) and hybrid positron emission tomography-magnetic resonance (PET/MR) in 10 healthy elders, 10 patients with subjective cognitive decline (SCD), 10 patients with amnestic mild cognitive impairment (aMCI) and 10 patients with typical Alzheimer's disease with dementia (AD). The fast activation/deactivation state dynamics of RSNs were assessed using hidden Markov modeling (HMM) of power envelope fluctuations at rest measured with MEG. Correlations were sought between temporal properties of HMM states and participants' cognitive test scores, whole hippocampal grey matter volume and regional brain glucose metabolism. The posterior default-mode network (DMN) was less often activated and for shorter durations in AD patients than matched healthy elders. No significant difference was found in patients with SCD or aMCI. The time spent by participants in the activated posterior DMN state did not correlate significantly with cognitive scores, nor with the whole hippocampal volume. However, it correlated positively with the regional glucose consumption in the right dorsolateral prefrontal cortex (DLPFC). AD patients present alterations of posterior DMN power activation dynamics at rest that identify an additional electrophysiological correlate of AD-related synaptic and neural dysfunction. The right DLPFC may play a causal role in the activation of the posterior DMN, possibly linked to the occurrence of mind wandering episodes. As such, these data might suggest a neural correlate of the decrease in mind wandering episodes reported in pathological aging.

Brain dysconnectivity relates to disability and cognitive impairment in multiple sclerosis

The pathophysiology of cognitive dysfunction in multiple sclerosis (MS) is still unclear. This magnetoencephalography (MEG) study investigates the impact of MS on brain resting-state functional connectivity (rsFC) and its relationship to disability and cognitive impairment. We investigated rsFC based on power envelope correlation within and between different frequency bands, in a large cohort of participants consisting of 99 MS patients and 47 healthy subjects. Correlations were investigated between rsFC and outcomes on disability, disease duration and 7 neuropsychological scores within each group, while stringently correcting for multiple comparisons and possible confounding factors. Specific dysconnections correlating with MS-induced physical disability and disease duration were found within the sensorimotor and language networks, respectively. Global network-level reductions in within- and cross-network rsFC were observed in the default-mode network. Healthy subjects and patients significantly differed in their scores on cognitive fatigue and verbal fluency. Healthy subjects and patients showed different correlation patterns between rsFC and cognitive fatigue or verbal fluency, both of which involved a shift in patients from the posterior default-mode network to the language network. Introducing electrophysiological rsFC in a regression model of verbal fluency and cognitive fatigue in MS patients significantly increased the explained variance compared to a regression limited to structural MRI markers (relative thalamic volume and lesion load). This MEG study demonstrates that MS induces distinct changes in the resting-state functional brain architecture that relate to disability, disease duration and specific cognitive functioning alterations. It highlights the potential value of electrophysiological intrinsic rsFC for monitoring the cognitive impairment in patients with MS.

Changes in electrophysiological static and dynamic human brain functional architecture from childhood to late adulthood

This magnetoencephalography study aimed at characterizing age-related changes in resting-state functional brain organization from mid-childhood to late adulthood. We investigated neuromagnetic brain activity at rest in 105 participants divided into three age groups: children (6-9 years), young adults (18-34 years) and healthy elders (53-78 years). The effects of age on static resting-state functional brain integration were assessed using band-limited power envelope correlation, whereas those on transient functional brain dynamics were disclosed using hidden Markov modeling of power envelope activity. Brain development from childhood to adulthood came with (1) a strengthening of functional integration within and between resting-state networks and (2) an increased temporal stability of transient (100-300 ms lifetime) and recurrent states of network activation or deactivation mainly encompassing lateral or medial associative neocortical areas. Healthy aging was characterized by decreased static resting-state functional integration and dynamic stability within the primary visual network. These results based on electrophysiological measurements free of neurovascular biases suggest that functional brain integration mainly evolves during brain development, with limited changes in healthy aging. These novel electrophysiological insights into human brain functional architecture across the lifespan pave the way for future clinical studies investigating how brain disorders affect brain development or healthy aging.

Consistency of local activation parameters at sensor- and source-level in neural signals

OBJECTIVE

Although magnetoencephalography and electroencephalography (M/EEG) signals at sensor level are robust and reliable, they suffer from different degrees of distortion due to changes in brain tissue conductivities, known as field spread and volume conduction effects. To estimate original neural generators from M/EEG activity acquired at sensor level, diverse source localisation algorithms have been proposed; however, they are not exempt from limitations and usually involve time-consuming procedures. Connectivity and network-based M/EEG analyses have been found to be affected by field spread and volume conduction effects; nevertheless, the influence of the aforementioned effects on widely used local activation parameters has not been assessed yet. The goal of this study is to evaluate the consistency of various local activation parameters when they are computed at sensor- and source-level.

APPROACH

Six spectral (relative power, median frequency, and individual alpha frequency) and non-linear parameters (Lempel-Ziv complexity, sample entropy, and central tendency measure) are computed from M/EEG signals at sensor- and source-level using four source inversion methods: weighted minimum norm estimate (wMNE), standardised low-resolution brain electromagnetic tomography (sLORETA), linear constrained minimum variance (LCMV), and dynamical statistical parametric mapping (dSPM).

MAIN RESULTS

Our results show that the spectral and non-linear parameters yield similar results at sensor- and source-level, showing high correlation values between them for all the source inversion methods evaluated and both modalities of signal, EEG and MEG. Furthermore, the correlation values remain high when performing coarse-grained spatial analyses.

SIGNIFICANCE

To the best of our knowledge, this is the first study analysing how field spread and volume conduction effects impact on local activation parameters computed from resting-state neural activity. Our findings evidence that local activation parameters are robust against field spread and volume conduction effects and provide equivalent information at sensor- and source-level even when performing regional analyses.

Prediction and classification of sleep quality based on phase synchronization related whole-brain dynamic connectivity using resting state fMRI

Recently, functional network connectivity (FNC) has been extended from static to dynamic analysis to explore the time-varying functional organization of brain networks. Nowadays, a majority of dynamic FNC (dFNC) analysis frameworks identified recurring FNC patterns with linear correlations based on the amplitude of fMRI time series. However, the brain is a complex dynamical system and phase synchronization provides more informative measures. This paper proposes a novel framework for the prediction/classification of behaviors and cognitions based on the dFNCs derived from phase locking value. When applying to the analysis of fMRI data from Human Connectome Project (HCP), four dFNC states are identified for the study of sleep quality. State 1 exhibits most intense phase synchronization across the whole brain. States 2 and 3 have low and weak connections, respectively. State 4 exhibits strong phase synchronization in intra and inter-connections of somatomotor, visual and cognitive control networks. Through the two-sample t-test, we reveal that for the group with bad sleep quality, state 4 shows decreased phase synchronization within and between networks such as subcortical, auditory, somatomotor and visual, but increased phase synchronization within cognitive control network, and between this network and somatomotor/visual/default-mode/cerebellar networks. The networks with increased phase synchronization in state 4 behave oppositely in state 2. Group differences are absent in state 3, and weak in state 1. We establish a prediction model by linear regression of FNC against sleep quality, and adopt a support vector machine approach for the classification. We compare the performance between conventional FNC and PLV-based dFNC with cross-validation. Results show that the PLV-based dFNC significantly outperforms the conventional FNC in terms of both predictive power and classification accuracy. We also observe that combining static and dynamic features does not significantly improve the classification over using dFNC features alone. Overall, the proposed approach provides a novel means to assess dFNC, which can be used as brain fingerprints to facilitate prediction and classification.

Frequency-Dependent Intrinsic Electrophysiological Functional Architecture of the Human Verbal Language Network

Functional magnetic resonance imaging (fMRI) allowed the spatial characterization of the resting-state verbal language network (vLN). While other resting-state networks (RSNs) were matched with their electrophysiological equivalents at rest and could be spectrally defined, such correspondence is lacking for the vLN. This magnetoencephalography (MEG) study aimed at defining the spatio-spectral characteristics of the neuromagnetic intrinsic functional architecture of the vLN. Neuromagnetic activity was recorded at rest in 100 right-handed healthy adults (age range: 18-41 years). Band-limited power envelope correlations were performed within and across frequency bands (θ, α, β, and low γ) from a seed region placed in the left Broca's area, using static orthogonalization as leakage correction. K-means clustering was used to segregate spatio-spectral clusters of resting-state functional connectivity (rsFC). Remarkably, unlike other RSNs, within-frequency long-range rsFC from the left Broca's area was not driven by one main carrying frequency but was characterized by a specific spatio-spectral pattern segregated along the ventral (predominantly θ and α) and dorsal (β and low-γ bands) vLN streams. In contrast, spatial patterns of cross-frequency vLN functional integration were spectrally more widespread and involved multiple frequency bands. Moreover, the static intrinsic functional architecture of the neuromagnetic human vLN involved clearly left-hemisphere-dominant vLN interactions as well as cross-network interactions with the executive control network and postero-medial nodes of the DMN. Overall, this study highlighted the involvement of multiple modes of within and cross-frequency power envelope couplings at the basis of long-range electrophysiological vLN functional integration. As such, it lays the foundation for future works aimed at understanding the pathophysiology of language-related disorders.

Multi-channel whole-head OPM-MEG: Helmet design and a comparison with a conventional system

Magnetoencephalography (MEG) is a powerful technique for functional neuroimaging, offering a non-invasive window on brain electrophysiology. MEG systems have traditionally been based on cryogenic sensors which detect the small extracranial magnetic fields generated by synchronised current in neuronal assemblies, however, such systems have fundamental limitations. In recent years, non-cryogenic quantum-enabled sensors, called optically-pumped magnetometers (OPMs), in combination with novel techniques for accurate background magnetic field control, have promised to lift those restrictions offering an adaptable, motion-robust MEG system, with improved data quality, at reduced cost. However, OPM-MEG remains a nascent technology, and whilst viable systems exist, most employ small numbers of sensors sited above targeted brain regions. Here, building on previous work, we construct a wearable OPM-MEG system with ‘whole-head’ coverage based upon commercially available OPMs, and test its capabilities to measure alpha, beta and gamma oscillations. We design two methods for OPM mounting; a flexible (EEG-like) cap and rigid (additively-manufactured) helmet. Whilst both designs allow for high quality data to be collected, we argue that the rigid helmet offers a more robust option with significant advantages for reconstruction of field data into 3D images of changes in neuronal current. Using repeat measurements in two participants, we show signal detection for our device to be highly robust. Moreover, via application of source-space modelling, we show that, despite having 5 times fewer sensors, our system exhibits comparable performance to an established cryogenic MEG device. While significant challenges still remain, these developments provide further evidence that OPM-MEG is likely to facilitate a step change for functional neuroimaging.