Effects of Spatial Pattern Scale of Brain Activity on the Sensitivity of DOT,fMRI, EEG and MEG

The objective of this work is to quantify how patterns of cortical activity at different spatial scales are measured by noninvasive functional neuroimaging sensors. We simulated cortical activation patterns at nine different spatial scales in a realistic head model and propagated this activity to magnetoencephalography (MEG), electroencephalography (EEG), diffuse optical tomography (DOT), and functional magnetic resonance imaging (fMRI) sensors in arrangements that are typically used in functional neuroimaging studies. We estimated contrast transfer functions (CTF), correlation distances in sensor space, and the minimum resolvable spatial scale of cortical activity for each modality. We found that CTF decreases as the spatial extent of cortical activity decreases, and that correlations between nearby sensors depend on the spatial extent of cortical activity. For cortical activity on the intermediate spatial scale of 6.7 cm², the correlation distances (r>0.5) were 1.0 cm for fMRI, 2.0 cm for DOT, 12.8 for EEG, 9.5 cm for MEG magnetometers and 9.7 cm for MEG gradiometers. The resolvable spatial pattern scale was found to be 1.43 cm² for MEG magnetometers, 0.88 cm² for MEG gradiometers, 376 cm² for EEG, 0.75 cm² for DOT, and 0.072 cm² for fMRI. These findings show that sensitivity to cortical activity varies substantially as a function of spatial scale within and between the different imaging modalities. This information should be taken into account when interpreting neuroimaging data and when choosing the number of nodes for network analyses in sensor space.     

Modelling and detecting deep brain activity with MEG and EEG

We introduce an anatomical and electrophysiological model of deep brain structures dedicated to magnetoencephalography (MEG) and electroencephalography (EEG) source imaging. So far, most imaging inverse models considered that MEG/EEG surface signals were predominantly produced by cortical, hence superficial, neural currents. Here we question whether crucial deep brain structures such as the basal ganglia and the hippocampus may also contribute to distant, scalp MEG and EEG measurements. We first design a realistic anatomical and electrophysiological model of these structures and subsequently run Monte-Carlo experiments to evaluate the respective sensitivity of the MEG and EEG to signals from deeper origins. Results indicate that MEG/EEG may indeed localize these deeper generators, which is confirmed here from experimental MEG data reporting on the modulation of alpha (10–12 Hz) brain waves.     

Divergent Cortical Generators of MEG and EEG during Human Sleep Spindles Suggested by Distributed Source Modeling

Background

Sleep spindles are ∼1-second bursts of 10–15 Hz activity, occurring during normal stage 2 sleep. In animals, sleep spindles can be synchronous across multiple cortical and thalamic locations, suggesting a distributed stable phase-locked generating system. The high synchrony of spindles across scalp EEG sites suggests that this may also be true in humans. However, prior MEG studies suggest multiple and varying generators.

Methodology/Principal Findings

We recorded 306 channels of MEG simultaneously with 60 channels of EEG during naturally occurring spindles of stage 2 sleep in 7 healthy subjects. High-resolution structural MRI was obtained in each subject, to define the shells for a boundary element forward solution and to reconstruct the cortex providing the solution space for a noise-normalized minimum norm source estimation procedure. Integrated across the entire duration of all spindles, sources estimated from EEG and MEG are similar, diffuse and widespread, including all lobes from both hemispheres. However, the locations, phase and amplitude of sources simultaneously estimated from MEG versus EEG are highly distinct during the same spindles. Specifically, the sources estimated from EEG are highly synchronous across the cortex, whereas those from MEG rapidly shift in phase, hemisphere, and the location within the hemisphere.

Conclusions/Significance

The heterogeneity of MEG sources implies that multiple generators are active during human sleep spindles. If the source modeling is correct, then EEG spindles are generated by a different, diffusely synchronous system. Animal studies have identified two thalamo-cortical systems, core and matrix, that produce focal or diffuse activation and thus could underlie MEG and EEG spindles, respectively. Alternatively, EEG spindles could reflect overlap at the sensors of the same sources as are seen from the MEG. Although our results generally match human intracranial recordings, additional improvements are possible and simultaneous intra- and extra-cranial measures are needed to test their accuracy.     

Replicability of MEG and EEG measures of the auditory N1/N1m-response

We investigated the replicability of the source location, amplitude and latency measures of the auditory evoked N1 (EEG) and N1m (MEG) responses. Each of the 5 subjects was measured 6 times in two recording sessions. Responses to monaural stimuli were recorded from 122 MEG and 64 EEG channels simultaneously. The EEG data were modeled with a symmetrically-located dipole pair. For the MEG data, one dipole in each hemisphere was located independently using a subset of channels. Standard deviation (SD) was used as a measure for replicability. The average SD of the x, y and z coordinates of the contralateral N1m dipole was about 2 mm, whereas the corresponding figures for the ipsilateral N1m and the contra- and ipsilateral N1 were about twice as large. The SDs of the dipole amplitudes and latencies were almost equal with MEG and EEG. The amplitude and latency measures of the MEG field gradient waveforms were almost as replicable as those of the dipole models. The results suggest that both MEG and EEG can be used for investigating the simultaneous activity of the left and right auditory cortices independently, MEG being superior in certain experimental setups.     

Accounting for Linear Transformations of EEG and MEG Data in Source Analysis

Analyses of electro- and magnetoencephalography (EEG, MEG) data often involve a linear modification of signals at the sensor level. Examples include re-referencing of the EEG, computation of synthetic gradiometer in MEG, or the removal of artifactual independent components to clean EEG and MEG data. A question of practical relevance is, if such modifications must be accounted for by adapting the physical forward model (leadfield) before subsequent source analysis. Here, we show that two scenarios need to be differentiated. In the first scenario, which corresponds to re-referencing the EEG and synthetic gradiometer computation in MEG, the leadfield must be adapted before source analysis. In the second scenario, which corresponds to removing artifactual components to ‘clean’ the data, the leadfield must not be changed. We demonstrate and discuss the consequences of wrongly modifying the leadfield in the latter case for an example. Future EEG and MEG studies employing source analyses should carefully consider whether and, if so, how the leadfield must be modified as explicated here.     

A review of EEG and MEG for brainnetome research

The majority of brain activities are performed by functionally integrating separate regions of the brain. Therefore, the synchronous operation of the brain’s multiple regions or neuronal assemblies can be represented as a network with nodes that are interconnected by links. Because of the complexity of brain interactions and their varying effects at different levels of complexity, one of the corresponding authors of this paper recently proposed the brainnetome as a new –ome to explore and integrate the brain network at different scales. Because electroencephalography (EEG) and magnetoencephalography (MEG) are noninvasive and have outstanding temporal resolution and because they are the primary clinical techniques used to capture the dynamics of neuronal connections, they lend themselves to the analysis of the neural networks comprising the brainnetome. Because of EEG/MEG’s applicability to brainnetome analyses, the aim of this review is to identify the procedures that can be used to form a network using EEG/MEG data in sensor or source space and to promote EEG/MEG network analysis for either neuroscience or clinical applications. To accomplish this aim, we show the relationship of the brainnetome to brain networks at the macroscale and provide a systematic review of network construction using EEG and MEG. Some potential applications of the EEG/MEG brainnetome are to use newly developed methods to associate the properties of a brainnetome with indices of cognition or disease conditions. Associations based on EEG/MEG brainnetome analysis may improve the comprehension of the functioning of the brain in neuroscience research or the recognition of abnormal patterns in neurological disease.     

Harmony: EEG/MEG Linear Inverse Source Reconstruction in the Anatomical Basis of Spherical Harmonics

EEG/MEG source localization based on a “distributed solution” is severely underdetermined, because the number of sources is much larger than the number of measurements. In particular, this makes the solution strongly affected by sensor noise. A new way to constrain the problem is presented. By using the anatomical basis of spherical harmonics (or spherical splines) instead of single dipoles the dimensionality of the inverse solution is greatly reduced without sacrificing the quality of the data fit. The smoothness of the resulting solution reduces the surface bias and scatter of the sources (incoherency) compared to the popular minimum-norm algorithms where single-dipole basis is used (MNE, depth-weighted MNE, dSPM, sLORETA, LORETA, IBF) and allows to efficiently reduce the effect of sensor noise. This approach, termed Harmony, performed well when applied to experimental data (two exemplars of early evoked potentials) and showed better localization precision and solution coherence than the other tested algorithms when applied to realistically simulated data.     

EEG-MEG Integration Enhances the Characterization of Functional and Effective Connectivity in the Resting State Network

At the sensor level many aspects, such as spectral power, functional and effective connectivity as well as relative-power-ratio ratio (RPR) and spatial resolution have been comprehensively investigated through both electroencephalography (EEG) and magnetoencephalography (MEG). Despite this, differences between both modalities have not yet been systematically studied by direct comparison. It remains an open question as to whether the integration of EEG and MEG data would improve the information obtained from the above mentioned parameters. Here, EEG (64-channel system) and MEG (275 sensor system) were recorded simultaneously in conditions with eyes open (EO) and eyes closed (EC) in 29 healthy adults. Spectral power, functional and effective connectivity, RPR, and spatial resolution were analyzed at five different frequency bands (delta, theta, alpha, beta and gamma). Networks of functional and effective connectivity were described using a spatial filter approach called the dynamic imaging of coherent sources (DICS) followed by the renormalized partial directed coherence (RPDC). Absolute mean power at the sensor level was significantly higher in EEG than in MEG data in both EO and EC conditions. At the source level, there was a trend towards a better performance of the combined EEG+MEG analysis compared with separate EEG or MEG analyses for the source mean power, functional correlation, effective connectivity for both EO and EC. The network of coherent sources and the spatial resolution were similar for both the EEG and MEG data if they were analyzed separately. Results indicate that the combined approach has several advantages over the separate analyses of both EEG and MEG. Moreover, by a direct comparison of EEG and MEG, EEG was characterized by significantly higher values in all measured parameters in both sensor and source level. All the above conclusions are specific to the resting state task and the specific analysis used in this study to have general conclusion multi-center studies would be helpful.     

Modelling of cortical and thalamic 600 Hz activity by means of oscillatory networks

The purpose of this study is to investigate information processing in the primary somatosensory system with the help of oscillatory network modelling. Specifically, we consider interactions in the oscillatory 600 Hz activity between the thalamus and the cortical Brodmann areas 3b and 1. This type of cortical activity occurs after electrical stimulation of peripheral nerves such as the median nerve. Our measurements consist of simultaneous 31-channel MEG and 32-channel EEG recordings and individual 3D MRI data. We perform source localization by means of a multi-dipole model. The dipole activation time courses are then modelled by a set of coupled oscillators, described by linear second-order ordinary delay differential equations (DDEs). In particular, a new model for the thalamic activity is included in the oscillatory network. The parameters of the DDE system are successfully fitted to the data by a nonlinear evolutionary optimization method. To activate the oscillatory network, an individual input function is used, based on measurements of the propagated stimulation signal at the biceps. A significant feedback from the cortex to the thalamus could be detected by comparing the network modelling with and without feedback connections. Our finding in humans is supported by earlier animal studies. We conclude that this type of rhythmic brain activity can be modelled by oscillatory networks in order to disentangle feed forward and feedback information transfer.     

Measuring spectrally-resolved information transfer

Information transfer, measured by transfer entropy, is a key component of distributed computation. It is therefore important to understand the pattern of information transfer in order to unravel the distributed computational algorithms of a system. Since in many natural systems distributed computation is thought to rely on rhythmic processes a frequency resolved measure of information transfer is highly desirable. Here, we present a novel algorithm, and its efficient implementation, to identify separately frequencies sending and receiving information in a network. Our approach relies on the invertible maximum overlap discrete wavelet transform (MODWT) for the creation of surrogate data in the computation of transfer entropy and entirely avoids filtering of the original signals. The approach thereby avoids well-known problems due to phase shifts or the ineffectiveness of filtering in the information theoretic setting. We also show that measuring frequency-resolved information transfer is a partial information decomposition problem that cannot be fully resolved to date and discuss the implications of this issue. Last, we evaluate the performance of our algorithm on simulated data and apply it to human magnetoencephalography (MEG) recordings and to local field potential recordings in the ferret. In human MEG we demonstrate top-down information flow in temporal cortex from very high frequencies (above 100Hz) to both similarly high frequencies and to frequencies around 20Hz, i.e. a complex spectral configuration of cortical information transmission that has not been described before. In the ferret we show that the prefrontal cortex sends information at low frequencies (4-8 Hz) to early visual cortex (V1), while V1 receives the information at high frequencies (> 125 Hz).     

Attention modulates earliest responses in the primary auditory and visual cortices

A fundamental question about the neural correlates of attention concerns the earliest sensory processing stage that it can affect. We addressed this issue by recording magnetoencephalography (MEG) signals while subjects performed detection tasks, which required employment of spatial or nonspatial attention, in auditory or visual modality. Using distributed source analysis of MEG signals, we found that, contrary to previous studies that used equivalent current dipole (ECD) analysis, spatial attention enhanced the initial feedforward response in the primary visual cortex (V1) at 55-90 ms. We also found attentional modulation of the putative primary auditory cortex (A1) activity at 30-50 ms. Furthermore, we reproduced our findings using ECD modeling guided by the results of distributed source analysis and suggest a reason why earlier studies using ECD analysis failed to identify the modulation of earliest V1 activity.     

High-frequency (600 Hz) SEP activities originating in the subcortical and cortical human somatosensory system

Digitally high-pass filtered median nerve SEP show an oscillatory burst of low-amplitude high-frequency (600 Hz) wavelets superimposed on the N20 component which itself is generated by excitatory postsynaptic potentials of area 3b pyramidal cells. Prior studies using magnetoencephalography (MEG) localized one wavelet generator close to the primary somatosensory hand cortex. Since MEG recordings are biased towards tangentially oriented and superficial generators, a dipole source analysis of 32-channel electric SEP recordings was employed here to test for the possibility of deep and/or radially oriented burst generators: in 10 normal subjects low noise (16 000 averages) median nerve SEP were evaluated using dipole source analysis before and after applying a digital 475 Hz high-pass filter. Two main oscillatory 600 Hz burst sources were modeled; (i) a deep burst source close to the thalamus, most active in a time window between the brain-stem P14 and the cortical N20 sources, detectable in 7 of 10 subjects; most probably, this activity originates from deep axon segments of thalamocortical fibers; and (ii) a subsequent burst source timed around the N20 and located in the vicinity of the primary somatosensory hand cortex in all subjects, which was already known from MEG data. This superficial oscillatory source may be dominated by repetitive activity conducted in the terminal segments of the thalamocortical projection fibers initiated by the thalamic burst generator.     

Revealing instances of coordination among multiple cortical areas

Cognitive functions must involve interactions between several (perhaps many) cortical regions. The instances of such interactions may not be tightly time locked to any external cue. Thus averaging over repeated trials of brain activity or its spectrograms may miss these instances. Here, coordinated activity among multiple cortical locations is revealed in ongoing activity with millisecond accuracy without the need for averaging over time or frequencies. This is based on reconstructions of the cortical current dipole amplitudes at multiple points from MEG recordings. In these current dipole traces, instances of brief activity undulations (BAUs) are automatically detected and used to reveal where and when cortical points interact. The article shows that these BAUs truly represent the reorganization of activity at the cortex and are strongly connected to behavior.     

Large-Scale Functional Networks Identified from Resting-State EEG Using Spatial ICA

Several methods have been applied to EEG or MEG signals to detect functional networks. In recent works using MEG/EEG and fMRI data, temporal ICA analysis has been used to extract spatial maps of resting-state networks with or without an atlas-based parcellation of the cortex. Since the links between the fMRI signal and the electromagnetic signals are not fully established, and to avoid any bias, we examined whether EEG alone was able to derive the spatial distribution and temporal characteristics of functional networks. To do so, we propose a two-step original method: 1) An individual multi-frequency data analysis including EEG-based source localisation and spatial independent component analysis, which allowed us to characterize the resting-state networks. 2) A group-level analysis involving a hierarchical clustering procedure to identify reproducible large-scale networks across the population. Compared with large-scale resting-state networks obtained with fMRI, the proposed EEG-based analysis revealed smaller independent networks thanks to the high temporal resolution of EEG, hence hierarchical organization of networks. The comparison showed a substantial overlap between EEG and fMRI networks in motor, premotor, sensory, frontal, and parietal areas. However, there were mismatches between EEG-based and fMRI-based networks in temporal areas, presumably resulting from a poor sensitivity of fMRI in these regions or artefacts in the EEG signals. The proposed method opens the way for studying the high temporal dynamics of networks at the source level thanks to the high temporal resolution of EEG. It would then become possible to study detailed measures of the dynamics of connectivity.     

Cortical transients preceding voluntary movement

The process of initiating a voluntary muscular movement evidently involves a focusing of diffuse brain activity onto a highly specific location in the primary motor cortex. Even the very simple stereotypic movements used to study the ‘contingent negative variation’ and the ‘readiness potential’ begin with EEG indicative of widely distributed brain activity. In natural settings the involvement of diffuse cortical networks is undoubtedly even more important. Eventually, however, activity must coalesce onto specific neurons for the intended movement to ensue. Here we examine that focusing process from a mathematical point of view. Using a digital simulation, we solve the global equations for cortical dynamics and model the flow from diffuse onset to localized spike. From this perspective the interplay between global and local effects is seen as a necessary consequence of a basic cortical architecture which supports wave propagation. Watching the process evolve over time allows us to estimate some characteristic amplitudes and delays.     

Electro/magnetoencephalography beamforming with spatial restrictions based on sparsity

We present a source localization method for electroencephalographic (EEG) and magnetoencephalographic (MEG) data which is based on an estimate of the sparsity obtained through the eigencanceler (EIG), which is a spatial filter whose weights are constrained to lie in the noise subspace. The EIG provides rejection of directional interferences while minimizing noise contributions and maintaining specified beam pattern constraints. In our case, the EIG is used to estimate the sparsity of the signal as a function of the position, then we use this information to spatially restrict the neural sources to locations out of the sparsity maxima. As proof of the concept, we incorporate this restriction in the “classical” linearly constrained minimum variance (LCMV) source localization approach in order to enhance its performance. We present numerical examples to evaluate the proposed method using realistically simulated EEG/MEG data for different signal-to-noise (SNR) conditions and various levels of correlation between sources, as well as real EEG/MEG measurements of median nerve stimulation. Our results show that the proposed method has the potential of reducing the bias on the search of neural sources in the classical approach, as well as making it more effective in localizing correlated sources.     

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic–haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

What is the relationship between functional network architectures inferred from electromagnetic and haemodynamic data? This study shows that superposition of electromagnetic networks at canonical frequency bands manifests as highly structured patterns of haemodynamic functional connectivity in the human brain, reflecting systematic variation in cytoarchitecture.     

Comparative power spectral analysis of simultaneous elecroencephalographic and magnetoencephalographic recordings in humans suggests non-resistive extracellular media

The resistive or non-resistive nature of the extracellular space in the brain is still debated, and is an important issue for correctly modeling extracellular potentials. Here, we first show theoretically that if the medium is resistive, the frequency scaling should be the same for electroencephalogram (EEG) and magnetoencephalogram (MEG) signals at low frequencies (<10 Hz). To test this prediction, we analyzed the spectrum of simultaneous EEG and MEG measurements in four human subjects. The frequency scaling of EEG displays coherent variations across the brain, in general between 1/f and 1/f ², and tends to be smaller in parietal/temporal regions. In a given region, although the variability of the frequency scaling exponent was higher for MEG compared to EEG, both signals consistently scale with a different exponent. In some cases, the scaling was similar, but only when the signal-to-noise ratio of the MEG was low. Several methods of noise correction for environmental and instrumental noise were tested, and they all increased the difference between EEG and MEG scaling. In conclusion, there is a significant difference in frequency scaling between EEG and MEG, which can be explained if the extracellular medium (including other layers such as dura matter and skull) is globally non-resistive.     

Frequency Change Detection in Human Auditory Cortex

We offer a model of how human cortex detects changes in the auditory environment. Auditory change detection has recently been the object of intense investigation via the mismatch negativity (MMN). MMN is a preattentive response to sudden changes in stimulation, measured noninvasively in the electroencephalogram (EEG) and the magnetoencephalogram (MEG). It is elicited in the oddball paradigm, where infrequent deviant tones intersperse a series of repetitive standard tones. However, little apart from the participation of tonotopically organized auditory cortex is known about the neural mechanisms underlying change detection and the MMN. In the present study, we investigate how poststimulus inhibition might account for MMN and compare the effects of adaptation with those of lateral inhibition in a model describing tonotopically organized cortex. To test the predictions of our model, we performed MEG and EEG measurements on human subjects and used both small- (<1/3 octave) and large- (>5 octaves) frequency differences between the standard and deviant tones. The experimental results bear out the prediction that MMN is due to both adaptation and lateral inhibition. Finally, we suggest that MMN might serve as a probe of what stimulus features are mapped by human auditory cortex.