Linear Inverse Source Estimates Research Articles

Transition from cortical slow oscillations of sleep to spike-wave seizures

In cats, spike-wave (SW) seizures have been observed to emerge from cortical slow oscillations (CSOs) and spindles in sleep. We examined a patient's generalized SW seizures that appeared to emerge from CSOs in sleep. Dense array (256-channel) electroencephalographic (dEEG) data were recorded during long-term monitoring for localizing seizure onset for neurosurgical planning. CSOs were identified on the basis of topographic criteria from dEEG studies of normal human sleep. SW discharges were identified in the surface dEEG and examined with distributed linear inverse neural source estimation. This patient's SW discharges appeared to emerge from a series of CSOs that engaged the left frontal pole. Furthermore, 12Hz sleep spindles in this patient were often synchronized in time with the CSOs, and in some instances they shared a similar localization over the left frontal pole at the point of seizure onset. The initial discharges of each seizure engaged both the left frontal pole, which was the primary neural source of this patient's CSOs, and the left temporal lobe, which seemed critical for the evolution of this patient's seizures. The correlation of SW seizures with CSOs in animal studies may be observed in humans as well, providing clues to the pathology of arousal regulation in some cases of nocturnal epilepsy.

W17 Single and paired pulse TMS: Comparisons between natural and magnetic activation and between monophasic and biphasic TMS

ized linear inverse source estimates of cortical current density. Determination of the priors in the resolution of the linear inverse problem was performed with the use of information from the hemodynamic responses of the cortical areas as revealed by block-designed (strength of activated voxels) fMRI. We estimate functional cortical connectivity by computing the directed transfer function (DTF) and the Partial Directed Coherence (PDC) on the estimated cortical current density waveforms in regions of interest (ROIs) on the modeled cortical mantle. Furthermore, this method allows to detect changes in the time course of information flow between cortical regions in different frequency bands. Results: The reliability of these techniques was further demonstrated by simulation studies and the elaboration of high-resolution EEG and fMRI signals collected during different cognitive and motor tasks in healthy subjects. Discussion: Multimodal integration of high resolution EEG and fMRI data offers a significant aid for the study of the ‘‘when’’ , ‘‘where’’, and ‘‘how’’ in the cortical brain regions.

Filling-in in Schizophrenia: a High-density Electrical Mapping and Source-analysis Investigation of Illusory Contour Processing

Evidence is accumulating that patients with schizophrenia exhibit relatively severe deficits in early visual sensory processing within the dorsal stream, while processing within the ventral stream appears to be relatively more intact. Here, illusory contour (IC) processing was investigated in a cohort of schizophrenia patients and age-matched healthy controls using high-density visual evoked potentials (VEPs), spatiotemporal topographic analyses and the Local Auto-Regressive Average distributed linear inverse source estimation. IC processing was assessed because it is now known to be an excellent metric of early processing within regions of the ventral visual stream. Results in the present study show that IC processing (106-194 ms) is spared in patients with schizophrenia, providing strong evidence that early ventral stream processing is essentially normal. This is so despite equally strong evidence that early dorsal stream processing is severely impaired in this population, as indexed by a robust decrement in amplitude of the P1 component in patients and a large topographic difference between groups for this component (54-104 ms). Source analysis confirmed that the flow of activity into the dorsal stream was substantially decreased in patients. As such, these results suggest that some aspects of early ventral processing are not entirely reliant on intact inputs from the dorsal stream. Lastly, we show that later phases of visual processing (240-400 ms) also rely on the activity of different brain networks in controls and patients, with the latter recruiting strong frontal activity perhaps as compensation for impaired ventral stream processing during this period. We interpret the present findings in the context of a two-stage processing model. Under this model, it is suggested that the second stage of ventral stream processing is dependent on the fidelity of inputs from the dorsal visual stream and that impairment of this critical modulatory input may underlie the failure of 'higher-level' ventral stream processes in this population.

Estimation of the effective and functional human cortical connectivity with structural equation modeling and directed transfer function applied to high-resolution EEG

Different brain imaging devices are presently available to provide images of the human functional cortical activity, based on hemodynamic, metabolic or electromagnetic measurements. However, static images of brain regions activated during particular tasks do not convey the information of how these regions are interconnected. The concept of brain connectivity plays a central role in the neuroscience, and different definitions of connectivity, functional and effective, have been adopted in literature. While the functional connectivity is defined as the temporal coherence among the activities of different brain areas, the effective connectivity is defined as the simplest brain circuit that would produce the same temporal relationship as observed experimentally among cortical sites. The structural equation modeling (SEM) is the most used method to estimate effective connectivity in neuroscience, and its typical application is on data related to brain hemodynamic behavior tested by functional magnetic resonance imaging (fMRI), whereas the directed transfer function (DTF) method is a frequency-domain approach based on both a multivariate autoregressive (MVAR) modeling of time series and on the concept of Granger causality. This study presents advanced methods for the estimation of cortical connectivity by applying SEM and DTF on the cortical signals estimated from high-resolution electroencephalography (EEG) recordings, since these signals exhibit a higher spatial resolution than conventional cerebral electromagnetic measures. To estimate correctly the cortical signals, we used a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from individual MRI, a distributed source model and a regularized linear inverse source estimates of cortical current density. Before the application of SEM and DTF methodology to the cortical waveforms estimated from high-resolution EEG data, we performed a simulation study, in which different main factors (signal-to-noise ratio, SNR, and simulated cortical activity duration, LENGTH) were systematically manipulated in the generation of test signals, and the errors in the estimated connectivity were evaluated by the analysis of variance (ANOVA). The statistical analysis returned that during simulations, both SEM and DTF estimators were able to correctly estimate the imposed connectivity patterns under reasonable operative conditions, that is, when data exhibit an SNR of at least 3 and a LENGTH of at least 75 s of nonconsecutive EEG recordings at 64 Hz of sampling rate. Hence, effective and functional connectivity patterns of cortical activity can be effectively estimated under general conditions met in any practical EEG recordings, by combining high-resolution EEG techniques and linear inverse estimation with SEM or DTF methods. We conclude that the estimation of cortical connectivity can be performed not only with hemodynamic measurements, but also with EEG signals treated with advanced computational techniques.

Multimodal integration of EEG, MEG and fMRI data for the solution of the neuroimage puzzle

In this paper, advanced methods for the modeling of human cortical activity from combined high-resolution electroencephalography (EEG), magnetoencephalography (MEG) and functional magnetic resonance imaging (fMRI) data are presented. These methods include a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from magnetic resonance images, multidipole source model and regularized linear inverse source estimates of cortical current density. Determination of the priors in the resolution of the linear inverse problem was performed with the use of information from the hemodynamic responses of the cortical areas as revealed by block-designed (strength of activated voxels) fMRI. Examples of the application of these methods to the estimation of the time varying cortical current density activity in selected region of interest (ROI) are presented for movement-related high-resolution EEG data.

Estimation of the cortical functional connectivity with the multimodal integration of high-resolution EEG and fMRI data by directed transfer function

Nowadays, several types of brain imaging device are available to provide images of the functional activity of the cerebral cortex based on hemodynamic, metabolic, or electromagnetic measurements. However, static images of brain regions activated during particular tasks do not convey the information of how these regions communicate with each other. In this study, advanced methods for the estimation of cortical connectivity from combined high-resolution electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) data are presented. These methods include a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from individual magnetic resonance images, multidipole source model, and regularized linear inverse source estimates of cortical current density. Determination of the priors in the resolution of the linear inverse problem was performed with the use of information from the hemodynamic responses of the cortical areas as revealed by block-designed (strength of activated voxels) fMRI. We estimate functional cortical connectivity by computing the directed transfer function (DTF) on the estimated cortical current density waveforms in regions of interest (ROIs) on the modeled cortical mantle. The proposed method was able to unveil the direction of the information flow between the cortical regions of interest, as it is directional in nature. Furthermore, this method allows to detect changes in the time course of information flow between cortical regions in different frequency bands. The reliability of these techniques was further demonstrated by elaboration of high-resolution EEG and fMRI signals collected during visually triggered finger movements in four healthy subjects. Connectivity patterns estimated for this task reveal an involvement of right parietal and bilateral premotor and prefrontal cortical areas. This cortical region involvement resembles that revealed in previous studies where visually triggered finger movements were analyzed with the use of separate EEG or fMRI measurements.

Setting boundaries: brain dynamics of modal and amodal illusory shape completion in humans.

Normal visual perception requires differentiating foreground from background objects. Differences in physical attributes sometimes determine this relationship. Often such differences must instead be inferred, as when two objects or their parts have the same luminance. Modal completion refers to such perceptual "filling-in" of object borders that are accompanied by concurrent brightness enhancement, in turn termed illusory contours (ICs). Amodal completion is filling-in without concurrent brightness enhancement. Presently there are controversies regarding whether both completion processes use a common neural mechanism and whether perceptual filling-in is a bottom-up, feedforward process initiating at the lowest levels of the cortical visual pathway or commences at higher-tier regions. We previously examined modal completion (Murray et al., 2002) and provided evidence that the earliest modal IC sensitivity occurs within higher-tier object recognition areas of the lateral occipital complex (LOC). We further proposed that previous observations of IC sensitivity in lower-tier regions likely reflect feedback modulation from the LOC. The present study tested these proposals, examining the commonality between modal and amodal completion mechanisms with high-density electrical mapping, spatiotemporal topographic analyses, and the local autoregressive average distributed linear inverse source estimation. A common initial mechanism for both types of completion processes (140 msec) that manifested as a modulation in response strength within higher-tier visual areas, including the LOC and parietal structures, is demonstrated, whereas differential mechanisms were evident only at a subsequent time period (240 msec), with amodal completion relying on continued strong responses in these structures.

Solving the “neuroimaging puzzle” with the multimodal integration of EEG and functional magnetic resonance recordings

In this paper, advanced methods for the modeling of human cortical activity from combined high-resolution electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) data are reviewed. These methods include a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from magnetic resonance images, multidipole source model, and regularized linear inverse source estimates. Determination of the priors in the resolution of the linear inverse problem was performed with the use of information from the hemodynamic responses of the cortical areas as revealed by block-designed fMRI.

FMRI Priors for the Linear Inverse Estimation of EEG Cortical Sources

In this paper, advanced methods for the modeling of human cortical activity from combined high-resolution electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) data are reviewed. These methods include a subject's multicompartment head model (scalp, skull, dura mater, cortex) constructed from magnetic resonance images, multidipole source model, and regularized linear inverse source estimate based on boundary element mathematics. Furthermore, determination of the priors in the resolution of the linear inverse problem was performed with the use of information from the hemodynamic responses of the cortical areas as revealed by block-designed (strength of activated voxels) and event-related (coupling of activated voxels) fMRI. Linear inverse source estimates of cortical activity were regularized by taking into account the covariance of background EEG sensor noise. As an example, these methods were applied to EEG (128 electrodes) and fMRI data, which were recorded in separate sessions while normal subjects executed voluntary right one-digit movements.

Linear inverse source estimate of combined EEG and MEG data related to voluntary movements.

A method for the modeling of human movement-related cortical activity from combined electroencephalography (EEG) and magnetoencephalography (MEG) data is proposed. This method includes a subject's multi-compartment head model (scalp, skull, dura mater, cortex) constructed from magnetic resonance images, multi-dipole source model, and a regularized linear inverse source estimate based on boundary element mathematics. Linear inverse source estimates of cortical activity were regularized by taking into account the covariance of background EG and MEG sensor noise. EEG (121 sensors) and MEG (43 sensors) data were recorded in separate sessions whereas normal subjects executed voluntary right one-digit movements. Linear inverse source solution of EEG, MEG, and EEG-MEG data were quantitatively evaluated by using three performance indexes. The first two indexes (Dipole Localization Error [DLE] and Spatial Dispersion [SDis]) were used to compute the localization power for the source solutions obtained. Such indexes were based on the information provided by the column of the resolution matrix (i.e., impulse response). Ideal DLE values tend to zero (the source current was correctly retrieved by the procedure). In contrast, high DLE values suggest severe mislocalization in the source reconstruction. A high value of SDis at a source space point mean that such a source will be retrieved by a large area with the linear inverse source estimation. The remaining performance index assessed the quality of the source solution based on the information provided by the rows of the resolution matrix R, i.e., resolution kernels. The i-th resolution kernels of the matrix R describe how the estimation of the i-th source is distorted by the concomitant activity of all other sources. A statistically significant lower dipole localization error was observed and lower spatial dispersion in source solutions produced by combined EEG-MEG data than from EEG and MEG data considered separately (P < 0.05). These effects were not due to an increased number of sensors in the combined EEG-MEG solutions. They result from the independence of source information conveyed by the multimodal measurements. From a physiological point of view, the linear inverse source solution of EEG-MEG data suggested a contralaterally preponderant bilateral activation of primary sensorimotor cortex from the preparation to the execution of the movement. This activation was associated with that of the supplementary motor area. The activation of bilateral primary sensorimotor cortical areas was greater during the processing of afferent information related to the ongoing movement than in the preparation for the motor act. In conclusion, the linear inverse source estimate of combined MEG and EEG data improves the estimate of movement-related cortical activity.

High-resolution electro-encephalogram: source estimates of Laplacian-transformed somatosensory-evoked potentials using a realistic subject head model constructed from magnetic resonance images.

A novel high-resolution electro-encephalographic (EEG) procedure is proposed, including high spatial sampling (128 channels), a realistic magnetic resonance-constructed subject head model, a multi-dipole cortical source model and regularised weighted minimum-norm linear inverse source estimation (WMN). As an innovation, EEG potentials (two healthy subjects; median-nerve, short-latency somatosensory-evoked potentials (SEPs)) are preliminarily Laplacian-transformed (LT) to remove brain electrical activity generated by subcortical sources (i.e. not represented in the source model). LT-WMN estimates are mathematically evaluated by figures of merit (WMN estimates as a reference). Results show higher dipole identifiability (0.69; 0.088), lower dipole localisation error (0.6 mm; 7.8 mm) and lower spatial dispersion (8.6 mm; 24 mm) in LT-WMN than in WMN estimates (Bonferroni corrected p < 0.001). These estimates are presented on the subject modelled cortical surface to highlight the increased spatial information content in LT-WMN compared with WMN estimates. The proposed high-resolution EEG technique is useful for the study of somatosensory functions in basic research and clinical applications.

Integration of High Resolution EEG and Functional Magnetic Resonance in the Study of Human Movement-Related Potentials

Cortical sources of human movement-related potentials (i.e. unilateral finger extension) were modeled using functional magnetic resonance imaging (fMR) data as a constraint of a linear inverse source estimation from highly sampled (128 channels) EEG data. Remarkably, this estimation was performed within realistic subject's MR-constructed head models by boundary element techniques. An appropriate figure of merit served to set the optimal amount of fMR constraints. With respect to standard linear inverse source estimates, fMR-constrained ones presented increased spatial detail and provided a more reliable timing of activation in bilateral sensorimotor cortical regions of interest.