Analysis of Primate IBN Spike Trains Using System Identification Techniques. II. Relationship to Gaze, Eye, and Head Movement Dynamics During Head-Free Gaze Shifts
Kathleen E. Cullen and
Daniel Guitton
Aerospace Medical Research Unit and the Montreal Neurological Institute, McGill University, Montreal, Quebec H3G 1Y6, Canada
 |
ABSTRACT |
Cullen, Kathleen E. and Daniel Guitton. Analysis of primate IBN spike trains using system identification techniques. II. Relationship to gaze, eye, and head movement dynamics during head-free gaze shifts. J. Neurophysiol. 78: 3283-3306, 1997. We have investigated the relationships among the firing frequency B(t) of inhibitory burst neurons (IBNs) and the metrics and dynamics of the eye, head, and gaze (eye + head) movements generated during voluntary combined eye-head gaze shifts in monkey. The same IBNs were characterized during head-fixed saccades in our first of three companion papers. In head-free gaze shifts, the number of spikes (NOS) in a burst was, for 82% of the neurons, better correlated with gaze amplitude than with the amplitude of either the eye or head components of the gaze shift. A multiple regression analysis confirmed that NOS was well correlated to the sum of head and eye amplitudes during head-free gaze shifts. Furthermore, the mean slope of the relationship between NOS and gaze amplitude was significantly less for head-free gaze shifts than for head-fixed saccades. NOS is a global parameter. To refine we used system identification techniques to evaluate a series of dynamic models in which IBN spike trains were related to gaze or eye movements. We found that gaze- and eye-based models predicted the discharges of IBNs equally well. However, the bias values required by gaze-based models were comparable to those required in our head-fixed models whereas those required by eye-based models were significantly larger. The difference in biases between gaze- and eye-based models was very strongly correlated to the mean head velocity (
) during gaze shifts [R =
0.93 ± 0.15 (SD)]. This result suggested that the increased bias required by the eye-based models reflected an unmodeled
input onto these cells. To pursue this argument further we investigated a series of dynamic models that included both eye velocity (
) and
terms and this confirmed the importance of these two terms. As in our head-fixed analysis of companion paper I, the most valuable model formulation also included an eye saccade amplitude term (
E) and was given by B(t) = r0 + r1
E + b1
+ g1
where r0, r1, b1, and g1 are constants. The amplitude of the head velocity coefficient was significantly less than that of the eye velocity coefficient. Furthermore, in our population long-lead IBNs tended to have a smaller head velocity coefficients than short-lead IBNs. We conclude that during head-free gaze shifts, the head velocity signal carried to the abducens nucleus by primate excitatory burst neurons (EBNs; if EBNs and IBNs carry similar signals) must be offset by other premotor cells.
 |
INTRODUCTION |
Inhibitory burst neurons (IBNs) in the paramedian pontine reticular formation, caudal to the abducens motor nucleus, project to the contralateral abducens nucleus (ABD) and inhibit its activity during saccades (Hikosaka and Kawakami 1977
; Hikosaka et al. 1978
, 1980
; Scudder et al. 1988
; Strassman et al. 1986
; Yoshida et al. 1982
). In the preceding companion paper (Cullen and Guitton 1997a
) we analyzed, by using system identification techniques, the discharge characteristics of IBNs in the head-fixed primate. In this paper we are concerned with the discharge of these cells in the primate whose head is free to move. We define gaze = eye-in-space = eye-in-head + head-in-space.
IBNs are the inhibitory component of the saccade burst generator; the excitatory drive onto motoneurons (MNs) being provided by excitatory burst neurons (EBNs). EBNs project to IBNs and both cell types are thought to carry the same signals (see reviews by Fuchs et al. 1985
; Hepp et al. 1989
). Burst neurons (BNs) have been further subdivided into short- and long-lead bursters depending on the magnitude of the time interval by which the discharge leads movement. IBNs are thought to be gated by inhibitory projections from omnipause neurons (OPNs) (Evinger et al. 1982
; Furuya and Markham 1982
; Keller 1974
; King and Fuchs 1977
; Nakao et al. 1980
). Investigations in the cat demonstrated that OPNs pause for the entire duration of head-free gaze shifts, rather than for the duration of the saccadic eye movement's contribution (Paré and Guitton 1990
). The latter duration is usually shorter, especially in large gaze shifts that are bigger than the oculomotor range and for which the eye frequently remains immobile, at an eccentric orbital position, while gaze motion is carried by the head. There is also evidence that the pause in activity of primate OPNs may also be gaze related (Phillips 1993
). Because the OPN pause disinhibits IBNs, it follows that IBNs should discharge during the entire duration of gaze shifts in both cat and monkey. We have seen this in the head-free cat; the duration of an IBN's burst was better correlated with gaze than with saccade duration (Cullen et al. 1993
).
The first studies of BNs in the head-free primate classified short-lead BNs as "gaze related" or "saccade related" (Whittington et al. 1984
). It was not known whether the neurons were IBNs or EBNs. Gaze-related BNs were defined as cells for which the total number of spikes (NOS) in a burst is the same for a given displacement of the visual axis irrespective of whether the head is fixed or free. Saccade-related BNs were defined as cells for which the total NOS in a burst is the same for eye movements of equal amplitude irrespective of whether the head is fixed or free. Gaze-related bursters were thought to be converted to saccade-related bursters by subtracting from the former a signal encoding head velocity (
). Our studies regarding the cat (Cullen et al. 1993
) demonstrated that the discharge of cat IBNs during combined eye-head gaze shifts differed from that of the prototypical saccade-related and gaze-related BNs. In the head-free condition, the total NOS in a burst was well correlated with gaze amplitude and generally poorly correlated with saccade amplitude. This is a property of gaze-related neurons, not saccade-related ones. However, in the head-fixed condition the slope of the relationship between the NOS in a burst and eye saccade amplitude was significantly higher than that between NOS and gaze amplitude measured during head-free gaze shifts. This is not a property of gaze-related BNs. In summary, IBNs in cats were neither the saccade- nor gaze-related neurons defined by Whittington et al. (1984)
; they were related to the movement of the visual axis in different ways depending on whether the cat's head was fixed or free. Put another way, the discharge characteristics of cat IBNs could not be extrapolated from head-fixed data to predict accurately either eye or gaze saccade amplitude in the head-free condition.
Recently a number of studies in primates have suggested that the nature of the signals carried by their IBNs may also be more complex than initially proposed. For example, Phillips et al. (1995)
proposed what they called command, head, eye, and gaze BNs. By comparison, Tomlinson and Brance (1991)
recorded saccade-related bursters in the IBN area of the head-free primate but surprisingly, these neurons also discharged throughout the duration of large gaze shifts when the eye was relatively immobile at an eccentric position in the orbit. The results of these analyses have not permitted a clear quantitative determination of whether BN discharges encode gaze, eye, or head. To address this issue we carried out dynamic analyses of IBN discharges in the totally head-free macaque monkey, using the objective system identification algorithms described in Cullen et al. (1996)
, and applied it to IBNs in the head-fixed primate as described in companion paper I (Cullen and Guitton 1997a
).
Figure 1 shows an extension of the local feedback system for controlling head-fixed saccades (described in companion paper I Cullen and Guitton 1997a
), to the control of gaze shifts in the head-free condition. As in companion paper I Fig. 1 demonstrates that the BN signal can be considered from either the upstream or downstream perspectives. This paper deals with the latter; the subsequent companion paper III (Cullen and Guitton 1997b
) will consider the former. The schematic proposes that gaze shifts are controlled by a signal encoding gaze motor error (
g) produced by subtracting the angle the eye (
E*) and head (
H*) have rotated from the desired angular rotation (
T); i.e.,
g =
T
E*
H*. There is considerable agreement on this conceptual model (Fuller et al. 1983
; Guitton and Volle 1987
; Guitton et al. 1984
, 1990
; Laurutis and Robinson 1986
; Pélisson and Prablanc 1986
; Pélisson et al. 1988
, 1989
, 1995
; Roucoux et al. 1980
; Tomlinson and Bahra 1986a
,b
; Tomlinson 1990
; see review by Guitton 1992
). In our simplified schema illustrated, the BN pool (B) is shown in an ambiguous state, driven by
g but itself only driving eye MNs. Furthermore, neck MNs are shown as being driven by an unknown source. We will show that the organization in Fig. 1 is inadequate. This has already been demonstrated for the cat where a modified version of this model was proposed by Galiana predicts movement trajectories and burst generator firing frequency profiles (Cullen et al. 1993
; Galiana and Guitton 1992
; Guitton et al. 1990
). This model assumed that the cat burst generator is strongly coupled to both eye and head motor circuits (as suggested by the results of Grantyn and Berthoz 1987
and Guitton et al. 1990
), but because this is uncertain in the monkey we have not included a link between B and neck MNs in Fig. 1. Our present analysis of monkey IBNs reveals strong analogies with our previous findings in the cat (Cullen et al. 1993
).

View larger version (16K):
[in this window]
[in a new window]
| FIG. 1.
Classic local feedback model for saccade generation (shown in Fig. 2 of companion paper I, Cullen and Guitton 1997a ) extended to describe the control of head-free gaze shifts. E*, angular rotation of the eye since onset of gaze shift; obtained by integrating the burst neuron (BN) output here assumed (incorrectly as our analysis will show) to encode eye velocity. H*, angular rotation of the head in space (or relative to body) since onset of gaze shift; obtained by integrating the canal input (or from propioceptive feedback from the head plant). T, desired angular rotation of visual axis in space (or relative to body). Difference between ( E* + H*) and T provides an estimate of gaze motor error [ g(t)]. E(t) and H(t), actual eye and head movements generated as a result of motoneuron (MN) signals passing through plant dynamics of the eye and head plants, respectively. BNs are thought to be gated by inhibitory projections from omnipause neurons (OPNs). As discussed in companion paper I (Cullen and Guitton 1997a ), BNs can be analyzed from upstream and downstream perspectives. In this study we focused exclusively on predicting inhibitory BN (IBN) firing frequency based on downstream signals (indicated by the dotted box). We have investigated a series of models in which the spike train dynamics of the IBNs were described in terms of eye, head, and gaze (eye + head) trajectories. B, BNs; , neural integrator; ?, unknown origin.
|
|
An important conclusion of our analysis is that IBNs in the head-free primate carry both eye velocity (
) and head velocity (
) signals to the ABD and that vestibular projections to this motor nucleus are necessary to offset the
signal. A further conclusion of our analysis is that head-fixed models are not adequate for predicting the discharges of primate IBNs during head-free gaze shifts.
 |
METHODS |
The neurons described in this study were obtained from the same two monkeys (Macaca fascicularis) whose IBNs were studied in the head-fixed condition and described in companion paper I (Cullen and Guitton 1997a
). The surgical preparation of the animals and the methods used for obtaining extracellular recordings from IBNs were identical to those previously described in that paper. The same cells are analyzed in the present study, and as described in companion paper I, these neurons were categorized as IBNs on the basis of their physiological responses during head-fixed saccades, vestibular nystagmus, smooth pursuit, and their location in the IBN area.
Data collection and analysis
After we collected data from a neuron during saccadic eye movements made by a monkey with its head restrained, we carefully released the animal's head and attempted to continue recording the activity of the same cell during active orienting gaze shifts made with the head completely unrestrained. If the neuron was lost at this stage in the experiment (which occurred regularly), we again restrained the head and attempted to locate another IBN in the head-fixed condition. For neurons that remained well isolated during the head-free recording period, data were recorded while the monkey oriented to 1) a target light that was stepped between positions ±10, 20, and 35° relative to the straight ahead position and 2) a fruit target that appeared unexpectedly on either side of an opaque screen placed in front of the monkey. Gaze and head positions were recorded by using the magnetic search coil technique and stored on digital audio tape (DAT) tape with the recorded unit activity sampled at 20 kHz. We corrected for the nonlinearity in recorded eye movements, which was inherent to the measurement system, during off-line analysis. In head-free recordings the action potential amplitude frequently varied during a movement, probably because of slight movement of the electrode tip relative to the cell. Hence, it was important to review the data off-line to carefully monitor the shape of the unit spike wave form to guarantee an uncontaminated single unit recording. Off-line, the level of a windowing circuit was set by visual inspection at an appropriate level to generate a pulse coincident with the rising phase of the action potential. This pulse was sent to the event channel of the intelligent CED1401 (Cambridge Electronic Design) peripheral device, with the result that the event times of the spikes were logged by the computer along with the position signals, low-pass filtered (250 Hz, 8 pole Bessel), and sampled at 1,000 Hz.
During off-line analysis using Matlab, the sampled gaze- and head-movement traces were digitally filtered at 125 Hz (Cullen et al. 1996
). A spike density function, in which a Gaussian function with a width of 5 ms was convolved with the spike train, was used to represent the neural discharge. This procedure assured that the frequency content of the movement and spike density traces were similar (see Cullen et al. 1996
). Eye position was calculated from the difference between the recorded gaze and head position signals. Gaze, head, and eye velocities were derived digitally from position data. We considered only those gaze shifts in which the vertical position component was less than one-third of the entire amplitude; this assured the analysis of only nearly horizontal gaze movements. The onset and offset of the ocular saccade and head movement components of the gaze shift and the gaze shift itself were defined by using a 20°/s eye velocity criterion. In cases where the head continued to move slowly, long after the end of the gaze shift (>100 ms; in all cases the neuron had ceased firing), the head movement offset was set to the time 100 ms after the termination of the gaze shift. Burst duration was defined as the time between the onset and offset of the burst. Burst onset was defined as when the first spike of the burst occurred. Because these neurons frequently generated a few additional spikes after a saccade or gaze shift, burst offset was defined as when 95% of the spikes in the burst had occurred.
Models for BN firing rate
The goal of this analysis was to determine whether IBN discharges were best linked to gaze, eye, or head trajectories or a combination thereof. We employed system identification techniques (Cullen et al. 1996
) to objectively analyze different models that relate the dynamics of spike trains to quantities describing gaze, eye, or head trajectories (e.g., velocity). The specific models are considered in RESULTS. The methods for parameter estimation were the same as those used in companion paper I (Cullen and Guitton 1997a
). The gaze, eye, and head movements that occurred between the onset and offset of the gaze shift were used in our analysis. Model fits were carried out over the entire duration of the gaze shift, because IBNs continued to discharge throughout the gaze shift, even during the end portion where the eye had stopped moving and/or had reversed direction (Figs. 2 and 3). For each model, optimal fits were made to an ensemble of ~40 gaze shifts of different amplitudes between 10 and 70°. Models were ranked by noting whether a model provided an increase in the variance accounted for (VAF) as well as a simultaneous decrease in a cost index (the Bayesian information criterion, BIC) whose value will decrease for an increasingly complex model only when the addition of parameters is warranted.

View larger version (22K):
[in this window]
[in a new window]
| FIG. 2.
Behavior of typical short-lead IBN (SLIBN; L0702) during 3 combined eye-head gaze shifts in the ipsilateral direction. Gaze position (G) = visual-axis-in-space = eye-in-head (E) + head-in-space (H). Note that the neuron continued firing throughout entire duration of gaze shift, even during periods in which eye (E) was stationary in orbit and/or had reversed direction. Dotted vertical lines, onset and offset of gaze shift; SD, spike density; , , and , gaze, eye and head velocity, respectively.
|
|

View larger version (26K):
[in this window]
[in a new window]
| FIG. 3.
Behavior of typical long-lead IBN (LLIBN; H0409) during three combined eye-head gaze shifts in ipsilateral direction. LLIBNs, like SLIBNs, continued to discharge throughout entire duration of a gaze shift. Abbreviations are the same as those used in Fig. 2.
|
|
 |
RESULTS |
In this study we consider the same population of 28 IBNs whose activity during saccadic eye movements made by the head-fixed monkey was presented in companion paper I (Cullen and Guitton 1997a
). Here we specifically investigate the nature of the relationship among the discharge of these cells and the eye, head, and gaze movements generated during voluntary combined eye-head gaze shifts.
General discharge characteristics
We retained the same classification of short-lead IBNs (SLIBNs; n = 16) and long-lead IBNs (LLIBNs; n = 12) used in companion paper I (Cullen and Guitton 1997a
), so as to be compatible with our analysis in that paper and with Scudder et al. (1988)
. Recall that the mean period between the onset of the first spike and the onset of eye velocity during head-fixed saccades was taken as
15 ms for SLIBNs and >15 ms for LLIBNs. In this paper we will use the same example neurons as in companion paper I.
Figure 2 demonstrates the discharge of our typical SLIBN (L0702) that accompanied gaze shifts made with coordinated eye and head movements by the head-free monkey. As determined in companion paper I (Cullen and Guitton 1997a
), during head-fixed saccades the first spike in this neuron's discharge led ipsilateral saccadic eye movements on average by 15 ms; the duration of its burst was tightly correlated with saccade duration (slope = 1.0; correlation coefficient R = 0.94); and the total NOS generated during a burst was well correlated with the amplitude of the horizontal component of eye saccades (slope = 0.95; R = 0.79). During head-free gaze shifts this cell's discharge led ipsilateral gaze movements on average by 14 ms, and the duration of the burst was better correlated with the duration of the entire gaze shift than with the duration of the ocular saccade(R = 0.97 vs. 0.84). The strong correlation between burst and gaze duration is particularly evident in the middle column of Fig. 2. The cell continued to fire until the gaze shift had ended even though the eye movement reached a position limit (plateau) where the eye velocity was near zero (i.e., the rapid ocular saccade had ended) and the final portion of the gaze shift was generated by the head movement alone. In other gaze shifts this cell continued firing until the end of the gaze shift, despite the fact that the eye movement actually reversed its direction before gaze stabilized.
Figure 3 illustrates the behavior of our typical LLIBN (H0409) during three example head-free gaze shifts. In the head-fixed condition, this IBN generated a burst of firing whose first spike led ipsilateral saccadic eye movements by 51 ms; the total NOS generated during a burst was proportional to the horizontal component of the saccade (slope = 1.0; R = 0.85) and the duration of its burst was correlated with the duration of the saccade (slope = 0.88; R = 0.57). During head-free gaze shifts, the first spike in this cell's discharge led ipsilateral gaze movements on average by 54 ms, and the duration of the entire discharge (prelude + burst) was as well correlated with the duration of the entire gaze shift as with the duration of the ocular saccade component (R = 0.77 vs. 0.76). The inability of discharge duration to specify whether gaze or eye is being encoded by this LLIBN is typical of these cells and, as will be seen later, is due to "noise" created by the prelude discharge. In the three examples illustrated in Fig. 3, the cell continued to fire until the gaze shift ended even in the plateau phase when eye velocity was near zero. This result was similar to the SLIBN described previously. Indeed, all the IBNs we recorded discharged during the plateau phase of the eye movement and throughout the terminal portion of the gaze shift where the eye movement actually reversed its direction.
As noted in companion paper I (Cullen and Guitton 1997a
), the IBNs we recorded burst most robustly for ipsilateral "ON-direction" (OND) saccades and discharged far less, if at all, during pure vertical and contralateral "OFF-direction" (OFFD) saccades. Figure 4 compares, by using the same scale, the gaze OND and OFFD responses for neurons L0702 and H0409, introduced in Figs. 2 and 3, and neuron H0925, to be considered in a subsequent section. The former two cells had negligible OFFD discharges in both the head-fixed and -free conditions. The latter neuron was unusual in that its direction tuning curve was located well above zero during contralaterally directed saccades. This cell had the best OFFD responses in our population for both the head-fixed and head-free conditions. Note in Fig. 4 the clear correlation for this cell between firing frequency and gaze velocity in the ONDs and OFFDs. We will consider this phenomenon more quantitatively in a subsequent section.

View larger version (18K):
[in this window]
[in a new window]
| FIG. 4.
IBN firing rate during ON-direction (OND) and OFF-direction (OFFD) gaze shifts for 3 example IBNs (SLIBN L0702 and LLIBNs H0409 and H0925). A and B: units L0702 and H0409 were typical of 1/2 of the cells in our population. There was a strong similarity between trajectories of neurons' firing frequency and horizontal gaze velocity during gaze shifts directed ipsilateral to the recording site in the cell's OND. In contrast, these cells discharged few, if any, spikes during contralateral gaze shifts in cell's OFFD. C: by comparison, unit H0925 demonstrated robust OFFD response that resembled an attenuated version of OND response. Response of this cell was a particularly striking example of behavior of the other1/2 of cells in our population that demonstrated a substantial OFFD response. Abbreviations are the same as in Fig. 2. Vertical calibration bar denotes 400 spikes/s, 400°/s, and 40°.
|
|
We employed the same technique, as those used in our head-fixed analysis, to obtain an accurate measure of each IBN's optimal direction with respect to the direction of both the ocular component of the saccade and the entire gaze shift (Fig. 5). Each neuron's optimal direction was calculated by fitting a nonlinear function to relate the NOS in the unit discharge to polar angle of the gaze and eye movements. The gaze shift amplitudes used to determine the preferred direction ranged from 30 to 40°. Example fits are illustrated for LLIBN H0409 in the inserts in Fig. 5, A and B (eye and gaze, respectively). As in the case of head-fixed saccades, the IBNs fired most strongly in relation to eye movements directed ipsilaterally and closely aligned with the horizontal plane (Fig. 5A). Likewise, IBNs fired preferentially for gaze shifts aligned with the horizontal plane (Fig. 5B). The mean OND of the IBNs was 3.5 ± 21.1° (SD) and
1.9 ± 18.8° with respect to eye and gaze movements, respectively. This result can be compared with the value
1.0 ± 15.2° found in the head-fixed condition. In the present analysis the calculated preferred directions (relative to either gaze or eye movements) varied somewhat more than in the head-fixed analysis (gaze: 51° upward to 68° downward); however, this additional variance was accounted for by only two neurons (7% of our sample). As in the previous study (Cullen and Guitton 1997a
) the population average of preferred directions approximated well the ipsilateral horizontal; consequently, the subsequent analyses in this study were confined to gaze, eye, and head movements directed horizontally or nearly so (see METHODS).

View larger version (34K):
[in this window]
[in a new window]

View larger version (36K):
[in this window]
[in a new window]
| FIG. 5.
Preferred direction of each of the IBNs in our population is illustrated in polar coordinates; the length of the vector indicates the sensitivity of each IBN measured in spikes/deg as indicated by scale to left of each panel. A: average preferred direction (heavy arrow) of our sample of IBNs was aligned with ipsilateral eye movements that were directed 4° above horizontal. Inset: direction tuning of example IBN (LLIBN H0409) is illustrated by a plot of the number of spikes (NOS) vs. eye movement direction. B: average preferred direction (heavy arrow) of the entire sample of IBNs was aligned with 2° of ipsilateral gaze displacement. Inset: direction tuning of LLIBN H0409 is illustrated by a plot of the NOS vs. gaze direction. Gaze shifts included in analysis ranged from 30 to 40°. Optimal direction for this cell was determined by fitting data to a sum of 2 sinusoids by using nonlinear least squares.
|
|
Estimation of lead time
The time by which a BN's discharge led a gaze shift was calculated by the two methods used in companion paper I (Cullen and Guitton 1997a
). Briefly, lead time was determined by 1) calculating the period between the onset of the first spike and the onset of gaze velocity and 2) shifting the unit discharge in time td until an optimal model fit (i.e., the maximal VAF) was obtained for the simple dynamic model
|
(1)
|
where r is a bias term, b1 is a gain term, and
is the input that can either be eye, head, or gaze velocity.
Figure 6 shows the results of applying the latter procedure to cell L0702. The relationship among VAFs obtained by Eq. 1 and lead times used relative to td are illustrated. The values of optimal dynamic latency td for this cell were 12, 11, and 13 ms for the eye-, head-, and gaze-based model fits, respectively. The figure shows that a clear maximum VAF can be identified for eye-, head-, and gaze-based versions of Eq. 1. It is important to note that correlations between the entire profile of the temporally shifted burst and the gaze or eye movement trajectories, not only the beginning of the burst, were used to determine the optimal value of td. Lead times on each trace are indicated that would be coincident with the onset of the burst relative to the indicated movement (
).

View larger version (17K):
[in this window]
[in a new window]
| FIG. 6.
Decreasing variances accounted for (VAFs) that arise when burst is shifted by times other than the optimal. Data are illustrated for cell L0702. Optimal dynamic lead time td is defined as value for which maximal VAF is obtained using model 2d (Eq. 1) to fit the firing frequency profile. Values of td for this cell were 12, 11, and 13 ms for eye-, head-, and gaze-based models 2d, respectively. A clear optimal lead time value (latency = 0 on the abscissa) can be identified for each eye, head, and gaze model. Arrow, lead times that would correspond to the burst not being shifted in time relative to each component.
|
|
The onset of gaze and eye movements are coincident in most gaze shifts because the early portion of the gaze shift is dominated by the eye movement. Therefore the onset of gaze and eye velocity can be considered to be equivalent. The distributions of lead times determined from the difference in timing between the onset of the first spike and the onset of gaze (or eye) velocity are shown in Fig. 7A. The distributions of the dynamic lead times that were estimated by using either eye or gaze velocity, respectively, as input to Eq. 1 are illustrated in Fig. 7, B and C. A comparison of the histograms in Fig. 7, C and B, with that in Fig. 7A reveals that the dynamic lead times have a much narrower distribution than those based on the first spike.

View larger version (21K):
[in this window]
[in a new window]
| FIG. 7.
Histograms illustrating lead time by which IBN discharges led head-free gaze shifts. Lead times were obtained via 2 methods. A: time between onset of 1st spike and onset of gaze shift. By using this method, the onset of 1st spike led the onset of gaze shift by 21.0 ± 11.4 ms. B and C: to obtain an estimate of IBNs' dynamic lead time, unit discharge wasshifted in time until an optimal model fit was obtained for model 2d, [B(t) =R + b1 ], where = eye velocity in B and gaze velocity in C. By usingthis method, the BN discharge led gaze velocity and eye velocity by 12.6 ±2.9 ms and 13.7 ± 3.6 ms, respectively. Measures of dynamic lead time in (B) and (C) provided, in general, lower values of lead time than conventional method that utilized the onset of 1st spike (A). N, total number of neurons.
|
|
Table 1 shows a comparison of the average lead times for SLIBNs and LLIBNs obtained by the first-spike and dynamic lead time methods. These numbers can be compared either along rows (for each neuron group, the first-spike method vs. the dynamic method) or along columns (for each method the difference between SLIBNs and LLIBNs). Let us first consider the rows in Table 1. In our accompanying head-fixed study, lead times for LLIBNs determined by the onset of the first spike were significantly longer than dynamic lead times. The values were comparable for SLIBNs. For head-free gaze shifts, the dynamic lead times (both gaze and eye velocity based) were significantly less (P < 0.005) for both SLIBNs and LLIBNs than mean lead times determined from the onset of the first spike. For the combined SLIBN and LLIBN groups, the first spike in the discharge led head-free gaze shifts on average by 21.0 ± 11.4 ms, (see also Fig. 7A) whereas the dynamic lead time analysis gave 12.6 ± 2.9 ms (Fig. 7B) and 13.7 ± 3.6 ms (Fig. 7C) for eye and gaze inputs, respectively. For the columns in Table 1, as in the case of head-fixed saccades, the estimated dynamic lead times (either gaze or eye velocity based) for SLIBNs were significantly less (P < 0.005) than those estimated for LLIBNs (gaze input, 11.7 ± 2.4 vs. 16.5 ± 3.1; eye input, 11.1 ± 2.7 vs. 14.8 ± 2.1).
For head-free gaze shifts, the relationships among the dynamic lead times that were calculated on the basis of gaze and eye velocity and those determined by the onset of the first spike are compared for each neuron in Fig. 8. The dynamic lead time was better correlated with the lead time of the first spike when the dynamic model was gaze velocity based (R = 0.87) than when it was eye velocity based(R = 0.78). Furthermore, the slope of the linear regression line was <1 because, as also shown in Fig. 7, the range of dynamic lead times was much less than that based on the onset of first spike.

View larger version (26K):
[in this window]
[in a new window]
| FIG. 8.
Relationship between lead time values obtained by 2 separate methods used in this study (see text). Dynamic lead times calculated by using both eye- ( ) and gaze- ( ) based models were well correlated with lead times determined with onset of 1st spike (R = 0.78 and 0.87, respectively). This result indicates that neurons that began to fire just before gaze shifts tended to have shorter dynamic lead times than those IBNs with longer preludes.
|
|
When the onset of the first spike was used to determine lead times, the values determined for head-free gaze shifts were well correlated (R = 0.89) but significantly larger (P < 0.005) than those determined for head-fixed saccades across our sample of neurons (Table 1). However, the significantly greater first-spike lead time during head-free versus head-fixed gaze shifts was not accompanied by corresponding differences in the estimated dynamic lead times. A comparison of the lead times calculated in head-fixed and -free conditions for each unit further revealed that the head-free dynamic lead time was not particularly well correlated with the head-fixed value(R = 0.54, gaze-based model; R = 0.66, eye-based model). This observation probably reflects the fact that the estimated values of dynamic lead times (gaze- and eye-based) for all neurons were rather similar and consistently covered a small range (~12 ms) compared with the broad range (~52 ms) determined by using the onset of the first spike.
Relationships between IBN activity and saccade metrics
For the head-fixed condition we reported significant relationships between NOS in the IBN burst and saccade amplitude and between burst duration and saccade duration. We now consider analogous relationships in the head-free condition. These relationships are summarized in Table 2 and Fig. 9 and to facilitate comparison, data from the head-fixed analysis of the same population of neurons are included in the table.

View larger version (37K):
[in this window]
[in a new window]
| FIG. 9.
Relationships among burst duration and duration of ocular component of gaze shift (A), duration of head movement (B), and duration of entire gaze shift for our population of SLIBNs (C; see Table 1). Heavy lines, average regression fit for all SLIBNs. Inset: data for example SLIBN L0702. Relationship among burst duration and duration of ocular component of gaze shift (D), duration of head movement (E), and duration of entire gaze shift for our population of LLIBNs (F). Heavy lines, average regression fits for all LLIBNs (see Table 1). Inset: data for example LLIBN H0409.
|
|
BURST DURATION VERSUS EYE, HEAD, OR GAZE DURATION.
Figure 9, A-C shows regression plots of burst duration versus ocular saccade duration, head movement duration, and gaze shift duration for each SLIBN. Insets show data for SLIBN L0702. Table 2 shows average values. In general, burst duration (DB) was well correlated with both gaze (Dg) and (De) eye saccade duration (mean R = 0.90, range = 0.73-0.98; R = 0.90, range = 0.61-0.99, respectively) whereas the correlation between head (Dh) and burst duration was significantly poorer (R = 0.63, range = 0.39-0.86). The mean intercept values calculated in the head-free analysis of burst duration versus either gaze (19.8) or eye saccade (10.0) duration were not significantly different from the intercept value (10.6) calculated in the head-fixed analysis of burst versus eye saccade duration. Furthermore, the average slope of the SLIBN burst-gaze saccade duration relationship (0.81) was not significantly different from that observed during head-fixed saccades (0.87). By comparison, the mean slope of the relationship between burst duration and head-free saccadic eye movement duration (0.95) was significantly higher (P < 0.005) than that observed during head-fixed saccades. In summary, according to these tests the head-free relation between burst and gaze shift duration was statistically similar to the head-fixed SLIBN burst eye saccade duration relationship.
The results obtained for LLIBNs are shown in Table 2 and Fig. 9, D-F. Insets show data for cell H0409. As was the case for SLIBNs, LLIBN burst duration was well correlated with both gaze and eye saccade duration (mean R = 0.85, range = 0.68-0.97; R = 0.85, range = 0.71-0.95, respectively), whereas the correlation between head movement and burst duration was significantly poorer (R = 0.58). The mean slope of the burst-gaze saccade duration relationship (0.80) was significantly smaller (contrary to that found for SLIBNs), and the mean slope of the burst eye saccade duration relationship (1.07) was significantly larger than that observed during head-fixed saccades (0.97; P < 0.005). The intercept calculated in the head-free analysis of burst versus eye saccade duration (26.2) was closest to the head-fixed intercept (26.5). In summary, for LLIBNs no head-free relations were statistically similar to the head-fixed relations.
When the entire population of IBNs was considered, the slope of the relationship between burst and eye saccade duration head-fixed (0.91) was halfway between the head-free gaze (0.80) and eye (1.00) slopes. The population head-free eye intercept was the most similar to the head-fixed value. We conclude, from our analysis of burst duration, that this quantity does not indicate clearly whether IBNs as a population discharge for the duration of either saccadic gaze or eye motion. This ambiguity is due to the noise attributable to the preamble discharge in LLIBNs and the very similar durations of the eye and gaze saccades. However, when a clear distinction is available between the end of the eye saccade and gaze saccade, as in the middle panel of Fig. 2, then the burst clearly continues for the duration of the gaze shift and exceeds saccade duration.
NOS VERSUS EYE, HEAD, OR GAZE AMPLITUDE.
In our head-fixed analysis of IBNs we reported a robust correlation between NOS in the discharge and the horizontal amplitude of the ocular saccade. Figure 10 illustrates for the head-free condition the relationship between NOS and the amplitudeof the eye (ocular; top), head (middle), and gaze (bottom) components of gaze shifts for SLIBNs (Fig. 10, A-C) and LLIBNs (Fig. 10, D-F). Insets show data for the examples SLIBN (L0702) and LLIBN (H0409). Table 2 gives mean values. For SLIBNs, the NOS was best correlated with the amplitude of the entire gaze shift (mean R = 0.73) than with the amplitude of the eye (R = 0.62) and head (R = 0.59) components. Indeed, the NOS in a burst was better related to gaze than eye amplitude for 82% of the cells in this study. In addition, it is interesting to note that the NOS generated by SLIBNs was nearly as well related to gaze amplitude during head-free gaze shifts as to saccade amplitude during head-fixed saccades (mean R = 0.73 vs. 0.79, respectively). Note that the mean slope (0.89 NOS/deg) of the relationship between NOS and gaze amplitude was significantly less than that observed for saccades made head-fixed (1.1 NOS/deg, P < 0.005). By comparison the mean slope of the relationship between NOS and the amplitude of the ocular component of a head-free gaze shift (1.3 NOS/deg) was significantly larger than that observed during head-fixed saccades. The intercept values calculated in the head-free analysis of NOS versus gaze amplitude and NOS versus eye amplitude were similar and significantly greater (P > 0.005) than that calculated in the comparable head-fixed analysis.

View larger version (43K):
[in this window]
[in a new window]
| FIG. 10.
A-C: SLIBNs. Relationships among NOS in a burst and amplitude of ocular component of gaze shift (A), amplitude of head movement (B), and amplitude of gaze shift (C). Heavy lines, average regression fit for all SLIBNs (see Table 1). Inset: data for example SLIBN (unit L0702). D-F: LLIBNs. Relationships among NOS in a burst and amplitude of ocular component of gaze shift (D), amplitude of head movement (E), and amplitude of gaze shift (F). Heavy lines, average regression fit for all LLIBNs (see Table 1). Inset: data for example LLIBN (unit H0409).
|
|
Results obtained from the analysis of LLIBNs (Fig. 10, D-F) yielded similar results. The NOS in a burst was best correlated with the amplitude of the entire gaze shift (R = 0.68) and more weakly correlated with the amplitude of the eye and head components (R = 0.60 and 0.53, respectively). The NOS gaze amplitude correlation was not as good as that obtained head-fixed. The slope of the relationship between NOS and gaze amplitude (0.65 NOS/deg) was significantly less than that between NOS and the amplitude of saccades made head-fixed (0.95 NOS/deg, P < 0.05). Unlike SLIBNs, however, there was no significant difference in the slopes of the relationship between NOS versus ocular saccade amplitude during head-free gaze shifts and head-fixed saccades.
When data from the entire population of IBNs were combined, NOS was best related to gaze amplitude (R = 0.71, range = 0.26-0.86) and more weakly correlated with eye amplitude (R = 0.61, range = 0.15-0.87) and head amplitude (R = 0.57, range = 0.21-0.84). Furthermore, the slope of the relationship NOS versus amplitude was significantly lower for gaze (0.79 NOS/deg) and higher for eye amplitude (1.14 NOS/deg) than it was in the corresponding head-fixed relationship (1.03 NOS/deg). We conclude from our analysis of NOS that this quantity is best correlated to head-free gaze shift amplitude, and furthermore the correlation coefficients for head-free gaze and head-fixed eye saccades are nearly the same. These results suggest that NOS encodes gaze amplitude, but note that the slope NOS/deg for head-free gaze is less than that for head-fixed eye.
NOS VERSUS EYE AND HEAD AMPLITUDE.
In general, the results of the foregoing analysis of the relationships between primate IBN discharges and the metrics of head-fixed saccades and head-free gaze shifts bear a striking resemblance to our previously published results in the cat (Cullen et al. 1993
). We can probe further the question of whether the discharges of our IBNs were gaze related by carrying out a multiple regression analysis. For this analysis we fit the following equation
where a, b, and g are constants and
E and
H represent the amplitude of the eye and head components at the end of the gaze shift, respectively. The results of this analysis are illustrated in Table 2. If the IBNs in our study were gaze related, then the values of b and g should be nearly equal. We found that this was, in fact, the case; there was no significant difference between the values of b and g estimated for either SLIBNs (1.00 ± 0.58 vs. 0.78 ± 0.39), LLIBNs (0.77 ± 0.31 vs. 0.54 ± 0.44), or the entire population of cells (0.90 ± 0.49 vs. 0.68 ± 0.42). In addition, there was only a slight improvement in the correlation coefficient (Table 2, mean R = 0.75) versus the value obtained in our prior analysis of the relationship between NOS and gaze amplitude (Table 2, mean R = 0.71). Finally, there was no significant difference between the values of the y-intercept obtained from this multiple regression analysis (Table 2, 14.0 ± 9.3) and that obtained from the comparable head-fixed relationship, NOS versus saccade amplitude (Table 2, 10.4 ± 4.2). Taken together, these findings lend support to the hypothesis that the discharges of IBNs are gaze related during head-free gaze shifts. This hypothesis is examined in more detail in the next section and the implications of these findings are addressed in the DISCUSSION.
Dynamic models of IBN discharge during head-free gaze shifts
GAZE- VERSUS EYE-BASED DYNAMIC MODELS.
We investigated the same downstream models as those studied in companion paper I (Cullen and Guitton 1997a
), but now IBN spike train dynamics were estimated by using either gaze, eye, or head movement trajectories as inputs. We will show in a later section that head movement parameters were poor (as might be concluded from the poor relationship between head movement amplitude and NOS shown in the previous analysis) at predicting burst dynamics, so we focused on the link between burst and either eye or gaze parameters. Three representative models along with the average parameter values associated with either gaze or eye motion as inputs are listed in Tables 3 (models 1d, 2d, and 7d). The numbering of each model corresponds to that of the comparable model tested in companion paper I (Cullen and Guitton 1997a
). The parameter estimation methods were described in detail elsewhere (Cullen et al. 1996
). In the analysis of all models that were tested we found that there were no notable differences between the averaged VAF and BIC values for the subpopulations of SLIBNs and LLIBNs. Therefore in our description of model fits we will consider average values for the entire population of IBNs.
The simplest model: is firing frequency proportional to eye or gaze velocity? The results of applying the most simple model in our series of downstream models (Table 3, model 1d) to the analysis of head-free IBN discharges are illustrated in Fig. 11 for cell L0702. For this neuron and all others in our study model 1d was much better at predicting the dynamics of IBN firing when gaze velocity, rather than eye velocity, was the model input. The IBN continued firing throughout the entire gaze shift (onset and offset indicated by vertical dotted lines) although eye velocity sometimes approached 0°/s
the so-called plateau phase
more than 100 ms before the gaze shift ended. This is particularly evident in the middle column where eye velocity also reversed direction just before the gaze shift ended. The firing rate predicted from eye velocity was much too low, starting particularly from the middle of the gaze shift onward.

View larger version (28K):
[in this window]
[in a new window]
| FIG. 11.
Examples of model fits to SLIBN (L0702) firing frequency during head-free gaze shifts. Shaded trajectories, unit firing frequency profiles; superimposed curves, fit produced by each model. Top panel: simple modelB(t) = b1 produced significantly better fit (heavy line) of the actual firing rate (shaded curve) than the model B(t) = b1 (2nd panel). Values of parameters, estimated by using >40 gaze shifts, provided below each model fit for this neuron. Bottom traces: accompanying gaze and eye velocity trajectories shifted in time by estimated optimal dynamic lead time td.
|
|
The VAF and BIC for model 1d are listed in Table 4 for short-lead, long-lead, and the combined population of IBNs. (The values for the parameters in Table 3 will be considered in the DISCUSSION where we will evaluate their physiological significance as well as their relation to the head-fixed data.) A comparison of the VAF values indicates, not surprisingly after viewing Fig. 11, that the gaze velocity-based model 1d provided a substantially better prediction of IBN firing than the eye velocity-based model (population meansVAF =
0.04 vs.
0.35, respectively). Note, however, the negative VAFs. Thus in spite of the better gaze prediction, as in the case of head-fixed saccades, this simple model was objectively poor at predicting the firing behavior of IBNs. Indeed, the negative VAF, for both SLIBNs and LLIBNs, when either gaze or eye velocity were the model inputs, indicates that the fit was worse than that about a constant mean value.
Importance of a bias term to gaze and eye velocity-based models. The addition of a bias term (r) to model 1d resulted in model 2d (Table 3) and greatly improved the prediction of IBN firing in both gaze and eye velocity-based models. This result is the same as that found in our preceding analyses of these neurons during head-fixed saccades. Comparison of Fig. 12 with Fig. 11 illustrates for cell L0702 the improvement in the fit of this SLIBN's firing rate that accompanied the addition of a bias term when either gaze or eye velocity was the input to the model. (Note the different bias values in these 2 cases.)

View larger version (28K):
[in this window]
[in a new window]
| FIG. 12.
Model fits to SLIBN (L0702) activity for same 3 example gaze shifts illustrated in Fig. 11. When bias term was added to model illustrated in Fig. 11 there was a significant improvement between goodness of the model fit when either gaze (top panel) or eye (2nd panel) velocity were used as inputs to model. Values of the parameters, estimated by using >40 gaze shifts, provided below each model fit. Gaze- and eye-based models provided comparable fits to data; however, bias term estimated when gaze velocity was input was similar to that estimated for the head-fixed data (for this example neuron, mean 184 vs. 175) (see companion paper I Cullen and Guitton 1997a ), whereas that estimated when eye velocity was input was significantly larger (255).
|
|
The VAF increased substantially from
0.04 to 0.20 for the gaze velocity-based model 2d, and from
0.35 to 0.22 for the eye velocity-based model (Table 4). There was no significant difference between the ability of the gaze and eye velocity-based models to predict the firing discharge of IBNs. In effect, the addition of a bias to the model effectively permitted the velocity term in the eye velocity-based model to contribute to predicting the modulation of the IBNs discharge once the eye had reversed direction toward the end of the gaze shift (Fig. 12). Model 2d provided a fit of the IBN spike trains that was as good, on average, as that described by a correlation coefficient (R) of 0.46 in a linear regression analysis. The accompanying relatively important decrease in the average BIC value between models 1d and 2d [population means 10.07 vs. 9.54 (gaze) and 10.56 vs. 9.50 (eye)] states that the addition of a bias term to both the eye and gaze models was warranted. It should be noted, however, that whereas the bias value estimated for the gaze model was comparable to that determined for head-fixed saccades (222.9 ± 96.6 vs. 222.6 ± 64.6, respectively), the bias value estimated for the eye model was significantly larger (279.3 ± 107.9; P < 0.005). The reasons for this are due to the influence of head velocity and will become evident in the DISCUSSION.
Importance of acceleration, higher-order nonlinear and position terms. As was the case in the dynamic analysis of IBN discharges during head-fixed saccades, the addition of acceleration, higher-order velocity, and position terms to model 2d had little effect on increasing the VAF (or decreasing the BIC) of either gaze or eye velocity-based head-free models for the population of IBNs. The addition of an acceleration term to model 2d resulted in no change in the VAF (0.20 for gaze-based models; 0.22 for eye-based models). Furthermore, when higher-order terms (squared and cubed velocity terms) as well as an acceleration term were added to model 2d (to approximate the nonlinearity proposed by Van Gisbergen et al. (1981)
the VAF only increased very slightly (0.22 for gaze-based models; 0.24 for eye-based models). Finally, the addition of a position term to model 2d did not produce a marked improvement in model fits for either gaze- or eye-based models (VAF = 0.21 and 0.23, respectively). Consequently, as in the case of the head-fixed analysis, in downstream head-free models it appears that the bias term is much more critical than either acceleration, higher-order velocity, or position terms for representing IBN discharges.
Estimation of pole term, initial conditions, and variable bias. The addition of a pole term (the derivative of the IBN firing rate) and acceleration term to model 2d improved the model fits slightly when initial conditions (ICs) were taken directly from the data parameters (VAF = 0.25 for gaze-based models; 0.27 for eye-based models). However, the ability of this model to estimate IBN discharges improved considerably when ICs were estimated as parameters(VAF = 0.43 for gaze-based models; 0.43 for eye-based models). Such a model is comparable to model 6d in our investigation of head-fixed models (companion paper I, Cullen and Guitton 1997a
). It is important to note that this model required a large number of additional parameters, because a separate IC was estimated for each gaze shift. As was the case in our accompanying investigation of head-fixed model 6d, an important limitation of the pole-term model was that the estimated bias was highly variable and most often a large negative value (mean bias =
2,147 ± 4,314 vs.
1,166 ± 3,145 for gaze- and eye-based models, respectively). As discussed in companion paper I, it is unlikely that such a large negative bias has any physiological significance.
The value of model 6d in companion paper I (Cullen and Guitton 1997a
) was that it led to model 7d and ultimately to 8d. Accordingly, we also tested model 7d using eye or gaze inputs (Table 3). For this model the bias term was estimated separately for each gaze shift. The method of parameter estimation was less complex than it was for models that include a pole term for which the ICs are estimated as parameters (see Cullen et al. 1996
), although the number of parameters is comparable. Model 7d provided significantly better fits than model 2d (Table 4: VAF = 0.40 for gaze-based models; 0.40 for eye-based models). Comparison of Fig. 13 with Figs. 11 and 12 illustrates for cell L0702 the impressive fit of IBN firing rate that accompanied the addition of a variable bias by using model 7d. As was the case for model 2d, the mean bias value estimated for the gaze-based model was comparable to that determined for head-fixed saccades (209.2 ± 91.1 vs. 222.5 ± 65.3, respectively), whereas the bias value estimated for the eye-based model was significantly larger (274.5 ± 109.4; P < 0.005).

View larger version (27K):
[in this window]
[in a new window]
| FIG. 13.
Model fits to SLIBN (L0702) activity for same 3 example gaze shifts illustrated in Figs. 11 and 12. When model included velocity and variable bias terms there was significant improvement between goodness of the model fit when either gaze (top panel) and eye (2nd panel) velocity was input to model. Gaze- and eye-based models provided comparable fits to data. Values of parameters, estimated by using >40 gaze shifts, provided below each model fit. Values of rk listed represent the bias estimated for each of 3 gaze shifts. Mean of the variable bias term estimated when gaze velocity was input was similar to that estimated from the head-fixed data (for this example neuron mean 181 vs. 169) (see companion paper I, Cullen and Guitton 1997a ), whereas mean bias estimated when eye velocity was input was significantly larger (244).
|
|
Dependence of ICs and/or biases on gaze shift metrics. In our head-fixed investigation we demonstrated that for 50% of our cells the estimated ICs in model 6d (Table 2, companion paper I, Cullen and Guitton 1997a
) and the estimated bias terms in model 7d (Table 2, companion paper I) were correlated with the metrics of saccades. In the head-free analysis of these neurons we found a similar phenomenon; the estimated ICs and/or biases were correlated with the metrics of eye and gaze saccades, for gaze-based models in 5 of 16 SLIBNs and 9 of 12 LLIBNs (
50% of the whole population). Seventy-one percent of the cells (5 of 9 SLIBNs and 5 of 5 LLIBNs) that demonstrated head-fixed correlations between saccade metrics and ICs (and/or variable biases) also demonstrated correlations between eye and gaze saccade metrics and these parameters in the head-free models 6d and 7d. The best correlations were between the bias term of the gaze-based model (rkg) and eye and gaze saccade amplitudes. Significant correlations were also found with peak eye and gaze saccade velocity, but it is well known that eye-gaze amplitude and velocity are correlated during gaze shifts. By comparison, the estimated ICs or biases were correlated with the metrics of ocular and gaze saccades for eye-based head-free models 6d and 7d in only 2 of 16 SLIBNs and 4 of 12 LLIBNs (21% of the population).
Some of these relationships are illustrated in Fig. 14 for the same example IBN (H0925) shown in companion paper I (Cullen and Guitton 1997a
). For this neuron the estimated bias term (rkg) in gaze velocity-based model 7d was well correlated with eye and gaze saccade amplitudes (Figs. 14, A and B; 
, amplitudes >10°; R =
0.54 and
0.57, respectively) but much less to peak gaze velocity (Fig. 14C; R =
0.39). By comparison, the bias term (rke) in eye velocity-based model 7d was only weakly correlated with eye amplitude (Fig. 14D; R =
0.33) and even more poorly related to peak eye velocity (not shown, R =
0.13). The filled symbols in Fig. 14, A and B are for gaze shifts accompanied by some head motion(
), generally for gaze amplitudes >10°. We also analyzedsmaller gaze shifts (<10°) for which the head, although completely unimpeded, did not move. These points (open symbols) lie near the head-fixed line [- - -; taken from our analysis in companion paper I (Cullen and Guitton 1997a
)] for this cell, suggesting a nonlinear bias versus gaze amplitude relationship related to the nonlinear contribution of head motion to gaze amplitude. This finding is further addressed in the following section.

View larger version (28K):
[in this window]
[in a new window]
| FIG. 14.
Analysis of relationship between multiple biases (in model 7d) for LLIBN H0925 during head-free gaze shifts. Estimated bias term (rkg) in gaze velocity-based model 7d was correlated with eye (A) and gaze (B) amplitudes but much less to peak gaze velocity (C). In contrast, the bias term (rke) in eye velocity-based model 7d was only weakly correlated with eye amplitude (D). Filled squares and circles, bias values for gaze shifts accompanied by some head motion; open square, bias values for small gaze shifts in which the head did not move; solid line, fit through filled squares and circles; dashed line, taken from Fig. 13C of companion paper I (Cullen and Guitton 1997a ) and shows bias values in head-fixed condition. Note that the open square lie near this line. Majority of cells (71%) that demonstrated head-fixed correlations among saccade metrics and variable biases also demonstrated correlations between eye and gaze saccade metrics and biases in head-free condition.
|
|
We also investigated the relationships between the estimated ICs and/or biases and the position of either gaze or eye immediately before and after the saccade. Only two IBNs demonstrated a relationship; this was between ICs or biases (models 6d and 7d, gaze and eye velocity-based) and the initial position of the eye before the onset of the gaze shift.
EYE- AND HEAD-BASED DYNAMIC MODELS.
We showed in the previous section that only model 1d showed a significant difference in VAF between gaze-based and eye-based models (Table 4; VAF =
0.04 vs.
0.35, respectively). However, the mean VAF values that were obtained using model 1d were negative, indicating that neither gaze- nor eye-based models provided a useful description of IBN discharges. In models that included a bias term, including models 2d and 7d, the bias value estimated for gaze-based models was comparable to that in our head-fixed analysis, whereas that estimated for the eye-based models was significantly larger. One possible explanation for this result is that the increased bias required in the eye-based models reflected an unmodeled head velocity input onto these neurons. To further investigate this possibility, the difference in bias obtained in gaze- and eye-based models (model 7d) was compared with the mean head velocity during each gaze shift. We found that, for our sample of cells, the difference in biases was very strongly correlated with the mean head velocity during each gaze shift (R =
0.93 ± 0.15). This relationship is demonstrated for a typical neuron in Fig. 15. Based on these results we investigated a second series of models that include linear combinations of separate eye and head velocity terms (Table 5). These models address objectively the issue of whether IBNs carry both eye and head velocity signals during head-free gaze shifts.

View larger version (13K):
[in this window]
[in a new window]
| FIG. 15.
Difference between biases (rk, model 7d) estimated in gaze- and eye-based models plotted as function of mean head velocity during gaze shift for LLIBN H0925. This cell was typical of all cells in the study in that difference between estimated biases was very well correlated with the mean head velocity (R = 0.93).
|
|
Simplest model: is firing frequency described by sum of independent eye and head velocity inputs? The results of applying the most simple model, based on the sum of independent eye and head velocity inputs (model 1d), is listed in Tables 5 and 6. In this model the estimated head velocity coefficient (1.37 ± 0.52) was significantly larger than the estimated eye velocity coefficient (0.96 ± 0.28; Table 5). In addition, the estimated head velocity coefficient was more variable across the cell population than was the estimated eye velocity coefficient. Although model 1d predicted better the dynamics of IBN firing when a linear combination of eye and head velocity were used (Table 6; VAF = 0.03) than when either gaze or eye velocity was the only input to the model (Table 4; VAF =
0.04 vs.
0.35, respectively), the mean VAF for model 1d with independent eye and head velocity terms was nearly zero. This result indicated that the fit of model 1d was not much more effective than a fit about a constant mean value.
Importance of bias term to models with independent eye and head velocity inputs. The addition of a bias term (r) to model 1d resulted in 2d (Table 5) and greatly improved the ability of models that include linear combinations of independent eye and head velocity terms to predict IBN discharges (Table 6; VAF = 0.23). This result is consistent with our preceding analyses of these neurons during head-fixed saccades (companion paper I Cullen and Guitton 1997a
) and during head-free gaze shifts using gaze-based and eye-based models (see GAZE- VERSUS EYE-BASED DYNAMIC MODELS). In contrast to model 1d, the estimated head velocity coefficient (0.17 ± 0.42) was significantly less and more variable than the estimated eye velocity coefficient (0.52 ± 0.17; Table 5). Model 2d predicted the dynamics of IBN firing only slightly better when a linear combination of eye and head velocity were used (Table 6, VAF = 0.23) than when either gaze or eye velocity was the only input to the model (Table 4, VAF = 0.20 vs. 0.22, respectively). This small increase in VAF was accompanied by a small decrease in the average BIC value (Table 6, 9.47) relative to the gaze-based and eye-based model 2d that had one less parameter [Table 4, 9.54 (gaze) vs. 9.50 (eye)]. This decrease in BIC suggests that the addition of independent eye and head velocity terms to the model was warranted. The bias value estimated for the eye or head model was not significantly different from that determined for head-fixed saccades (255.8 ± 95.8 vs. 222.6 ± 64.6, respectively).
Inclusion of gaze shift and ocular saccade amplitude terms. Our finding that estimated ICs and biases were often correlated with saccade amplitude (companion paper I Cullen and Guitton 1997a
) and gaze shift amplitude (see GAZE- VERSUS EYE-BASED DYNAMIC MODELS) led us to construct a model that included an amplitude-dependent term. This model (Table 5, model 8d) is a simple extension of model 2d, to which a term proportional to either gaze amplitude (r1
G) or ocular saccade amplitude (r1
E) was added. This modification resulted in a small improvement in VAF for the population of neurons [Table 6, VAF = 0.23 in model 2d vs. 0.25 in 8d (gaze) and 0.25 in 8d (eye)]. The increase in VAF was accompanied by a small decrease in the average BIC value (Table 6, BIC = 9.44) in both gaze- and eye amplitude-dependent models relative to model 2d (BIC = 9.47), which suggested that the addition of an amplitude term to the latter model was warranted. As in model 2d, the estimated head velocity coefficient g1 [Table 5, 0.26 ± 0.34 (gaze-amplitude model) and 0.25 ± 0.36 (eye-amplitude model)] was significantly less than the estimated eye velocity coefficient b1 [Table 5, 0.54 ± 0.17 (gaze-amplitude model) and 0.54 ± 0.17 (eye-amplitude model)]. Similarly, g1 was much more variable across the cell population than b1 and furthermore g1 was less for LLIBNs than for SLIBNs. Note also that the sign of the amplitude-dependent term was negative, similar to the head-fixed analysis, indicating that this component of the firing rate was inversely correlated with gaze shift and/or ocular saccade amplitude.
OFFD discharges
Figure 3 showed the OFFD of the three neurons we have extensively illustrated in this paper and in companion paper I (Cullen and Guitton 1997a
). The OFFD of SLIBN L0702 and LLIBN H0409 consisted of one or two spikes and could not be analyzed. Fifty percent of our IBNs gave similar weak responses. Of the remaining cells, H0925 had the best model fits and we will consider its characteristics in detail.
For cell H0925 in the head-free condition, the dynamic lead time in the OFFD was 12 ms compared with 18 ms in the OND. Recall from companion paper I (Cullen and Guitton 1997a
) that these values were 11 and 17 ms, respectively, in the head-fixed condition. The relation between NOS and gaze amplitude in the OFFD had a correlation coefficient of 0.84 compared with 0.80 when NOS was related to eye amplitude. The excellent correlation to gaze compares well with that to the eye in the head-fixed condition.
We calculated the fits for all models listed in Table 5. Table 7 compares for cell H0925 the characteristics of the model fits for models 1d, 2d, and 8d that contain independent eye and head velocity terms. Surprisingly, the OFFD fits were better than those for the OND for each model. Table 7 also shows that in model 2d the absolute value of the velocity gain was similar in OFFD (eye = 0.64; head = 0.32) as in OND (eye = 0.78; head = 0.44) and, as in the head-fixed analysis of this cell, the bias term was much less and essentially zero in the OFFD compared with the OND. In model 8d (gaze and eye) the eye and head velocity gains in the OFFD were similar to their values in model 2d, but unlike the OND for this neuron, the bias was small and no relation was found between the estimated bias values and gaze shift amplitude.
View this table:
[in this window]
[in a new window]
|
TABLE 7.
Parameter estimations for downstream models with independent eye and head velocity terms;
ON- vs. OFF-direction response (unit H0925)
|
|
 |
DISCUSSION |
In this paper we have applied the system identification methods described in Cullen et al. (1996)
to analyze the discharge of both SLIBNs and LLIBNs in the head-free monkey. Our approach has permitted an objective comparison between the discharges of these neurons in the head-fixed (see companion paper I Cullen and Guitton 1997a
) and head-free animals. Such a comparison is important to understand how eye and head motor signals are processed by the brain stem to generate gaze shifts. Our most important conclusion is that IBNs carry head as well as eye velocity-related signals and we will consider this point in the last summary section of this DISCUSSION.
Determining a neuron's lead time
As in our previous head-fixed analysis, we calculated the lead time for each neuron by determining the time interval (dynamic lead time) by which the unit discharge should be shifted relative to either eye velocity or gaze velocity to obtain an optimal fit by using an accurate dynamic model. Theoretically, the dynamic lead time should provide a better global estimate of the true physiological delay between a BN's discharge and the resultant movement, because it takes into account the entirety of a cell's discharge rather than just the first action potential. Accordingly, the range of the estimated head-free gaze and eye dynamic lead times was greatly reduced (8-21 ms) from that determined from the onset of the first spike (5-55 ms) (Fig. 7).
The first-spike lead time, measured head-free, was well correlated for each cell with the head-fixed value, but the mean onset of the first spike in both SLIBNs and LLIBNs was systematically greater, by ~3 ms, for head-free gaze shifts than for head-fixed saccades (Table 1). Scudder et al. (1988)
reported that, in the head-fixed monkey, onset lead times vary with saccade size and velocity. Accordingly, the slightly longer first-spike lead times in our head-free condition might have been due to the larger range of gaze shift amplitudes that naturally arise in this condition; on average the gaze shifts were 10° larger head-free compared with head-fixed. However, our dynamic lead time calculations show that, unlike the first-spike method, there was no significant difference between the head-fixed and head-free mean dynamic lead times, using either eye or gaze movement parameters as model inputs (Table 1). Interestingly, in spite of similar global means, on a cell by cell basis the head-free and head-fixed dynamic latencies were only weakly correlated in either the gaze- or eye-based models. This effect may be due to the small range in dynamic lead times.
Metrics
Our previous studies in the head-free cat (Cullen et al. 1993
, 1994
) demonstrated that the duration of IBN discharges are better correlated with gaze than eye saccade duration. In the cat the range of ocular motility is small and one frequently observes gaze shifts whose motion is carried by the head, the eye being immobile at a small position saturation in the orbit (Guitton et al. 1984
, 1990
). In the monkey some large gaze shifts demonstrated this pattern, but for most gaze shifts it was hard to measure precisely a difference between the end of rapid eye and gaze motion. Indeed, gaze shifts like those of the middle and right columns of Figs. 2 and 3, respectively, where the eye is immobile in the orbit during a portion of the gaze shift, were a minority. Hence the difference in correlation seen in the cat was not observed in monkey IBNs; the duration of primate IBN bursts during gaze shifts was equally well correlated with both gaze and eye saccade duration. By comparison, the correlation between head movement and burst durations was significantly poorer.
In the head-free cat (Cullen et al. 1993
) the total NOS in a burst was well correlated with gaze amplitude and generally poorly correlated with eye saccade amplitude. In agreement with these findings (Table 2), the NOS in the burst of head-free primate IBNs was better correlated with the amplitude of the entire gaze shift (all neurons R = 0.71) than with the amplitude of the eye movement (R = 0.61) or head movement (R = 0.57). The better correlation between gaze amplitude and NOS was found in 82% of the neurons. Furthermore, NOS was nearly as well correlated with the amplitude of head-free gaze shifts (R = 0.71) for all cells as it was with the amplitude of head-fixed saccades (R = 0.79) (companion paper I, Cullen and Guitton 1997a
). This suggests that monkey IBNs encode displacements of the visual axis in both the head-fixed and head-free conditions and is in line with our finding that dynamic models require
terms and also with the predictions of the model proposed by Galiana and Guitton (1992)
.
Interestingly, in both the present study and in our analysis of cat IBNs (Cullen et al. 1993
) we found that the slope of the relationship between NOS and gaze amplitude was significantly lower head-free (0.79) than head-fixed (1.03). This result states that a given NOS predicts, with excellent fidelity, a larger displacement of the visual axis when the animal's head is free to move compared with when it is fixed. This is due to the head velocity signal carried by IBNs and is discussed in the last section of this paper. This effect may also be related to a close coupling between bursts driving eye and head motion, evidence for which exists in the cat (Grantyn and Berthoz 1987
; Guitton et al. 1990
) and that we have modeled elsewhere (Galiana and Guitton 1992
; Guitton et al. 1990
).
In a previous study by Whittington et al. (1984)
two types of primate BNs were described: 1) gaze BNs (cells that generate the same NOS in a burst for a given displacement of the visual axis irrespective of whether the head is fixed or free) and 2) eye BNs (cells that generate the same NOS in a burst for eye movement of equal amplitude irrespective of whether the head is fixed or free). In their model, gaze-burst cells were converted to eye-burst cells by subtracting from the former a signal encoding
. In this schema gaze motion is identical in head-fixed and head-free conditions because the head's contribution to a gaze shift is effectively nulled by the subtraction of head velocity at the level of the eye-burst cells. This is generally not the case. Several behavioral studies have demonstrated that the gain of the vestibuloocular reflex, the source of the
signal, is clearly reduced during large gaze shifts; furthermore head-free gaze shifts do not follow the same trajectory as head-fixed saccades (Fuller et al. 1983
; Guitton and Volle 1987
; Laurutis and Robinson 1986
; Pélisson et al. 1988
; Tomlinson 1990
; Tomlinson and Bahra 1986a
,b
).
Our neurons did not belong to either of the categories proposed by Whittington et al. (1984)
. The relationships between NOS and either gaze or eye metrics head-free are clearly different from the head-fixed values. Recent studies by Tomlinson and Brance (1991)
of monkey gaze shifts during passive whole-body rotations corroborate the more complex nature of the signals carried by BNs in the primate IBN area. They reported that the relationship between the NOS in the burst and saccadic eye movement amplitude was similar in the head-fixed and head-free conditions, thereby labeling these cells as eye saccade-related BNs. However, in an additional observation difficult to reconcile with the former, they reported that the duration of the burst in these neurons was best related to gaze movement duration. The activity of our IBNs also differed from this more recent categorization of burst cells, because we found that the relationships between NOS and eye metrics was significantly different during head-free gaze shifts and head-fixed saccades. In particular, the results of our multiple regression analysis (Table 5) indicate that the activity of IBNs, as a population, is related to head as well as eye movements during head-free gaze shifts.
Significance of estimated parameters in dynamic models
Based on the NOS analysis, IBNs seem to better encode gaze than eye metrics in the head-free cat and monkey. In our analysis of cat IBNs, we noted that their dynamics appeared to be better related to the trajectory of gaze-in-space than to eye-in-head (Cullen et al. 1993
). In the present study, a visual inspection of some primate IBN firing frequency profiles (Figs. 2 and 3) also suggests that the discharge of these cells is better related to gaze; for example, IBNs discharge during gaze shifts when eye velocity is zero. However, an objective analysis of firing frequency using our optimization algorithms gives a more complex picture. Indeed, a comparison of the VAF and BIC values in Table 3 indicates that both gaze and eye velocity-based models 2d and 7d were equally good, albeit with different parameters, at predicting IBN activity for both SLIBNs and LLIBNs. Therefore the question arises that if gaze- and eye-based dynamic models are equally good at predicting firing frequency, then what criteria can be used to establish how abducens MNs "read" the incoming burst signal. We will consider this in the following sections.
BIAS TERM.
The addition of a bias term to model 1d led to 2d and considerably improved the fit of both gaze and eye velocity-based models. However, whereas the mean bias term estimated in 2d for the gaze velocity-based models was similar to that estimated for head-fixed saccades in the same population of cells (Table 3; 222.9 head-free vs. 222.6 head-fixed), the bias term in the eye velocity-based head-free model was significantly larger (279.3, P < 0.005). This same result was observed for the mean bias values estimated for each cell by model 7d.
But what is the significance of the bias term? In these models it is that part of the discharge that is not modulated by gaze or eye velocity. It could reflect the following three mechanisms: 1) the intercept of a piece-wise linear fit to a nonlinear relation between firing frequency and velocity, 2) the "resting" discharge that would result when the inhibition of BNs by brain stem OPNs is removed during saccades (and gaze shifts), and 3) the presence of unmodeled inputs onto eye muscle MNs that contribute unknown signals to these cells that influence eye motion in addition to the burst signal. Interpretation 1 has some merit because the head-fixed nonlinear model 4d produces a lower bias than 3d (companion paper I Cullen and Guitton 1997a
), but unfortunately in both the head-free and head-fixed versions of model 3d the nonlinearity did not improve the VAF. Interpretation 2's major weakness is that the bias is different in the ONDs and OFFDs of each neuron (Fig. 7) which would not be the case if it simply reflected a resting discharge uncovered by an omnidirectional pause in OPN activity. We suggest therefore that the bias reflects unmodeled inputs that are further considered in a subsequent section.
VELOCITY GAIN TERM.
By convention, the positive OND velocity gain term for all IBNs indicates that they increase their firing rate for gaze shifts (and the ocular component of gaze shifts) in the ipsilateral direction. As seen in Table 3, the gaze velocity gain coefficients determined for models 2d and 7d during head-free gaze shifts were significantly lower than those determined with the same model for head-fixed saccades. Recall that there was no accompanying difference in the estimated bias terms. The eye velocity gain coefficients determined for these models during head-free gaze shifts were also significantly lower than those determined for head-fixed saccades (but note, with a significant accompanying increase in the estimated bias term). These results show that the velocity gain parameter optimized from head-fixed data cannot be used to predict head-free IBN discharges on the basis of either eye or gaze dynamics. We will consider again this point in the next paragraphs, in relation to Fig. 17.

View larger version (28K):
[in this window]
[in a new window]
| FIG. 17.
Prediction of head-free discharges by using head-fixed datafor SLIBN (L0702). Firing frequency calculated by using head-fixedmodel (companion paper I, Cullen and Guitton 1997a ). B(t) = 176 +0.71 (t td) is shown by the heavy lines when gaze (top panel) and eye (2nd panel) are inputs ( ) to model. Using gaze velocity as input to head-fixed model results in overestimation of unit discharge, whereas using eye velocity as input to model results in underestimation of unit discharge. VAFs when gaze and eye were model inputs were 0.17 and 0.04, respectively. Conventions similar to Fig. 11.
|
|
AMPLITUDE-DEPENDENT TERMS.
In companion paper I (Cullen and Guitton 1997a
) we described a component of IBN firing that varied inversely with the amplitude of the saccadic eye movement, and we characterized this phenomena by using the estimated amplitude-dependent coefficient (r1) in model 8d (Table 2 of companion paper I, Cullen and Guitton 1997a
). During head-free gaze shifts the coefficient was also negative, but its absolute value was significantly less than that estimated in our head-fixed analysis of the same cells (Table 5, model 8d). The significance of this term is considered in companion papers I and III (Cullen and Guitton 1997a
,b
).
SIGNIFICANCE OF ESTIMATED PARAMETERS IN MODELS WITH INDEPENDENT EYE AND HEAD TERMS.
As emphasized previously, the bias value estimated for gaze-based models was comparable to that in our head-fixed analysis, whereas that estimated for the eye-based models was significantly larger. We found that the difference in bias between the gaze- and eye-based versions of model 7d was very strongly correlated with the mean head velocity during each gaze shift (Fig. 15). This result suggests that the increased bias in the eye-based models reflects an unmodeled head velocity input onto these cells. To pursue this argument further, a series of models that included independently estimated eye and head velocity terms was investigated (Table 5). These models show that IBNs as a population carry head as well as eye velocity signals during head-free gaze shifts.
A comparison between the estimated eye velocity coefficients (b1) obtained in our analysis of models with independent eye and head velocity terms (Table 5) versus those obtained in our head-fixed analysis (Table 2 of companion paper I, Cullen and Guitton 1997a
) reveals that the value estimated in the head-free analysis was significantly lower. This was true for each model that was tested. Figure 16A illustrates the relationship between the value of b1 for each unit during head-fixed saccades and head-free gaze-shifts for model 8d. The data from all but one cell lie below the plotted line of unity slope, indicating that the eye velocity sensitivity of IBNs is less during head-free gaze shifts than during head-fixed saccades.

View larger version (13K):
[in this window]
[in a new window]
| FIG. 16.
A: comparison of eye velocity gain term b1 estimated in head-free model 8d by using independent eye and head velocity terms with value of b1 estimated in our analysis of head-fixed saccades (companion paper I, Cullen and Guitton 1997a ) also using model 8d. IBNs generally carried a smaller eye velocity signal during head-free gaze shifts. B: relationship between the eye velocity gain b1 and head velocity gain g1 estimated by using independent eye and head velocity terms in head-free model 8d. Most cells carried a larger eye velocity signal than head velocity signal during head-free gaze shifts; some cells had no, or negative, head velocity gains. C: comparison of value of amplitude-related term r1 estimated by using independent eye and head velocity terms in model 8d, with the value of r1 estimated in our head-fixed analysis (companion paper I, Cullen and Guitton 1997a ) also for model 8d. Value of r1 was generally lower during head-free gaze shifts.
|
|

View larger version (27K):
[in this window]
[in a new window]
| FIG. 18.
A: extension of Fig. 18A of companion paper I (Cullen and Guitton 1997a ) that illustrates inputs to abducens nuclei on the left and right (ABDR and ABDL, respectively) that are active during rightward head-free gaze shifts. Specific contralateral inputs to ABDR arise from VNL-PHL (left vestibular-prepositus nucleus) and include I-PVPL (type I-position-vestibular-pause cells), RBPL (burst position), RPL, (position cells) and LPL. The 1st subscript denotes the OND of each cell type's response; final subscript denotes side of the brainstem where cell is located. Inputs from the ipsilateral (right, in this example) VN/PH include projections from I-PVPR, LBPR, LPR, LE/HR (eye/head-velocity cells), RBPR, and RPR. Excitatory and inhibitory projections are indicated by gray and black circles, respectively; SCCL and SCCR, left and right semicircular canals, respectively; NIC, nucleus intercalatus. Equations describing the ON [B(t)]- and OFF [B(t)]-direction responses of the BN population are based on results of our dynamic analysis. For additional details see text. B: simple schematic that illustrates necessity of removing from the ocular MNs the head velocity term carried by BNs during head-free gaze shifts. In addition, it is suggested that signal processing at the level of the MNs includes 1) removal of amplitude-dependent term E if transfer function between MNs and plant dynamics is in fact linear and 2) additional source of eye velocity signal during head-free gaze shifts to match head-fixed transfer function between IBNs and MNs.
|
|
We next compared the estimated eye and head velocity coefficients. In both models 2d and 8d the estimated
coefficients were significantly less than the estimated
coefficients. The relationship between the eye and head velocity coefficients is illustrated in Fig. 16B for model 8d. If our neurons encoded gaze then one would expect the data points to lie along the 45° equi-gain line. Some do, but as seen in Table 5, as a population our IBNs generally encoded a larger
signal than
signal. It is interesting that LLIBNs had a significantly lower head velocity sensitivity than the SLIBNs (P = 0.025). Figure 16B shows that five LLIBNs in our sample had an
sensitivity and nearly no, or a negative,
sensitivity. Except for lead-time, this is the only difference that we have seen between these two cell populations in all of our analyses (this paper; Cullen and Guitton 1997a
,b
). The standard error bars of the parameter estimates of b1 and g1 lie within the area of the points themselves in Fig. 16B. It follows that, for the vast majority of cells, the estimated head velocity sensitivity coefficient g1 was significantly different from zero. On average, the 95% confidence interval for g1 spanned approximately ±0.06 on either side of the parameter estimate. Hence, within our population of 28 cells, 26 had a significant head velocity signal (Fig. 16B). Furthermore, on average the 95% confidence interval for b1 spanned approximately ±0.04 on either side of the parameter estimate. Hence, the eye velocity sensitivity of all but one IBN was significantly less during head-free gaze shifts than during head-fixed saccades (Fig. 16A).
The findings illustrated in Fig. 16, A and B, further support our hypothesis that IBNs cannot be considered gaze or eye BNs in the classical sense. According to the schema proposed by Whittington et al. (1984)
, a BN described by the following equation during head-fixed saccades
|
(2)
|
where b1x is the head-fixed eye velocity gain, would be an eye BN if it were described by this same equation during head-free gaze shifts. If instead the head-free discharge of the neuron in Eq. 2 is described by the equation
|
(3a)
|
that can be rewritten as
|
(3b)
|
it would be categorized as a gaze BN. However, our analysis indicated that the best description of IBN activity was
|
(4)
|
where b1f and g1f are the eye and head velocity gains, respectively.
Our findings that b1f < b1x and g1f > 0 indicate that the IBNs in this study were not eye BNs. Furthermore, the findings that b1f > g1f (or put another way that b1f
g1f) and b1f
b1x indicate that the IBNs in this study were also not gaze BNs. This point is taken up again in our discussion of Fig. 17.
Comparison of the value of the amplitude-related term r1, estimated by using independent eye and head velocity terms in model 8d, with the value of r1, estimated in our head-fixed analysis (companion paper I Cullen and Guitton 1997a
) also for model 8d, demonstrated that r1 was generally lower during head-free gaze shifts (Fig. 16C).
Head-fixed models cannot be used to predict head-free discharges
We consider again the important question of whether head-fixed models can be used to predict head-free discharges. As we reviewed in the previous section, the classic view of the gaze control network is that a BN characterized head-fixed can become head-free, either a gaze burster that predicts gaze motion or, when head velocity is subtracted from its discharge, an eye burster that predicts eye motion. We concluded from our metric and dynamic analyses that IBNs cannot be classified according to this scheme.
To more clearly illustrate this point, the parameters calculated for the head-fixed model 2d were used to predict the discharge of IBNs during head-free gaze shifts, when either the corresponding gaze or eye velocity was taken as the input to the model. Figure 17 shows the result of this analysis for the example IBN L0702 that can be compared with the optimal model 2d fit in Fig. 12. When gaze velocity was the input to model 2d, the resultant fit using the head-fixed parameters specific to this cell provided an overestimation of the unit discharge (top panel), whereas when eye velocity was the input to the model the resultant fit was an underestimation of the unit discharge (2nd panel). This result applied to all the neurons in our study; head-free gaze velocity inputs to the head-fixed model consistently overestimated head-free IBN discharges, whereas head-free eye velocity inputs to the head-fixed model underestimated the discharge. These differences arose because b1f < b1x (preceding section) and also because of the lower bias value for head-fixed versus head-free eye-based models. These results corroborate the conclusions of our metric analysis showing that the slope of the relationship between the total NOS and gaze amplitude during head-free gaze shifts was significantly lower than that observed during head-fixed saccades. In summary, these observations taken together lead to the important conclusion that head-fixed models cannot be applied to the analysis of IBN discharges during head-free gaze shifts.
Conclusions: implications of IBN head-free model fits
In companion paper I (Cullen and Guitton 1997a
) we proposed that model 8d combined simplicity, physiological soundness, and good predictive capabilities to provide a VAF of 0.34 in its prediction of burst frequency. Accordingly, we concluded that the net IBN signal arriving at the ABD and vestibular nuclei-prepositus-hypoglossi nuclei (VN-PH) complex in the head-fixed monkey was
|
(5)
|
In this paper we have characterized the influence of the IBN signal on the ABD and the VN-PH complex during head-free gaze shifts. To obtain the net signal arising from the burst generator it is necessary, as we did in companion paper I (Cullen and Guitton 1997a
), to take account of both the OND and the OFFD discharges (Van Gisbergen et al. 1981
). It is important to recall that we found a difference between SLIBN and LLIBN sensitivity to
, the latter being weaker. We have not included this observation in our summary schema of Fig. 18 because it complicates matters considerably. In this section we have lumped together the properties of the two IBN subtypes. In the head-free condition, model 8d provided the following descriptions of the OND discharges of IBN
|
(6)
|
and
|
(7)
|
The VAFs given by these two representations are equivalent, but we will use Eq. 7 because the
E term is conceptually easier to consider given that IBNs are known to project to extraocular MNs.
As seen in companion paper I (Cullen and Guitton 1997a
), evaluating the OFFD response is a more difficult process because of a zero OFFD in ~50% of the cells. We proposed in companion paper I a method to estimate the mean OFFD by taking a proportion of the best OFFD response found in cell H0925. In the head-free condition, model 8d (gaze) for this cell gives the following equation
|
(8)
|
and model 8d (eye) gives the following equation
|
(9)
|
Note that we have not inverted the signs as in the Tables, because Fig. 18 shows the response of left IBNs to a positive (rightward) movement. Again, the VAFs given by these two representations are equivalent, and for the purposes of this discussion we will use Eq. 9.
With the use of arguments presented in the last section of companion paper I (Cullen and Guitton 1997a
), the postulated average OFFD r0 =
9/4 =
2; r1 = 0.8/4 = 0.2; b1 = 0.63/4 = 0.16; and g1 = 0.25/4 = 0.06. The population average OFFD response simplifies to
|
(10)
|
In Fig. 18 this small bias is neglected. Figure 18A is an extension of companion paper I's (Cullen and Guitton 1997a
) Fig. 16A to the head-free condition. If we assume that IBNs and EBNs carry the same signal, then the net excitation on ABDR and VN-PHR, or net inhibition on ABDL and VN-PHL, will be the difference between the mean OND and OFFD response (Van Gisbergen et al. 1981
)
|
(11)
|
or rounded off (Fig. 18)
|
(12)
|
A comparison of Eqs. 12 and 5 reveals that the bias input to the ABD was nearly identical in the head-fixed and -free condition, and that a component of IBN firing rate remains inversely related to the amplitude of the eye movement in both conditions. However, there were important differences that we observed: 1) a striking reduction in the eye amplitude and eye velocity gain-related signals that reach the motor nuclei during head-free gaze shifts and 2) an
-related signal contained in head-free IBN discharges (Fig. 18B).
In the head-fixed condition the net inhibitory IBN input to ABDL assures the silence of these antagonist MNs during rightward saccades. In the head-free condition the interesting conclusion arises that ABDL is also inhibited but now during the entire duration of a rightward gaze shift, even during a large gaze shift when the eye is immobile in the plateau phase and the change in gaze is accomplished by the movement of the head. (Recall that normally when the eye is not moving there is tonic activity in all MNs). On the agonist side, ABDR in this example, the question arises as to how the net excitatory BN signal is transformed into the eye motor command seen on MNs, given that the net excitatory signal on ABDR contains a signal proportional to
. Because the firing frequency of MNs in ABDR during right saccades is roughly proportional to
(Robinson 1970
) it follows that the
signal must be removed at MN. This is shown generally in Fig. 18B. Specifically, in Fig. 18A, a plausible mechanism for this is via the I-PVPL cells. During right saccades these cells, on average, do not pause completely (Scudder and Fuchs 1992
). We postulate that, in this context, I-PVPL will carry a signal from the left semicircular canal (SCCL) proportional to
that decreases with right
. A decreasing excitatory
signal in ABDR is theoretically equivalent to an increasing inhibitory
signal that is needed to cancel the
signal carried by the burst generation to MN.
The reader well versed in the gaze-control literature will certainly be troubled by these arguments. Indeed, an important property of the control of gaze shifts of amplitude
25° is that the SCC signal impinging on MN is thought to have a low gain (Pelisson et al. 1988; Tomlinson 1990
; reviewed in Guitton 1992
). For example, the experiments of Laurutis and Robinson (1986)
and Guitton and Volle (1987)
indicated that a sudden mechanically imposed reduction of
during a large gaze shift did not result in a simultaneous increase in
.
We have an explanation for this phenomenon that is compatible with a canal signal impinging on both B and MN: IBN firing frequency in the monkey and the cat is well related to
during passive (i.e., experimenter generated) head-on-body rotations (Cullen and Guitton 1995a
,b
). Most importantly, when passively generated head motion in the cat is mechanically braked, burst frequency decreases markedly with a lag time of 10-15 ms. These observations imply that, during a rightward gaze shift that is suddenly halted, the I-PVPL signal will increase to resting rate, but the burst generator signal will simultaneously decrease. The two signals may cancel at the MNs and eye motion will be minimally affected.
To summarize, our use of system identification methods to analyze IBN discharges shows that both SLIBNs and LLIBNs discharge with a frequency approximately proportional to the sum of
and
signals (LLIBNs may have a lower
sensitivity than SLIBNs), and this leads to the prediction that I-PVPL will carry a signal proportional to
during rightward gaze shifts. Furthermore, the
and
signals on BNs also get sent to the VN-PH complex. If this structure acts as a neural mathematical integrator (Cannon and Robinson 1987
; Cheron et al. 1986a
,b
; Cheron and Godaux 1987
; Escudero et al. 1992
; Lopez-Barneo et al. 1982
; McFarland and Fuchs 1992
; Mettens et al. 1994
) or plant model (Galiana and Guitton 1992
), then one might expect some cells in the PH of the head-free monkey to carry head- as well as eye-position signals.
 |
ACKNOWLEDGEMENTS |
We are grateful to H. L. Galiana for many helpful discussions and to C. G. Rey for contribution to the development of the analysis methodology and software. We thank Dr. J.A.M. Van Gisbergen for a number of useful suggestions that have helped improve this manuscript greatly; we also thank Dr. W. M. King for critically reading this manuscript.
This study was supported by the Medical Research Council of Canada, the National Institutes of Health, and the Human Frontiers Science Organization.
 |
FOOTNOTES |
Address for reprint requests: K. E. Cullen, Aerospace Medical Research Unit, 3655 Drummond St., Montreal, Quebec H3G 1Y6, Canada.
Received 13 August 1996; accepted in final form 2 July 1997.
 |
REFERENCES |
-
CANNON, S. C.,
ROBINSON, D. A.
Loss of the neural integrator of the oculomotor system from brain stem lesions in monkey.
J. Neurophysiol.
57: 1383-1409, 1987.[Abstract/Free Full Text]
-
CHERON, G.,
GILLES, P.,
GODAUX, E.
Lesions in the cat prepositus-vestibular complex: effects on the optokinetic system.
J. Physiol. (Lond.)
372: 95-111, 1986a.[Abstract]
-
CHERON, G.,
GODAUX, E.
Disabling of the oculomotor neural integrator by kainic acid injections in the prepositus-vestibular complex of the cat.
J. Physiol. (Lond.)
394: 267-290, 1987.[Abstract]
-
CHERON, G.,
GODAUX, E.,
LAUNE, J. M.,
VANDERKELEN, B.
Lesions in the cat prepositus-vestibular complex: effects on the vestibular ocular reflex and saccades.
J. Physiol. (Lond).
372: 75-94, 1986b.[Abstract]
-
CULLEN, K. E. AND GUITTON, D. IBN Activity Encodes Gaze, Not Eye, Metrics and Dynamics During Head on Body Rotations: Evidence That Vestibular Signals Supplement Visual Information in the Control of Gaze Shifts. New York: NY Acad. Sci., 1995a.
-
CULLEN, K. E.,
GUITTON, D.
Analysis of primate IBN spike train dynamics: gaze versus eye based models during combined eye-head gaze shifts.
Soc. Neurosci. Abstr.
21: 1915, 1995b.
-
CULLEN, K. E.,
GUITTON, D.
Analysis of primate IBN spike trains using system identification techniques. I. Influence of eye-movement dynamics during head-fixed saccades.
J. Neurophysiol.
78: 3259-3282, 1997a.[Abstract/Free Full Text]
-
CULLEN, K. E.,
GUITTON, D.
Analysis of primate inhibitory burst neuron spike trains using system identification techniques. III. Influence of motor error during head-fixed saccades and combined eye-head gaze shifts.
J. Neurophysiol.
78: 3307-3322, 1997b.[Abstract/Free Full Text]
-
CULLEN, K. E.,
GUITTON, D.,
JIANG, W.,
REY, C. G.
Gaze-related activity of putative inhibitory burst neurons in the head-free cat.
J. Neurophysiol.
70: 2678-2683, 1993.[Abstract/Free Full Text]
-
CULLEN, K. E.,
REY, C. G.,
GUITTON, D.
Dynamic models of primate and cat burst neuron activity during head-fixed saccades and head-free gaze shifts.
In: Contemporary Ocular Motor and Vestibular Research,
edited by
A. Fuchs,
T. Brandt,
U. Buttner,
and D. Zee
. Stuttgart, Germany: Thieme, 1994, p. 274-281
-
CULLEN, K. E.,
REY, C. G.,
GUITTON, D.,
GALIANA, H. L.
The use of system identification techniques in the analysis of oculomotor burst neuron spike train dynamics.
J. Comput. Neurosci.
3: 347-368, 1996.[Medline]
-
ESCUDERO, M.,
DE LA CRUZ, R. R.,
DELGADO-GARCIA, J. M. A
comparative neurophysiological study of prepositus hypoglossi and vestibular neurons projecting to the abducens nucleus in the alert cat.
J. Physiol. (Lond.)
458: 539-560, 1992.[Abstract]
-
EVINGER, C.,
KANEKO, C.R.S.,
FUCHS, A. F.
The activity of omnipause neurons in alert cats during saccadic eye movements and visual stimuli.
J. Neurophysiol.
47: 827-844, 1982.[Abstract/Free Full Text]
-
FUCHS, A. F.,
KANEKO, C.R.S.,
SCUDDER, C. A.
Brainstem control of saccadic eye movements.
Annu. Rev. Neurosci.
8: 307-337, 1985.[Medline]
-
FULLER, J. H.,
MALDONADO, H.,
SCHLAG, J.
Vestibular-oculomotor interaction in cat eye-head movements.
Exp. Brain Res.
271: 241-250, 1983.
-
FURUYA, M.,
MARKHAM, C.
Direct inhibitory synaptic linkage of pause neurons with burst inhibitory neurons.
Brain Res.
245: 139-143, 1982.[Medline]
-
GALIANA, H. L. AND GUITTON, D. Central organization and modeling of eye-head coordination during orienting gaze shifts. In: Sensing and Controlling Motion, edited by B. Cohen, D. L. Tomko, and F. Guedry. New York: NY Acad. Sci., 1992, p. 452-471.
-
GRANTYN, A.,
BERTHOZ, A.
Reticulo-spinal neurons participating in the control of synergic eye and head movements during orienting in the cat: I. Behavioral properties.
Exp. Brain Res.
66: 339-354, 1987.[Medline]
-
GUITTON, D.
Control of eye-head coordination during orienting gaze shifts.
Trends Neurosci.
15: 174-179, 1992.[Medline]
-
GUITTON, D.,
DOUGLAS, R. M.,
VOLLE, M.
Eye-head coordination in cats.
J. Neurophysiol.
52: 427-459, 1984.
-
GUITTON, D.,
MUNOZ, D. P.,
GALIANA, H. L.
Gaze control in the cat: studies and modeling of the coupling between orienting eye and head movements in different behavioral tasks.
J. Neurophysiol.
64: 509-531, 1990.[Abstract/Free Full Text]
-
GUITTON, D.,
VOLLE, M.
Gaze control in humans: eye-head coordination during orienting movements to targets within and beyond the oculomotor range.
J. Neurophysiol.
58: 427-459, 1987.[Abstract/Free Full Text]
-
HEPP, K.,
HENN, V.,
VILIS, T.,
COHEN, B.
Brainstem regions related to saccade generation.
In: The Neurobiology of Saccadic Eye Movements,
edited by
R. H. Wurtz,
and M. E. Goldberg
. Amsterdam: Elsevier/North-Holland, 1989, p. 105-211
-
HIKOSAKA, O.,
IGUSA, Y.,
IMAI, H.
Inhibitory connections of nystagmus-related reticular burst neurons with neurons in the abducens, prepositus hypoglossi and vestibular nuclei in the cat.
Exp. Brain Res.
39: 301-312, 1980.[Medline]
-
HIKOSAKA, O.,
IGUSA, Y.,
NAKAO, S.,
SHIMAZU, H.
Direct inhibitory synaptic linkage of pontomedullary reticular burst neurons with abducens motoneurons in the cat.
Exp. Brain Res.
33: 337-352, 1978.[Medline]
-
HIKOSAKA, O.,
KAWAKAMI, T.
Inhibitory neurons related to the quick phase of vestibular nystagmus their location and projection.
Exp. Brain Res.
27: 377-396, 1977.[Medline]
-
KELLER, E. L.
Participation of the medial pontine reticular formation in eye movement generation in monkey.
J. Neurophysiol.
37: 316-332, 1974.[Free Full Text]
-
KING, W. M.,
FUCHS, A. F.
Neural activity in the mesencephalon related to vertical eye movements.
In: Control of Gaze by Brain Stem Neurons,
edited by
R. Baker,
and A. Berthoz
. New York: Elsevier, 1977, p. 327-336
-
LAURUTIS, V. P.,
ROBINSON, D. A.
The vestibulo-ocular reflex during human saccadic eye movements.
J. Physiol. (Lond.)
373: 209-233, 1986.[Abstract]
-
LJUNG, L.
In: System Identification for the User.,
. Englewood Cliffs, NJ: Prentice-Hall, 1987
-
LOPEZ-BARNEO, J.,
DARLOT, C.,
BERTOZ, A.,
BAKER, R.
Neural activity in prepositus nucleus correlated with eye movements in the alert cat.
J. Neurophysiol.
47: 329-352, 1982.[Abstract/Free Full Text]
-
MCFARLAND, J. L.,
FUCHS, A. F.
Discharge patterns of nucleus prepositus hypoglossi and adjacent vestibular nucleus during horizontal eye movement in behaving macaques.
J. Neurophysiol.
68: 319-332, 1992.[Abstract/Free Full Text]
-
METTENS, P.,
GODAUX, E.,
CHERON, G.,
GALIANA, H. L.
Effect of muscimol microinjections into the prepositus hypoglossi and medial vestibular nuclei on cat eye movements.
J. Neurophysiol.
72: 785-802, 1994.[Abstract/Free Full Text]
-
NAKAO, S.,
CUTHROYS, I.,
MARKHAM, C.
Direct inhibitory projections of pause neurons to nystagmus-related pontomedullary reticular burst neurons in the cat.
Exp. Brain Res.
40: 283-293, 1980.[Medline]
-
PARÉ, M.,
GUITTON, D.
Gaze-related activity of brainstem omnipause neurons during combined eye-head gaze shifts in the alert cat.
Exp Brain Res.
83: 210-214, 1990.[Medline]
-
PÉLISSON, D.,
GUITTON, D.,
GOFFART, L.
On-line compensation of gaze shifts perturbed by micro-stimulation of the superior colliculus in the cat with unrestrained head.
Exp. Brain Res.
106: 196-204, 1995.[Medline]
-
PÉLISSON, D.,
GUITTON, D.,
MUNOZ, D. P.
Compensatory eye and head movements generated by the cat following stimulation-induced perturbations in gaze position.
Exp. Brain Res.
78: 654-658, 1989.[Medline]
-
PÉLISSON, D.,
PRABLANC, C.
Vestibulo-ocular reflex (VOR) induced by passive head rotation and goal-directed saccadic eye movements do not simply add in man.
Brain Res.
380: 397-400, 1986.[Medline]
-
PÉLISSON, D.,
PRABLANC, C.,
URQUIZAR, C.
Vestibulo-ocular reflex inhibition and gaze saccade control characteristics during eye-head orientation in humans.
J. Neurophysiol.
59: 997-1013, 1988.[Abstract/Free Full Text]
-
PHILLIPS, J. O. The Brainstem Control of Coordinated Eye-Head Movement: Neuroanatomical, Behavioral and Neurophysiological Studies of the Control of Gaze Behavior in the Monkey (PhD thesis). Seattle, WA: Univ. of Washington, 1993.
-
PHILLIPS, J. O.,
LING, L.,
FUCHS, A. F.,
SIEBOLD, C.,
PLORDE, J. J.
Rapid horizontal gaze movement in the monkey.
J. Neurophysiol.
73: 1632-1652, 1995.[Abstract/Free Full Text]
-
ROBINSON, D. A.
Oculomotor unit behavior in the monkey.
J. Neurophysiol.
33: 393-404, 1970.[Free Full Text]
-
ROUCOUX, A.,
GUITTON, D.,
CROMMELINCK, M.
Stimulation of the superior colliculus in the alert cat. II. Eye and head movements evoked when the head is unrestrained.
Exp. Brain Res.
39: 75-85, 1980.[Medline]
-
SCUDDER, C. A.,
FUCHS, A. F.
Physiological and behavioral indentification of vestibular nucleus neurons mediating the horizontal vestibuloocular reflex in trained rhesus monkeys.
J. Neurophysiol.
69: 244-264, 1992.
-
SCUDDER, C. A.,
FUCHS, A. F.,
LANGER, T. P.
Characteristics and functional identification of saccadic inhibitory burst neurons in the alert monkey.
J. Neurophysiol.
59: 1430-1454, 1988.[Abstract/Free Full Text]
-
STRASSMAN, A.,
HIGHSTEIN, S. M.,
MCCREA, R. A.
Anatomy and physiology of saccadic burst neurons in the alert squirrel monkey. II. Inhibitory burst neurons.
J. Comp. Neurol.
249: 358-380, 1986.[Medline]
-
TOMLINSON, D. R.
Combined eye-head gaze shifts in the primate. III. Contributions to the accuracy of gaze saccades.
J. Neurophysiol.
64: 1873-1891, 1990.[Abstract/Free Full Text]
-
TOMLINSON, R. D.,
BAHRA, P. S.
Combined eye-head gaze shifts in the primate. I. Metrics.
J. Neurophysiol.
56: 1542-1557, 1986a.[Abstract/Free Full Text]
-
TOMLINSON, R. D.,
BAHRA, P. S.
Combined eye-head gaze shifts in the primate. II. Interactions between saccades and the vestibuloocular reflex.
J. Neurophysiol.
56: 1558-1570, 1986b.[Abstract/Free Full Text]
-
TOMLINSON, D. R.,
BRANCE, M.
Brainstem control of coordinated eye-head gaze shifts.
In: The Head-Neck Sensory Motor System,
edited by
A. Bertoz,
P. P. Vidal,
and W. Graf
. Oxford, UK: Oxford Univ. Press, 1991, p. 356-361
-
VAN GISBERGEN, J.A.M.,
ROBINSON, D. A.,
GIELEN, S. A
quantitative analysis of generation of saccadic eye movements by burst neurons.
J. Neurophysiol.
45: 417-442, 1981.[Free Full Text]
-
WHITTINGTON, D. A.,
LESTIENNE, F.,
BIZZI, E.
Behavior of preocular burst neurons during eye-head coordination.
Exp. Brain Res.
55: 215-222, 1984.[Medline]
-
YOSHIDA, K.,
BERTHOZ, A.,
VIDAL, P. P.,
MCCREA, R. A.
Morphological and physiological characteristics of inhibitory burst neurons controlling rapid eye movements on the alert cat.
J. Neurophysiol.
48: 761-784, 1982.[Free Full Text]