Dipole separability in a neuromagnetic source analysis

Explaining event-related fields by a mechanistic model encapsulating the anatomical structure of auditory cortex

Event-related fields of the magnetoencephalogram are triggered by sensory stimuli and appear as a series of waves extending hundreds of milliseconds after stimulus onset. They reflect the processing of the stimulus in cortex and have a highly subject-specific morphology. However, we still have an incomplete picture of how event-related fields are generated, what the various waves signify, and why they are so subject-specific. Here, we focus on this problem through the lens of a computational model which describes auditory cortex in terms of interconnected cortical columns as part of hierarchically placed fields of the core, belt, and parabelt areas. We develop an analytical approach arriving at solutions to the system dynamics in terms of normal modes: damped harmonic oscillators emerging out of the coupled excitation and inhibition in the system. Each normal mode is a global feature which depends on the anatomical structure of the entire auditory cortex. Further, normal modes are fundamental dynamical building blocks, in that the activity of each cortical column represents a combination of all normal modes. This approach allows us to replicate a typical auditory event-related response as a weighted sum of the single-column activities. Our work offers an alternative to the view that the event-related field arises out of spatially discrete, local generators. Rather, there is only a single generator process distributed over the entire network of the auditory cortex. We present predictions for testing to what degree subject-specificity is due to cross-subject variations in dynamical parameters rather than in the cortical surface morphology.

From evoked potentials to cortical currents: Resolving V1 and V2 components using retinotopy constrained source estimation without fMRI

Despite evoked potentials' (EP) ubiquity in research and clinical medicine, insights are limited to gross brain dynamics as it remains challenging to map surface potentials to their sources in specific cortical regions. Multiple sources cancellation due to cortical folding and cross-talk obscures close sources, e.g. between visual areas V1 and V2. Recently retinotopic functional magnetic resonance imaging (fMRI) responses were used to constrain source locations to assist separating close sources and to determine cortical current generators. However, an fMRI is largely infeasible for routine EP investigation. We developed a novel method that replaces the fMRI derived retinotopic layout (RL) by an approach where the retinotopy and current estimates are generated from EEG or MEG signals and a standard clinical T1-weighted anatomical MRI. Using the EEG-RL, sources were localized to within 2 mm of the fMRI-RL constrained localized sources. The EEG-RL also produced V1 and V2 current waveforms that closely matched the fMRI-RL's (n = 2) r(1,198)  = 0.99, P < 0.0001. Applying the method to subjects without fMRI (n = 4) demonstrates it generates waveforms that agree closely with the literature. Our advance allows investigators with their current EEG or MEG systems to create a library of brain models tuned to individual subjects' cortical folding in retinotopic maps, and should be applicable to auditory and somatosensory maps. The novel method developed expands EP's ability to study specific brain areas, revitalizing this well-worn technique. Hum Brain Mapp 37:1696-1709, 2016. © 2016 Wiley Periodicals, Inc.

Source cancellation profiles of electroencephalography and magnetoencephalography

Recorded electric potentials and magnetic fields due to cortical electrical activity have spatial spread even if their underlying brain sources are focal. Consequently, as a result of source cancellation, loss in signal amplitude and reduction in the effective signal-to-noise ratio can be expected when distributed sources are active simultaneously. Here we investigate the cancellation effects of EEG and MEG through the use of an anatomically correct forward model based on structural MRI acquired from 7 healthy adults. A boundary element model (BEM) with four compartments (brain, cerebrospinal fluid, skull and scalp) and highly accurate cortical meshes (~300,000 vertices) were generated. Distributed source activations were simulated using contiguous patches of active dipoles. To investigate cancellation effects in both EEG and MEG, quantitative indices were defined (source enhancement, cortical orientation disparity) and computed for varying values of the patch radius as well as for automatically parcellated gyri and sulci. Results were calculated for each cortical location, averaged over all subjects using a probabilistic atlas, and quantitatively compared between MEG and EEG. As expected, MEG sensors were found to be maximally sensitive to signals due to sources tangential to the scalp, and minimally sensitive to radial sources. Compared to EEG, however, MEG was found to be much more sensitive to signals generated antero-medially, notably in the anterior cingulate gyrus. Given that sources of activation cancel each other according to the orientation disparity of the cortex, this study provides useful methods and results for quantifying the effect of source orientation disparity upon source cancellation.

Face activated neurodynamic cortical networks

Previous neuroimaging studies have shown that complex visual stimuli, such as faces, activate multiple brain regions, yet little is known on the dynamics and complexity of the activated cortical networks during the entire measurable evoked response. In this study, we used simulated and face-evoked empirical MEG data from an oddball study to investigate the feasibility of accurate, efficient, and reliable spatio-temporal tracking of cortical pathways over prolonged time intervals. We applied a data-driven, semiautomated approach to spatio-temporal source localization with no prior assumptions on active cortical regions to explore non-invasively face-processing dynamics and their modulation by task. Simulations demonstrated that the use of multi-start downhill simplex and data-driven selections of time intervals submitted to the Calibrated Start Spatio-Temporal (CSST) algorithm resulted in improved accuracy of the source localization and the estimation of the onset of their activity. Locations and dynamics of the identified sources indicated a distributed cortical network involved in face processing whose complexity was task dependent. This MEG study provided the first non-invasive demonstration, agreeing with intracranial recordings, of an early onset of the activity in the fusiform face gyrus (FFG), and that frontal activation preceded parietal for responses elicited by target faces.

The folding fingerprint of visual cortex reveals the timing of human V1 and V2

Primate neocortex contains over 30 visual areas. Recent techniques such as functional magnetic resonance imaging (fMRI) have successfully identified many of these areas in the human brain, but have been of limited value for revealing the temporal dynamics between visual areas. The electroencephalogram (EEG) provides information with high temporal precision, but has had limited success separating out the signals from individual neighboring cortical areas. Consequently, controversies exist over the temporal dynamics across cortical areas. In order to address this problem we developed a new method to identify the sources of the EEG. An individual's unique cortical pattern of sulci and gyri along with a visual area's functional retinotopic layout provides a folding fingerprint that predicts specific scalp topographies for stimuli presented in different parts of the visual field. Using this folding fingerprint with a 96 or 192 location stimulus severely constrains the solution space making it relatively easy to extract the temporal response of multiple visual areas to multiple stimulus locations. The large number of stimuli also provides a means to validate the waveforms by comparing across stimulus sets, an important feature not present in most EEG source identification procedures. Using this method our data reveal that both V1 and V2 waveforms have similar onset latencies, and their temporal dynamics provide new information regarding the response latencies of these areas in humans. Our method enables the previously unattainable separation of EEG responses from neighboring brain areas. While we applied the method to the first two cortical visual areas, V1 and V2, this method is also applicable to somatosensory areas that have defined mappings. This method provides a means to study the rapid information flow in the human brain to reveal top-down and bottom-up cognitive processes.

Effects of the task of categorizing FM direction on auditory evoked magnetic fields in the human auditory cortex

We examined effects of the task of categorizing linear frequency-modulated (FM) sweeps into rising and falling on auditory evoked magnetic fields (AEFs) from the human auditory cortex, recorded by means of whole-head magnetoencephalography. AEFs in this task condition were compared with those in a passive condition where subjects had been asked to just passively listen to the same stimulus material. We found that the M100-peak latency was significantly shorter for the task condition than for the passive condition in the left but not in the right hemisphere. Furthermore, the M100-peak latency was significantly shorter in the right than in the left hemisphere for the passive and the task conditions. In contrast, the M100-peak amplitude did not differ significantly between conditions, nor between hemispheres. We also analyzed the activation strength derived from the integral of the absolute magnetic field over constant time windows between stimulus onset and 260 ms. We isolated an early, narrow time range between about 60 ms and 80 ms that showed larger values in the task condition, most prominently in the right hemisphere. These results add to other imaging and lesion studies which suggest a specific role of the right auditory cortex in identifying FM sweep direction and thus in categorizing FM sweeps into rising and falling.

Effects of language experience: Neural commitment to language-specific auditory patterns

Linguistic experience alters an individual's perception of speech. We here provide evidence of the effects of language experience at the neural level from two magnetoencephalography (MEG) studies that compare adult American and Japanese listeners' phonetic processing. The experimental stimuli were American English /ra/ and /la/ syllables, phonemic in English but not in Japanese. In Experiment 1, the control stimuli were /ba/ and /wa/ syllables, phonemic in both languages; in Experiment 2, they were non-speech replicas of /ra/ and /la/. The behavioral and neuromagnetic results showed that Japanese listeners were less sensitive to the phonemic /r-l/ difference than American listeners. Furthermore, processing non-native speech sounds recruited significantly greater brain resources in both hemispheres and required a significantly longer period of brain activation in two regions, the superior temporal area and the inferior parietal area. The control stimuli showed no significant differences except that the duration effect in the superior temporal cortex also applied to the non-speech replicas. We argue that early exposure to a particular language produces a "neural commitment" to the acoustic properties of that language and that this neural commitment interferes with foreign language processing, making it less efficient.

Neuromagnetic localization of rhythmic activity in the human brain: a comparison of three methods

Cortical rhythmic activity is increasingly employed for characterizing human brain function. Using MEG, it is possible to localize the generators of these rhythms. Traditionally, the source locations have been estimated using sequential dipole modeling. Recently, two new methods for localizing rhythmic activity have been developed, Dynamic Imaging of Coherent Sources (DICS) and Frequency-Domain Minimum Current Estimation (MCE(FD)). With new analysis methods emerging, the researcher faces the problem of choosing an appropriate strategy. The aim of this study was to compare the performance and reliability of these three methods. The evaluation was performed using measured data from four healthy subjects, as well as with simulations of rhythmic activity. We found that the methods gave comparable results, and that all three approaches localized the principal sources of oscillatory activity very well. Dipole modeling is a very powerful tool once appropriate subsets of sensors have been selected. MCE(FD) provides simultaneous localization of sources and was found to give a good overview of the data. With DICS, it was possible to separate close-by sources that were not retrieved by the other two methods.

Model Selection in Spatio-Temporal Electromagnetic Source Analysis

Several methods [model selection procedures (MSPs)] to determine the number of sources in electroencephalogram (EEG) and magnetoencphalogram (MEG) data have previously been investigated in an instantaneous analysis. In this paper, these MSPs are extended to a spatio-temporal analysis if possible. It is seen that the residual variance (RV) tends to overestimate the number of sources. The Akaike information criterion (AIC) and the Wald test on amplitudes (WA) and the Wald test on locations (WL) have the highest probabilities of selecting the correct number of sources. The WA has the advantage that it offers the opportunity to test which source is active at which time sample.

Magnetoencephalography and its Achilles' heel

Magnetoencephalography (MEG) has practically unlimited temporal resolution. Fundamental physical reasons, however, restrict the capability of MEG to separate simultaneously active sources. After a brief tutorial introduction into MEG, various aspects of spatial resolution are reviewed with the help of examples. First the estimation of a single current dipole is examined. A consideration of the resolution field shows that the spatial selectivity of the estimated dipole moment is highly dependent on methodological issues. A subsequent consideration of various two-dipole configurations illustrates how the topography of the magnetic field depends on the distance between the two dipoles and their relative orientations. The resolution fields associated with the estimation of the dipole moments reveal a strong interference for closely spaced dipoles. A simple model suggests that the standard deviations of the estimated moments are inversely proportional to the distance of the dipoles. Spatial information provided by techniques like functional magnetic resonance imaging (fMRI) could help to overcome problems resulting from the limited spatial resolution of MEG (multimodal integration). But a straightforward synthesis, according to the principle that fMRI provides the spatial structure of the sources and MEG adds the temporal information, is probably doomed to failure in many situations. A serious dilemma, among other problems, is that the fMRI signal generally represents a temporal integral over several seconds: The knowledge that a certain brain region was active sometime or other is not necessarily helpful for disentangling the MEG activity within a specified short time window. An intriguing fact is that the spatio-temporal pattern of the MEG signals can be considered as a signature of the brain which is suitable for hypothesis testing with high temporal and spatial resolution.

Model selection in electromagnetic source analysis with an application to VEFs

In electromagnetic source analysis, it is necessary to determine how many sources are required to describe the electroencephalogram or magnetoencephalogram adequately. Model selection procedures (MSPs) or goodness of fit procedures give an estimate of the required number of sources. Existing and new MSPs are evaluated in different source and noise settings: two sources which are close or distant and noise which is uncorrelated or correlated. The commonly used MSP residual variance is seen to be ineffective, that is it often selects too many sources. Alternatives like the adjusted Hotelling's test, Bayes information criterion and the Wald test on source amplitudes are seen to be effective. The adjusted Hotelling's test is recommended if a conservative approach is taken and MSPs such as Bayes information criterion or the Wald test on source amplitudes are recommended if a more liberal approach is desirable. The MSPs are applied to empirical data (visual evoked fields).

Comparison of minimum current estimate and dipole modeling in the analysis of simulated activity in the human visual cortices

Magnetoencephalographic(MEG) data are typically interpreted using source models because of the nonunique inverse problem. Although single current dipoles, adequately representing local active areas, can be identified accurately, multiple and overlapping sources form a challenge for MEG modeling. We tested the performances of multidipole modeling and minimum current estimate (MCE) in the analysis of complicated source configurations. Simulated current sources were placed to physiologically meaningful areas of the human visual cortices. Ten volunteers from the laboratory staff analyzed four different simulations with both dipole modeling and MCE without prior information of the sources. In general, the same sources were found using both modeling methods. The subjects tended to report more false sources with MCE than with dipole model, in part due to their inexperience with the method. Dipole model was more accurate than MCE both in time and space for nonsimultaneous sources but both methods performed similarly when sources overlapped in time. For all source configurations, considerably smaller source amplitudes were reported with MCE than with dipole model.

A quantitative assessment of the sensitivity of whole-head MEG to activity in the adult human cortex

MagnetoEncephaloGraphy (MEG) relies on the detection of cortical current flow by measurement of the associated magnetic field outside the head. The amplitude of this magnetic field depends strongly on the depth of the electrical brain activity. Additionally, radially orientated sources are magnetically silent in a concentrically homogeneous volume conductor, giving rise to the anecdotal assumptions that MEG is insensitive to both deep and gyral sources. Utilising cortical surfaces extracted from Magnetic Resonance Images (MRIs) of two adult brains we constructed all possible single source elements and examined the proportion of active neocortex that is actually detectable with a whole-head MEG system. We identified those electrically active regions to which MEG is maximally sensitive by analytically computing the probability of detecting a source within a specified confidence volume. Our findings show that source depth, and not orientation, is the main factor that compromises the sensitivity of MEG to activity in the adult human cortex. There are thin strips (approximately 2 mm wide) of poor resolvability at the crests of gyri; however, these strips account for only a relatively small proportion of the cortical area and are abutted by elements with nominal tangential component yet high resolvability due to their proximity to the sensor array. Finally, we varied the extent of the patches of cortical activity, showing that small patches have a small net-current moment and are therefore less visible whereas large patches have a strong net-current moment, are generally more visible to the MEG system, yet are less appropriately modelled as single dipoles.

Latency of auditory evoked field deflection N100m ruled by pitch or spectrum?

The auditory evoked field (AEF) in response to pure tones of 250 and 1000 Hz and a complex tone with a periodicity of 4 ms (composed of the frequencies 1000, 1250, 1500, 1750, and 2000 Hz), corresponding to a pitch of 250 Hz, was recorded with a 37-channel neuromagnetometer system. The intensity was 60 dB sensation level (SL). Two different stimulus durations were examined in 12 subjects: 500 ms (long tones) and 100 ms (short tones). The stimulus onset asynchrony (SOA) was uniformly distributed between 3 and 4 s for the long tones and between 0.8 and 1.2 s for the short tones. Each subject was investigated four times, to assess the intraindividual variability. The mean latency of the AEF deflection N100m turned out to be similar for the long and the short tones: about 98 and 87 ms for the pure tones of 250 Hz and 1000 Hz, respectively, and 95 ms for the complex tone with a pitch of 250 Hz. However, a great interindividual variability was observed, exhibiting no consistent relationship between the N100m latencies for the three different tones, except that the response to the pure tone of 1000 Hz generally occurred earlier. In conclusion, this study does not support the proposal that the N100m latency represents a code for pitch, although a low pitch appears to be a factor favoring a longer N100m latency.

Electric potential produced by a dipole in a homogeneous conducting sphere

The potential produced by a dipole in a homogeneous conducting sphere is useful in simulation study, and the current available solutions still suffer from some shortcomings. In this communication, a closed solution is developed for the precise calculation of the potential anywhere in the spherical model.

Supplementary motor area activation preceding voluntary movement is detectable with a whole-scalp magnetoencephalography system

Despite the fact that the knowledge about the structure and the function of the supplementary motor area (SMA) is steadily increasing, the role of the SMA in the human brain, e.g., the contribution of the SMA to the Bereitschaftspotential, still remains unclear and controversial. The goal of this study was to contribute further to this discussion by taking advantage of the increased spatial information of a whole-scalp magnetoencephalography (MEG) system enabling us to record the magnetic equivalent of the Bereitschaftspotential 1, the Bereitschaftsfeld 1 (BF 1) or readiness field 1. Five subjects performed a complex, and one subject a simple, finger-tapping task. It was possible to record the BF 1 for all subjects. The first appearance of the BF 1 was in the range of -1.9 to -1.7 s prior to movement onset, except for the subject performing the simple task (-1 s). Analysis of the development of the magnetic field distribution and the channel waveforms showed the beginning of the Bereitschaftsfeld 2 (BF 2) or readiness field 2 at about -0.5 s prior to movement onset. In the time range of BF 1, dipole source analysis localized the source in the SMA only, whereas dipole source analysis containing also the time range of BF 2 resulted in dipole models, including dipoles in the primary motor area. In summary, with a whole-head MEG system, it was possible for the first time to detect SMA activity in healthy subjects with MEG.

Visualization of magnetoencephalographic data using minimum current estimates

The locations of active brain areas can be estimated from the magnetic field the neural current sources produce. In this work we study a visualization method of magnetoencephalographic data that is based on minimum[symbol: see text] (1)-norm estimates. The method can represent several local or distributed sources and does not need explicit a priori information. We evaluated the performance of the method using simulation studies. In a situation resembling typical magnetoencephalographic measurement, the mean estimated source strength exceeded baseline level up to 2 cm from the simulated point-like source. The method can also visualize several sources, activated simultaneously or in a sequence, which we demonstrated by analyzing magnetic responses associated with sensory stimulation and a picture naming task.