JN AJP: Gastrointestinal and Liver Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 93: 1633-1646, 2005. First published November 3, 2004; doi:10.1152/jn.00383.2004
0022-3077/05 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
93/3/1633    most recent
00383.2004v1
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 Google Scholar
Google Scholar
Right arrow Articles by Hess, B. J. M.
Right arrow Articles by Misslisch, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hess, B. J. M.
Right arrow Articles by Misslisch, H.

Canal–Otolith Interactions After Off-Vertical Axis Rotations. II. Spatiotemporal Properties of Roll and Pitch Postrotatory Vestibuloocular Reflexes

Bernhard J. M. Hess, Karin Jaggi-Schwarz and Hubert Misslisch

Department of Neurology, University of Zurich, Zurich, Switzerland

Submitted 14 April 2004; accepted in final form 31 October 2004


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We have examined the spatiotemporal characteristics of postrotatory eye velocity after roll and pitch off-vertical axis rotations (OVAR). Three rhesus monkeys were placed in one of 3 orientations on a 3-dimensional (3D) turntable: upright (90° roll or pitch OVAR), 45° nose-up (45° roll OVAR), and 45° left ear-down (45° pitch OVAR). Subjects were then rotated at ±60°/s around the naso-occipital or interaural axis and stopped after 10 turns, in one of 7 final head orientations, each separated by 30°. We found that postrotatory eye velocity showed horizontal–vertical components after roll OVAR and horizontal–torsional components after pitch OVAR that varied systematically as a function of final head orientation. The quantitative analysis suggests that, in contrast to the analogous yaw OVAR paradigm, a system of up to 3 real, gravity-dependent eigenvectors and eigenvalues determines the spatiotemporal characteristics of the residual eye velocities after roll and pitch OVAR. One of these eigenvectors closely aligned with gravity, whereas the other 2 determined the orientation of the earth horizontal plane. We propose that the spatial characteristics of eye velocity after roll and pitch OVAR follow the physical constraints of stationary orientation in a gravitational field and reflect the brain’s best estimate of head-in-space orientation within an internal representation of 3D space.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
When we move in 3-dimensional (3D) space, the brain uses afferent information from semicircular canals and otolith organs to generate and constantly update an internal representation of current head-in-space orientation with respect to gravity. Most head movements occur in planes that activate both the otolith and semicircular canal organs. Detailed knowledge about the central processing of these signals and their interactions with other self-motion signals is essential for a better understanding of the brain’s programming of motor behavior underlying spatial orientation, navigation, orienting, and postural control.

One fruitful way to explore the interaction of semicircular canal and otolith information involves testing unusual vestibular stimulus configurations. For example, tilting the head after prolonged constant-velocity rotation away from the former spinning axis creates a strong conflict between the canal-borne head rotation and the otolith-borne stationary head orientation signals. The resulting central vestibular signal configuration, i.e., strong rotational velocity cues in the presence of stationary graviceptive information, is physically compatible only with a rotation of the head about an axis parallel to earth vertical. Accordingly, the postrotatory vestibuloocular reflex (VOR) determined in such paradigms shows an eye velocity vector that reorients toward alignment with the direction of gravity (Angelaki and Hess 1994Go; Merfeld et al. 1993Go). Another slightly different conflict paradigm consists in evoking a sense of rotation by full-field optic flow stimuli in the head horizontal plane while being tilted relative to gravity. In this case, the velocity vector of the induced optokinetic response reorients toward earth vertical as soon as the dominant optokinetic visual drive is switched off by putting the subject in complete darkness. In contrast to the head-horizontal system, torsional or vertical optic flow stimuli create apparently little conflicting sense of rotation, presumably because they occur more commonly in everyday locomotion in the same plane as a relatively stationary gravity vector (Dai et al. 1991Go; Gizzi et al. 1994Go; Raphan and Sturm 1991Go). In line with this, little or no reorientation of torsional and vertical optokinetic afternystagmus (OKAN) has been observed in static tilt paradigms (Dai et al. 1991Go).

We have earlier proposed that after yaw head rotations or optokinetic stimulation, the spatial properties of the postrotatory responses can be understood in terms of an active tendency to rotate the head velocity cues toward alignment with the earth vertical. Accordingly, the horizontal postrotatory velocity tends to counter tilt whenever the head is tilted or comes to a halt in a tilted orientation relative to earth vertical (rotation hypothesis). To test this hypothesis, we have more recently investigated the spatial properties of horizontal vestibuloocular responses in a static vestibuloocular paradigm, which is comparable to the static optokinetic nystagmus (OKN)/OKAN paradigm. Subjects were rotated in the head yaw plane around off-vertical axes (OVAR) and successively stopped in various positions with respect to the direction of gravity (Jaggi-Schwarz et al. 2000Go). A detailed quantitative analysis of these 3D postrotatory yaw OVAR responses corroborated previous findings using a similar paradigm in humans (Harris and Barnes 1987Go), monkey (Raphan et al. 1992Go), and cat (Harris 1987Go) and showed mathematically that the underlying mechanism can indeed be described by an active rotation of the head velocity signals. This same rotation mechanism, however, cannot readily explain the spatial properties of roll or pitch VOR or OKN responses that have been described so far for the dynamic postrotatory tilt paradigm only. First of all, there is little or no reorientation of torsional or vertical optokinetic afterresponses (Dai et al. 1991Go) unless when using a dynamic postrotatory tilt paradigm (Hess and Angelaki 1995Go). Second, the spatial properties of postrotatory roll or pitch VOR or OKN that are observed in the dynamic tilt paradigm appear to result from facilitation of components along the earth vertical and disfacilitation of components that are oriented in the earth horizontal plane. That is, torsional or vertical components of postrotatory velocity perpendicular to gravity are suppressed, suggesting a cosine tuning of these responses with respect to gravity (projection or cosine tuning hypothesis; Angelaki and Hess 1994Go; Hess and Angelaki 1995Go).

Rather little is known about the spatial characteristics of postrotatory roll or pitch VOR as a function of static head orientations in space, i.e., without applying conflicting postrotatory tilts from supine or ear-down positions toward upright. In contrast to the dynamic tilt paradigm, a spatially stationary orientation of the head rotation axis ensures that during and immediately after stop of rotation the semicircular canal signals and the phasic otolith signals provide congruent information about the axis of ongoing or immediately preceding head rotation. There is thus no physical conflict between the centrally estimated orientation of axis of head rotation in space and the estimated actual head orientation at rotation stop. Because of the opposite directions of the angular velocity information from the semicircular canals and the centrally generated otolith-borne steady-state velocity, both represented in head-fixed coordinates, most of the vestibular head velocity signal will cancel at stop of rotation. Based on the residual eye velocity signals, however, this OVAR paradigm provides the unique opportunity to probe the postrotatory roll and pitch VOR in terms of 1) which coordinate system is used to represent residual torsional and vertical VOR velocity and 2) which mechanisms, a rotation or a projection if any, underlies the transformation of postrotatory roll and pitch VOR from head-fixed to space-fixed coordinates (Fig. 1).



View larger version (17K):
[in this window]
[in a new window]
 
FIG. 1. Projection vs. rotation model applied to postrotatory roll vestibuloocular reflex (VOR). A: projection. B: rotation of postrotatory torsional eye velocity onto/toward the direction of gravity in 45° backward tilted position. Top panels: horizontal ({Omega}hor), vertical ({Omega}ver), and torsional ({Omega}tor) postrotatory slow-phase velocity. Bottom panels: horizontal vs. torsional velocity (side view), and horizontal vs. vertical velocity (front view). Gray and black curves illustrate the effect of increasing coupling strengths (feedback) from 0.4 (gray) to 0.8 (black). G, direction of gravity; dashed straight/curved arrows, direction of projection/rotation of postrotatory velocity; {vartheta}, rotation angle about axis (R-axis). Note that the postrotatory response trajectory is determined in A by the estimated direction of gravity (specified by 2 angles, not shown), in B by the orientation of the rotation axis (R-axis) that is perpendicular to both gravity and the head rotation axis (specified by one angle, not shown), and a rotation angle ({vartheta} angle).

 
The answer to these questions cannot be deduced from dynamic tilt experiments as studies of optokinetic afterresponses show (Dai et al. 1991Go; Hess and Angelaki 1995Go). The postrotatory OVAR paradigm, in which otolith and semicircular canal cues interact in different static head orientations, prove to be useful to investigate the coordinate system, in which the central estimation of torsional and vertical head-in-space orientation is represented. It will be shown that, in contrast to the rotation model, the projection model leads to a high correlation between the geometric properties of postrotatory responses and the actual head-in-space orientation.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Experiments were performed on 3 rhesus monkeys (Macaca mulatta). Details of animal preparation, measurement, and calibration of 3D eye position were described in Jaggi-Schwarz et al. (2000)Go. All procedures were in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and approved by the Veterinary Office of the Canton of Zurich.

Experimental setup and protocol

In this study, animals were rotated around the naso-occipital (x-axis) or the interaural (y-axis; Fig. 2, A and B) in 2 different experimental configurations with a constant velocity of ±60°/s (positive/negative roll = clockwise/counterclockwise; positive/negative pitch = downward/upward) for 10 cycles (initial acceleration 180°/s2): in one configuration the axis of rotation was oriented 90° and in the other 45° off-vertical. Animals were stopped (decelerations 180°/s2) at one out of 7 end positions in space-fixed coordinates: ±0°, ±30°, ±60°, ±90° (positive, i.e., clockwise roll or downward pitch OVAR; negative, i.e., counterclockwise roll or upward pitch OVAR).



View larger version (41K):
[in this window]
[in a new window]
 
FIG. 2. Schematic of vestibular stimulation types and coordinate systems. A: 90° roll or pitch off-vertical axis rotation (OVAR): angle ({vartheta}) between rotation axis ({omega}) and gravity (G) is 90°. B: 45° roll or pitch OVAR: angle ({vartheta}) between rotation axis ({omega}) and gravity (G) is 45°. RP, rotation plane; RRP, resultant response plane. C: roll OVAR: front view of a subject after rotation in the roll plane to a final head orientation of –30°, expressed in the space-fixed coordinates yE and zE. Vector of gravity rotated by {varphi} = 30° in head-fixed coordinates (y, z). D: pitch OVAR: side view of a subject after rotation in the pitch plane to a final head orientation of 30°, expressed in space-fixed coordinates yE and zE. Vector of gravity rotated by {varphi} = –30° in head-fixed coordinates (x, z).

 
Data analysis

For data analyses and data presentation, we express the orientation of the animal relative to gravity as equivalent orientation of the gravity vector in a head-fixed coordinate system. Thus a final body orientation of, for instance, –30° roll or 30° pitch expressed relative to space (xEyEzE coordinates in Fig. 2, C and D) corresponded to a rotation of the gravity vector by +30° in roll and –30° in pitch when expressed in the head-fixed coordinate system (xyz coordinates in Fig. 2, C and D).

TEMPORAL DOMAIN ANALYSIS. Details of computation and desaccading of eye angular velocity ({Omega}) as well as the empirical analysis of the spatial orientation of postrotatory VOR velocity have been described previously (Jaggi-Schwarz et al. 2000Go). Briefly, we first computed postrotatory eye velocity components in the rotation plane, i.e., the plane orthogonal to the rotation axis (see gray discs and dashed lines labeled RP in Fig. 2). After roll OVAR, these orthogonal velocity components represent a trajectory, {Omega}orth = {Omega}ver + {Omega}hor, that lay in the yz plane (Fig. 2C), whereas following pitch OVAR, the analogous trajectory, {Omega}orth = {Omega}tor + {Omega}hor, lay in the xz plane (Fig. 2D). Second, we fitted (in a least-squares sense) the sum of 2 exponential functions to each of the 3 components of postrotatory eye velocity as well as to the orthogonal velocity trajectory ({Omega}orth). Third, we determined the peak of the exponential function fitted to the orthogonal response vector. This peak was used to delimit the final response portion of the postrotatory response, defined as the response from peak velocity to zero or to a nonzero offset value.

SPATIAL DOMAIN ANALYSIS. The remaining analyses were done in the spatial response domain, using the time domain fits of the postrotatory response. By plotting the horizontal, vertical, and torsional components against each other, we analyzed the projection of the 3D velocity trajectories onto the different head planes as follows. In a first empirical approach, we estimated the angle that the head-vertical z-axis and the orthogonal response trajectory subtended in the rotation plane by measuring the slope of a line fitted to the orthogonal postrotatory trajectory (i.e., the vertical–horizontal trajectory after roll OVAR and the torsional–horizontal trajectory after pitch OVAR; see Fig. 3, A and B). We call the plane formed by this line and the rotation axis the resultant response plane (RRP) (see gray rectangles labeled "RRP" in Fig. 2, A and B and gray rectangles in Fig. 3). The angle describes the slope of a straight line fitted to the postrotatory eye velocity trajectory in the resultant response plane relative to the rotation axis. It characterizes how far the final part of the eye velocity vector tilted (in the resultant plane) away from the rotation axis toward alignment with the gravity vector (Fig. 3, A and B). This angle was difficult to assess empirically because of the generally small magnitude and sluggishness of the peak of the postrotatory signals.



View larger version (22K):
[in this window]
[in a new window]
 
FIG. 3. 3D orientation of postrotatory eye velocity ({Omega}tor, {Omega}ver, {Omega}hor). After roll or pitch OVAR, the head stops in an orientation relative to gravity as indicated by the thick blue lines in A and B. A tangent line (dashed black line) to the trajectory of postrotatory eye velocity (red curve starting at location indicated by a dot) and the rotation axis (x-axis in A, y-axis in B) specify the resultant response plane (RRP). Angle specifies how far the resultant response plane is rotated relative to the xz plane (after roll OVAR in A) or relative to the yz plane (after pitch OVAR in B). Angle specifies the slope of the tangent line in the resultant response plane. Responses from animal JU in A after stop at 60°, = 55°, = 73°; from animal RO in B after stop at –240°, = –246°, = 68°.

 
The angular orientation of both, the tangent line and the resultant response plane, approximated the 3D spatial orientation of the final time course of postrotatory eye velocity (see angles and in Fig. 3).

In a parallel theoretical treatment, we fitted the postrotatory response with a system of differential equations that implemented a gravity-dependent 3D projection mechanism (cosine tuning) upstream to a leaky integrator within a feedback system. More specifically, we fitted the following function to the postrotatory VOR responses

(1)
This function mathematically describes the postrotatory VOR velocity as a sum of a signal coding head velocity ({omega}oto) at stop of rotation ("bias velocity") and a spatially transformed signal, called inertial velocity ({omega}inertial), based on afferent inputs from the 4 vertical semicircular canals ({omega}rlv: lumped vertical semicircular canal signal generated at stop of rotation) and otolith information about 3D head orientation in space (Fig. 4).



View larger version (10K):
[in this window]
[in a new window]
 
FIG. 4. Head velocity signals from the vertical semicircular canals ({omega}rlv) and otolith-borne velocity signals ({omega}oto) interact in a positive feedback loop with inertial signals ({omega}inertial) that are subject to a spatial (P) and temporal ({int}) filtering. The spatial filter depends on the current estimate of head-in-space orientation relative to gravity (g). K, diagonal matrix controlling the feedback strengths.

 
The rationale for this model is in short as follows: As far as the dynamics is concerned Eq. 1 represents the time domain solution of a 3D generalization of the feedback model proposed by Robinson (1977)Go to explain the horizontal postrotatory VOR. Translated from the frequency (Laplace) into the time domain this model reads as follows: + a(1 – k){omega} = rlv + a{omega}rlv, where "a" is the reciprocal time constant of a leaky integrator and "k" is a feedback constant. In words: head velocity and acceleration signals are driven in a central feedback network by the semicircular canal input that conveys information about head acceleration and velocity in a certain ratio (1:a). Note that the central signals represent the input one-to-one if there is no feedback (k = 0). To extend this model to 3 dimensions we replaced the reciprocal time constant a = 1/T and the feedback constant k with diagonal matrices A (all diagonal elements aii = a, off-diagonal elements aik = 0, for i != k) and K (3 diagonal elements: k1, k2, k3, off-diagonal elements zero). To implement the spatial response properties, we introduced a gravity-dependent spatial filter upstream to the integrator in the feedback loop, which modifies the velocity according to its orientation relative to gravity (Fig. 4). The effect of this filter can be described by a cosine tuning of the velocity signals ({omega}) in the feedback loop with respect to gravity (unit vector g, expressed in head coordinates), i.e., P{omega} = (g · {omega})g, where the dot" · " designs the scalar product between g and {omega}. Altogether, these changes modify the 1D Robinson equation as follows: + A(1 – KP){omega} = rlv + A{omega}rlv. The eigenvectors and eigenvalues of the resulting matrix F = A(1 – KP) on the left-hand side of this 1st-order differential equation fully characterize the spatiotemporal characteristics of this feedback system (for details see APPENDIX). Specifically, the eigenvectors determine the response directions of Eq. 1 as a function of the initial response vector {omega}oto and the centrally estimated direction of gravity g, whereas the eigenvalues characterize the dynamics of the response components along each of the 3 eigenvectors (see Table 1 in RESULTS). Depending on head orientation the response depends on 2 or 3 distinct eigenvectors (Fig. 5). Alternatively, we also used an analogous form of Eq. 1 to test the rotation hypothesis. For this, we replaced the cosine tuning in the feedback loop by a rotation operator, i.e., the matrix P in the feedback matrix F by a rotation matrix R [i.e., F = A(1 – KR)]. This system now reflected the spatial properties of the matrix R, which has only one real eigenvector, that is, the axis of rotation, in contrast to the matrix P, which in general has 3 distinct real eigenvectors (for details see APPENDIX).


View this table:
[in this window]
[in a new window]
 
TABLE 1. Estimated response parameters

 


View larger version (25K):
[in this window]
[in a new window]
 
FIG. 5. Feedback matrix F = A(1 – KP) exhibits 3 real eigenvectors v1, v2, and v3. Depending on head orientation relative to gravity, these vectors specify 3 (A) or only 2 distinct directions in space (B and C). Eigenvector v3 parallels the direction of gravity in the head-fixed coordinates x, y, and z (see monkey sketch), depending on the feedback strengths k1, k2, and k3, whereas eigenvectors v1 and v2 define a plane (A) or a line (B and C) perpendicular to gravity.

 
To further characterize the spatial properties as a function of head orientation we computed the volume of the parallelepiped that is spanned by 3 normalized eigenvectors (i.e., unit vectors 1 = v1/||v1||, 2 = v2/||v2||, and 3 = v3/||v3||), using the equation V = 1 {wedge} 2 · 3, where the wedge ("{wedge}") designs the cross-product. As shown in the APPENDIX, V depends on the feedback constants and gravity as follows: V = –k2g2g12 {sum}i kigi2 (ki ≥ 0). Accordingly, the eigenvectors form a right- (left-) handed triple depending on whether the head rolls right (g2 < 0) or left side down (g2 > 0). Technically, a value V = 1 or V = –1 means that the vectors are mutually orthogonal, whereas a V = 0 means that the 3 vectors lie in one plane (i.e., they are linearly dependent; Fig. 5).

SPATIAL ORIENTATION OF THE EIGENVECTOR v3 RELATIVE TO THE HEAD. The estimated spatial orientation of the head relative to gravity is described by the 3D orientation of the eigenvector v3 relative to the head. It was obtained by calculating the polar angles and of v3 as follows: First we calculated the orientation of the eigenvector v3 relative to the axis of rotation (ea: unit vector directed along the x-axis for roll OVAR or along the y-axis for pitch OVAR)

(2)
where 3 = v3/||v3|| is the normalized eigenvector v3 (3 = unit vector, ||v3|| = length of v3). The angular orientation of the eigenvector v3 in the rotation plane was obtained by (Fig. 3, A and B)

(3)

ESTIMATION OF POSTROTATORY EYE VELOCITY PARAMETERS AS A FUNCTION OF HEAD ORIENTATION. The solution (Eq. 1) of the feedback system (Fig. 4) provides an expression of postrotatory eye velocity as a function of a few relevant physiological parameters. These parameters constitute the angular orientation of eye velocity relative to gravity, the dominant response time constants, and the feedback gains in the horizontal, vertical, and torsional eye movement channel. To obtain a handle on these parameters, we represented the dynamics of the semicircular canals ({omega}rlv in Eq. 1) according to the model introduced by Steinhausen (1933)Go, with time constants fixed to 0.003 and 5 s (Wilson and Melvill Jones 1979Go). To estimate the otolith-driven head velocity signal, we empirically measured in each trial the amplitude of 3D eye velocity at stop of head rotation [{omega}oto = {omega}VOR(t = 0) in Eq. 1], which corresponded to the steady-state velocity that was maintained throughout the 10 rotations. Its main component was vertical for pitch OVAR and torsional for roll OVAR. Finally, we analytically solved the integral on the right-hand side of Eq. 1 by evaluating the matrix exponential etF, with F = A(1 – KP) being a unique function of final head orientation. According to the head-in-space orientation hypothesis, the spatial orientation of the postrotatory response depends ultimately on estimated head orientation relative to gravity, as described by the 2 polar angles and that characterize the orientation of the eigenvector v3 in head coordinates (see Fig. 3; see Eqs. A11 and A12 in the APPENDIX). Thus the nonlinear least-squares fit of Eq. 1 to the postrotatory response yielded an estimate of head orientation relative to gravity at stop of rotation (parameters and ), the time constant T (see Eq. A4), the 3 feedback gains k1, k2, and k3 (see Eqs. A4 and A5 in the APPENDIX), and the amplitude of 3D eye velocity, elicited by activating the vertical semicircular canals at stop of head rotation. A detailed description of Eq. 1 as a function of these parameters is presented in the APPENDIX.

To illustrate the fitted curves of postrotatory eye velocity in the spatial domain, we used normalized postrotatory curves, pn(t), which we computed as follows: If p(t) = [ptor(t), pver(t), phor(t)]T was the fitted postrotatory curve, the normalized curve was defined as pn(t) = p(t)/||p0||, where was the amplitude of the curve p(t) at stop of rotation (time t = 0: onset of postrotatory decay).

STATISTICS. Statistical analyses were performed with the Matlab statistics toolbox (The MathWorks, Natick, MA). Repeated-measures Friedman’s ANOVA was used to determine the significance of differences in the time and the feedback constants as a function of the tilt angle and rotation direction. To compute the variances of volume V and the time constant T3 we considered the propagation of errors of the 6 independent variables gi (i = 1 to 3) and ki (i = 1 to 3) (Bevington and Robinson 1992Go).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Figure 6A shows a representative example of perrotatory VOR velocity during the last 2 cycles of 90° roll OVAR at –60°/s (counterclockwise) and the subsequent postrotatory VOR velocity after the stop of rotation (vertical dashed line) in 60° left ear-down position. This corresponds to a tilt of gravity of +60° because we expressed gravity in head-fixed coordinates in a range of ±90° (Fig. 2A). Superimposed on the torsional, vertical, and horizontal components of eye velocity (gray data points in top, middle, and bottom panels), we plotted the exponential functions fitted to each component of eye velocity (black lines), starting at stop of rotation. The perrotatory velocity traces exhibit the well-known steady-state torsional bias velocity, in combination with the head-position–dependent modulation of torsional, vertical, and horizontal eye velocity (Angelaki and Hess 1996a, bGo; Hess and Dieringer 1990Go). After stop of rotation, the VOR response consisted of a negative torsional (counterclockwise; see METHODS) eye velocity, which gradually declined to a zero value.



View larger version (18K):
[in this window]
[in a new window]
 
FIG. 6. Three-dimensional VOR velocity during and after roll and pitch OVAR. A: when a subject is stopped in 60° left ear-down position after 90° roll OVAR, postrotatory eye velocity does not decline along the former (torsional) stimulation axis. Rather, nontorsional eye velocity components first build up before gradually declining to zero. Monkey JU. B: when a subject is stopped in 60° nose-down position after 90° pitch OVAR, torsional and horizontal VOR velocity components emerge and decay in addition to the vertical velocity induced by activation of the semicircular canals at the stop of rotation. Monkey SU. {Omega}tor, torsional; {Omega}ver, vertical; {Omega}hor, horizontal slow phase eye velocity (gray data); black curves, exponential fits to individual eye velocity components; Hroll/pitch, roll/pitch head position (potentiometer output reset each 360°).

 
Note that postrotatory eye velocity did not simply decrease along the former (roll) stimulation axis, as indicated by the additional negative vertical (upward) eye velocity, which built up and then slowly declined toward zero, and the positive progressively decreasing horizontal eye velocity. Also postrotatory vertical eye velocity did not decline to zero, an observation that was occasionally found in all animals, for both the 45° and 90° OVAR paradigms.

Figure 6B, using the same conventions as in Fig. 6A, illustrates an example of reorientation of postrotatory eye velocity after 90° pitch OVAR at 60°/s with a final head orientation of 60° nose down. This corresponds to a tilt of gravity of –60° in head coordinates (Fig. 2D). As indicated by the evolution of horizontal and torsional velocity components, postrotatory eye velocity did not decline along the pitch stimulation axis. The emergence of nonprincipal velocity components after stop of head rotation indicates that the postrotatory eye velocity vector reorients such that it aligns with the direction of gravity.

A better way to illustrate this finding, using the same data as shown in Fig. 6, is to plot the postrotatory eye velocity in the spatial domain (Fig. 7).



View larger version (26K):
[in this window]
[in a new window]
 
FIG. 7. Examples of spatial orientation of postrotatory eye velocity. AC: 90° roll OVAR paradigm. A (3D-view): increase and subsequent decline of horizontal ({Omega}hor), vertical ({Omega}ver), and torsional ({Omega}tor) eye velocity, represented by the normalized best-fit exponential curves to the same data as shown in Fig. 6A; the dot specifies the onset of postrotatory response. Note that eye velocity decreases along a line (dashed line) that is relatively close to the true direction of the gravity vector (blue line) in the rotation plane (= yz plane). B (front view): Evolution of vertical and horizontal (= orthogonal) eye velocity in the rotation plane. Angle specifies the orientation of the line of intersection of the resultant response plane and the rotation plane. C (top view): evolution of postrotatory velocity in the resultant plane, characterized by the angle of a tangent line to it (dashed line). DE: 90° pitch OVAR paradigm. D (3D-view): evolution of overshooting vertical, horizontal and torsional eye velocity, represented by the best-fit 3D exponential curves to the same data as shown in Fig. 6B; the dot specifies the onset of postrotatory response. Note that eye velocity decreases ultimately along a line (dashed line) that is close to the true direction of the gravity vector (blue line) in the rotation plane (= xz plane). E (side view): decline of horizontal and torsionial (= orthogonal) eye velocity in the rotation plane specifying the angle . F (top view): evolution of postrotatory response in the resultant plane specifying the angle .

 
For this, we used the normalized postrotatory curves, which we computed from the fitted postrotatory curves as described in METHODS. Inspection of the normalized exponential fits to postrotatory eye velocity after 90° roll OVAR in Fig. 7, AC (solid red lines) reveals that from the initial velocity (solid circle) both vertical and horizontal velocity components build up (Fig. 7B) and then decrease toward zero (the solid circle specifies the onset of the velocity vector’s trajectory). In other words, postrotatory eye velocity at first reorients and subsequently declines along a line that is in approximate alignment with gravity (compare in Fig. 7, A and B the dashed lines representing the line fit on eye velocity with the heavy blue line representing the true direction of gravity in head-fixed coordinates). Note that the emergence of nontorsional eye velocity components means that postrotatory eye velocity does not precisely decrease along the roll (x-) axis. We quantified the alignment of postrotatory eye velocity within the resultant response plane (plane labeled RRP in Fig. 7A) by computing the angle between the line fitted to the orthogonal response and the negative z-axis (Fig. 7B). In this example, = 55°, underestimating by 5° the true value of 60°. For a complete characterization of the eye velocity vector’s 3D orientation, we quantified the slope of its decay in the resultant response plane by determining the angle of the tangent line, fitted to the decreasing part of the postrotatory velocity (Fig. 7C). Here, = 73°, deviating by 17° from the orientation of gravity relative to the rotation axis of 90°.

In other cases the alignment of postrotatory velocity with the direction of gravity was achieved by changing merely the ratio between components rather than by a buildup of nonprincipal velocity components. The postrotatory eye velocity after 90° pitch OVAR illustrated in Fig. 6B provides an example of this mechanism. Plots of the respective exponential fits in the spatial domain (Fig. 7, DF, solid red lines) show that the alignment with the direction of gravity occurs—after an initial overshooting in vertical direction—by smoothly adjusting the torsional–horizontal (Fig. 7E) and the torsional–vertical postrotatory response components (Fig. 7F). In this example, the angle yielded a value of –59°, closely matching the true head orientation (–60°).

The slope of the decreasing part of postrotatory orthogonal eye velocity, here the vector sum of horizontal and torsional eye velocity (i.e., {Omega}orth = {Omega}tor + {Omega}hor) was 88°, underestimating the true inclination of the head relative to gravity in the resultant response plane by 2° (see Fig. 7F).

Spatial orientation of postrotatory eye velocity after roll and pitch OVAR

The reorientation of postrotatory eye velocity seen in the examples of Figs. 6 and 7 was observed in all final body orientations independent of the stimulus configuration, i.e., after 90° roll or pitch OVAR (Fig. 8) as well as after 45° roll or pitch OVAR (Fig. 9). More precisely, the diagrams in the top and bottom panels of Fig. 8, A and B and Fig. 9, A and B plot the polar angles and of eye velocity (solid circle) as a function of the head orientation relative to gravity. The open squares in the bottom panels represent empirically evaluated values, averaged over both rotational directions (clockwise/counterclockwise roll, forward/backward pitch OVAR) and all 3 subjects (means ± 1SD). The solid circles represent the simultaneously determined and angles based on least-squares fits of Eq. 1 to postrotatory responses.



View larger version (25K):
[in this window]
[in a new window]
 
FIG. 8. Spatial orientation of postrotatory eye velocity after 90° roll (A) and pitch (B) OVAR. Top panels: solid black circles represent mean values (±SD) of the slope of the tangent line to postrotatory velocity in the resultant response plane. Note that after roll OVAR is very close to the expected 90° independent of head orientation. After pitch OVAR exceeds on average the expected 90° in 30° and 60° nose-up position. Bottom panels: close relationship between average orientations of eye velocity (solid circles and open squares) and gravity (dashed lines) in the rotation plane represented by the angle . Data points (open squares) are averages over both OVAR directions and all 3 subjects (±SD) for all final body orientations obtained from the empirical analysis. Solid black circles are obtained from model-based least-squares fits to postrotatory responses yielding simultaneous estimates (bottom panels) and (top panels) (±SD).

 


View larger version (26K):
[in this window]
[in a new window]
 
FIG. 9. Spatial orientation of postrotatory eye velocity after 45° roll (A) and pitch (B) OVAR. Top panels: solid black circles represent mean values (±SD) of the slope of the tangent line to postrotatory velocity in the resultant response plane. Note that after roll OVAR undershoots the expected 135° on average by about 25°. After pitch OVAR it is on average close to the expected value of 45°. Bottom panels: close relationship between average orientations of eye velocity (open squares and solid circles) and gravity (dashed line) in the rotation plane represented by the angle . Data points (open squares) are averages over both OVAR directions and all 3 subjects (±SD) for all final body orientations obtained from the empirical analysis. Solid black circles are obtained from model-based least-squares fits to postrotatory responses yielding simultaneous estimates of (bottom panels) and (top panels) (±SD).

 
Because there was no significant difference between clockwise and counterclockwise roll with the exception of in the 45° roll OVAR paradigm (90° roll OVAR P{varphi} > 0.5; P{vartheta} > 0.7; 45° roll OVAR P{varphi} > 0.5; P{vartheta} < 0.01) or forward and backward pitch rotations (pitch OVAR paradigms: P{varphi} > 0.7; P{vartheta} > 0.2), and angles were averaged over the 2 directions and all 3 subjects (means ± 1SD). Thereby we did not consider the differences in between clockwise and counterclockwise roll in the 45° roll OVAR paradigm, which was attributed to a response asymmetry in right ear-down position (45° roll OVAR: P{vartheta} > 0.02 after exclusion of right ear-down tilts). Final body orientations of 0°, 90°, and –90° correspond to upright, left ear-down, and right ear-down (roll) or upright, supine, and prone (pitch). The predictions of a perfectly reorienting VOR velocity in the rotation plane are indicated by the dashed diagonal in Figs. 8 and 9. For instance, the estimated should ideally be 60° when animals are stopped in a final position such that the gravity vector’s orientation in the pitch or roll plane is 60°. This statement applies to both 90° and 45° roll or pitch OVAR because the {varphi} angle was always measured in the rotation plane. In general, we found a close correlation between the estimated values and the direction of gravity in the rotation plane, as supported by the high r2 values from linear regression analyses: 90° or 45° roll OVAR r2 = 0.96 or 0.94, 90° or 45° pitch OVAR r2 = 0.98 or 1.0. Linear regression of the angle yielded a slope of 0.75 ± 0.17 and an offset = –0.8 ± 5.4° for 90° roll OVAR compared with 0.86 ± 0.13 and 1.5 ± 8.4° for 90° pitch OVAR. Thus the direction of postrotatory responses after both roll or pitch OVAR rather tightly followed the direction of gravity in the tilt plane. At larger tilt angles there was a tendency to undershoot the veridical head tilt, in particular in the 90° roll OVAR paradigm.

For a complete description of the spatial orientation we also quantified how much gravity determined the trajectory of the postrotatory eye velocity response in a plane orthogonal to the rotation plan. In contrast to the angle that described the orientation of eye velocity in the pitch or roll plane, the angle determined the orientation in the orthogonal plane set up by the estimated direction of eye velocity in the rotation plane (i.e., the line fitted to the decreasing part of orthogonal eye velocity) and the rotation axis (see resultant response plane, Figs. 2, 3, and 7). More specifically, we defined as the angle subtended by a tangent line to the final portion of the response trajectory and the rotation axis. The average angle after 90° roll and pitch OVAR was = 79.7 ± 7.5° and 93.0 ± 24.1°, respectively (in both paradigms the veridical angle is 90°; Figs. 2A and 8, A and B). Whereas the {vartheta} angle in the roll paradigm was estimated close to veridical, at least for tilt angles between +60° and –60° (Fig. 8A) the estimated values in the 90° pitch paradigm deviated from the true values by an average of 36.9° when the animals were stopped in 30° and 60° nose-up position (Fig. 8B). Note also that the estimates showed much larger SDs in the pitch than in the roll OVAR paradigm. After 45° roll OVAR, was 110.4 ± 3.2° (the veridical angle {vartheta} = 135°; Figs. 2B and 9A) and after 45° pitch OVAR was 48.4 ± 11.6° (the veridical angle {vartheta} = 45°; Figs. 2B and 9B).

Properties of feedback constants, roll and pitch eigenvectors, and time constants

The feedback strengths in the different eye movement channels (see diagonal matrix K in Fig. 4) ultimately determine, together with the central estimate of head-in-space orientation, g, as a function of and (see Fig. 3), the spatial orientation as well as the associated dynamics of the postrotatory response (see Eqs. A3 and A4 in the APPENDIX). In fact, the response parallels the estimated head-in-space-orientation only if the ratios of the feedback strengths in the relevant (i.e., nonzero) feedback channels are approximately equal. In the 90° roll OVAR paradigm, the feedback strength of the horizontal (k3) and vertical channel (k2), which control the spatiotemporal characteristics of the orthogonal postrotatory response in the roll plane, did not differ significantly from each other (Friedman’s ANOVA, P > 0.5). We therefore pooled these values over the 3 animals and the 2 directions of rotations and compared these averages with the average feedback in the torsional channel k1 (Fig. 10A). The torsional feedback was significantly smaller than the horizontal or the mean of the horizontal and vertical feedback (Friedman’s ANOVA, P < 0.01 for both comparisons). This means that most of the response ultimately decayed in the roll plane as illustrated by the average angles in Fig. 8A. Similarly, one would expect the feedback strength in the vertical channel to be relatively small in the 90° pitch paradigm, such that responses ultimately decayed in the pitch plane. However, this was not observed, consistent with the large variances and the deviations of the average angles from the pitch plane in Fig. 8B. Because there were no significant differences in strength of the three feedback channels in the remaining paradigms (Friedman’s ANOVA, P > 0.1), we pooled the values in each of the 90° pitch, 45° roll, and 45° pitch OVAR paradigms over the 3 animals and the 2 directions (Fig. 10, BD).



View larger version (30K):
[in this window]
[in a new window]
 
FIG. 10. Properties of feedback constants as a function of head orientation after 90° roll (A), 90° pitch (B), 45° roll (C), and 45° pitch (D) OVAR. In A, solid black circles represent the feedback averaged over the horizontal and vertical eye movement channels and gray circles represent the average torsional feedback; in B to D solid black circles (±SD) represent the feedback averaged over all 3 eye movement channels.

 
Because the eigenvectors of Eq. 1 in general change orientation as a function of head-in-space orientation, it is of interest to see whether the orientation of the 3 eigenvectors, which ultimately determine the direction of postrotatory eye velocity, was independent of how final head orientations were reached in different paradigms. As a consequence of Eq. 1 only the postrotatory responses after 45° roll/pitch OVAR were characterized by 3 linearly independent eigenvectors (see Fig. 5), as indicated by the nonzero volume V of the parallelepiped that these vectors encompassed throughout the tilt range (Fig. 11). Although linearly independent, the eigenvectors were not mutually orthogonal, as demonstrated by the fact that the magnitude of the volume V was less than unity. Furthermore, the sign of V systematically differed in 45° roll responses: In right ear-down/nose-up positions, the orientation of the eigenvectors resulted in a positive V (Fig. 11A), whereas V was always negative in left ear-down/nose-up positions. Technically, this observation means that the eigenvectors switched from a right-handed to a left-handed triple when the y-component of gravity changed from positive to negative (Fig. 11A). In the 45° pitch OVAR paradigm, the orientation of the eigenvectors did not switch, as expected (Fig. 11B). Comparison of the orientations of the eigenvectors across the 2 paradigms shows that they were similar in the same final head orientations (led/nu positions in Fig. 11, A and B).



View larger version (22K):
[in this window]
[in a new window]
 
FIG. 11. Properties of eigenvectors as a function of head orientation after 45° roll (A) and 45° pitch OVAR (B). Three eigenvectors encompass a positive volume V, indicating a right-handed system in position with the right ear down (A), and a negative volume V in positions with the left ear down (A, B), indicating a left-handed system. Solid circles represent means (±SD).

 
Associated with each eigenvector is an eigenvalue (= reciprocal time constant) that characterizes the dynamics of the postrotatory response component along this eigenvector. The time constants associated with the first and second eigenvector were equal (see Eq.A5 in the APPENDIX). They were similar for each stimulus configuration and did not differ significantly for clockwise versus counterclockwise and forward versus backward rotation (T1: Friedman’s ANOVA, P > 0.25 in 90° OVAR, P > 0.05 in 45° OVAR paradigms). They assumed average values between 4 and 6 s (see Table 1). In contrast, the time constants associated with response components along the third eigenvector were on average about 5 times longer and also did not differ significantly for clockwise versus counterclockwise and forward versus backward rotation (T3: Friedman’s ANOVA, P > 0.25 in 90° OVAR, P > 0.05 in 45° OVAR paradigms).

In addition to the estimated {varphi} and {vartheta}, the feedback constants, and the time constants, the least-squares fit of Eq. 1 yielded also estimates of the amplitude (v0) of the semicircular canal–induced response generated at stop of rotation. The averages (±1SD) across the 2 rotation directions and all tilt angles of all 3 animals are summarized for the 90° and 45° roll and pitch OVAR paradigms in Table 1.

We also fitted equation 1 to all postrotatory responses after roll and pitch OVAR in all 3 animals after replacing the projection P in the feedback matrix F with a rotation R (for details see APPENDIX). Although the fitted curves were often indistinguishable from those obtained with the projection model, the corresponding angular head orientation parameters did not correlate well with the final head orientation at stop of rotation. Table 2 summarizes the variances accounted for (VAF) of the estimated angles for both the projection and rotation model. It also shows the average estimates of the angles for the 2 models.


View this table:
[in this window]
[in a new window]
 
TABLE 2. Comparison of estimated {varphi} and {vartheta} obtained with the projection and rotation model

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Roll and pitch movements of the head from upright activate both the vertical semicircular canals as well as the otolith organs, which code information about head acceleration, velocity, and position. The brain uses these vestibular and other proprioceptive signals to generate an internal space-centered reference system, which enables orientation of gaze in space independent of self-motion. Our analysis of postrotatory responses after roll and pitch OVAR suggests that the spatial properties of residual 3D eye velocity reflect the brain’s best estimate of head-in-space orientation. The underlying mechanisms can conceptually be thought of as a projection of the vector of residual postrotatory eye velocity onto a central representation of gravity within a feedback system. Alternatively, postrotatory eye velocity could actively rotate toward alignment with an internal estimate of gravity, as described earlier for responses after yaw OVAR (Jaggi-Schwarz et al. 2000Go). For comparision, simulations of these 2 fundamentally different mechanisms have been illustrated in Fig. 1 for different feedback strengths. In contrast to the yaw OVAR paradigm, we found that after roll and pitch OVAR not only the centrally estimated direction of gravity but also 2 orthogonal directions, which depended on the particular head orientation relative to gravity, determined the geometric properties of postrotatory velocity.

Residual 3D postrotatory eye velocity reflect 3D head-in-space orientation

In our paradigms the subjects were stopped after several cycles of rotations in complete darkness, well after the rotation cues from the semicircular canals have subsided. In this situation self-orientation in space faces 2 problems, summarized here as the head-in-space orientation problem: First, to estimate the final head and body orientation relative to gravity as quickly and accurately as possible and second, based on this estimation, to further estimate whether there is any head translation by comparing the total afferent otolith input with the centrally estimated head orientation. This second step is necessary because the otolith afferents carry information about the sum of gravitational and inertial head accelerations.

After a sudden stop head velocity signals, generated by activation of the vertical semicircular canals and otolith signals, reflecting head-in-space velocity and instantaneous orientation relative to gravity (Angelaki 1992a, bGo; Hain 1986Go; Schnabolk and Raphan 1992Go), are thought to interact within feedback loops with central signals. Immediately after the stop most of the persisting angular eye velocity ({omega}oto in Fig. 4) is cancelled by vertical canal signals ({omega}rlv in Fig. 4), a process that is virtually perfect under visual conditions.

To quantitatively analyze the characteristics of the residual velocity signals when rotation stops in the dark, we adopted the following assumptions: First, we assume that the inertial velocity signals conform to the physical constraints that exist between the rate of change in head orientation, the instantaneous head orientation, and head angular velocity. These constraints can be quantitatively expressed by the vector equation daoto/dt = {omega}head {wedge} aoto + d{alpha}oto/dt u (see Goldstein 1980Go), which captures the spatiotemporal relation between the time derivative of otolith signals (denoted by daoto/dt) that encodes the rate of change in head orientation (i.e., rotation and/or translation), and a rotation of the head relative to gravity (denoted by whead {wedge} aoto), and/or a translation of the head. Note that a pure translation results in a change in acceleration magnitude but not in direction (thus: daoto/dt = d{alpha}oto/dt u, with aoto = {alpha}otou, where u = unit direction of translation). If the brain’s best estimate is that head-in-space orientation (denoted in Fig. 4 by the estimated orientation of gravity, g, in head coordinates) is stationary at stop of rotation (i.e., g {approx} aoto), it follows from the physical stationary condition for otolith signals that the residual head angular velocity must parallel the estimated head-in-space orientation (i.e., parallel vectors: {omega}head ~ g). An implicit assumption of our approach is that the centrally estimated head orientation relative to gravity does not much change over the time course of postrotatory nystagmus. Recent human psychophysical studies appear to support this notion (Jaggi-Schwarz and Hess 2003Go; Jaggi-Schwarz et al. 2003Go), although some earlier studies found considerable variability in estimates of the subjective visual vertical over time, particularly in roll tilt positions beyond 90° (Udo de Haes and Schöne 1970Go; von Holst and Grisebach 1951Go). Second, we assume that the mechanisms, which underlie the transformation of residual head velocity signals toward alignment with the estimated head-in-space orientation, are based on cosine tuning, also referred to as projection, because only the vector component along the estimated direction of gravity survives.

Finally, we account for the well-documented central enhancement of the time constants of the postrotatory responses (velocity storage: Cohen et al. 1977Go; Raphan et al. 1977Go, 1979Go) by a leaky integration within a positive feedback loop (Robinson 1977Go). The stationary hypothesis (vector parallelism, {omega}head ~ g, if the brain estimates dg/dt = 0) finds ample support in respective conflict paradigms (Angelaki and Hess 1994Go; Dai et al. 1991Go; Gizzi et al. 1994Go; Harris 1987Go; Harris and Barnes 1987Go; Hess and Angelaki 1995Go; Jaggi-Schwarz et al. 2000Go; Merfeld et al. 1993Go, 1999Go, 2001Go; Raphan and Sturm 1991Go; Raphan et al. 1992Go), where it predicts that postrotatory head velocity signals in the VOR ultimately must align with the brain’s best estimate of direction of gravity. This is so because it represents the only physically possible configuration of a postrotatory head velocity signal and a stationary head-in-space-orientation. That the head must be stationary in these conflict experiments is the brain’s best estimate based on otolith and proprioceptive somatosensory cues. Mechanisms underlying the spatial transformations of head angular velocity signals toward alignment with gravity have been suggested in a number of studies (Angelaki and Hess 1994Go, 1995Go; Angelaki et al. 1995Go; Hess and Angelaki 1995Go; Merfeld 1995Go; Merfeld and Zupan 2002Go; Mergner and Glasauer 1999Go; Raphan and Sturm 1991Go; Zupan et al. 2002Go). Specifically, Angelaki and Hess (1994)Go showed evidence that the spatial transformation of postrotatrory responses show different characteristics in the horizontal versus the torsional and vertical eye movement pathways: Whereas horizontal postrotatory responses after dynamic head tilts in pitch or roll showed characteristics compatible with a rotation, torsional and vertical responses after tilts in the yaw plane showed characteristics of a projection. Recently we have extended this finding by quantitatively showing that the reorientation of angular velocity signals in the yaw OVAR paradigm is compatible with the geometric properties of a rotation (Jaggi-Schwarz et al. 2000Go). In the following we discuss the experimental evidence for a projection (cosine tuning) mechanism in the torsional and vertical VOR pathways and compare the utility of rotation and projection mechanisms.

Characteristics of head-in-space orientation estimates after roll or pitch OVAR

Our analysis of the orientation properties of postrotatory responses after roll and pitch OVAR results in reasonably good predictions of the actual head orientation of the monkey in space. In contrast to the yaw OVAR paradigm (Jaggi-Schwarz et al. 2000Go), head-in-space-orientation could not be predicted from the postrotatory responses after roll and pitch OVAR by implementing a rotation mechanism (rotation matrix R instead of projection P in Fig. 4). As earlier outlined, the 2 transformations have different spatial properties, which entail different dependencies of the postrotatory response profiles on head-in-space orientation. In the projection model, the feedback loop enforces spatial alignment of residual eye velocity with the direction of gravity as the feedback constants tend toward equal strengths in the relevant eye movement channels. In the rotation model, the feedback works differently in that it mainly controls the direction but not as efficiently the magnitude of rotation. The magnitude of rotation depends in a highly nonlinear fashion on the parameters of head orientation in space. Simulations of this model show that the feedback needs to be limited in strength to avoid overshooting responses (see example in Fig. 1B). We found that least-squares fits of the rotation model to postrotatory responses yielded angles and that poorly correlated with the final head-in-space orientation in both the roll and pitch OVAR paradigm (see Table 2). This was true even though these fits were often not distinguishable from those of the projection model.

In all OVAR paradigms monkeys estimated head-in-space orientation significantly better in the tilt plane than in a plane perpendicular to it (called "resultant response plane"; see Fig. 2). This is true for the earlier studied yaw OVAR paradigm (Jaggi-Schwarz et al. 2000Go) as well as for the roll and pitch paradigms. Specifically, we observe negligible differences between the average estimated and true head orientation in the 90° roll or pitch OVAR paradigm if the tilt range is restricted to ±60° around upright or to upright and forward tilted positions, respectively. Nose-up orientations are often under- or overestimated with respect to the resultant plane ({vartheta} angle) or even systematically underestimated, such as after the animal has been rotated about a 45° nose-up tilted axis (Fig. 9A).

To further characterize the otolith-dependent orientation signals, we quantified the mutual orientation of the eigenvectors as a function of head orientation. This analysis demonstrates that, independent of the paradigm, the spatial trajectory of postrotatory responses are attributed to similar orientations of the eigenvectors at the same final head orientation (compare led/nu position in the 45° roll and pitch paradigms, Fig. 11). A comparison of the estimated orientation in the resultant response plane ({vartheta} angle), however, shows that this orientation parameter is not consistently estimated in the same final head positions in the roll and pitch OVAR paradigm. A possible explanation for this discrepancy is that estimation of the {vartheta} parameter depends not only on the final head orientation but also on the previous rotation.

Tilt plane and head position dependent deficits in 3D estimates of head-in-space orientation

What is the origin of the observed head orientation–dependent estimation errors? One possibility is that it reflects a systematic limitation of the projection or cosine-tuning hypothesis. Raphan and Sturm (1991)Go assumed in their velocity storage model that only the yaw eigenvector reorients toward earth vertical, whereas the roll and pitch eigenvectors remain head fixed. Although these assumptions are perfectly sufficient to explain the reorientation of postrotatory eye velocity in yaw OVAR paradigms they are incompatible with a reorientation of postrotatory responses after pitch and roll OVAR. The projection hypothesis, on the other hand, has been successfully applied in explaining the spatial orientation of torsional and vertical postrotatory eye velocity in dynamic tilt paradigms (Angelaki and Hess 1994Go; Angelaki et al. 1995Go). Furthermore it has been used explicitly in some observer-type models to explain the interaction of eye velocity and gravity (Droulez and Darlot 1990Go; Zupan et al. 2002Go) but not in others (Merfeld 1995Go; Merfeld and Zupan 2002Go; Merfeld et al. 1993Go). This is not surprising because, to our knowledge, none of these models has explicitly addressed the head-in-space orientation problem after rotations in the roll and pitch plane.

Other reasons that could limit the accuracy of estimating head-in-space orientation in nose-up positions could be rooted in the way of how backward directed otolithic shear forces are encoded and processed. Anatomically it is known that the anterior portion of the utricles is bent upward in primates (Lindenman 1973Go; Spoendlin 1966Go), possibly to improve the sensitivity for head orientation cues in a wider range of forward tilted positions. From clinical examinations it is known that static imbalance, such as in one-leg stance, increases when the head is extended, suggesting that the otolith signal processing is less effective in this posture (Brandt 1999Go).

Consequences of misjudging head-in-space orientation

Erroneous interpretation of gravitoinertial cues are the source of a number of illusions (for a review see Bos and Bles 2002Go; Young 2002Go). In a recent study in humans, Merfeld and colleagues (1999)Go reported the perception of illusionary tilt in a postrotatory tilt paradigm involving strong conflicting rotational cues. Along these lines, the head-in-space orientation errors in the nose-up positions may suggest that the monkeys tended to misjudge their orientation in a direction perpendicular to the tilt plane and might have experienced an illusionary perception of forward motion at stop of rotation. To our knowledge, there are no quantitative studies on the perception of self-orientation after roll tilt in nose-up orientations.

Differences in the spatial transformation properties of yaw versus pitch and roll postrotatory OVAR

Why are there fundamental differences in the spatial representation of yaw versus roll or pitch head velocity signals? One reason could be that these differences reflect the geometric constraints on vestibular balance control while performing gaze orientating head and body movements during locomotion and foraging. To maintain balance in a gravitational environment it is advantageous to keep the angular head and body velocity approximately aligned with gravity to comply as much as possible to the dynamic equilibrium condition {omega} {wedge} g = 0 (i.e., keep head and body angular velocity as much as possible aligned with gravity). Clearly, rotations about roll and/or pitch body axes are much more threatening for balance control than yaw movements because of the nearly orthogonal orientation of these movements relative to gravity. To maintain balance these movements must be limited in amplitude. Along these lines, we propose that the projection mechanism observed in head-in-space orientation signals reflects the signal processing of a central vestibulomotor balance control system. In contrast to roll and pitch, yaw head and body movements are much less threatening because the rotation axis is usually close to earth vertical. To establish a space-fixed relation between self-motion signals and gravity, it is more useful to use a rotation in this case. The rotation is optimally suited to transform velocity signals from head- to space-fixed coordinates because it is symmetric with respect to the longitudinal head and body axis, and therefore also with respect to the usual orientation relative to the earth horizontal plane.


    APPENDIX
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Interaction of otolith and vertical semicircular canal signals

In the present experiments, monkeys were rotated either in roll (about the naso-occipital or x-axis) or in pitch (about the interaural or y-axis). These rotations were applied with the axis in the earth-horizontal plane (90° OVAR), or tilted nose up or right ear up by 45° relative to the earth-horizontal plane (45° OVAR). Altogether the following 4 paradigms were tested: 1) 90° roll OVAR (see Fig. 1A, {vartheta} = 0), 2) 45° roll OVAR (Fig. 1B, {vartheta} = 45°), 3) 90° pitch OVAR (Fig. 1C, {vartheta} = 0), and 4) 45° pitch OVAR (Fig. 1D, {vartheta} = 45°). In the following we describe the central interaction of vertical semicircular canal and otolith signals when the rotation stops in different head orientations.

The integration of otolith and semicircular canal-borne head velocity signal can be characterized by a set of linear 1st-order differential equations. Generalizing Robinson’s 1D velocity storage model (Robinson 1977Go) to 3D we obtain the following set of differential equations in the time domain

(A1)
The system matrix F is obtained from a diagonal matrix A of 3 equal time constants [A]ii = 1/T (i = 1 ... 3) and a diagonal matrix K with 3 different feedback constants [K]ii = ki (i = 1 ... 3). The driving inputs to this leaky integration network are the lumped angular velocity signals from the vertical semicircular canals, represented by the vector {omega}rlv. Note that the input {omega}rlv is not modified, that is, {omega} = {omega}rlv, if the feedback is turned off.

PROJECTION MODEL. We have tested the hypothesis that the spatial transformation in this network is realized by a gravity-dependent projection mechanism (cosine tuning). Such a mechanism modifies the velocity signals in the loop according to

(A2)
where g is the estimated vector of gravity g. Thus the output of the spatial filter P is a vector directed along gravity with magnitude proportional to the cosine of the angle subtended by the input vector and gravity. The effect of this mechanism is to facilitate central integration of vertical canal signals according to their alignment with gravity.

By introducing the projection mechanism (Eq.A2) into the feedback system (Eq.A1) we obtain the following system matrix

(A3)
The eigenvectors and eigenvalues of the matrix F describe the spatiotemporal characteristics of A1

(A4)

(A5)
Because of the algebraic multiplicity of the eigenvalue {lambda}1, the eigenvectors (Eq.A4) must be replaced by generalized eigenvectors to correctly capture the spatial properties of the system in Eq.A1. Using the Jordan decomposition we find that eigenvector v1 needs to be replaced by

(A4')
Depending on the feedback, the eigenvector v3 tends to align with gravity, whereas the other 2 are always perpendicular to gravity (v1 · g = 0; v2 · g = 0). The time constant of the response component along the eigenvector v3, T3 = 1/{lambda}3 (Eq. A4) depends on both the feedback in the system and the system’s orientation relative to gravity, whereas the time constants along the other 2 eigenvectors v1 and v2 are invariable and equal. As the feedback constants approach unity, the system enhances responses along the eigenvector v3 by the factor 1/(1 – {sum}i kigi2), whereas responses along orthogonal directions are not affected.

ROTATION MODEL. Along the same lines, we tested the alternative hypothesis that the spatial transformation is a gravity-dependent rotation of postrotatory eye velocity. In this case, the system matrix F

(A6)
represents a rotation defined by the axis of rotation a = [a1 a2 a3]T (||a|| = 1) and the angle {vartheta} of rotation about this axis. The eigenvectors and eigenvalues of matrix F are (if a1 and a2 are not both equal zero)

(A7)

(A8)
The spatiotemporal properties of this system are characterized by one real eigenvector (v3 = a) and eigenvalue ({lambda}3) and two complex conjugate eigenvectors and eigenvalues (v1, {lambda}1 and v2, {lambda}2). We assume that the rotation axis a aligns with the line of intersection of the earth horizontal plane and the plane of head rotation (see Jaggi-Schwarz et al. 2000Go). Accordingly, we aligned the eigenvector v3 in the roll OVAR paradigm with the y-axis; that is, v3 = a = [0 cos ({varphi}a) sin ({varphi}a)]T with {varphi}a = 0 in the upright position. Similarly in the pitch OVAR paradigm we aligned it with the x-axis; that is, a = [cos ({varphi}a) 0 –sin ({varphi}a)]T with {varphi}a = 0 in the upright position. Thus when the subject stops after roll OVAR in 60° left ear-down position, the rotation axis, about which eye velocity is assumed to rotate through an angle {vartheta}a toward alignment with gravity, would be: a = [0 cos (60°) sin (60°)]T. The complete reorientation of postrotatory eye velocity toward gravity is described by the angles {varphi}a and {vartheta}a.

The solution of Eq.A1 with the gravity-dependent system matrix F (Eq.A3 or Eq.A6) yields for an initial response vector {omega}0 (={omega}oto) (Kailath 1980Go)

(A9)
This solution constitutes a purely otolith-dependent and semicircular canaldependent component that in turn also depends on gravity by virtue of the gravity-dependent matrix exponential etF (Eq.A3). It is easily seen that this matrix exponential has the same eigenvalues and eigenvectors as the matrix F because of the relation

(A10)
where {lambda}i are the eigenvalues and V is the matrix of eigenvectors of F.

For a parametric evaluation of the matrix elements of F we note that during roll and pitch OVAR, the vector of gravity rotates about the subject as follows

(A11)

(A12)
The angle {varphi} describes the angle of rotation of the vector of gravity about the naso-occipital or interaural axis (positive for leftward or backward rotation), whereas the angle {vartheta} represents the tilt of the vector of gravity relative to the respective rotation axis (0 ≤ {vartheta} ≤ 180°). In the 90° roll or pitch OVAR paradigm {vartheta} = 90°, whereas in the 45° roll or pitch OVAR {vartheta} = 135° and 45°. In the projection model, Eqs.A11 and A12 determine the matrix elements of the projection operator P (Eq.A3) at stop of rotation as functions of {varphi} and {vartheta}. In the rotation model, the same angles {varphi} and {vartheta} describe the orientation of the axis a (i.e., {varphi}a = {varphi}) and the amount of rotation of postrotatory eye velocity about this axis, away from the head rotation axis (i.e., {vartheta}a = {vartheta}). This is so because the vectors g and a have been defined such that they are always perpendicular to each other.

For the experimental evaluation of Eq.A9 we approximated the postrotatory vertical semicircular canal signals at stop of rotation by {omega}rlv(t) = v0(e–t/T1 – e–t/T0)rlv with T1 = 5 s, T0 = 0.003 s, and unit vector rlv directed along the naso-occipital and interaural axis for roll and pitch OVAR, respectively. With this information at hand, we analytically solved the time integral on the right-hand side of Eq.A9. For a detailed quantitative evaluation, we fitted the solution of Eq.A9 with nonlinear least-squares methods to the postrotatory response. The estimated parameters were the tilt angles {vartheta} and {varphi} (gi in Eqs.A11 and A12); the inertial response time constant T (Eq.A5); the feedback constants k1, k2, and k3; and the sensitivity of the semicircular canal response (v0).


    GRANTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This work was supported by the Swiss National Science Foundation Grant 31-47 287.96.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank B. Disler and E. Buffone for excellent technical assistance.


    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: B.J.M. Hess, Frauenklinikstrasse 26, CH-8091 Zürich, Switzerland (E-mail: bhess{at}neurol.unizh.ch)


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Angelaki DE. Detection of rotating gravity signals. Biol Cybern 67: 523–533, 1992a.[CrossRef][Web of Science][Medline]

Angelaki DE. Spatio-temporal convergence (STC) in otolith neurons. Biol Cybern 67: 83–96, 1992b.[CrossRef][Web of Science][Medline]

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 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. I. Linear acceleration responses during off-vertical axis rotation. J Neurophysiol 75: 2405–2424, 1996a.[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, 1996b.[Abstract/Free Full Text]

Angelaki DE, Hess BJM, and Suzuki J. 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]

Bevington PR and Robinson DK. Data Reduction and Error Analysis for the Physical Sciences. Boston, MA: McGraw-Hill, 1992.

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

Brandt T. Vertigo: Its Multisensory Syndromes. London: Springer-Verlag, 1999.

Cohen B, Matsuo V, and Raphan T. Quantitative analysis of the velocity characteristics of optokinetic nystagmus and optokinetic after-nystagmus. J Physiol 270: 321–344, 1977.[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]

Droulez J and Darlot C. The geometric and dynamic implications of the coherence constraints in three-dimensional sensorimotor interactions. In: Attention and Performance XIII. Motor Representation and Control, edited by Jeannerod M. Hillsdale, NJ: Erlbaum, 1990, p. 495–526.

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, 1976.[Abstract/Free Full Text]

Gizzi M, Raphan T, Rudolph S, and Cohen B. Orientation of human optokinetic nystagmus to gravity: a model-based approach. Exp Brain Res 99: 347–360, 1994.[Web of Science][Medline]

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

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

Green AM and Angelaki DE. An integrative neural network for detecting inertial motion and head orientation. J Neurophysiol March 31, 2004; 10.1152/jn.01234.2003b.

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]

Harris LR and Barnes GR. Orientation of vestibular nystagmus is modified by head tilt. In: The Vestibular System: Neurophysiologic and Clinical Research, edited by Graham MD and Kemink JL. New York: Raven Press, 1987, p. 539–548.

Hess BJM and Angelaki DE. Inertial representation of visual and vestibular self-motion signals. In: Multisensory Control of Posture, edited by Mergner T and Hlavacka F. New York: Plenum Press, 1995, p. 183–190.

Hess BJM and Dieringer N. Spatial organization of the maculo-ocular reflex of the rat: responses during off-vertical axis rotation. Eur J Neurosci 2: 909–919, 1990.[CrossRef][Web of Science][Medline]

Jaggi-Schwarz K and Hess BJM. Influence of dynamic tilts on the perception of earth-vertical. Exp Brain Res 149: 340–350, 2003.[Web of Science][Medline]

Jaggi-Schwarz K, Misslisch H, and Hess BJM. Canal–otolith interactions after off-vertical axis rotations. I. Spatial reorientation of horizontal vestibuloocular reflex. J Neurophysiol 83: 1522–1535, 2000.[Abstract/Free Full Text]

Jaggi-Schwarz K, Ortega M, and Hess BJM. Reciprocal error behavior in estimated body position and subjective visual vertical. Exp Brain Res 150: 122–125, 2003.[Web of Science][Medline]

Kailath T. Linear Systems. Englewood Cliffs, NJ: Prentice Hall, 1980.

Lindeman HH. Anatomy of the otolith organs. Adv Otorhinolaryngol 20: 405–433, 1973.[Medline]

Loe PR, Tomko DL, and Werner G. The neural signal of angular head position in primary afferent vestibular nerve axons. J Physiol 230: 29–50, 1973.[Abstract/Free Full Text]

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, Young LR, Paige GD, and Tomko DL. Three dimensional eye movements of squirrel monkeys following postrotatory tilt. J Vestib Res 3: 123–139, 1993.[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 Gifford CA. Neural processing of gravito-inertial cues in humans. II. Influence of the semicircular canals during eccentric rotation. J Neurophysiol 85: 1648–1660, 2001.[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]

Moler C and Van Loan C. Nineteen dubious ways to compute the exponential of a matrix. SIAM Rev 20: 801–837, 1978.[CrossRef]

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

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. Velocity storage in the vestibulo-ocular reflex arc (VOR). Exp Brain Res 35: 229–248, 1979.[Web of Science][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]

Robinson DA. Vestibular and optokinetic symbiosis: an example of explaining by modeling. In: Control of Gaze by Brain Stem Neurons. Developments in Neuroscience, edited by Baker H and Berthoz A. Amsterdam: Elsevier/North-Holland Biomedical Press, 1977, p. 49–58.

Schnabolk C and Raphan T. Modeling 3-D slow phase velocity estimation during off-vertical-axis rotation (OVAR). J Vestib Res 2: 1–14, 1992.[Medline]

Spoendlin HH. The ultrastructure of the vestibular sense organ. In: The Vestibular System and Its Diseases, edited by Wolfson RJ. Philadelphia, PA: Univ. of Pennsylvania Press, 1966, p. 39–68.

Steinhausen W. Über die Beobachtung der Cupula in den Bogengangsampullen des Labyrinths des lebenden Hechts. Pfluegers Arch 232: 500–512, 1933.[CrossRef]

Udo de Haes HA and Schone H. Interaction between statolith organs and semicircular canals on apparent vertical and nystagmus. Investigations on the effectiveness of the statolith organs. Acta Otolaryngol 69: 25–31, 1970.[Medline]

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.

von Holst E and Grisebach E. Einfluss des Bogengangssystems auf die "subjektive Lotrechte" beim Menschen. Naturwissenschaften 38: 67–68, 1951.[CrossRef][Web of Science]

Wilson VJ and Melvill-Jones JG. Mammalian Vestibular Physiology. New York: Plenum, 1979.

Young LR. Spatial orientation. In: Principles and Practice of Aviation Psychology (Volume in the Human Factors in Transportation Series), edited by Tsang PS and Vidulich MA. Mahwah, NJ: Erlbaum, 2002, p. 1–41.

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 1. Biol Cybern 86: 209–230, 2002.[CrossRef][Web of Science][Medline]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
93/3/1633    most recent
00383.2004v1
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 Google Scholar
Google Scholar
Right arrow Articles by Hess, B. J. M.
Right arrow Articles by Misslisch, H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hess, B. J. M.
Right arrow Articles by Misslisch, H.


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