Functional Brain Mapping Using Intracranial Source Imagingengold.ui.ac.ir/~mahnam/IBME/Research/Functional Brain Mapping.pdf · Functional Brain Mapping Using Intracranial Source

1140_book.fm Page 1 Wednesday, July 9, 2003 8:46 AM

12Functional Brain

Mapping UsingIntracranial Source

Imaging

CONTENTS12.1 Introduction12.2 Neurophysiological Signals

Origin of Activity • Localization of the Sources

12.3 Mathematical ModelingSource Models • Volume Conductor Models • Source Estimation Procedure

12.4 Neurophysiological RecordingsElectroencephalography — EEG • Magnetoencephalography — MEG

12.5 Data AnalysisSurface Maps • Source Imaging Procedure • Mapping-Based on EEG or MEG

12.6 Application Examples Auditory Evoked Responses • Somatosensory Evoked Responses • Visual Evoked Responses • Epileptogenic Spike Localization • Language Mapping

12.7 Concluding RemarksReferences

12.1 Introduction

The treatment of certain complex neurological diseases, such as pharmacologically intractable epilepsy,brain tumors, and arteriovenous malformations, often includes surgical intervention. In all cases, accuratelocalization of the lesion to be resected is of paramount importance, and the margin between therapeutictreatment and debilitating side effects due to surgery is very narrow. Several additional factors contributeto the success of surgery, chief among which is the accurate identification of the cortical areas that areresponsible for certain brain functions, such as sensation, movement, and speech. This latter procedure,known as functional brain mapping, is critical, because resection of such vital brain areas can havedevastating results.

George ZouridakisUniversity of Houston

Darshan IyerUniversity of Houston

© 2004 by CRC Press LLC


A number of noninvasive functional imaging modalities, including functional magnetic resonanceimaging (fMRI) (Binder et al., 1997; Cuenod et al., 1995), positron emission tomography (PET) (O’Learyet al., 1996; Peterson et al., 1988), regional cerebral blood flow (rCBF) (Friberg, 1993), and single photonemission computed tomography (SPECT) (Gomez-Tortosa et al., 1994), can map brain functions withvarying degrees of success. However, the most reliable approach to map the human brain still relies ondirect electrical stimulation of the exposed cortex (Ojemann et al., 1989), a procedure that is highlyinvasive and unpleasant, and it is performed mostly intraoperatively on an awake patient.

Recently, however, completely noninvasive procedures have been successfully used for brain mapping(Breier et al., 2000; Ebersole and Wade, 1990; Gallen et al., 1994; Hamalainen et al., 1993; Papanicolaou etal., 1999; Zouridakis et al., 1998b). These procedures rely on the fact that performance of certain brainfunctions activates only a small population of cortical neurons whose activation gives rise to electromagneticsignals that can be recorded externally. The electrical aspects of brain activation are recorded in an electro-encephalogram (EEG) by placing a set of electrodes on the scalp, while the corresponding magnetic aspectscan be captured in a magnetoencephalogram (MEG) by placing an array of coils close to the head.

Electrical or magnetic source imaging (SI) refers to the localization of the intracranial sources thatgive rise to externally recorded electrical or magnetic signals, respectively. This procedure is a form offunctional imaging and combines the neurophysiological EEG or MEG data with structural MRI scans.That is, while MRI alone provides information on the brain’s anatomy, SI provides information aboutits function. Moreover, SI provides very high temporal and spatial resolution, which makes it far superiorthan any other functional imaging modality in that respect.

In this chapter we introduce the basic principles that make EEG- and MEG-based SI procedurespossible, and the equipment used to accomplish it. Additionally, we provide a few concrete examples thatdemonstrate the usefulness of SI as a valuable clinical tool for functional brain mapping.

12.2 Neurophysiological Signals

12.2.1 Origin of Activity

When considering surface neurophysiological activity, two major aspects must be distinguished. Thefirst one is the characterization of the sources underlying the observed signals. The human brain iscomposed of vast numbers of electrically active neurons and other supporting cells (e.g., glial cells)that are assembled in functional groups. In particular, the outer surface of the brain, the cerebralcortex, consists of a thin and highly compact intricate network of cells arranged in layers (Kandel etal., 1991). The average thickness of the cortex is 2.5 mm, and its density reaches an impressive 105

cells per millimeter which form approximately 1015 synapses. Conceptually, the building block of thecerebral cortex is a column, a structure of neurons arranged radially to the head surface, with a diameterof about 0.5 to 3 mm (Kelly, 1991). A specific type of neurons in a column, called pyramidal cells,have a linear structure with dendrites that are arranged parallel to each other. In the resting state, eachof these neurons is negatively polarized relative to the surrounding electrolyte solution, because of theactivity of the cell membrane. When a cell is in the active state, large quantities of positive and negativeions — namely, sodium (Na+), potassium (K+), and chloride (Cl–) — cross the cell membrane, movingfrom the intracellular to the extracellular fluid, and vice versa. For all practical purposes, this ionmovement is equivalent to a current flow, and it is responsible for all electrophysiological signalsrecorded externally. In particular, the electrical potentials and the magnetic fields recorded on thescalp as the EEG and MEG, respectively, result mainly from the temporal and spatial summation ofpostsynaptic activity generated along the apical dendrites of the pyramidal neurons. These fieldsrepresent the synchronous activation of a large number of neurons, estimated to be between 104 and105 cells (Okada et al., 1992; Williamson and Kaufman, 1987).

The second aspect of surface neurophysiological activity is the influence of the tissues surroundingthe active neurons, the so-called volume conductor effects. A detailed analysis of the relationship betweenthe intracranial sources and the extracranially recorded signals is beyond the scope of this chapter, but



excellent reviews can be found elsewhere (Hamalainen et al., 1993; Williamson and Kaufman, 1987).Briefly, neural activity can be represented mathematically as a current density in a closed volume of finiteconductivity. Outside this volume, the conductivity and current density are zero. The intracellularcurrents give rise to an electric field outside the cells which, in turn, results in currents that flow passivelythrough the rest of the conducting medium. Thus, the total current density can be divided into twocomponents, primary (intracellular) and secondary (extracellular). The primary currents are more inter-esting in neurophysiology, because they are associated directly with cell activation. Under specific con-ditions (Lewine, 1990), that are fairly well approximated in many practical applications, it can be assumedthat the electromagnetic fields recorded on the scalp as the familiar EEG and MEG are associated mainlywith the primary (intracellular) currents, and thus they represent cell activation.

The relationship between the intracranial sources and the extracranially recorded biological signalscan be computed as a function of the distance between the observation point and the source location.Both the magnetic field and the potential are inversely proportional to the square of this distance.

12.2.2 Localization of the Sources

Determining the extracranial potentials and magnetic fields that result from the intracranial currentsources is known as the forward problem. Conversely, the inverse problem requires localization of theintracranial current sources that give rise to the externally recorded electrical potentials (EEG) andmagnetic fields (MEG). Unfortunately, the inverse problem has no unique solution, as the correspondingmathematical equations can be satisfied by an infinite number of possible source configurations. Toovercome this limitation, it is necessary to solve a “constrained inverse problem” and make assumptionsfor both the sources and the surrounding head tissues (Morse and Feshbach, 1963). These assumptionsmust be neurophysiologically based, otherwise the results computed and their interpretation may becompletely erroneous (Jayakar et al, 1991; Nunez et al., 1991). A typical approach involves the so-called“quasi-static” approximation (Plonsey and Heppner, 1967) under which all potentials and currents atany given instant in time are determined by the properties of the sources at that time only. Moreover,the various brain tissues are assumed to be linear and completely characterized by their electrical con-ductivity, which is frequency independent. The medium can be inhomogeneous, if its conductivity isdifferent for different compartments, and anisotropic, if the conductivity of the tissues is different alongdifferent directions in the three-dimensional space.

12.3 Mathematical Modeling

12.3.1 Source Models

Over the past several years, many different approaches have been developed to identify the sources thatexplain the fields recorded on the head surface. The most common approach has been to describe theactivation of a well-localized (Nunez et al., 1991) small population of neurons at a particular point intime as a single equivalent current dipole. Indeed, when a neuron is activated, ion currents flow throughit and this “electrical generator” can be modeled as a small current dipole, with a moment of about 10–13

A m. A schematic diagram of such an arrangement is shown in Figure 12.1.However, to identify the sources that account for a complex field distribution over a time interval, a

series of single dipoles is required (Supek and Aine, 1993). These dipoles can move in position andorientation over time, and when displayed on a single image, they appear as a moving dipole.15,16 Anotherapproach is to consider current distributions in the brain (Ioannides et al., 1990; Morse and Feshbach,1963; Okada et al., 1992; Wang et al., 1992), and estimate several fixed dipoles at discrete points withina region of interest. Spatiotemporal source modeling (Baumgartner et al., 1991; Scherg, 1992; Schergand von Cramon, 1985) is yet another approach that takes into consideration overlapping activity ofmultiple generators. In this case, dipoles are fixed in location and orientation, but they can vary overtime in strength and polarity to explain the temporal evolution of the recorded fields.



In general, solutions for multiple-dipole sources are much less reliable than those for a single dipole(Cuffin, 1998), because in the former case the solutions provided are very sensitive to noise. Thus, mostclinical applications still use single-dipole approaches.

12.3.2 Volume Conductor Models

Volume conductor models have evolved from extremely simple (Rush and Driscoll, 1978) to very realisticand computationally intensive ones (Fender, 1991). The first and simplest model was that of an infinite,homogeneous, isotropic medium. The next sophistication level considered a homogeneous but boundedmedium (specifically, a sphere) as a head model (Wilson and Bayley, 1950). Later on, to account for thevarious tissue layers within the head, several concentric-sphere models were developed (Ruth and Driscoll,1968; Schneider, 1972; Witner et al., 1972) that included three concentric spheres to account for the scalp,skull, and the brain. More recently, however, additional approaches, such as the eccentric-sphere model(Mejis and Peters, 1987) and the four-shell model, have been proposed. The latter uses a cerebrospinalfluid layer underneath the skull, and it is especially useful for hydrocephalic patients (Nishijo et al., 1996).

However, more realistic head models can be obtained if, instead of using arbitrary spheres, MRIs areemployed to determine the curvature of a particular subject’s head and fit the spheres of the multicom-partment model (Lopes da Silva et al., 1991). More realistic models computed entirely from MRIs (Mejisand Peters, 1987) are implemented using numerical techniques, such as the boundary element method(BEM) (Mejis et al., 1989) and the finite element method (FEM). Finally, realistic head models thataccount also for the effects of anisotropy in the electrical conductivity have been developed (de Munck,1988; Zhou and van Oosterom, 1991).

12.3.3 Source Estimation Procedure

Initially, inverse source localization procedures were implemented on spherical models (Rush andDriscoll, 1978; Smith et al., 1993) and they were later extended to work on realistic geometry. In thecase of a single-current dipole, to identify the intracranial sources that explain the fields recorded onthe head surface, the procedure takes the measured values from all sensors at a given instant in timeand searches, using iterative minimization techniques, for an equivalent dipole within the head thatcould generate such fields (Henderson et al., 1975; Kavanagh et al., 1978; Schneider, 1972; Sidman et

FIGURE 12.1 A well-localized small population of neurons, with pyramidal cells oriented radially to the headsurface. The sum of currents (arrows on the cells) flowing in the apical dendrites can be modeled as a single equivalentcurrent dipole.



al., 1978). Techniques such as the brain electrical source analysis (BESA) (Jirsa et al., 2000) and CURRY(Neuroscan, El Paso, TX) are iterative with automatic picking of dipole sources and define one dipolefor each point in time.

In more general cases, the localization procedure starts with some initial estimate of the distributedsources and then recursively enhances the strength of some of the elements, while decreasing the strengthof the rest of the elements until they become zero. In the end, only a small number of elements willremain nonzero, yielding a localized solution. This method is implemented, for example, in the FOCUSS(Gorodnitsky et al., 1995) and LORETA (Pascual-Marqui et al., 1994) algorithms.

Another approach to source localization employs a spatiotemporal model, under the assumption thatthere are several dipolar sources that maintain their position and orientation fixed, while they vary onlytheir amplitude as a function of time. Rather than fitting dipoles to measurements from one instant intime, dipoles are fitted by minimizing the residual least-square error over the entire time interval (Schergand von Cramon, 1985). More advanced spatiotemporal approaches, such as the multiple signal classi-fication algorithm (MUSIC) (Mosher et al., 1992), have also been developed, whereby the algorithmsearches for a single dipole through the three-dimensional head volume. To localize the source, the usermust search the head volume for local peaks. Extensions of this algorithm automate this search througha recursive approach (Mosher and Leahy, 1963).

12.4 Neurophysiological Recordings

12.4.1 Electroencephalography — EEG

The EEG is a record of the electrical aspects of neurophysiological activity. More specifically, a one-channelEEG represents the potential difference between two sites on the head, and requires two electrodes connectedto the input of a differential amplifier — one active and one for reference. These electrodes are typicallygold-plated or sintered Ag-AgCl, and they are attached to the scalp with a special glue (collodion) orelectrolyte gel. Typically, the active electrode is placed close to the structures of interest, whereas the referenceelectrode is placed at a distant location (usually on the left or linked mastoids, behind the ears). This typeof connection is called monopolar, or referential. When both the active and the reference electrodes areplaced close to the structures of interest, the placement is called bipolar or differential. The advantage ofthis latter connection is that far-field activity common to both electrodes is cancelled and thus sharplocalizations of events can be obtained. Alternatively, the reference electrode can be connected to a circuitthat measures the average activity of all electrodes to obtain a recording with an average reference. Thecommon electrical connection for the electrical circuitry of the recording equipment is called the ground,and all voltages are measured with reference to it. The patient is also connected to this point through aseparate electrode typically placed on the patient’s forehead. The potential differences in each channel areamplified by high-gain, differential amplifiers, and then it is either displayed on a monitor, saved on a harddisk, or simply plotted on paper, depending on the technology used.

For clinical recordings, electrodes are arrange on the head in a montage that follows the 10–20 inter-national placement system and its extensions. Routine clinical equipment employs about 21 electrodes,whereas more recent systems specialized primarily for research may include hundreds of electrodes. Inthis case, arrays of electrodes are arranged on a cap for convenience, and these densely sampled recordingsprovide the so-called dense-array EEG, or dEEG. A system that uses 256 channels is described next.

12.4.1.1 EEG Recording Device

The latest dEEG systems incorporate several unique features: they offer high spatial sampling with up to256 recording channels; they are portable, and often the complete set of 256 amplifiers is battery operatedand is about the size of a book; they can use active electrodes for noise cancellation, in which case eachelectrode incorporates a micropreamplifier that preconditions the EEG signals before transmission tothe main amplifiers through the electrode leads, and this eliminates the need for a shielded room; andfinally, they provide very high sampling rates of up to 5 kHz per channel.



The active electrode is a sensor with a very low output impedance. By integrating the first amplifierstage with a regular electrode, extremely low-noise measurements are now possible without any skinpreparation. The noise levels achieved can be as low as the thermal noise of the electrode impedance,which is the theoretically minimum level. Active electrodes eliminate all problems associated with highelectrode impedances, cable shielding, capacitive coupling between the cables, and sources of interference,as well as any artifacts due to cable or connector movements. The amplifiers of a dEEG system (ActiveTwo,BioSemi, the Netherlands) available in our lab, along with an active electrode and a cap with the electrodeholders, are shown in Figure 12.2.

12.4.2 Magnetoencephalography — MEG

The MEG measures the extracranial magnetic fields produced by intracranial electrical currents. Neuro-magnetic signals are many orders of magnitude weaker than the ambient magnetic noise, which is dueto the earth’s field and to the presence of ferromagnetic objects and electrical instrumentation. Typicalscalp-recorded magnetic fields have a peak amplitude of about 100 fT (Gomez-Tortosa et al, 1994; Lewine,1990), whereas environmental electromagnetic noise in a hospital (power lines, elevators, MRI magnets,etc.) may be as high as 1 T (1015 fT) in extreme cases (Lewine, 1990). Therefore, to detect this kind ofbiological activity, it is necessary to use highly sensitive instrumentation and, at the same time, attemptto eliminate extraneous magnetic fields.

12.4.2.1 MEG Recording Device

MEG measurements were practically impossible before the introduction of superconductive instru-mentation. The latest generation of biomagnetometers are composed of large arrays of pick-up coils,each of which is connected to a SQUID (superconducting quantum interference device) that acts asa very low noise, ultrahigh gain, current-to-voltage converter. The SQUIDs and induction coils areimmersed in liquid helium to maintain a superconducting state. This type of device can detect evenvery small changes in magnetic flux, such as the one resulting from neurophysiological activity. Figure12.3 shows a schematic diagram of a multisensor MEG system along with a detection coil and a SQUIDof a single channel.

FIGURE 12.2 Active electrodes, electrode cap, and amplifiers of a portable dEEG system with 256 channels.



To reduce the amount of extraneous magnetic noise, MEG systems are operated in specially designedmagnetically shielded rooms (Gallen et al., 1994). Additional improvements in the quality of the MEGsignals are obtained by selecting an appropriate type for the detection coils. Coils can be single-loopmagnetometers that measure the magnetic field directly, and first-, second-, or third-order gradiometersthat consist of two or more loops that measure spatial changes of magnetic field (Flynn, 1994). Variousexamples of coils for both planar and axial geometry are shown in Figure 12.4.

A gradiometer (Figure 12.4d) can selectively cancel out environmental magnetic noise, depending on its“baseline,” i.e., the distance between the two coils; because the magnetic fields decrease very rapidly withdistance from the source, fields generated by a neural source close to the gradiometer will be detected withdifferent strengths in each loop, and therefore they will induce currents of different magnitude in the twoloops, thus resulting in a nonzero net current in the gradiometer. On the other hand, fields generated faraway from the gradiometer will be practically uniform in space and will result in zero net current.

Because of the sensitivity of MEG recordings to noise, patients implanted with electrically activemedical devices such as cardiac pacemakers, neurostimulators, and infusion pumps cannot be studiedwith MEG, due to the electrical interference from the implanted devices. Also, equipment entering theroom must be screened for large metallic components or electromagnetic activity that could introduceartifactual signals.

The cost and complexity of instrumentation led to the initial development of MEG systems containingonly a few channels. Recently, however, whole-head systems with large arrays comprising 248 magne-tometers have been developed (model 360WH, 4DNeuroimaging, San Diego, CA). An example of a 148-channel system is shown in Figure 12.5. The need for cryogenics, and for a shielded room and a rigidhelmet-type sensor limit the scope of applications of the MEG-based mapping procedures.

FIGURE 12.3 Schematic diagram of a multisensor MEG system (left) along with a detection coil and SQUID in asingle channel (right).

FIGURE 12.4 MEG detection coils include planar (a) magnetometers and (b) gradiometers, and axial (c) magne-tometers and (d) gradiometers. In a gradiometer, a uniform field B creates currents J and –J of opposite polarity aproduces a zero net signal.

Squid

Coil

Single SensorDewar

Liquid Helium

Skull



12.5 Data Analysis

12.5.1 Surface Maps

Since both EEG and MEG activity have a common basis, their analyses can be carried out using similarmethods. During the last 30 years, new techniques, such as color mapping of activity distribution overthe scalp (Desmedt et al., 1987; Duffy et al., 1979), have shifted the focus from an analysis of waveforms,derived at selected scalp sites, to the spatial distribution over the entire scalp, at particular points in time.The recording sensors typically span a two-dimensional surface over the skull, and each sensor measuresa time series. When these measurements are interpolated and color-coded to provide a spatially contin-uous representation of the measurements, a so-called topographic scalp map is obtained at each timepoint. Figure 12.6 shows the location of the sensors on the head and a three-dimensional map of theelectrical fields corresponding to the peak of the N1 auditory component.

Topographic analysis is either done qualitatively, by interpreting focal maxima and minima in themap, or quantitatively by dipole localization methods (Brandt, 1992). Alternatively, a combined spa-tiotemporal approach, where topographic information over different instances in time, called a space-time series, is taken for source localization.

Recently, much effort has been made in the development of source imaging techniques based on high-resolution or dense-array EEG (dEEG). For example, techniques such as scalp Laplacian mapping andcortical imaging attempt to restore the high-frequency content of brain electrical activity that is smearedand distorted by the low-conductivity skull (Babiloni et al., 1996; He et al., 1995; Hjorth, 1975; Le et al.,1994). These techniques can enhance the spatial resolution of the EEG by deconvolving the low-passspatial filtering effect of the head volume conduction (Freeman, 1980; He et al., 1996; Sidman et al.,1992). Since a surface cortical map lies closer to the actual sources, its potential distribution may providemore direct view of the neural sources (Fender, 1991; Jayakar et al., 1991).

FIGURE 12.5 A whole-head MEG system with 148 recording channels operated in a magnetically shielded room.



12.5.2 Source Imaging Procedure

The SI procedure consists of several steps that culminate in the display of functional information obtainedusing MEG or EEG onto high-resolution anatomic images typically obtained by MRI. For this final step,it is necessary to translate EEG/MEG locations into MRI coordinates.

The precise location of the measurement points on the scalp is determined electronically withreference to a Cartesian coordinate system anchored on three landmarks (fiducial points) on eachsubject’s head: two external ear canal points and the nasion. The line passing through the two pre-auricular points defines the y-axis of the system. The line perpendicular to the y-axis passing throughthe nasion defines the x-axis, and the line perpendicular to the x-y plane passing through the x-yorigin defines the z-axis.

The EEG/MEG anatomic reference frame is established electronically: a set of three receivers trian-gulate the signal from a stylus-type transmitter placed successively at several reference points on thesubject’s head (typically the three fiducial points, the vertex, and the inion). Additionally, small vitaminE-containing capsules visible on MRI are placed on the fiducials. The locations of the common markerson the MR images and the EEG/MEG measurements serve for translating EEG/MEG locations intothe MRI reference frame. Validation studies have shown that the accuracy of this approach can be ashigh as a few millimeters (Gallen et al., 1994). The same stylus transmitter can be used to define thecurvature of the head by tracing its surface.

After the anatomic reference frame has been established, the fiducial points registered, and the patient’shead digitized, the actual recordings are obtained. The latter consist of time-varying measurements ateach detector position. The data undergo further analyses that allow the sources underlying the recordedactivity to be localized.

The equivalent current dipole (ECD) model yields a description of the instantaneous current dipolein terms of its location, strength, and orientation, along with an estimate of its reliability (confidencevolume). Even though there is some variation across labs regarding acceptance criteria for the dipolesolutions, a high correlation between data and model and a small confidence volume are always desirable.It is, therefore, possible to select a set of “best-fitting dipoles” to describe the data.

Once the best fitting dipoles have been identified, they are co-registered onto the subject’s MRI scan.The resulting images are then printed on MRI films with different symbols and colors, each correspondingto a distinct type of neurophysiological activity.

FIGURE 12.6 Example of three-dimensional topographic map. (Left) Arrangement of sensors on the head; (right)a three-dimensional map of the electrical fields obtained at the peak of the N1 auditory component.



12.5.3 Mapping-Based on EEG or MEG

Both theory and experiment suggest that the MEG offers no significant advantage over the EEG (Cohenand Cuffin, 1991). MEG systems, however, are very expensive (total cost about $3 million); they requirespecial cryogenic equipment, a magnetically shielded room, daily monitoring and maintenance, and theyare available only in a handful of places around the world. Thus, the clinical usefulness of MEG can bevery limited. On the contrary, in addition to the unique features mentioned earlier, the latest EEG systemsincorporate several advantages: they are readily available in practically all clinical settings, and even themost sophisticated systems are much less expensive than MEG (total cost about $150,000). Therefore,successful brain mapping based on EEG can have a significant impact on patient care.

12.6 Application Examples

In general, the neurophysiological signals recorded on the surface of the head can be separated into twocategories, spontaneous and elicited. For instance, epileptogenic discharges, such as interictal spikes, areevents of the former type. The latter type of activity, known as evoked potentials in the case of EEG andevoked fields in the case of MEG, results from external stimulation of a specific sensory pathway, suchas the auditory, the somatosensory, or the visual.

In the next sections we give examples from our studies during the past few years that illustrate theprocess of measuring surface activity and localizing the underlying sources.

12.6.1 Auditory Evoked Responses

The most prominent component of an evoked response obtained from transient auditory stimulationis the N1, which occurs at approximately 100 msec after stimulus onset. Due to the magnitude andconsistency of the response, the sources underlying the N1 can be easily localized. Typically, neuro-physiological signals are collected in “trials,” each of which consists of a few hundred milliseconds ofactivity prior to and following the stimulus onset. Single-trial responses are contaminated by noise,and to improve the signal-to-noise ratio, many trials are collected and averaged (typically between 100and 500), using the stimulus onset as a time reference. The averaging approach relies on the assumptionthat the neuronal responses to all stimuli are identical. Additionally, to eliminate artifacts from periodicevents, such as arterial pulsations or the 60-Hz power line interference, the interstimulus interval israndomly varied within a predefined range.

A variety of stimuli can elicit this component. We used pure tones of 1-kHz frequency and 50 msecduration (10-msec rise/fall and 30-msec plateau) that were delivered binaurally at an intensity of 80 dBnHL, and a mean stimulus rate of 0.4/sec (Iyer et al., 2002). Figure 12.7 shows an example of the resultingevoked response. Each of the superimposed tracings represents the activity recorded on one of the 256channels, while the highlighted area denotes the peak of the N1 component. The three-dimensionalsurface maps superimposed on the subject’s MRI were reconstructed at the N1 peak, while the two dipolesindicate the cortical areas activated at that time.

The estimated sources underlying the N1 component were consistently found to be on the floorof the Sylvian fissure, i.e., in the area of the primary auditory cortex, as shown in Figure 12.8. Thesources depicted correspond to activity during the highlighted interval around the N1 peak, asshown in Figure 12.7.

12.6.2 Somatosensory Evoked Responses

Localization of the central sulcus and the adjacent precentral and postcentral gyri is a fundamentalobjective in most mapping studies involving patients, when brain lesions are located in the parietal orfrontal cortex. In a recent MEG study (Zouridakis et al., 1999), we used tactile stimulation to map thesomatosensory cortex. Stimuli, i.e., bursts of compressed air, were delivered to a plastic diaphragm clippedto the patient’s fingertip, lip, or toe (Benzel et al., 1993). Approximately 500 single trials were required,



each composed of 100 msec of prestimulus and 100 msec of poststimulus activity. The interstimulusinterval was approximately 500 msec. The average responses were digitally filtered with a bandpass filterbetween 2 and 40 Hz to eliminate high frequency noise, low frequency artifacts, and baseline drifts. Ateach time point between 30 and 75 msec following stimulus onset, dipole coordinates were obtainedusing a single-dipole model, and the ones with the best correlation between the measured and predictedfields were selected as the source location. Following this approach detailed mapping of the sensoryhomunculus can be achieved. In the case of EEG-based mapping, evoked responses are typically obtainedusing electrical stimulation. Figure 12.9 shows an example mapping, where left cortical areas are activatedafter stimulation of the subject’s right toe, little finger, index finger, thumb, and lower lip. The topographicarrangement of the areas identified resembles the “homunculus” typically seen in textbooks.

FIGURE 12.7 The N1 evoked potential resulting from auditory stimulation as recorded from 256 channels aroundthe head (left), surface distribution of potentials at the N1 peak, and cortical areas activated during the highlightedportion of the N1 (right).

FIGURE 12.8 Localization of the auditory N1 component in the primary auditory area resulting from 1-kHz tonestimuli.



12.6.3 Visual Evoked Responses

In cases in which a lesion involves the occipital cortex, visual evoked responses can be used to assesspreoperatively potential postoperative complications. Figure 12.10 shows an example of the cortical areasactivated after stimulation of a patient’s left hemifield of view with a checkerboard pattern (check size1.4∞, pattern reversal rate 1/sec). The maximum magnetic responses were obtained in the right hemi-sphere, at a latency of approximately 130 msec; this is indicated by a square in the figure.

12.6.4 Epileptogenic Spike Localization

Ictal events, such as spikes, are very important for localizing epileptogenic regions in the brain. During thepast several years, MEG has been used in clinical settings as a noninvasive method for localizing the sourcesof ictal activity (Sutherling and Barth, 1989), and more recently, due to the availability of large-array sensors,

FIGURE 12.9 The cortical areas activated after stimulation of a patient’s right toe (1), little finger (2), index finger(3), thumb (4), and lower lip (5) follow the well-known topography of a textbook “homunculus.”

FIGURE 12.10 Cortical areas activated after stimulation of a patient’s left hemifield of view with a checkerboardpattern stimulus.



the sources of interictal spikes (Aung et al., 1995; Papanicolaou et al., 1999). Several studies have suggestedthat noninvasive source imaging may reduce dependence on invasive electrophysiology for localization ofepileptic foci (Papanicolaou et al., 1999; Smith et al., 1993; Zouridakis et al., 1999). Most MEG epilepsystudies report localizations of interictal activity only, because the large artifacts from muscle activity typicallyassociated with seizures make localization of ictal sources extremely difficult. As an example, 2 sec ofsimultaneous MEG and EEG activity is shown in Figure 12.11 where, in addition to the artifacts resultingfrom the electrical activity of the heart, an interictal epileptic spike can be seen in both modalities.

Source localizations of several spikes obtained during this study are overlaid on MRI scans in Figure12.12, where they are shown as triangles. In this case, a focal area of interictal epileptic activity over thelateral surface of the left posterior temporal lobe can be seen.

12.6.5 Language Mapping

Identifying brain regions involved in language functions is often an integral part of the evaluation thatbrain surgery patients. At present, this task is accomplished mostly through invasive techniques (Ojemannet al., 1989), but recently noninvasive techniques based on MEG have also been used. For example, usinga task for continuous recognition of either printed or spoken single words, evoked responses wererecorded for 1 sec after the onset of a word and the corresponding intracranial generators were modeledas single-current dipoles. The number of dipoles obtained in each hemisphere was used to quantify theextent of cerebral activation (Breier et al., 2000; Zouridakis et al., 1998a). The areas thus identifiedincluded the posterior part of the superior temporal gyrus, and the supramarginal and angular gyri. Acharacteristic MRI slice with the localized sources superimposed is shown in Figure 12.13, where thecircles and triangles correspond to two repetitions of the same experimental task.

12.7 Concluding Remarks

The previous paragraphs give some examples of how SI is used for research and clinical purposes. BothMEG and EEG large-array systems have received Food and Drug Administration (FDA) approval for

FIGURE 12.11 An epileptic interictal spike detected simultaneously on the MEG and EEG records, along withECG artifacts.

MEG spike

EEG spike

ECG

I



clinical use, and they are gradually becoming available in many centers around the world. As newtechnological and computational advances come into existence, the clinical relevance of this techniquebecomes more apparent. In particular, in the case of epilepsy surgery, the usefulness of SI as a noninvasivetool to preoperatively delineate the extent of a lesion (to be resected) and of the eloquent cortex (to bepreserved) has already been recognized (Breier et al., 2000; Gallen et al., 1994; Hamalainen et al., 1993;Papanicolaou et al., 1999; Zouridakis et al., 1998a, 1998b).

With the development of larger sensor arrays and of general source modeling algorithms, it is possiblethat future clinical applications of source imaging may be extended to include brain injury and strokeassessment, dementias, developmental disorders, as well as further characterization of higher cortical

FIGURE 12.12 A cluster of epileptic interictal spikes localized in the left posterior temporal lobe.

FIGURE 12.13 Cortical areas involved with receptive language. The circles and triangles correspond to two repeti-tions of the same experimental task.



areas involved in attention, memory, and cognition. As a result, these procedures may lead to safer, faster,and more cost-effective clinical interventions.

References

Aung, M., Sobel, D.F., Gallen, C.C., et al., Potential contribution of bilateral magnetic source imaging tothe evaluation of epilepsy surgery candidates, Neurosurgery, 37, 1113–1121, 1995.

Babiloni, F., Babiloni, C., Carducci, F., Fattorini, L., Onorati, P., and Urbano, A., Spline Laplacian estimateof EEG potentials over a realistic magnetic resonance-constructed scalp surface model, Electroen-ceph. Clin. Neurophysiol., 98(4), 363–373, 1996.

Baumgartner, C., Sutherling, W., Di, S., and Barth, D.S., Spatiotemporal modeling of cerebral evokedmagnetic fields to median nerve stimulation, Electroenceph. Clin. Neurophysiol., 79, 27–35, 1991.

Benzel, E.C., Lewine, J.D., Bucholz, R.D., and Orrison, W.W., Jr., Magnetic source imaging: a review ofthe Magnes system of biomagnetic technologies incorporated, Neurosurgery, 33(2), 252–259, 1993.

Binder, J.R., Frost, J.A., Hammeke, T.A., Cox, R.W., Rao, S.M., and Prieto, T., Human brain areas identifiedby functional magnetic resonance imaging, J. Neurosci., 17, 353–362, 1997.

Binder, J.R., Rao, S.M., Hammeke, T.A., Frost, J.A., Bandettini, P.A., Jesmanowicz, A., and Hyde, J.S.,Lateralized human brain systems demonstrated by task subtraction functional magnetic resonanceimaging, Arch. Neurol., 52, 593–601, 1995.

Brandt, M.E., Topographic mapping of brain electromagnetic signals: a review of current technology,Am. J. Physiol. Imaging, 7(3–4), 160–174, 1992.

Breier, J.I., Simos, P.G., Zouridakis, G., and Papanicolaou, A.C., Lateralization of activity associated withlanguage function using magnetoencephalography: a reliability study, J. Clin. Neurophysiol., 17(5),503–510, 2000.

Cohen, D. and Cuffin, B.N., EEG versus MEG localization accuracy: theory and experiment, Brain Topogr.,4(2), 95–103, 1991.

Cuenod, C.A., Brookheimer, S.Y., Hertz-Pannier, L., Zeffiro, T.A., Theodore, W.H., and Le Bihan, D.,Functional MRI during word generation using conventional equipment: a potential tool for lan-guage localization in the clinical environment, Neurology, 45, 1821–1827, 1995.

Cuffin, B.N., EEG dipole source localization, IEEE Eng. Med. Biol. Mag., 17(5):118–122, 1998.de Munck, J., The potential distribution in a layered anisotropic spherical volume conductor, J. Appl.

Phys., 64, 464–470, 1988.Desmedt, J.E., Nguyen, T.H., and Bourguet, M., Bit-mapped color imaging of human evoked potentials

with reference to the N20, P22, P27 and N30 somatosensory response, Electroenceph. Clin. Neu-rophysiol., 68, 1–19, 1987.

Duffy, F.H., Burchfiel, J.L., and Lomroso, C.T., Brain electrical activity mapping (BEAM): a method forextending the clinical utility of EEG and evoked potential data, Ann. Neurol., 5, 309–321, 1979.

Ebersole, J.S., EEG dipole modeling in complex partial epilepsy, Brain Topogr., 4, 113–123, 1991.Ebersole, J.S. and Wade, P.B., Spike voltage topography and equivalent dipole localization in complex

partial epilepsy, Brain Topogr., 3, 21–34, 1990.Fender, D.H., Models of the human brain and the surrounding media: their influence on the reliability

of source localization, J. Clin. Neurophysiol., 8(4), 381–390, 1991.Flynn, E.R., Factors which affect spatial resolving power in large array biomagnetic sensors, Rev. Sci.

Instrum., 65, 922–935, 1994.Freeman, W., Use of spatial deconvolution to compensate for distortion of EEG by volume conductor,

IEEE Trans. Biomed. Eng., 27, 421–429, 1980.Friberg, L., Brain mapping in thinking and language function, Acta Neurochir. Suppl., 56, 34–39, 1993.Gallen, C.C., Schwartz, B., Rieke, K., Pantev, C., Sobel, D., Hirschkoff, E., and Bloom, F.E., Intrasubject

reliability and validity of somatosensory source localization using a large array biomagnetometer,Electroencephal. Clin. Neurophysiol., 90, 145–156, 1994.



Gomez-Tortosa, E., Martin, E.M., Sychra, J.J., and Dujovny, M., Language-activated single-photon emis-sion tomography imaging in the evaluation of language lateralization-evidence from a case ofcrossed aphasia: case report, Neurosurgery, 35, 515–519, 1994.

Gorodnitsky, I.F., George, J.S., and Rao, B.D., Neuromagnetic source imaging with FOCUSS: a recursiveweighted minimum norm algorithm, Electroencephalogr. Clin. Neurophysiol., 95(4), 231–251, 1995

Hamalainen, M., Hari, R., IImoniemi, R., Knuutila, J., and Lounasmaa, O., Magnetoencephalography —theory, instrumentation, and applications to noninvasive studies of the working human brain, Rev.Mod. Phys., 65, 413–497, 1993.

Hari, R., Pelizzone, M., Makela, J.P., Hallstrom, J., Leinonen, L., and Lounasmaa, O.V., Neuromagneticresponses of the human auditory cortex to on- and off-sets of noise bursts, Audiology, 26, 31–43,1987.

He, B., Chernyak, Y., and Cohen, R.J., An equivalent body surface charge model representative three-dimensional bioelectrical activity, IEEE Trans. Biomed. Eng., 42(7), 637–646, 1995.

He, B., Wang, Y., Pak, S., and Ling, Y., Cortical source imaging from scalp electroencephalograms, Med.Biol. Eng. Comput., 34(Suppl. 2), 257–258, 1996.

Henderson, C.J., Butler, S.R., and Glass, A., The localization of equivalent dipoles of EEG sources by theapplication of electrical field theory, Electroenceph. Clin. Neurophysiol., 39, 117–130, 1975.

Hjorth, B., An on-line transformation of EEG scalp potentials into orthogonal source derivations, Elec-troenceph. Clin. Neurophysiol., 39, 526–530, 1975.

Ioannides, A.A., Bolton, J.P.R., and Clarke, C.J.S., Continuous probabilistic solutions to the biomagneticinverse problem, Inverse Problems, 6, 523–542, 1990.

Iyer, D., Boutros, N.N., and Zouridakis, G., Independent Component Analysis of Multichannel AuditoryEvoked Potentials, 2nd Joint EMBS/BMES Conference, Houston, October 2002.

Jayakar, P., Duchowny, M., Resnick, T., and Alvarez, L.. Localization of seizure foci: pitfalls and caveats,J. Clin. Neurophysiol., 8, 414–431, 1991.

Jirsa, V.K., Fink, P., Foo, P., and Kelso, J.A.S., Parametric stabilization of biological coordination: atheoretical model, J. Biol. Phys., 6, 85–112, 2000.

Kandel, E.R., Schwartz, J.H., and Jessell, T.M., Eds., Principles of Neural Science, 3rd ed., Elsevier, NewYork, 1991.

Kavanagh, R.N., Darcey, T.M., Lehmann, D., and Fender, D.H., Evaluation of methods for three-dimen-sional localization of electrical sources in the human brain, IEEE Trans. Biomed. Eng., 25, 421–429,1978.

Kelly, J.P., The neural basis of perception and movement, in Principles of Neural Science, 3rd ed., Kandel,E.R., Schwartz, J. H., and Jessell, T.M., Eds., Elsevier Science, New York, 1991.

Le, J., Menon, V., and Gevins, A., Local estimate of surface Laplacian derivation on a realistically shapedscalp surface and its performance on noisy data, Electroenceph. Clin. Neurophysiol., 90, 433–441,1994.

Lewine, J.D., Neuromagnetic techniques for the noninvasive analysis of brain function, in NoninvasiveTechniques in Biology and Medicine, Freeman, S.E., Fukushima, E., and Greene, E.R., Eds., SanFrancisco Press, San Francisco, 1990.

Lopes da Silva, F.H., Wieringa, H.J., and Peters, M.J., Source localization of EEG versus MEG: empiricalcomparison using visually evoked responses and theoretical considerations, Brain Topogr., 4(2),133–142, 1991.

Mejis, J.W.H. and Peters, M.J., The EED, EEG and MEG, using a model of eccentric spheres to modelthe head, IEEE Trans. Biomed. Eng., BME-34, 913–920, 1987.

Mejis, J.W.H., Weier, O.W., Peters, M.J., and van Oosterom, A., On the numerical accuracy of theboundary element method, IEEE Trans. Biomed. Eng., BME-36, 1038–1049, 1989.

Morse, P.M. and Feshbach, H., Methods of Theoretical Physics I and II, McGraw Hill, New York, 1963.Mosher, J.C. and Leahy, R.M., Recursive MUSIC: a framework for EEG and MEG source localization,

IEEE Trans. Biomed.. Eng., 45(11), 1342–1354, 1998.



Mosher, J.C., Lewis, P.S., and Leahy, R., Multiple dipole modeling and localization from spatiotemporalMEG data, IEEE Trans. Biomed. Eng., 39, 541–557, 1992.

Nishijo, H., Ikeda, H., Miyamoto, H., Endo, S., and Ono, T., Localization of an ictal onset zone using arealistic 4-shell head model of scalp, skull, liquor, and brain, Soc. Neurosci. Abst., 22, 185, 1996.

Nunez, P. and Pilgreen, K., The spline Laplacian in clinical neurophysiology: a method to improve theEEG spatial resolution, J. Clin. Neurophysiol., 38, 356–368, 1975.

Nunez, P.L., Pilgreen, K.L., Westdorp, A.F., Law, S.K., and Nelson, A.V., A visual study of surface potentialsand Laplacians due to distributed neocortical sources: computer simulations and evoked potentials,Brain Topogr., 4(2), 151–168, 1991.

Ojemann, G., Ojemann, J., Lettich, E., and Burger, M., Cortical language localization in left, dominanthemisphere: an electrical stimulation mapping investigation in 117 patients, J. Neurosurg., 71,316–326, 1989.

Okada, Y., Huang, J.C., and Xu, C., A hierarchical minimum norm estimation method for reconstructingcurrent densities in the brain from remotely measured magnetic fields, in Biomagnetism, ClinicalAspects, Hoke, M., Ernd, S.N., Okada, Y.C., and Romani, G.L., Eds., Excerpta Medica, Amsterdam,1992, pp. 729–734.

O’Leary, D.S., Andreason, N.C., Hurtig, R.R., Hichwa, R.D., Watkins, G.L., Ponto, L.L., Rogers, M., andKirchner, P.T., A positron emission tomography study of binaurally and dichotically presentedstimuli: effects of level of language and directed attention, Brain Lang., 53, 20–39, 1996.

Papanicolaou, A.C., Simos, P.G., Breier, J.I., Zouridakis, G., Willmore, L.J., Wheless, J.W., Constantinou,J.E., Maggio, W.W., and Gormley, W.B., Magnetoencephalographic mapping of the language-specific cortex, J. Neurosurg., 90(1), 85–93, 1999.

Pascual-Marqui, R.D., Michel, C.M., and Lehmann, D., Low resolution electromagnetic tomography: anew method for localizing electrical activity in the brain, Int. J. Psychophysiol., 18(1), 49–65, 1994.

Peterson, S.E., Fox, P.T., Posner, M.I., Mintun, M., and Raichle, M.E., Positron emission tomographicstudies of cortical anatomy of single-word processing, Nature, 1988, 331, 585–589, 1988.

Plonsey, R., The nature of sources of bioelectric and biomagnetic fields, Biophys. J., 39, 309–312, 1982.Plonsey, R. and Heppner, D., Considerations of quasistationarity in electrophysiological systems, Bull.

Math. Biophys., 29, 657–664, 1967.Rush, S. and Driscoll, D., Current distribution in the brain from surface electrodes, Anesth. Analg., 47,

717–721, 1978.Scherg, M., Functional imaging and localization of electromagnetic brain activity, Brain Topogr., 5,

103–112, 1992.Scherg, M. and von Cramon, D., Two bilateral sources of the late AEP as identified by a spatio-temporal

dipole model, Electroenceph. Clin. Neurophysiol., 62, 32–44, 1985.Schneider, M.R., A multistage process for computing virtual dipolar sources of EEG discharges from

surface information, IEEE Trans. Biomed. Eng., 19, 1–12, 1972.Sidman, R., Vincent, D., Smith, D., and Lu, L., Experimental test of the cortical imaging technique-

applications to the response to median nerve stimulation and the localization of epileptiformdischarges, IEEE Trans. Biomed. Eng., 27, 437–444, 1992.

Sidman, R.D., Giambalvo, V., Allison, T., and Bergey, P., A method for localization of sources of humancerebral potentials evoked by sensory stimuli, Sensory Processes, 2, 116–129, 1978.

Smith, D., Baker, G., Davies, G., Dewey, M., and Chadwick, D.W., Outcomes of add-on treatment withlamotrigine in partial epilepsy, Epilepsia, 34(2), 312–322, 1993.

Supek, S. and Aine, C.J., Simulation studies of multiple dipole neuromagnetic source localization: modelorder and limits of source resolution, IEEE Trans. Biomed. Eng., 40, 529–540, 1993.

Sutherling, W.W. and Barth, D.S., Neocortical propagation in temporal lobe spike foci on magnetoen-cephalography and electroencephalography, Ann. Neurol., 25(4), 373–381, 1989.

Wang, J., Williamson, S.J., and Kaufman, L., Magnetic source images determined by a lead-field analysis:the unique minimum-norm least squares estimation, IEEE Trans. Biomed. Eng., 39, 665–675, 1992.



Williamson, S.J. and Kaufman, L., Analysis of neuromagnetic signals, in Handbook of Electroencephalog-raphy and Clinical Neurophysiology, Vol. 1, Gevins, A.S. and Remond, A., Eds., Elsevier, New York,1987, pp. 405–448.

Wilson, F.N. and Bayley, R.H., The electric field of an eccentric dipole in an homogeneous sphericalconducting medium, Circulation, 1, 84–92, 1950.

Witwer, J., Trezek, G., Jewett, D.L., and Yamashita, Y., The effect of media inhomogeneities upon intrac-ranial electrical fields, IEEE Trans. Biomed. Eng., BME-19, 352–362, 1972.

Yeh, G.C.K. and Martinek, K., The potential of an general dipole model in a homogeneous conductingprolate spheroid, Ann. N.Y. Acad. Sci., 65, 1003–1006, 1957.

Zhou, H. and van Oosterom, A., Computation of the potential distribution in a four-layer anisotropicconcentric spherical volume conductor, IEEE Trans. Biomed. Eng., BME-38, 154–158, 1991.

Zouridakis, G., Simos, P.G., and Papanicolaou, A.C., Multiple bilaterally asymmetric cortical sourcesaccount for the auditory N1m component, Brain Topogr., 10, 183–189, 1998a.

Zouridakis, G., Simos, P.G., Breier, J.I., and Papanicolaou, A.C., Functional hemispheric asymmetryassessment in a visual language task using MEG, Brain Topogr., 11, 57–65, 1998b.

Zouridakis, G., Simos, P.G., Papanicolaou, A.C., and Breier, J.I., Magnetic source imaging: introductionand application examples, J. Clin. Eng., 24, 51–61, 1999.


Functional Brain Mapping Using Intracranial Source Imagingengold.ui.ac.ir/~mahnam/IBME/Research/Functional Brain Mapping.pdf · Functional Brain Mapping Using Intracranial Source

Documents