JN Watch the video to see how APS reaches out to developing nations.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 92: 905-925, 2004. First published March 31, 2004; doi:10.1152/jn.01234.2003
0022-3077/04 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
92/2/905    most recent
01234.2003v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (25)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Green, A. M.
Right arrow Articles by Angelaki, D. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Green, A. M.
Right arrow Articles by Angelaki, D. E.

An Integrative Neural Network for Detecting Inertial Motion and Head Orientation

Andrea M. Green and Dora E. Angelaki

Department of Anatomy and Neurobiology, Washington University School of Medicine, St. Louis, Missouri 63110

Submitted 19 December 2003; accepted in final form 23 March 2004


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DEVELOPMENT
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The ability to navigate in the world and execute appropriate behavioral responses depends critically on the contribution of the vestibular system to the detection of motion and spatial orientation. A complicating factor is that otolith afferents equivalently encode inertial and gravitational accelerations. Recent studies have demonstrated that the brain can resolve this sensory ambiguity by combining signals from both the otoliths and semicircular canal sensors, although it remains unknown how the brain integrates these sensory contributions to perform the nonlinear vector computations required to accurately detect head movement in space. Here, we illustrate how a physiologically relevant, nonlinear integrative neural network could be used to perform the required computations for inertial motion detection along the interaural head axis. The proposed model not only can simulate recent behavioral observations, including a translational vestibuloocular reflex driven by the semicircular canals, but also accounts for several previously unexplained characteristics of central neural responses such as complex otolith–canal convergence patterns and the prevalence of dynamically processed otolith signals. A key model prediction, implied by the required computations for tilt–translation discrimination, is a coordinate transformation of canal signals from a head-fixed to a spatial reference frame. As a result, cell responses may reflect canal signal contributions that cannot be easily detected or distinguished from otolith signals. New experimental protocols are proposed to characterize these cells and identify their contributions to spatial motion estimation. The proposed theoretical framework makes an essential first link between the computations for inertial acceleration detection derived from the physical laws of motion and the neural response properties predicted in a physiologically realistic network implementation.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DEVELOPMENT
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The vestibular system plays a critical role in spatial perception and motor compensation for head movement. As in any inertial guidance system, vestibular information about the angular and linear components of head motion is provided by distinct types of sensors. The semicircular canals act as integrating angular accelerometers to provide the brain with an estimate of rotational head velocity (Fernandez and Goldberg 1971Go), whereas the otolith organs behave as linear accelerometers that sense net gravito-inertial force (Angelaki and Dickman 2000Go; Fernandez and Goldberg 1976aGo,1976bGo). Because everyday activities such as locomotion give rise to both translations and rotations within a gravitational field, integration of both sets of vestibular signals is vital for motion detection.

Among the most important and well-studied functions of the vestibular system is its contribution to gaze stabilization. Vestibular stimulation elicits short-latency compensatory ocular responses to head motion known as the vestibuloocular reflexes (VORs) that ensure the ability to maintain ocular stability and thus high visual acuity while moving. Early investigations of the VOR focused mainly on the sensorimotor transformations associated with rotational motion [rotational VOR (RVOR)]. More recently, compensatory responses to translation [translational VOR (TVOR)] have been investigated (Angelaki 1998Go; Busettini et al. 1994Go; Paige and Tomko 1991aGo,1991bGo; Paige et al. 1998Go; Schwarz and Miles 1991Go; Schwarz et al. 1989Go; Telford et al. 1997Go).

Linear acceleration information is provided to the brain by primary otolith afferents. However, linear accelerometers (including the otoliths) respond similarly to inertial and gravitational accelerations (Einstein's equivalence principle; Einstein 1908Go). Thus, otolith afferents provide inherently ambiguous sensory information, given that the encoded acceleration could have been generated during either actual translation or a head reorientation relative to gravity (Angelaki and Dickman 2000Go; Fernandez and Goldberg 1976aGo,1976bGo). Yet, behavioral responses to tilts and translation are different. In the oculomotor system, for example, lateral translation elicits horizontal eye movements (Angelaki 1998Go; Paige and Tomko 1991aGo; Schwarz et al. 1991Go; Telford et al. 1997Go), whereas roll tilt generates mainly ocular torsion (Angelaki and Hess, 1996Go; Crawford and Vilis 1991Go; Haslwanter et al. 1992Go; Seidman et al. 1995Go). The integration of available sensory information to ensure that otolith signals are correctly processed to generate appropriate perceptual or motor responses thus represents an essential task for the central nervous system.

It has long been proposed that the brain integrates information from both otolith and semicircular canal afferents to differentiate translation from tilt (Guedry 1974Go; Mayne 1974Go; Young 1984Go). Theoretically, the canals should then also contribute to driving the TVOR (Glasauer and Merfeld 1997Go; Merfeld 1995Go; Merfeld and Zupan 2002Go; Mergner and Glasauer 1999Go; Zupan et al. 2002Go). This has been confirmed experimentally by examining oculomotor responses to simultaneous roll tilt and translation stimuli, carefully matched to ensure that the gravitational and translational components of acceleration along the interaural head axis cancelled one another out. Despite the absence of a net lateral acceleration stimulus to the otoliths, horizontal ocular responses appropriately directed to compensate for the translational component of motion were nevertheless elicited (Angelaki et al. 1999Go; Green and Angelaki 2003Go). The contribution of semicircular canal cues to the generation of these horizontal eye movements was directly demonstrated by the fact that they were no longer evoked in canal-plugged animals (Angelaki et al. 1999Go). Recently, it has been shown that these canal-driven responses represent an extra-otolithic TVOR that exhibits dynamic properties and a dependency on viewing distance similar to those of the purely otolith-driven reflex (Green and Angelaki 2003Go). Quantitative analyses demonstrated that the horizontal eye velocity profile associated with this extra-otolith driven TVOR was best correlated with angular head position, suggesting that angular velocity signals from the semicircular canals are processed by an additional neural integrator in the TVOR pathways. It was proposed that the integrative network known as the "velocity storage integrator" might perform this function (Green and Angelaki 2003Go).

These experimental results are consistent with the predictions of several theoretical studies that propose that the brain explicitly constructs internal estimates of gravity and translational acceleration (Glasauer et al. 1997Go; Merfeld 1995Go; Merfeld and Zupan 2002Go; Merfeld et al. 1993bGo; Mergner and Glasauer 1999Go; Zupan et al. 2002Go). To do so, the brain must effectively solve a vector differential equation that relies on an estimate of head velocity to calculate the rate of change of gravitational acceleration in a head-fixed reference frame. The solution of any differential equation requires the process of temporal integration. Thus, calculation of the instantaneous gravity vector implies a central neural integration of head angular velocity information, in agreement with experimental observations (Green and Angelaki 2003Go). This gravity estimate can then be combined with the net gravito-inertial acceleration sensed by the otoliths to extract the translational component of acceleration. Although such models have provided computationally rigorous solutions to the problem that are consistent with many experimentally observed behaviors, they are difficult to relate directly to the response properties of individual neurons. Specifically, these models use 3-component vector representations to perform the calculations required to compute head orientation in 3-dimensions (e.g., vector cross-products), whereas the instantaneous firing rate of an individual neuron is a scalar quantity. Thus, although significant progress has been made in outlining the computational requirements for resolving the tilt–translation ambiguity problem, few predictions have been made regarding the types of neural responses expected from a network that effectively implements these calculations. Specifically, how could cells involved in these nonlinear vector computations be identified and what types of experimental and analytical approaches should be used to interpret their responses?

The goals of the current investigation were twofold: 1) to illustrate how an integrative network within the traditional VOR circuitry (Green and Angelaki 2003Go) could implement these abstract vector computations to distinguish head translations from reorientations relative to gravity; 2) to investigate the predictions of such a structure at the neural level, with the goal of elucidating experimental paradigms appropriate for identifying and characterizing the physiological response properties of neurons involved in inertial motion detection. Preliminary versions of these results have been presented in abstract form (Green and Angelaki 2002Go; Green et al. 2002Go).


    MODEL DEVELOPMENT
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DEVELOPMENT
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
General framework

The general framework for the current theoretical investigation was previously proposed (Green and Angelaki 2003Go) and is summarized in the schematic of Fig. 1, which illustrates a feedforward model for the VOR. The classical model for sensorimotor transformations in the RVOR (Robinson 1981Go; Skavenski and Robinson 1973Go) consists of a parallel set of pathways (labeled as RVOR in Fig. 1) that convey angular head velocity signals (), sensed by the semicircular canals, C(s), to extraocular motoneurons (Mn) both "directly" (i.e., via the short latency 3–neuron–arc VOR pathways) and "indirectly" via the oculomotor neural integrator NI2 (e.g., Cannon and Robinson 1987Go). The bottom projection (labeled as TVOR in Fig. 1) represents a proposal for the dynamic processing of linear acceleration signals () sensed by the otoliths, O(s), during translation (Angelaki et al. 2001aGo; Green and Galiana 1998Go; Musallam and Tomlinson 1999Go). A second neural integrator, NI1, accounts for the recently established contribution of integrated semicircular canal signals to the TVOR (Angelaki et al. 1999Go; Green and Angelaki 2003Go). This integrative network, which could be the so-called velocity storage integrator (Raphan et al. 1977Go, 1979Go), computes a dynamic estimate of head orientation relative to gravity that, when combined with otolith sensory signals, could be used to extract the component of linear acceleration due to translation (i.e., on cell VOT) (Green and Angelaki 2003Go).



View larger version (15K):
[in this window]
[in a new window]
 
FIG. 1. Overview of sensorimotor processing in the vestibuloocular reflexes. Schematic illustrates a classical feedforward diagram for the rotational vestibuloocular reflex (RVOR) (Skavenski and Robinson 1973Go) extended to incorporate the contributions of both otolith and semicircular canal signals to the translational VOR (TVOR) (Green and Angelaki 2003Go). Inputs to the model are angular head velocity , sensed by the semicircular canals C(s), and linear acceleration , sensed by the otolith organs O(s). Output of the model is eye position . Boxes in the schematic are dynamic elements that represent either a sensor [C(s) and O(s)], the motor plant [P(s)], or a neural filtering process [NI1(s), NI2(s)]. Two central filtering processes include the "oculomotor integrator" [NI2(s)] and an integrative network for canal signals [NI1(s)] that could be the "velocity storage integrator" (Raphan et al. 1979Go). Circles are summing junctions used to represent particular cell populations including vestibular-only (VO) cells, VOR and VOT, that mediate sensory signal flow in the RVOR and TVOR pathways, respectively, premotor eye-movement–sensitive (EM) neurons and motoneurons (Mn). Circle labeled with an X in the schematic indicates a proposed head-orientation–dependent modulation in the coupling between semicircular canal sensory signals and the integrative network NI1.

 
In general, the semicircular canal signals that must be combined with otolith information to extract translational acceleration depend on current head orientation. For example, both roll rotation from upright and yaw rotation from supine positions cause a reorientation relative to gravity that stimulates the otoliths along the head's interaural axis. The extra-otolith signals required to ensure that reorientations relative to gravity are not incorrectly interpreted as translation must therefore arise mainly from vertical canals in the upright orientation, yet from the horizontal canals in the supine orientation. Thus, the circle labeled with an X at the input of NI1 illustrates a required head-orientation–dependent modulation in the coupling between semicircular canal sensory signals [i.e., at the output of C(s)] and this integrative network. The result of such a nonlinear (multiplicative) coupling of semicircular canal signals to integrator NI1 is the construction of an internal dynamic estimate of head reorientation relative to gravity (or equivalently a dynamic estimate of the gravity vector in a head-fixed reference frame), as described below.

Mathematics of tilt–translation discrimination

Several theoretical studies have proposed models for tilt–translation discrimination that combine otolith and canal sensory information in 3 dimensions (3D) to extract the component of linear acceleration that is due to translation from that due to head reorientation relative to gravity (Angelaki et al. 1999Go; Glasauer and Merfeld 1997Go; Merfeld 1995Go; Merfeld and Zupan 2002Go; Mergner and Glasauer 1999Go; Zupan et al. 2002Go). All are based on the premise that canal information about angular velocity is used to keep track of changes in orientation of the gravity vector relative to the head (in which the vestibular sensors are fixed), as described by the first-order differential equation (Goldstein 1980Go)

(1)
where and are vector representions of gravity and angular velocity, respectively, in head-fixed coordinates and x denotes a vector cross-product. Using the additional information that the net acceleration sensed by the otoliths is the vectorial difference of translational () and gravitational () components

(2)
the translational acceleration component can easily be computed by solving the vector differential equation resulting from the combination of Eq. 1 and 2 (Angelaki et al. 1999Go; Hess and Angelaki 1997Go; Viéville and Faugeras 1990Go)

(3)
Equivalently, the same basic equations have been used to postulate that the otolith measurement of gravito-inertial force is explicitly resolved into central estimates of gravitational and translational accelerations (e.g., Glasauer and Merfeld 1997Go; Merfeld and Zupan 2002Go). Specifically, an internal 3D estimate of the gravity vector in head coordinates is obtained from Eq. 1 (i.e., = –{int} x, assuming a known set of initial conditions). The translational acceleration component can be subsequently obtained from Eq. 2. These implementations are mathematically equivalent (Angelaki et al. 2001bGo). Both involve a central (nonlinear) neural integration of angular velocity, a theoretical prediction that is consistent with recent experimental observations (Green and Angelaki 2003Go).

Although vector Eqs 1 and 2 can be used to discriminate tilt from translation in 3D, a key focus of the current study is to predict the responses of central neurons whose firing rates are scalar quantities. Thus, to simplify the problem we have chosen to examine tilt–translation discrimination along only the interaural head axis. In particular, we can expand the vector cross-product in Eq. 1 into components as

(4)
where , , and are unit vectors in a right-handed coordinate system along the x [nasooccipital (NO)], y [interaural (IA)], and z [dorsoventral (DV)] axes, respectively. Integration of each vector component in Eq. 4 yields the gravitational acceleration along the x, y, and z axes, as illustrated in Fig. 2. Because we restrict consideration to tilt–translation discrimination along the interaural axis (i.e., y-axis associated with unit vector ) we focus on the calculation of gy (shaded region in Fig. 2)

(5)
According to Eq. 5, calculation of the gravity component, along the interaural axis (y-axis) requires estimates of the components of gravitational acceleration along the other 2 axes, illustrating the 3D nature of the problem (i.e., coupled integrative networks in Fig. 2). However, the solution can be simplified if we decouple computation of gy from its dependency on gx and gz by limiting consideration to conditions where x- and z-axis translations are small (i.e., fx = fz {approx} 0). Under these conditions, where gx {approx}{alpha}x and gz {approx}{alpha}z, gy can be approximated as

(6)
The preceding equation indicates that, for the restricted conditions we consider here, the component of gravity along the interaural head axis can be computed by integrating estimates of yaw and roll head velocities that have been premultiplied by the net instantaneous accelerations sensed by otolith afferents along the nasooccipital and dorsoventral axes, respectively.



View larger version (33K):
[in this window]
[in a new window]
 
FIG. 2. Schematic illustration of the computations necessary to calculate the gravitational component of acceleration in a head-fixed reference frame as the head reorients relative to gravity (i.e., = –{int}x; also see Eq. 4). Computations involve multiplications between the vector components (i.e., x-, y-, and z-axis components) of angular velocity signals (solid lines) and gravitational accelerations (dashed lines), followed by integration. Notice that calculation of the gravitational acceleration along any particular axis depends on the components of gravitational acceleration along the 2 other axes, illustrating the requirement for coupled integrative networks in 3D. Here we consider tilt–translation discrimination along only the interaural axis (i.e., y-axis) and thus focus on the computation of gy (shaded region). To simplify the problem we decouple computation of gy from its dependency on gx and gz by assuming that x- and z-axis translational accelerations are small (i.e., fx = fz {approx} 0) such that gx {approx}{alpha}x and gz {approx}{alpha}z (see MODEL DEVELOPMENT).

 
Proposed network for tilt–translation discrimination

General description.    Figure 3A illustrates one of many possible integrative networks (representing NI1 expanded from Fig. 1) that could perform the computations of Eq. 2 and 6. Circles in the schematic represent summing junctions that are used to represent different vestibular-only (VO; i.e., eye-movement–insensitive) cell populations, whereas boxes represent dynamic elements or filters. These include first-order dynamic approximations of the semicircular canals, C(s) = Tcs/(Tcs + 1), (Fernandez and Goldberg 1971Go) and the otolith organs, O(s) = 1/(Tos + 1) (Fernandez and Goldberg 1976bGo) as well as the neural filter, CLP(s), which represents a low-pass internal model of the semicircular canals [CLP(s) = 1/(Tcs + 1)].



View larger version (29K):
[in this window]
[in a new window]
 
FIG. 3. Proposed model for detecting head orientation and translation. A: implementation of the integrative network NI1 (illustrated in Fig. 1). Circles labeled VO1 through VO5 represent different populations of VO cells. Box labeled CLP(s) represents a low-pass internal model of the semicircular canals [CLP(s) = 1/(Tcs + 1)]. Inputs to the model include yaw and roll velocities, {omega}z and {omega}x, sensed by the horizontal and vertical canals [C(s) = Tcs/(Tcs + 1)], respectively, and interaural acceleration {alpha}y, sensed by the otolith organs [O(s) = 1/(Tos + 1)]. X's indicate a multiplicative modulation in the strengths of the projections from cells VO1 and VO2 onto cell VO4 by nasooccipital and dorsoventral accelerations, {alpha}x and {alpha}z, respectively. Note that these multiplications implement the head-orientation–dependent modulation illustrated in Fig. 1 at the input to NI1. Here, this nonlinear signal interaction is incorporated as part of the network. Parameters associated with different pathways represent the strength or weight of the projection. B: previously proposed implementation of the premotor network for the RVOR and TVOR consisting of a feedback loop around an internal model [F(s) = Kf/(Tfs + 1)] of the eye plant [P(s) = Kp/(Tps + 1)] (Angelaki et al. 2001aGo; Green and Galiana 1998Go). Semicircular canal signals project directly onto eye-contra–sensitive premotor vestibular neurons (EMC) that represent the key interneurons in the shortest-latency RVOR pathways. Eye-ipsi–sensitive premotor neurons (EMI) are proposed to mediate the shortest-latency TVOR pathways (Angelaki et al. 2001aGo) and receive projections from cell VO3 in A. Separate but structurally identical premotor networks were assumed to simulate horizontal and torsional eye movements. Model parameters are: pdh = 1, ph = 2.568, pd{nu} = 0.56, p{nu} = 2.568, q1 = 0.25, q2 = 0.22, q3 = 0.1975, K1o = 0, K2o = 0, K1a = 1, K2a = 1, K1b = 0.1, K2b = 0.1, K3 = 0.061, K3b = 1.6, K3c = 6.25, a = 0.19, b = 0.75, d1 = 0.21, d2 = 1.1, Tc = 6, To = 0.0159, Kf = 2.81, Kp = 1, Tf = 0.25, Tp = 0.25. Parameters qoh = 0.39, eh = 0.03, and qot = 0.04, et = 0.15 are qo and e values for the horizontal and torsional premotor systems, respectively. Inputs {omega}z and {omega}x are in units of deg/s, {alpha}y is in units of cm/s2, and {alpha}z and {alpha}x are in units of g.

 
For the simplified case of tilt–translation discrimination along the interaural axis considered here (Eq. 6), inputs to the model include interaural acceleration, {alpha}y, sensed by the otoliths (mainly the utricles), and yaw and roll head velocities, {omega}z and {omega}x, sensed mainly by the horizontal and vertical semicircular canals, respectively. The 2 orthogonal accelerations, {alpha}x and {alpha}z, are proposed to modulate the strengths of semicircular canal projections onto the network. Specifically, by multiplying the yaw and roll head velocity projections onto VO4 by either {alpha}x or {alpha}z, the network effectively implements Eq. 6. Accordingly, cell VO5 encodes gy. The network output arises from cell VO3, which performs the addition implied by Eq. 2 to extract translational acceleration (i.e., fy = {alpha}y + gy). VO3 then projects directly into the downstream premotor TVOR pathways.

We focus here on an integrative network of VO neurons for two reasons: 1) populations of VO neurons have recently been observed that code mainly for either translation (vestibular and fastigial nuclei; Angelaki et al. 2003Go) or tilt (vestibular nuclei; Zhou et al. 2000Go); 2) the ability to distinguish head tilts and translations is important for both perceptual and motor responses. Both observations suggest that a network for distinguishing tilt versus translation occurs upstream of the premotor oculomotor networks of eye-movement–sensitive neurons (i.e., network in Fig 3B). In keeping with a previous proposal that the required integrative network could represent that known as the velocity storage integrator (Green and Angelaki 2003Go), we assume a model structure that is based on a feedback implementation of this integrative network originally proposed by Robinson (1977)Go. Accordingly, neurons that receive direct vestibular sensory projections (i.e., cells VO1, VO2, and VO3) are each interconnected in a feedback loop with an assumed common internal low-pass canal model, CLP(s), to form a distributed integrative neural network. Potentially, many other model structures could be used to implement the requirements for tilt–translation discrimination described by Eq. 2 and 6. All such networks, however, will be common in the requirements for 1) performing a central neural integration of canal signals and 2) implementing a head-orientation–dependent coupling between canal and otolith-derived sensory information. The key arguments to be made in this study focus on the implications of these common requirements and therefore are largely independent of the particular model structure. In the example network proposed here, the required integration is a distributed network property implemented by positive feedback loops through the low-pass filter, CLP(s), whereas the head-orientation–dependent sensory coupling is implemented by the multiplicative interactions, denoted by X's, in Fig. 3A.

The premotor oculomotor network used to simulate compensatory VOR responses is illustrated in Fig. 3B. It represents a feedback implementation of the "neural integration and eye plant compensation network" of Fig. 1 and was previously described in detail (Angelaki et al. 2001aGo; Green 2000Go; Green and Galiana 1998Go). For simplicity, we assume separate, but structurally identical, premotor networks for driving horizontal and torsional eye movements.

Dynamic processing of sensory signals.    The goal in this section is to illustrate the relationship between the dynamic computations performed by the model in Fig. 3A and the necessary computations for tilt–translation discrimination presented above. In particular, we will demonstrate that the network can perform the computations described by Eq. 6 to construct an internal estimate of dynamic head orientation relative to gravity on cell VO5. For descriptive purposes we will express this cell's response in the Laplace domain. Note, however, that because canal-related projection weights onto cell VO4 (i.e., projections from cells VO1 and VO2) modulate as a function of {alpha}x(t) or {alpha}z(t), the system is nonlinear. Laplace domain descriptions thus cannot generally be used here to predict response trajectories over time. However, they can approximate the dynamic characteristics of the system for small movements about an average static head orientation (i.e., a given operating point) when the system exhibits close to linear performance. In particular, for small head movements around a given pitch orientation we can assume static approximations to the linear accelerations along the NO and DV axes (i.e., {alpha}x {approx} –sin {phi} and {alpha}z {approx} cos {phi} in units of g) where angle {phi} describes the pitch angle from upright. At mid-high frequencies (>0.1 Hz), where semicircular canal cues were previously confirmed to play a crucial role in resolving ambiguous otolith sensory information (Angelaki et al. 1999Go), the response of cell VO5 can then be approximated as

(7)
where GVOH5 and GVOV5 are static gain terms (GVOH5 = phK1aK3b/Tc and GVOV5 = p{nu}K2aK3b/Tc). Equation 7 represents a high-frequency, small-angle approximation to the general expression for cell VO5 that assumes model parameters chosen to ensure close to ideal tilt–translation discrimination (see APPENDIX).

Because 1/s is the Laplace domain description of an integrator, Eq. 7 illustrates that the network integrates angular velocity signals, {omega}z and {omega}x, that have been multiplied by {alpha}x or {alpha}z, as required to construct an internal estimate of the gravitational acceleration component along the interaural axis (i.e., compare Eq. 7 with Eq. 6). Given an internal (scaled) estimate of gy on cell VO5, Eq. 2 predicts that the translational component of the acceleration, fy, can be obtained by combining this estimate with the net interaural acceleration signal, {alpha}y. In the model this occurs on cell VO3 [i.e., VO3(s) = q1{alpha}y(s) + K3cVO5(s) {approx} q1[{alpha}y(s) + gy(s)] {approx} q1fy(s)]. Note that, although the analytical expressions presented in this section are valid only for small dynamic head reorientations relative to gravity, the model can simulate appropriate responses even for large tilts in all pitch head orientations, assuming that any concurrent translational acceleration is mainly directed along the y-axis of the head (i.e., fx, fz {approx} 0; see above). Further details of the Laplace domain descriptions of cell and motor responses are provided in the APPENDIX.

Model simulations

The proposed model was implemented using the MATLAB simulation toolbox SIMULINK (MathWorks, Natick, MA). All simulations were performed using a fixed-step Runge–Kutta numerical integration routine (ode4 in SIMULINK) with time steps fixed at 0.01 s. The model parameters provided in the caption of Fig. 3 were chosen to satisfy the criteria outlined in the APPENDIX.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DEVELOPMENT
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The basic predictions of the model in Fig. 3 will be explored during pure rotations and translations in different head orientations as well as for combinations of rotational and translation motion stimuli that have recently been used to investigate tilt and translation responses. To begin, we will illustrate the ability of the model to generate appropriate behavioral responses. We will then use the model to predict the neural response properties expected from the cell populations involved in estimating spatial orientation and inertial motion.

Behavioral responses

Frequency response predictions.    Figure 4 illustrates predicted horizontal and torsional ocular responses during yaw rotation, interaural translation, and small angle roll rotation (e.g., <30°) from upright orientation. Both yaw rotation and interaural translation are predicted to elicit large horizontal eye movements (Fig. 4, A and C; solid black curves), whereas head roll generates torsional ocular responses (Fig. 4B, solid gray curves), as required for gaze stabilization. Small torsional responses are also elicited in response to head translation, as observed experimentally (Fig. 4C, gray traces; Angelaki 1998Go; Paige and Tomko 1991aGo).



View larger version (28K):
[in this window]
[in a new window]
 
FIG. 4. Behavioral predictions of the model. Horizontal (black) and torsional (gray) ocular gains and phases are plotted as a function of frequency during (A) yaw rotation, (B) roll rotation, and (C) interaural translation in the upright orientation. Horizontal eye movement sensitivities to the interaural acceleration sensed by the otoliths during translation (black) and roll tilt (gray) are compared in D (in similar units of deg·s–1·g–1). Model predictions in B and C are superimposed on experimentally observed horizontal (black circles) and torsional (gray triangles) ocular responses in rhesus monkeys (data in B: ±22°, 0.1–0.5 Hz, and ±5° at 1-Hz peak roll amplitudes; Angelaki et al. 1999Go; data in C: 0.1 g at 0.16, 0.2 Hz, and 0.4 g, 0.3- to 5-Hz peak translational accelerations; Angelaki 1998Go). Illustrated frequency response plots represent analytical predictions that are equivalent to the simulated model behavior for all amplitudes of head movement in A and C. Theoretically predicted responses in B and D (solid black and gray curves) represent linearized approximations valid for small head reorientations relative to gravity (e.g., <30° where static approximations to {alpha}x(t) and {alpha}z(t) can be assumed; see MODEL DEVELOPMENT and APPENDIX). Gray squares and dash-dot curves in B describe response gains and phases obtained by model simulation for large head roll angles (±90°). Dashed black trace in B represents the horizontal gain predicted if the canals were perfect transducers of head velocity at all frequencies.

 
Although interaural acceleration-sensitive otolith afferents modulate similarly during both lateral translation and roll tilt, the horizontal ocular responses to these motions are quite different. Ambiguous sensory otolith information has therefore been appropriately combined with integrated semicircular canal signals in the model to ensure that dynamic head tilts do not give rise to a TVOR. This is explicitly illustrated in Fig. 4D where tilt and translation responses are both expressed relative to the acceleration stimulus; the horizontal eye velocity sensitivity to interaural linear acceleration during roll tilt is clearly attenuated compared to that during translation across all frequencies (Fig. 4D). Notice, however, that at low frequencies the predicted horizontal ocular responses to small amplitude roll rotations are nevertheless significantly larger than those reported experimentally in rhesus monkeys (compare black traces with black circles in Fig. 4B; Angelaki et al. 1999Go). The difference between model predictions and experimental data is much smaller when large-amplitude roll rotations are simulated (e.g., ±90° peak roll deviations, gray squares and dash–dot lines in Fig. 4B). This is the case because of the inherent nonlinearities associated with large reorientations relative to gravity.

The deterioration in model performance with decreasing frequency occurs because the semicircular canals cease to provide perfect estimates of head velocity. In fact, when the canals are assumed to be perfect transducers of head velocity [i.e., canal transfer function C(s) = 1] negligible horizontal responses to tilt are predicted across all frequencies (Fig. 4B, black dashed lines). In the absence of an accurate canal estimate of head velocity in the real physiological system, other strategies for distinguishing tilt from translation are required to achieve appropriate behavior at low frequencies (e.g., Mergner and Glasauer 1999Go; Paige and Tomko 1991aGo; Telford et al. 1997Go). In the following, we will focus on simulations of cell responses at frequencies above 0.1 Hz, where the semicircular canals provide reliable estimates of head velocity and the network appropriately discriminates tilts and translations.

Simulated behavioral responses to tilt–translation combinations.    Novel combinations of translational and roll tilt movements have recently been used to investigate semicircular canal and otolith contributions to oculomotor responses (Angelaki et al. 1999Go; Green and Angelaki 2003Go). Similar stimulus combinations were used to simulate the performance of the model. Four protocols are illustrated at the top of Fig. 5A that consist of either lateral translation (Translation only), roll tilt (Roll tilt only), or combined lateral translation and roll tilt motion stimuli (Roll tilt + Translation motion and Roll tilt Translation motion). Because the interaural acceleration ({alpha}y) stimulus to the otoliths was matched for Translation and Roll tilt motions (each set to a peak of 0.2 g at 0.5 Hz; Fig. 5A, bottom row), combined motion stimuli result either in a doubling of the interaural acceleration stimulus (Roll tilt + Translation) or zero acceleration (Roll tiltTranslation), depending on the relative directions of the two stimuli.



View larger version (40K):
[in this window]
[in a new window]
 
FIG. 5. Simulated behavioral responses. A: horizontal and torsional eye velocity responses (hor, solid black; tor, dashed gray) to 4 tilt–translation stimulus paradigms at 0.5 Hz (i.e., Translation only, Roll tilt only, Roll tilt + Translation, and Roll tiltTranslation; Angelaki et al. 1999Go; Green and Angelaki 2003Go). B: ocular responses during Supine Roll, Supine Translation, Supine Yaw, and Upright Yaw motions. Bottom 2 traces in A and B illustrate the stimuli: roll and horizontal angular velocities ({omega}x, solid black; {omega}z, dashed gray lines); net and translational components of the interaural acceleration ({alpha}y, solid black; fy, dashed gray lines).

 
Simulations of the model output for the sinusoidal tilt–translation protocols are illustrated in Fig. 5A. Horizontal eye velocities (TVOR) of similar amplitude are predicted during all movements that include a translational component (Fig. 5A, solid black traces, top row). In addition, torsional eye velocity (RVOR) is elicited whenever the movement includes a roll component (Fig. 5A, dashed gray traces, top row). In contrast, no horizontal TVOR is predicted during pure roll tilt, despite the presence of an identical interaural acceleration stimulus to the otoliths (Fig. 5A, compare columns 1 and 2). Notice that an appropriate TVOR is elicited even in the absence of an interaural acceleration stimulus to the otoliths during Roll tiltTranslation motion. These predictions are compatible with recent observations in rhesus monkeys (Angelaki et al. 1999Go; Green and Angelaki 2003Go).

The model can also predict appropriate eye movement responses with the head in supine orientation (Fig. 5B). For example, roll rotation and lateral translation elicit torsional and horizontal responses, respectively (Fig. 5B, columns 1 and 2). More interesting is the case of supine yaw rotation that dynamically stimulates the otoliths along the interaural head axis (Fig. 5B, column 3). This stimulus condition presents a similar ambiguity problem to the case of roll tilt from upright. Specifically, because the otoliths are dynamically stimulated during supine yaw rotation, horizontal eye velocity could potentially reflect a combination of TVOR and RVOR response components. If otolith signals were not appropriately interpreted by the brain, supine yaw rotation would elicit significantly smaller horizontal responses than upright yaw rotation, whereas larger horizontal eye velocities would be predicted in the prone orientation. However, model simulations result in identical horizontal eye movements during both upright and supine yaw rotations (Fig. 5B; compare columns 3 and 4). Thus, just as in the case of roll tilt from upright, interaural accelerations attributed to a head reorientation relative to gravity during supine yaw rotation are appropriately distinguished from translation.

Neural response predictions and simulations

Given a model that predicts appropriate behavioral responses, we may now address the predicted response properties of different average neural populations within such a network. First, we will illustrate the frequency response predictions for each cell type during pure roll tilts and translations, as well as their responses to earth-vertical–axis rotations in different pitch head orientations when the canals are stimulated in isolation. We will show that in the integrative network proposed here, traditional interpretations of the responses to these stimuli embed assumptions that can lead to incorrect conclusions with respect to the signals encoded by central neurons. The goals of this section will be to illustrate this point and then to consider experimental protocols appropriate to reveal the underlying properties of the network.

Frequency response and earth-vertical–axis rotation predictions.    Neurons VO1 and VO2 exhibit the expected characteristics of semicircular-canal–sensitive cells, coding for head rotation in head coordinates (Fig. 6). In both upright and supine orientations, cell VO1 modulates in phase with angular yaw velocity but does not respond during roll rotation (Fig. 6, B and C, black traces), whereas cell VO2 responds exclusively to roll rotations (Fig. 6, B and C, gray traces). Neither cell group modulates during a pure translational stimulus (Fig. 6A). The canal afferent-like behavior of cells VO1 and VO2 also holds for earth-vertical–axis rotations at different static pitch orientations (Fig. 6, D and E). Cell VO1 responds maximally during earth-vertical–axis rotation in the upright orientation (i.e., during yaw rotation for a pitch angle = 0°), exhibiting a sensitivity that drops off with the cosine of pitch angle from upright, as predicted for a dominantly horizontal canal-sensitive cell. Similarly, cell VO2 exhibits no modulation during earth-vertical–axis rotation with the head upright and maximal responses in the prone and supine positions (i.e., pitch angles of ±90°) consistent with this cell receiving mainly vertical semicircular canal inputs (Fig. 6E). Accordingly, the properties of VO1 and VO2 are consistent with cells described as "canal-only" neurons (Dickman and Angelaki 2002Go).



View larger version (32K):
[in this window]
[in a new window]
 
FIG. 6. Predicted responses for cell populations VO1 and VO2. AC: frequency response predictions for cells VO1 (black) and VO2 (gray) during (A) translation, (B) roll rotation, and (C) yaw rotation in upright (solid curves) and supine (dashed curves) orientations. Cell response sensitivities and phases are expressed relative to head velocity. DE: predicted responses for cells (D) VO1 and (E) VO2 during earth-vertical–axis (EVA) rotations at 0.5 Hz in different static pitch orientations (0°, upright; 90°, prone; –90°, supine).

 
Cells VO3 and VO4 reflect a more complicated processing of sensory information. These neural populations exhibit strong modulations during high frequency translation, in phase with either ipsilaterally or contralaterally directed linear acceleration in both upright and supine head orientations (i.e., phase of ±90° relative to head velocity; Fig. 7A). Their activities thus reflect high-pass responses to translational velocity. Although the robust translational responses of cells VO3 and VO4 clearly demonstrate an otolithic sensory contribution to their activities, there is no evidence for a canal contribution during earth-vertical–axis rotations (Fig. 7, D and E). Nonetheless, cell VO4 exhibits a robust response at high frequencies during roll tilt from upright or yaw tilt from supine orientations (solid gray curves in Fig. 7B and dashed gray curves in Fig. 7C, respectively) that is in phase with angular head velocity (i.e., phase of –180° during upright roll and 0° during supine yaw). In contrast, cell VO3 exhibits little modulation during high frequency head reorientations relative to gravity despite the presence of a gravitational stimulus to the otoliths. Thus, as predicted, cell VO3 distinguishes between inertial and gravitational accelerations, coding selectively for translational accelerations. The response of this model cell population is compatible with those of cells in the vestibular and fastigial nuclei that have been shown to preferentially encode the translational component of linear acceleration (Angelaki et al. 2003Go).



View larger version (34K):
[in this window]
[in a new window]
 
FIG. 7. Predicted responses for cell populations VO3 and VO4. AC: frequency response predictions for cells VO3 (black) and VO4 (gray) during (A) translation, (B) roll rotation, and (C) yaw rotation in upright (solid curves) and supine (dashed curves) orientations. Cell response sensitivities and phases are expressed relative to stimulus velocity. DE: predicted responses for cells (D) VO3 and (E) VO4 during EVA rotations at 0.5 Hz in different static pitch orientations (0°, upright; 90°, prone; –90°, supine).

 
In contrast to cells VO3 and VO4, VO5 does not respond during translation (Fig. 8A). It also exhibits no response during earth-vertical–axis rotations (not shown) but modulates in phase with angular position (–180° or 0° phase) during roll or yaw rotations when the rotation involves a dynamic head reorientation relative to gravity (Fig. 8, B and C). Thus in agreement with theoretical predictions (see MODEL DEVELOPMENT), cell VO5 encodes angular head position relative to gravity or tilt (more generally, gy; see APPENDIX). Although we focus here on responses at mid-high frequencies (>0.1 Hz) it is relevant to note that VO5 is also predicted to respond to static head tilt. Cells that respond robustly to both static and dynamic head tilts from upright but not to high-frequency translation have been observed in the vestibular nuclei (Zhou et al. 2000Go).



View larger version (27K):
[in this window]
[in a new window]
 
FIG. 8. Predicted responses for cell population VO5. AC: frequency response predictions during (A) translation, (B) roll rotation, and (C) yaw rotation in upright (solid curves) and supine (dashed curves) head orientations. Cell response sensitivities are expressed relative to translational velocity (A) or angular head position (BC).

 
In summary, the idealized model predicts several average cell populations that are compatible with those observed experimentally, including cells that respond exclusively to canal stimulation (i.e., cells VO1 and VO2) as well as cells that preferentially encode either head translation (cell VO3) or tilt (cell VO5). Yet, using the stimuli described in Fig. 6 through 8, incorrect conclusions would likely be reached regarding the sensory contributions to their responses. In particular, because cells VO3, VO4, and VO5 fail to modulate during earth-vertical–axis rotations, it might be concluded that their activities do not reflect canal signal contributions (a conclusion that is incorrect). Neither can their responses be directly correlated with the sensory otolith signal because none of these cell groups encodes the net gravito-inertial acceleration. As a result, their responses during earth-horizontal–axis rotations (i.e., combined otolith and canal stimulation) do not reflect a simple summation of their activities during translation (estimated otolith contribution) and earth-vertical–axis rotations (estimated canal contribution).

This observation has in fact been made for many central translation-sensitive vestibular neurons (Dickman and Angelaki 2002Go) and was used to suggest that the activities of such neurons might reflect intermediate processing stages in distinguishing between head tilt and translation. Yet the question of how this distinction is made remained unanswered. The model proposed here provides an explanation. Specifically, the responses of cells VO3, VO4, and VO5 differ from those of sensory otolith signals because their activities do indeed reflect the contribution of semicircular canal inputs (Fig. 3A). However, these cells do not code for rotation in the head-fixed reference frame of the canal sensors. Because of the multiplicative canal–otolith interactions at the input to cell VO4, the activities of these cells instead reflect spatially referenced canal signal contributions aligned with the earth-horizontal axis (e.g., vertical canal signals during tilts from upright and horizontal signals during tilts from supine head orientations). They thus encode only the component of rotation orthogonal to gravity that is not observed during earth-vertical–axis rotations. Such a postulated canal-derived signal would nevertheless be difficult to detect because it contributes only under conditions that typically simultaneously stimulate the otoliths. In the following section, we will address how the proposed "hidden" semicircular canal signal contribution to the responses of these neurons can be isolated.

Simulated responses to combined tilt–translation stimuli.    To isolate the contribution of semicircular canal signals to the neural activities of cells VO3, VO4, and VO5 we may investigate the simulated responses of the model neurons for the same set of stimulus combinations employed to examine behavioral responses (Fig. 9). As expected based on the analytical predictions in Figures 6 through 8, cells VO1 and VO2 code for horizontal and roll head velocity, respectively, regardless of head orientation (Fig. 9, A and B). In contrast, VO3 exhibits negligible responses to roll and yaw rotations, in both upright and supine head orientations (e.g., Fig. 9C, columns 2 and 5) but encodes the translational component of acceleration (Fig. 9C, similar responses in columns 1, 3, and 4). Cell VO5, on the other hand, codes specifically for head tilt relative to gravity, responding during either upright roll (Fig. 9E, columns 2, 3, and 4) or supine yaw (Fig. 9E, column 5) but not during pure translation. Finally, cell VO4 exhibits a more complex behavior, being sensitive to both translation and head reorientation relative to gravity.



View larger version (41K):
[in this window]
[in a new window]
 
FIG. 9. Simulations of model cell responses. AE: firing rate modulation of cells VO1 through VO5 during the 4 tilt–translation protocols in the upright orientation (columns 14; see Fig. 5 legend and text) and for supine yaw rotation (column 5). F: angular roll and yaw velocities ({omega}x, solid black; {omega}z, dashed gray lines). G: net and translational interaural accelerations ({alpha}y, solid black; fy, dashed gray lines).

 
Although earth-vertical–axis rotations do not reveal the contribution of semicircular canal signal inputs to cells VO3, VO4, and VO5 (e.g., Fig. 7, D and E), notice that all 3 cell populations modulate with sinusoidal responses during Roll tilt Translation motion (Fig. 9, column 4). This observation is relevant because during this unique protocol there is no net linear acceleration along the interaural head axis and thus no sensory stimulus along this axis to the otoliths (Fig. 9G, column 4). Therefore, the modulation in the firing rates of these cells must be attributed to semicircular-canal–derived signals. Specifically, during the upright Roll tilt Translation protocol, VO4 cell responses are entirely driven by vertical semicircular canal signals. These signals are subsequently integrated by the network to produce the head tilt position responses of cell VO5. It is a scaled version of this same integrated canal signal that appears on cell VO3 (Fig. 9C, column 4) and is subsequently used to drive compensatory responses to translation during Roll tiltTranslation motion, even in the absence of an appropriate dynamic otolith stimulus. Thus, the Roll tilt Translation motion paradigm unmasks a semicircular canal contribution to the activities of these cells that might otherwise remain "hidden" during earth-vertical–axis rotations in all head orientations.

Variability in individual cell responses

Using a particular model structure and parameter set chosen to achieve ideal tilt–translation discrimination we have illustrated the response properties of several average cell populations within this network. However, in contrast to these idealized model cells, the majority of experimentally observed translation-sensitive neurons do in fact exhibit responses to earth-vertical–axis rotations that reflect the sensory contribution of signals from multiple orthogonal canals (Dickman and Angelaki 2002Go). Furthermore, relatively few cells have been isolated that respond exclusively to tilt or translation (Angelaki et al. 2003Go; Zhou et al. 2000Go), suggesting that these variables may be encoded mainly as population averages (i.e., as represented by our average model neurons).

In the following sections we will explore the effects of varying particular parametric assumptions in the model with two key goals: 1) to examine a more realistic representation of the properties of individual neurons that contribute to the average population responses; 2) to further explore the implications of the most fundamental properties of the proposed model that must be shared by any neural network that effectively implements Eq. 4. These include its function as a neural integrator and the requirement for a coordinate transformation of canal signal contributions. To illustrate these issues we will consider how changes in the coupling of semicircular canal and otolith signals onto the network impact on the expected properties of individual neurons and their compatibility with experimental observations.

Variable semicircular canal-related projection weights.    To achieve close to ideal tilt–translation discrimination in the proposed model we assumed that the net projections from each canal onto cell VO4 (i.e., indirect projections from cells VO1 and VO2) were equal in strength and entirely head-orientation–dependent (i.e., K1a = K2a and K1o = K2o = 0 in Fig. 10A and Fig. 3A; see Fig. 3 legend). Average cell populations VO3, VO4, and VO5 were then predicted to exhibit no response during earth-vertical–axis rotations (e.g., solid black trace with zero gain in Fig. 10B; see also Fig. 7, D and E). Using the VO4 cell population as an example, we will now examine the effect of relaxing these parametric assumptions to explore the expected range of responses from individual neurons within this population.



View larger version (22K):
[in this window]
[in a new window]
 
FIG. 10. Predicted EVA rotation responses for cell population VO4, allowing for potential variability in individual cell characteristics. A: schematic of a cell (part of the VO4 population, Fig. 3A) that receives both head-orientation–dependent and head-orientation–independent horizontal and vertical canal inputs (weights K1a, K2a and K1o, K2o, respectively). Head-orientation–dependent inputs are shown to be multiplied by nasoocciptal and dorsoventral accelerations, {alpha}x and {alpha}z. B: cell VO4 exhibits no response during EVA rotations (solid black line, zero gain) when K1o = K2o = 0 and K1a = K2a. When K1o = K2o = 0.2, the cell exhibits a response consistent with equal sensitivities to horizontal and vertical canal stimuli (solid gray line). For comparison, the response predicted if all inputs to the cell were instead head-orientation–invariant (e.g., K1a = K2a = 0, K1o = 1.2, K2o = 1.2) is illustrated with the dotted line. C: responses predicted when horizontal and vertical canal head-orientation–dependent projections have unequal strengths (solid black line: K1a = 2, K2a = 1 or K1a = 1, K2a = 2 and K1o = K2o = 0; solid gray line: K1a = 3, K2a = 1 or K1a = 1, K2a = 3 and K1o = K2o = 0). For comparison, the response predicted from a neuron with a head-orientation–independent sensitivity only to vertical canal input (K2o = 1, K1o = K1a = K2a = 0) is illustrated with the dotted line. D: response patterns observed for different combinations of head-orientation–independent and -dependent projection strengths. Solid black: K1a = 2, K2a = 1, K1o = K2o = 0; solid gray: K1a = 2, K2a = 1, K1o = 0.4, K2o = 0; dashed black: K1a = 1, K2a = 2, K1o = 0.2, K2o = 0; dashed gray: K1a = 1, K2a = 2, K1o = 0, K2o = –1; dotted gray: K1a = 2, K2a = 1, K1o = 2, K2o = 0. Phases (not shown) are always either 0° or 180°.

 
We may begin by examining what happens if there is a head-orientation–independent canal component to an individual VO4 cell's response (i.e., K1o != 0 and/or K2o != 0). For example, in addition to head-orientation–dependent projections (still assumed to be equal at this point), a VO4 neuron could also receive small head-orientation–invariant horizontal and vertical canal contributions (K1o = 0.2, K2o = 0.2). Under these conditions, the cell would exhibit the cosine-type tuning during earth-vertical–axis rotations expected for a neuron equally sensitive to horizontal and vertical canal inputs, thus demonstrating evidence for orthogonal canal–canal convergence (Fig. 10B, solid gray curve). The tuning exhibited by the cell is nonetheless clearly different from the responses that would be observed if all sensory canal projections were invariant to head orientation (Fig. 10B; compare dotted black and solid gray curves). Although the cell would no longer be classified as an "otolith-only" neuron, a significant component of the canal contribution to its response would remain hidden unless somehow explicitly unmasked (e.g., during Roll tilt Translation motion).

More fundamental to the arguments here are the predictions made when the assumption of equal orientation-dependent horizontal and vertical canal-related projections is relaxed (i.e., K1a != K2a). In this case, the predicted responses no longer reflect simple cosine tuning patterns, but rather exhibit second harmonic spatial tuning properties (e.g., solid black and gray curves in Fig. 10C) indicative of head-orientation–dependent canal signal contributions (see APPENDIX for details). Notably, such spatial tuning might not be apparent for small pitch angles, given that these curves could appear similar to those of a neuron that receives exclusively orientation-independent vertical canal signals (Fig. 10C, compare solid black and dotted traces). Examination of cell response properties over a large range of static tilt angles (e.g., ±90°) is therefore likely to be necessary to reveal the presence of head-orientation–dependent rotational sensitivities.

More generally, individual neurons are likely to receive different combinations of orientation-dependent and -independent canal projections, giving rise to a range of response patterns during earth-vertical–axis rotations that reflect different degrees of convergence from multiple orthogonal canal sensors, as observed experimentally (Fig. 10D; Dickman and Angelaki 2002Go; Siebold et al. 2001Go). When examined over a larger range of head orientations than those typically used (≤30°; e.g., Dickman and Angelaki 2002Go; Siebold et al. 2001Go), however, their responses are expected to differ considerably from the simple cosine-tuned behavior that has traditionally been assumed for vestibular neurons. Such response patterns do not simply reflect the particular model structure chosen here but would be expected in any network in which head-orientation–dependent canal- and otolith-derived signals converge to distinguish tilts and translations according to the requirements implied by Eq. 2 and 4. Hence, although these complex tuning properties have yet to be observed experimentally, they represent a fundamental model prediction that remains to be tested.

Variable otolith signal projection weights.    Otolith signals couple onto the model network not only through cell VO3, but also through sensory projections directly onto cell population VO4 (projection with weight q2) and into the neural filter CLP(s) (projection with weight q3). The weights of these parameters were chosen both to set particular cell sensitivities to translation and to ensure that cell VO3 exhibits high-pass filtered responses to otolith stimulation with a minimal response to static head tilt (see APPENDIX). With the chosen parameter set the average neural populations responsive to translation (i.e., cells VO3 and VO4) were predicted to modulate closely in phase with head acceleration (Fig. 7). However, a consistent, yet so far unexplained, experimental observation is that eye-movement–insensitive vestibular neurons exhibit dynamic responses to translation that are highly variable and often quite different from those of otolith afferents (Angelaki and Dickman 2000Go; Chen-Huang and Peterson 2002Go; Dickman and Angelaki, 2002Go; Musallam and Tomlinson 2002Go; Tomlinson et al. 1996Go). Next we will illustrate that by changing a single parameter, weight q3, individual neurons with a wide range of dynamic responses to head translation are predicted in the proposed model. Changes in this single parameter are sufficient to illustrate a large range of response gains and phases during translation without affecting the integrative properties of the network.

Figure 11 illustrates distributions of the predicted gains and phases (relative to translational acceleration) of neural populations VO3, VO4, and VO5 at 0.5 Hz when weight q3 was randomly varied according to a Gaussian distribution about its nominal value (see Fig. 11 legend). Notice that both populations VO3 and VO4 can exhibit a broad distribution of response phases reflecting either lags relative to ipsilaterally directed or leads relative to contralaterally directed acceleration (Fig. 11, A and B). This is particularly evident in the case of the VO3 population that exhibits a phase distribution peak more closely aligned with translational velocity (±90°) than acceleration (0° or –180°). In addition, whereas with the nominal parameter set cell VO5 was predicted to respond only to head tilt, with changes in weight q3 the cell may also exhibit responses to translation. Notice, however, that these translational responses are always in phase with either ipsilaterally or contralaterally directed translational velocity (±90°; Fig. 11C). These points are further emphasized in Fig. 11D where the response gains and phases for cells VO3, VO4, and VO5 are plotted as a function of frequency for an example case where q3 = 1. Each cell population now exhibits distinct dynamic characteristics such that at 0.5 Hz a range of response phases relative to translational acceleration are evident, consistent with experimental observations (Dickman and Angelaki, 2002Go). Notably, cell VO5 modulates in phase with translational velocity, exhibiting low-pass response characteristics relative to translational acceleration. This observation is consistent with the dynamic properties of cells shown to predominantly encode head tilt at high frequencies (Zhou et al. 2000Go).



View larger version (32K):
[in this window]
[in a new window]
 
FIG. 11. Predicted variability in cell response characteristics during translation. AC: histograms illustrating the distributions of 0.5-Hz response gain and phase predicted for individual cells contributing to populations (A) VO3, (B) VO4, and (C) VO5 when weight q3 was varied about its nominal mean value of 0.1975 according to a normal distribution with SD equal to 10 times this parameter (i.e., {approx}68% of the randomly chosen q3 values fell within ±10 times the mean value of 0.1975). Each histogram displays a total of 1,000 gains or phases, associated with different values of q3, that are grouped into 100 bins. Gray stars denote the mean gain and phase values of each cell type predicted for the nominal parameter set used in all simulations (see Fig. 3 legend). Black stars indicate the population means (calculated vectorially; e.g., Wei and Angelaki 2001Go) when q3 was varied about its nominal value. Gray and black stars superimpose in C, indicating no difference in the population and nominal parameter set mean values. Note, however, that even though q3 was varied according to a normal distribution about its nominal value, large changes in average cell responses are nonetheless apparent in A and B because gains and phases in the network are generally nonlinear functions of the model parameters. D: frequency response predictions for cells VO3 (dashed), VO4 (dotted), and VO5 (solid) for q3 = 1. In contrast to Figs. 68, gains and phases are plotted here relative to translational acceleration.

 
It should be noted that because the characteristics of cell population VO3 are quite sensitive to changes in q3, varying this parameter in isolation alters the VO3 population response and hence the average model behavior. This is indicated by the stars in Fig. 11A, which compare the mean gains and phases of cell population VO3 using the nominally chosen parameter set (gray stars) with the average response characteristics after varying parameter q3 (black stars). Behavioral performance identical to that illustrated previously is nonetheless possible if inputs to the premotor TVOR pathways (i.e., onto cell EMI in Fig. 3B) are presumed to include a weighted combination of direct sensory otolith projections in addition to inputs from cell VO3 (assuming also appropriate adjustments in canal projection weights onto network NI1). Indeed, a contribution of direct sensory projections to the TVOR is compatible with the neuroanatomical observation that the shortest latency, although weak, sensory otolith projections to the abducens are mono- and disynaptic (Uchino et al. 1994Go).

Most relevant to the investigation here is the fact that the framework provides a functional rationale for the wide range of central response dynamics observed during translation (Angelaki and Dickman 2000Go; Chen-Huang and Peterson 2002Go; Dickman and Angelaki 2002Go; Musallam and Tomlinson 2002Go; Tomlinson et al. 1996Go; Zhou et al. 2000Go, 2001Go). Although such dynamic characteristics could appear consistent with a low-pass versus high-pass filtering strategy for distinguishing tilts and translations (Paige et al. 1991aGo; Telford et al. 1997Go), they arise here as the result of coupling otolith sensory signals to an integrative network that combines otolith and canal sensory information to solve the differential equations necessary to detect head orientation in space.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DEVELOPMENT
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The detection of motion in a gravitational field poses particular computational demands on the central nervous system because inertial and gravitational linear accelerations are sensed equivalently by the otoliths. To accurately distinguish between the two requires an estimate of changes in head orientation, implying the requirement for an integration of information from multiple sensory organs (Guedry 1974Go; Mayne 1974Go; Young 1984Go). Indeed, recent experimental studies of perception and oculomotor behavior have confirmed not only that tilts and translations are appropriately discriminated in the frequency range where the canals provide accurate estimates of head velocity, but that information from these angular sensors (in the absence of nonvestibular cues) is essential to resolve the ambiguity (Angelaki et al. 1999Go; Green and Angelaki 2003Go; Merfeld et al. 2003Go; also see Merfeld et al. 1999Go). Furthermore, both experimental observations (Green and Angelaki 2003Go) and theoretical predictions (Angelaki et al. 1999Go; Bos and Bles 2002Go; Glasauer and Merfeld 1997Go; Hess and Angelaki 1997Go; Merfeld 1995Go; Merfeld and Zupan 2002Go; Mergner and Glasauer 1999Go; Viéville and Faugeras 1990Go; Zupan et al. 2002Go) imply that resolution of the problem involves a temporal integration of angular velocity information.

Here, we have illustrated an integrative network that can compute an internal estimate of head orientation relative to gravity along the interaural head axis at mid-high frequencies (e.g., >0.1 Hz). In conjunction with net linear acceleration signals provided by the otoliths, the translational acceleration component can then be extracted to provide an appropriate drive to motor and perceptual systems. Our investigation here is novel in that we explore the problem in a physiologically relevant (albeit limited) context, where the scalar firing rate activities of different populations of neurons can be predicted. We illustrate that the computations for discriminating tilt and translation can be performed by an integrative network in which the strengths of canal signal projections modulate nonlinearly (i.e., in a multiplicative sense) as a function of current head orientation. The most important feature of the approach is that it implies the existence of cells that respond to rotational stimuli from the canals differently depending on head orientation. As described below, this observation has important implications both for the characterization of cell responses and for identifying their contributions to spatial motion estimation.

Comparison of model predictions with experimentally observed neural response properties

The properties of neural populations likely to be involved in inertial motion estimation have been characterized during rotation and/or translation in several brain areas including the vestibular nuclei (Angelaki and Dickman 2000Go; Angelaki et al. 1993Go; Baker et al. 1984aGo,1984bGo; Brettler and Baker 2001Go; Bush et al. 1993Go; Chen-Huang and Peterson 2002Go; Dickman and Angelaki 2002Go; Endo et al. 1995Go; Fukushima et al. 1990Go; Graf et al. 1993Go; Iwamoto et al. 1996Go; Kasper et al. 1988Go; Musallam and Tomlinson 2002Go; Perlmutter et al. 1998Go, 1999Go; Schor et al. 1984Go, 1985Go; Wilson et al. 1990Go, 1996Go; Yakushin et al. 1999Go) and the rostral fastigial nucleus of the cerebellum (Siebold et al. 1997Go, 1999Go, 2001Go; Zhou et al. 2001Go). However, despite extensive characterization of the spatial tuning characteristics of such neurons, in many studies otolith and canal signal contributions were not adequately distinguished (see below). Furthermore, the neural correlates for the computations underlying tilt–translation discrimination are only now beginning to be investigated. Thus, the data with which our model predictions can be directly compared are currently limited. Using one potential model structure that can perform the computations implied by Eq. 2 and 6 to distinguish tilts and translations, we have nonetheless illustrated the ability to predict distinct average cell populations that are compatible with those observed experimentally to date. These include neurons that respond mainly to head rotation ("canal-only" cells) with spatial response selectivities closely aligned with a particular canal plane (e.g., cells VO1 and VO2; Dickman and Angelaki 2002Go) as well as neurons that code explicitly for internal dynamic estimates of either translation (cell VO3; Angelaki et al. 2003Go) or tilt (cell VO5; Zhou et al. 2000Go).

Most central translation-sensitive neurons, however, encode neither exclusively translation nor tilt (Angelaki et al. 2003Go). In fact, the most consistent experimental observation is a wide distribution of response properties that are difficult to interpret (Angelaki and Dickman 2000Go; Chen-Huang and Peterson 2002Go; Dickman and Angelaki 2002Go). As will be elaborated on below, two particular aspects of the framework explored here are not only compatible with general characteristics of experimentally observed responses, but lend particular insight as to the strategies likely used by the brain to estimate spatial orientation and motion.

A first general experimental observation is that the dynamic characteristics of central responses to translation are highly variable and typically quite different from those of sensory otolith afferents (Angelaki and Dickman 2000Go; Chen-Huang and Peterson 2002Go; Dickman and Angelaki 2002Go; Musallam and Tomlinson 2002Go; Tomlinson et al. 1996Go; Zhou et al. 2001Go). Cells with a wide range of response phases relative to the sensory stimulus have been observed during 0.5-Hz translations, and many modulated more closely in phase with translational velocity than acceleration (Angelaki and Dickman 2000Go; Dickman and Angelaki 2002Go). Although the contribution of an integrative element to their response dynamics was apparent (Angelaki and Dickman, 2000Go; Dickman and Angelaki 2002Go; also see Musallam and Tomlinson 2002Go) its origin and functional role remained unexplained. The current model provides a potential functional rationale for this observation. Both a variability in response phase and a tendency toward encoding head velocity at midrange frequencies (e.g., 0.5 Hz) are not only predicted but in fact almost inevitable consequences of coupling otolith sensory signals to an integrative network for detecting head orientation (e.g., Fig. 11).

A second general observation is that neurons sensitive to dynamic otolith stimulation in the vestibular and fastigial nuclei exhibit either no response during earth-vertical–axis rotations or complex spatial tuning properties with characteristics of orthogonal canal–canal convergence (Dickman and Angelaki 2002Go; Siebold et al. 2001Go). A functional significance for this experimental observation is also provided by the framework proposed here. We have illustrated that the ability to distinguish tilts and translations at mid-high frequencies implies a convergence of otolith information with sensory signals from multiple sets of canals (e.g., Figs. 3A and 10). Furthermore, because the problem of tilt–translation discrimination also implies a head-orientation–dependent coupling between canal and otolith signals, many central vestibular neurons may not encode rotational signals in the head-fixed reference frame of the canal sensors. As a result, any component of their responses aligned with an earth-horizontal axis would not be isolated during earth-vertical–axis rotations. This can account for the experimental observation that their responses during earth-horizontal–axis rotations (i.e., combined otolith and canal stimulation) often cannot be predicted by a vector summation of their activities during translation (estimated otolith contribution) and earth-vertical–axis rotations (estimated canal contribution; Dickman and Angelaki 2002Go).

In the following section, we consider how these two key features of the proposed model (i.e., its integrative properties and multiplicative canal–otolith interactions) expected to be common to any network that effectively solves Eq. 2 and 4, have an impact on the interpretation of neural activities and the protocols required to identify their functional contributions to the problem of distinguishing inertial and gravitational accelerations.

Implications for investigating neural response properties

The requirements for discriminating tilt and translation are relatively straightforward from a theoretical standpoint. However, the task of identifying how the brain actually uses the available sensory information to perform the required computations represents a much greater challenge because rotations that reorient the head relative to gravity stimulate both the canals and the otoliths simultaneously. Key to addressing how the required processing takes place is the question of how the contributions from each sensor should be identified.

Most studies to date have characterized cell responses only during purely rotational movements. The extent of canal–otolith convergence was often then estimated based on the assumption that the contributions from each sensor would mimic the dynamic characteristics of the corresponding afferent population. Specifically, whereas canal afferents encode angular velocity, otolith afferents modulate in phase with linear acceleration, a signal that appears proportional to angular head position during small angular reorientations relative to gravity (e.g., Fig. 12A, top). Experimental observations of large changes in response phase, including spatiotemporal convergence (STC) properties, have thus typically been attributed to a convergence of canal and otolith signals (Baker 1984bGo; Endo et al. 1995Go; Kasper et al. 1988Go; Perlmutter et al. 1998Go, 1999Go; Siebold et al. 1997Go, 1999Go; Wilson et al. 1990Go). Notably, however, studies in which the otoliths were stimulated in isolation have illustrated that such an assumption is invalid (Angelaki and Dickman 2000Go; Bush et al. 1993Go; Dickman and Angelaki 2002Go; Schor et al. 1985Go). This is also suggested by the theoretical investigation here.



View larger version (24K):
[in this window]
[in a new window]
 
FIG. 12. Implications of the proposed framework for distinguishing canal and otolith sensory contributions to neural responses. A: in the absence of dynamic processing, canal and otolith contributions to central neural responses mimic the dynamic characteristics of the corresponding afferent population (top). Canal contributions are in phase with angular velocity {omega}, whereas otolith contributions are in phase with linear acceleration {alpha}, or equivalently angular head position {theta}, for small tilt angles {theta}. If sensory signals are dynamically processed by an integrative network, the dynamic contributions of the respective sensors will be reversed (bottom): integrated canal signals modulate in phase with angular position {theta}. Integrated otolith signals will appear in phase with linear velocity v, during translation or in phase with angular velocity {omega}, at a given frequency during dynamic tilt. B: central vestibular neurons are typically assumed to encode angular velocity in the head-fixed coordinate frame of the semicircular canals ({omega}HEAD-FIXED; top). The multiplicative canal–otolith interactions implied by Eq. 6 to estimate dynamic head orientation relative to gravity imply a spatially referenced internal estimate of angular velocity aligned with the earth-horizontal axis ({omega}SPACE-FIXED,EARTH-HORIZONTAL; middle). More generally, with imbalances in the strengths of projections derived from different canals (c1 != c2), individual neurons may encode a representation of angular velocity that is neither entirely space-fixed, nor head-fixed ({omega}INTERMEDIATE; bottom). This is predicted to be reflected in second-harmonic tuning characteristics during EVA rotations (Fig. 10).

 
A key feature of the proposed model is that it performs a temporal integration of semicircular canal signals to construct an internal estimate of dynamic head tilt. If both the otoliths and canals converge on such a network, the respective contributions from each sensor may reflect signals in phase with either head position or velocity at various points within the circuit. As summarized in Fig. 12A, whereas the otolith- and canal-related signal contributions of cells at the input of such an integrative network would closely reflect the dynamic characteristics of the respective sensory signals, cells at the output end of such an integrator (e.g., cell VO5) would reflect apparently "opposite" dynamic characteristics (e.g., Fig. 12A, bottom). For example, the tilt responses of cell VO5 in phase with angular head position at mid-high frequencies arise almost exclusively from semicircular canal signal contributions, whereas any otolith sensory contribution to cell population VO5 is predicted to modulate in phase with translational velocity (e.g., Fig. 12A, bottom, Fig. 11, C and D). Other cells may reflect sensory signals that have been dynamically processed to different extents. Clearly, therefore, response dynamics should not be used to dissociate otolith versus canal sensory contributions at central stages of vestibular processing.

Recent studies of central vestibular neurons have alternatively examined cell responses during earth-vertical–axis rotations to isolate the canal contributions to cell responses in the absence of a dynamic stimulus to the otoliths (Dickman and Angelaki 2002Go; Siebold et al. 2001Go). A key assumption in employing this approach, however, is that central sensitivities to canal signals are invariant to head orientation (e.g., Fig. 12B, top). Again, the present investigation emphasizes that this approach may be invalid. The vector cross-product computations required to solve Eq. 4 imply that the brain constructs spatially referenced central estimates of head velocity aligned with an earth-horizontal axis (e.g., Fig. 12B, middle). In the proposed model this transformation occurs on cell VO4 with the result that particular cell populations (i.e., neurons VO3, VO4, and VO5) code mainly for the component of rotation orthogonal to gravity. Thus, their activities reflect a canal-derived signal that can be observed only under conditions that typically simultaneously activate the otolith organs. This implies that previous characterizations of the 3D central neural sensitivities to canal stimulation during earth-vertical–axis rotations may have been inaccurate, potentially leading to inaccurate estimates of dynamic otolith contributions.

Validation of the hypothesized existence of cells that encode head velocity at least partially in a spatially referenced coordinate frame is fundamental to understanding how the brain constructs central estimates of head orientation and motion. Here we have illustrated that a "hidden" canal component to cell activities can be unmasked when canal and otolith sensory contributions are combined to cancel the net interaural acceleration stimulus to the otoliths (i.e., during Roll tiltTranslation motion; Angelaki et al. 1999Go; Green and Angelaki 2003Go). However, to verify the existence of a head-orientation–dependent coding of canal-derived signals requires an examination of central responses in different head orientations. Specifically, a cell that explicitly uses a combination of otolith and spatially referenced canal-derived rotational signals to construct an internal estimate of translation must respond to Roll tilt Translation motion when upright but Yaw tiltTranslation motion in supine or prone orientations. Alternatively, because individual cells involved in the required computations are unlikely to exclusively encode the component of rotation orthogonal to gravity (i.e., earth-horizontal component) an examination of central responses during earth-vertical–axis rotations at different static head orientations may be sufficient to reveal head-orientation–dependent canal signal contributions. In particular, we have illustrated that these could be identified by a significant second harmonic component in the cell's spatial tuning properties during earth-vertical–axis rotations (Fig. 12B, bottom; Fig. 10; also see APPENDIX). However, to validate this model prediction, cell responses must be characterized over a broad range of head orientations (e.g., ±90°; Yakushin et al. 1999Go) much larger than those that have typically been used to date (≤30°; e.g., Dickman and Angelaki 2002Go; Siebold et al. 2001Go).

Relationship to velocity storage

We previously speculated that the network involved in constructing internal estimates of head orientation and motion might be that known as the "velocity storage integrator" (Green and Angelaki 2003Go). In particular, both the problem of tilt–translation discrimination and "velocity storage" behavior require the presence of central integrative networks. At least a subset of the central cells in each network are thus expected to reflect these integrative properties by exhibiting improved responses to low frequency rotations compared to those predicted from the dynamic characteristics of the semicircular canals (e.g., Fig. 13). Such extended low frequency rotational dynamics are prevalent in the activities of vestibular nuclei neurons (Buettner et al. 1978Go; Dickman and Angelaki 2003Go; Reisine and Raphan 1992Go; Waespe and Henn 1977Go, 1978Go, 1979Go). However, it remains to be addressed whether subsets of these same neurons participate in the computations surrounding tilt–translation discrimination.



View larger version (19K):
[in this window]
[in a new window]
 
FIG. 13. "Velocity-storage" type properties of cell VO4. Simulated response of cell VO4 (solid line) to a 100 deg/s step in head velocity (dotted line) when the cell is assumed to receive a small head-orientation–independent canal signal (K1o = 0.1). Cell's response reflects extended low-frequency dynamics compared to those of the canal afferent population (dashed line). Head velocity and simulated canal afferent traces have been scaled for ease of comparison with the predicted cell response.

 
A second link between the proposed integrative network and that traditionally associated with velocity storage relates to their 3-dimensional properties. It has been shown that the ocular nystagmus elicited by constant velocity rotation of either the subject or the visual surround tends to align with gravity or the net gravito-inertial vector (Angelaki and Hess 1994Go; Dai et al. 1991Go; Harris 1987Go; Merfeld and Young 1995Go; Merfeld et al. 1993aGo; Raphan et al. 1992Go; Wearne et al. 1999Go). This observation has led to the proposal that the velocity storage network constructs an estimate of rotational motion in a space-fixed (inertial) reference frame (Angelaki and Hess 1994Go, 1995Go; Angelaki et al. 1995Go). Yet, this is exactly what we have argued here is required to distinguish tilt from translation. Further insight as to the relationship between the two networks may be provided by examining how modifications of the spatial and/or dynamic characteristics of velocity storage, associated with lesions of the nodulus and ventral uvula (Angelaki and Hess 1995Go; Waespe et al. 1985Go; Wearne et al. 1998Go) or sectioning the commissural fibers in the region of the rostral medial vestibular nuclei (Katz et al. 1991Go; Wearne et al. 1997Go), impact on the ability to discriminate tilts and translations. Intimately related to this issue is the extension of the model to address the problem of inertial motion detection at low frequencies where the semicircular canals do not provide accurate estimates of head velocity and other strategies including frequency segregation (Paige and Tomko 1991aGo; Telford et al. 1997Go) and/or the use of otolith-based estimates of head rotation (Angelaki 1992aGo,1992bGo; Angelaki and Hess 1996Go; Correia and Money 1970Go; Hain 1986Go; Harris 1987Go; Kushiro et al. 2002Go; Raphan and Schnabolk 1988Go) may make a contribution. As illustrated in Fig. 2, extension of the framework to 3 dimensions implies coupling of interdependent integrative networks that is expected to lend further insight as to the relationship between tilt-translation discrimination and the complex 3D spatial characteristics of velocity storage (Angelaki and Hess 1994Go, 1995Go; Dai et al. 1991Go; Hess and Angelaki 1997Go; Raphan and Sturm 1991Go; Wearne et al. 1998Go, 1999Go).

In conclusion, estimation of head orientation and inertial motion in space is a fundamental task of the vestibular system, yet one that remains to be understood at the neural level. Here we propose, for the first time, a physiologically relevant implementation of the theoretical relationships derived from the physical laws of motion, necessary to distinguish tilts and translations. Average neural populations that encode internal neural estimates of dynamic translational and gravitational accelerations are shown to emerge within such a network. Yet, the high degree of variability and complex patterns of otolith–canal convergence observed in central neurons responsive to translation are not only difficult to interpret, but suggest a distributed central coding of such internal estimates. Here we illustrate that these seemingly complex properties are consistent with those of an integrative network that implements the nonlinear canal–otolith interactions required to distinguish inertial and gravitational accelerations. A key prediction of this model is the existence of central neurons that encode semicircular canal signals at least partially in a space-fixed coordinate frame, rather than in the traditionally assumed head-fixed reference frame defined by the peripheral vestibular system. New experimental paradigms and alternative data interpretations that account both for the implied coordinate transformation and the integrative properties of the system will be required if the contributions of sensory vestibular signals to constructing internal estimates of head motion and orientation are to be elucidated.


    APPENDIX
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DEVELOPMENT
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Analytical descriptions of cell and ocular responses predicted by the proposed model (Fig. 3) are summarized here. To describe the basic dynamic characteristics of the network, Laplace domain descriptions were obtained by treating the schematics in Fig. 3 as signal flow graphs in which circles are summing junctions that denote different cell populations and boxes are dynamic elements that represent either a sensor [canals, C(s) = Tcs/(Tcs + 1); otoliths, O(s) = 1/(Tos + 1)], the eye plant [P(s) = Kp/(Tps + 1)], or a neural filtering process [CLP(s) = 1/(Tcs + 1); F(s) = Kf/(Tfs + 1)]. Model parameters denote the strength of a given projection. Because projection weights from cells VO1 and VO2 onto cell VO4 vary as a function of the linear acceleration sensed either along the nasooccipital or the dorsoventral head axis, respectively, the system is nonlinear. Hence, the Laplace domain descriptions presented below can be used to approximate the dynamic characteristics of the network only for small head movements about a given operating point (i.e., a given average head orientation) when the system exhibits close to linear performance.

VO cell responses

General expressions.    Assuming small head deviations for a given static pitch angle {phi}, the responses of the VO cells in Fig. 3A to yaw and roll head velocity signals ({omega}z and {omega}x), sensed by the canals, and to interaural linear acceleration signals ({alpha}y), sensed by the otoliths, can be described by the general Laplace domain equation

(A1)
where s represents the complex Laplace variable and n is used to denote the cell number. The common gain term GNI1 can be written in terms of model parameters as

(A2)
The network time constant TNI1 is given by

(A3)
Model parameters were adjusted to set time constant TNI1 to a large value of 20 s in the upright orientation such that the network behaves as a leaky neural integrator. Notice that all response terms in Eq. A1 reflect this long time constant. VO cell responses are therefore predicted to embed the integrative properties of the network such that they may exhibit "velocity-storage" type characteristics (e.g., see Fig. 13).

Expression for gains GVOH, GVOV, and GVOL and time constants TVOH, TVOV, and TVOL as functions of the model parameters are summarized in Tables A1 and A2, respectively, for each cell type. Note that most gain and time constant terms depend on some combination of weights K1 and K2 that vary as a function of {alpha}x(t) or {alpha}z(t), respectively [i.e., K1(t) = K1o + K1a{alpha}x(t); K2(t) = K2oK2a{alpha}z(t)]. Thus these terms can be evaluated only for small deviations about a particular operating point where static approximations to weights K1 and K2 can be assumed. Specifically, the accelerations along the nasooccipital and dorsoventral head axes [{alpha}x(t) and {alpha}z(t)] can be approximated as static values {alpha}x {approx} –sin {phi} and {alpha}z {approx} cos {phi}, for a given head pitch angle {phi}. Parameters K1 and K2 can then be approximated as static weights K1 {approx} K1oK1asin {phi} and K2 {approx} K2oK2a cos {phi}. To calculate the bode plots in Figs. 6 through 8 estimates of K1 and K2 for a given average pitch angle were substituted into the appropriate gain and time constant expressions (above and in Tables A1 and A2) and used to evaluate Eq. A1 for each VO cell population.


View this table:
[in this window]
[in a new window]
 
TABLE A1. Gains in Eq. A1 expressed as a function of model parameters

 

View this table:
[in this window]
[in a new window]
 
TABLE A2. Time constants in Eq. A1 expressed as a function of model parameters

 
Choice of model parameters.    Model parameters in the network of Fig. 3A were chosen to satisfy several criteria for ideal tilt–translation discrimination at mid-high frequencies as well as to replicate particular properties of experimentally observed cell types. These conditions included: 1) network time constant TNI1 of about 20 s in the upright head orientation (see Eq. A2); 2) "canal-only" type responses on cells VO1 and VO2 that differ only in the axis of rotation for which each cell population is maximally stimulated (small values for K1b = K2b of 0.1, ph = p{nu}); 3) translation-selective responses on cell VO3 at mid-high frequencies (GVOH3 = GVOV3 = q1981{pi}Tc/180); 4) high-pass response to sensory otolith signals on cell VO3 (i.e., close to zero static tilt sensitivity, –1/TVOL3 {approx} 0); 5) sensitivities to acceleration during 0.5-Hz translation on cells VO3 and VO4 of 250 and 180 spikes·s–1·g–1, respectively (Dickman and Angelaki, 2002Go); 6) close to ideal estimate of head orientation relative to gravity (i.e., gy) on cell VO5 (K1o = K2o = 0; K1a = K2a; GVOL5 small, {approx}0 in units of cm/s2); 7) VO5 cell sensitivity to both high frequency dynamic and static tilts of approximately 0.7 spikes·s–1·deg–1, similar to that observed experimentally for cells sensitive mainly to tilt (Zhou et al. 2000Go).

Because of condition 6 above, the most general expression for cell VO5, given by Eq. A1, simplifies at high frequencies and for small head movements around a given static pitch orientation, (i.e., {alpha}x {approx} –sin {phi} and {alpha}z {approx} cos {phi} in units of g) to that presented in Eq 7 (see MODEL DEVELOPMENT)

(A4)
Note that this equation approximates gy only for small dynamic head reorientations relative to gravity as it predicts that cell VO5 encodes angular tilt position. For example, during head roll from upright, cell VO5 is predicted to encode a scaled estimate of angular roll position, {psi}, even though the actual interaural gravitational acceleration during roll is gy = –sin({psi}) (in units of g). This occurs because {alpha}x and {alpha}z are assumed static (i.e., constant during a given head movement) for ease of a linear systems exposition. In the actual model, however, the signal inputs to cell VO4 are multiplied by real dynamic signals {alpha}x(t) and {alpha}z(t). As a result, the simulated response of cell VO5 closely approximates an internal scaled estimate of gy that is accurate even for large angles during rotation. In addition, notice that Eq. A4 implies that canal signals make the dominant contribution to the response of cell VO5 at high frequencies. In contrast, at very low frequencies (e.g., <0.01 Hz) otolith signals dominate such that cell VO5 is also predicted to respond to static head tilts.

Earth-vertical–axis rotations.    During earth-vertical–axis rotations only the semicircular canals are stimulated. For an earth-vertical–axis rotation velocity {omega}e{nu} and pitch angle {phi}, the yaw and roll components of velocity can be expressed as {omega}z(s) = {omega}e{nu}(s)cos {phi} and {omega}x(s) = –{omega}e{nu}(s)sin {phi}, respectively. VO cell responses can be then described by the first two terms in Eq. A1 as

(A5)
To illustrate why cells VO3, VO4, and VO5 fail to respond in the model during earth-vertical–axis rotations we may examine the example case of cell VO4. For a given pitch head orientation, K1 = K1oK1asin {phi} and K2 = K2oK2acos {phi}. The cell's response can then be approximated at mid-high frequencies as

(A6)
When phK1a = p{nu}K2a and K1o = K2o = 0 (as for the chosen parameter set; see Fig. 3 legend) Eq. A6 illustrates that cell VO4 will not respond during earth-vertical–axis rotations. It can be shown that in this case the cell population only exhibits a sensitivity to canal signals along an axis orthogonal to gravity (i.e., earth-horizontal axis). However, more generally, if phK1a != p{nu}K2a on an individual VO4 cell, the neuron will indeed respond during earth-vertical–axis rotations but is predicted to respond with second-harmonic tuning (i.e., associated with the 1st two terms in Eq. A6; e.g., see Fig. 10C). Any head-orientation–invariant canal contribution to the cell's input (i.e., last two terms in Eq. A6 when K1o != 0 and/or K2o != 0) will contribute a first-harmonic component to the cell's tuning. In general, on any given cell, both first- and second-harmonic tuning components are predicted to be observed (e.g., Fig. 10D). The same considerations apply to model cells VO3 and VO5.

Ocular responses

Horizontal and torsional ocular responses (for small head movements about an average pitch setpoint) can be expressed in the Laplace domain as

(A7)
where

(A8)
is the oculomotor neural integrator time constant (i.e., long time constant set to 20 s; Angelaki et al., 2001aGo; Robinson, 1981Go) and

(A9)
is a gain term common to all response components. Gains GRVOR and GTVOR and time constant TTVOR can be expressed in terms of model parameters as

(A10)

(A11)

(A12)
where weights p (pdh, pd{nu}), e (eh, et), and qo (qoh, qot) differ slightly for horizontal versus torsional premotor eye movement networks (i.e., see Fig. 3 legend) to reflect differences in horizontal versus torsional RVOR gains and responses to interaural translation. Note that to achieve eye plant compensation in RVOR responses, Eq. A7 assumes that the neural filter, F(s) = Kf/(Tfs + 1), in Fig. 3B represents a scaled internal model of the eye plant, P(s) = Kp/(Tps + 1), where Tf = Tp (Galiana and Outerbridge 1984Go; Green and Galiana 1998Go) and that d2a = d1 (Angelaki et al. 2001aGo).

A more intuitive feel for the response predicted by Eq. A7 may be obtained by considering an approximation valid at mid-high frequencies [f >> 1/(2{pi}TNI1), f >> 1/(2{pi}TNI2), f >> 1/(2{pi}Tc)] that incorporates a simplified expression for cell VO3, valid for the chosen parameter set and for small head movements about a given static pitch orientation, {phi}

(A13)
where GVOH3 = phK1aK3bK3c/Tc and GVOV3 = p{nu}K2aK3bK3c/Tc. The first term in Eq. A13 predicts an eye velocity response proportional to angular head velocity. This corresponds to the RVOR component of ocular responses (either horizontal or torsional) driven by direct canal projections (associated with weights pdh and pd{nu}) onto premotor cells sensitive to contralaterally directed eye movements (i.e., EMC cells in Fig. 3B). The second term reflects canal and otolith contributions to the TVOR component of the ocular response. These are conveyed from cell VO3 onto the premotor horizontal and torsional networks by projections onto cells sensitive to ipsilaterally directed eye movements (i.e., projection with weight qo onto EMI cells in Fig. 3B).

Details of the differences in dynamic processing that occur in the RVOR versus TVOR premotor pathways of the proposed model have been described previously (Angelaki et al. 2001aGo; Green and Galiana 1998Go). Model parameters a, b, d1, d2, Kp, Kf, Tp, Tf, p = pdh, and e = eh in Fig. 3B are identical to those previously published (Angelaki et al. 2001aGo). Parameters qoh, qot, and et were adjusted here in conjunction with weight ql (i.e., otolith sensory input to cell VO3) to approximate experimentally observed horizontal versus torsional ocular responses to interaural translation (Angelaki 1998Go; see data in Fig. 4C). In addition, weight pdv was chosen to simulate a torsional RVOR gain during head roll of 0.76, slightly lower than the simulated horizontal RVOR gain of 0.88 (e.g., Angelaki and Hess 1996Go).


    GRANTS
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DEVELOPMENT
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This work was supported by National Institutes of Health Grants F32-DC-05271, EY-12814, and DC-0460 and National Aeronautics and Space Administration Grants NAG2-1493 and NNAO4CC77G.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DEVELOPMENT
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank B. Hess, D. Dickman, T. Hullar, E. Klier, A. Haque, and A. Shaikh for comments on the manuscript.


    FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Address for reprint requests and other correspondence: A. Green, Dept. of Anatomy and Neurobiology, Box 8108, Washington University School of Medicine, 660 South Euclid Ave., St. Louis, MO 63110 (E-mail: agreen{at}pcg.wustl.edu).


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DEVELOPMENT
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Angelaki DE. Two-dimensional coding of linear acceleration and the angular velocity sensitivity of the otolith system. Biol Cybern 67: 511–521, 1992a.[CrossRef][Web of Science][Medline]

Angelaki DE. Detection of rotating gravity signals. Biol Cybern 67: 523–533, 1992b.[CrossRef][Web of Science][Medline]

Angelaki DE. Three-dimensional organization of otolith-ocular reflexes in rhesus monkeys. III. Responses to translation. J Neurophysiol 80: 680–695, 1998.[Abstract/Free Full Text]

Angelaki DE, Bush GA, and Perachio AA. Two-dimensional spatiotemporal coding of linear acceleration in vestibular nuclei neurons. J Neurosci 13: 1403–1417, 1993.[Abstract]

Angelaki DE and Dickman JD. Spatiotemporal processing of linear acceleration: primary afferent and central vestibular neuron responses. J Neurophysiol 84: 2113–2132, 2000.[Abstract/Free Full Text]

Angelaki DE, Green AM, and Dickman JD. Differential sensorimotor processing of vestibulo-ocular signals during rotation and translation. J Neurosci 21: 3968–3985, 2001a.[Abstract/Free Full Text]

Angelaki DE and Hess BJM. Inertial representation of angular motion in the vestibular system of rhesus monkeys. I. Vestibuloocular reflex. J Neurophysiol 71: 1222–1249, 1994.[Abstract/Free Full Text]

Angelaki DE and Hess BJM. Inertial representation of angular motion in the vestibular system of rhesus monkeys. II. Otolith-controlled transformation that depends on an intact cerebellar nodulus. J Neurophysiol 73: 1729–1751, 1995.[Abstract/Free Full Text]

Angelaki DE and Hess BJM. Three-dimensional organization of otolith-ocular reflexes in rhesus monkeys. II. Inertial detection of angular velocity. J Neurophysiol 75: 2425–2440, 1996.[Abstract/Free Full Text]

Angelaki DE, Hess BJM, and Suzuki J-I. Differential processing of semicircular canal signals in the vestibulo-ocular reflex. J Neurosci 15: 7201–7216, 1995.[Abstract]

Angelaki DE, McHenry MQ, Dickman JD, Newlands SD, and Hess BJM. Computation of inertial motion: neural strategies to resolve ambiguous otolith information. J Neurosci 19: 316–327, 1999.[Abstract/Free Full Text]

Angelaki DE, Shaikh A, and Dickman JD. Responses of vestibular nuclei neurons during combinations of tilt and translation. Program No. 593.8 2003 Abstract Viewer/Itinerary Planner. Washington, DC: Society for Neuroscience, 2003. Online.

Angelaki DE, Wei M, and Merfeld DM. Vestibular discrimination of gravity and translational acceleration. Ann NY Acad Sci 942: 114–127, 2001b.[Web of Science][Medline]

Baker J, Goldberg J, Hermann G, and Peterson B. Optimal response planes and canal convergence in secondary neurons in vestibular nuclei of alert cats. Brain Res 294: 133–137, 1984a.[CrossRef][Web of Science][Medline]

Baker J, Goldberg J, Hermann G, and Peterson B. Spatial and temporal response properties of secondary neurons that receive convergent input in vestibular nuclei of alert cats. Brain Res 294: 138–143, 1984b.[CrossRef][Web of Science][Medline]

Bos JE and Bles W. Theoretical considerations on canal-otolith interaction and an observer model. Biol Cybern 86: 191–207, 2002.[CrossRef][Web of Science][Medline]

Brettler SC and Baker JF. Directional sensitivity of anterior, posterior, and horizontal canal vestibulo-ocular neurons in the cat. Exp Brain Res 140: 432–442, 2001.[CrossRef][Web of Science][Medline]

Buettner UW, Büttner U, and Henn V. Transfer characteristics of neurons in vestibular nuclei of the alert monkey. J Neurophysiol 41: 1614–1628, 1978.[Abstract/Free Full Text]

Busettini C, Miles FA, Schwarz U, and Carl JR. Human ocular responses to translation of the observer and of the scene: dependence on viewing distance. Exp Brain Res 100: 484–494, 1994.[Web of Science][Medline]

Bush GA, Perachio AA, and Angelaki DE. Encoding of head acceleration in vestibular neurons. I. Spatiotemporal response properties to linear acceleration. J Neurophysiol 69: 2039–2055, 1993.[Abstract/Free Full Text]

Cannon SC and Robinson DA. Loss of the neural integrator of the oculomotor system from brain stem lesions in monkey. J Neurophysiol 57: 1383–1409, 1987.[Abstract/Free Full Text]

Chen-Huang C and Peterson BW. Three-dimensional rotational and translational responses of vestibular neurons in squirrel monkeys. Program No. 168.11 2002 Abstract Viewer/Itinerary Planner. Washington, DC: Society for Neuroscience, 2002. Online.

Correia MJ and Money KE. The effect of blockage of all six semicircular canal ducts on nystagmus produced by dynamic linear acceleration in the cat. Acta Otolaryngol 69: 7–16, 1970.[Medline]

Crawford JD and Vilis T. Axes of eye rotation and Listing's law during rotations of the head. J Neurophysiol 65: 407–423, 1991.[Abstract/Free Full Text]

Dai MJ, Raphan T, and Cohen B. Spatial orientation of the vestibular system: dependence of optokinetic after-nystagmus on gravity. J Neurophysiol 66: 1422–1439, 1991.[Abstract/Free Full Text]

Dickman JD and Angelaki DE. Vestibular convergence patterns in vestibular nuclei neurons of alert primates. J Neurophysiol 88: 3518–3533, 2002.[Abstract/Free Full Text]

Dickman JD and Angelaki DE. Dynamics of vestibular neurons during rotational motion in alert rhesus monkeys. Exp Brain Res (Epub ahead of print), 2003

Einstein A. Uber das Relativitätsprinzip und die aus demselben gezogenen Folgerungen. Jahrb Radioakt 4: 411–462, 1908.

Endo K, Thomson DB, Wilson VJ, Yamaguchi T, and Yates BJ. Vertical vestibular input to and projections from the caudal parts of the vestibular nuclei of the decerebrate cat. J Neurophysiol 74: 428–436, 1995.[Abstract/Free Full Text]

Fernandez C and Goldberg JM. Physiology of peripheral neurons innervating semicircular canals of the squirrel monkey. II. Response to sinusoidal stimulation and dynamics of peripheral vestibular system. J Neurophysiol 34: 661–675, 1971.[Free Full Text]

Fernandez C and Goldberg JM. Physiology of peripheral neurons innervating otolith organs of the squirrel monkey. I. Response to static tilts and to long-duration centrifugal force. J Neurophysiol 39: 970–984, 1976a.[Abstract/Free Full Text]

Fernandez C and Goldberg JM. Physiology of peripheral neurons innervating otolith organs of the squirrel monkey. III. Response dynamics. J Neurophysiol 39: 996–1008, 1976b.[Abstract/Free Full Text]

Fukushima K, Perlmutter SI, Baker JF, and Peterson BW. Spatial properties of second-order vestibulo-ocular relay neurons in the alert cat. Exp Brain Res 81: 462–478, 1990.[Web of Science][Medline]

Galiana HL and Outerbridge JS. A bilateral model for central neural pathways in the vestibuloocular reflex. J Neurophysiol 51: 210–241, 1984.[Abstract/Free Full Text]

Glasauer S and Merfeld DM. Modelling three-dimensional vestibular responses during complex motion stimulation. In: Three-Dimensional Kinematics of Eye, Head and Limb Movements, edited by Fetter M, Haslwanter T, Misslisch H, and Tweed D. Amsterdam: Harwood Academic, 1997, p. 387–398.

Goldstein H. Classical Mechanics. Reading, MA: Addison–Wesley, 1980.

Graf W, Baker J, and Peterson BW. Sensorimotor transformation in the cat's vestibuloocular reflex system. I. Neuronal signals coding spatial coordination of compensatory eye movements. J Neurophysiol 70: 2425–2441, 1993.[Abstract/Free Full Text]

Green AM. Visual–Vestibular Interaction in a Bilateral Model of the Rotational and Translational Vestibulo-Ocular Reflexes: An Investigation of Viewing-Context-Dependent Reflex Performance (PhD thesis). Montreal, Canada: McGill University, 2000.

Green AM and Angelaki DE. A model to explore the relationship between tilt/translation discrimination and velocity storage. Vestibular Influences on Movement satellite of the 22nd Meeting of the Barany Society, Orcas Island, Sep. 22–26, 2002.

Green AM and Angelaki DE. Resolution of sensory ambiguities for gaze stabilization requires a second neural integrator. J Neurosci 23: 9265–9275, 2003.[Abstract/Free Full Text]

Green AM, Dickman JD, and Angelaki DE. Investigation of the relationship between tilt/translation discrimination and velocity storage. Program No. 565.8. 2002 Abstract Viewer/Itinerary Planner. Washington, DC: Society for Neuroscience, 2002. Online.

Green AM and Galiana HL. Hypothesis for shared central processing of canal and otolith signals. J Neurophysiol 80: 2222–2228, 1998.[Abstract/Free Full Text]

Guedry FE. Psychophysics of vestibular sensation. In: Handbook of Sensory Physiology: The Vestibular System, Pt 2, Psychophysics, Applied Aspects and General Interpretations, edited by Guedry FE Jr. Berlin: Springer-Verlag, 1974.

Hain TC. A model of the nystagmus induced by off vertical axis rotation. Biol Cybern 54: 337–350, 1986.[CrossRef][Web of Science][Medline]

Harris LR. Vestibular and optokinetic eye movements evoked in the cat by rotation about a tilted axis. Exp Brain Res 66: 522–532, 1987.[CrossRef][Web of Science][Medline]

Haslwanter T, Straumann D, Hess BJM, and Henn V. Static roll and pitch in the monkey: shift and rotation of Listing's plane. Vision Res 32: 1341–1348, 1992.[CrossRef][Web of Science][Medline]

Hess BJ and Angelaki DE. Inertial vestibular coding of motion: concepts and evidence. Curr Opin Neurobiol 7: 860–866, 1997.[CrossRef][Web of Science][Medline]

Iwamoto Y, Perlmutter SI, Baker JF, and Peterson BW. Spatial coordination by descending vestibular signals. II. Response properties of medial and lateral vestibulospinal tract neurons in alert and decerebrate cats. Exp Brain Res 108: 85–100, 1996.[Web of Science][Medline]

Kasper J, Schor RH, and Wilson VJ. Response of vestibular neurons to head rotations in vertical planes. I. Response to vestibular stimulation. J Neurophysiol 60: 1753–1764, 1988.[Abstract/Free Full Text]

Katz E, Vianney de Jong JM, Buettner-Ennever JA, and Cohen B. Effects of midline medullary lesions on velocity storage and the vestibulo-ocular reflex. Exp Brain Res 87: 505–520, 1991.[Web of Science][Medline]

Kushiro K, Dai M, Kunin M, Yakushin SB, Cohen B, and Raphan T. Compensatory and orienting eye movements induced by off-vertical axis rotation (OVAR) in monkeys. J Neurophysiol 88: 2445–2462, 2002.[Abstract/Free Full Text]

Mayne RA. A systems concept of the vestibular organs. In: Handbook of Sensory Physiology: The Vestibular System, edited by Kornhuber HH. New York: Springer-Verlag, 1974, p. 493–580.

Merfeld DM. Modeling the vestibulo-ocular reflex of the squirrel monkey during eccentric rotation and roll tilt. Exp Brain Res 106: 123–134, 1995.[Web of Science][Medline]

Merfeld DM, Park S, Gianna-Poulin C, Black F, and Wood S. Neural strategies to resolve ambiguous otolith cues: influence of rotational cues on human perception of tilt and translation. Program No. 268.18. 2003 Abstract Viewer/Itinerary Planner. Washington, DC: Society for Neuroscience, 2003. Online.

Merfeld DM and Young LR. The vestibulo-ocular reflex of the squirrel monkey during eccentric rotation and roll tilt. Exp Brain Res 106: 111–122, 1995.[Web of Science][Medline]

Merfeld DM, Young LR, Oman CM, and Shelhamer MJ. A multidimensional model of the effect of gravity on the spatial orientation of the monkey. J Vest Res 3: 141–161, 1993b.[Medline]

Merfeld DM, Young LR, Paige GD, and Tomko DL. Three dimensional eye movements of squirrel monkeys following postrotatory tilt. J Vest Res 3: 123–139, 1993a.[Medline]

Merfeld DM and Zupan LH. Neural processing of gravitoinertial cues in humans. III. Modeling tilt and translation responses. J Neurophysiol 87: 819–833, 2002.[Abstract/Free Full Text]

Merfeld DM, Zupan LH, and Peterka RJ. Humans use internal models to estimate gravity and linear acceleration. Nature 398: 615–618, 1999.[CrossRef][Medline]

Mergner T and Glasauer S. A simple model of vestibular canal-otolith signal fusion. Ann NY Acad Sci 871: 430–434, 1999.[CrossRef][Web of Science][Medline]

Musallam S and Tomlinson RD. Asymmetric integration recorded from vestbular-only cells in response to position transients. J Neurophysiol 88: 2104–2113, 2002.[Abstract/Free Full Text]

Musallam WS and Tomlinson RD. Model for the translational vestibuloocular reflex (VOR). J Neurophysiol 82: 2010–2014, 1999.[Abstract/Free Full Text]

Paige GD, Telford L, Seidman SH, and Barnes GR. Human vestibuloocular reflex and its interactions with vision and fixation distance during linear and angular head movement. J Neurophysiol 80: 2391–2404, 1998.[Abstract/Free Full Text]

Paige GD and Tomko DL. Eye movement responses to linear head motion in the squirrel monkey. I. Basic characteristics. J Neurophysiol 65: 1170–1182, 1991a.[Abstract/Free Full Text]

Paige GD and Tomko DL. Eye movement responses to linear head motion in the squirrel monkey. II. Visual-vestibular interactions and kinematic considerations. J Neurophysiol 65: 1183–1196, 1991b.[Abstract/Free Full Text]

Perlmutter SI, Iwamoto Y, Baker JF, and Peterson BW. Interdependence of spatial properties and projection patterns of medial vestibulospinal tract neurons in the cat. J Neurophysiol 79: 270–284, 1998.[Abstract/Free Full Text]

Perlmutter SI, Iwamoto Y, Baker JF, and Peterson BW. Spatial alignment of rotational and static tilt responses of vestibulospinal neurons in the cat. J Neurophysiol 82: 855–862, 1999.[Abstract/Free Full Text]

Raphan T, Dai M, and Cohen B. Spatial orientation of the vestibular system. Ann NY Acad Sci 656: 140–157, 1992.[Web of Science][Medline]

Raphan T, Matsuo V, and Cohen B. A velocity storage mechanism responsible for optokinetic nystagmus (OKN), optokinetic after-nystagmus (OKAN) and vestibular nystagmus. In: Control of Gaze by Brain Stem Neurons, edited by Baker R and Berthoz A. Amsterdam: Elsevier, 1977, p. 37–47.

Raphan T, Matsuo V, and Cohen B. Velocity storage in the velocity-ocular reflex arc (VOR). Exp Brain Res 35: 229–248, 1979.[Web of Science][Medline]

Raphan T and Schnabolk C. Modeling slow phase velocity generation during off-vertical axis rotation. Ann NY Acad Sci 545: 29–50, 1988.[Medline]

Raphan T and Sturm D. Modeling the spatiotemporal organization of velocity storage in the vestibuloocular reflex by optokinetic studies. J Neurophysiol 66: 1410–1421, 1991.[Abstract/Free Full Text]

Reisine H and Raphan T. Neural basis for eye velocity generation in the vestibular nuclei of alert monkeys during off-vertical axis rotation. Exp Brain Res 92: 209–226, 1992.[Web of Science][Medline]

Robinson DA. Linear addition of optokinetic and vestibular signals in the vestibular nucleus. Exp Brain Res 30: 447–450, 1977.[Web of Science][Medline]

Robinson DA. The use of control systems analysis in the neurophysiology of eye movements. Annu Rev Neurosci 4: 463–503, 1981.[CrossRef][Web of Science][Medline]

Schor RH, Miller AD, Timerick SBJ, and Tomko DL. Responses to head tilt in cat central vestibular neurons. II. Frequency dependence of neural response vectors. J Neurophysiol 53: 1444–1452, 1985.[Abstract/Free Full Text]

Schor RH, Miller AD, and Tomko DL. Responses to head tilt in cat central vestibular neurons. I. Direction of maximum sensitivity. J Neurophysiol 51: 136–146, 1984.[Abstract/Free Full Text]

Schwarz U, Busettini C, and Miles FA. Ocular responses to linear motion are inversely proportional to viewing distance. Science 245: 1394–1396, 1989.[Abstract/Free Full Text]

Schwarz U and Miles FA. Ocular responses to translation and their dependence on viewing distance. I. Motion of the observer. J Neurophysiol 66: 851–864, 1991.[Abstract/Free Full Text]

Seidman SH, Telford L, and Paige GD. Vertical, horizontal and torsional eye movement responses to head roll in the squirrel monkey. Exp Brain Res 104: 218–226, 1995.[Web of Science][Medline]

Siebold C, Anagnostou E, Glasauer S, Glonti L, Kleine JF, Tchelidze T, and Büttner U. Canal–otolith interaction in the fastigial nucleus of the alert monkey. Exp Brain Res 136: 169–178, 2001.[CrossRef][Web of Science][Medline]

Siebold C, Glonti L, Glasauer S, and Büttner U. Rostral fastigial nucleus activity in the alert monkey during three-dimensional passive head movements. J Neurophysiol 77: 1432–1446, 1997.[Abstract/Free Full Text]

Siebold C, Kleine JF, Glonti L, Tchelidze T, and Büttner U. Fastigial nucleus activity during different frequencies and orientations of vertical vestibular stimulation in the monkey. J Neurophysiol 82: 34–41, 1999.[Abstract/Free Full Text]

Skavenski AA and Robinson DA. Role of abducens neurons in vestibuloocular reflex. J Neurophysiol 36: 724–738, 1973.[Free Full Text]

Telford L, Seidman SH, and Paige GD. Dynamics of squirrel monkey linear vestibuloocular reflex and interactions with fixation distance. J Neurophysiol 78: 1775–1790, 1997.[Abstract/Free Full Text]

Tomlinson RD, McConville KMV, and Na E-Q. Behavior of cells without eye movement sensitivity in the vestibular nuclei during combined rotational and translational stimuli. J Vest Res 6: 145–158, 1996.[CrossRef][Web of Science][Medline]

Uchino Y, Ikegami H, Sasaki M, Endo K, Imagawa M, and Isu N. Monosynaptic and disynaptic connections in the utriculo-ocular reflex arc of the cat. J Neurophysiol 71: 950–958, 1994.[Abstract/Free Full Text]

Viéville T and Faugeras OD. Cooperation of the inertial and visual systems. In: Traditional and Non-Traditional Robotic Sensors, edited by Henderson TC. Berlin: Springer-Verlag, 1990.

Waespe W, Cohen B, and Raphan T. Dynamic modification of the vestibulo-ocular reflex by the nodulus and uvula. Science 228: 199–202, 1985.[Abstract/Free Full Text]

Waespe W and Henn V. Neuronal activity in the vestibular nuclei of the alert monkey during vestibular and optokinetic stimulation. Exp Brain Res 27: 523–538, 1977.[Web of Science][Medline]

Waespe W and Henn V. Conflicting visual-vestibular stimulation and vestibular nucleus activity in alert monkeys. Exp Brain Res 33: 203–211, 1978.[Web of Science][Medline]

Waespe W and Henn V. The velocity response of vestibular nucleus neurons during vestibular, visual, and combined angular acceleration. Exp Brain Res 37: 337–347, 1979.[Web of Science][Medline]

Wearne S, Raphan T, and Cohen B. Contribution of vestibular commissural pathways to spatial orientation of the angular vestibuloocular reflex. J Neurophysiol 78: 1193–1197, 1997.[Abstract/Free Full Text]

Wearne S, Raphan T, and Cohen B. Control of spatial orientation of the angular vestibuloocular reflex by the nodulus and uvula. J Neurophysiol 79: 2690–2715, 1998.[Abstract/Free Full Text]

Wearne S, Raphan T, and Cohen B. Effects of tilt of the gravito-inertial acceleration vector on the angular vestibuloocular reflex during centrifugation. J Neurophysiol 81: 2175–2190, 1999.[Abstract/Free Full Text]

Wei M and Angelaki DE. Cross-axis adaptation of the translational vestibulo-ocular reflex. Exp Brain Res 138: 304–312, 2001.[CrossRef][Web of Science][Medline]

Wilson VJ, Ikegami H, Schor RH, and Thomson DB. Tilt responses of neurons in the caudal descending nucleus of the decerebrate cat: influence of the caudal cerebellar vermis and of neck receptors. J Neurophysiol 75: 1242–1249, 1996.[Abstract/Free Full Text]

Wilson VJ, Yamagata Y, Yates BJ, Schor RH, and Nonaka S. Response of vestibular neurons to head rotations in vertical planes. III. Response of vestibulocollic neurons to vestibular and neck stimulation. J Neurophysiol 64: 1695–1703, 1990.[Abstract/Free Full Text]

Yakushin SB, Raphan T, and Cohen B. Spatial properties of otolith units recorded in the vestibular nuclei. Ann NY Acad Sci 871: 458–462, 1999.[CrossRef][Web of Science][Medline]

Young LR. Perceptions of the body in space: mechanisms. In: Handbook of Physiology. The Nervous System. Sensory Processes. Bethesda, MD: Am. Physiol. Soc., 1984, sect. 1, vol. III, pt. 2, p. 1023–1066.

Zhou W, Tang B-F, and King WM. Vestibular neurons encode linear translation and head tilt with respect to gravity. Soc Neurosci Abstr 26: 1491, 2000.

Zhou W, Tang B-F, and King WM. Responses of rostral fastigial neurons to linear acceleration in an alert monkey. Exp Brain Res 139: 111–115, 2001.[CrossRef][Web of Science][Medline]

Zupan LH, Merfeld DM, and Darlot C. Using sensory weighting to model the influence of canal, otolith and visual cues on spatial orientation and eye movements. Biol Cybern 86: 209–230, 2002.[CrossRef][Web of Science][Medline]




This article has been cited by other articles:


Home page
IOVSHome page
N. Vuillerme and R. Cuisinier
Sensory Supplementation through Tongue Electrotactile Stimulation to Preserve Head Stabilization in Space in the Absence of Vision
Invest. Ophthalmol. Vis. Sci., January 1, 2009; 50(1): 476 - 481.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
P. R. MacNeilage, N. Ganesan, and D. E. Angelaki
Computational Approaches to Spatial Orientation: From Transfer Functions to Dynamic Bayesian Inference
J Neurophysiol, December 1, 2008; 100(6): 2981 - 2996.
[Abstract] [Full Text] [PDF]


Home page
J. Neurosci.Home page
T. Yakusheva, P. M. Blazquez, and D. E. Angelaki
Frequency-Selective Coding of Translation and Tilt in Macaque Cerebellar Nodulus and Uvula
J. Neurosci., October 1, 2008; 28(40): 9997 - 10009.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
K. L. McArthur and J. D. Dickman
Canal and Otolith Contributions to Compensatory Tilt Responses in Pigeons
J Neurophysiol, September 1, 2008; 100(3): 1488 - 1497.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
R. F. Lewis, C. Haburcakova, and D. M. Merfeld
Roll Tilt Psychophysics in Rhesus Monkeys During Vestibular and Visual Stimulation
J Neurophysiol, July 1, 2008; 100(1): 140 - 153.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
H. Zhu, J. R. Jordan, S. P. G. Hardy, B. Fulcher, C. Childress, C. Varner, B. Windham, B. Jeffcoat, R. W. Rockhold, and W. Zhou
Linear acceleration-evoked cardiovascular responses in awake rats
J Appl Physiol, August 1, 2007; 103(2): 646 - 654.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
W. Zhou, Y. Xu, I. Simpson, and Y. Cai
Multiplicative Computation in the Vestibulo-Ocular Reflex (VOR)
J Neurophysiol, April 1, 2007; 97(4): 2780 - 2789.
[Abstract] [Full Text] [PDF]


Home page
J. Neurosci.Home page
A. M. Green, H. Meng, and D. E. Angelaki
A Reevaluation of the Inverse Dynamic Model for Eye Movements
J. Neurosci., February 7, 2007; 27(6): 1346 - 1355.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
W. Zhou, B. F. Tang, S. D. Newlands, and W. M. King
Responses of Monkey Vestibular-Only Neurons to Translation and Angular Rotation
J Neurophysiol, December 1, 2006; 96(6): 2915 - 2930.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
S. B. Yakushin, T. Raphan, and B. Cohen
Spatial Properties of Central Vestibular Neurons
J Neurophysiol, January 1, 2006; 95(1): 464 - 478.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
L. L. Chen and M. M. G. Walton
Head Movement Evoked By Electrical Stimulation in the Supplementary Eye Field of the Rhesus Monkey
J Neurophysiol, December 1, 2005; 94(6): 4502 - 4519.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
A. G. Shaikh, F. F. Ghasia, J. D. Dickman, and D. E. Angelaki
Properties of Cerebellar Fastigial Neurons During Translation, Rotation, and Eye Movements
J Neurophysiol, February 1, 2005; 93(2): 853 - 863.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
92/2/905    most recent
01234.2003v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (25)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Green, A. M.
Right arrow Articles by Angelaki, D. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Green, A. M.
Right arrow Articles by Angelaki, D. E.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2004 by the The American Physiological Society.