Fast realistic modeling in bioelectromagnetism using lead-field interpolation.

The practical use of realistic models in bioelectromagnetism is limited by the time-consuming amount of numerical calculations. We propose a method leading to much higher speed than currently available, and compatible with any kind of numerical methods (boundary elements (BEM), finite elements, finite differences). Illustrated with the BEM for EEG and MEG, it applies to ECG and MCG as well. The principle is two-fold. First, a Lead-Field matrix is calculated (once for all) for a grid of dipoles covering the brain volume. Second, any forward solution is interpolated from the pre-calculated Lead-Fields corresponding to grid dipoles near the source. Extrapolation is used for shallow sources falling outside the grid. Three interpolation techniques were tested: trilinear, second-order Bézier (Bernstein polynomials), and 3D spline. The trilinear interpolation yielded the highest speed gain, with factors better than x10,000 for a 9,000-triangle BEM model. More accurate results could be obtained with the Bézier interpolation (speed gain approximately 1,000), which, combined with a 8-mm step grid, lead to intrinsic localization and orientation errors of only 0.2 mm and 0.2 degrees. Further improvements in MEG could be obtained by interpolating only the contribution of secondary currents. Cropping grids by removing shallow points lead to a much better estimation of the dipole orientation in EEG than when solving the forward problem classically, providing an efficient alternative to locally refined models. This method would show special usefulness when combining realistic models with stochastic inverse procedures (simulated annealing, genetic algorithms) requiring many forward calculations.

Early signs of visual categorization for biological and non-biological stimuli in humans.

In a previous experiment aimed at studying gender processing from faces, we had found unexpected early ERP differences (45-85 ms) in task-irrelevant stimuli between a condition in which the stimuli of each gender were delivered in separate runs, and a condition in which the stimuli of both genders were mixed. Similar effects were observed with hand stimuli. These early ERP differences were tentatively related to incidental categorization processes between male and female stimuli. The present study was designed to test the robustness of these early effects for faces, and to examine whether similar effects can also be generated between two classes of non-biological stimuli. We replicated the previous findings for faces, and found similar early differential effects (50-65 ms) for non-biological stimuli (grey and hatched geometrical shapes) only, however, when the two shape categories were separated by conspicuous visual characteristics. While these results can partly be explained by phenomena related to neuronal habituation in the visual cortex, they may also suggest the existence of coarse and automatic categorization processes for rapid distinction between two wide classes of stimuli with strong psychosocial significance for humans.

Neurophysiological mechanisms of auditory selective attention in humans.

This chapter reviews the main data on the physiological substrates of auditory selective attention and their contribution to theoretical models of cognitive psychology. While event-related potentials, magnetoencephalography, and more recently neuroimaging techniques have provided fundamental information on the neural correlates of attention in the central cortical system, measurements of the frequency-following responses in the brainstem and evoked otoacoustic emissions at the cochlea strongly suggest attentional phenomena at the auditory periphery. We propose an adaptive filtering mechanism for selective auditory attention that can be flexibly and dynamically tuned depending on the attentional demand.

Neurophysiological correlates of face gender processing in humans.

Event-related potentials (ERPs) were recorded while subjects were involved in three gender-processing tasks based on human faces and on human hands. In one condition all stimuli were only of one gender, preventing any gender discrimination. In a second condition, faces (or hands) of men and women were intermixed but the gender was irrelevant for the subject's task; hence gender discrimination was assumed to be incidental. In the third condition, the task required explicit gender discrimination; gender processing was therefore assumed to be intentional. Gender processing had no effect on the occipito-temporal negative potential at approximately 170 ms after stimulation (N170 component of the ERP), suggesting that the neural mechanisms involved in the structural encoding of faces are different from those involved in the extraction of gender-related facial features. In contrast, incidental and intentional processing of face (but not hand) gender affected the ERPs between 145 and 185 ms from stimulus onset at more anterior scalp locations. This effect was interpreted as evidence for the direct visual processing of faces as described in Bruce and Young's model [Bruce, V. & Young, A. (1986) Br. J. Psychol., 77, 305-327]. Additional gender discrimination effects were observed for both faces and hands at mid-parietal sites around 45-85 ms latency, in the incidental task only. This difference was tentatively assumed to reflect an early mechanism of coarse visual categorization. Finally, intentional (but not incidental) gender processing affected the ERPs during a later epoch starting from approximately 200 ms and ending at approximately 250 ms for faces, and approximately 350 ms for hands. This later effect might be related to attention-based gender categorization or to a more general categorization activity.

ERP manifestations of processing printed words at different psycholinguistic levels: time course and scalp distribution.

The aim of the present study was to examine the time course and scalp distribution of electrophysiological manifestations of the visual word recognition mechanism. Event-related potentials (ERPs) elicited by visually presented lists of words were recorded while subjects were involved in a series of oddball tasks. The distinction between the designated target and nontarget stimuli was manipulated to induce a different level of processing in each session (visual, phonological/phonetic, phonological/lexical, and semantic). The ERPs of main interest in this study were those elicited by nontarget stimuli. In the visual task the targets were twice as big as the nontargets. Words, pseudowords, strings of consonants, strings of alphanumeric symbols, and strings of forms elicited a sharp negative peak at 170 msec (N170); their distribution was limited to the occipito-temporal sites. For the left hemisphere electrode sites, the N170 was larger for orthographic than for nonorthographic stimuli and vice versa for the right hemisphere. The ERPs elicited by all orthographic stimuli formed a clearly distinct cluster that was different from the ERPs elicited by nonorthographic stimuli. In the phonological/phonetic decision task the targets were words and pseudowords rhyming with the French word vitrail, whereas the nontargets were words, pseudowords, and strings of consonants that did not rhyme with vitrail. The most conspicuous potential was a negative peak at 320 msec, which was similarly elicited by pronounceable stimuli but not by nonpronounceable stimuli. The N320 was bilaterally distributed over the middle temporal lobe and was significantly larger over the left than over the right hemisphere. In the phonological/lexical processing task we compared the ERPs elicited by strings of consonants (among which words were selected), pseudowords (among which words were selected), and by words (among which pseudowords were selected). The most conspicuous potential in these tasks was a negative potential peaking at 350 msec (N350) elicited by phonologically legal but not by phonologically illegal stimuli. The distribution of the N350 was similar to that of the N320, but it was broader and including temporo-parietal areas that were not activated in the "rhyme" task. Finally, in the semantic task the targets were abstract words, and the nontargets were concrete words, pseudowords, and strings of consonants. The negative potential in this task peaked at 450 msec. Unlike the lexical decision, the negative peak in this task significantly distinguished not only between phonologically legal and illegal words but also between meaningful (words) and meaningless (pseudowords) phonologically legal structures. The distribution of the N450 included the areas activated in the lexical decision task but also areas in the fronto-central regions. The present data corroborated the functional neuroanatomy of word recognition systems suggested by other neuroimaging methods and described their timecourse, supporting a cascade-type process that involves different but interconnected neural modules, each responsible for a different level of processing word-related information.

A ring-shaped distribution of dipoles as a source model of induced gamma-band activity.

As opposed to slow waves, spontaneous and stimulus-induced oscillations in the gamma-band show no polarity reversal in cortical depth, which cannot be explained by the classical equivalent current dipole model usually proposed as a model of pyramidal cell synaptic activity. Here we propose a ring-shaped distribution of dipoles as a source model for these fast oscillations. This distribution generates a field potential that does not reverse through cortical depth. Such a geometry could correspond to horizontally oriented dendritic fields. Moreover, this distribution generates a potential field, but no, or weak, magnetic field on the scalp surface, which corresponds to the observation that visually-induced gamma-band oscillations are detectable in EEG data, but not in simultaneously recorded MEG data.

An evaluation of dipole reconstruction accuracy with spherical and realistic head models in MEG.

MEG forward problem has been solved for about 2000 dipoles placed on the brain surface using a very fine 3-layer realistic model of the head and the boundary element method (BEM). For each dipole, spherical models, one-layer realistic BEM models and coarser 3-layer realistic BEM models, were used to reconstruct the dipole. It was found that the localization bias induced by using a spherical model of the head increased from 2.5 mm in the upper part of the head to 12 mm in the lower part, on average. It was also found that, for the same computing time, a 3-layer model of the head gave on average 2 mm better localization errors than a one-layer model of the head. Orientation errors of less than 20 degrees could only be retrieved with a 3-layer realistic model. Localization and orientation errors highly depended on the dipole position in the brain.

Induced gamma-band activity during the delay of a visual short-term memory task in humans.

It has been hypothesized that visual objects could be represented in the brain by a distributed cell assembly synchronized on an oscillatory mode in the gamma-band (20-80 Hz). If this hypothesis is correct, then oscillatory gamma-band activity should appear in any task requiring the activation of an object representation, and in particular when an object representation is held active in short-term memory: sustained gamma-band activity is thus expected during the delay of a delayed-matching-to-sample task. EEG was recorded while subjects performed such a task. Induced (e.g., appearing with a jitter in latency from one trial to the next) gamma-band activity was observed during the delay. In a control task, in which no memorization was required, this activity disappeared. Furthermore, this gamma-band activity during the rehearsal of the first stimulus representation in short-term memory peaked at both occipitotemporal and frontal electrodes. This topography fits with the idea of a synchronized cortical network centered on prefrontal and ventral visual areas. Activities in the alpha band, in the 15-20 Hz band, and in the averaged evoked potential were also analyzed. The gamma-band activity during the delay can be distinguished from all of these other components of the response, on the basis of either its variations or its topography. It thus seems to be a specific functional component of the response that could correspond to the rehearsal of an object representation in short-term memory.

Human cortical responses evoked by dichotically presented tones of different frequencies.

Behavioral and patient studies have suggested that during dichotic listening the ipsilateral auditory pathways are strongly inhibited, so that each hemisphere is treats the sound coming to the contralateral ear. We analysed the auditory N100m neuromagnetic evoked response following passive listening of dichotic tones of different frequencies. We found that the N100m in each hemisphere depended on both ipsilateral and contralateral stimuli, revealing no strong inhibition of ipsilateral pathways. The N100m increased with the interaural frequency disparity and was reduced as both ears received identical stimuli. The results can be explained by the existence of a frequency-dependent excitatory/inhibitory organization of the auditory cortex, as has been described in the cat. We suggest that the N100m might also reflect automatic processes involved in multiple-stream perception.

Dynamics of MLAEP changes in midazolam-induced sedation.

This study aimed at assessing the effects of midazolam (MDZ) sedation on auditory brainstem (BAEP) and middle latency (MLAEP) evoked potentials in intensive care conditions. Ten ventilated comatose patients were receiving an intravenous MDZ bolus dose (0.2 mg/kg) followed by a 2 h continuous infusion (0.1 mg/kg/h). MLAEPs and BAEPs elicited by clicks (90 dB HL + masking) were simultaneously and continuously monitored during the first 6 h and for 30 min the next morning. We found no effect of MDZ sedation on BAEPs. Only MLAEP components were modified. However, none of the patients presented any total abolition of the MLAEPs. Bolus injection led to very early alteration of cortical responses, beginning after 5 min and lasting almost 1 h (maximum Pa latency increase, 3.1 ms; maximum Pa-Nb amplitude decrease, 46%). During continuous infusion, MLAEPs remained slightly, although significantly, altered (Pa latency, +1.3 ms; Pa-Nb amplitude, 27%). The Nb wave seemed to be modified earlier and to return to normality later than the Pa wave. These findings incite a careful interpretation of MLAEP tracings acquired during the first hour following MDZ bolus injection. If possible, MDZ should be administered as continuous infusion for reliable interpretation of evoked potential changes in intensive care unit, or during surgery.

Effects of altitude hypoxia on middle latency auditory evoked potentials in humans.

BACKGROUND: This paper presents an investigation of auditory evoked responses in humans subjected to high altitude hypoxic conditions. METHODS: Middle latency (MLAEPs) as well as short latency (BAEPs) evoked potentials were recorded in 10 healthy subjects, first at sea level (N), then 24 h (H1) and 72 h (H3) after their arrival at an altitude of 4350 m. At the same time, arterial blood parameters (PaO2, PaCO2 and pH) were measured and the clinical status of the subjects was assessed. RESULTS: In altitude conditions, the amplitude of BAEP peak V decreased (-17%). The MLAEP waves showed variations in the shape of their latest waveforms. Their amplitudes, however, were not affected. The Pa-Nb interpeak latency significantly decreased (-2.2 ms) between N and H1, and remained stable during the stay at high altitude. CONCLUSION: A correlation was found between the relative decrease of PaCO2 and the shortening of Nb wave latency, suggesting that the variations in MLAEPs could be preferentially related to the ventilatory response of the subjects in hypoxic conditions. However, no correlation was found between the clinical status of the subjects (Acute Mountain Sickness score) and the parameters of the waves.

A systematic evaluation of the spherical model accuracy in EEG dipole localization.

This paper presents a study of the intrinsic localization error bias due to the use of a spherical geometry model on EEG simulated data obtained from realistically shaped models. About 2000 dipoles were randomly chosen on the segmented cortex surface of a particular subject. Forward calculations were performed using a uniformly meshed model for each dipole located at a depth greater than 20 mm below the brain surface, and locally refined models were used for shallower dipoles. Inverse calculations were performed using four different spherical models and another uniformly meshed model. It was found that the best spherical model lead to localization errors of 5-6 mm in the upper part of the head, and of 15-25 mm in the lower part. The influence of the number of electrodes upon this intrinsic bias was also studied. It was found that using 32 electrodes instead of 19 improves the localization by 2.7 mm on average, while using 63 instead of 32 electrodes lead to improvements of less than 1 mm. Finally, simulations involving two simultaneously active dipoles (one in the vicinity of each auditory cortex) show localization errors increasing by about 2-3 mm.

Matching of digitised brain atlas to magnetic resonance images.

A method has been developed to match a standard digitised brain atlas onto MR images for identification of cerebral structures in anatomical images. This method uses, first, a three-dimensional crude registration based on the proportional system of Talairach. Then, a two-dimensional refined registration is performed using a deformation function based on a set of homologous landmarks on both images (MR and atlas). Displacements vectors are computed between each corresponding landmark. These vectors are interpolated by thin-plate splines, generating an unwarping function defined on the whole image. This function can then be applied on any structure of the atlas. An evaluation of the matching procedure has been performed. First, the influence of the choice of the landmarks has been evaluated for the fine registration method. The latter has been then compared to the crude registration method considered as a classical reference method. These results show the advantages of the fine registration approach.

Improved dipole localization using local mesh refinement of realistic head geometries: an EEG simulation study.

A systematic evaluation of dipole localization accuracy using the boundary element method is presented. EEG simulations are carried out with dipoles located in the right parietal and temporal regions of the head. Uniformly meshed and locally refined head models are considered in both spherical and realistic geometries. An initial study determines the influence upon the localization accuracy of the dipole depth below the brain surface, of its orientation (radial and tangential), and of the global and local mesh densities. Simulated potential data are computed analytically in the spherical case, and numerically using a very fine (locally refined) model in the realistic case. Results in both geometries show that in order to get localization errors of about 2-4 mm, uniformly meshed models may be used for dipoles located at depths greater than 20 mm, whereas locally refined models should be used for shallower dipoles. Two other studies show how the localization accuracy depends upon the size of the local refinement area and upon the number of electrodes (19, 32, 63). Results show that a large number of electrodes brings significant improvements, especially for shallow and tangential dipoles.

Stimulus specificity of phase-locked and non-phase-locked 40 Hz visual responses in human.

Considerable interest has been raised by non-phase-locked episodes of synchronization in the gamma-band (30-60 Hz). One of their putative roles in the visual modality is feature-binding. We tested the stimulus specificity of high-frequency oscillations in humans using three types of visual stimuli: two coherent stimuli (a Kanizsa and a real triangle) and a noncoherent stimulus ("no-triangle stimulus"). The task of the subject was to count the occurrences of a curved illusory triangle. A time-frequency analysis of single-trial EEG data recorded from eight human subjects was performed to characterize phase-locked as well as non-phase-locked high-frequency activities. We found in early phase-locked 40 Hz component, maximal at electrodes Cz-C4, which does not vary with stimulation type. We describe a second 40 Hz component, appearing around 280 msec, that is not phase-locked to stimulus onset. This component is stronger in response to a coherent triangle, whether real or illusory: it could reflect, therefore, a mechanism of feature binding based on high-frequency synchronization. Because both the illusory and the real triangle are more target-like, it could also correspond to an oscillatory mechanism for testing the match between stimulus and target. At the same latencies, the low-frequency evoked response components phase-locked to stimulus onset behave differently, suggesting that low- and high-frequency activities have different functional roles.

Improved forward EEG calculations using local mesh refinement of realistic head geometries.

A method for semi-automatically constructing realistic surface meshes of 3 head structures--scalp, skull and brain--from a stack of MR images is described. Then an evaluation is given for both spherical and realistic dipolar models, using the boundary element method (BEM). In both cases, locally refined models were considered. Two characteristic mesh parameters were defined: the global and the local mesh densities (in triangles per cm2). In spherical geometries, numerical and analytical solutions were compared, and in the realistic case, all models were compared to a highly refined one, considered as a reference. Both geometries gave comparable results. It was found that for "deep dipoles" located at more than 20-30 mm under the brain surface, meshes with a global density of 0.5 tri/cm2 gave "acceptable" results, whereas for more superficial dipoles (2-3 mm < depth < 20-30 mm), it was necessary to locally refine meshes near the source location up to a local density of about 5-8 tri/cm2, to get comparable results.

Gamma-range activity evoked by coherent visual stimuli in humans.

We tested the hypothesis of a role of gamma-range synchronized oscillatory activity in visual feature binding by recording evoked potentials from 12 subjects to three stimuli: two coherent ones (a Kanizsa triangle and a real triangle) and a non-coherent one (a Kanizsa triangle in which the inducing disks had been rotated so that no triangle could be perceived). The evoked potentials were analysed by convoluting the signal for each subject and each stimulation type by Gabor wavelets centred from 28 up to 46 Hz, providing a continuous measure of frequency-specific power over time. A first peak of activity was found around 38 Hz and 100 ms with a maximum at electrode Cz in each experimental condition. A second peak of activity occurred around 30 Hz and 230 ms, with a maximum at O1 in response to the real triangle and a maximum at Cz in the case of the illusory triangle. At 100 ms we did not find any variations of the gamma-band component of the evoked potential with stimulation type, but the power of the 30 Hz component of the evoked potential between 210 and 290 ms differed from noise only in the case of a coherent triangle, no matter whether real or illusory. We thus found a 30 Hz component whose power correlates with stimulus coherency, which supports the hypothesis of a functional role of high-frequency synchronization in feature binding.(ABSTRACT TRUNCATED AT 250 WORDS)

Tonotopic organization of the human auditory cortex: N100 topography and multiple dipole model analysis.

The tonotopic organization of the human auditory cortex has been investigated by means of scalp potential mapping and dipole modelling of the evoked response occurring around 100 msec after the stimulus onset. The major characteristics of the topographical changes observed with increasing stimulus frequency were statistically demonstrated. Using a 3-concentric sphere head model, the scalp potential distributions can be explained in first approximation by two equivalent current dipoles, located in the supratemporal plane and mimicking the activity of both auditory cortices. To take into account the temporal aspects of the brain activities, 3 time-varying dipole strategies were tested. Frequency dependence of the dipole orientation has been evidenced in both hemispheres with the 3 models, whereas no significant change in dipole position was found. The tilt in dipole orientation could be related to the folding geometry of Heschl's gyrus, which varies with depth. In agreement with previous MEG findings, this brings new evidence for a tonotopic organization of the auditory cortical area involved in the N100 wave generation. Moreover, distinct frequency dependences of the equivalent current dipoles were observed in the early and the late parts of the N100. This study demonstrates that simple dipolar models, applied on electrical data, make it possible to reveal functionally distinct cortical areas.

Separate representation of stimulus frequency, intensity, and duration in auditory sensory memory: an event-related potential and dipole-model analysis.

Abstract The present study analyzed the neural correlates of acoustic stimulus representation in echoic sensory memory. The neural traces of auditory sensory memory were indirectly studied by using the mismatch negativity (MMN), an event-related potential component elicited by a change in a repetitive sound. The MMN is assumed to reflect change detection in a comparison process between the sensory input from a deviant stimulus and the neural representation of repetitive stimuli in echoic memory. The scalp topographies of the MMNs elicited by pure tones deviating from standard tones by either frequency, intensity, or duration varied according to the type of stimulus deviance, indicating that the MMNs for different attributes originate, at least in part, from distinct neural populations in the auditory cortex. This result was supported by dipole-model analysis. If the MMN generator process occurs where the stimulus information is stored, these findings strongly suggest that the frequency, intensity, and duration of acoustic stimuli have a separate neural representation in sensory memory.

Two separate frontal components in the N1 wave of the human auditory evoked response.

Scalp current density analysis of the auditory evoked response to 1-kHz tone bursts delivered at various interstimulus intervals (ISIs) (from 1 s to 2 min in separate runs) shows that two different frontal components can be observed and functionally dissociated in the N1 time range: one is elicited for all ISIs, peaks at about 95 ms poststimulus, and has a full recovery time below 8 s; the second is elicited only by infrequent stimuli (ISIs > 4 s), peaks around 140 ms, and significantly increases in amplitude with increasing ISIs. The first component can be considered a new obligatory component in N1 elicited simultaneously with the responses in auditory cortex; the later component could correspond to the orienting Component III of Näätänen and Picton (1987).