Skip Navigation

This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (31)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Halgren, E.
Right arrow Articles by Chauvel, P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Halgren, E.
Right arrow Articles by Chauvel, P.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Cerebral Cortex, Vol. 12, No. 7, 710-728, July 2002
© 2002 Oxford University Press

Rapid Distributed Fronto-parieto-occipital Processing Stages During Working Memory in Humans

E. Halgren,1,2, C. Boujon,3, J. Clarke, C. Wang,1 and P. Chauvel,2

Institut National de la Santé et de la Recherche Médicale, CJF90-12, Neurologie, CHU Pontchaillou, F-35033 Rennes, France

E. Halgren, MGH NMR Center, 149 13th Street, Charlestown, MA 02129, USA. Email: halgren{at}nmr.mgh.harvard.edu.


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Cortical potentials were recorded from implanted electrodes during a difficult working memory task requiring rapid storage, modification and retrieval of multiple memoranda. Synchronous event-related potentials were generated in distributed occipital, parietal, Rolandic and prefrontal sites beginning ~130 ms after stimulus onset and continuing for >500 ms. Coherent phase-locked, event-related oscillations supported interaction between these dorsal stream structures throughout the task period. The Rolandic structures generated early as well as sustained potentials to sensory stimuli in the absence of movement. Activation peaks and phase lags between synaptic populations suggested that perceptual processing occurred exclusively in the visual association cortex from ~90 to 130 ms, with its results projected to fronto-parietal areas for interpretation from ~130 to 280 ms. The direction of interaction then appeared to reverse from ~300 to 400 ms, consistent with mental arithmetic being performed by fronto-parietal areas operating upon a visual scratch pad in the dorsolateral occipital cortex. A second reversal, from ~420 to 600 ms, may have represented an updating of memoranda stored in fronto-parietal sites. Lateralized perisylvian oscillations suggested an articulatory loop. Anterior cingulate activity was evoked by feedback signals indicating errors. These results indicate how a fronto-centro-parietal ‘central executive’ might interact with an occipital visual scratch pad, perisylvian articulatory loop and limbic monitor to implement the sequential stages of a complex mental operation.


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Working memory maintains information for use in on-line cognition. As conceptualized by Baddeley, working memory is composed of a central executive controlling ‘slave systems’, primarily a ‘visuospatial scratch pad’ and an ‘articulatory loop’ (Baddeley, 1986Go). Working memory thus combines a mental workspace with short-term buffers and executive processes for interpreting incoming stimuli, updating actively stored items accordingly, and maintaining the intermediate results of processing for further calculations. Working memory has become a key concept in many models of higher cognitive activities such as reading (Just et al., 1996Go). The amount of resources in working memory and the efficiency of their dynamic allocation are seen as the main factors underlying limitations in fundamental intellectual abilities (Kyllonen and Christal, 1990Go).

In some working memory tasks, a single stimulus is presented, then retained as presented for ~2–30 s and finally retrieved. Such tasks probe what has been termed primary memory (James, 1890Go) and are amenable to neurophysiological analysis because the input, retention and retrieval stages are separated temporally. In contrast, other working memory tasks require a rapid interaction between perception, interpretation and memory. Previous studies have used positron emission tomography (PET) or functional magnetic resonance imaging (fMRI) for measuring haemodynamic or metabolic activation during working memory. Both varieties of working memory task activate an extended network of the prefrontal, parietal and medial frontal cortices (Swartz et al., 1995Go; Fiez et al., 1996Go; Smith et al., 1998Go; Ungerleider et al., 1998Go).

Haemodynamic and neuropsychological studies suggest that it may be possible to distinguish the contributions that these different areas make to working memory. Consistent with a long tradition in neuropsychology, the prefrontal cortex has been associated with executive function (Baddeley, 1986Go). In some haemodynamic studies, the more ventral prefrontal areas have been associated with the simple retention of material in primary memory, whereas the more dorsolateral areas have been engaged when the task requires a higher level of executive control (Petrides et al., 1993Go; D'Esposito et al., 1995Go; Cohen et al., 1997Go; Owen, 1997Go; Smith et al., 1998Go). Other studies have found the location of prefrontal activation to depend upon the nature of the material to be retained, with verbal/non-verbal material mapping to the dominant/non-dominant hemisphere and spatial/ non-spatial material mapping to the dorsal/ventral areas (Smith et al., 1998Go; Ungerleider et al., 1998Go). Similarly, lesion studies have suggested that the parietal cortex (in particular near the intraparietal sulcus) is specialized in spatial analysis (Critchley, 1953Go). Conversely, the left perisylvian activation is associated with the maintenance of information in the articulatory loop according to both neuropsychological (Shallice and Vallar, 1990Go) and haemodynamic evidence (Fiez et al., 1996Go). Finally, anterior cingulate activation has been associated with error monitoring and stimulus response mapping (Bush et al., 2000Go).

These studies support hypotheses associating different cortical areas with different working memory components. However, the task/area double dissociations predicted by these hypotheses are often difficult to obtain. In particular, parietal and anterior cingulate activations in working memory tasks tend to be correlated with prefrontal activation across task conditions (Carpenter et al., 2000Go; Diwadkar et al., 2000Go) and a strict correlation of working memory deficits with prefrontal lesions has not been found (Frisk and Milner, 1990Go). Co-activation of a parieto-fronto-cingulate network has also been found in other tasks, in particular those involving attention (Corbetta, 1998Go) and task difficulty/self-monitoring/response choice (Cohen et al., 2000Go). Similarly, neuroanatomical and neuropsychological evidence supports an integrated parieto-fronto-cingulate circuit in executive functions (Goldman-Rakic, 1988Go; Mesulam, 1990Go). Finally, efforts to isolate cortical areas that are exclusively activated by the necessity of controlled allocation of attentional resources in dual task situations have yielded inconsistent results, suggesting that high-level executive functions may arise through the interaction of multiple structures (Adcock et al., 2000Go).

One interpretation of these data is that, although the parietal and prefrontal cortices make distinct contributions to complex processing (as revealed by lesions), they are so interconnected that their activation patterns are nearly indistinguishable. Alternatively, it is possible that the different structures not only have distinct functions but that they are activated sequentially and this cannot be seen with haemodynamic measures because of their poor temporal resolution (>~1 s) (Buxton et al., 1998Go; Dale and Halgren, 2001Go). That is, although these areas are coactivated on a time-scale of ~1 s and belong to the same general system, it is conceivable that their activation is modular and sequential rather than continuously interactive and overlapping.

The spatiotemporal dynamics of cortical activation during working memory also have implications for whether the prefrontal cortex acts in working memory as a repository of primary memory, as opposed to as the central executive. As primary memory, specific and sustained firing of prefrontal neurons would encode significant memoranda as well as the cognitive context. This sustained firing would function not only as a memory store but also as a source of task-specific instructions and temporal markers for the posterior association cortex, thus freeing it from immediate sensory input (Goldman-Rakic, 1988Go; Fuster, 1989Go; Halgren, 1994Go). Alternatively, as central executive, the prefrontal cortex would provide goal-directed sequencing and control of the posterior association cortex. Initially, this control directs attention to a sensory channel by activating the associated cortical area (Desimone, 1996Go; Hillyard and Anllo-Vento, 1998Go; Mangun et al., 1998Go; Knight et al., 1999Go). Subsequent executive control can be conceptualized as an extension of this process, with the prefrontal cortex successively activating cortical areas concerned with perceptual processing, interpretation of the stimulus with respect to task instructions, consequent mental operations on active memoranda, memory updating and response selection. These central executive functions would imply phasic activation of prefrontal neurons corresponding to the successive processing stages.

A related question is whether primary memory is held in the prefrontal cortex or the posterior sensory association cortex or both. If it is only in the prefrontal cortex, then sensory association areas may complete their work soon after stimulus presentation, pass the results to anterior areas and then no longer be active. Conversely, the model mentioned above (Baddeley, 1986Go) suggests that sensory association areas are involved in later stages simultaneously with higher cortical areas, in particular in the prefrontal cortex.

In summary, an adequate neural model of working memory should specify not only which structures are involved, but also the dynamics of their interaction. The most extensive source of data regarding these dynamics is unit-recording studies in macaques during working memory (Goldman-Rakic, 1988Go). Sustained stimulus-specific firing by prefrontal neurons supports the hypothesized prefrontal contribution to primary memory. However, phasic firing patterns are also observed in prefrontal neurons and, conversely, the same categories of task-related firing are observed in temporal and parietal sites, with the time-course of activation in these areas entirely overlapping with prefrontal activation (Fuster, 1989Go; Chafee and Goldman-Rakic, 1998Go). Furthermore, the parietal and prefrontal firing patterns are modified but fundamentally maintained during cooling of either location, suggesting that they arise out of the coordinated activity of an extended network (Fuster et al., 1985Go; Chafee and Goldman-Rakic, 2000Go). The lack of stronger effects is interpreted as being due to the high degree of redundancy in connections in the generating circuit (Goldman-Rakic, 1988Go). However, the similarity of activation in different areas that seems to reflect their seamless integration leaves no leverage for explicating their distinct roles (Chafee and Goldman-Rakic, 2000Go). Furthermore, it is uncertain how to generalize these results to humans in more complicated working memory tasks that engage the central executive heavily.

Direct electrophysiological recordings in humans during working memory would thus be useful in understanding the temporal dynamics of engagement by different cortical areas. Such recordings with scalp electrodes have also been interpreted as indicating continuing interactions between the anterior and posterior cortex (Gevins et al., 1996Go). However, the cortical generators of these signals cannot be inferred unambiguously (Hamalainen et al., 1993Go; Dale and Halgren, 2001Go). Intracranial electroencephalography (iEEG) recordings have high temporal resolution and under some circumstances can also localize generators with high spatial resolution (Halgren et al., 1980Go). Such recordings can only be obtained from electrodes implanted for clinical purposes and, thus, are limited by possible contamination with pathological activity, as well as clinical limitations in spatial sampling. Nonetheless, they offer a view of the spatiotemporal activation patterns that underlie cognitive events, which is complementary to non-invasive measures.

In the current study, we report the results of intracranial recordings during Mental Counters, a working memory task that requires the rapid perception and evaluation of six successive stimuli interleaved with the maintenance and updating of active memoranda (Larson and Saccuzzo, 1989Go). The memoranda were digits or figures in different variants that were intended to probe different material-specific ‘scratch pads’. Following the input stimuli, a probe stimulus was presented and then a feedback tone following the response.

The results are most consistent with a sequential but not modular model of cortical function in working memory. Although early processing was briefly confined to the visual association cortex, activation soon spread to multiple occipital, parietal and frontal sites, which all remained active for the entire epoch. In particular, the co-activation of visual association with the fronto-parietal cortices suggests that it contributes to working memory beyond perceptual analysis. The results support the presence of at least four processing stages in working memory, with each lasting ~50–200 ms. Phase-locked oscillations from ~4 to 12 Hz were prominent in multiple structures including the prefrontal cortex, as would be expected if it contributes to the central executive. Single-trial spectral analysis of the dominant oscillations suggested organized interactions between the areas, with rapidly shifting directions of information flow. Event-related oscillations ranged from 4 to 50 Hz. They were commonly embedded in slower oscillations synchronized with the individual stimuli (~1 Hz), which in turn were embedded in sustained activation from the first stimulus to the probe (<0.1 Hz). The latter supports a primary memory contribution of the prefrontal cortex to working memory. Limbic areas, in particular the anterior cingulate gyrus, were preferentially activated by the probe and feedback stimuli. Overall, the results reveal multiple embedded oscillations that may reflect the temporal organization of processing and memory functions across multiple cortical areas during this difficult task.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Patients

Seventeen patients (10 females and seven males) suffering from intractable partial epilepsy consented to participate in this study while undergoing iEEG for possible neurosurgical treatment. All subjects except one were right-handed and that subject was left dominant for language according to the Wada test. Their mean age was 26.2 ± 7.7 years.

Behavioral Tasks

The subjects performed different variants of the Mental Counters paradigm of Larson, which were modified to be appropriate for neurophysiological recordings in patients (Larson and Saccuzzo, 1989Go). These tests all required the maintenance of different memoranda and their updating in response to rapidly presented stimuli. In the standard numeric paradigm with two memoranda presented in the visual modality (N2V), the location of the stimuli indicated which memorandum was to be operated on and which operation was to be performed on it (Fig. 1Go). Based on preliminary testing with a control group, NV2 was chosen as the basic task variant to be given to all patients. Of the 17 subjects, 16 performed the N2V variant. Some subjects also performed variants of the Mental Counters paradigm that differed from N2V in material (using figures or letters instead of numbers, i.e. F2V or L2V), modality (auditory rather than visual, i.e. N2A) or memory load (three rather than two memoranda, i.e. N3V). Note that, in all of the visual variants, identical stimuli were used and the tasks were different only in the nature of the memoranda and the operations to be performed upon them. Of the 17 subjects, 10 performed the F2V variant, four the L2V variant, three the N3V variant and only one the N2A variant. Each trial consisted of six stimuli, 36 trials were presented per block and one to four blocks of a variant were obtained in a particular subject. In addition to electrophysiological recordings, reaction times and response accuracy were obtained.



View larger version (24K):
[in this window]
[in a new window]
 
Figure 1. The Mental Counters task. Examples of two variants of mental counters are shown. Other variants using three locations or letters as memoranda were also used, as described in the text.

 
The study was described to the subject in the subject's hospital room, usually on the day prior to recordings and informed consent was requested. The subjects then participated in a training session of 45 min during which each variant was initially presented with inter-stimuls intervals (ISIs) of 2500 ms. This was gradually decreased to 1240 ms and finally 960 (nine subjects) or 840 ms (eight subjects) if the patient was capable. Approximately half of the patients receiving training were unable to perform the task at an acceptable speed, presumably due to psychomotor slowing induced by anticonvulsant medications. Physiological recordings were only obtained from the remaining 17 subjects.

Numbers–Two Memoranda–Visual (N2V)

At the beginning of each trial the initial values of the memoranda (three zeros) were presented on the left, right and middle of a 16 in colour monitor. When the subject pressed a key, three zeros were replaced in the same positions by three points. Then at ISIs of 840 ms (240 ms stimulus exposure) a short line appeared. If the line was horizontal (42%) or vertical (42%) then the subject made a subtraction or addition to either the left or right memorandum, according to the location of the target. If the line was above a dot then the operation was addition of one and if the line was below a dot then the operation was subtraction. The new result was maintained in memory to be used as a basis for the next line. Diagonal lines (16%) indicated that no operation was to be performed (i.e. ‘catch’ trials). Six lines were presented in each trial. The location of the lines and their orientation varied randomly on successive presentations with the exception that the correct value of any given memorandum was maintained between –3 and +3. At the end of these six presentations a proposed solution consisting of three numbers was provided on the screen, one at each of the original dot locations. Fifty percent of the proposed solutions matched the correct result. The subject had to press the left key in <3 s if the proposed solution was correct and the right key if it was not. A feedback tone indicated whether the subject's answer was correct (high tone) or not (low tone). Throughout the entire trial the subject maintained fixation on the centre of the monitor, which was occupied by either a dot or a number. The response and feedback procedures as well as all time-intervals and stimulus-type proportions were identical across all variants.

Numbers–Three Memoranda–Visual (N3V)

N3V was identical to N2V, except that the lines could appear above or below the middle point on the monitor, i.e. three rather than two memoranda needed to be maintained in memory and incremented or decremented accordingly.

Numbers–Two Memoranda–Auditory (N2A)

N2A used sounds as the stimuli. As in N2V each trial started with the presentation of three zeros on the left, right and centre of the monitor. The dots then remained on the screen until they were replaced with the response probe. Throughout this period the subject maintained fixation on the centre dot. Adding one to the memorandum was signalled by pure high-pitched sounds, subtraction by low-pitched sounds and catch trials by complex sounds. In order to indicate that the operation was to be carried out on the right-side memorandum, the sound was presented in the right ear and vice versa. After six sounds a solution was proposed on the screen.

Letters–Two Memoranda–Visual (L2V)

L2V used letters as the memoranda. Each trial began with three Ls on the left and right and in the centre of the monitor. These were replaced by dots and then short lines started appearing above or below the left or right dot, exactly as in N2V. Horizontal or vertical lines above dots signalled going ‘up’ a letter (or backwards) in the alphabet (e.g. L to K) and targets below dots signalled going ‘down’ a letter (e.g. L to M). Again, diagonal lines signalled no operation. After six lines were processed a response probe consisting of three letters was presented.

Figures–Two Memoranda–Visual (F2V)

F2V used two figures (left or right) as the memoranda, each consisting of a horizontal line, a vertical line or both (i.e. a plus). Each trial began with three vertical lines on the left, right and centre of the monitor. These were replaced by dots, and then short lines started appearing above or below the left or right dot, exactly as in N2V. The appearance of a horizontal or vertical line above a dot signalled that the corresponding element should be added to the corresponding figure (Fig. 1Go). Similarly, the appearance of a horizontal or vertical line below a dot signalled that the corresponding element should be removed. Diagonal lines signalled no operation. After six lines were processed a response probe consisting of three figures was presented.

Stereoencephalographic Recordings and Stimulation

Stereoencephalographic (SEEG) recordings were obtained from a total of five to eight multicontact probes surgically inserted into the depths of one or both cerebral hemispheres of each patient (Fig. 2Go). The anatomical location of these probes was decided entirely on clinical grounds for localizing the origins of the patients' epileptic seizures (Chauvel et al., 1996aGo). Typically, the probes remained implanted for 4–7 days in order to obtain recordings during seizure activity. The recordings in the present study were obtained during seizure-free periods. Each probe contained between five and 15 stainless steel contacts, thereby allowing for recordings from different medial to lateral sites aligned perpendicular to the sagittal plane (Fig. 2Go). The contacts were 2.0 mm in length and were separated by 1.5 mm. Recordings were obtained across all 17 subjects from 667 left hemisphere and 701 right hemisphere sites. The contacts were localized using human stereotaxic statistical data (Talairach et al., 1967Go; Talairach and Tournoux, 1988Go) and confirmed with stereoscopic stereotaxic angiography (Szikla et al., 1977Go) and, in most cases, stereotaxic MRI (Musolino et al., 1990Go). The probes are lettered (primed in the left hemisphere). The contacts are numbered according to their distance from the midline in millimetres. The contact locations are also indicated in Talairach's coordinates (x = distance to the right of the midline, y = distance forward from the anterior commissure or AC and z = height above the anterior commissure–posterior commissure or AC–PC line) (Talairach and Tournoux, 1988Go). The actual coordinates are listed first, then normalized in italics. Since the y and z coordinates are constant for all contacts of a given probe, they are omitted where they are clear from the context.



View larger version (40K):
[in this window]
[in a new window]
 
Figure 2. Electrode locations. The probe entry points in the left and right hemispheres (LH and RH) are normalized following an earlier study (Talairach and Tournoux, 1988Go). Recording contacts are located every 3.5 mm (centre to centre) along the probe's lateral to medial trajectory orthogonal to the sagittal plane.

 
In many cases, the contacts were stimulated electrically in order to delineate the origin and spread of seizure discharges better (Chauvel et al., 1996aGo). Unless otherwise noted, these stimulations were bipolar between adjacent electrode contacts, with pulses of 1 ms delivered in a train of 5 s at 50 pulses/s at stimulus intensities below that necessary for evoking a local after-discharge (as monitored by simultaneous SEEG recordings). Those stimulations with effects that provided confirmation of functional localizations are presented below.

The contacts were located in the following regions (the counts include contacts located in the underlying white matter): the occipital cortex (n = 62), precuneus and superior parietal lobule (n = 15), posterior cingulate gyrus (n = 60), supramarginal gyrus (n = 128), lingual and fusiform gyri at the occipitotemporal junction (n = 30), middle temporal gyrus at the occipitotemporal junction (n = 48), posterior hippocampal formation (n = 70), posterior middle temporal gyrus (n = 105), anterior hippocampal formation (n = 80), middle level of the middle temporal gyrus (n = 117), amygdala region (n = 75), anterior level of the middle temporal gyrus (n = 107), posterior superior temporal gyrus (n = 80), anterior superior temporal gyrus (n = 18), ventromedial prefrontal cortex (n = 71), ventrolateral prefrontal cortex (n = 157), dorsomedial prefrontal cortex (n = 15), ventrolateral prefrontal cortex (n = 24), paracentral gyrus (n = 19) and Rolandic cortex (n = 87) (including the precentral and postcentral gyri and sulci). Scalp recordings were usually obtained at Cz or Pz. The specifics of localization are presented in the context of descriptions of individual responses.

Recordings were digitized on-line with an accuracy of 12 bits. In the initial eight subjects, recordings were made as separate epochs during the six input (I1–I6), probe and feedback stimuli at 6 ms per sample with a bandpass of 0.1–35 Hz. The epochs began 40 ms before the stimulus and ended 794 or 914 ms after the stimulus (for the 840 or 960 ms stimulus-to-stimulus onset intervals, ISIs, respectively) for the I1–I6 and probe stimuli. The epochs began 190 ms before the stimulus and ended 644 or 764 ms after the stimulus for the feedback stimuli. Data were acquired continuously in the last 10 subjects during the entire task at 200 Hz and a bandpass of 0.5–50 Hz. This bandpass may have reduced sustained potentials related to primary memory. Eye movement artefacts (>75 µV) were monitored with an electrode placed at the outer canthus of the right eye. Unipolar recordings to a relatively inactive reference (the nose) were employed in order to understand the field distributions in terms of both gradients in amplitude as well as absolute levels. The potentials recorded at the nose relative to a distant (i.e. non-cephalic) reference were negligible in amplitude compared to the amplitude of intracranial potentials and, thus, it can be considered an inactive reference in the current study.

Analysis

The recordings were screened for artefacts (e.g. epileptiform activity and eye or body movements) and then averaged within each task variant in order to yield event-related potentials (ERPs) of input stimuli indicating addition, subtraction or no operation (‘catch’), of true and false probes, and of correct versus incorrect responses. In addition, ERPs were made of the entire trial in order to identify longer duration processes. These whole-trial ERPs in the initial eight subjects were made by concatenating the ERPs to all conditions of the I1–I6, probe and feedback epochs. Consequently, the ERPs began 40 ms before the I1 epoch and ended 644 or 764 ms after the feedback epoch, with small gaps between the I6 and probe epochs and between the probe and feedback epochs. Data were acquired continuously in the last 10 subjects, thereby permitting a continuous epoch from before the I1 epoch until after the feedback epoch.

Previous studies have suggested that most task-related EEG activity is not phase locked to the evoking stimulus and thus does not contribute to the ERP (Klopp et al., 1999Go; Tallon-Baudry and Bertrand, 1999Go). Furthermore, averaging obscures the temporal relations between EEGs in different structures. In order to obtain a more comprehensive estimate of task-related responses, the individual trial spectral power, phase locking and phase lag were calculated using a wavelet function for different frequencies and latencies with respect to the stimulus. The resulting time–frequency plots were averaged across trials and normalized with respect to the pre-stimulus baseline. For these measures, the single trial signal si(t) for each channel i was convoluted with complex Morlet's wavelets w(t,f0):

The time–frequency representation TFi(t,f0) has both amplitude [Ai(t,f0)] and phase [{phi}i(t,f0)] information. The estimated spectral amplitude is the average of Ai(t,f0) across all trials. The phase-locking factor is defined (Lachaux et al., 1999Go) as

The phase lag factor is defined as

Each Morlet's wavelet (Kronland-Martinet et al., 1987Go) used here had a Gaussian shape both in the time domain (SD {sigma}t) and in the frequency domain (SD {sigma}f) around its central frequency f0:

with {sigma}f = 1/2{pi}{sigma}t. Wavelets are normalized so that their total energy is unity, the normalization factor A being equal to ({sigma}t{pi}0.5)-0.5. Relatively constant temporal and spatial resolution across target frequencies were obtained by adjusting the wavelet widths (f0/{sigma}f ratio) according to the target frequency. In the initial eight subjects, the wavelet widths increased linearly from 1.5 to 10 as the frequency increased from 4 to 40 Hz, resulting in a {sigma}t value of 60–40 ms and a {sigma}f value of 2.7–4 Hz. In the last 10 subjects, the wavelet widths increased from 2 to 13 as the frequency went from 4 to 50 Hz, resulting in a {sigma}t value of 79–41 ms and a {sigma}f value of 2.0–3.8 Hz. Tests with simulated data confirmed that the methods used here can show clear phase-locking and lag patterns, even at 4 Hz. The findings were confirmed with standard spectral power and coherence measures on successive overlapping windows as described previously (Klopp et al., 1999Go, 2000Go).

Like extracranial EEG and magnetoencephalography (MEG), SEEG has high temporal resolution (limited by the sampling interval) but, unlike them, SEEG also has fine anatomical resolution (limited by the spacing between recording sites). Although recordings are only obtained from epileptic patients, recording locations and epochs may be selected where the EEG appears to be normal. Indeed, when it is possible to compare the results from SEEG with those from non-invasive measures in normal subjects, very similar results are obtained. The main limitation of SEEG is incomplete sampling: in order to sample the entire brain at a spacing of 3.5 mm over 10 000 contacts would be needed, but only ~100 contacts are sampled in each patient. Furthermore, in the tissue that is actively generating the EEG response the amplitude and polarity change rapidly from contact to contact. Thus, it is very difficult to average responses across subjects, even if they are nominally recorded from the same region.

These considerations lead to a data analysis strategy where every recording is examined for rapid changes in amplitude and/or polarity over short distances. Electrode contacts with proof of locally generated ERP components are then further intensively examined for the cognitive correlates of those components (by comparison across other tasks), for the precise location of the contact (anatomically by comparison with MRI and stereotactic arteriography and functionally by comparison with the effects of local electrical stimulation) and for the presence of activation in the spectral domain. When multiple ERP generators in different sites are simultaneously recorded, then evidence for the presence and direction of communication between these sites is sought using phase-locking and phase lag measures. Initial analysis thus treats the data as a series of single-case studies, thereby proving local generation of particular ERP components in particular locations in particular people. The consistency of the data is evaluated by identifying replications across subjects. The certain localization and high spatiotemporal resolution of SEEG can then be used to help in interpreting particular activations identified by fMRI, PET, MEG and/or EEG studies.


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Behaviour

Behavioural performance for the two most common task variants is shown in Table 1Go. On average, the subjects responded correctly in 69% of the trials with a reaction time of 1064 ms to the probe. Performance was slightly better (t-test, P < 0.05) on the version with numerical memoranda (N2V) than that with figural memoranda (F2V), but the reaction times did not differ.


View this table:
[in this window]
[in a new window]
 
Table 1 Behavioural performance
 
Overview

The task evoked a sustained multiphasic response that was distributed across widespread cortical regions. This is illustrated in a single subject with probes sampling the occipital, parietal, Rolandic and medial regions (Figs 3–5GoGoGo). Despite the lack of somatosensory stimuli or behavioural responses during the input stimuli, they evoked reliable oscillations from ~0.05 to 35 Hz in the peri-Rolandic area (Figs 6–8GoGoGo). In contrast, limbic sites, in particular the anterior cingulate gyrus, tended to be activated by the feedback tones indicating incorrect responses (Figs 9 and 10GoGo). Single-trial, phase-locking measures (Figs 5 and 6GoGo) confirmed the similarity of the activity in widespread occipital, parietal, central, medial and prefrontal sites, and phase lag measures suggested shifting patterns of information flow across the multiple processing phases.



View larger version (52K):
[in this window]
[in a new window]
 
Figure 3. Sustained co-activation of the occipital, parietal and Rolandic cortices. Each set of waveforms shows activity recorded from one of 26 different contacts located on five probes recorded simultaneously from a subject performing the Mental Counters task. Each recording begins 40 ms before and ends 800 ms after stimulus onset. Activity is maintained simultaneously in occipital, parietal, central and medial sites for much of the trial. In the upper traces, multiple inversions in polarity are recorded from ~100 to 800 ms by two anterior occipital probes (O and D) (the numbers indicate the distance in millimetres from the centre of the recording contact to the cerebral midline). Several components peaking from 105 to 760 ms post-stimulus in probe D (upper right column) gradually increase in amplitude as the probe passes from the deep white matter (contact D28) into the fundus of the anterior occipital sulcus (aOs) (contact D35). Large potentials with multiple inversions are seen in contacts D38–D52 as the probe passes first through the anterior bank of the sulcus (i.e. the most posterior part of the middle temporal gyrus or mTg) and then the posterior bank (i.e. the most anterior part of the middle occipital gyrus or mOg). The earlier peaks at ~105 and 180 ms invert more medially ({circ}) than do later peaks at ~320, 500 and 680 ms ({bigtriangleup}). The internal contacts of probe O (left column) are in the right superior parietal lobule (sPl) at its junction with the superior occipital gyrus and the external contacts are in the inferior parietal lobule (iPl) at its junction with the middle occipital gyrus. Multiple large peaks are seen to polarity invert between contacts O25, O29 and O32, for example at ~200 ms ({bullet}). These contacts lie in the fundus of the intraparietal sulcus (iPs) at its limit with the superior occipital sulcus. The rapid oscillations in contacts D35–42 and O25–29 are superimposed on a slower oscillation with a period approximately equal to the ISI (, m). Localization of the responsive contacts in visual motion areas was confirmed by the subjective illusions of head movement evoked by electrical stimulation between contacts D38 and D42 and between contacts O22 and O25 (1 s and 2.5 mA). In the lower traces, potentials in the same latency range are shown from probes passing through the parietal, cingulate and Rolandic cortices. The deepest contacts (P4–P11) in the cingulate sulcus (Cis) record a small but polarity-inverting N200 ({diamond}). The probe then passes though the white matter of the superior parietal lobule (contact P28) where it records slowly changing N200–P275–N360–P440 components ({bigtriangleup}, H), which change morphology and then polarity invert as the probe enters the fundus of the most anterior part of the intraparietal sulcus (contacts P32–P39). The recordings again stabilize as the electrode passes through the grey matter of the inferior parietal lobule (contacts P42–P46), with a final low-amplitude inversion at the most lateral contact (P53) in the external grey matter of the superior limit of the inferior parietal lobule ({square}). Polarity inversions of what appear to be the same components ({diamond}) are also recorded by contacts S53–S60 in the supramarginal gyrus (sMg) (the anterior bank of the ascending ramus of the lateral fissure; see the coronal MRI in Fig. 5Go) and similar potentials are recorded in the central sulcus (Ces) (G55; see Fig. 8Go for MRI), as well as by other contacts in the cingulate sulcus (S18 and G13). The three waveforms in each set represent the activation evoked by Mental Counters' variants differing in the material memorized (numbers, letters or figures). No significant material-specific differences are visible. x/y/z coordinates (Talairach and Tournoux, 1988Go) of the probes (actual followed by normalized): probe D (28–45/–71/14 and 27–44/–75/15), probe O (15–36/–79/30 and 15–35/–84/32), probe P (4–53/–39/54 and 4–52/–41/58), probe S (18–60/–29/35 and 18–59/–31/37) and probe G (55/–10/39 and 57/–11/42). See Figure 5Go for electrode localization. Pt. 2.

 


View larger version (55K):
[in this window]
[in a new window]
 
Figure 4. Effect of stimulus modality and required calculation. Waveforms from the most active contacts in Figure 3Go are shown in response to additional task conditions. In the left column no significant differences are seen in the responses recorded by any of these contacts to visual stimuli that required adding versus subtracting from the memoranda. However, there may be subtle differences in the late components to catch trials (e.g. S53) ({square}). In the right column the responses are summed across operations and compared with the corresponding potentials evoked by auditory input stimuli. As expected, no significant response to auditory stimuli is recorded in the visual association cortices (aOs and p-iPs). In contrast, auditory stimuli evoke responses in parietal and central contacts (sMg, iPs, sPs and CeS), which are similar to the responses evoked by visual stimuli except that the auditory responses are 60–100 ms shorter in latency. Most striking in these sites is a negative component at 100 ms to auditory stimuli ({bullet}) that seems to correspond in distribution to the visual N200 ({bigtriangleup}). See Figure 5Go for the electrode localizations and Figure 3Go for the abbreviations and Talairach coordinates (Talairach and Tournoux, 1988Go). Pt. 2.

 


View larger version (130K):
[in this window]
[in a new window]
 
Figure 5. Single-trial spectral analysis of occipito-parieto-Rolandic interactions. Interactions were calculated between seven sites in occipital, parietal, central and medial cortex, the average ERPs of which are shown in Figures 3 and 4GoGo. Each coloured box plots z-scores comparing spectral measures for each frequency (y-axis), at each latency (x-axis), for every trial, to those calculated in the baseline period. Spectral power is plotted in the boxes on the diagonal, phase locking in the boxes at the upper right and phase lag in the boxes at the lower left. In the phase lag plots red indicates that the site listed above the plot leads the site listed to the side. Event-related spectral changes were found in all frequencies and latencies. Across sites the most consistent changes were in the theta and alpha bands (displayed in the bottom of each square). These can be considered in three epochs, which are indicated by vertical tan, grey and pink stripes. Shifts in the phase relations between sites in the alpha and theta bands are visible as alternating red and blue vertical stripes in the lower parts of the time–frequency plots. The first epoch, from ~150 to 280 ms, is characterized by a phase lead of posterior over anterior sites. For example, the ‘early’ occipital leads D38 and O29 (the sites with the earliest polarity-inverting ERPs; see Fig. 3Go) consistently lead the later occipital site (D52), as well as parietal (S53 and P39), cingulate (S18) and central (G55) sites ({bigtriangleup}) (left two columns). The parietal sites in turn lead the central site ({bigtriangleup}). Phase locking decreases at early latencies between several occipital, parietal, medial and central sites (downward open arrows) and then increases (upward open arrows). All sites except O29 showed increased spectral power in the theta range (jagged open circle), with extensions into the gamma band in the early occipital site D38 ({diamond}) and into the alpha band in anterior sites ({diamond}). The second epoch, from ~300 to 400 ms, is characterized by a reversal of phase lags, so that anterior cortex now tends to lead posterior cortex. Specifically, parietal cortex leads occipital and central cortex ({circ}) and central cortex leads occipital cortex ({bullet}). Phase locking may show a second dip, at the transition between different phase lag directions (downward filled arrows). Spectral power tends to decrease, particularly in the occipital sites in the alpha range. The third epoch, from ~40 to –530 ms, is characterized by a restoration of the original phase lag directions, with early occipital cortex again leading central and parietal cortex ({square}). Phase locking is generally high (upward filled arrows). Gamma range increases are present in the central sulcus (H). Data are combined from all input stimuli (I1–I6) of three task variants differing in material (number, letter or figure). Electrode locations are indicated below, on the schema traced from intraoperative teleradiography integrated with MRI and the surface-rendered MRI, as well as sagittal and coronal MRI sections. See Figure 3Go for the Talairach coordinates (Talairach and Tournoux, 1988Go). Pt. 2.

 


View larger version (108K):
[in this window]
[in a new window]
 
Figure 6. Sustained activations and interactions in the peri-Rolandic cortex. Pt. 9 spectra (top). Event-related spectral time–frequency plots were calculated on individual trials and then summed and normalized with respect to the baseline period, for sites in the primary motor (R), premotor (M) and prefrontal (K and F) cortices. Spectral power (left) is plotted continuously over the entire trial (upper left), consisting of six input stimuli (I1–I6), the answer probe (P) and response/feedback (F), as well as collapsed across the input stimuli (I1–I6 at top middle). The most consistent spectral changes are in the theta and low alpha bands from ~300 to 400 ms (vertical grey stripes). These changes include increases in spectral power ({square}) and phase locking ({bigtriangleup} at upper right). Activity in the prefrontal cortex (K49) tends to lead other sites during this epoch (upward open arrow). In addition, a high frequency response (45–50 Hz) is present in K49 ({circ}). During the preceding epoch (~150–280 ms) (vertical tan stripes), spectral activation tends to be in the high alpha and beta range ({diamond}) and other sites tend to lead the premotor cortex (M35) (downward open arrows). Pt. 9 waveforms (middle left). ERPs from ~130 to 600 ms are recorded from two probes (K and M) passing through the precentral sulcus (preCs). Multiple inversions are seen between adjacent contacts separated by 1.5 mm: an N160/N210 complex in K46 inverts to a P170/P220 in K49 ({circ}-arrows) and an N230 in M28 inverts to a P230 in M31 ({bullet}). Additional small non-inverting negativities in the 400–500 ms range are also seen. Note that, as in the postcentral cortex, these responses are small (~25 µV), but the inversions are very clear. Similar responses are seen to stimuli that evoke addition, subtraction or no operation, and that lead to updating numerical versus figural memoranda (N2I versus F2I). Electrical stimulation (2–5 mA, 20 pulses/s and 2.5 s) of M28–M31 and M31–M35 provoked a discrete fall of the left arm and flexion of digits in the left hand. Together with the anatomical location, these results indicate that M28–M35 are located in the premotor cortex. Stimulation (1 pulse/s) of R28–31 evoked contractions of the left hand and fingers, and of R31–34 evoked contractions of the eyelids, confirming localization in the primary motor cortex. Stimulation of K42–K49 and of F23–F30 provoked no response. Probe coordinates (actual followed by normalized): probe M (28–35/3/62 and 28–35/3/70), probe K (46–49/12/36 and 46–49/12/40), probe F (26/32/31 and 26/32/35) and probe R (31/–7/41 and 31/–7/46). Pt. 13 (lower left). A series of potentials in the same latency range polarity inverting over short distances, this time in the postcentral sulcus (posCs). Components at ~210 and ~265 ms invert between S33 and S36 ({bigtriangleup}), again between S40 and S43 ({square}) and possibly again between S47 and S50 (). An earlier potential (peaking at ~180 ms) as well as later potentials (at ~360 and 500 ms), which do not invert, are also seen. Stimulation of S22–S26 (21–24) (1.5 mA) provoked a feeling of a painful electrical discharge in the left hand, stimulation of S33–S36 (1.5 mA) provoked shivers in both palms and stimulation of S40–S43 (1.5 mA) provoked no subjective response. Interictal EEG recorded by probe S was within normal limits and the patient's seizures were found to arise in the temporal lobe, without spread to electrode S. Probe coordinates (actual followed by normalized): probe S (33–50/–23/35 and 30–46/–25/34). Pt. 14 (right). Recordings in the left postcentral sulcus show an extended activation, with peaks at ~200 and 450 ms. The latter peak is large (60 µV) and inverts polarity between P'52 and P'59 (H-arrows). This response is larger to stimuli that direct the updating of numerical ({diamond}) (N2I) as compared to figural ({bullet}) (F2I) memoranda. Similar, but smaller responses are recorded to response probe stimuli (h-arrows) (N2P). Feedback stimuli (h-arrows) (N2F) also evoke similar waveforms, but at ~70–100 ms shorter latencies. Probe coordinates (actual followed by normalized): probe P' (–49 to 59/–14/27 and –47 to 57/–14/26).

 


View larger version (31K):
[in this window]
[in a new window]
 
Figure 7. Widespread embedded fronto-parieto-temporal oscillations. Pt. 14 (left). The potentials recorded across an entire trial consisting of warning (W), six input stimuli (I1–I6), response probe (P) and feedback (F). Strong oscillations begin with the warning stimulus in the left inferior postcentral sulcus (P'52) and continue through the input stimuli (h, {bigtriangleup}). Weaker oscillatory responses are also present in Broca's area (posterior inferior frontal gyrus or piFg) (R'60) where the response is more prominent to the response probe than to the warning. The oscillatory response is weaker in the probe anterior to Broca's area (I'42). Both the prefrontal and parietal responses are larger to stimuli signalling change to numerical memoranda (N2) than identical stimuli signalling updates of figural memoranda (F2). Stimulation (1.5 mA) of R'57–R'60 and adjacent contacts caused an immediate arrest of reading and/or paraphasias. Stimulation of I'42–I'45 (2.5 mA) only evoked paraphasia when it also provoked a seizure that spread to R'60. Thus, R'60 lies in Broca's area according to stimulation and I'42 probably lies immediately anterior to it. Note that, since the time from the probe to the feedback was not constant, ‘FB’ is written at the average time. Probe coordinates (actual followed by normalized): probe R' (–60/11/15 and –57/11/14) and probe I' (–42/26/7 and –39/26/7). Pt. 10 (right upper). Slow oscillations recorded below the left superior temporal sulcus (sTs) (C'56). Probe coordinates (actual followed by normalized): probe C' (–56/–27/–8.5 and –56/–26/–9). Pt. 13 (right lower). Slow oscillations recorded below the right superior temporal sulcus (D61) and, in particular, in the supramarginal gyrus (sMg) (P52) (H, {bigtriangleup}). Like the postcentral sulcus response, the supramarginal gyrus response clearly begins to the warning stimulus. Probe coordinates (actual followed by normalized): probe D (61/–62/14 and 61/–68/14) and probe P (52/–40/31 and 48/–44/31).

 


View larger version (60K):
[in this window]
[in a new window]
 
Figure 8. Widespread CNV-like potentials and embedded oscillations. Pt. 2 (left). Oscillating activations at different time-scales are seen in different simultaneously recorded contacts in occipital, parietal and central sites. Phasic responses lasting 50–100 ms are seen in the supramarginal gyrus ({diamond}; S53), anterior intraparietal sulcus ({bigtriangleup}; P39–P42), posterior intraparietal sulcus ({bullet}; O25–O29) and central sulcus (G41–G48). These responses often habituate to successive input stimuli ({diamond}, {bigtriangleup}, {circ}). For example, in contact S53 the N210–P280 declines from ~85 µV to I1 ({diamond}), ~65 µV to I2 ({diamond}), ~55 µV to I3 and ~45 µV to I4–I6. Similarly, the P280 in P39 declines from ~65 µV to I1 ({bigtriangleup}), to ~45 µV to I2 ({bigtriangleup}) and ~25 µV to subsequent stimuli. A later component at ~500 ms also declines from ~40 to ~20 µV over the same period. Finally, in G48 the N210–P290 peak to peak measure strongly declines in amplitude from 50 µV to I1 to 10–20 µV to succeeding stimuli. A weaker habituation of ~40% occurs in contacts G41 and G44. A dual-peaked response with an overall periodicity equal to the input stimulus ISI of 960 ms inverts between O25 and O29 in the posterior intraparietal sulcus. Polarity inversion is particularly clear for the earlier peak at ~220 ms (h). Habituation across successive stimuli in a trial is prominent in the later peak at ~350 ms, with baseline to peak values declining from ~110 µV to I1 ({bullet}) to ~65 µV to I3–I6 ({circ}). Finally, although sustained responses lasting ~6 s are seen in several leads (e.g. P39), it is largest and polarity inverts between the grey matter of the precentral gyrus (G44) and the central sulcus (G48) (H), according to the MRI (lower right). Direct electrical stimulation between contacts G41–G44 evoked left mouth and tongue twitches. See Figure 5Go for additional electrode locations. Pt. 7 (right upper). Long-duration potentials in another patient are recorded in the left anterior cingulate sulcus (R'21), where the potential is negative and in the vicinity of the H-shaped orbital sulcus (O'15), where it is positive. Note the rapid habituation of the P290 in the cingulate sulcus ({square}, {square}). Probe coordinates (actual followed by normalized): probe R' (–21/36/–20 and –20/41/–20) and probe O' (–15/–10/47 and –15/–11/46).

 


View larger version (58K):
[in this window]
[in a new window]
 
Figure 9. Late occipital, parietal and hippocampal potentials to catch and probe stimuli. Small, late (peak latency ~500 ms), negative potentials are recorded to ‘catch’ input stimuli (N2I) in the supramarginal gyrus (P47) ({bullet}) and the most posterior portion of the middle temporal gyrus (p-mTg) (O41 and O44) ({square}), anterior to the usual location of area MT+. The more medial contacts on the O probe (which electrical stimulation indicated were located in the retinotopic cortex) do not show this response. Potentials with the same latencies and polarity are recorded in the same sites as false probes (N2P) ({bigtriangleup}). The same potentials appear to be evoked at an earlier latency to feedback indicating incorrect responses (N2F) ({diamond}). An earlier latency is expected because the feedback signals are in the auditory modality. Simultaneously recorded hippocampal formation (HCF) contacts (B40) are selectively responsive to the feedback tones (h). Scalp potentials at Pz do not closely resemble any of the depth sites and are usually inverted in polarity from the posterior portion of the middle temporal gyrus and supramarginal gyrus recordings. Probe coordinates (actual followed by normalized): probe P (47/–40/30 and 48/–28/29), probe O (34–44/–59/20 and 34–45/–56/19) and probe B (40/–18/–16 and 41/–17/–16). Pt. 5.

 


View larger version (54K):
[in this window]
[in a new window]
 
Figure 10. Late potentials in anterior cingulate gyrus and ventral temporofrontal sites to feedback and catch stimuli. Pt. 3 (upper left). A negative potential at ~400 ms to catch trials in the posterior inferior frontal gyrus (p-iFg) (G52) ({bullet}) is not recorded in the adjacent contact, suggesting local generation. Probe coordinates (actual followed by normalized): probe G (52/27/1 and 54/25/1). Pt. 6 (middle left). Negative potentials at ~300 ms are recorded in the anterior cingulate gyrus (aCg) (O3) ({square}) to feedback tones indicating wrong responses (L2F). Probe coordinates (actual followed by normalized): probe O (3/38/–5 and 3/38/–5). Pt. 10 (lower left). Negative potentials at ~300 ms are again recorded in the anterior cingulate gyrus (O'7) ({bigtriangleup}) to feedback indicating wrong responses (N2F). Probe coordinates (actual followed by normalized): O' (7/45/–13 and –7/44/–12). Pt. 1 (right). Again, a potential to incorrect feedback tones is recorded in the anterior cingulate gyrus (G12) ({diamond}), as well as in other limbic sites including the posterior hippocampus (p-HC) (C26) ({square}) and the posterior parahippocampal gyrus (p-pHCg) (D24) ({circ}). These sites all generate a positivity at ~300 ms that resembles the simultaneous P300 recorded at the scalp (Pz) ({bigtriangleup}). The same sites showed little or no response to the input stimuli, except for the catch trials (shown only for the anterior cingulate gyrus) (G9) ({diamond}). Probe coordinates (actual followed by normalized): probe D (24/–45/11 and 24/–44/11), probe C (26/26/–39/–3 and 26/–38/–4) and probe G (9–12/36/2 and 9–12/36/2).

 
Visual Association Cortices

Multiple inverting evoked potentials were recorded in the dorsal visual pathway, with peaks as early as ~100 ms and continuing until ~700 ms, i.e. virtually the end of the recording epoch. Local potential inversions and gradients demonstrated that these potentials were locally generated. The precise anatomical location of the most responsive cortical recording sites indicated that they were in the visual association cortex of the dorsal stream, an inference confirmed in some cases by the effects of local electrical stimulation.

Probes sampling two such sites are shown at the top of Figure 3Go. Multiple large peaks separated by ~100 ms inverted repeatedly between successive contacts of probes D and O. The different components inverted in different locations, indicating that the earlier components were generated in more medial cortex. These rapid oscillations were superimposed on a slower oscillation with a period approximately equal to the ISI, i.e. 840 ms. The slower oscillations also showed multiple inversions, large amplitudes and steep voltage gradients, and thus were locally generated (Figs 3 and 8GoGo). Leads passing through more anterior association cortex also recorded similar responses. Such leads could be located anterior to area MT+ (probe O, Fig. 9Go) or just below the posterior superior temporal sulcus (contacts C'56 and D61, Fig. 7Go).

Local electrical stimulation of probe D in patient 2 evoked illusory head movements. The location of probe D (visualized directly on his MRI, Fig. 5Go) corresponded to that of human visual motion area MT+. Earlier studies have found specific responses in this area to coherent as compared to incoherent movement in humans (Watson et al., 1993Go; Tootell et al., 1995Go; Ulbert et al., 2001Go) and macaques (Lappe and Duffy, 1999Go; Andersen et al., 2000Go), leading to theories that assign a role to area MT+ (and in particular to MSTd) in visual flow phenomena, inferring the direction and speed of head motion from coherent visual motion. The current finding that stimulation in this area resulted in a sensation of head motion is consistent with these theories. This contrasts with the experience of moving visual stimuli evoked by stimulation of another visual motion area in the medial parieto-occipital sulcus, which appears to correspond to SPO or V7 (Richer et al., 1991Go). Illusory head movements were also evoked by local electrical stimulation of probe O, which appeared to be located in an area that has also been associated with visual motion processing with fMRI, and termed LIP (Brandt et al., 1999Go). The neural correlates of activation in this structure and its primate homologue are poorly understood.

Occipital responses typically did not change with different stimulus materials (numbers, letters or figures; Fig. 3Go) or different processing requirements (adding, subtracting or no operation; Fig. 4Go). However, responses to auditory stimuli in the homologous task were absent (Fig. 4Go). The responses to the probe stimuli were similar to those evoked by the I1–I6 stimuli (Fig. 7Go). There were exceptions, for example the most external contacts in probe O of Figure 9Go, which showed an enhanced late response to the catch stimuli, peaking at ~500 ms. This response was seen in the same leads to the feedback tones (in particular to those that indicated an incorrect response) and to the probe stimuli (in particular those that indicated a mismatch to the numbers held in working memory). In some cases, habituation was observed across successive input stimuli in a trial (e.g. contact O29, Fig. 8Go).

Parietal Lobe

A series of potentials from ~130 to 500 ms were recorded in multiple parietal sites, including the intraparietal sulcus, inferior and superior parietal lobules, posterior cingulate and supramarginal and postcentral regions (Figs 3–9GoGoGoGoGoGoGo). Typically, parietal waveforms exhibited negative–positive–negative–positive peaks at ~200–280–350–430 ms. Frequent polarity inversions, most prominently in lateral inferior parietal sites, indicated local generation, for example components in Figure 3Go inverted at approximately the same latency in the most anterior part of the intraparietal sulcus (a-iPs) and the supramarginal gyrus (sMg) (see also probe P, Fig. 9Go). Similar potentials were also recorded in several locations near the cingulate sulcus (Cis; probes P, S and G of Fig. 3Go) and in the postcentral sulcus (posCs; probe S of Fig. 6Go), where they again polarity inverted over short distances. Although relatively low amplitude, the double inversion of peaks at ~210–265 ms between successive contacts of the probe confirmed local generation. Stimulation of the probe evoked somatosensory responses that were typical of non-primary cortex (Penfield and Jasper, 1954Go).

Later potentials were also found to invert in the postcentral sulcus, for example at 450 ms (probe P', Fig. 6Go). When viewed over a longer time-scale, this late potential was seen to form a slow oscillation with approximately the same periodicity as the ISI (Fig. 7Go). Superimposed on these slow oscillations are the rapid oscillations reflected in the early peaks. Even longer duration modulatory responses (~6 s) were rare in the parietal lobe recordings, but an example of a potential slowly developing over the course of I1–I6 stimuli was seen at the anterior limit of the intraparietal sulcus (Fig. 8Go).

Responses in the supramarginal gyrus from ~200 to 500 ms commonly habituated over the course of I1–I6 stimuli. The parietal responses usually did not significantly differ to I1–I6 stimuli signalling addition versus subtraction versus no action, nor in relation to the material of the stimuli. However, there were rare exceptions. For example, the longest latency component to no-operation trials was somewhat larger than to the subtraction or addition trials in the supramarginal contacts shown in Figures 4 and 9GoGo. Similarly, the response to numbers was clearly larger than that to figures in the left postcentral sulcus (P'52) site shown in Figure 6Go (N2I versus F2I). Oscillatory potentials were also observed to auditory stimuli when the homologous task was administered (Fig. 4Go). These components appeared to be similar to those evoked by the visual stimuli, except for an ~90 ms earlier latency in the auditory modality. This is consistent with other SEEG studies finding both visual and auditory responses in this area (Halgren et al., 1995aGo).

Parietal sites also generally responded to the probe and feedback stimuli. Often it appeared that fewer oscillations might be present to these stimuli, in particular to the feedback stimuli. However, this was difficult to confirm because of signal:noise considerations: parietal responses are often small and require substantial averaging, but six times as many trials were available by combining across input stimuli I1–I6 as compared to the probe or feedback stimuli. For example, in the left inferior post-central sulcus (probe P') (Fig. 6Go), focal and/or polarity-inverting components from ~200 to 500 ms were evoked by the probe stimuli. These resembled those evoked by the input stimuli, but had smaller amplitudes. The feedback stimuli also evoked similar waveforms, except they were smaller and of earlier latency, corresponding to the shorter latencies often observed in the auditory modality. Similarly, the supramarginal site shown in Figure 9Go showed an enhanced late potential to the wrong feedback, catch input trials and probe stimuli.

Frontal Lobe

The parietal pattern of rapid oscillations (~12 Hz) from ~130 to 600 ms sometimes superimposed on slower waves with the period of an entire stimulus or trial was also observed in frontal recordings. For example, the latencies, polarities and waveforms of the components recorded by probe G in the precentral area in Figure 3Go are essentially identical to those recorded by probe S in the supramarginal gyrus of the same patient. Direct visualization of the probe's path with MRI (Fig. 8Go), as well as the effects of direct electrical stimulation, confirmed that this probe was located in the primary motor cortex. Although no inversions were seen in this probe, other nearby probes did record them. Two such probes (K and M) in the precentral sulcus are shown in Figure 6Go. Note that, as in the postcentral cortex, these responses were small (20–30 µV), but the multiple inversions between contacts separated by 1.5 mm provide clear evidence of local generation. Stimulation-induced movements indicated that the superior probe is located in the premotor cortex (Chauvel et al., 1996bGo). Habituation of frontal responses was common across stimuli I1–I6 (Fig. 8Go).

Cycling activation with periodicity of the input stimuli (~1.2 Hz) was also recorded in and near Broca's area (see, for example, contacts I'42 and R'60 in Fig. 7Go). Note that the onset of each cycle of this activation appeared to precede the stimulus onset. The slower oscillations in the frontal cortex were focal, but frank polarity inversions were not observed.

Long duration responses (~6 s) were observed in several frontal contacts. A contingent negative variation (CNV) that was negative within the precentral gyrus is shown to polarity invert at its surface in Figure 8Go. Electrical stimulation of these contacts identified the generating cortex as pri