Sensitivity of MEG and EEG to source orientation

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.

Identifying epileptogenic abnormality by decomposing intracranial EEG and MEG power spectra

BACKGROUND

Accurate identification of abnormal electroencephalographic (EEG) activity is pivotal for diagnosing and treating epilepsy. Recent studies indicate that decomposing brain activity into periodic (oscillatory) and aperiodic (trend across all frequencies) components can illuminate the drivers of spectral activity changes.

NEW METHODS

We analysed intracranial EEG (iEEG) data from 234 subjects, creating a normative map. This map was compared to a cohort of 63 patients with refractory focal epilepsy under consideration for neurosurgery. The normative map was computed using three approaches: (i) relative complete band power, (ii) relative band power with the aperiodic component removed, and (iii) the aperiodic exponent. Abnormalities were calculated for each approach in the patient cohort. We evaluated the spatial profiles, assessed their ability to localize abnormalities, and replicated the findings using magnetoencephalography (MEG).

RESULTS

Normative maps of relative complete band power and relative periodic band power exhibited similar spatial profiles, while the aperiodic normative map revealed higher exponent values in the temporal lobe. Abnormalities estimated through complete band power effectively distinguished between good and bad outcome patients. Combining periodic and aperiodic abnormalities enhanced performance, like the complete band power approach.

COMPARISON WITH EXISTING METHODS AND CONCLUSIONS

Sparing cerebral tissue with abnormalities in both periodic and aperiodic activity may result in poor surgical outcomes. Both periodic and aperiodic components do not carry sufficient information in isolation. The relative complete band power solution proved to be the most reliable method for this purpose. Future studies could investigate how cerebral location or pathology influences periodic or aperiodic abnormalities.

Learning Spatiotemporal Brain Dynamics in Adolescents via Multimodal MEG and fMRI Data Fusion Using Joint Tensor/Matrix Decomposition

OBJECTIVE

Brain function is understood to be regulated by complex spatiotemporal dynamics, and can be characterized by a combination of observed brain response patterns in time and space. Magnetoencephalography (MEG), with its high temporal resolution, and functional magnetic resonance imaging (fMRI), with its high spatial resolution, are complementary imaging techniques with great potential to reveal information about spatiotemporal brain dynamics. Hence, the complementary nature of these imaging techniques holds much promise to study brain function in time and space, especially when the two data types are allowed to fully interact.

METHODS

We employed coupled tensor/matrix factorization (CMTF) to extract joint latent components in the form of unique spatiotemporal brain patterns that can be used to study brain development and function on a millisecond scale.

RESULTS

Using the CMTF model, we extracted distinct brain patterns that revealed fine-grained spatiotemporal brain dynamics and typical sensory processing pathways informative of high-level cognitive functions in healthy adolescents. The components extracted from multimodal tensor fusion possessed better discriminative ability between high- and low-performance subjects than single-modality data-driven models.

CONCLUSION

Multimodal tensor fusion successfully identified spatiotemporal brain dynamics of brain function and produced unique components with high discriminatory power.

SIGNIFICANCE

The CMTF model is a promising tool for high-order, multimodal data fusion that exploits the functional resolution of MEG and fMRI, and provides a comprehensive picture of the developing brain in time and space.

Single neuron contributions to the auditory brainstem EEG

The auditory brainstem response (ABR) is an acoustically evoked EEG potential that is an important diagnostic tool for hearing loss, especially in newborns. The ABR originates from the response sequence of auditory brainstem nuclei, and a click-evoked ABR typically shows three positive peaks ('waves') within the first six milliseconds. However, an assignment of the waves of the ABR to specific sources is difficult, and a quantification of contributions to the ABR waves is not available. Here, we exploit the large size and physical separation of the barn owl first-order cochlear nucleus magnocellularis (NM) to estimate single-cell contributions to the ABR. We simultaneously recorded NM neurons' spikes and the EEG, and found that > 5,000 spontaneous single-cell spikes are necessary to isolate a significant spike-triggered average response at the EEG electrode. An average single-neuron contribution to the ABR was predicted by convolving the spike-triggered average with the cell's peri-stimulus time histogram. Amplitudes of predicted contributions of single NM cells typically reached 32.9 nV (mean, range: 2.5 - 162.7 nV), or 0.07% (median, range: 0.01 - 4.0%) of the ABR amplitude. The time of the predicted peak coincided best with the peak of the ABR wave II, and this coincidence was independent of the click sound level. Our results suggest that wave II of the ABR is shaped by a small fraction of NM units.

Simultaneous MEG and EEG source imaging of electrophysiological activity in response to acute transcranial photobiomodulation

Introduction

Transcranial photobiomodulation (tPBM) is a non-invasive neuromodulation technique that improves human cognition. The effects of tPBM of the right forehead on neurophysiological activity have been previously investigated using EEG in sensor space. However, the spatial resolution of these studies is limited. Magnetoencephalography (MEG) is known to facilitate a higher spatial resolution of brain source images. This study aimed to image post-tPBM effects in brain space based on both MEG and EEG measurements across the entire human brain.

Methods

MEG and EEG scans were concurrently acquired for 6 min before and after 8-min of tPBM delivered using a 1,064-nm laser on the right forehead of 25 healthy participants. Group-level changes in both the MEG and EEG power spectral density with respect to the baseline (pre-tPBM) were quantified and averaged within each frequency band in the sensor space. Constrained modeling was used to generate MEG and EEG source images of post-tPBM, followed by cluster-based permutation analysis for family wise error correction (p < 0.05).

Results

The 8-min tPBM enabled significant increases in alpha (8-12 Hz) and beta (13-30 Hz) powers across multiple cortical regions, as confirmed by MEG and EEG source images. Moreover, tPBM-enhanced oscillations in the beta band were located not only near the stimulation site but also in remote cerebral regions, including the frontal, parietal, and occipital regions, particularly on the ipsilateral side.

Discussion

MEG and EEG results shown in this study demonstrated that tPBM modulates neurophysiological activity locally and in distant cortical areas. The EEG topographies reported in this study were consistent with previous observations. This study is the first to present MEG and EEG evidence of the electrophysiological effects of tPBM in the brain space, supporting the potential utility of tPBM in treating neurological diseases through the modulation of brain oscillations.

A Multivariate analysis on evoked components of Chinese semantic congruity: an OP-MEG study with EEG

The application of wearable magnetoencephalography using optically-pumped magnetometers has drawn extensive attention in the field of neuroscience. Electroencephalogram system can cover the whole head and reflect the overall activity of a large number of neurons. The efficacy of optically-pumped magnetometer in detecting event-related components can be validated through electroencephalogram results. Multivariate pattern analysis is capable of tracking the evolution of neurocognitive processes over time. In this paper, we adopted a classical Chinese semantic congruity paradigm and separately collected electroencephalogram and optically-pumped magnetometer signals. Then, we verified the consistency of optically-pumped magnetometer and electroencephalogram in detecting N400 using mutual information index. Multivariate pattern analysis revealed the difference in decoding performance of these two modalities, which can be further validated by dynamic/stable coding analysis on the temporal generalization matrix. The results from searchlight analysis provided a neural basis for this dissimilarity at the magnetoencephalography source level and the electroencephalogram sensor level. This study opens a new avenue for investigating the brain's coding patterns using wearable magnetoencephalography and reveals the differences in sensitivity between the two modalities in reflecting neuron representation patterns.

Comparing the performance of beamformer algorithms in estimating orientations of neural sources

The efficacy of transcranial electric stimulation (tES) to effectively modulate neuronal activity depends critically on the spatial orientation of the targeted neuronal population. Therefore, precise estimation of target orientation is of utmost importance. Different beamforming algorithms provide orientation estimates; however, a systematic analysis of their performance is still lacking. For fixed brain locations, EEG and MEG data from sources with randomized orientations were simulated. The orientation was then estimated (1) with an EEG and (2) with a combined EEG-MEG approach. Three commonly used beamformer algorithms were evaluated with respect to their abilities to estimate the correct orientation: Unit-Gain (UG), Unit-Noise-Gain (UNG), and Array-Gain (AG) beamformer. Performance depends on the signal-to-noise ratios for the modalities and on the chosen beamformer. Overall, the UNG and AG beamformers appear as the most reliable. With increasing noise, the UG estimate converges to a vector determined by the leadfield, thus leading to insufficient orientation estimates.

The neurophysiology of closed-loop auditory stimulation in sleep: A magnetoencephalography study

Closed-loop auditory stimulation (CLAS) is a brain modulation technique in which sounds are timed to enhance or disrupt endogenous neurophysiological events. CLAS of slow oscillation up-states in sleep is becoming a popular tool to study and enhance sleep's functions, as it increases slow oscillations, evokes sleep spindles and enhances memory consolidation of certain tasks. However, few studies have examined the specific neurophysiological mechanisms involved in CLAS, in part because of practical limitations to available tools. To evaluate evidence for possible models of how sound stimulation during brain up-states alters brain activity, we simultaneously recorded electro- and magnetoencephalography in human participants who received auditory stimulation across sleep stages. We conducted a series of analyses that test different models of pathways through which CLAS of slow oscillations may affect widespread neural activity that have been suggested in literature, using spatial information, timing and phase relationships in the source-localized magnetoencephalography data. The results suggest that auditory information reaches ventral frontal lobe areas via non-lemniscal pathways. From there, a slow oscillation is created and propagated. We demonstrate that while the state of excitability of tissue in auditory cortex and frontal ventral regions shows some synchrony with the electroencephalography (EEG)-recorded up-states that are commonly used for CLAS, it is the state of ventral frontal regions that is most critical for slow oscillation generation. Our findings advance models of how CLAS leads to enhancement of slow oscillations, sleep spindles and associated cognitive benefits and offer insight into how the effectiveness of brain stimulation techniques can be improved.

Solving large-scale MEG/EEG source localisation and functional connectivity problems simultaneously using state-space models

State-space models are widely employed across various research disciplines to study unobserved dynamics. Conventional estimation techniques, such as Kalman filtering and expectation maximisation, offer valuable insights but incur high computational costs in large-scale analyses. Sparse inverse covariance estimators can mitigate these costs, but at the expense of a trade-off between enforced sparsity and increased estimation bias, necessitating careful assessment in low signal-to-noise ratio (SNR) situations. To address these challenges, we propose a three-fold solution: (1) Introducing multiple penalised state-space (MPSS) models that leverage data-driven regularisation; (2) Developing novel algorithms derived from backpropagation, gradient descent, and alternating least squares to solve MPSS models; (3) Presenting a K-fold cross-validation extension for evaluating regularisation parameters. We validate this MPSS regularisation framework through lower and more complex simulations under varying SNR conditions, including a large-scale synthetic magneto- and electro-encephalography (MEG/EEG) data analysis. In addition, we apply MPSS models to concurrently solve brain source localisation and functional connectivity problems for real event-related MEG/EEG data, encompassing thousands of sources on the cortical surface. The proposed methodology overcomes the limitations of existing approaches, such as constraints to small-scale and region-of-interest analyses. Thus, it may enable a more accurate and detailed exploration of cognitive brain functions.

Cleaning up the Brickyard: How Theory and Methodology Shape Experiments in Cognitive Neuroscience of Language

The capacity for language is a defining property of our species, yet despite decades of research, evidence on its neural basis is still mixed and a generalized consensus is difficult to achieve. We suggest that this is partly caused by researchers defining "language" in different ways, with focus on a wide range of phenomena, properties, and levels of investigation. Accordingly, there is very little agreement among cognitive neuroscientists of language on the operationalization of fundamental concepts to be investigated in neuroscientific experiments. Here, we review chains of derivation in the cognitive neuroscience of language, focusing on how the hypothesis under consideration is defined by a combination of theoretical and methodological assumptions. We first attempt to disentangle the complex relationship between linguistics, psychology, and neuroscience in the field. Next, we focus on how conclusions that can be drawn from any experiment are inherently constrained by auxiliary assumptions, both theoretical and methodological, on which the validity of conclusions drawn rests. These issues are discussed in the context of classical experimental manipulations as well as study designs that employ novel approaches such as naturalistic stimuli and computational modeling. We conclude by proposing that a highly interdisciplinary field such as the cognitive neuroscience of language requires researchers to form explicit statements concerning the theoretical definitions, methodological choices, and other constraining factors involved in their work.

Both stronger and weaker cerebro-cerebellar functional connectivity patterns during processing of spoken sentences in autism spectrum disorder

Cerebellar differences have long been documented in autism spectrum disorder (ASD), yet the extent to which such differences might impact language processing in ASD remains unknown. To investigate this, we recorded brain activity with magnetoencephalography (MEG) while ASD and age-matched typically developing (TD) children passively processed spoken meaningful English and meaningless Jabberwocky sentences. Using a novel source localization approach that allows higher resolution MEG source localization of cerebellar activity, we found that, unlike TD children, ASD children showed no difference between evoked responses to meaningful versus meaningless sentences in right cerebellar lobule VI. ASD children also had atypically weak functional connectivity in the meaningful versus meaningless speech condition between right cerebellar lobule VI and several left-hemisphere sensorimotor and language regions in later time windows. In contrast, ASD children had atypically strong functional connectivity for in the meaningful versus meaningless speech condition between right cerebellar lobule VI and primary auditory cortical areas in an earlier time window. The atypical functional connectivity patterns in ASD correlated with ASD severity and the ability to inhibit involuntary attention. These findings align with a model where cerebro-cerebellar speech processing mechanisms in ASD are impacted by aberrant stimulus-driven attention, which could result from atypical temporal information and predictions of auditory sensory events by right cerebellar lobule VI.

Multi-objective optimization via evolutionary algorithm (MOVEA) for high-definition transcranial electrical stimulation of the human brain

Designing a transcranial electrical stimulation (tES) strategy requires considering multiple objectives, such as intensity in the target area, focality, stimulation depth, and avoidance zone. These objectives are often mutually exclusive. In this paper, we propose a general framework, called multi-objective optimization via evolutionary algorithm (MOVEA), which solves the non-convex optimization problem in designing tES strategies without a predefined direction. MOVEA enables simultaneous optimization of multiple targets through Pareto optimization, generating a Pareto front after a single run without manual weight adjustment and allowing easy expansion to more targets. This Pareto front consists of optimal solutions that meet various requirements while respecting trade-off relationships between conflicting objectives such as intensity and focality. MOVEA is versatile and suitable for both transcranial alternating current stimulation (tACS) and transcranial temporal interference stimulation (tTIS) based on high definition (HD) and two-pair systems. We comprehensively compared tACS and tTIS in terms of intensity, focality, and steerability for targets at different depths. Our findings reveal that tTIS enhances focality by reducing activated volume outside the target by 60%. HD-tTIS and HD-tDCS can achieve equivalent maximum intensities, surpassing those of two-pair tTIS, such as 0.51 V/m under HD-tACS/HD-tTIS and 0.42 V/m under two-pair tTIS for the motor area as a target. Analysis of variance in eight subjects highlights individual differences in both optimal stimulation policies and outcomes for tACS and tTIS, emphasizing the need for personalized stimulation protocols. These findings provide guidance for designing appropriate stimulation strategies for tACS and tTIS. MOVEA facilitates the optimization of tES based on specific objectives and constraints, advancing tTIS and tACS-based neuromodulation in understanding the causal relationship between brain regions and cognitive functions and treating diseases. The code for MOVEA is available at https://github.com/ncclabsustech/MOVEA.

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.

A high-density EEG and structural MRI source analysis of the frequency following response to missing fundamental stimuli reveals subcortical and cortical activation to low and high frequency stimuli

Pitch is a perceptual rather than physical phenomenon, important for spoken language use, musical communication, and other aspects of everyday life. Auditory stimuli can be designed to probe the relationship between perception and physiological responses to pitch-evoking stimuli. One technique for measuring physiological responses to pitch-evoking stimuli is the frequency following response (FFR). The FFR is an electroencephalographic (EEG) response to periodic auditory stimuli. The FFR contains nonlinearities not present in the stimuli, including correlates of the amplitude envelope of the stimulus; however, these nonlinearities remain undercharacterized. The FFR is a composite response reflecting multiple neural and peripheral generators, and their contributions to the scalp-recorded FFR vary in ill-understood ways depending on the electrode montage, stimulus, and imaging technique. The FFR is typically assumed to be generated in the auditory brainstem; there is also evidence both for and against a cortical contribution to the FFR. Here a methodology is used to examine the FFR correlates of pitch and the generators of the FFR to stimuli with different pitches. Stimuli were designed to tease apart biological correlates of pitch and amplitude envelope. FFRs were recorded with 256-electrode EEG nets, in contrast to a typical FFR setup which only contains a single active electrode. Structural MRI scans were obtained for each participant to co-register with the electrode locations and constrain a source localization algorithm. The results of this localization shed light on the generating mechanisms of the FFR, including providing evidence for both cortical and subcortical auditory sources.

A neural marker of the human face identity familiarity effect

Human adults associate different views of an identity much better for familiar than for unfamiliar faces. However, a robust and consistent neural index of this behavioral face identity familiarity effect (FIFE)-not found in non-human primate species-is lacking. Here we provide such a neural FIFE index, measured implicitly and with one fixation per face. Fourteen participants viewed 70 s stimulation sequences of a large set (n = 40) of widely variable natural images of a face identity at a rate of 6 images/second (6 Hz). Different face identities appeared every 5th image (1.2 Hz). In a sequence, face images were either familiar (i.e., famous) or unfamiliar, participants performing a non-periodic task unrelated to face recognition. The face identity recognition response identified at 1.2 Hz over occipital-temporal regions in the frequency-domain electroencephalogram was 3.4 times larger for familiar than unfamiliar faces. The neural response to familiar faces-which emerged at about 180 ms following face onset-was significant in each individual but a case of prosopdysgnosia. Besides potential clinical and forensic applications to implicitly measure one's knowledge of a face identity, these findings open new perspectives to clarify the neurofunctional source of the FIFE and understand the nature of human face identity recognition.

Electromagnetic source imaging predicts surgical outcome in children with focal cortical dysplasia

OBJECTIVE

To evaluate the diagnostic accuracy of electromagnetic source imaging (EMSI) in localizing spikes and predict surgical outcome in children with drug resistant epilepsy (DRE) due to focal cortical dysplasia (FCD).

METHODS

We retrospectively analyzed magnetoencephalography (MEG) and high-density (HD-EEG) data from 23 children with FCD-associated DRE who underwent intracranial EEG and surgery. We localized spikes using equivalent current dipole (ECD) fitting, dipole clustering, and dynamical statistical parametric mapping (dSPM) on EMSI, electric source imaging (ESI), and magnetic source imaging (MSI). We calculated the distance from the seizure onset zone (DSOZ) and resection (DRES). We estimated receiver operating characteristic (ROC) curves with Youden's index (J) to predict outcome.

RESULTS

EMSI presented shorter DSOZ (15.18 ± 9.06 mm) and DRES (8.56 ± 6.24 mm) compared to ESI (DSOZ: 25.04 ± 16.20 mm, p < 0.009; DRES: 18.88 ± 17.30 mm, p < 0.03) and MSI (DSOZ: 23.37 ± 8.98 mm, p < 0.03; DRES: 15.51 ± 10.11 mm, p < 0.02) for clustering in patients with good outcome. Clustering showed shorter DSOZ and DRES compared to ECD fitting and dSPM (p < 0.05). EMSI had higher performance as outcome predictor (J = 70.63%) compared to ESI (J = 41.27%) and MSI (J = 33.33%) for clustering.

CONCLUSIONS

EMSI provides superior localization and improved predictive performance than individual modalities.

SIGNIFICANCE

EMSI can help the surgical planning and facilitate the localization of epileptogenic foci.

Multi-modal and multi-model interrogation of large-scale functional brain networks

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

Large-scale brain networks and intra-axial tumor surgery: a narrative review of functional mapping techniques, critical needs, and scientific opportunities

In recent years, a paradigm shift in neuroscience has been occurring from "localizationism," or the idea that the brain is organized into separately functioning modules, toward "connectomics," or the idea that interconnected nodes form networks as the underlying substrates of behavior and thought. Accordingly, our understanding of mechanisms of neurological function, dysfunction, and recovery has evolved to include connections, disconnections, and reconnections. Brain tumors provide a unique opportunity to probe large-scale neural networks with focal and sometimes reversible lesions, allowing neuroscientists the unique opportunity to directly test newly formed hypotheses about underlying brain structural-functional relationships and network properties. Moreover, if a more complete model of neurological dysfunction is to be defined as a "disconnectome," potential avenues for recovery might be mapped through a "reconnectome." Such insight may open the door to novel therapeutic approaches where previous attempts have failed. In this review, we briefly delve into the most clinically relevant neural networks and brain mapping techniques, and we examine how they are being applied to modern neurosurgical brain tumor practices. We then explore how brain tumors might teach us more about mechanisms of global brain dysfunction and recovery through pre- and postoperative longitudinal connectomic and behavioral analyses.

Spatiotemporal Alterations in Working Memory-Related Beta Band Neuromagnetic Activity of Patients With Schizophrenia On and Off Antipsychotic Medication: Investigation With MEG

BACKGROUND AND HYPOTHESIS

We used the uniquely high combined spatial and temporal resolution of magnetoencephalography to characterize working memory (WM)-related modulation of beta band activity in neuroleptic-free patients with schizophrenia in comparison to a large sample of performance-matched healthy controls. We also tested for effects of antipsychotic medication on identified differences in these same patients.

STUDY DESIGN

Inpatients with schizophrenia (n = 21) or psychotic disorder not otherwise specified (n = 4) completed N-back and control tasks during magnetoencephalography while on placebo and during antipsychotic medication treatment, in a blinded, randomized, counterbalanced manner. Healthy, performance-matched controls (N = 100) completed the same tasks. WM-related neural activation was estimated as beta band (14-30 Hz) desynchronization throughout the brain in successive 400 ms time windows. Voxel-wise statistical comparisons were performed between controls and patients while off-medication at each time window. Significant clusters resulting from this between-groups analysis were then used as regions-of-interest, the activations of which were compared between on- and off-medication conditions in patients.

STUDY RESULTS

Controls showed beta-band desynchronization (activation) of a fronto-parietal network immediately preceding correct button press responses-the time associated with WM updating and task execution. Altered activation in medication-free patients occurred largely during this time, in prefrontal, parietal, and visual cortices. Medication altered patients' neural responses such that the activation time courses in these regions-of-interest more closely resembled those of controls.

CONCLUSIONS

These findings demonstrate that WM-related beta band alterations in schizophrenia are time-specific and associated with neural systems targeted by antipsychotic medications. Future studies may investigate this association by examining its potential neurochemical basis.

Inter-individual single-trial classification of MEG data using M-CCA

Neuroscientific studies often involve some form of group analysis over multiple participants. This requires alignment of recordings across participants. A naive solution is to assume that participants' recordings can be aligned anatomically in sensor space. However, this assumption is likely violated due to anatomical and functional differences between individual brains. In magnetoencephalography (MEG) recordings the problem of inter-subject alignment is exacerbated by the susceptibility of MEG to individual cortical folding patterns as well as the inter-subject variability of sensor locations over the brain due to the use of a fixed helmet. Hence, an approach to combine MEG data over individual brains should relax the assumptions that a) brain anatomy and function are tightly linked and b) that the same sensors capture functionally comparable brain activation across individuals. Here we use multiset canonical correlation analysis (M-CCA) to find a common representation of MEG activations recorded from 15 participants performing a grasping task. The M-CCA algorithm was applied to transform the data of a set of multiple participants into a common space with maximum correlation between participants. Importantly, we derive a method to transform data from a new, previously unseen participant into this common representation. This makes it useful for applications that require transfer of models derived from a group of individuals to new individuals. We demonstrate the usefulness and superiority of the approach with respect to previously used approaches. Finally, we show that our approach requires only a small number of labeled data from the new participant. The proposed method demonstrates that functionally motivated common spaces have potential applications in reducing training time of online brain-computer interfaces, where models can be pre-trained on previous participants/sessions. Moreover, inter-subject alignment via M-CCA has the potential for combining data of different participants and could become helpful in future endeavors on large open datasets.

Predictive coding across the left fronto-temporal hierarchy during language comprehension

We used magnetoencephalography (MEG) and event-related potentials (ERPs) to track the time-course and localization of evoked activity produced by expected, unexpected plausible, and implausible words during incremental language comprehension. We suggest that the full pattern of results can be explained within a hierarchical predictive coding framework in which increased evoked activity reflects the activation of residual information that was not already represented at a given level of the fronto-temporal hierarchy ("error" activity). Between 300 and 500 ms, the three conditions produced progressively larger responses within left temporal cortex (lexico-semantic prediction error), whereas implausible inputs produced a selectively enhanced response within inferior frontal cortex (prediction error at the level of the event model). Between 600 and 1,000 ms, unexpected plausible words activated left inferior frontal and middle temporal cortices (feedback activity that produced top-down error), whereas highly implausible inputs activated left inferior frontal cortex, posterior fusiform (unsuppressed orthographic prediction error/reprocessing), and medial temporal cortex (possibly supporting new learning). Therefore, predictive coding may provide a unifying theory that links language comprehension to other domains of cognition.

Brain-Computer Interfaces in Visualized Medicine

The brain-computer interface (BCI), also known as a brain-machine interface (BMI), has attracted extensive attention in biomedical applications. More importantly, BCI technologies have substantially revolutionized early predictions, diagnostic techniques, and rehabilitation strategies addressing acute diseases because of BCI's innovations and clinical translations. Therefore, in this chapter, a comprehensive description of the basic concepts of BCI will be exhibited, and various visualization techniques employed in BCI's medical applications will be discussed.

Network alterations in eating epilepsy during resting state and during eating using Magnetoencephalography

OBJECTIVE

Eating epilepsy presents various imaging and electrophysiological features along with various seizure triggers. As such, network changes in eating epilepsy have not been comprehensively explored. This study was conducted to illustrate resting state network changes in eating epilepsy and to study the changes in network configurations during eating.

METHODS

Magnetoencephalography recordings of nineteen patients with drug-resistant eating epilepsy were compared with healthy controls during resting state. A subgroup of nine patients and 12 controls had MEG recordings during eating. Network changes were analyzed using phase lag index across 5 frequency bands [delta, theta, alpha, beta, and gamma] using clustering coefficient (CC), betweenness centrality (BC), path length (PL), modularity (Q), and small worldness (SW).

RESULTS

During the resting state, PL was decreased in patients with epilepsy in the delta, theta, and gamma band. Q was lower in patients with epilepsy in the beta and gamma bands. During eating, in patients with epilepsy, PL and SW were increased in all frequency bands, and Q was decreased in the beta band and increased in the rest of the frequency bands. Patients with mixed types of seizures showed higher PL in all bands except alpha, higher Q in all bands, and higher SW in the alpha and beta bands. Node-wise changes in CC and BC implicated changes in DMN and 'eating' networks.

CONCLUSION

Reflex Eating epilepsy presents with a hyperconnected network that exacerbates during eating. The cause of seizure onset and loss of consciousness in eating epilepsy might be due to aberrant network interaction between the regions of the brain involved with eating, such as the sensorimotor cortex, lateral parietal cortex, and insula with the limbic cortex and default mode network across multiple frequency bands.

Wearable OPM-MEG: A changing landscape for epilepsy

Magnetoencephalography with optically pumped magnometers (OPM-MEG) is an emerging and novel, cost-effective wearable system that can simultaneously record neuronal activity with high temporal resolution ("when" neuronal activity occurs) and spatial resolution ("where" neuronal activity occurs). This paper will first outline recent methodological advances in OPM-MEG compared to conventional superconducting quantum interference device (SQUID)-MEG before discussing how OPM-MEG can become a valuable and noninvasive clinical support tool in epilepsy surgery evaluation. Although OPM-MEG and SQUID-MEG share similar data features, OPM-MEG is a wearable design that fits children and adults, and it is also robust to head motion within a magnetically shielded room. This means that OPM-MEG can potentially extend the application of MEG into the neurobiology of severe childhood epilepsies with intellectual disabilities (e.g., epileptic encephalopathies) without sedation. It is worth noting that most OPM-MEG sensors are heated, which may become an issue with large OPM sensor arrays (OPM-MEG currently has fewer sensors than SQUID-MEG). Future implementation of triaxial sensors may alleviate the need for large OPM sensor arrays. OPM-MEG designs allowing both awake and sleep recording are essential for potential long-term epilepsy monitoring.

Analysis of the robustness and dynamics of spin-locking preparations for the detection of oscillatory magnetic fields

Extracting quantitative information of neuronal signals by non-invasive imaging is an outstanding challenge for understanding brain function and pathology. However, state-of-the-art techniques offer low sensitivity to deep electrical sources. Stimulus induced rotary saturation is a recently proposed magnetic resonance imaging sequence that detects oscillatory magnetic fields using a spin-lock preparation. Phantom experiments and simulations proved its efficiency and sensitivity, but the susceptibility of the method to field inhomogeneities is still not well understood. In this study, we simulated and analyzed the dynamic of three spin-lock preparations and their response to field inhomogeneities in the presence of a resonant oscillating field. We show that the composite spin-lock preparation is more robust against field variations within the double resonance effect. In addition, we tested the capability of the chosen composite spin-lock preparation to recover information about the spectral components of a composite signal. This study sets the bases to move one step further towards the clinical application of MR-based neuronal current imaging.

Neuroimaging Assessment of Pain

Pain is an unpleasant sensory and emotional experience. Understanding the neural mechanisms of acute and chronic pain and the brain changes affecting pain factors is important for finding pain treatment methods. The emergence and progress of non-invasive neuroimaging technology can help us better understand pain at the neural level. Recent developments in identifying brain-based biomarkers of pain through advances in advanced imaging can provide some foundations for predicting and detecting pain. For example, a neurologic pain signature (involving brain regions that receive nociceptive afferents) and a stimulus intensity-independent pain signature (involving brain regions that do not show increased activity in proportion to noxious stimulus intensity) were developed based on multivariate modeling to identify processes related to the pain experience. However, an accurate and comprehensive review of common neuroimaging techniques for evaluating pain is lacking. This paper reviews the mechanism, clinical application, reliability, strengths, and limitations of common neuroimaging techniques for assessing pain to promote our further understanding of pain.

Spatiotemporal Signatures of Surprise Captured by Magnetoencephalography

Surprise and social influence are linked through several neuropsychological mechanisms. By garnering attention, causing arousal, and motivating engagement, surprise provides a context for effective or durable social influence. Attention to a surprising event motivates the formation of an explanation or updating of models, while high arousal experiences due to surprise promote memory formation. They both encourage engagement with the surprising event through efforts aimed at understanding the situation. By affecting the behavior of the individual or a social group via setting an attractive engagement context, surprise plays an important role in shaping personal and social change. Surprise is an outcome of the brain’s function in constantly anticipating the future of sensory inputs based on past experiences. When new sensory data is different from the brain’s predictions shaped by recent trends, distinct neural signals are generated to report this surprise. As a quantitative approach to modeling the generation of brain surprise, input stimuli containing surprising elements are employed in experiments such as oddball tasks during which brain activity is recorded. Although surprise has been well characterized in many studies, an information-theoretical model to describe and predict the surprise level of an external stimulus in the recorded MEG data has not been reported to date, and setting forth such a model is the main objective of this paper. Through mining trial-by-trial MEG data in an oddball task according to theoretical definitions of surprise, the proposed surprise decoding model employs the entire epoch of the brain response to a stimulus to measure surprise and assesses which collection of temporal/spatial components in the recorded data can provide optimal power for describing the brain’s surprise. We considered three different theoretical formulations for surprise assuming the brain acts as an ideal observer that calculates transition probabilities to estimate the generative distribution of the input. We found that middle temporal components and the right and left fronto-central regions offer the strongest power for decoding surprise. Our findings provide a practical and rigorous method for measuring the brain’s surprise, which can be employed in conjunction with behavioral data to evaluate the interactive and social effects of surprising events.

Variability and task-responsiveness of electrophysiological dynamics: scale-free stability and oscillatory flexibility

Cortical oscillations and scale-free neural activity are thought to influence a variety of cognitive functions, but their differential relationships to neural stability and flexibility has never been investigated. Based on the existing literature, we hypothesize that scale-free and oscillatory processes in the brain exhibit different trade-offs between stability and flexibility; specifically, cortical oscillations may reflect variable, task-responsive aspects of brain activity, while scale-free activity is proposed to reflect a more stable and task-unresponsive aspect. We test this hypothesis using data from two large-scale MEG studies (HCP: n = 89; CamCAN: n = 195), operationalizing stability and flexibility by task-responsiveness and spontaneous intra-subject variability in resting state. We demonstrate that the power-law exponent of scale-free activity is a highly stable parameter, which responds little to external cognitive demands and shows minimal spontaneous fluctuations over time. In contrast, oscillatory power, particularly in the alpha range (8-13 Hz), responds strongly to tasks and exhibits comparatively large spontaneous fluctuations over time. In sum, our data support differential roles for oscillatory and scale-free activity in the brain with respect to neural stability and flexibility. This result carries implications for criticality-based theories of scale-free activity, state-trait models of variability, and homeostatic views of the brain with regulated variables vs. effectors.

Performance Analysis of Optically Pumped 4He Magnetometers vs. Conventional SQUIDs: From Adult to Infant Head Models

Optically pumped magnetometers (OPMs) are new, room-temperature alternatives to superconducting quantum interference devices (SQUIDs) for measuring the brain's magnetic fields. The most used OPM in MagnetoEncephaloGraphy (MEG) are based on alkali atoms operating in the spin-exchange relaxation-free (SERF) regime. These sensors do not require cooling but have to be heated. Another kind of OPM, based on the parametric resonance of ⁴He atoms are operated at room temperature, suppressing the heat dissipation issue. They also have an advantageous bandwidth and dynamic range more suitable for MEG recordings. We quantitatively assessed the improvement (relative to a SQUID magnetometers array) in recording the magnetic field with a wearable ⁴He OPM-MEG system through data simulations. The OPM array and magnetoencephalography forward models were based on anatomical MRI data from an adult, a nine-year-old child, and 10 infants aged between one month and two years. Our simulations showed that a ⁴He OPMs array offers markedly better spatial specificity than a SQUID magnetometers array in various key performance areas (e.g., signal power, information content, and spatial resolution). Our results are also discussed regarding previous simulation results obtained for alkali OPM.

Ear-EEG compares well to cap-EEG in recording auditory ERPs: a quantification of signal loss

Objective:Ear-EEG (Electroencephalography) allows to record brain activity using only a few electrodes located close to the ear. Ear-EEG is comfortable and easy to apply, facilitating beyond-the-lab EEG recordings in everyday life. With the unobtrusive setup, a person wearing it can blend in, allowing unhindered EEG recordings in social situations. However, compared to classical cap-EEG, only a small part of the head is covered with electrodes. Most scalp positions that are known from established EEG research are not covered by ear-EEG electrodes, making the comparison between the two approaches difficult and might hinder the transition from cap-based lab studies to ear-based beyond-the-lab studies.Approach:We here provide a reference data-set comparing ear-EEG and cap-EEG directly for four different auditory event-related potentials (ERPs): N100, MMN, P300 and N400. We show how the ERPs are reflected when using only electrodes around the ears.Main results:We find that significant condition differences for all ERP-components could be recorded using only ear-electrodes. The effect sizes were moderate to high on the single subject level. Morphology and temporal evolution of signals recorded from around-the-ear resemble highly those from standard scalp-EEG positions. We found a reduction in effect size (signal loss) for the ear-EEG electrodes compared to cap-EEG of 21-44%. The amount of signal loss depended on the ERP-component; we observed the lowest percentage signal loss for the N400 and the highest percentage signal loss for the N100. Our analysis further shows that no single channel position around the ear is optimal for recording all ERP-components or all participants, speaking in favor of multi-channel ear-EEG solutions.Significance:Our study provides reference results for future studies employing ear-EEG.

Age-related differences in bottom-up and top-down attention: Insights from EEG and MEG

Attention operates through top‐down and bottom‐up processes, and a balance between these processes is crucial for daily tasks. Imperilling such balance could explain ageing‐associated attentional problems such as exacerbated distractibility. In this study, we aimed to characterize this enhanced distractibility by investigating the impact of ageing upon event‐related components associated with top‐down and bottom‐up attentional processes. MEG and EEG data were acquired from 14 older and 14 younger healthy adults while performing a task that conjointly evaluates top‐down and bottom‐up attention. Event‐related components were analysed on sensor and source levels. In comparison with the younger group, the older mainly displayed (1) reduced target anticipation processes (reduced CMV), (2) increased early target processing (larger P50 but smaller N1) and (3) increased processing of early distracting sounds (larger N1 but reduced P3a), followed by a (4) prolonged reorientation towards the main task (larger RON). Taken together, our results suggest that the enhanced distractibility in ageing could stem from top‐down deficits, in particular from reduced inhibitory and reorientation processes.

People with High Empathy Show Increased Cortical Activity around the Left Medial Parieto-Occipital Sulcus after Watching Social Interaction of On-Screen Characters

People with high empathy interpret others’ mental states in daily social interactions. To investigate their characteristics of social cognitive processing, we compared neuromagnetic activities between 20 males with high empathy and 23 males with low empathy while watching social interactions between two characters. Twenty stories of four-panel comic strips were presented; the first three panels described social interactions, and the last panel described empathic/nonempathic behaviors. People with high empathy exhibited increased cortical activity in the right occipital region, medial part of the bilateral superior frontal gyri, and right posterior insula while watching social interaction scenes, which suggests that they paid attention to others’ faces and bodies, and inferred others’ mental states. They also exhibited increased cortical activity in the left superior frontal gyrus while watching empathic behaviors. Moreover, they exhibited increased cortical activity in the region around the left medial parieto-occipital sulcus, which is related to self-projection, while passively watching both empathic and nonempathic endings. Taken together, these results suggest that people with high empathy pay attention to others and actively infer others’ mental states while watching social interactions and that they reconstruct others’ mental states and intentions through self-projection after watching a sequence of others’ behaviors.

MEG correlates of temporal regularity relevant to pitch perception in human auditory cortex

We recorded neural responses in human participants to three types of pitch-evoking regular stimuli at rates below and above the lower limit of pitch using magnetoencephalography (MEG). These bandpass filtered (1–4 kHz) stimuli were harmonic complex tones (HC), click trains (CT), and regular interval noise (RIN). Trials consisted of noise-regular-noise (NRN) or regular-noise-regular (RNR) segments in which the repetition rate (or fundamental frequency F0) was either above (250 Hz) or below (20 Hz) the lower limit of pitch. Neural activation was estimated and compared at the senor and source levels.

The pitch-relevant regular stimuli (F0 = 250 Hz) were all associated with marked evoked responses at around 140 ms after noise-to-regular transitions at both sensor and source levels. In particular, greater evoked responses to pitch-relevant stimuli than pitch-irrelevant stimuli (F0 = 20 Hz) were localized along the Heschl's sulcus around 140 ms. The regularity-onset responses for RIN were much weaker than for the other types of regular stimuli (HC, CT). This effect was localized over planum temporale, planum polare, and lateral Heschl's gyrus. Importantly, the effect of pitch did not interact with the stimulus type. That is, we did not find evidence to support different responses for different types of regular stimuli from the spatiotemporal cluster of the pitch effect (∼140 ms).

The current data demonstrate cortical sensitivity to temporal regularity relevant to pitch that is consistently present across different pitch-relevant stimuli in the Heschl's sulcus between Heschl's gyrus and planum temporale, both of which have been identified as a “pitch center” based on different modalities.

Magnetoencephalography and electroencephalography can both detect differences in cortical responses to vibrotactile stimuli in individuals on the autism spectrum

Autism Spectrum (AS) is defined primarily by differences in social interactions, with impairments in sensory processing also characterizing the condition. In the search for neurophysiological biomarkers associated with traits relevant to the condition, focusing on sensory processing offers a path that is likely to be translatable across populations with different degrees of ability, as well as into animal models and across imaging modalities. In a prior study, a somatosensory neurophysiological signature of AS was identified using magnetoencephalography (MEG). Specifically, source estimation results showed differences between AS and neurotypically developing (NTD) subjects in the brain response to 25-Hz vibrotactile stimulation of the right fingertips, with lower inter-trial coherence (ITC) observed in the AS group. Here, we examined whether these group differences can be detected without source estimation using scalp electroencephalography (EEG), which is more commonly available in clinical settings than MEG, and therefore offers a greater potential for clinical translation. To that end, we recorded simultaneous whole-head MEG and EEG in 14 AS and 10 NTD subjects (age 15-28 years) using the same vibrotactile paradigm. Based on the scalp topographies, small sets of left hemisphere MEG and EEG sensors showing the maximum overall ITC were selected for group comparisons. Significant differences between the AS and NTD groups in ITC at 25 Hz as well as at 50 Hz were recorded in both MEG and EEG sensor data. For each measure, the mean ITC was lower in the AS than in the NTD group. EEG ITC values correlated with behaviorally assessed somatosensory sensation avoiding scores. The results show that information about ITC from MEG and EEG signals have substantial overlap, and thus EEG sensor-based ITC measures of the AS somatosensory processing biomarker previously identified using source localized MEG data have a potential to be developed into clinical use in AS, thanks to the higher accessibility to EEG in clinical settings.

Imaging the extent and location of spatiotemporally distributed epileptiform sources from MEG measurements

Non-invasive MEG/EEG source imaging provides valuable information about the epileptogenic brain areas which can be used to aid presurgical planning in focal epilepsy patients suffering from drug-resistant seizures. However, the source extent estimation for electrophysiological source imaging remains to be a challenge and is usually largely dependent on subjective choice. Our recently developed algorithm, fast spatiotemporal iteratively reweighted edge sparsity minimization (FAST-IRES) strategy, has been shown to objectively estimate extended sources from EEG recording, while it has not been applied to MEG recordings. In this work, through extensive numerical experiments and real data analysis in a group of focal drug-resistant epilepsy patients' interictal spikes, we demonstrated the ability of FAST-IRES algorithm to image the location and extent of underlying epilepsy sources from MEG measurements. Our results indicate the merits of FAST-IRES in imaging the location and extent of epilepsy sources for pre-surgical evaluation from MEG measurements.

Face-selective responses in combined EEG/MEG recordings with fast periodic visual stimulation (FPVS)

Fast periodic visual stimulation (FPVS) allows the recording of objective brain responses of human face categorization (i.e., generalizable face-selective responses) with high signal-to-noise ratio. This approach has been successfully employed in a number of scalp electroencephalography (EEG) studies but has not been used with magnetoencephalography (MEG) yet, let alone with combined MEG/EEG recordings and distributed source estimation. Here, we presented various natural images of faces periodically (1.2 Hz) among natural images of objects (base frequency 6 Hz) whilst recording simultaneous EEG and MEG in 15 participants. Both measurement modalities showed face-selective responses at 1.2 Hz and harmonics across participants, with high and comparable signal-to-noise ratio (SNR) in about 3 min of stimulation. The correlation of face categorization responses between EEG and two MEG sensor types was lower than between the two MEG sensor types, indicating that the two sensor modalities provide independent information about the sources of face-selective responses. Face-selective EEG responses were right-lateralized as reported previously, and were numerically but non-significantly right-lateralized in MEG data. Distributed source estimation based on combined EEG/MEG signals confirmed a more bilateral face-selective response in visual brain regions located anteriorly to the common response to all stimuli at 6 Hz and harmonics. Conventional sensor and source space analyses of evoked responses in the time domain further corroborated this result. Our results demonstrate that FPVS in combination with simultaneously recorded EEG and MEG may serve as an efficient localizer paradigm for human face categorization.

Trait sensation seeking is associated with heightened beta-band oscillatory dynamics over left ventrolateral prefrontal cortex during reward expectancy

BACKGROUND

Sensation Seeking, the proclivity toward novel and stimulating experiences, is associated with greater left ventrolateral prefrontal cortex (vlPFC) activity during uncertain reward expectancy. Here, we examined relationships between sensation seeking and vlPFC oscillatory dynamics using electroencephalography (EEG).

METHODS

In 26 adolescents/young adults (16 female; 22.3 ± 1.7yrs), EEG was measured during uncertain reward expectancy. Event-related spectral perturbations (ERSP) from 15-80 Hz (beta/gamma bands) were compared as a function of uncertain reward expected value and assessed for relationships with feedback-related negativity (FRN) response to outcome feedback and response tendency measures of risk for BD.

RESULTS

Event-related synchronization (ERS) between 15-25 Hz (beta) over left vlPFC was sensitive to the expected value of uncertain reward (rho=0.46; p = 0.048), and correlated with sensation seeking (r = 0.49, p < 0.01) and feedback-related negativity (FRN), where greater beta ERS was related to larger FRN (r = -0.39, p = 0.047). FRN was also related to behavioral inhibition (r = 0.49, p < 0.01).

LIMITATIONS

It is unknown whether results may extrapolate to clinical populations, given the healthy sample used here. Further, although we have confidence that the beta-band signal we measure in this study arises from left prefrontal cortex, we largely infer a left vlPFC source.

CONCLUSIONS

These findings highlight the role of left vlPFC in evaluation of immediate rewards. We now provide a link between reward expectancy-related left vlPFC activity and the well-characterized FRN, with a known role in attentive processing. These findings can guide treatment development for mania/hypomania at-risk individuals, including transcranial alternating current stimulation.

Dynamic relationships between spontaneous and evoked electrophysiological activity

Spontaneous neural activity fluctuations have been shown to influence trial-by-trial variation in perceptual, cognitive, and behavioral outcomes. However, the complex electrophysiological mechanisms by which these fluctuations shape stimulus-evoked neural activity remain largely to be explored. Employing a large-scale magnetoencephalographic dataset and an electroencephalographic replication dataset, we investigate the relationship between spontaneous and evoked neural activity across a range of electrophysiological variables. We observe that for high-frequency activity, high pre-stimulus amplitudes lead to greater evoked desynchronization, while for low frequencies, high pre-stimulus amplitudes induce larger degrees of event-related synchronization. We further decompose electrophysiological power into oscillatory and scale-free components, demonstrating different patterns of spontaneous-evoked correlation for each component. Finally, we find correlations between spontaneous and evoked time-domain electrophysiological signals. Overall, we demonstrate that the dynamics of multiple electrophysiological variables exhibit distinct relationships between their spontaneous and evoked activity, a result which carries implications for experimental design and analysis in non-invasive electrophysiology.

Statistical power: Implications for planning MEG studies

Statistical power is key for robust, replicable science. Here, we systematically explored how numbers of trials and subjects affect statistical power in MEG sensor-level data. More specifically, we simulated "experiments" using the MEG resting-state dataset of the Human Connectome Project (HCP). We divided the data in two conditions, injected a dipolar source at a known anatomical location in the "signal condition", but not in the "noise condition", and detected significant differences at sensor level with classical paired t-tests across subjects, using amplitude, squared amplitude, and global field power (GFP) measures. Group-level detectability of these simulated effects varied drastically with anatomical origin. We thus examined in detail which spatial properties of the sources affected detectability, looking specifically at the distance from closest sensor and orientation of the source, and at the variability of these parameters across subjects. In line with previous single-subject studies, we found that the most detectable effects originate from source locations that are closest to the sensors and oriented tangentially with respect to the head surface. In addition, cross-subject variability in orientation also affected group-level detectability, boosting detection in regions where this variability was small and hindering detection in regions where it was large. Incidentally, we observed a considerable covariation of source position, orientation, and their cross-subject variability in individual brain anatomical space, making it difficult to assess the impact of each of these variables independently of one another. We thus also performed simulations where we controlled spatial properties independently of individual anatomy. These additional simulations confirmed the strong impact of distance and orientation and further showed that orientation variability across subjects affects detectability, whereas position variability does not. Importantly, our study indicates that strict unequivocal recommendations as to the ideal number of trials and subjects for any experiment cannot be realistically provided for neurophysiological studies and should be adapted according to the brain regions under study.

A comprehensive study on electroencephalography and magnetoencephalography sensitivity to cortical and subcortical sources

Signal-to-noise ratio (SNR) maps are a good way to visualize electroencephalography (EEG) and magnetoencephalography (MEG) sensitivity. SNR maps extend the knowledge about the modulation of EEG and MEG signals by source locations and orientations and can therefore help to better understand and interpret measured signals as well as source reconstruction results thereof. Our work has two main objectives. First, we investigated the accuracy and reliability of EEG and MEG finite element method (FEM)-based sensitivity maps for three different head models, namely an isotropic three and four-compartment and an anisotropic six-compartment head model. As a result, we found that ignoring the cerebrospinal fluid leads to an overestimation of EEG SNR values. Second, we examined and compared EEG and MEG SNR mappings for both cortical and subcortical sources and their modulation by source location and orientation. Our results for cortical sources show that EEG sensitivity is higher for radial and deep sources and MEG for tangential ones, which are the majority of sources. As to the subcortical sources, we found that deep sources with sufficient tangential source orientation are recordable by the MEG. Our work, which represents the first comprehensive study where cortical and subcortical sources are considered in highly detailed FEM-based EEG and MEG SNR mappings, sheds a new light on the sensitivity of EEG and MEG and might influence the decision of brain researchers or clinicians in their choice of the best modality for their experiment or diagnostics, respectively.

Magnetoencephalography: physics, techniques, and applications in the basic and clinical neurosciences

Magnetoencephalography (MEG) is a technique used to measure the magnetic fields generated from neuronal activity in the brain. MEG has a high temporal resolution on the order of milliseconds and provides a more direct measure of brain activity when compared with hemodynamic-based neuroimaging methods such as magnetic resonance imaging and positron emission tomography. The current review focuses on basic features of MEG such as the instrumentation and the physics that are integral to the signals that can be measured, and the principles of source localization techniques, particularly the physics of beamforming and the techniques that are used to localize the signal of interest. In addition, we review several metrics that can be used to assess functional coupling in MEG and describe the advantages and disadvantages of each approach. Lastly, we discuss the current and future applications of MEG.

The 10 Common Evidence-Supported Indications for MEG in Epilepsy Surgery: An Illustrated Compendium

Unfamiliarity with the indications for and benefits of magnetoencephalography (MEG) persists, even in the epilepsy community, and hinders its acceptance to clinical practice, despite the evidence. The wide treatment gap for patients with drug-resistant epilepsy and immense underutilization of epilepsy surgery had similar effects. Thus, educating referring physicians (epileptologists, neurologists, and neurosurgeons) both about the value of epilepsy surgery and about the potential benefits of MEG can achieve synergy and greatly improve the process of selecting surgical candidates. As a practical step toward a comprehensive educational process to benefit potential MEG users, current MEG referrers, and newcomers to MEG, the authors have elected to provide an illustrated guide to 10 everyday situations where MEG can help in the evaluation of people with drug-resistant epilepsy. They are as follows: (1) lacking or imprecise hypothesis regarding a seizure onset; (2) negative MRI with a mesial temporal onset suspected; (3) multiple lesions on MRI; (4) large lesion on MRI; (5) diagnostic or therapeutic reoperation; (6) ambiguous EEG findings suggestive of "bilateral" or "generalized" pattern; (7) intrasylvian onset suspected; (8) interhemispheric onset suspected; (9) insular onset suspected; and (10) negative (i.e., spikeless) EEG. Only their practical implementation and furtherance of personal and collective education will lead to the potentially impactful synergy of the two-MEG and epilepsy surgery. Thus, while fulfilling our mission as physicians, we must not forget that ignoring the wealth of evidence about the vast underutilization of epilepsy surgery - and about the usefulness and value of MEG in selecting surgical candidates - is far from benign neglect.

Separating vascular and neuronal effects of age on fMRI BOLD signals

Accurate identification of brain function is necessary to understand the neurobiology of cognitive ageing, and thereby promote well-being across the lifespan. A common tool used to investigate neurocognitive ageing is functional magnetic resonance imaging (fMRI). However, although fMRI data are often interpreted in terms of neuronal activity, the blood oxygenation level-dependent (BOLD) signal measured by fMRI includes contributions of both vascular and neuronal factors, which change differentially with age. While some studies investigate vascular ageing factors, the results of these studies are not well known within the field of neurocognitive ageing and therefore vascular confounds in neurocognitive fMRI studies are common. Despite over 10 000 BOLD-fMRI papers on ageing, fewer than 20 have applied techniques to correct for vascular effects. However, neurovascular ageing is not only a confound in fMRI, but an important feature in its own right, to be assessed alongside measures of neuronal ageing. We review current approaches to dissociate neuronal and vascular components of BOLD-fMRI of regional activity and functional connectivity. We highlight emerging evidence that vascular mechanisms in the brain do not simply control blood flow to support the metabolic needs of neurons, but form complex neurovascular interactions that influence neuronal function in health and disease. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Influence of unfused cranial bones on magnetoencephalography signals in human infants

OBJECTIVE

To clarify the effects of unfused cranial bones on magnetoencephalography (MEG) signals during early development.

METHODS

In a simulation study, we compared the MEG signals over a spherical head model with a circular hole mimicking the anterior fontanel to those over the same head model without the fontanel for different head and fontanel sizes with varying skull thickness and conductivity.

RESULTS

The fontanel had small effects according to three indices. The sum of differences in signal over a sensor array due to a fontanel, for example, was < 6% of the sum without the fontanel. However, the fontanel effects were extensive for dipole sources deep in the brain or outside the fontanel for larger fontanels. The effects were comparable in magnitude for tangential and radial sources. Skull thickness significantly increased the effect, while skull conductivity had minor effects.

CONCLUSION

MEG signal is weakly affected by a fontanel. However, the effects can be extensive and significant for radial sources, thicker skull and large fontanels. The fontanel effects can be intuitively explained by the concept of secondary sources at the fontanel wall.

SIGNIFICANCE

The minor influence of unfused cranial bones simplifies MEG analysis, but it should be considered for quantitative analysis.

Modified covariance beamformer for solving MEG inverse problem in the environment with correlated sources

Magnetoencephalography (MEG) is a neuroimaging method ideally suited for non-invasive studies of brain dynamics. MEG's spatial resolution critically depends on the approach used to solve the ill-posed inverse problem in order to transform sensor signals into cortical activation maps. Over recent years non-globally optimized solutions based on the use of adaptive beamformers (BF) gained popularity. When operating in the environment with a small number of uncorrelated sources the BFs perform optimally and yield high spatial resolution. However, the BFs are known to fail when dealing with correlated sources acting like poorly tuned spatial filters with low signal-to-noise ratio (SNR) of the output timeseries and often meaningless cortical maps of power distribution. This fact poses a serious limitation on the broader use of this promising technique especially since fundamental mechanisms of brain functioning, its inherent symmetry and task-based experimental paradigms result into a great deal of correlation in the activity of cortical sources. To cope with this problem, we developed a novel data covariance modification approach that allows for building beamformers that maintain high spatial resolution when operating in the environments with correlated sources. At the core of our method is a projection operation applied to the vectorized sensor-space covariance matrix. This projection does not remove the activity of the correlated sources from the sensor-space covariance matrix but rather selectively handles their contributions to the covariance matrix and creates a sufficiently accurate approximation of an ideal data covariance that could hypothetically be observed should these sources be uncorrelated. Since the projection operation is reciprocal to the PSIICOS method developed by us earlier (Ossadtchi et al., 2018) we refer to the family of algorithms presented here as ReciPSIICOS. We assess the performance of the novel approach using realistically simulated MEG data and show its superior performance in comparison to the classical BF approaches and well established MNE as a method immune to source synchrony by design. We have also applied our approach to the MEG datasets from the two experiments involving two different auditory tasks. The analysis of experimental MEG datasets showed that beamformers from ReciPSIICOS family, but not the classical BF, discovered the expected bilateral focal sources in the primary auditory cortex and detected motor cortex activity associated with the audio-motor task. In most cases MNE managed well but as expected produced more spatially diffuse source distributions. Notably, ReciPSIICOS beamformers yielded cortical activity estimates with SNR several times higher than that obtained with the classical BF, which may indirectly indicate the severeness of the signal cancellation problem when applying classical beamformers to MEG signals generated by synchronous sources.

Recent advances in the noninvasive detection of high-frequency oscillations in the human brain

In recent decades, a significant body of evidence based on invasive clinical research has showed that high-frequency oscillations (HFOs) are a promising biomarker for localization of the seizure onset zone (SOZ), and therefore, have the potential to improve postsurgical outcomes in patients with epilepsy. Emerging clinical literature has demonstrated that HFOs can be recorded noninvasively using methods such as scalp electroencephalography (EEG) and magnetoencephalography (MEG). Not only are HFOs considered to be a useful biomarker of the SOZ, they also have the potential to gauge disease severity, monitor treatment, and evaluate prognostic outcomes. In this article, we review recent clinical research on noninvasively detected HFOs in the human brain, with a focus on epilepsy. Noninvasively detected scalp HFOs have been investigated in various types of epilepsy. HFOs have also been studied noninvasively in other pathologic brain disorders, such as migraine and autism. Herein, we discuss the challenges reported in noninvasive HFO studies, including the scarcity of MEG and high-density EEG equipment in clinical settings, low signal-to-noise ratio, lack of clinically approved automated detection methods, and the difficulty in differentiating between physiologic and pathologic HFOs. Additional studies on noninvasive recording methods for HFOs are needed, especially prospective multicenter studies. Further research is fundamental, and extensive work is needed before HFOs can routinely be assessed in clinical settings; however, the future appears promising.

Parametrizing the Conditionally Gaussian Prior Model for Source Localization with Reference to the P20/N20 Component of Median Nerve SEP/SEF

In this article, we focused on developing the conditionally Gaussian hierarchical Bayesian model (CG-HBM), which forms a superclass of several inversion methods for source localization of brain activity using somatosensory evoked potential (SEP) and field (SEF) measurements. The goal of this proof-of-concept study was to improve the applicability of the CG-HBM as a superclass by proposing a robust approach for the parametrization of focal source scenarios. We aimed at a parametrization that is invariant with respect to altering the noise level and the source space size. The posterior difference between the gamma and inverse gamma hyperprior was minimized by optimizing the shape parameter, while a suitable range for the scale parameter can be obtained via the prior-over-measurement signal-to-noise ratio, which we introduce as a new concept in this study. In the source localization experiments, the primary generator of the P20/N20 component was detected in the Brodmann area 3b using the CG-HBM approach and a parameter range derived from the existing knowledge of the Tikhonov-regularized minimum norm estimate, i.e., the classical Gaussian prior model. Moreover, it seems that the detection of deep thalamic activity simultaneously with the P20/N20 component with the gamma hyperprior can be enhanced while using a close-to-optimal shape parameter value.