The Saccadic Re-Centering Bias is Associated with Activity Changes in the Human Superior Colliculus

Being able to effectively explore our visual world is of fundamental importance, and it has been suggested that the straight-ahead gaze (primary position) might play a special role in this context. We employed fMRI in humans to investigate how neural activity might be modulated for saccades relative to this putative default position. Using an endogenous cueing paradigm, saccade direction and orbital starting position were systematically manipulated, resulting in saccades toward primary position (centripetal) and away from primary position (centrifugal) that were matched in amplitude, directional predictability, as well as orbital starting position. In accord with earlier research, we found that fMRI activity in the superior colliculus (SC), as well as in the frontal eye fields and the intraparietal sulcus, was enhanced contralateral to saccade direction across all saccade conditions. Furthermore, the SC exhibited a relative activity decrease during re-centering relative to centrifugal saccades, a pattern that was paralleled by faster saccadic reaction times. In contrast, activity within the cortical eye fields was not significantly modulated during re-centering saccades as compared to other saccade types, suggesting that the re-centering bias is predominantly implemented at a subcortical rather than cortical processing stage. Such a modulation might reflect a special coding bias facilitating the return of gaze to a default position in the gaze space in which retinotopic and egocentric reference frames are aligned and from which the visual world can be effectively explored.

1996; Luna et al., 1998;Koyama et al., 2004). Numerous studies in animals and humans have indicated that the IPS is part of the direct visual pathway for the generation of reflexive saccades and seems to be involved in the transformation of saccade vector coordinates that will finally be transmitted to the SC (e.g., Duhamel et al., 1992;Gaymard et al., 1998). In turn, FEF has been specifically associated with target selection during the generation of voluntary saccades (e.g., Sommer and Tehovnik, 1997;Mort et al., 2003). Beyond the generation of saccadic eye movements, both the FEF and the IPS have also been associated with the covert orienting of attention in a wide range of paradigms, consistent with the view that these two functions are supported by overlapping neural systems (Corbetta and Shulman, 2002).
One important property that appears to be shared by saccadeprocessing regions, including the SC, the FEF, and the IPS, is the contralateral representation of saccade targets based on the underlying retinotopic (eye-centered) reference frame, thus resulting in enhanced activity in the hemisphere contralateral to the targeted visual field. This contralateral predominance (also referred to as contralateral bias) has been demonstrated in a large number of animal studies using single-unit recordings and micro-stimulation measures targeting the SC (Goldberg and Wurtz, 1972;Robinson, 1972;Schiller and Stryker, 1972;Kustov and Robinson, 1996;Moschovakis, 1996;Schiller and Tehovnik, 2001;Ignashchenkova

IntroductIon
The execution of eye movements to objects or locations of interest is carried out by a complex network of cortical and subcortical brain regions enabling a precise alignment of the observer's fovea with the object of interest (for a review see Wurtz and Albano, 1980;Sparks, 2002). One of the crucial structures both for the reflexive orienting toward salient stimuli and for the initiation of voluntary saccadic eye movements is the superior colliculus (SC), a small layered structure that forms part of the tectum of the midbrain. The SC receives multiple cortical and subcortical inputs and represents the final relay structure for eye-movement execution (Fries, 1984;Sparks et al., 2000;Schiller and Tehovnik, 2001;Sparks, 2002;Sommer and Wurtz, 2008). Moreover, the SC has been shown to be involved in the selection of future gaze locations based on a visual salience map (e.g., Fecteau and Munoz, 2006;Boehnke and Munoz, 2008), as well as in covert shifts of attention in the absence of eye movements (e.g., Ignashchenkova et al., 2004).
Among its wide range of connections, the SC receives prominent input from two cortical regions that are linked to both overt saccades and covert attentional shifts: the frontal eye fields (FEF) located in the premotor cortex (Paus, 1996;Blanke et al., 2000;Koyama et al., 2004) and an area within the intraparietal sulcus (IPS; encompassing BA7 and BA40) considered to be the human equivalent of the primate lateral intraparietal area (LIP, Muri et al., et al., 2004;Field et al., 2008) and the FEF (Bruce and Goldberg, 1985;Schall, 1995;Tehovnik and Slocum, 2000;Crapse and Sommer, 2009). In the monkey LIP, although topographically organized, most neurons seem to exhibit a relatively small contralateral bias (Blatt et al., 1990;Platt and Glimcher, 1997). In humans, the contralateral representation of stimulus input and saccade targets has been mostly demonstrated at the cortical level, including FEF, occipital cortex, and IPS (Sereno et al., 2001;Astafiev et al., 2003;Schluppeck et al., 2005;Silver et al., 2005;Kastner et al., 2007;Serences and Yantis, 2007;Curtis and Connolly, 2008;Ikkai and Curtis, 2008;Konen and Kastner, 2008). With the increasing quality of neuroimaging techniques, however, studies have started to investigate the functional contralateral organization of the human SC during visual stimulation and covert attention shifts (DuBois and Cohen, 2000;Kastner, 2005, 2009;Wall et al., 2009) and during saccadic eye movements (Krebs et al., 2010).
The retinotopic reference frame underlying this contralateral representation is separated into the left and right visual field by the vertical meridian that intersects visual fixation. Considering this relationship within an eye-centered reference frame, it might be inferred that the representation of saccade targets is completely independent of the eye position relative to the head, and thus, given a head-restrained setting, that saccade vectors and parameters can be fully described in retinotopic coordinates. However, there is evidence that saccade processing is not independent of the egocentric (head-centered) reference frame. A number of behavioral studies in humans (Laurutis and Robinson, 1986;Zambarbieri et al., 1995;Fuller, 1996), as well as studies in animals (Albano and Wurtz, 1982;McIlwain, 1986;Pare and Munoz, 1996), have shown that saccades toward the center of straight-ahead gaze (primary position) are facilitated as compared to saccades to other gaze locations, as measured by saccadic reaction times (SRT). It has been proposed that this "re-centering bias" (also referred to as centripetal bias) might subserve the efficient coordination of eye-and head-movements by providing a common default position (Fuller, 1996;Pare and Munoz, 2001;Tatler, 2007). This primary position, in which the retinotopic and egocentric reference frame are aligned, is considered to be a particularly efficient starting point for exploring the visual environment in an energy-efficient way (Vitu et al., 2004;Tatler, 2007;Kardamakis and Moschovakis, 2009;Durand et al., 2010).
The neural processes underlying the behavioral findings of a re-centering bias, however, are not well understood. In addition to the well described retinotopic saccade-target representation, such a modulation would require information about the current eye position relative to the current head-position being available to the relevant saccade-controlling brain regions. Indeed, based on studies investigating saccade-related SC activity in animals, there is evidence that the SC might be involved in the processing of eyeposition signals (McIlwain, 1986;Van Opstal et al., 1995;Pare and Munoz, 2001;Campos et al., 2006) and that such a signal might contribute to the facilitation of re-centering the eyes to the primary position. More specifically, it has been proposed that saccades of identical direction and amplitude require different levels of innervation at the motor plant depending on the orbital deviation from primary position and that this might be reflected in a dynamic adaptation of the neuronal activity in the SC (Van Opstal et al., 1995;Campos et al., 2006). With regard to the role of human SC in this context, we recently observed differential activity changes in the SC for saccades toward (centripetal) and away (centrifugal) from primary position (Krebs et al., 2010). However, in our previous paradigm, centripetal (i.e., re-centering) and centrifugal saccades differed in their directional predictability. In addition, since the fMRI acquisition in that study was limited to a small volume optimized for the brainstem, the data provided no information about cortical eye fields, i.e., FEF and IPS, and their possible modulation during re-centering saccades.
In the present study we employed fMRI during a paradigm in which subjects performed endogenously cued 7°-saccades to the left and right between five target positions along the horizontal meridian (see Figure 1). Saccade vectors were oriented toward (centripetal) or away from (centrifugal) the primary central position (orientation factor), starting from positions that were either proximal or distal to the primary gaze position (proximity factor). In addition to the contralateral coding of saccade direction (contralateral predominance factor), this design allowed us to probe the processing level at which the facilitation of saccades toward the primary position primarily occurs.
On the basis of earlier work in animals and humans, we hypothesized that the re-centering bias would be specifically reflected by activity changes within the SC, due to its role in saccade programming and its capacity to incorporate both eye-and head-position information (Van Opstal et al., 1995;Campos et al., 2006). Importantly, by using endogenous cues, the present paradigm enabled us to isolate neural responses in the absence of salient peripheral targets which strongly modulate SC activity (e.g., Boehnke and Munoz, 2008), and may have tended to dilute the neural processing activity related to the re-centering bias. Furthermore, we expected relatively little modulation related to the re-centering bias in cortical areas (FEF and IPS). More specifically, given that all saccades in the current paradigm were endogenously cued, including the re-centering saccades, they would be likely to engage similar levels of top-down cortical control during cue interpretation and voluntary target selection before the saccade-initiation command is passed on to the SC. At this later processing stage, the SC could then efficiently facilitate saccades toward the primary position by incorporating eye-and head-position information.

MaterIals and Methods fMrI subjects and paradIgM
Sixteen healthy right-handed subjects participated in the fMRI study (mean age ± standard deviation (SD): 25 ± 4.6, eight female). Four of initially 20 subjects had to be excluded based on their performance in the saccade task (two subjects) or due to technical difficulties (two subjects). Participants gave written informed consent before the experiment in accordance with the Duke Medical Center Institutional Review Board for human subjects.
We used an event-related fMRI saccade paradigm that manipulated saccade direction (left vs. right), the orientation of the saccade vector (centrifugal vs. centripetal), and the proximity of the saccade vector (proximal vs. distal) relative to the primary position. At the beginning of each experimental run, subjects saw an instruction screen for 6 s that displayed the assignment of the color cues (e.g., green = saccade to the left, blue = saccade to the right, red = maintain fixation). The color-direction mapping was consistent throughout minimize blinking. The timing of cue onsets was pseudo-randomly varied with a stimulus onset asynchrony (SOA) of 1.5-6 s (average SOA = 2.5 s) to allow for effective event-related blood-oxygen level-dependent (BOLD) response estimation (Hinrichs et al., 2000). The event-related design resulted in eight saccade conditions ( Figure 1B), along with the cued fixation condition that served as a baseline in which an instructional cue was presented but no eye movement was performed.
Importantly, the employed design allowed us to compare centripetal saccades and centrifugal saccades (see Figure 1B: proximal centripetal 1-to-0 vs. distal centrifugal 1-to-2) with identical shift amplitude, identical directional probability (i.e., the direction of the saccade was not predictable), and, importantly, with an equivalent pre-saccadic deviation from the primary position (as illustrated in Figure 1C). These saccades differed particularly in the orientation of their movement with respect to the primary position (centripetal vs. centrifugal) and were thus particularly well suited for investigating the neural correlates associated with the re-centering bias. The additional saccade conditions (see Figure 1B: proximal centrifugal 0-to-1 and distal centripetal 2-to-1) were mainly included to provide equal directional probabilities throughout the task. the experiment for each subject, but was counterbalanced across subjects. The instructions were followed by the default stimulus display consisting of five white squares (each 0.5° × 0.5° of visual angle) on a black background, with the central square located at the center of the screen and with an even between-square spacing of 7° ( Figure 1A). In regards to this layout, we will refer to the central position as being at position 0, the adjacent ones as being at position 1 on the left or right, and most eccentric ones as being at position 2 on the left or right.
At the beginning of each run, subjects fixated the center square until a colored cue (duration 300 ms) appeared at that position, with its color signaling to either execute the first 7°-saccade to the adjacent square in the indicated direction or to maintain fixation. After the eye movement was performed, the subjects' gaze remained at the new square until the next cue appeared at the newly foveated position. In fixation trials (e.g., as indicated by, for example, a red cue), the subjects' gaze remained at the same square until the following cue signaled a new eye movement. Saccade direction from all positions was randomized throughout each run, except from the ±14°-position, which necessarily required a return to the neighboring more medial position. Subjects were asked to execute each saccade as quickly and as accurately as possible and to FIguRe 1 | Trial structure and study design. (A) Subjects performed saccades from square to square along the horizontal meridian. Instructional color cues (300 ms duration) at the currently foveated square signaled the direction of the next saccade (e.g., "blue = right" vs. "green = left"; here represented by horizontal vs. vertical lines) or indicated to maintain fixation (e.g., "red = stay"). After each saccade, the subjects' gaze remained at the new square until the next saccade was triggered by a color cue. The stimulus onset asynchrony (SOA) was jittered between 1.5 and 6 s for all saccade and fixation cues. (B) Saccade vectors were oriented either toward (centripetal) or away from (centrifugal) the primary position (defined as head-centered straight gaze) and also differed in their proximity to the primary position (proximal vs. distal). The design thus resulted in four saccade trial types represented by vectors with specific start and end positions (start-to-end), i.e., proximal centrifugal (0-to-1), distal centrifugal (1-to-2), proximal centripetal (1-to-0), and distal centripetal (2-to-1), each of which could be performed to the left or right side. (C) Depending on the initial deviation of the eyes relative to primary position, rightward saccades of identical amplitude either re-center the eyes or move them further away from primary position. a random effects analysis using one-sample T-tests for voxel-wise comparisons (significance threshold p < 0.001, voxel-extent threshold k > 15). At the voxel-wise level a contrast for "all saccades vs. fixation" was calculated, which was used to establish the general saccade-related activity.
To investigate whether saccade-related regions were sensitive to the direction, the orientation, and the proximity of the saccade vector relative to the primary position, we performed a region of interest (ROI) analysis using the MarsBar analysis toolbox (Ollinger et al., 2001;Brett et al., 2002). ROIs were defined bilaterally as cubes centered at the local activity maxima derived from the conditionindependent voxel-wise comparison "all saccades vs. fixation" within FEF (8 mm × 8 mm × 8 mm) and IPS (8 mm × 8 mm × 8 mm), as well as in the SC (4 mm × 4 mm × 4 mm). From these ROIs the mean parameter estimates (beta values) for each condition of interest, i.e., left and right for proximal centrifugal (0-to-1), proximal centripetal (1-to-0), distal centrifugal (1-to-2), and distal centripetal (2-to-1), were extracted. The contrast "all saccades vs. fixation" was chosen to allow for a condition-independent selection of the areas within SC, FEF, and IPS that were associated with general saccade generation.
For all ROIs (FEF, IPS, and SC), the extracted parameter estimates were collapsed across hemispheres, yielding activity that was contralateral to saccade direction (i.e., activity within left hemisphere during saccades to the right and vice versa) and ipsilateral to saccade direction (i.e., activity within left hemisphere during saccades to the left and vice versa), respectively. Data were then analyzed via repeated-measures analyses of variance (rANOVA), with the factors contralateral predominance (contralateral vs. ipsilateral), orientation (centrifugal vs. centripetal), and center proximity (proximal vs. distal). In order to evaluate the actual region-specific shape of the BOLD response, we extracted the activity time course for each ROI across all saccade conditions based on a shape-assumption-free finite-impulse-response (FIR) model. This indicated that, in accord with earlier reports (Ikkai and Curtis, 2008), the non-standard HRF with a somewhat earlier peak was in fact also effective for characterizing activity within FEF and IPS.
Based on our previous findings regarding differential saccaderelated activity in a frontal region associated with the default-mode network of attention (Krebs et al., 2010), we performed an additional ROI analysis in bilateral medial prefrontal cortex (mPFC). The center coordinates for the cubical ROIs (4 mm × 4 mm × 4 mm; left: x, y, z = −6, 44, 10; right: x, y, z = 12, 38, −12) were derived from the contrast "fixation vs. all saccades" based on the standard HRF, yielding an unbiased ROI for assessing differential activation in the default-mode network . The extracted parameter estimates (beta values) were collapsed across left and right mPFC and analyzed via rANOVA with factors orientation (centrifugal vs. centripetal) and center proximity (proximal vs. distal).

electro-oculograM study
Although the eye-tracking signal acquired during the fMRI session can be used to quantify general task performance (e.g., the number of correct and erroneous saccades), online measures during fMRI scanning are often of limited use. Particularly, various sources of artifacts (e.g., eye-tracking disturbances caused by the onset of the fMRI gradients) lead to relatively low number of analyzable trials. Moreover, the relatively low temporal resolution of fMrI data acquIsItIon Prior to actual scanning, subjects performed a 10-min training session to get familiarized with the task. Inside the scanner, subjects performed seven experimental runs, each of 7-min duration, resulting in a total of 112 trials in each saccade condition: left and right proximal centrifugal (0-to-1), left and right proximal centripetal back to primary position (1-to-0), left and right distal centrifugal (1-to-2), left and right distal centripetal (2-to-1), and a total of 70 fixation trials (i.e., 14 trials at each of the five squares). fMRI images were acquired using a 3 Tesla GE EXCITE HD scanner with an eightchannel head-coil array. Each functional run consisted of 290 images that were acquired in an axial slice orientation (30 slices, thickness 3 mm) using an interleaved scanning order (inward spiral sequence with SENSE acceleration factor of 2, TR = 1500 ms, TE = 25 ms, FoV = 192 mm, matrix size of 64 × 64 yielding an in-plane resolution of 3 mm × 3 mm). The first five time points of each run were discarded to allow a steady magnetization to be reached.
For each subject, a T1-weighted high-resolution whole-brain anatomical scan (3D FSPGR sequence, FoV = 256 mm, yielding a voxel size of 1 mm × 1 mm × 1 mm) was acquired to enable coregistration and normalization to a template brain. During the entire experiment, subjects were instructed to keep their head-position straight and remain still, which was monitored online and facilitated by passive padding inside the head coil. In order to control for the correct execution of the saccade task, eye movements were monitored online throughout all runs using an MR-compatible eyetracking system (Viewpoint, Arrington Research, Scottsdale, AZ, USA). Eye position was monitored at 30 Hz via an infrared camera that was mounted above the subject's right eye and calibrated at the beginning of the scanning session.

fMrI data analysIs
Images were preprocessed and analyzed using the Statistical Parametric Mapping software package (SPM5; Wellcome Department of Imaging Neuroscience, University College, London, UK). Anatomical images were co-registered to the SPM template and spatially normalized. Functional images were corrected for acquisition delay, spatially realigned, and co-registered to the original T1-weighted image. After spatial normalization to a final voxel size of 2 mm × 2 mm × 2 mm, functional images were smoothed with an isotropic 4-mm full-width half-maximum Gaussian kernel. Before model estimation, a high-pass temporal filter of 128 s was applied (Ashburner and Friston, 1999).
A two-stage model was used for statistical analysis (Friston et al., 1995). In the first stage, BOLD responses were modeled by delta functions at the stimulus onsets for the nine event types (left and right centrifugal 0-to-1 saccades, left and right centripetal 1-to-0 saccades, left and right centrifugal 1-to-2 saccades, left and right centripetal 2-to-1 saccades, and fixation trials), which were then all convolved with a canonical hemodynamic response function (HRF). In order to optimize the signal estimation within the SC, as suggested by Wall et al. (2009) and Krebs et al. (2010), we used a version of the canonical HRF that peaked at 4.5 s instead of at the standard 6 s. The resulting estimates together with the corresponding temporal and dispersion derivatives as well as six realignment parameters for each run formed covariates of a general linear model (GLM, Friston et al., 1995). Individual subjects' contrast images were entered into vs. right), orientation (centrifugal vs. centripetal), and proximity (proximal vs. distal). Only the EOG data, used for deriving saccade latency data in a corresponding experiment, will be reported here. The EEG data from this corresponding experiment will be reported elsewhere.

results behavIoral results fMrI study
The horizontal eye-position trace extracted from 24 s of a single subject's fMRI session is depicted in Figure 2A and includes all eight saccade types. The eye-position data of three of the subjects could not be analyzed due to technical problems. After rejecting an average of 10% of all trials due to blinks and other artifacts (e.g., eye-tracking signal disturbances caused by the fMRI gradients) in the remaining 13 data sets, the average saccade-task performance in the included trials was well above 95% (rates of correctly executed saccades: centrifugal 0-to-1 = 0.97; centripetal 1-to-0 = 0.966; centrifugal 1-to-2 = 0.978; centripetal 2-to-1 = 0.979). However, the general noise level and the restricted temporal resolution of 30 Hz did not allow for an accurate analysis of SRT differences in the scanner. The training prior to the scanning session, the online video monitoring, and the high percentage of correct trials indicated that the task was performed accurately. For a more precise analysis, the fMRI eye-tracker (sampling rate of 30 Hz) is not sufficient to provide a sensitive measure for small differences in SRTs between the conditions of interest. To compensate for these limitations, we derived saccade onset times by analyzing the electro-oculogram (EOG) channels recorded during an equivalent electroencephalographic (EEG) study that was performed outside the scanner in a sparsely lit electrically shielded chamber. Here, 16 subjects (mean age ± SD: 22 ± 2.6, six female) performed the exact same task used in the fMRI study, with the only differences being a higher number of trials per condition (160) and shorter SOAs (randomly jittered between 1 and 1.5 s). The EOG channels were continuously recorded with a band-pass filter of 0.01-100 Hz at a sampling rate of 500 Hz (SynAmps amplifiers Neuroscan; http:// www.neuroscan.com).
After saccade-locked segmentation of the horizontal EOG channel, trials containing artifacts (e.g., blinks) were removed from the data while clean trials were separated in correctly executed and erroneous saccades using the Brain Vision Analyzer 1.05 software (Brainproducts Co, Munich, Germany). Saccade onset latencies of correctly executed saccades were averaged for each experimental condition (left and right centrifugal 0-to-1, left and right centripetal 1-to-0, left and right centrifugal 1-to-2, left and right centripetal 2-to-1) and analyzed via rANOVA with the factors direction (left

FIguRe 2 | eye-position traces and saccadic reaction times (SRTs). (A)
A single subject's horizontal eye-position trace from the fMRI scanning session is depicted for a sequence of 24 s (16 TRs) including saccades to all five possible positions at 0°, ±7°, and ±14° of visual angle. According saccade direction is schematically depicted below, i.e., R = rightward and L = leftward. (B) The same trial sequence is shown as a single subject's horizontal EOG channel depicting rightward and leftward saccades, respectively (upward vs. downward trace shifts). (C) SRTs extracted from the EOG data and averaged across leftward and rightward saccades for each condition (arrow legend exemplary depicts rightward shifts). SRTs were shorter for centripetal as compared to centrifugal saccades. Moreover, the SRT advantage for centripetal saccades increased with greater distance to the primary position. Error bars depict the standard error (SE) of the mean across subjects; asterisks indicate significant T-tests.

sIgnal tIMe course wIthIn selected roIs
In order to estimate the actual shape of the BOLD signal, the time course was extracted from the selected ROIs based on a shape-assumption-free FIR model. The extracted values averaged across all saccade conditions and hemispheres for the selected ROIs are displayed in Figure 3B. The BOLD response within all three areas extracted from this shape-assumption-free FIR model peaked at 4.5 s, in accordance with our use of the non-standard canonical HRF.

roI-based analysIs of paraMeter estIMates
In order to investigate effects of the three main factors (contralateral predominance, orientation, and proximity) in an orthogonal fashion, parameter estimates were extracted from the conditionindependent-selected ROIs (i.e., all saccade trials vs. fixation trials; Figure 3A) and analyzed via rANOVA. The analyses confirmed the contralateral predominance across all saccade conditions within the selected bilateral-pair ROIs (bilateral SC, FEF, and IPS; Figure 3C). The activity enhancement in the hemisphere contralateral to saccade direction was most robust within the IPS (main effect of contralateral predominance: F (1,15) = 56.8, p < 0.0001), followed by the SC [F (1,15) = 15.7, p = 0.001] and the FEF [F (1,15) = 5.1, p = 0.039].
With regard to the saccade orientation and proximity relative to the primary position, we observed a main effect of orientation within the SC [F (1,15) = 7.3, p = 0.016], confirming that centrifugal saccades were associated with enhanced collicular activity as compared to centripetal saccades (Figure 4). In the cortical ROIs, however, no main effects of saccade orientation or proximity to primary position were observed. However, the contralateral enhancement in IPS was significantly pronounced during centrifugal saccades [interaction of contralateral predominance × orientation: F (1,15) = 7.4, p = 0.016]. In addition, IPS displayed a tendency for a greater contralateral enhancement with increasing distance to the primary position, as indicated by a trend in the interaction term [contralateral predominance × proximity: F (1,15) = 4.19, p = 0.059]. An analogous analysis within the FEF revealed a significant interaction between orientation and proximity [F (1,15) = 5.9, p = 0.028]. This interaction reflected an inversion of activity differences between centripetal and centrifugal saccades dependent on the proximity to the primary position (i.e., proximal: centripetal > centrifugal; distal: centripetal < centrifugal). No additional main or interaction effects were observed in the ROI analyses (all p-values >0.1).
Post hoc specific comparisons using two-tailed paired T-tests were performed focusing on activity within the predominant SC (i.e., contralateral to saccade direction) to further investigate the observed differences due to orientation and proximity within the respective predominant hemisphere. Activity in the predominant SC was significantly higher during distal centrifugal saccades (i.e., 1-to-2) than during both proximal (1-to-0) and distal (2-to-1) centripetal saccades (proximal T (15) = 3.08, p = 0.008; distal T (15) = 2.24, p = 0.04). Furthermore, centripetal saccades from the most eccentric position were associated with diminished SC activity as compared to centrifugal saccades from primary position [T (15) = 3.99, p = 0.001].
Regarding activity in the predominant hemispheres of FEF and IPS, the only significant difference observed was between distal centrifugal (1-to-2) and distal centripetal saccades (2-to-1), with the particularly with respect to SRTs, the saccadic task performance in the equivalent EOG study, with much higher temporal resolution capability, was analyzed. behavIoral results eog study Figure 2B depicts a single subject's EOG trace of the identical trial sequence as the fMRI study (albeit with shorter SOAs). On average, no more than 6% of the trials in each condition had to be excluded due to artifacts in the EOG (e.g., blinks). Regarding the remaining trials, subjects consistently performed saccades in the cued direction as indicated by a mean percentage of correct saccades well above 95% (rates of correctly executed saccades: centrifugal 0-to-1 = 0.993; centripetal 1-to-0 = 0.994; centrifugal 1-to-2 = 0.989; centripetal 2-to-1 = 0.999). Thus, the least number of errors was made during centripetal saccades starting from the most eccentric square, presumably reflecting the predictability of saccade direction at this position.
In regards to the SRTs, rightward saccades tended to be slightly faster than leftward ones, although this difference did not quite reach significance [direction F (1,15) = 4.2, p = 0.058]. Regardless, in that there were also no significant interactions of any of the factors with left-vs. right saccade direction, the SRTs were collapsed across leftward and rightward saccades. Statistical analysis of the collapsed SRTs revealed shorter latencies for centripetal (proximal 413 ms; distal 382 ms) as compared to centrifugal (proximal 419 ms; distal 421 ms) saccades as indicated by a main effect of orientation [F (1,15) = 47.03, p < 0.0001]. Moreover, there was a main effect of proximity [F (1,15) = 15.22, p = 0.001] which was accompanied by a significant interaction of orientation by proximity [F (1,15) = 28, p < 0.0001], resulting from the SRT advantage for centripetal saccades increasing with greater distance from the primary position ( Figure 2C). Post hoc T-tests revealed that centripetal saccades from the most eccentric position yielded faster SRTs as compared to all other saccades (all p-values <0.0001), again most likely reflecting the high directional predictability. Furthermore, centripetal saccades toward primary position (re-centering saccades) yielded shorter SRTs as compared to centrifugal saccades (both proximal and distal conditions: p < 0.01). No significant SRT difference was observed between proximal and distal centrifugal saccades (p = 0.379). Figure 3A displays representative slices of the comparison between all saccade trials and fixation trials, a contrast that subtracts out the sensory processing of the stimulus change at the foveated location. The voxel-wise contrast revealed a network of regions involved in the execution of endogenously cued saccades that were located along the horizontal meridian, including large portions of bilateral primary visual cortex along the calcarine fissure, lateral occipital cortex, FEF, the supplementary eye field (SEF), IPS, and SC (a detailed description of activation clusters is provided in Table 1; T-values refer to the local activity maxima). Note that the displayed activation clusters are all based on the modified HRF model peaking at 4.5 s, which was more effective for the BOLD signal within the SC (see also Wall et al., 2009). For comparison, an analogous voxel-wise analysis based on the standard HRF peaking at 6 s can be found in Table 2.

voxel-wIse fMrI actIvatIons durIng saccade executIon
It is important to note that the high directional predictability of distal centripetal saccades (2-to-1) imposes limitations on their use for direct comparisons with the remaining saccade types. In contrast, proximal centripetal (1-to-0) and distal centrifugal (1-to-2) saccades were closely matched regarding shift amplitude, directional probability, and the starting gaze-deviation from primary position, thus yielding the "purest" comparison regarding the investigation of the re-centering bias.
In summary, activity within the SC, as well as in FEF and IPS, was enhanced contralateral to saccade direction (contralateral predominance, see Figure 3C). When focusing on the predominant hemisphere, activity within the SC was diminished during centripetal as compared to centrifugal saccades (Figure 4). This was especially true for saccades starting from the most eccentric position, but the effect was still present for the more proximal former eliciting higher activity in the predominant FEF [T (15) = 2.68, p = 0.017]. However, there was a tendency for higher activity within the predominant FEF for proximal centrifugal saccades starting from primary position (0-to-1) as compared to distal centripetal saccades [2-to-1; T (15) = 2, p = 0.064]. No further post hoc T-tests reached significance (all p-values >0.1). An additional ROI analysis within the mPFC bilaterally that was conducted to investigate differential default-mode network activity during the different saccade types revealed a significant effect of orientation [F (1,15) = 10.49, p = 0.006], with increased deactivation levels for centrifugal as compared to centripetal saccades. Post hoc specific T-tests verified that this main effect was mainly driven by significant differences between distal centrifugal (1-to-2) and proximal centripetal (1-to-0) saccades [T (15) = 2.39, p = 0.030] and between proximal centrifugal (0-to-1) and distal centripetal (2-to-1) saccades [T (15) = 3.37, p = 0.004].

FIguRe 3 | Saccade-related BOLD activity. (A)
Saccade trials, as compared to fixation trials, were associated with robust activity in a broad network that included bilateral IPS, FEF, SEF, SC, as well as large parts of the visual cortex. This contrast was based on an HRF that has been optimized for subcortical activity, with a peak at 4.5 s. Activity is displayed on a standard T1-weighted image (display cut-off T = 3; extent threshold k = 15). (B) Mean BOLD time courses for the selected ROIs, namely the SC (depicted in red), FEF (blue), and IPS (yellow) are displayed averaged across all saccade trials and both hemispheres. (C) The collapsed ROI-based parameter estimates (beta) for leftward and rightward saccades with respect to the hemisphere contra-and ipsilateral to saccade direction revealed a relative activity enhancement in the contralateral (gray bars) as compared to the ipsilateral (white bars) hemisphere within SC, FEF, and IPS.
In line with the notion of a facilitation for saccades toward the primary position (re-centering bias), we observed shorter SRTs for centripetal saccades, as has previously been reported in both animals (Albano and Wurtz, 1982;McIlwain, 1986;Pare and Munoz, 1996) and humans (Laurutis and Robinson, 1986;Fuller, 1996). This latency bias has been initially associated with the general dynamics of eye-head coordination in natural head-free gaze shifts (combined eye-and head-movement), but appears to be also conserved in a headrestrained setting. Importantly, the re-centering bias was observed not only for the saccades starting from the most eccentric positions, which differ in their directional predictability from all other saccades, but also when comparing centripetal and centrifugal saccades starting at position "1" (1-to-0 vs. 1-to-2; see Figure 1C). This latter comparison is considerably more compelling since these saccade types differ only in their centripetal (i.e., re-centering) vs. centrifugal orientation relative to the primary position while being matched regarding both starting position and directional predictability.

saccade-related actIvIty and contralateralIty
The comparison of all saccades to fixation trials provided a general activity map for saccade generation, including the primary areas of interest in this study -namely, the SC, IPS, and FEF. We also observed activity clusters in various regions of visual cortex, as well as in some additional subcortical regions related to stimulus and saccade processing (including medial thalamus, lateral geniculate nucleus, pulvinar, and putamen). Although several of these areas are implicated in saccade processing (e.g., primary visual cortex, see Bodis-Wollner et al., 1997;Sylvester et al., 2005;Rieger et al., 2008), we focused on the two main cortical areas (i.e., FEF and IPS) and the key subcortical area (i.e., SC) that are central to the generation of saccadic eye movements.
ones. Importantly, this re-centering activity modulation in the SC closely matched the pattern of the corresponding SRTs (Figure 2C). At the cortical level, the only robust orientation effect within paired comparisons was the increased activity in the predominant FEF for distal centrifugal saccades (1-to-2) as compared to centripetal saccades starting at the most eccentric position (2-to-1).

dIscussIon behavIoral perforMance
The general performance level in the saccade task was well above 95% in both the fMRI study and the analogous EOG one. Based on the highly comparable performance, we based the SRT analysis on the EOG data, due to its considerably higher temporal resolution, as well as the higher trial numbers and superior data quality. Mean SRTs varied between 380 and 420 ms depending on the condition, consistent with the normal SRT range of centrally cued saccades (Walker et al., 2000).  lie in the visual field ipsilateral to the saccade direction, and thus a robust contribution of the visual input should tend to produce a relatively increased contralateral predominance for centripetal saccades close to the primary position. None of these asymmetries was observed in the present data. Furthermore, a recent study comparing a covert and overt version of a spatial attention task found highly similar activity levels in the cortical eye fields, despite the fact that only overt trials entailed a change in eye position and thus a change of the current visual input (Ikkai and Curtis, 2008).

dIfferentIal actIvIty durIng re-centerIng saccades
In keeping with the behavioral analysis, the best controlled comparison for investigating the underlying neural pattern of the recentering bias is the comparison of centripetal and centrifugal saccades starting at the same orbital position (1-to-0 vs. 1-to-2). Our results indeed show that activity within the predominant SC side (contralateral to saccade direction) was significantly diminished during centripetal re-centering saccades as compared to centrifugal saccades from the equivalent starting position (Figure 4), a pattern that was paralleled by the SRT differences ( Figure 2C). Our observation of reduced evoked SC activity in humans during re-centering as compared to centrifugal saccades is consistent with the animals studies described above, and might be related to dynamic gain field modulations in collicular movement neurons The ROI-based analysis of the BOLD activity revealed that activity in the SC was significantly enhanced contralateral to the direction of the saccade target, a pattern that has been demonstrated in numerous studies in animals (e.g., Goldberg and Wurtz, 1972;Field et al., 2008) and recently in humans (Krebs et al., 2010). A contralateral bias was also observed in the cortical eye fields (FEF and IPS), matching the results of previous fMRI studies using memory-guided saccade paradigms (Schluppeck et al., 2005;Curtis and Connolly, 2008). Importantly, the contralateral predominance appeared to be independent of the additional experimental factors, i.e., orientation and proximity relative to the primary position, except for a relative increase in contralateral IPS activity during centrifugal as compared to centripetal saccades, which might be due to overlapping topographic representations of left and right IPS (Kastner et al., 2007;Konen and Kastner, 2008).
It is important to consider whether the observed contralateral activity enhancement in saccade-related regions could have been simply induced by differential sensory input. First, configurational aspects argue against this possibility. For centrifugal saccades, the majority of both the pre-and post-saccadic visual input from the stimulus array would lie in the visual field contralateral to the saccade direction, which would thus tend to lead to enhanced activity in the ipsilateral hemisphere. For centripetal saccades, on the other hand, the majority of pre-and post-saccadic visual input would FIguRe 4 | Condition-specific BOLD response. ROI-based parameter estimates are collapsed for leftward and rightward saccades with respect to the hemisphere contralateral (gray bars) and ipsilateral (white bars) to saccade direction. Activity in the predominant SC (i.e., contralateral) was modulated by the orientation of the saccade relative to primary position. In particular, the generation of re-centering saccades (1-to-0) was associated with significantly decreased activity as compared to centrifugal saccades (1-to-2) of identical amplitude, directional probability, and equivalent starting position (arrow legend exemplary depicts rightward shifts). Activity in the contralateral FEF and IPS was not significantly modulated by the saccade orientation relative to primary position. Error bars depict the standard error (SE) of the activity means across subjects; asterisks indicate significant T-tests focusing on contralateral activity only. the current paradigm required the voluntary initiation of saccades following central instructional cues, we expected a similar degree of engagement of these regions irrespective of saccade orientation. The lack of significant cortical activity changes during re-centering saccades is thus consistent with the view that the re-centering bias is a fairly automatized process subsequent to cue interpretation and target selection and is thus mainly implemented at the subcortical level, where current eye-and head-positions as well as the future gaze position are integrated (e.g., Campos et al., 2006). It remains possible, however, that potential cortical contributions to the re-centering bias may have been diluted by neural processes that are common to all saccade types, i.e., cue interpretation and target selection, which can not be easily separated due to the low temporal resolution of the BOLD signal. Such a possibility might explain why there appeared to be a trend in both FEF and IPS for reduced activity during 1-to-0 as compared to 1-to-2 saccades (see Figure 4) which did not reach significance. Along similar lines, it is furthermore possible that this trend might have arisen from a corollary signal from SC to these cortical regions Wurtz, 2004, 2008) conveying a copy of the re-centering-related motor signal.
While no significant activity modulation was found in cortical areas that are directly implicated in saccade processing, the present data still provide some evidence for subcortical-cortical interactions. In particular, and consistent with our previous work (Krebs et al., 2010) the mPFC, a region associated with the defaultmode network Raichle et al., 2001), exhibited differential activity during centrifugal and centripetal saccades with considerably more pronounced deactivation for the former. Based on the notion that less demanding tasks are associated with reduced deactivations in this area, this pattern might be related to the diminished neural activity that we observed in the SC during the execution of re-centering saccades.

conclusIons
The current study investigated the neural basis of the re-centering bias in humans that has been demonstrated behaviorally. We systematically varied saccade direction and orbital starting position while keeping saccade amplitude and directional predictability constant across the conditions of interest. The observed reduced BOLD response in the SC during re-centering saccades is in line with the notion that this subcortical brain region integrates information about future gaze positions as well as eye-and head-position signals, thus engendering effective exploration of the visual environment.

acknowledgMents
We thank Ken Roberts and Michael Scholz for technical support. This work was funded by NIH grants R01-MH060415 and R01-NS051048 to Marty G. Woldorff. (Van Opstal et al., 1995;Salinas and Thier, 2000). In particular, the firing activity of a large portion of collicular neurons are modulated by the orbital position (Van Opstal et al., 1995;Pare and Munoz, 2001;Campos et al., 2006), and the SC thus appears to be in the position to effectively integrate eye-and head-position information (Sparks, 2002;Kardamakis and Moschovakis, 2009). It has been proposed that for eye positions that deviate from the primary position, stronger muscular elastic forces need to be overcome to initiate a saccade to an even more centrifugal position as compared to a centripetal saccade. For example, the initiation of rightward saccades with identical amplitude would require differential levels of activity in the left SC (contralateral to saccade direction) depending on the starting point relative to the primary position (see Figure 1C). Alternatively, or in addition, there might be an involvement of specific neurons in the posterior SC that have been shown to re-center the eyes independently of the initial position (Guitton et al., 1980), a mechanism that could result in less requirements for movement neurons that code a specific saccade vector. A recent study has shown evidence for a preferential coding of the primary position in early visual cortex (V1) of monkeys (Durand et al., 2010). Therefore it seems likely that the re-centering bias relies on projections between visual cortex and the SC (Fries, 1984), which is consistent with a recent study reporting egocentric referenceframe properties of human visual cortex using retinotopic mapping (Rieger et al., 2008).
While the SC was clearly sensitive to the orientation of the saccade vector relative to the primary position, neither of the cortical eye fields (i.e., FEF and IPS) was systematically influenced by this factor. Although the IPS is generally involved in the continuous updating of visual space by transforming coordinates of the eye-and head-centered reference frames (Duhamel et al., 1992;Colby et al., 1995;Medendorp et al., 2003), the present data reveals no re-centering-specific activity modulation in this region, which is consistent with earlier findings (Rieger et al., 2008). Regarding the FEF, it has been reported that, under some conditions, different orbital starting positions change the resulting saccade vector elicited by micro-stimulation in humans (Blanke and Seeck, 2003). However, these modulations were restricted to differences in starting positions along the vertical meridian while variable starting positions along the horizontal meridian did not influence direction or amplitude of the elicited saccades.
It is well established that activity in the cortical eye fields is generally increased during tasks requiring higher levels of top-down processing both in overt and covert attention paradigms in humans (reviewed in Corbetta and Shulman, 2002). More specifically, activity in both IPS and FEF is increased during the execution of endogenously cued (voluntary) saccades, as compared to exogenously cued (reflexive) saccades, an effect that has been interpreted as likely reflecting the higher processing demands during cue interpretation and saccade preparation for the former (Mort et al., 2003). Since