Extracting Visual Evoked Potentials from EEG Data Recorded During fMRI-guided Transcranial Magnetic Stimulation

Summary

This paper describes a method for collecting and analyzing electroencephalography (EEG) data during concurrent transcranial magnetic stimulation (TMS) guided by activations revealed with functional magnetic resonance imaging (fMRI). A method for TMS artifact removal and extraction of event related potentials is described as well as considerations in paradigm design and experimental setup.

Abstract

Transcranial Magnetic Stimulation (TMS) is an effective method for establishing a causal link between a cortical area and cognitive/neurophysiological effects. Specifically, by creating a transient interference with the normal activity of a target region and measuring changes in an electrophysiological signal, we can establish a causal link between the stimulated brain area or network and the electrophysiological signal that we record. If target brain areas are functionally defined with prior fMRI scan, TMS could be used to link the fMRI activations with evoked potentials recorded. However, conducting such experiments presents significant technical challenges given the high amplitude artifacts introduced into the EEG signal by the magnetic pulse, and the difficulty to successfully target areas that were functionally defined by fMRI. Here we describe a methodology for combining these three common tools: TMS, EEG, and fMRI. We explain how to guide the stimulator's coil to the desired target area using anatomical or functional MRI data, how to record EEG during concurrent TMS, how to design an ERP study suitable for EEG-TMS combination and how to extract reliable ERP from the recorded data. We will provide representative results from a previously published study, in which fMRI-guided TMS was used concurrently with EEG to show that the face-selective N1 and the body-selective N1 component of the ERP are associated with distinct neural networks in extrastriate cortex. This method allows us to combine the high spatial resolution of fMRI with the high temporal resolution of TMS and EEG and therefore obtain a comprehensive understanding of the neural basis of various cognitive processes.

Introduction

Transcranial Magnetic Stimulation (TMS) generates momentary interference to the normal neural activity in target areas of the brain. By creating this transitory neural interference and measuring a behavioral or physiological change, we can draw a causal link between the target area and the measured experimental effect (for a review see Pascual-Leone et al. and Taylor et al.^1,2). Such an experimental effect may be, for example, a performance on a cognitive task or a change in electrophysiological (EEG) activity. Indeed, in recent years researchers have started using TMS in combination with EEG to directly relate cortical areas with event-related potentials (ERP) or oscillatory activity patterns (e.g.^2-7). In this methodological paper we will describe a particular and useful framework for combining TMS and EEG: fMRI-guided TMS during an ERP experiment. First, we will detail how to apply TMS to areas predefined by fMRI, while recording EEG data. We will then describe an experimental design that allows extraction of reliable ERP. The goal of such an experiment is to causally link brain areas revealed with functional MRI to ERP components of interest. Lastly, we will give a specific example of a study relating face and body selective ERPs with face and body selective areas that are revealed with fMRI.

What is the benefit of linking EEG signals with fMRI activations? EEG and fMRI are commonly used tools to measure cortical responses to visual input. For example, category-selectivity in the visual pathway was assessed for different visual object categories such as faces, body parts, and written words, both by means of ERP extracted from EEG data^8,9, and functional MRI^10-12. The signals measured by these two common research tools are, however, of fundamentally different nature. EEG carries information about neural electrical activity with great temporal precision, but very low spatial resolution and may reflect a mixture of many separate underlying sources. The fMRI provides an indirect measure of neuronal activity relying on the slow hemodynamic changes occurring during stimulus presentation or/and task execution, but presents this activity with a higher spatial resolution. Establishing a correlation between the two measures can thus be of great interest, but is limited in that it does not imply a causal link between the scalp-recorded electrophysiological response and the areas revealed with functional MRI. Even when measured simultaneously (e.g.^13-15), a directional causal relation between EEG and activity in functionally defined cortical areas cannot be determined. TMS is a tool that can assist achieving the establishment of such a causal relationship.

A simultaneous EEG-TMS study is methodologically challenging, mostly due to the high voltage artifact introduced to the EEG signal by the magnetic stimulation (see Figure 1, for a review see Ilmoniemi et al.¹⁶). This artifact consists of a transient short living pulse-related disturbance, often followed by a slower secondary (or residual) artifact that may last a few hundreds of milliseconds after the pulse is delivered Figure 2A, thus overriding most ERP components of interest. This secondary artifact may include mechanical sources such as currents induced by the magnetic pulse into the wiring and the slow decay of these currents in the skin, and physiological sources such as muscular activity over the scalp and auditory or somatosensory evoked potentials elicited by the operation of the coil^17-20. Although the mechanical sources of interference probably produce larger amplitude artifacts than the physiological ones, these different artifacts cannot be separated, and the existence of any of them in the signal can confound the results. One possible solution is the application of repetitive TMS pulses prior to EEG recording ("offline TMS"), as opposed to simultaneous EEG-TMS. The inhibitory effect of such a protocol on cortical activity persists for several minutes (and up to half an hour) after the stimulation, and EEG can be measured during this effective time window and compared with baseline, pre-TMS, EEG data. Repetitive stimulation, however, is by definition lacking the high temporal resolution that online TMS may offer, where pulses can be administered at a precise timing relative to trial onset at the millisecond resolution. The effect of repetitive stimulation may also propagate via cortical connections across a wider area than desired and therefore significantly reduce the spatial resolution as well.

To take advantage of both the spatial and temporal resolution that TMS can provide, a simultaneous EEG-TMS combination can be applied. However, this requires methods for removal of artifacts generated by the magnetic stimulation on the EEG signal. Very few offline mathematical solutions for TMS artifact removal have been proposed^16,21,22, although no method is agreed upon, and no one method may be optimal for all experimental designs. A "clipping" system, consisting of a sample-and-hold circuitry, was also developed to momentarily stop EEG acquisition during TMS pulse delivery 20. This technique not only requires specialized hardware, but may not completely remove the residual TMS artifact. In this paper we will describe an adaptation of an EEG-TMS methodology developed by Thut and colleagues¹⁹, particularly suitable for ERP studies. This technique allows reliable extraction of ERP while eliminating all the residual noise components caused by the TMS pulse Figure 2. We will further provide general guidance towards a successful EEG-TMS experimental setup.

Another challenge in TMS studies addressed in this methodological paper is finding the best coil position and angle for an accurate targeting of the desired cortical area. We will describe the use of a stereotactic navigation system to coregister the subject's head with the pre-acquired functional MRI images. Although the navigation system can be used to localize anatomically-defined brain structures, an fMRI-guided targeting is particularly useful since for many functions and experimental effects the precise location of activation cannot be inferred from anatomical markers alone. For such functional regions of interest (ROI), the definition of an area is made for each participant individually.

To illustrate all of the above, we will provide an example of a study we conducted previously, in which EEG was recorded concurrently with TMS guided by fMRI activations⁷. In this study, a double dissociation was made between face-selective and body-selective ERPs: although face and body ERPs peak around the same latency and electrode sites, targeting individually defined face-selective and body-selective areas in the lateral occipital lobe enabled us to dissociate the neural networks underlying each ERP response. Finally, we will try to give more general advise for optimizing EEG recording during TMS application.

Protocol

The experiment is conducted in two separate sessions. During the first session a functional MRI experiment (e.g. a functional localizer) is carried out in order to define the desired TMS target areas on an individual subject basis. The fMRI results are then fed into a stereotactic navigation system for an accurate TMS targeting. The second session is held following the analysis of the fMRI data, during which EEG is recorded concurrently with TMS. The protocol described here was approved by the ethics committee of the Tel-Aviv Sourasky Medical Center.

In the example given in this paper, data were analyzed with MATLAB version 7.7 (R2008b).  The Statistical Parametric Mapping (SPM 5) software for MATLAB and MarsBar toolbox for SPM 23 were used for fMRI data processing.

1. fMRI Session and fMRI Data Analysis

1. Run a functional MRI task using an echo planar imaging (EPI) sequence to determine the desired foci of activation to be targeted with TMS. For a better dissociation between adjacent areas, such as the OFA and EBA in the example below, high-resolution scanning is recommended. Voxels of 3 mm³ or smaller, obtained with an MRI head-coil of 8 channels or more, are sufficient for delineating adjacent regions.
2. Run a T1-weighted structural scan to obtain neuroanatomical data. Make sure that the face of the participant is entirely included in the field of view of this scan, since external markers on the image (e.g. the tip of the nose) will be later used to coregister the subject's head with its scan.
3. After data collection, use the MarsBar toolbox for SPM to define the desired brain regions of interest based on contrasts between the experimental conditions. Use the contrasts face>objects to define the Occipital Face Area (OFA), and bodies>objects to define the Extrastirate Body Area (EBA). To further ensure that the two cortical targets are functionally distinct, use "conjunction" (in MarsBar) to mask from each ROI any voxels responding to the other experimental condition (exclude face-voxels from the EBA, and body-voxels from the OFA).
4. Corregister the structural T1 images with the functional scans, using SPM.
5. Copy the files of the structural scan, as well as the relevant functional contrast files, to a portable drive in order to upload to the navigation system.

2. Preparing a Paradigm for the EEG-TMS Experiment that will allow ERP Extraction

Described in the section below is a method for collecting EEG data during TMS application in a way that allows extraction of reliable and replicable ERPs¹⁹. The advantage of this technique is that it easily handles the secondary, long lasting, TMS artifact, and is robust enough even to allow restoration of data at electrodes located right under the TMS coil, where the artifact is of the highest voltage and longest duration.

1. Organization of the paradigm
   1. Run the different TMS conditions (the different target brain areas, as well as a no-TMS condition) in separate blocks.
   2. Within each block present the participant with all stimulus conditions (e.g. faces, objects, scenes and so forth) randomly in an event-related design.
   3. For better quality of ERP and TMS-noise template (below) make sure to have at least 50 trials per condition.
2. Set the timing of the TMS pulse/pulses to the desired latency after image onset. This is done via writing to a parallel port, from which a cable goes to the TMS stimulator. This function is available in most software for psychological experiments, such as Psychtoolbox (versions 2 or 3) for MATLAB 24 or E-Prime (see Materials Table). Jitter the inter stimulus interval (ISI) to reduce stimulus (and pulse) predictability (e.g. add a random value between 0 to 500 msec at each ISI).
3. Prepare an additional blank-screen condition:
   1. Prepare trials during which TMS will be applied at the same intensity, but with no stimulus presentation on the screen. These blank-screen TMS trials will be served to calculate a TMS artifact template in the absence of visual stimulation.
   2. Set the number of repetitions of the blank trials to be identical to the number of repetitions of any one of the experimental conditions within the block.
   3. For an accurate representation of the shape of the TMS residual artifact, randomize the blank trials throughout the whole block rather than presenting them all at the beginning or at the end.

3. Setting up the EEG and the Neuronavigation System, and Conducting the Experiment

Accurate TMS targeting of individually defined ROIs is possible with the use of a stereotactic navigation system, comprised of an infrared camera, infrared sensors mounted on the participant's head, and a specialized software.

1. Screen participants based on TMS safety criteria. Exclude from participation subjects with a self or family history of epilepsy, subjects with other neurological conditions or with frequent migraines, and subjects on psychoactive drugs. Although not typically screened for, subjects with suspicion for autonomic disorders such as vasovagal syncopes (mostly manifested as a tendency to faint easily) may also be excluded. Instruct participants to avoid alcoholic drinks starting the prior evening, and caffeinated drinks at least 2 hr before the experiment.  For further instructions and discussion of safety see Rossi et al²⁵, and the Magstim safety review (http://joedevlin.psychol.ucl.ac.uk/tms/docs/magstim_safety.pdf).
2. Prepare the neuronavigation system:
   1. Before the session starts, feed the structural scan files into the navigation system's software.
   2. Overlay the functional MRI results (the contrasts) on the structural images.
   3. Using the Neuronavigation software, mark the desired targets on the images, as well as the external anatomical markers that will serve for coregistration: the tip of the nose, the deepest part of the nose bridge often referred to as nasion, and the tragus of each ear.
3. Mount the EEG cap on the participant's head and connect the electrodes:
   1. Try to keep electrode impedance not higher than 5 kΩ.
   2. To avoid TMS-related heating of the electrodes, use as little gel as possible. To achieve good impedance with little amount of gel perform a thorough skin preparation. Optionally, ask participants to wash their hair before coming to the experiment.
   3. Make sure that the electrode wires do not cross each other and are oriented away from the coil's location. Avoid loops in the wires.
   4. Use a high sampling rate for a better representation of the noise artifact. It is recommended to use 1 kHz or above, as most previous studies using this method have done^7,26-28.
   5. Place the reference and ground electrodes as far from the coil as possible. In this example, areas in the occipital cortex were targeted using a nose reference and an Fz ground⁷. For other examples see^3,4,27,29,30. Note that data can be re-referenced offline to a new reference as needed, such as the common average.
      Note: For a review about optimization of TMS-EEG setup, see Veniero et al³¹.
4. Coregister the subject's head with the scan, as follows:
   1. Mount the infrared detectors on the participant's head.
   2. Coregister the head location with the navigation system using the predefined markers (tip of the nose, etc. See Figure 3). It is advised to repeat the coregistration between blocks to ensure accurate coil placement at all stages.
5. Locate the target areas:
   1. Have the subject seated with the chin resting on a chinrest at the desired distance from the screen.
   2. Make sure the participants are comfortable at their chair, as they are asked to refrain from movements during the experimental blocks (important for an accurate noise-template measurement).
   3. Choose a TMS target from the navigation system Figure 3.
   4. Using the pointer tool (see Materials Table), let the navigation system guide the user to the optimal coil location and mark it with a small sticker on the electrode cap. It is important to hold the pointer perpendicular to the head. Repeat this stage prior to each block. Note that it is not recommended to use online navigation during the block (navigating the coil itself while holding it) since any coil movement should be avoided for a best measurement of the TMS artifact template. It was found that predefining and marking the TMS location is the optimal way to achieve stable stimulation.
   5. Guide the coil's precise center to the marker, while being held by a holder. Make sure it is tangential to the head.
6. Set the TMS intensity to the desired value. Administer one test pulse for the participant's approval.
7. Run the experimental block.
8. No-TMS condition: if a specialized sham TMS coil is not available, place the TMS coil next to the subject's head and tilt it in 90°. Run the block as usual, including the blank trials.
   The reader may also refer to the JoVE video paper by Andon and Zatorre³² for further demonstration of the navigation system.

4. Analyzing the EEG Data and Computing ERP

1. Remove the immediate pulse artifact, as follows:
   1. If a clipping device is not available (see above), the first step in the EEG data processing would be to cut out the immediate TMS pulse artifact itself from the data. Note that this step can be skipped if filters are not desired. But if filters are applied, the sharp-edge shape of the artifact will create distortions in the data. A narrow time-window of 10 to 15 msec after pulse onset should suffice, but be sure to verify this by visual inspection of the data.
   2. Connect the two cut ends created after pulse removal. Two main methods for accomplishing this have been suggested in previous reports: 1. simply join together the two ends remaining after pulse removal (see Fuggetta et al.²⁶ and Figure 1); 2. interpolate a line between the two cut ends by generating equally spaced values between them⁷. The equation used for this interpolation is as follows: for each missing data point y at sample x, compute y=y₀+((y₁-y₀)*(x-x₀))/(x₁-x₀) where x₀ and y₀ are the coordinates for the last data point before the cut segment, and x₁ and y₁ are the coordinates for the first data point after the cut segment. Both techniques assure that filters can be applied on the data without producing ripples due to sharp voltage steps, as demonstrated in Figure 1. Refer to the work by Reichenbach and colleagues²⁷ for a more complex 3^rd order polynomial interpolation.
2. Apply the subtraction method:
   1. For each experimental block, including a no-TMS control condition if present, calculate an averaged ERP to the blank screen trials by time-locking them to trial start (as if an image had been presented).
   2. Subtract this averaged template from each trial of all other stimulus conditions. If several blocks were run for the same stimulation site, do this separately for each block, as the template will be slightly different between the blocks.
3. All other preprocessing and processing steps are carried out as in any other ERP experiment.

Results

A concurrent EEG-TMS investigation was used to reveal whether the ERP responses to faces and bodies recorded over the occipito-temporal scalp are dissociated. When visual stimuli are presented, a prominent N1 component is recorded at posterior-lateral electrode sites. In particular, the N1 component is typically larger for faces and bodies than to other stimulus categories^8,33. By assessing the effect of stimulation on face and body-selective brain areas defined with fMRI on their respective face and body N1 c...

Discussion

Having the unique ability to momentarily disrupt normal neuronal activity in selected cortical areas, in precise time points and with a relatively good spatial accuracy, TMS allows to causally link a stimulated brain area with a behavioral or a neurophysiological measure. In this paper we described a method for measuring EEG during concurrent TMS application, targeting functionally defined cortical areas, and applying an analysis that enables a reliable measuring of ERP responses. We gave an example from the literature i...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
3.0 T Signa MRI scanner	General Electric
BrainAmp amplifier	Brain Products GmbH	BP-01300
Electrode input box	Brain Products GmbH		Optional
PowerPack - battery for amplifier	Brain Products GmbH	BP-02615
BrainCap - 32 flat electrodes on a flexible cap	Brain Products GmbH	BP-0300MR	Flat electrodes should be used to assure a shorter distance beween coil and scalp. If larger (e.g. pin type) electrodes are used, remove the ones under the coil
TMS Super Rapid2 stimulator	Magstim
50 mm double coil	Magstim
Coil holder			Any mechanical arm or tripod that can hold the coil, be adjusted to the right angle and location, and keep the coil steady during stimulation
Chinrest
Polaris infrared camera	Rogue Research Inc
Polaris trackers and pointer tool	Rogue Research Inc
BrainSight workstation and software	Rogue Research Inc
BrainVision Recorder software	Brain Products GmbH	BP-00010
MATLAB software	The MathWorks Inc
SPM for Matlab	Wellcome Department of Imaging Neuroscience, London, UK
MarsBar region of interest toolbox for SPM
Psychtoolbox for MATLAB			This toolbox and the E-prime software (below) are examples for stimulus presentation software capable of delivering commands to the TMS stimulator and to the EEG recorder with reliable timing
E-Prime software	Psychology Software Tools, Inc.