Analysis of Primate IBN Spike Trains Using System Identification Techniques. I. Relationship to Eye Movement Dynamics During Head-Fixed Saccades

Kathleen E. Cullen and Daniel Guitton

Aerospace Medical Research Unit and the Montreal Neurological Institute, McGill University, Montreal, Quebec H3G 1Y6, Canada

    ABSTRACT
Abstract
Introduction
Methods
Results
Discussion
References

Cullen, Kathleen E. and Daniel Guitton. Analysis of primate IBN spike trains using system identification techniques. I. Relationship to eye movement dynamics during head-fixed saccades. J. Neurophysiol. 78: 3259-3282, 1997. The dynamic behavior of primate (Macaca fascicularis) inhibitory burst neurons (IBNs) during head-fixed saccades was analyzed by using system identification techniques. Neurons were categorized as IBNs on the basis of their anatomic location as well as by their activity during horizontal head-fixed saccadic and smooth pursuit eye movements and vestibular nystagmus. Each IBN's latency or "dynamic lead time" (td) was determined by shifting the unit discharge in time until an optimal fit to the firing rate frequency B(t) profile was obtained by using the simple model based on eye movement dynamics,B(t) = r b1Ė(t); where Ė is eye velocity. For the population of IBNs, the dynamic estimate of lead time provided a significantly lower value than a method that used the onset of the first spike. We then compared the relative abilities of different eye movement-based models to predict B(t) by using objective optimization algorithms. The most important terms for predicting B(t) were eye velocity gain (b1) and bias terms (r) mentioned above. The contributions of higher-order velocity, acceleration, and/or eye position terms to model fits were found to be negligible. The addition of a pole term [the time derivative of B(t)] in conjunction with an acceleration term significantly improved model fits to IBN spike trains, particularly when the firing rates at the beginning of each saccade [initial conditions (ICs)] were estimated as parameters. Such a model fit the data well (a fit comparable to a linear regression analysis with a R2 value of 0.5, or equivalently, a correlation coefficient of 0.74). A simplified version of this model [B(t) = rk + b1Ė(t)], which did not contain a pole term, but in which the bias term (rk) was estimated separately for each saccade, provided nearly equivalent fits of the data. However, models in which ICs or rks were estimated separately for each saccade contained too many parameters to be considered as useful models of IBN discharges. We discovered that estimated ICs and rks were correlated with saccade amplitude for the majority of short-lead IBNs (SLIBNs; 56%) and many long-lead IBNs (LLIBNs; 42%). This observation led us to construct a more simple model that included a term that was inversely related to the amplitude of the saccade, in addition to eye velocity and constant bias terms. Such a model better described neuron discharges than more complex models based on a third-order nonlinear function of eye velocity. Given the small number of parameters required by this model (only 3) and its ability to fit the data, we suggest that it provides the most valuable description of IBN discharges. This model emphasizes that the IBN discharges are dependent on saccade amplitude and implies further that a mechanism must exist, at the motoneuron (MN) level, to offset the effect of the bias and amplitude-dependent terms. In addition, we did not find a significant difference in the variance accounted for by any of the downstream models tested for SLIBNs versus LLIBNs. Therefore we conclude that the eye movement signals encoded dynamically by SLIBNs and LLIBNs are similar in nature. Put another way, SLIBNs are not closer, dynamically, to MNs than LLIBNs.

    INTRODUCTION
Abstract
Introduction
Methods
Results
Discussion
References

To generate a saccadic eye movement, motoneurons (MNs) are required to provide signals to the extraocular muscles that overcome the restraining forces of the orbital tissues and drive the eye quickly to its new position and hold it there. Robinson (1964) proposed that the following fourth-order differential equation provides an accurate description of MN discharge in this process
<IT>MN</IT>(<IT>t</IT>) = <IT>r + b</IT><SUB>0</SUB><IT>E + b</IT><SUB>1</SUB><IT><A><AC>E</AC><AC>˙</AC></A></IT>+ <IT>b</IT><SUB>2</SUB><IT><A><AC>E</AC><AC>¨</AC></A></IT>+ <IT>b</IT><SUB>3</SUB><IT><A><AC>E</AC><AC>&cjs1167;</AC></A></IT>+ <IT>b</IT><SUB>4</SUB><IT><A><AC>E</AC><AC>⃜</AC></A></IT>+ <IT>c<A><AC>M</AC><AC>˙</AC></A>N</IT>(<IT>t</IT>) (1)
where MN(t) is MN firing rate, E is the time-varying eye position, letters with an overdot are derivatives (e.g., Ė is eye velocity), and c, r, b0, b1, b2, b3, b4 are constants. The key predictions of this formulation were as follows: 1) MNs are required to generate a brief burst (a "pulse") of action potentials to overcome the viscous drag of the orbital tissues (b1Ė) to drive a saccade; 2) a lower frequency tonic discharge (a "step") is required to balance the elastic restoring forces of the orbital tissues to hold the eye steady in its new position (b0E); and 3) between the pulse and the step a slow decay or "slide" component exists, lasting until ~250 ms after the saccade, which just offsets the relaxation of the slow viscoelastic elements in the orbit [c&Mdot;N(t)]. Subsequent extracellular recordings have confirmed these predictions and have demonstrated the pulse-slide-step nature of the MN signal (Delgado-Garcia et al. 1986; Fuchs and Luschei 1970; Fuchs et al. 1988; Goldstein and Robinson 1984; Stahl and Simpson 1995).

A number of investigators have attempted to simplify Eq. 1 by retaining only the most significant terms. For example, Keller (1973) and Van Gisbergen et al. (1981) proposed a second-order model including eye position (E), velocity (Ė), and acceleration (Ë) terms. Goldstein and Robinson (1984) and Optican and Miles (1985) also proposed a second-order form but included as well a pole term [proportional to the derivative of the neuron's firing rate, &Mdot;N(t)] to explain the exponential decay of MN firing following a saccadic eye movement (the "postsaccadic slide"). More recently, Fuchs et al. (1988) demonstrated that a third-order model, with only inertia of the globe neglected in the derivation of the model Eq. 1, provides an adequate description of MN(t). However, Stahl and Simpson (1995) noted that the rejection of the second-order approximation of Eq. 1 may be premature because of questionable constraints placed on parameters by Fuchs et al. (1988).

The most prominent simplification to Eq. 1 is the first-order model (Robinson 1970)
<IT>MN</IT>(<IT>t</IT>) = <IT>r + b</IT><SUB>0</SUB><IT>E + b</IT><SUB>1</SUB><IT><A><AC>E</AC><AC>˙</AC></A></IT> (2)
In this "pulse-step" model, the burst from the MN pool (the dominant signal during saccades being proportional to eye velocity) drives the eye to the new position, whereas the tonic activity of the neuron holds the eye steady in its new position; there is no slide component [i.e., the model does not include the pole term c&Mdot;N(t)]. Although this model is not accurate, its wide-spread use prompted the hypothesis that the position and velocity terms are provided by separate inputs onto MNs, namely the integrator and burst generator signals, respectively. The finding that higher-order models are necessary for describing the downstream relationship between MN(t) and eye movement trajectory begs the question as to whether it is still useful to consider each term in Eq. 1 as a separate neural signal. Mounting evidence suggests that this is not the case. Indeed, neurons in the vestibular nucleus (VN) that project to MNs carry both position and velocity information and perhaps even the pole term (discussed in Fuchs et al. 1988). It is natural to ask whether other inputs to MNs also carry mixed signals. In this paper we consider the burst generator in this context.

For horizontal saccades, excitatory burst neurons (EBNs) project to the ipsilateral abducens nucleus (ABD) to activate it during saccades and act together with a group of inhibitory burst neurons (IBNs) that silence antagonist abducens MNs (Hikosaka and Kawakami 1977; Hikosaka et al. 1978; Igusa et al. 1980; Sasaki and Shimazu 1981; Scudder et al. 1988; Strassman et al. 1986a,b; Yoshida et al. 1982). IBNs and EBNs discharge vigorously or "burst" during ipsilateral saccadic eye movements and are thought to carry similar signals, because IBNs are driven by EBNs (Sasaki and Shimazu 1981; Strassman et al. 1986a) and EBNs appear to be a primary target of the contralateral projections of IBNs (Strassman et al. 1986b).

In recent models of the head-fixed saccadic system (Fuchs et al. 1985; Jürgens et al. 1981; Robinson 1975; Scudder 1988), a saccade is controlled by a feedback loop wherein an eye motor-error signal (desired change in eye position minus actual change in eye position) is the "upstream" signal that activates burst neurons (BNs) whereas the "downstream" output of these neurons drives MNs (Fig. 2B). From the downstream perspective and within the context of the simple model Eq. 2, BNs provide the eye velocity signal (for review see Fuchs et al. 1985). Accordingly, most prior analyses of BN activity have focused largely on either the relationship between peak eye velocity and peak discharge frequency or global characteristics of the discharge such as the relationship between number of spikes (NOS) and horizontal component of saccade amplitude (Kaneko and Fuchs 1981; Scudder et al. 1988; Strassman et al. 1986a,b; Yoshida et al. 1982). These are static relationships that do not provide information concerning the time-varying relationship between BN spike trains and resultant eye movement trajectories. Only the study by Van Gisbergen et al. (1981) considered time-varying or dynamic relationships. They asked whether the difference in the ON-direction (OND) and OFF-direction (OFFD) responses of BN signals could provide some of the terms in Eq. 1 other than velocity and examined, with phase-plane trajectories, the influence of acceleration. Using a method that required removing the acceleration term by trial and error, they suggested that the burst discharge is dynamically related to both eye velocity and acceleration.


View larger version (32K):
[in this window]
[in a new window]
 
FIG. 2. A: location of recorded IBNs was verified with the aid of a marking lesion in monkey L. Lesion was made in area located just dorsal to the region where putative IBN neurons were recorded. Dotted vertical line, path of electrode; filled circles, locations of cells previously recorded at the same location. B: classic local feedback model for saccade generation proposed by Jürgens et al. (1981) as a modification to Robinson's (1975) model. Delta E*, eye motor efference signal giving angular rotation of the eye that has occurred since onset of saccade is obtained by integrating burst neuron (BN; B) output that encodes eye velocity. Delta T, desired angular rotation of the eye. Difference between Delta T and Delta E* provides an estimate of the eye motor error [epsilon (t)] relative to desired change in eye rotation. Saccade is triggered when omnipause neurons (OPNs) stop firing, so their inhibitory influence on the BN pool is removed. Activity of B can be studied from two perspectives: influence of B on the eye movement trajectory (downstream analysis) and the response of B to eyemotor-error signal driving it (upstream analysis). epsilon , summing junction; E(t), resultant eye position.

This paper is the first in a series of three companion papers wherein we analyze in detail the discharges of one class of BNs, IBNs. We have focused on IBNs because 1) their anatomic location is easily defined, which facilitates microelectrode access; 2) they show an excellent static relationship to saccade metrics; and 3) they are driven by EBNs and EBNs appear to be a primary target of the contralateral projections of IBNs (Sasaki and Shimazu 1981; Strassman et al. 1986a,b). In this paper we employ objective computer-based optimization algorithms to characterize the discharges of BNs in the context of Eq. 1. Our objective is to demonstrate which of the terms in Eq. 1 are important for describing the discharge of this class of premotor neurons. Using algorithms founded in system identification theory (Ljung 1987)we sought to construct the mathematical model that best describes the input-output relationship between an IBN's firing rate and properties of the saccadic eye movement trajectory. In the second companion paper of this three paper series (Cullen and Guitton 1997a) we apply these techniques to the analysis of IBN discharges in the head-free monkey to discover whether burst discharges are best correlated to eye, head, or gaze trajectories. In the third companion paper of this series (Cullen and Guitton 1997b) we model IBN firing rate as a function of the upstream motor-error signal (Fig. 2B) and compare upstream and downstream models in terms of how well they account for the neural discharge pattern.

Our analysis method accomplished the following: 1) provided an objective method for calculating lead time and for comparing the goodness-of-fit of different models; 2) provided an indication of whether it is warranted to increase the goodness-of-fit by using more complex models that have more free parameters; and 3) permitted a quantitative evaluation of the relationship between estimated initial firing rates and/or biases and movement characteristics such as initial and final eye position, saccade amplitude, and peak eye velocity. The results of our analysis indicate that 1) contrary to the current assumption, primate short-lead IBNs (SLIBNs) and long-lead IBNs (LLIBNs) encode eye movement dynamics equally well, and 2) eye velocity-based models that included a bias term (r), which is inversely related to the amplitude of the saccade, provide a good description of BN spike trains.

    METHODS
Abstract
Introduction
Methods
Results
Discussion
References

Two monkeys (Macaca fascicularis) were prepared for chronic extracellular recording. The surgical preparation of the animals and the methods used for obtaining the extracellular recordings from brain stem neurons were similar to those previously described by Cullen and McCrea (1993). The IBNs described in this study were identified on the basis of their characteristic burst discharge during rapid orienting horizontal eye movements in the ipsilateral direction and their location in the IBN area caudal and ventral to the ABD (see descriptions for cat in Hikosaka and Kawakami 1977; Kaneko and Fuchs 1981; Yoshida et al. 1982; and monkey in Scudder et al. 1988). Figure 1 illustrates the anatomic distribution of the entire population of IBNs recorded in this study. The location of these neurons was determined by using reconstructions of electrode tracts, stereotaxic locations relative to the ABD, and later histological confirmation. Figure 2A illustrates the marking lesion made in monkey L that was made just dorsal to the location where two IBNs were recorded. Cells were also characterized as BNs on the basis of their burst responses during quick phases of vestibular nystagmus evoked during horizontal rotation of the vestibular turntable and also by their lack of response during slow phases of nystagmus and smooth pursuit eye movements.


View larger version (15K):
[in this window]
[in a new window]
 
FIG. 1. Reconstructed locations of the entire population of inhibitory burst neurons (IBNs) analyzed in this study are illustrated on stereotaxic cross sections of the macaque brain stem. Sections were drawn at 400-µm intervals. Filled circles, short-lead IBNs (SLIBNs); open circles, long-lead IBNs (LLIBNs); ABN, abducens nucleus; PHN, nucleus prepositus hypoglossi; nVII, 7th nerve.

Experimental paradigms and data collection

During the experiments the head-fixed monkey was seated in a primate chair that was fixed to the suprastructure of a vestibular turntable. For a juice reward the monkey was trained to orient to a target light that was projected onto a cylindrical screen surrounding the turntable. The spot of light was generated by a HeNe laser that was projected onto the screen via a system of two mirrors, mounted above the monkey's head, whose positions were controlled by galvanometers (General Scanning). The screen onto which the visual target was projected was monochromatic and evenly illuminated. The target was stepped between positions ±5, 10, 20, and 35° horizontal relative to the straight ahead position. The monkeys were also trained to orient to a food target that appeared unexpectedly on any side of an opaque screen facing the monkey (the "barrier paradigm," see Guitton et al. 1984 for details). This procedure provided us with a rich variety of saccade vectors, starting positions, and, for saccades of the same vector, a variety of velocity profiles.

On average, the saccades in this study were slower than those reported previously in the monkey (e.g., Fuchs 1967). Figure 3A illustrates the relationship between saccade duration and saccade amplitude for the data set obtained from one animal (monkey H). As expected, saccade duration was related to saccade amplitude; however, there was considerable variation across the data set. For example, the duration of a 30° saccade ranged from 60 to 160 ms. This variability was a result of the wide range of velocities that the monkey utilized to generate a saccade of a given amplitude (Fig. 3B). The peak eye velocity generated during a 30° saccade ranged from 250 to 800°/s. In general, the saccadic eye velocities in our study were highly variable and often similar to those that were reported in the literature for humans (Robinson 1964). The difference between the metrics of the saccades in our study and other primate studies may have arisen as a result of two factors: 1) the more natural behavioral paradigms employed in our study resulted in monkeys producing less "stereotyped" saccade trajectories to orient to targets (approx 70% of the saccades in each neuron's data set) and 2) M. fascicularis were used in this study, whereas most previous studies have used M. mulatta. Although these two types of monkeys have a similar oculomotor range, it is more difficult to obtain consistent behavioral performance in M. fascicularis monkeys (Cullen, unpublished observation).


View larger version (19K):
[in this window]
[in a new window]
 
FIG. 3. Relationships between (A) saccade duration and saccade amplitude and (B) peak saccade velocity and amplitude. Data points represent all the saccades that were analyzed for monkey H.

For this paper BN activity was recorded during rapid voluntary orienting eye movements (saccades) made with the animal's head held fixed. Data for the head-free condition are reported in Cullen and Guiton 1997a. A 386 microcomputer equipped with a data translation DT2825 data acquisition board, a DATAQ waveform scroller board, and a CED1401 peripheral device (Cambridge Electronic Design) was used for controlling visual and vestibular stimuli, for acquiring data, and for on-line display and analysis of data. Eye movements were recorded by using the magnetic search coil technique. The recorded position signals were low-pass filtered (250 Hz, 8 pole Bessel) and sampled at 1,000 Hz by the computer. The low-pass filtering shifted the recorded position signals backward in time by 2 ms; the resultant delay was taken into account in our subsequent analysis.

Extracellular recordings were obtained with enamel insulated tungsten microelectrodes (7-10 MOmega impedance, Frederick Haer), which were inserted into the cerebellum through a 22-gauge guide tube. The guide tube was positioned with an X-Y stage (Narishige) that was attached to the monkey's recording well. The microelectrode was advanced along the Z axis into the brain stem with a miniature hydraulic drive (Narishige). Spike potentials were amplified and filtered (band-pass 400 Hz to 10 kHz). A unit's isolation was continuously monitored, so that the level of a windowing circuit was set at an appropriate level to generate a pulse coincident with the rising phase of a spike. This pulse was sent to the event channel of the intelligent peripheral (CED1401) device, logging the event times of the neuronal action potentials (see Cullen and McCrea 1993).

Position signals and the single unit activity were simultaneously recorded on Digital Audio Tape (DAT) tape (20 kHz maximum sampling frequency). The DAT recorded data were often replayed off-line, so that the unit data could be carefully retriggered [this option was particularly important for the head-free recordings presented in the companion papers (Cullen and Guitton 1997a,b)].

Analysis of BN discharges

During off-line analysis sampled eye movement traces were digitally filtered (Matlab) at 125 Hz, because Fourier analysis of ocular saccades has revealed little power above 50 Hz (for review, see Cullen et al. 1994, 1996). Eye velocity was derived digitally from position data. A nonlinear function relating the NOS in the unit discharge to polar angle of the saccade was employed to determine the preferred direction of each IBN (Scudder et al. 1988). The onset and offset of a saccade were defined using a 20°/s eye velocity criterion. 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 following a saccade or gaze shift, burst offset was defined as when 95% of the spikes in the burst had occurred. A spike density function, in which a Gaussian function [5 (SD) ms] was convolved with the spike train, was used to represent the neural discharge. Choosing a Gaussian of this width effectively low-pass filtered the neural discharge so that its frequency content was similar to that of the coincident eye movements (see Cullen et al. 1996 for details). In a previous study (Cullen et al. 1996) we had determined that decreasing the width of the Gaussian by as much as 50% had little effect on the results of our dynamic analysis (latencies and parameters were not markedly changed). Only those saccades in which the amplitude of the vertical position component was less than one-third of the horizontal component were studied; i.e., the amplitude of the horizontal component was between 95 and 100% of the overall movement amplitude. We considered these to be horizontal saccades.

Dynamic models of BN firing rate

We used system identification techniques (Ljung 1987) for the analysis of IBN spike train dynamics. The use of these algorithms permitted optimization across multiple eye trajectories simultaneously and provided a method to objectively determine 1) how well a given model predicts a neuron's discharge pattern; 2) whether increasing the complexity of a model is permissible; 3) what a BN's dynamic latency is; and 4) the relationship between initial conditions (ICs) and quantities describes saccade trajectories. Details of these techniques can be found in Cullen et al. (1996) wherein we describe system identification techniques as they relate to BN analysis. IBN discharges were evaluated in terms of the currently favored schema (Fig. 2B) (Jürgens et al. 1981) of the saccade control system that states that burst firing rate is related simultaneously to the upstream motor-error signal (Cullen and Guitton 1997b) and to the downstream dynamics of the saccade trajectory (this paper; Cullen and Guitton 1997a). In RESULTS of this paper we compare different linear and nonlinear dynamic models that predict IBN firing rate on the basis of downstream signals during head-fixed saccades.

METHODS FOR PARAMETER ESTIMATION. For downstream models, BN discharge was the output and a function of parameters describing saccade trajectory. Therefore these were inverse models because burst frequency was the output of the system, rather than the input as it is in Fig. 2B. This analysis procedure was necessitated by the estimation techniques that required that the noise be in the output of the model formulation; this is conveniently placed in the BN discharge (Cullen et al. 1996). By using the IBN firing rate as the output, we were able to make comparisons between the fits computed for different plausible models (Tables 2 and 3). For each model optimal fits were made to an ensemble of ~40 saccades of different amplitudes between 5 and 45° to either the left or right, depending on the neuron's OND. We chose to characterize the discharges of IBNs in terms of horizontal eye movements, because IBNs project to the MNs that innervate the muscles that are responsible for producing horizontal eye movements (i.e., the lateral and medial recti). Note, in a preliminary analysis we analyzed the responses of several cells along their optimal response direction (see RESULTS) and found that goodness-of-fit, relative importance of estimated parameters, and the optimal latencies estimated in our models were not markedly changed. Generally, the starting positions of the saccades in each data set were well distributed over a range of ±20° on either side of the center position. Optimal fits to the neuron's discharge during OFFD saccades were carried out separately.

 
View this table:
[in this window] [in a new window]
 
TABLE 2. Parameter estimations for the downstream models tested

 
View this table:
[in this window] [in a new window]
 
TABLE 3. VAF and BIC values for downstream models tested

ESTIMATION OF DYNAMIC LATENCY. We shifted the burst in time relative to the saccade by the latency (the optimal dynamic delay td) that provided the best fit of the BN firing rate when a simple but accurate downstream model was used (see section on lead time in RESULTS). It is important to note that the entire portion of the temporally shifted spike train that was coextensive with the saccade duration was used to fit the model. Thus every spike within this interval contributed to the latency calculation. When the BN was shifted by the appropriate td, the main portion of the burst was typically aligned with the duration of the saccade. Consequently, the action potentials discharged by the neuron before the main portion of the burst (the burst "prelude") were typically rejected by the optimization algorithms and were not included in the fit. The estimated delay (td) was then applied to the other more complex models (i.e., it was a fixed parameter) (see Cullen et al. 1996 for details).

COMPARISON OF MODEL FITS. The goodness-of-fit of each model was determined by considering the root mean square (RMS) of the fit errors and the variance accounted for (VAF). The VAF was defined as [1 - (var/SD)] where var is the variance of the fit's error and SD is the standard deviation of the data set about its mean. Using VAF allowed us to normalize the goodness of a model's fit across neurons. It should be noted that the VAF for a linear model is equivalent to the square of the correlation coefficient (R; i.e., a VAF of 0.5 indicates that 50% of the variability in the unit's firing rate is explained by the model and would correspond to R = 0.71). We tried models of different complexities. A more complex model has more adjustable parameters and hence has the potential for providing a better fit to the data. The question then arises as to whether increasing model complexity is warranted. As described in Cullen et al. (1996) we answered this question objectively by using the Bayesian information criterion (BIC) (Schwarz 1978); a "cost index" that indicates whether increasing the complexity of a model is justified when the accompanying decrease in the error of the fit is taken into account. We computed the BIC for each model estimation; a decrease in the BIC value indicated that an increase in model complexity was warranted. If this index did not decrease when a more complex model was used, then the model was considered to be no better at describing the data than the more simple model.

    RESULTS
Abstract
Introduction
Methods
Results
Discussion
References

We present data from twenty-eight IBNs that were recorded during saccadic eye movements in two head-fixed monkeys. A characterization of the activity of the cells in terms of the trajectory of head-fixed saccades is presented in this paper. In our companion papers we 1) analyze the activity of these neurons in the head-free condition and compare the results with the head-fixed data (Cullen and Guitton 1997a) and 2) consider IBN discharges in terms of their response to a motor-error signal (Cullen and Guitton 1997b).

General discharge characteristics

Neurons were identified as IBNs on the basis of their anatomic location (see METHODS) and activity during saccadic and horizontal smooth pursuit eye movements and horizontal vestibular nystagmus. The neurons were identified on the basis of their characteristic burst discharge during horizontal saccadic eye movements in the ipsilateral direction. We determined that the cells were not the "burst-driver" neurons that were reported in a similar region in the cat (Ohki et al. 1988) by verifying that they were inactive during slow phases of vestibular nystagmus in both directions. It was also verified that these neurons were inactive during smooth pursuit eye movements and therefore that none could be categorized as the reticular formation smooth pursuit cells previously described by Tomlinson and Brance (1991). The neurons in this study were categorized as either SLIBNs (n = 16) or LLIBNs (n = 12), depending on whether the mean period between the onset of the first spike and the onset of eye velocity was <= 15 ms or >15 ms, respectively; we used this arbitrary definition to be compatible with the prior IBN study of Scudder et al. (1988).

The discharge of a typical SLIBN (H0415) is illustrated in Fig. 4, A-C. This neuron generated a strong burst of firing in the OND (ipsilateral) and weak burst in the OFFD (contralateral), with the first spike leading ipsilateral saccadic eye movements on average by 11 ± 2.8 (SD) ms (Fig. 4A). In Figs. 4 and 5 and all subsequent figures showing spike frequency, the spike discharge was shifted to the right by a time equal to the calculated dynamic lead time of each neuron (see Estimation of lead times). The duration of this neuron's burst was tightly correlated with the duration of the saccade (slope = 0.9; R = 0.91) and the total NOS during a burst was proportional to the horizontal component of the amplitude of the saccadic eye movement (slope = 1.5; R = 0.92). During smooth pursuit eye movements this cell only discharged during saccades and did not discharge during slow-phase smooth pursuit movements in either the ipsilateral or contralateral directions (Fig. 4B). During vestibular stimulation this cell discharged only during quick phases and did not discharge during ipsilateral or contralateral slow phases (Fig. 4C).


View larger version (23K):
[in this window]
[in a new window]
 
FIG. 4. Behavior of typical SLIBN (H0415; A-C) and LLIBN (H0421, D-F) during saccadic eye movements (A and D), smooth pursuit (B and E), and vestibularly induced eye movements (C and F). A strong burst occurred whenever saccade, or vestibular quick phase, was directed toward the ipsilateral side. Duration of each neuron's burst was correlated with the duration of the saccade (R = 0.91 and R = 0.74, respectively) and total number of spikes (NOS) generated during a burst was proportional to horizontal component of the amplitude of saccadic eye movement (R = 0.92 and R = 0.62, respectively). Note that the IBNs did not discharge in response to slow-phase smooth pursuit or vestibularly induced eye movements in either the ipsilateral or contralateral direction. Top trace in each panel, output of spike threshold discriminator. SD, spike density function showing firing frequency distribution; Ėh, horizontal eye velocity; Eh, horizontal eye position; Ev vertical eye position; Hh, horizontal head position; Th, horizontal target position.


View larger version (18K):
[in this window]
[in a new window]
 
FIG. 5. IBN firing rate behavior during ON-direction (OND) and OFF-direction (OFFD) saccades for 3 example IBNs; SLIBN L0702 and LLIBNs H0409 and H0925. A, B: units L0702 and H0409 were typical of 1/2 the cells in our population. For these IBNs there was a striking similarity in the shape of the neural firing frequency and eye velocity trajectories during horizontal saccades directed ipsilateral to the recording site (cell's OND), whereas the cells discharged few, if any, spikes during contralateral saccadic eye movements in the cell's OFFD. C: in contrast, the other 1/2 of cells in our population demonstrated a substantial OFFD response. OFFD response of unit H0925 was particularly robust and resembled an attenuated version of the OND response. Vertical line, 400 spikes/s, 400°/s, and 40°; SD, spike density function; Ė, horizontal eye velocity; E, horizontal eye position.

Figure 4, D-F, illustrates the discharge of a typical LLIBN (H0421). This neuron generated a preamble of lower frequency discharge followed by a burst of firing, with the first spike in the preamble leading ipsilateral saccadic eye movements by 21 ± 7.9 ms (Fig. 4D). The duration ofthis neuron's spike train was correlated with the duration of the saccade (slope = 0.9; R = 0.74) and the total NOS generated during the train was proportional to the horizontal component of the amplitude of the saccadic eye movement (slope = 0.9; R = 0.62). As was the case for SLIBNs, this neuron did not discharge in response to slow-phase smooth pursuit or vestibularly induced eye movements in either the ipsilateral or contralateral directions (Fig. 4, E and F).

IBN preferred direction

The IBNs in this study burst vigorously for ipsilaterally directed saccades, and they typically discharged far fewer spikes during contralateral (OFFD) and pure vertical saccades. Figure 5 shows examples of discharges of three neurons we selected to illustrate our data throughout this paper and the companion papers (Cullen and Guitton 1997a,b): SLIBN L0702, LLIBN H0409, and LLIBN H0925. The latter neuron had the best OFFD discharge in our population. Note the clear correlation between burst frequency and saccade velocity for discharges in the ONDs and OFFDs for H0925 and in the OND for cells L0702 and H0409. A precise measure of a neuron's optimal saccade direction was calculated by fitting a nonlinear function to relate the NOS in a unit's discharge to the polar angle of the saccade (see METHODS). The saccade amplitudes used to determine the preferred direction of each cell in the study ranged from 15 to 25°. The preferred directions in our population (Fig. 6) varied from +29° (upward) to -30° (downward). The mean preferred direction was -1.0 ± 15.2° closely aligned with the horizontal stereotaxic plane (~15° upward from the plane of the horizontal canals). Our analysis of eye movements that were horizontal, or nearly so, is consistent with the mean preferred horizontal direction of IBNs (see METHODS).


View larger version (33K):
[in this window]
[in a new window]
 
FIG. 6. Preferred direction of each of the IBNs in our population is illustrated in polar coordinates; length of vector indicates sensitivity of each IBN measured in spikes/deg, according to the scale on left. Mean preferred direction of the entire sample of IBNs was -1° (range +29 to -30°) suggesting that the direction of the population response of these cells was very closely aligned with ipsilateral horizontal. Average population vector is represented by the heavy arrow. Saccades included in this analysis varied in amplitude over a range of 15-25°. Each unit's optimal direction was determined by fitting the data with the sum of 2 sinusoids by using a nonlinear least squares algorithm. Inset: direction tuning of an example IBN (SLIBN H0409) is illustrated by a plot of the NOS vs. saccade direction.

Estimation of lead times

The time by which a BN's discharge led eye saccades was calculated by two methods. Histograms illustrating the results of applying each method are shown in Fig. 7. In Fig. 7A the lead time was determined by calculating the period between the onset of the first spike and the onset of eye velocity. This method was used in many previous studies of IBNs and as a criterion to categorize the cells as either SLIBNs or LLIBNs (e.g., Scudder et al. 1988). In the second method (Fig. 7B), by using our optimization algorithms described in detail in Cullen et al. (1996), a neuron's dynamic lead time (td) was left as an unknown and determined by shifting progressively the unit discharge in time relative to the saccade until a best fit was obtained for models 2d and 6d (Table 2)
model 2d:<IT>B</IT>(<IT>t</IT>) = <IT>r + b</IT><SUB>1</SUB><IT><A><AC>E</AC><AC>˙</AC></A></IT>(<IT>t − t</IT><SUB><IT>d</IT></SUB>) (3)
and
model 6d:<IT>B</IT>(<IT>t</IT>) = <IT>r + b</IT><SUB>1</SUB><IT><A><AC>E</AC><AC>˙</AC></A></IT>(<IT>t − t</IT><SUB><IT>d</IT></SUB>) + <IT>b</IT><SUB>2</SUB><IT><A><AC>E</AC><AC>¨</AC></A></IT>(<IT>t − t</IT><SUB><IT>d</IT></SUB>) − <IT>c<A><AC>B</AC><AC>˙</AC></A></IT>(<IT>t</IT>) (4)
where r, the bias term; b1, the velocity gain; b2, the acceleration gain; c, the coefficient of the pole; and td are all varied until RMS error is minimized for the whole data set. These equations are accurate dynamic models of firing rate, model 2d being the simplest and 6d the most complex to analyze. Model 2d has three parameters to optimize, whereas 6d, having the pole term &Bdot;(t), has 44 (when 40 saccades are analyzed). For two neurons in which we studied lead time in detail, Eq. 3 provided a clear minimum error for a particular value of td. With Eq. 4 it was more difficult to optimize the value of td, because the influence of this quantity in the optimization process was diluted by the larger number of parameters (Cullen et al. 1996). Model 6d was less sharply "latency tuned" than model 2d. Nevertheless, we found no difference between the lead time estimates given by Eqs. 3 and 4. Accordingly, for all neurons we will present data obtained from Eq. 3. Figure 7B shows that the dynamic method provides latency estimates that have far less scatter than those obtained by using the onset of the first spike(Fig. 7A).


View larger version (16K):
[in this window]
[in a new window]
 
FIG. 7. Histograms illustrating lead time by which IBN discharges led head-fixed saccades. Lead time was obtained by 2 methods: (top) time between onset of 1st spike and onset of eye velocity; (bottom) dynamic lead time td obtained by shifting unit discharge in time until optimal fit was obtained for the model, B(t) = r + b1Ė(t - td). On average, the burst of monkey IBNs led ipsilateral head-fixed saccades by 17.0 ms (SD = 9.3) and 13.4 ms (SD = 3.1) by using onset of 1st-spike and dynamic estimation methods, respectively.

The relationship between the lead times calculated by using each method is plotted in Fig. 8. The values obtained from the two methods were not well correlated for the 28 primate neurons we analyzed (R = 0.51). However this relationship improved once a single outlier (the value within parentheses) was removed from the analysis (R = 0.67). A continuity can be seen between the dynamic lead times estimated for SLIBNs (black-diamond ) and LLIBNs (diamond ). The time-of-first-spike method yielded a mean lead time for SLIBNs that was statistically similar to the one using dynamic analysis with Eq. 3 (11.6 ± 3.8 ms vs. 11.8 ± 2.7 ms, respectively). However, for LLIBNs the first-spike method yielded significantly longer lead times (24.2 ± 9.7 ms vs. 15.5 ± 2.2 ms, respectively). For the entire population of neurons, including SLIBNs and LLIBNs, the mean lead times were 17.0 ± 9.3 ms and 13.4 ± 3.1 ms for the first-spike and dynamic methods, respectively. The difference between these values was significantly different (Student's t-test,P < 0.01). It is noteworthy that the mean lead time estimated dynamically for neurons classified as SLIBNs was significantly less (P < 0.001) than that estimated, using the same method, for neurons classified as LLIBNs [11.8 ± 2.7 ms vs. 15.5 ± 2.2 ms, respectively, a difference of 3.7 ms to be compared with the 8 ms estimated by Scudder et al. (1988)]. Only five neurons (4 LLIBNs and 1 SLIBN) had dynamic lead times longer than 15 ms.


View larger version (21K):
[in this window]
[in a new window]
 
FIG. 8. Relationship between lead time values obtained by 2 separate methods that were used in this study: the onset of 1st-spike and dynamic estimation method. The 2 lead time measures were not well correlated for total sample of IBNs (n = 28; R = 0.51). However, this relationship improved considerably (R = 0.67) when a single outlier (diamond ) was removed from the analysis. This analysis shows that SLIBNs (black-diamond ) tended to have shorter dynamic lead times than LLIBNs (diamond ). Distinction between SLIBNs and LLIBNs was made according to relative timing of 1st spike to onset of saccade.

Relationships between IBN activity and saccade metrics

Linear relationships between the duration of IBN discharges and ipsilaterally directed saccadic eye movements were reported in a number of previous studies and our data are comparable. Figure 9, A and B, illustrates the relationship between burst duration and saccade duration for neurons classified as SLIBNs and LLIBNs, respectively. The insets show data obtained from neurons that we selected as examples and whose data are illustrated throughout the paper. The heavy lines show the population means whose slopes, intercepts, and correlation coefficients are given in Table 1. The slopes for SLIBNs and LLIBNs were not significantly different (P = 0.11), but the intercepts were (P < 0.05). In Table 1 we also show data from Scudder et al. (1988). Both data sets concur; the average burst duration of SLIBNs closely approximates the duration of the associated saccade, whereas for LLIBNs the discharge preamble produces spike trains that are longer than saccades. The present analysis suggests that the duration of this preamble, for a given cell, is constant (because the slope is ~1) and independent of saccade duration (and amplitude).


View larger version (30K):
[in this window]
[in a new window]
 
FIG. 9. A: relationship between burst duration and saccade duration in the head-fixed condition for population of SLIBNs. Heavy line, average regression fit of this relationship (slope = 0.87, intercept = 10.6). Inset: data for our example SLIBN (unit L0702). B: relationship between burst duration and saccade duration in the head-fixed condition for population of LLIBNs. Heavy line, average regression fit of this relationship for LLIBNs (slope = 0.97, intercept = 26.5). Inset: data for our example LLIBN (unit H0409).

 
View this table:
[in this window] [in a new window]
 
TABLE 1. Linear correlations of IBN discharges with eye saccade metrics

A number of studies have reported a strong relationship between the NOS in an IBN's discharge and the horizontal amplitude of ipsilaterally directed saccades. To obtain the NOS in the portion of the burst that directly drives saccades, we shifted the spike train in time by the dynamic lead time (td, calculated by using model 2d of Table 2) and counted the NOS during the discharge period that was aligned with each saccade (e.g., between the vertical dashed lines in Fig. 4A). The NOS was proportional to the amplitude of the horizontal component of the saccadic eye movement for each neuron. (Because the saccades we analyzed were nearly horizontal, using the overall amplitude did not influence this conclusion.) This relationship is illustrated for SLIBNs and LLIBNs in Fig. 10, A and B, respectively. The heavy lines show the population means whose slopes, intercepts, and correlation coefficients are shown in Table 1, as well as the equations given by Scudder et al. (1988). The slope and intercept of the mean lines were not significantly different for SLIBNs and LLIBNs (P = 0.41 and 0.11, respectively) and are also similar to Scudder et al. (1988). We conclude that our method of defining the interval over which NOS were counted, by using an objective estimate of lead time, preserves the strong relationship that was observed in previous studies between total NOS in a burst and the horizontal amplitude of a saccade.


View larger version (31K):
[in this window]
[in a new window]
 
FIG. 10. A: relationship between total NOS in a discharge and saccade amplitude in the head-fixed condition for population of SLIBNs. Heavy line, average regression fit of this relationship (slope = 1.1, intercept = 9.3). Inset: data for our example SLIBN (unit L0702). B: relationship between total NOS and saccade amplitude in the head-fixed condition for population of LLIBNs. Heavy line, average regression fit of this relationship for LLIBNs (slope = 0.95, intercept = 11.9). Inset: data for example LLIBN (unit H0409).

Previous studies also reported a relationship between IBN peak firing frequency and peak saccade velocity (Scudder et al. 1988; Yoshida et al. 1982). We did not find this relationship for SLIBNs and LLIBNs (Table 1; mean R = 0.39) to be as robust as that in Scudder et al. (1988); it was significant for 38 and 8% of our SLIBNs and LLIBNs, respectively. This difference may be due to the less-stereotyped saccade profiles that were elicited during some of our protocols (see METHODS and DISCUSSION).

Dynamic models of IBN discharges during head-fixed saccades in OND

The foregoing analysis of IBN lead times and relationships with saccade metrics has 1) characterized our cells as typical IBNs by using criteria similar to those of other studies and 2) introduced and employed a dynamic estimate of discharge frequency lead time, which is further utilized in the following text to evaluate different models that relate IBN spike train dynamics to eye movement parameters during OND saccades. Note that the action potentials discharged by the neuron before the main portion of the burst (the burst prelude) and after (the tail) were typically excluded by the optimization algorithms (see METHODS). The resulting fits are illustrated throughout the paper for the example SLIBN (L0702) and LLIBN (H0409).

Table 2 shows the eight downstream models we investigated. It also provides the means and SDs of the parameters estimated for each model for SLIBNs, LLIBNs, and the entire population of IBNs. Table 3 gives the VAF and BIC produced by each of these models for each of the 28 neurons in this study as well as the means and SDs for SLIBNs, LLIBNs, and for the entire population.

We will now consider each model in turn beginning with the simplest.

BURST FREQUENCY PROPORTIONAL TO EYE VELOCITY. Model 1d in Table 2 is the simplest. This often used model addresses the question of whether instantaneous burst frequency is proportional to eye velocity. It is based on the well-known relationship (discussed previously) between peak discharge frequency and peak eye velocity (Kaneko and Fuchs 1981; Scudder et al. 1988; Strassman et al. 1986a,b; Yoshida et al. 1982).

The VAF and BIC values for model 1d are shown in Table 3. This model was a poor predictor of the data as indicated by the low VAFs. Indeed, a negative VAF in 12 of 28 of the neurons implies that the variance of this model's fit error was actually larger than the variance of the data about a constant mean level.

ADDITION OF BIAS TERM. Model 2d is an extended form of 1d and can be viewed as a first-order linear approximation to a nonlinear relationship between B(t) and Ė where r is the bias term.

When this bias term was added to the model (Table 2), the VAF for the entire population increased very significantly from 0.04 to 0.30 (Table 3). Note that this model actually provided a fit to the data that was on average as good as that described by a correlation coefficient of 0.55 in a linear regression analysis. The accompanying decrease in the average BIC value (from 10.23 to 9.60) indicates that the addition of a bias term was warranted.

Averaged measures of goodness-of-fit were also calculated separately for SLIBNs and LLIBNs for the models 1d [Table 3: VAF (short-lead) = 0.05; VAF (long-lead) = 0.03] and 2d [Table 3: VAF (short-lead) = 0.32; VAF (long-lead) = 0.28]. There was no significant difference in the VAF by each of these models nor by any other model we tested for the two subpopulations of cells; accordingly in Table 2 we pooled the results for both SLIBNs and LLIBNS and display the average value for each parameter in each model. The considerable improvement in the fit of IBN firing rate that occurred with the addition of the bias term is illustrated for three example saccades for SLIBN L0702 (Fig. 11A) and LLIBN H0409 (Fig. 11B). The goodness-of-fit to the SLIBN is particularly striking.


View larger version (46K):
[in this window]
[in a new window]
 
FIG. 11. Examples of model fits to IBN firing frequency profiles during head-fixed saccades for our 2 example IBNs illustrated in Figs. 5, 9, and 10, SLIBN L0702 (A) and LLIBN H0409 (B). The most simple downstream model tested [top row: B(t) = b1Ė(t - td)] generally provided a poor fit of IBN discharge during saccades. However, addition of bias term to the model greatly improved our ability to estimate IBN discharges (2nd row). In contrast, addition of acceleration term to this eye velocity-based model (3rd row) was not very effective in further decreasing model error of the fits to either short-lead BNs (SLBNs) or long-lead BNs (LLBNs) discharges. Values of parameters, estimated by using >40 saccades, are provided below each model fit. Three example saccades are illustrated for each IBN where shaded trajectories represent unit activity and superimposed curves represent the fits produced by each model. Bottom 2 traces: accompanying eye velocity and position trajectories that were shifted in time by estimated optimal dynamic latency td.

IMPORTANCE OF ACCELERATION TERM. The next downstream model we investigated utilized a simple approximation to the hypothesis of Van Gisbergen et al. (1981), that unit firing can be represented by a nonlinear function of eye velocity plus an eye acceleration term. As a first approximation to this proposal, the nonlinearity in velocity was (as in Eq. 3) simply represented by the bias term (r; Table 2: model 3d). For all our neurons the addition of the acceleration term (b2) had little effect on increasing the VAF (or decreasing the BIC); the addition of an acceleration term increased the mean VAF from 0.30 to 0.31 (Table 3). The same was true when the SLIBNs and LLIBNs were considered separately. Examples of the marginal improvement in the fit of the IBN firing rate that accompanied the addition of an acceleration term are illustrated in Fig. 11.

NONLINEAR VELOCITY TERMS. A more complex nonlinear function of eye velocity was tested that included a third-order nonlinearity in addition to an eye acceleration (Table 2: model 4d). The addition of the higher-order eye velocity terms (d1 and d2) decreased the error of the model slightly (Table 3: mean population VAF = 0.33). Taken together the results so far show that in downstream models a bias term is much more significant than either acceleration or nonlinear velocity terms for representing IBN discharges.

ESTIMATION OF POLE TERM. A model was also investigated which in addition to eye velocity, acceleration and bias terms allowed for the estimation of a pole term in the system (the derivative of the IBN firing rate; Table 2: models 5d and 6d). This corresponds to the slide component of the MN's response described in Eq. 1 by the term c&Mdot;N(t). The addition of this term provided a slightly better model fit than model 3d when the ICs in firing rate (necessitated by the pole term) were taken from the data (Table 3: VAF = 0.31 in model 3d vs. 0.37 in model 5d). However, when the ICs were estimated as parameters (Table 2: model 6d), the fit improved considerably (Table 3: VAF = 0.37 in 5d vs. 0.54 in 6d). For model 6d the fits to the data were as good, on average, as those described in a linear regression analysis by R values of 0.74 and 0.73 for SLIBNs (VAF = 0.55) and LLIBNs (VAF = 0.53), respectively. Examples of fits to firing rate are shown in Fig. 12 for the same example neurons and saccades as in Fig. 11. The improvement in fit when ICs were estimated as parameters (2nd panel from top) rather than taken from the data set (top panel) can be appreciated by visual inspection. In Fig. 12 the numbers above each firing frequency profile give the initial frequency B(t0).


View larger version (41K):
[in this window]
[in a new window]
 
FIG. 12. Model fits are illustrated for example SLIBN L0702 (A) and LLIBN H0409 (B) for the same saccades illustrated in Fig. 11. A model, which included eye velocity, eye acceleration, bias, and pole terms, produced particularly good fits of IBN discharges when firing rate at beginning of each saccade was estimated as a parameter (2nd row) rather than taken from data (top row). A model, without a pole term, for which bias term was estimated separately for each saccade (3rd row), provided fits that were nearly comparable to those of row 2. Parameters estimated from >40 saccades are given below each model fit. Values B(t0), listed separately for each of 3 saccades, represent actual values of firing rate at beginning of saccade (top row) and the values estimated by algorithm for initial firing rate when B(t0) was included as parameter in the model (2nd row). Values of rk listed separately for each of 3 saccades in the 3rd panel represent bias estimated for each saccade. Top 3 rows: shaded area shows actual firing frequency profiles; heavy solid line, estimated firing. Bottom 2 traces: accompanying eye velocity and position trajectories that were temporally shifted by estimated optimal dynamic latency td.

A model of the form 6d requires the estimation of many more parameters (equivalent to the number of saccades included in the analysis) than one in which the ICs were taken directly from the data (model 5d). To determine whether increasing the model complexity to this extent (from 4 to 44 parameters) was warranted, it was necessary to consider the BIC. Comparison of the BIC values for models 3d, 5d, and 6d (Table 3) shows that the addition of ICs as parameters to be evaluated was warranted (mean BIC = 9.57, 9.41, and 8.84, respectively). Model 6d, however, suffered from unrealistic values of the bias term and this will be considered in the DISCUSSION.

ESTIMATION OF VARIABLE BIAS. In model 6d, the average time constant of the decaying ICs is long relative to saccade duration and given by the value c = 2.4 s (see Table 2; DISCUSSION). Thus the ICs were effectively constants where the model estimated a different constant value for each saccade that was the sum of the estimated value of the bias and the constant IC. This observation led us to specifically test a model in which only one bias term was estimated separately for each saccade (Table 2: model 7d). This variable bias model had nearly as many parameters as did model 6d (42 vs. 44, respectively); however, the method of estimation was simpler for a model of this form than for models that include a pole term (see Cullen et al. 1996). Example fits obtained using model 7d are illustrated in Fig. 12 (3rd panel from top). This model provided fits that were nearly as good as those given by model 6d (Fig. 12; Table 3: VAF = 0.54 for model 6d vs. 0.52 for model 7d). The possible physiological implications of the variable bias are considered in the following section and in DISCUSSION.

ICs AND BIASES DEPEND ON SACCADE PARAMETERS. For 14 cells---the majority of SLIBNs (9/16) and many LLIBNs (5/12)---the values estimated for ICs in model 6d were inversely and significantly correlated with the metrics of saccades. This was not the case for the example neurons shown in Figs. 11 and 12, and so these relationships are illustrated for another LLIBN (H0925) in Fig. 13, A and B. The estimated ICs for this cell were well correlated with saccade amplitude (R = -0.86) and less strongly related to peak eye velocity (R = -0.62). For those neurons with significant correlations, ICs were better correlated with saccade amplitude (mean SLIBNs, R = -0.70; LLIBNs, R = -0.67) than with peak saccade velocity (mean SLIBNs,R = -0.63; LLIBNs, R = -0.59). Note that in all these neurons the value of the ICs varied inversely with amplitude and peak velocity (as in Fig. 13).


View larger version (31K):
[in this window]
[in a new window]
 
FIG. 13. Analysis of relationship between estimated initial conditions (ICs) and biases for LLIBN H0925. A and C: for this cell both ICs [IC = B(t0) + r; model 6d] and estimated biases (rk; model 7d) were negatively correlated with amplitude of the saccade (R = -0.86 and -0.89, respectively). In addition, estimated ICs and biases were also correlated with peak velocity of saccades for this neuron (B, D). However, this relationship was less robust. Insets: example of difference in the firing rate profiles generated by LLIBN H0925 during small vs. large saccades. Shaded area, actual firing frequency; superimposed solid line, estimated firing frequency using model 7d. Note that peak heights of burst firing frequency profiles are comparable in spite of the nearly 10-fold difference in amplitude between the 2 saccades. Best fit, by using model 7d, estimates a bias of approx 300 spikes/s for smaller saccade and approx 25 spikes/s for larger saccade, thereby producing a good fit of the unit discharge during wide range of saccade amplitudes.

For 9 of the 14 cells whose ICs were correlated to metrics, the bias term estimated for model 7d (Fig. 13, C and D) was also linearly correlated with the amplitude and/or peak velocity of the saccade. The estimated biases for cell H0925 were well correlated with saccade amplitude (R = -0.89) and less strongly related to peak eye velocity (R = -0.56). The bias term was largest for small saccades. This implies that model 2d, which has a fixed bias, averaged across all saccade amplitudes should be best at predicting medium-amplitude saccades. Furthermore, model 2d should underevaluate the discharge frequency of small saccades and overevaluate that of large ones. This is precisely what Fig. 14 shows for neuron H0925.


View larger version (28K):
[in this window]
[in a new window]
 
FIG. 14. Downstream model containing eye velocity and bias term (model 2d) was estimated to predict firing behavior of LLIBN H0925 during 40 saccades ranging in amplitude from 5 to 45° [variance accounted for (VAF) = 0.32]. Globally estimated model was then applied to 3 distinct subsets of original data set, each of which included only saccades within a specific amplitude range. In general, the model underestimated discharge of small saccades ranging from 5 to 10° (A) and overestimated discharge of large saccades ranging from 25 to 30° (C). The model fit was best for subset of saccades whose amplitudes were in the middle range (15-20°) of saccades contained in original data set (B). Accordingly, VAF when globally estimated model 2d was applied to subset of mid-range saccades was 0.49% (B), but only 0.15 and 0.25% when applied to the subsets of small and large saccades (A and C, respectively). Results of this analysis of LLIBN H0925 and that illustrated in Fig. 13 indicate that the downstream model in which bias term could be estimated separately for each saccade (model 7d) would better describe discharge of this neuron. For A-C, shaded area shows actual firing frequency, superimposed solid line represents estimated frequency; 2nd and 3rd traces, accompanying eye velocity and position trajectories temporally shifted by estimated dynamic latency td.

We also investigated whether a relationship existed among either the estimated biases and/or ICs and the position of the eye before and/or following a saccade. For two SLIBNs there was a significant relationship among the estimated ICs, as well as the estimated biases, and the position of the eye before the onset of the saccade; however, no relationship was observed with final position for any IBN in this study. No relationships were ever observed among the ICs when they were taken directly from the data (model 5d) and any of the measured saccade parameters (amplitude, peak velocity, or initial or final eye position).

EFFECT OF SACCADE AMPLITUDE ON IBN SPIKE TRAIN DYNAMICS. Our finding that estimated ICs and biases were often correlated with saccade amplitude led us to construct a more simple model that included an amplitude-dependent term. This model (Table 2: model 8d) is a simple extension of model 2d, to which a term that is proportional to saccade amplitude (r1Delta E) was added. The addition of the eye amplitude term to model 2d resulted in comparable improvements in VAFs of LLIBNs (Table 3: VAF = 0.28 in model 2d vs. 0.32 in model 8d) and SLIBNs (Table 3: VAF = 0.32 in 2d vs. 0.36 in 8d). On average, model 8d was slightly better at describing IBN discharges than the nonlinear function of eye velocity in model 4d (Table 3: VAF = 0.34 in 8d vs. 0.33 in 4d). However, model 8d did not fare as well as 7d, but it had 38 less parameters.

EFFECT OF EYE POSITION. The final two models that were tested addressed the issue of whether eye position was relevant in estimating IBN discharges. The simplest model included two terms, an eye position term and a bias. For each neuron this model was far less effective at describing the discharges of cells than the simple model of Eq. 3 that included an eye velocity term and a bias (VAF = 0.05 vs. 0.30 in model 2d). Furthermore, the addition of a position term to model 2d resulted in only a negligible increase in the VAF (VAF = 0.31).

Downstream model errors for SLIBNs vs. LLIBNs

We have shown that, as a population, SLIBNs have a shorter dynamic lead time than LLIBNs. This is compatible with a linear signal processing stream wherein it is possible that LLIBNs project to SLIBNs that in turn project to MNs (Fuchs et al. 1985). However, we argue here and in the companion paper (Cullen and Guitton 1997b) that this is not the case. To determine whether SLIBNs encode downstream-eye movement dynamics better than LLIBNs, we examined the relationship between lead time and the VAF values estimated for each downstream model. The VAF was not lower for those neurons that had shorter lead times; hence, SLIBNs do not appear closer than LLIBNs to the motor output for any of the eight models that we tested. The relationship between lead time and VAF for model 8d is illustrated in Fig. 15.


View larger version (14K):
[in this window]
[in a new window]
 
FIG. 15. Relationship between IBN lead times and VAF provided by downstream model 8d (Table 2). There was no clear relationship among lead times calculated by using either onset of 1st spike (triangle ) or the estimated dynamic lead time (bullet ) and VAF provided by this model or any other downstream model that was tested. This shows that there is no distinction between SLIBNs and LLIBNs regarding their link to motor output.

OFFD discharges

The OFFD discharge of IBNs is of considerable interest because it will subtract from the excitatory drive provided by EBNs on agonist MNs. It is therefore important to provide a quantitative evaluation of this effect if we are to understand how MN signals are generated (Van Gisbergen et al. 1981).

Figure 5 shows the OFFD of the three neurons that we used extensively as examples in this paper. The OFFD of SLIBN L0702 and LLIBN H0409 consisted of an occasional one or two spikes and could not be analyzed. It is important to note that 50% of each of our SLIBN and LLIBN populations gave similar weak OFFD discharges. The other 50% provided OFFD discharges that varied in the quality of their relationships to eye dynamics; of these, cell LLIBN H0925 exhibited the best model fits, as can be appreciated from Fig. 5. We will consider this neuron in detail.

For cell H0925, the dynamic lead time in the OFFD was 11 ms compared with 17 ms in the OND. The OND and OFFD data sets contained a similar distribution of saccade sizes. Consequently, the difference in OND and OFFD dynamic lead times does not appear to be the result of a difference in the amplitude of the saccades that were included in the analysis. As for the correlation coefficient in the linear relation between NOS and amplitude, the values were 0.86 and -0.71 in the ONDs and OFFDs, respectively. Table 4 compares for this neuron the characteristics of the model fits for each model we have tested. The VAFs in the OND and OFFD were generally comparable. The comments made in an earlier section on the relative accuracy of each model, as model characteristics and complexity were changed, apply to the OFFD. Thus models 6d and 7d best describe the OFFD response of cell H0925. Indeed, in some cases the OFFD fits were superior to the OND fits, as exemplified by model 6d where the VAF was 0.57 for the OFFD compared with 0.50 for the OND. The OFFD fit in this case is equivalent to a linear regression coefficient (R) of ~0.76.

 
View this table:
[in this window] [in a new window]
 
TABLE 4. ON- vs. OFF-direction response (unit H0925): parameter estimations for the downstream models

To compare parameters in the OND and OFFD models we will focus on models 2d, 7d, and 8d because of their similarity in structure and also because we later propose that 8d is the simplest, most accurate and physiologically plausible representation of the cell's OND discharges (see DISCUSSION). In model 2d the absolute value of the velocity gain, b1 = -0.83, was ~10% lower in the OFFD compared with the OND and, most strikingly, the bias was far less and essentially zero in the OFFD. This latter point is emphasized by the same VAF for the OFFD in models 1d and 2d. Similarly in model 8d, the absolute value of the velocity gain was ~20% lower in the OFFD (b1 = -0.86) compared with the OND, and the bias was nearly zero in the OFFD. The eye amplitude dependent gain factor (r1) had considerably less effect on the model fit in the OFFD (VAF = 0.33 in model 8d vs. 0.32 in model 2d) than in the OND (VAF = 0.45 in 8d vs. 0.32 in 2d). Accordingly, the amplitude gain r1 was considerably less in the OFFD (1.5) than in the OND (-7.5). In model 7d the velocity gain b1 = 0.90 was also ~10% lower in the OFFD, and unlike in the OND for this cell, no relation was found between bias values and saccade amplitude. This result is compatible with the relatively small value of r1 estimated for model 8d in the OFFD. Although the average bias in 7d was equal to that in 2d (i.e., essentially 0), the fact that the bias was allowed to vary slightly improved the VAF as noted above.

    DISCUSSION
Abstract
Introduction
Methods
Results
Discussion
References

Reproducibility of analysis

The main objective of this paper was to provide an accurate description of IBN burst discharges to better understand the signal carried by these cells and how it relates to MN discharge during saccade generation. Our analysis, by using system identification techniques, has provided an objective method for comparing the goodness-of-fit of different downstream models, as well as an indication as to whether increasing the complexity of a model by adding parameters is warranted.

Before going on to discuss the significance and implication of our results, it is important to demonstrate that a given model fit can be extended beyond the data set for which it was constructed. This is verified in Fig. 16 for models 8d and 7d applied to SLIBN H0923. Figure 16A shows three examples of each of the best and worst fits applied to the original data set. The average VAF by the three best fits was 0.50 and 0.60 in models 8d and 7d, respectively, compared with 0.28 and 0.42, respectively, for the three worst fits. Figure 16B shows the fits to a new data set from the same models as used in Fig. 16A and with the same parameters (r0, r1, and b1 in model 8d and b1 in model 7d). In the latter model the biases were fit as parameters to the new set of 40 saccades and are independent of those derived in the original data set. The VAF by models 8d and 7d for this new population was 0.35 and 0.50, respectively, compared with 0.37 and 0.59 in the original data set. The average VAF by the three best fits was 0.51 and 0.59 in models 8d and 7d compared with 0.11 and 0.32, respectively, for the worst fits. These data show that the models predict IBN firing in the new data set almost as well as in the set used to calculate the parameters.


View larger version (43K):
[in this window]
[in a new window]
 
FIG. 16. Downstream models containing eye velocity term, amplitude-dependent term, and a single bias that was estimated for all saccades (model 8d) and an eye velocity term and a bias that was estimated separately for each saccade (model 7d) were used to predict firing behavior of LLIBN H0923 over 40 saccades. These models provided VAFs of 0.37 and 0.59 for models 8d and 7d, respectively. A: left-hand and right-hand sides illustrate 3 saccades (in range of 20-30°) for which the fit of these 2 models was the best and the worst, respectively. B: models estimated from original data set (40 saccades) were applied to new data set containing 40 different saccades to validate estimated models. Means and SD of saccade amplitudes included in the "original" and "new" data sets were(mean = 19.7 ± 10.7 vs. 17.1 ± 11.0, respectively). VAFs were 0.35 and 0.50 for models 8d and 7d, respectively, indicating that original models were robust. As in (A), left-hand and right-hand sides of figure in (B) illustrate 3 saccades (range of 20-30°) for which the fit of 2 of these models was the best and the worst, respectively. A and B: shaded area, actual firing frequency; superimposed solid line, estimated frequency; bottom 2 traces, accompanying eye velocity and position trajectories temporally shifted by estimated dynamic latency td.

Comparison of results of our metric analysis with other studies

Our analysis of the relationship between saccade metrics and IBN discharges revealed that the cells in our sample were generally similar to the primate IBNs described in the study of Scudder et al. (1988). Each neuron demonstrated a monotonic increase in the NOS in a burst (excluding preludes) with increasing amplitude of the horizontal component of a saccade [R = 0.79 this study vs. 0.85 Scudder et al. (1988)]. Consequently, in the terminology of Hepp and Henn (1983), our SLIBNs as well as LLIBNs would have been classified as "directional" bursters rather than "vector" bursters. In addition, we found that the preferred direction of our population of IBNs was closely aligned with horizontal (OND was -1.0 ± 30°). This result is in agreement with that of Scudder and colleagues (1988) who reported a similar spread in the ranges of preferred directions for their sample of LLIBNs and SLIBNs. In both studies, the range of preferred directions is much smaller than that previously reported in squirrel monkeys (Strassman et al. 1986b).

Several previous studies have reported a strong relationship between IBN peak firing rate and peak eye velocity during saccades. Our use of a spike density function rather than reciprocal intervals to represent unit firing rate resulted in a filtering of the firing frequency profile; we expected to find even better correlations. We were surprised to find a poor relationship between IBN peak firing frequency and saccade peak velocity (mean R = 0.39). This difference could stem from the relatively nonstereotyped skewed saccade profiles (e.g., left saccade in bottom panel of Fig. 11B; right saccade in bottom right panel of Fig. 16B) elicited during one of our protocols (the barrier paradigm; see METHODS) where for a given amplitude the peak velocity was quite variable.

Significance of dynamic lead time

Previous studies of BNs determined lead times based on a variety of methods such as the time of first spike or transition from low to high frequency (see a review in Hepp et al. 1989). One advantage of the present study was that an objective dynamic lead time was calculated for each unit, which accounted for the average influence of each spike on the saccade trajectory. This dynamic estimate of lead time was then implemented in our analysis of different models to correct for latency. A neuron's dynamic lead time was determined by shifting the unit discharge in time (td) until an optimal model fit was obtained for the simple but accurate dynamic model of Eq. 3: B(t) = r + b1Ė(t - td). For LLIBNs the dynamic estimate provided a significantly lower value for lead time than that obtained by measuring the onset of the first spike. For SLIBNs, lead times determined by the two methods were the same. As seen in Fig. 8, the lead times calculated by using each method were moderately well correlated. Interestingly, neurons arbitrarily classified as SLIBNs on the basis of the timing of the onset of the first spike still had significantly smaller (by 3.7 ms) dynamic lead times than neurons classified by using the same criterion as LLIBNs. We will consider this in a subsequent section. As expected given this small 3.7-ms difference, there was a considerable amount of overlap among the dynamic lead times estimated for SLIBNs and LLIBNs (see Fig. 8); 8 of 12 neurons classified as LLIBNs with the first-spike method were now in the short-lead range.

Because of the mechanisms leading to cell depolarization the initial few spikes of an excitatory burst impinging on MNs should have a different influence than those spikes occurring in the middle of the burst. This phenomenon is also dependent on initial firing rate. In past studies the objective duration of the LLIBN prelude has been difficult to evaluate. Our analyses reveal that on average the time difference between the first-spike and the dynamic lead time was ~9 ms. Put another way, for LLIBNs 9 ms of activity does not carry information related to the dynamics of saccades, but is rather a preparatory or anticipatory response. As concerns the influence of spikes well within the burst, it is theoretically possible that td can vary throughout the saccade. We have not tested for this.

Evaluation of models and their estimated parameters

Table 2 gives mean values and SDs of the parameters that were estimated in each of the nine downstream models tested for SLIBNs and LLIBNs and for the entire population taken together. Specific values for our example neurons are provided in the model fits illustrated in Figs. 11 and 12. We have shown that the most influential terms for predicting IBN firing rate are eye velocity and a bias. A model (2d) comprised of these two terms provided a fit that corresponded to a correlation coefficient (R) of 0.55 in a linear regression analysis (VAF = 0.30). The sign of the estimated eye velocity term was positive (ipsilateral) for every neuron in our analysis, compatible with previous studies demonstrating increases in IBN firing rate accompanying increases in eye velocity in a direction ipsilateral to the neuron's location. The estimated bias term was positive in all but the model of Eq. 4 (model 6d) where ICs were estimated as parameters; this result will be considered later.

For the population of neurons, the estimated bias and eye velocity gain terms did not change significantly between models 2d and 3d. This is undoubtedly due to the small acceleration term in model 3d. Indeed, the mean estimated value of b2 in model 3d was very low (-0.0004 ± 0.0009) and was unexpectedly in the OFF-, or contralateral, direction. We have shown earlier (Cullen et al. 1996) that the magnitude of the acceleration term is critically dependent on the choice of lead time. For example a smaller than optimal lead time predicts a larger acceleration term because more of the initial burst frequency increase must be accounted for by early changes in eye trajectory. Our objective estimate of best lead time is responsible for minimizing b2. The addition of a position term to model 2d only very slightly improved the fit of IBN discharges. The mean estimated eye position gain was low (-0.04 ± 1.3) and also was in the contralateral direction.

The existence of a bias term around which IBN discharges are modulated during saccades was included in previous models of the brain stem burst generator (Scudder 1988; Van Gisbergen et al. 1981). The average value of the bias term in our estimations using model 2d was 223 spikes/s; relatively close to the value in the Scudder 1988 model (160 spikes/s).

The bias value was significantly lower in the fourth-order nonlinear representation of the Van Gisbergen et al. (1981) model (model 4d) than in the linearized version (model 3d) and the eye velocity gain term b1 was correspondingly higher in 4d than in 3d. These changes are compatible with an improved fit by the nonlinear model; but the improved VAF by model 4d with respect to 3d of ~7% (Table 3) is less than we expected based on the nonlinear relationships between B(t) and Ė(t) shown by Van Gisbergen et al. (1981). In fact model 8d, a simple extension of model 2d, which included an amplitude-dependent term, better described our IBN discharges than the nonlinear function of eye velocity in 4d. A number of these points are illustrated in Fig. 17, which shows phase-plane plots [B(t) vs. Ė(t - td)] for cell H0925. In Fig. 17A the dynamic latency was chosen as the optimal one (17 ms) for this cell and the phase-plane plot trajectories of each saccade are close, on average, to straight lines. Figure 17B shows the phase-plane plots, based on optimal latency, of three saccades of different amplitudes, all of which were included in Fig. 17A. Model 2d (------) approximates well the average phase-plane characteristics, but one can see from Fig. 17B that variations among saccades that are dependent on saccade amplitude could be fit better by including a saccade amplitude-dependent term in the dynamic model (i.e., model 8d) or by allowing for the estimation of variable biases (i.e., model 7d). The nonlinear model (model 4d) for this cell also is plotted in Fig. 17A (- - -) and illustrates well the small nonlinear contributions that we have observed for all cells over the velocity range studied.


View larger version (29K):
[in this window]
[in a new window]
 
FIG. 17. Phase-plane plots illustrating dynamic relationship between B(t) and Ė(t - td) for cell H0925. A: phase plane plot obtained when optimal dynamic latency estimated for this cell (td = 17 ms) was used. Fits to models 2d (------) and 4d (- - -) are shown for same neuron. B: phase-plane plots for 3 saccades of different sizes that were included in A. Fit to model 2d (------) is superimposed. Note the different phase plane trajectories depend on saccade amplitude.

When a pole term was added to model 3d to achieve model 5d, the bias increased significantly (ICs were taken from the data). However, the improved fit with 5d was small compared with that offered by 6d when ICs were fit as parameters. Model 6d was the most complex model tested and included velocity, acceleration, and a pole term; the ICs were estimated as parameters. For this model the bias value varied greatly between neurons and the mean was an extremely low value (Table 2: -1,795 ± 2,058). The values of the other parameters in this model were also extremely variable. A large negative bias resulted from the algorithm seeking to offset the large IC transient that is itself a bias, given that the average time constant of the ICs is long relative to saccade duration and given by the value c = 2.4 s. Thus negative biases, in this context, do not have an actual physiological significance; they add to the relatively constant (during a saccade) positive IC transient to produce an overall bias close to that given by model 7d.

To pursue these arguments, note in model 6d that c approx  b2. The significance of this can be appreciated by noting from model 2d that &Bdot;(tproportional to  Ë(t - td), approximately, and substituting c&Bdot;(t) = b2Ë(t - td) in 6d yields model 7d. Taken together, these arguments indicate that in 6d the best-fit algorithm has selected 1) c to null the acceleration term and 2) the bias to null the effect of c. Thus model 6d essentially reduces to 7d (which has neither acceleration nor pole terms) and that provides an excellent description of IBN firing rate. This further substantiates the validity of our algorithms and analysis procedure because 7d, with no pole term, and 6d, with a pole term, were analyzed in quite different ways (Cullen et al. 1996). Model 7d has no pole term, is of comparable complexity to 6d, and is easier to estimate. This formulation provided more meaningful estimates of bias (always positive). However, models 6d and 7d contained too many parameters to be considered as useful models of IBN discharges. Therefore the correlation, which we frequently observed between estimated ICs and biases, led us to construct a more simple model that included an amplitude-dependent term (model 8d). This model actually described IBN spike trains slightly better than a third-order nonlinear function of eye velocity (model 4d). The sign of the amplitude-dependent term in this model was negative for 26 of 28 neurons; that is, this component of the firing rate was inversely correlated with saccade amplitude. Given the small number of parameters in model 8d (only 3) and its ability to fit the data, we suggest that it provides the most valuable description of IBN discharges during saccades.

The amplitude-dependent component of the IBN burst that was revealed by the r1 term in model 8d and by the amplitude-dependent biases estimated in model 7d was not expected. As illustrated in Fig. 17B, this relationship does not result from piecewise linear fits to the nonlinear relation between B and Ė (where the bias is the intercept) because the nonlinear terms are very small. Both Van Gisbergen et al. (1981) and Scudder et al. (1988) found that more spikes are needed for a given amplitude change in small (<2°) saccades than for an equivalent amplitude change in large saccades. This result is compatible with a nonlinear relation between B(t) and Ė(t), which states that the effective viscosity is higher for small saccades. However, these arguments do not apply to our data set because we have analyzed saccades larger than 5°; Fig. 17B shows that there are still differences for cell H0925 between phase plane plots for 14, 21, and 36° saccades. We argue in companion paper III (Cullen and Guitton 1997b) that the amplitude-dependent component of the IBN burst may be due to the influence of a signal from the superior colliculus on IBNs.

Acceleration term does not improve fits of IBN discharges

Our finding that the addition of an eye acceleration term did not improve the fit to IBN discharges was unexpected in light of the previous results of Van Gisbergen et al. (1981). In their analysis, these investigators correlated movement trajectory with burst frequency profiles by considering phase-plane trajectories. In their analysis of downstream models they first plotted the average values of eye velocity versus firing rate for saccades of a given size and direction and noted that the phase plane plot was an open loop. They found that by subtracting the appropriate eye acceleration sensitivity from the BN discharge, the phase plane plot collapsed to a line. Consequently, they concluded that BN activity was related to both eye velocity and acceleration during saccadic eye movements. We believe that the acceleration term in their analysis is the result of two factors: 1) a different frequency content in the neural activity and eye movement signals and 2) an inappropriate choice of lead time.

Our approach differed from that of Van Gisbergen et al. (1981). We used a spike density frequency plot in which the spike train was convolved with a Gaussian pulse whose 5-ms width was chosen so that the neural activity was effectively filtered at a comparable frequency to that used for filtering the eye movement signals (Cullen et al. 1996). (To fit a transfer function in a linear system, it is necessary to have the same frequency content in the input and output.)In the latter paper we demonstrated that if the reciprocal interval was used, the firing frequency profile contained frequencies higher than those of the saccade, and this could result in an overestimation of the acceleration term in the downstream models. Furthermore, the reciprocal-interval method is nonlinear and sensitive to noise particularly at high frequencies (Richmond et al. 1987; Sanderson and Kobler 1976). In Van Gisbergen et al. (1981) the BN firing was calculated by using reciprocal intervals, but they averaged the firing rate of BNs during several amplitude matched saccades to reduce the high-frequency noise. It is possible that the significant acceleration sensitivity reported in their analysis could still reflect a difference in the frequency composition of the input eye velocity and output unit activity signals. Moreover, as discussed in Cullen et al. (1996), the choice of lead time strongly influences the estimate of the eye acceleration term. Whereas in our analysis the lead time was determined through optimization, Van Gisbergen and colleagues (1981) used an algorithm in which the lead time was simply defined as the time between the onset of the burst (when, arbitrarily, the discharge rate exceeded 100 spikes/s) and the saccade. As noted in Cullen et al. (1996), we could effectively produce a significant eye acceleration coefficient by shifting, by a short time interval, the unit discharge relative to the corresponding eye movement. A final point worthy of mention is that Van Gisbergen and colleagues (1981) analyzed the difference between the OND and OFFD responses of BNs, whereas we have chosen to analyze IBN OND and OFFD responses separately. However, we do not believe that this difference would have greatly effected our characterization of IBNs, because approximately one-half of our neurons generated negligible OFFD responses, whereas others were characterized equally well by OND and OFFD models such that their overall response could be obtained by linearly adding model terms.

Comparison between SLIBNs and LLIBNs

We have considered the discharge of our IBNs in the context of the local feedback model of saccade control (Fig. 2B) whereby the burst generator is driven by a signal representing motor error [to be analyzed in companion paper III (Cullen and Guitton 1997b)] and in turn drives MNs that move the eye against the resistive forces of the eye plant. In doing so we have explicitly assumed that IBNs and their excitatory counterpart EBNs are similar in their discharge characteristics. This is a common assumption (Fuchs et al. 1985; Hepp and Henn 1983; Luschei and Fuchs 1972; Scudder 1988) and is based on observations that 1) both EBN and IBN discharges encode the amplitude and velocity of saccades (Kaneko et al. 1981; Scudder et al. 1988; Strassman et al. 1986a,b; Van Gisbergen et al. 1981; Yoshida et al. 1982) and 2) EBNs project to the IBN area (Sasaki and Shimazu 1981; Strassman et al. 1986a). An analysis of EBNs, similar to what we have done here, is essential to test this assumption.

An important result of our analysis is that there was no significant difference between SLIBNs and LLIBNs in the VAF by each of the models we tested; i.e., the dynamic signals encoded by SLIBNs and LLIBNs were similar in nature. Hence, although we found a significantly shorter lead time (by 3.7 ms) for SLIBNs, this is not forceful evidence in favor of these cells being closer in their functional connectivity to MNs as is frequently assumed (for review see Fuchs et al. 1985).

Scudder et al. (1988) suggested that 1) LLIBNs might project to omnipause neurons (OPNs) so as to generate the gradual decline in discharge of OPNs before saccades, whereas it is the SLIBNs that provide the pause in antagonist MNs, and 2) LLIBNs may not be inhibitory at all but rather long-lead EBNs (LLEBNs) that project to SLIBNs in their vicinity thereby generating the burst in these latter cells. (Our 3.7-ms lead time of LLIBNs over SLIBNs is compatible with this scheme.) Our results showing a similarity between SLIBNs and LLIBNs in terms of their dynamic predictions of eye trajectory (Tables 2 and 3) as well as a continuity between the SLIBN and LLIBN populations in terms of their dynamic lead time(s) argue against this categorization of IBNs into different functional populations. Furthermore, there is no experimental evidence showing IBN projections into their own area; Strassman et al. (1986b) have shown, using intracellular staining techniques, that at least some LLIBNs project directly to the motor nucleus. However, the LLIBN projection to contralateral abducens MNs may be functionally weak (Scudder 1988) because LLIBN dynamic lead time begins l5.5 ms before a saccade (i.e., ~6.5 ms before the pause in the antagonist MNs) and this is very long for a monosynaptic connection. As Scudder et al. (1988) also pointed out, there is no evidence of a decline in MN tonic discharge beginning during the LLIBN prelude. Clearly, more work is needed to determine the role of LLIBNs and LLEBNs in oculomotor control and to distinguish these cells from other long-lead BNs (LLBNs) that are implicated in neck motor control (Grantyn and Berthoz 1985).

Conclusions: implications of IBN head-fixed model fits

Previous analyses of MN discharges during smooth sinusoidal eye movements have demonstrated that they have a resting discharge and "steady state" eye velocity and acceleration sensitivities (Fuchs et al. 1988; Keller 1973; Van Gisbergen et al. 1981). In addition, these neurons demonstrate a postsaccadic slide in discharge (required by the pole term in Eq. 1) (Fuchs et al. 1988; Goldstein and Robinson 1984; Stahl and Simpson 1995). However, the burst discharge of MNs during saccades was not carefully examined using techniques comparable to those that we used here. A thorough characterization of the dynamics of MN activity during saccades would provide an appropriate dynamic model of their discharge, which is important to evaluate the relative contributions of the many premotor inputs to MNs during saccades.


View larger version (27K):
[in this window]
[in a new window]
 
FIG. 18. A: schematic of inputs to abducens nuclei on left and right (ABDR and ABDL, respectively) that are active during rightward saccades. 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), and RPL (position cells). First subscript denotes OND of each cell type's response; final subscript denotes the side of the brain stem where the cell is located. Inputs from the ipsilateral (right, in this example) VN-PH include projections from I-PVPR, LBPR, LPR, RBPR, and RPR, as well as LE/HR (eye-head velocity cells). Excitatory and inhibitory projections are indicated by gray and solid circles, respectively. Equations describing OND and OFFD responses of the BN population B(t) are based on results of our dynamic analysis. For additional details see text. B: simple schematic illustrating signal processing required to remove bias and amplitude-dependent signals carried by BNs to ocular motoneurons (MNs) during saccadic eye movements. Removal of this component of BN discharges is required at level of the motor nucleus if transfer function between MNs and plant dynamics is linear. In addition, it is necessary to offset influence of this component of burst at the neural integrator.

To calculate the transfer function between a single BN (BNi) and the eye movement output, we have pooled the transfer functions together. What is the significance of pooling model parameters and is it a legitimate evaluation of the average BN signal? Van Gisbergen and Van Opstal (1989) have considered this problem with regard to the relationship between ocular MNs and plant time constants. To discuss this point they propose the thought experiment wherein an eye muscle motor unit is paralyzed while its MN is being recorded. The eye movement will surely not be altered by this very minor intervention and thus, they argue, the considerable differences in the time constant values (b1/b0 in Eq. 2) among individual MNs reflect more differences in inputs to, and intrinsic properties of, MNs rather than the existence of different time constants in the muscle motor units. It follows then that the average b1/b0 does not necessarily give the plant time constant.

The same thought experiment can be considered for BNs. Suppose we were to block the axon being studied. The eye movement would undoubtedly remain virtually unchanged because of the large number of BNs acting together. Hence, whereas the transfer function between the recorded cell and the output would remain the same, the neuron is no longer affecting the output. This is because the cell's transfer function is accounting more for the simultaneous input of many BNs onto a MN, and in turn many MNs driving a movement, rather than for its own effect on the movement. It follows that any transfer function is a sample within the IBN population variance itself. However, contrary to MNs, we have found only a small variability in the parameter and VAF values for each IBN and the distribution of parameter values is close to normal. It therefore appears legitimate in our case to take mean values of each parameter and give a mean transfer function for our population.

Figure 18A reviews inputs to the ABD on the right and left (ABDR and ABDL, respectively) that are implicated in the generation of rightward saccades. Although many cell types have been described in the VN and prepositus-hypoglossi nucleus (PH), only a few types have been identified as projecting to ABD. The major excitatory input to ABDR from VNL arises from type I position-vestibular-pause cells (I-PVPL) (cat: Escudero et al. 1992; monkey: Cullen 1991; McCrea et al. 1987; Scudder and Fuchs 1992). During saccades to the right (their OND) most of these cells have tonic activity, some burst, and only a minority pause (rhesus monkey: Scudder and Fuchs 1992). A minor projection arises contralaterally from burst position (RBPL) and position (RPL) cells that are also active during right saccades as indicated by the first subscript. Inhibitory projections from VN to ABDR arise in VNR. However, I-PVPR, LBPR, and LPR have their OND to the left and they pause for saccades to the right (- - -). Another cell type in the VNR, the eye-head velocity cells (LEHR), appear to be inhibitory to ABDR (Scudder and Fuchs 1992) but these also pause for rightward saccades.

With regards to the PH, the OND of its neurons is to the ipsilateral side. Excitatory inputs to the ABDR come from cells whose discharge is proportional to eye position to the right (RPR) (Escudero et al. 1992). RBPR in the right VN-PH appear to also send excitatory projections to the ABDR (McCrea et al. 1987). Inhibitory input comes from LPL in PHL, but these neurons pause for rightward saccades.

Let us now consider the projection of EBNs and IBNs to ABD and to the VN-PH, the latter being of considerable interest because of the key role VN-PH plays in the generation of the tonic signal needed to hold the eye at a fixed position in the orbit. As we saw above, for a rightward saccade PHR increases its tonic activity, which is projected to MNs, and PHL decreases its activity. The projection into the PHL of right IBNs (IBNR) is thought responsible for this decrease in tonic activity. Similarly, the projection into the PHR of right EBNs (EBNR) increases the tonic activity. In the simple model of MN discharge given by Eq. 2, it is assumed that the tonic signal is obtained by integrating the BN eye velocity signal (reviewed in Fuchs et al. 1985; Robinson 1975). VN-PH classically has been thought of as this neural integrator (Fig. 2B) (see 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). Our conclusion that IBN firing frequency is not simply proportional to Ė complicates the interpretation of how the integration process is performed.

The net signal of BNs onto ABD and VN-PH is the combination of the OND and OFFD discharges of EBNs and IBNs. For a saccade to the right, the right IBN population produces the signal
model 8d: <IT>B</IT><SUP>+</SUP>(<IT>t</IT>) = 282 − 4.1Δ<IT>E</IT>+ 0.84<IT><A><AC>E</AC><AC>˙</AC></A></IT> (5)
The left IBN population produces an OFFD for right saccades. It is difficult to give an average quantitative equation for this discharge because 50% of cells have no OFFD, whereas the remaining neurons produce, with model 8d, a VAF that varies from 0 in the worst cell to 0.33 in the best cell, H0925. If we assume that the velocity coefficient b1 varies linearly across the population of cells that do have an OFFD from 0 in the worst cell to 0.86 (Table 4) in the best cell (H0925), then b1 = -0.86/4 = -0.22 (b1 has a negative sign because we are now considering a left IBN discharging for a rightward (+) saccade). Because the values of r0 and r1 were negligible for the best cell (H0925), the average OFFD of left IBNs simplifies to
<IT>B</IT><SUP>−</SUP>(<IT>t</IT>) = −0.22<IT><A><AC>E</AC><AC>˙</AC></A></IT> (6)
If we assume that the average EBN and IBN discharges are similar, then as shown in Fig. 18A the average net inhibition on the left VN-PH and ABDL will be the sum between the inhibition from right IBNs and excitation from left EBNs
<IT>B</IT>(<IT>t</IT>) = 282 − 4.1Δ<IT>E</IT>+ 0.62<IT><A><AC>E</AC><AC>˙</AC></A></IT> (7)
and similarly, the average net excitation on the right VN-PH and ABDR will be
<IT>B</IT>(<IT>t</IT>) = 282 − 4.1Δ<IT>E</IT>+ 0.62<IT><A><AC>E</AC><AC>˙</AC></A></IT> (8)
The question now arises as to how this signal is utilized. This is summarized in Fig. 18B. If MNs require only an Ė term, then the bias and amplitude-dependent term must be offset at ABDR. However, as reviewed in the foregoing equations, our current knowledge of the inputs to ABDR shows no signal that could offset 282 - 4.1Delta E (because the inhibitory signals pause for right saccades). However, there is evidence that the plant has nonlinear properties with phase plane plots having characteristics similar to Fig. 17B (Van Gisbergen et al. 1981). This is compatible with the amplitude-dependent term diminishing with increasing saccade amplitude. If this is so, we predict that our optimization techniques applied to MNs will reveal in the MN discharge bias and amplitude-dependent terms. With regards to the VN-PH complex, if this circuit acts as an integrator then the bias and amplitude term should be offset (Fig. 18B). Alternatively, as proposed by Galiana (see Galiana and Guitton 1992) the VN-PH may contain a circuit that simulates the plant (a plant model) and, as was noted above for MNs, there may be no need to explicitly offset these terms at this level.

In conclusion, our analysis using rigorous system identification algorithms reveals that understanding the signal processing in even a "simple" structure like the ABD of the head-fixed monkey is a very complex and labor-intensive process. We cannot help but be struck by the awesome task of resolving the much more complex circuits involved in mediating higher brain functions.

    ACKNOWLEDGEMENTS

  We are grateful to C. G. Rey for useful discussions and contributions to the development of the analysis procedure and software and to H. L. Galiana for helpful discussions. We thank Dr. J.A.M. Van Gisbergen for a number of useful suggestions that helped to 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
Abstract
Introduction
Methods
Results
Discussion
References

0022-3077/97 $5.00 Copyright ©1997 The American Physiological Society