Skip to main content
Advertisement
  • Loading metrics

A triple distinction of cerebellar function for oculomotor learning and fatigue compensation

  • Jana Masselink ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    jana.masselink@uni-muenster.de

    Affiliation Institute for Psychology & Otto Creutzfeldt Center for Cognitive and Behavioral Neuroscience, University of Münster, Münster, Germany

  • Alexis Cheviet,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Validation, Writing – review & editing

    Affiliations IMPACT Team, Lyon Neuroscience Research Center, University Claude Bernard Lyon 1, Bron cedex, France, Department of Psychology, Durham University, South Road, Durham, United Kingdom

  • Caroline Froment-Tilikete,

    Roles Conceptualization, Data curation, Methodology, Validation, Writing – review & editing

    Affiliations IMPACT Team, Lyon Neuroscience Research Center, University Claude Bernard Lyon 1, Bron cedex, France, Hospices Civils de Lyon—Pierre-Wertheimer Hospital, Neuro-Ophtalmology Unit, Bron cedex, France

  • Denis Pélisson,

    Roles Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Validation, Writing – review & editing

    Affiliation IMPACT Team, Lyon Neuroscience Research Center, University Claude Bernard Lyon 1, Bron cedex, France

  • Markus Lappe

    Roles Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Validation, Writing – review & editing

    Affiliation Institute for Psychology & Otto Creutzfeldt Center for Cognitive and Behavioral Neuroscience, University of Münster, Münster, Germany

Abstract

The cerebellum implements error-based motor learning via synaptic gain adaptation of an inverse model, i.e. the mapping of a spatial movement goal onto a motor command. Recently, we modeled the motor and perceptual changes during learning of saccadic eye movements, showing that learning is actually a threefold process. Besides motor recalibration of (1) the inverse model, learning also comprises perceptual recalibration of (2) the visuospatial target map and (3) of a forward dynamics model that estimates the saccade size from corollary discharge. Yet, the site of perceptual recalibration remains unclear. Here we dissociate cerebellar contributions to the three stages of learning by modeling the learning data of eight cerebellar patients and eight healthy controls. Results showed that cerebellar pathology restrains short-term recalibration of the inverse model while the forward dynamics model is well informed about the reduced saccade change. Adaptation of the visuospatial target map trended in learning direction only in control subjects, yet without reaching significance. Moreover, some patients showed a tendency for uncompensated oculomotor fatigue caused by insufficient upregulation of saccade duration. According to our model, this could induce long-term perceptual compensation, consistent with the overestimation of target eccentricity found in the patients’ baseline data. We conclude that the cerebellum mediates short-term adaptation of the inverse model, especially by control of saccade duration, while the forward dynamics model was not affected by cerebellar pathology.

Author summary

Achieving a fine-grained understanding of how the cerebellum continuously recalibrates our movements is an ongoing challenge in sensorimotor neuroscience. Recently, we showed that recalibration of saccadic eye movements does not only operate in motor space, i.e. by adjusting the motor command, but also in external and internal visual space, i.e. by adjusting the spatial representation of the target and the internal saccade size. For this purpose, (1) we developed a paradigm that allowed us to monitor changes of the internal saccade size estimated from trans-saccadic target localizations, and (2) we unified the three learning processes in one computational modeling framework. Here we apply this approach to the saccade learning data of patients with a neurodegenerative cerebellar disease. First, we dissociate the cerebellar role in recalibration of the three sites of learning. Second, we show how learning is transposed to saccade kinematics. Third, we provide first insights into the perceptual consequences of cerebellar pathology that, according to our model, may be a mechanism to recover from disease-specific motor deficits. Our modeling framework may help to dissociate the contribution of specific sensorimotor areas to adaptive behavior as well as to improve the understanding of learning deficits and compensatory strategies in the clinical context.

Introduction

The cerebellum modulates supervised error-based learning [13]. If a movement is inaccurate, an error signal adapts the synaptic gains of cerebellar Purkinje cells [47] that implement an inverse model, i.e. the mapping of the visuospatial movement goal onto a motor command [4, 8, 9]. If e.g. the eye muscles fatigue, the synaptic gains are modulated such that the saccade does not fall short of the target [1012]. To compensate for continuous fluctuations of muscle dynamics, the cerebellum is constantly recalibrating the saccadic circuitry [1315].

However, the traditional view on the cerebellum as a pure motor hub has changed. The cerebellum is interconnected with almost the entire cortex [1618], including retinotopic areas [1921], and basal ganglia [18], and holds at least five visuospatial maps itself [22]. Moreover, it has been shown that saccadic motor learning also comprises changes in visuospatial representations. These changes occur, firstly, in the visuospatial target map when errors are assigned to an internal failure of spatial target representation [2327]. Secondly, these changes occur in the internal representation of saccade size as measured by trans-saccadic target localizations [25, 28, 29].

Changes in the internal saccade representation affect the forward dynamics model that transforms the corollary discharge, i.e. a copy of the motor command (CDM) [3032] into a visuospatial representation of saccade size (CDV) [4, 33, 34]. Intact CDV information is necessary to correctly bridge the spatial gap between pre- and post-saccadic visual input that supports our perception of a stable world [3537]. In Masselink et al. (2021) [38], we implemented the action-perception entanglement during learning in a computational model that collectively adapts the visuospatial target map, the inverse model and the forward dynamics model in response to postdictive motor error, i.e. motor error with respect to a postdicted target position.

Cerebellar perturbations restrain saccade changes during learning, as shown e.g. in lesioned monkeys [39, 40], in cerebellar patients [4143] and in healthy subjects with non-invasive cerebellar stimulation [4446]. Because of its circular learning architecture, its visuospatial maps and its CD projection pathway to frontal cortex [47], the cerebellum is an ideal candidate to also recalibrate the visuospatial target map and the forward dynamics model.

In Cheviet et al. (2022) [48] we found first empirical evidence that changes in visual target localizations during saccadic motor learning are impaired in cerebellar patients. These target localizations were performed (1) during fixation and (2) trans-saccadically after a saccade landed in the dark. The dataset includes eight patients with a neurodegenerative cerebellar disease and eight healthy control subjects who performed three conditions comprising (1) learning from peri-saccadic inward and (2) from peri-saccadic outward target steps, and (3) compensation of neural oculomotor fatigue, i.e. without target step. As learning relies on plasticity of the transformation processes between visual and motor representations (i.e. on how the target is scaled on the visuospatial map, on the inverse model and on the forward dynamics model), a quantification of the relative impairments of these transformations is essential to draw conclusions on cerebellar function.

Here we apply an expanded version of the visuomotor learning model of Masselink et al. (2021) [38] to the empirical data of Cheviet et al. (2022) [48] in order to quantify the recalibration of the visuospatial target map, of the inverse model and of the forward dynamics model in cerebellar patients. To achieve a deeper understanding of which control processes are perturbed, we expanded the model of Masselink et al. (2021) [38] with (1) how changes in the motor command are transposed to saccade kinematics (i.e. saccade peak velocity and saccade duration), and (2) how the visuomotor system compensates for a fatigue-induced decline in saccade peak velocity. Our results reveal that cerebellar pathology reduces short-term learning of the inverse model, especially when controlled by upregulation of saccade duration. Notably, the forward dynamics model and thus, the CDV signal, seemed correctly informed about the reduced saccade change. Learning of the visuospatial target map trended in learning direction only in control subjects, yet without significance. Moreover, we show that cerebellar pathology is accompanied by an overestimation of target eccentricity that, according to our model, could stem from error reduction to counteract oculomotor fatigue on a long timescale. We conclude that intact recalibration of the inverse model, particularly via control of saccade duration, depend on cerebellar integrity. The forward dynamics model, at least for learning from post-saccadic errors, is not affected by cerebellar pathology. Long-term learning at perceptual level may partially compensate for cerebellar motor deficits.

Results

To reveal how adaptation of the visuospatial target map, of the inverse model and of the forward dynamics model for saccade motor learning relies on cerebellar integrity, we fitted a visuomotor learning model to the saccade and visual target localization data of Cheviet et al. (2022) [48], including eight cerebellar patients and eight healthy control subjects. Table 1 lists the characteristics of the cerebellar patients, including disease type, disease duration and anatomical involvements.

thumbnail
Table 1. Disease types and anatomical involvements of cerebellar patients.

https://doi.org/10.1371/journal.pcbi.1011322.t001

Experimental data

Each subject participated in three conditions requiring saccades to a 20° rightward target as well as pre- and post-saccadic visual localizations of a shortly flashed bar presented around the usual target position (Fig 1A). In the pre-saccadic localization trials, subjects localized the flash with a mouse cursor while holding gaze at the fixation point. In the post-saccadic localization trials, subjects performed a saccade to a 20° rightward target and then localized the pre-saccadic flash that had appeared before target onset. The pre- and the post-exposure phase (phase 1 and 3) measured subjects’ state of saccade vector, pre- and post-saccadic visual localizations before and after exposure. The exposure phase (phase 2) induced either (1) saccade shortening in response to a 6° peri-saccadic inward target step (inward condition), (2) saccade lengthening in response to a 6° peri-saccadic outward target step (outward condition), or (3) tested the ability to maintain the saccade vector during repetitive, stereotyped saccades to a non-stepping target (no step condition). This task is known to induce oculomotor fatigue, i.e. a decline in saccade peak velocity. In healthy subjects, this is usually counteracted by upregulation of saccade duration such that the saccade vector stays stable [4955]. Compared to learning in response to an assumed fatigue of the eyes muscles, oculomotor fatigue is of neural origin and perhaps due to a loss of motivation or attention [55]. However, patients with cerebellar lesion of the vermis have shown to be impaired in compensating oculomotor fatigue [42].

thumbnail
Fig 1. Experimental tasks and model framework.

(A) In the saccade trials, subjects executed a saccade to a 20° rightward target. In the pre-exposure phase, the target was extinguished at saccade onset. In the exposure and post-exposure phase, the target either stepped 6° inward (inward condition), 6° outward (outward condition) or stayed at its initial position (no step condition). In the pre-saccadic localization trials, subjects localized a 12 ms flash with a mouse cursor while holding gaze at the fixation point. In the post-saccadic localization trials, subjects performed a saccade to a 20° rightward target and then localized the 12 ms pre-saccadic flash with a mouse cursor. Please note that the stimuli are not drawn to scale. The mouse cursor was a blue line pointer. The yellow circle illustrates gaze location but was not present at the stimulus display. (B) The target with the physical distance P1 is represented at the location V1 on the visuospatial map. An inverse model maps V1 onto a motor command M. Before saccade start, a forward dynamics model transforms a copy of the motor command CDM into visuospatial coordinates, i.e. into the computed displacement of visual space CDV. A forward outcome model then shifts retinal coordinates by CDV to predict the visual location of the post-saccadic target . For saccade execution, the motor command is transposed to saccade peak velocity κ and saccade duration λ, producing the saccade vector PM. The actual post-saccadic target appears at retinal position V2. A backward outcome model then postdicts the visual post-saccadic target back to pre-saccadic space . The visuomotor system evaluates the accuracy of the motor command with respect to the postdicted target position (Epost) in order to adapt its gains ωv, ωm and ωcd. A decrease of ωm (inward learning) is controlled by downregulating saccade peak velocity while an increase of ωm (outward learning) is controlled by upregulation of saccade duration. When saccades are performed repetitively to the same, non-stepping target, oculomotor fatigue occurs, i.e. peak velocity declines, which is usually compensated by an increase of saccade duration.

https://doi.org/10.1371/journal.pcbi.1011322.g001

Fig 2A shows the inward learning data of an example control subject and Fig 2B shows the inward learning data of an example patient. Amplitude adaptation is impaired in the patient compared to the control subject. Fig 2C shows that saccade variability was higher in patients than in control subjects. An increased variability of saccade endpoints has been reported before as a side-effect of cerebellar dysfunction [40, 42, 43, 56]. The standard error of saccade vectors in the pre-exposure phase averaged across the three conditions was 0.28 ± 0.03° in control subjects and 0.54 ± 0.06° in patients (difference between patients and control subjects t14 = -3.78, p = .002). Moreover, the mean saccade vector in the pre-exposure phase averaged across the three conditions differed more between patients than between control subjects (cross-subject standard error of saccade vectors, control subjects 0.37°, patients 0.57°). Fig 3 shows the individual saccade peak velocities and Fig 4 the individual saccade durations in dependence on saccade vector for control subjects (C1-C8) and patients (P1-P8) from early trials (blue) to late trials (yellow).

thumbnail
Fig 2. Example subject data for the inward condition and within-subject saccade variability in the pre-exposure phase.

Saccade vectors, pre- and post-saccadic localizations for the inward condition of (A) an example control subject and (B) an example patient. The pre-exposure phase measured subjects’ baseline state without target step. The exposure phase induced saccadic learning and the post-exposure phase measured adaptation-induced changes to the saccade vector, pre- and post-saccadic localizations. Similar to the pattern that we found at group level, the patient shows less learning in the saccade vector and higher saccade endpoint variability compared to the control subject. In the model, the predicted post-saccadic target is derived from post-saccadic localization relative to saccade landing. (C) Within-subject standard error of saccade vectors in the pre-exposure phase (averaged across the three conditions per subject) compared between control subjects and patients, *** p <.001, ** p <.01, * p <.05 and n.s. p ≥.05.

https://doi.org/10.1371/journal.pcbi.1011322.g002

thumbnail
Fig 3. Individual saccade peak velocities in dependence on saccade vectors from early trials (blue) to late trials (yellow).

Saccade peak velocities and saccade vectors are shown separately for each subject and condition (controls C1-C8, patients P1-P8 including disease type and disease duration in years). During inward learning of control subjects, the decrease of the saccade vector encompasses a decrease of saccade peak velocity. This is in line with saccade shortening being mainly controlled by downregulating peak velocity. During outward learning, the saccade vector increases while peak velocity stays stable. This is in line with saccade lengthening being mainly controlled by upregulating saccade duration (see Fig 4). In the no step condition, saccade peak velocity encompasses a small but substantial decline that, however, is not accompanied by a decline of the saccade vector. This is a sign of oculomotor fatigue being successfully compensated by healthy subjects. As known from cerebellar patients, the patients’ saccade vectors and peak velocities are more variable than in control subjects. Moreover, overall learning effects on saccade vector and peak velocity are smaller. In the no step condition, patients show substantial oculomotor fatigue, i.e. a decline in peak velocity that is not fully compensated by saccade duration.

https://doi.org/10.1371/journal.pcbi.1011322.g003

thumbnail
Fig 4. Individual saccade durations in dependence on saccade vectors from early trials (blue) to late trials (yellow).

Saccade durations and saccade vectors are shown separately for each subject (controls C1-C8, patients P1-P8) and condition. During inward learning of control subjects, saccade duration declines in subjects C3, C5, C7 and C8 (yet not significantly across subjects when compared between the pre- and the post-expose phase, t7 = -1.01, p = .346). The decrease of the saccade vector rather stems from a decline of saccade peak velocity (see Fig 3). In the outward condition, saccade duration increases together with the saccade vector. This is consistent with saccade lengthening being mainly controlled by upregulation of saccade duration. In the no step condition, saccade duration increases, thereby compensating for the peak velocity loss seen in Fig 3. In patients, small learning effects can be seen but saccade vectors and durations are more noisy. In the no step condition, saccade duration is increased only in patient P4 and P7. Overall, duration cannot sufficiently counteract the peak velocity loss.

https://doi.org/10.1371/journal.pcbi.1011322.g004

Please note that because the analysis of the empirical data has been presented in Cheviet et al. (2022) [48], this paper focuses on a model-based analysis.

Model fits

Fig 1B presents the framework of the expanded visuomotor learning model based on Masselink et al. (2021) [38] that we fitted to the saccade vectors and pre- and post-saccadic visual localizations, separately for each group (controls, patients) and each condition (inward, no step, outward). The visual pre-saccadic target V1 (fitted to the pre-saccadic visual localizations) is represented on the visuospatial map with the visual gain ωv (Fig 1B-1). An inverse model transforms V1 into a motor command M (fitted to the saccade vectors) with the motor gain ωm (Fig 1B-2). Before saccade execution, a copy of the motor command is routed into the CD pathway where a forward dynamics model estimates the visuospatial size CDV of the upcoming saccade based on the CD gain ωcd (Fig 1B-3). A forward outcome model then shifts retinal coordinates by CDV to predict the retinal position of the post-saccadic target (fitted to the post-saccadic visual localizations; Fig 1B-4). For saccade execution, the motor command is transposed to saccade peak velocity κ and saccade duration λ (Fig 1B-5). After the saccade is performed (Fig 1B-6), the actual post-saccadic target appears at retinal position V2 (Fig 1B-7). A backward outcome model then postdicts the visual post-saccadic target back to pre-saccadic space (Fig 1B-8). The postdictive motor error Epost evaluates the accuracy of the motor command with respect to the postdicted target position in order to adapt its gains ωv, ωm and ωcd (Fig 1B-9). A decrease of ωm (inward learning) is controlled by downregulating saccade peak velocity while an increase of ωm (outward learning) is controlled by upregulation of saccade duration (Fig 1B-5).

The no step condition in which saccades are performed repetitively to the same, non-stepping target should lead to oculomotor fatigue, i.e. a decline in saccade peak velocity which is usually compensated by an increase of saccade duration. We captured peak velocity and duration changes in the no step condition by fitting a peak velocity decay rate γκ ≥ 0 and a duration compensation rate γλ (with 0 ≤ γλ ≤ 1). The latter describes the within-trial percentage by which saccade duration compensates for the peak velocity loss.

Fig 5 presents the model fits to the mean data of control subjects (Fig 5A) and patients (Fig 5B) for the three conditions (inward, no step, outward). Dots with errorbars show the mean data of the pre- and post-exposure phases, respectively, jagged blue lines show the trial-by-trial saccade data (Fig 5A-1 and 5B-1: saccade vectors; Fig 5A-5 and 5B-5: saccade peak velocities; Fig 5A-6 and 5B-6: saccade durations) and smooth lines show the trial-by-trial model fit to the data. Fig 6 depicts the baseline state as well as pre- to post-exposure changes for each condition, compared between patients and control subjects. V1 corresponds to the pre-saccadic target localization, M is equivalent to the measured saccade vector and CDV relies on the interplay of post-saccadic target localization, pre-saccadic target localization and the measured saccade vector.

thumbnail
Fig 5. Model fits to the group data.

Pre- and post-exposure mean ± standard error (dots with error bars), trial-by-trial saccade data of the exposure phase (jagged blue lines) and model fits to the data (smooth lines) for (A) control subjects and (B) cerebellar patients, separately for each condition (inward, no step, outward). (1) Visual pre-saccadic target V1 (fitted to pre-saccadic localizations), motor command M (fitted to saccade vectors), computed displacement of visual space CDV, postdicted pre-saccadic target . (2) Predicted post-saccadic target (fitted to post-saccadic localizations with respect to saccade landing), visual post-saccadic target V2. (3) Visual gain ωv, motor gain ωm, CD gain ωcd. (4) Postdictive motor error Epost. (5) Saccade peak velocity κ. (6) Saccade duration λ. Asterisks indicate significant difference between the pre- and post-exposure phase with *** p <.001, ** p <.01, * p <.05 and n.s. p ≥.05.

https://doi.org/10.1371/journal.pcbi.1011322.g005

thumbnail
Fig 6. Baseline state and pre- to post-exposure changes.

(A) Baseline state of the visuomotor system averaged across the pre-exposure phases of the three conditions (mean ± standard error), separately for control subjects and patients. Group-centered asterisks indicate significant difference from 20° (V1, M, CDV), from 1 (ωv, ωm, ωcd) or from zero (Epost). (B) Pre- to post-exposure changes in the inward condition. (C) Pre- to post-exposure changes in the no step condition. (D) Pre- to post-exposure changes in the outward condition. For B-D group-centered asterisks indicate significant difference from zero. Panel-centered asterisks above a horizontal indicate significant difference between control subjects and patients with *** p <.001, ** p < .01, * p <. 05 and n.s. p ≥ .05.

https://doi.org/10.1371/journal.pcbi.1011322.g006

Table 2 gives an overview on the fitted parameters and residual standard errors as a measure of goodness of fit.

Differences in learning of the visuospatial target map between patients and controls

The fit of the visual pre-saccadic target position V1 to the pre-saccadic target localizations (Fig 5A-1 and 5B-1, green) captures how the visuospatial target map, i.e. the visual gain ωv (Fig 5A-3 and 5B-3, green), adapts on a short time-scale when errors are credited to an internal failure of visuospatial target representation. A small change of V1 in learning direction for control subjects could be observed, however this change did not reach significance (inward -0.7 ± 2.0°, t7 = -1.00, p = .351; outward 1.2 ± 1.8°, t7 = 1.86, p = .105; Fig 5A-1). As the change in pre-saccadic target localization across learning is a stable, substantial effect that, however, is usually rather small [2427, 38, 57], it seems not surprising that we could not replicate a significant effect with a sample size of eight subjects. Please note that in the analysis of Cheviet et al. (2022) [48] (which comprised more trials including those in which the flash was presented with a small distance to the 20° target position), these changes were subtly significant. Accordingly, the change of the visual gain ωv of control subjects was -0.04 ± 0.10 (t7 = -1.00, p = .351; Fig 6B) in the inward condition and 0.06 ± 0.09 (t7 = 1.86, p = .105; Fig 6D) in the outward condition.

The visual pre-saccadic target position V1 of cerebellar patients did not show any change in the direction of the target step. Small changes were observed in the opposite direction (inward 0.6 ± 1.3°, t7 = 1.35, p = .218; outward -0.9 ± 1.1°, t7 = -2.28, p = .057; Fig 6B and 6D) but were also not significant in the analysis of Cheviet et al. (2022) [48]. Adding a mechanism that allows the model to learn in the opposite direction of error reduction seemed not justified in this context. Accordingly, we fitted a stable V1 and a stable ωv to the patients’ data (Fig 5B-1 and 5B-3, green). In sum, patients did not show short-term learning of the visual pre-saccadic target position V1 as compared to the tendency observed in control subjects. The V1 change was significantly different between patients and control subjects in the outward condition (t14 = 2.78, p = .015, Fig 6D) but failed to reach significance in the inward condition (t14 = -1.57, p = 0.138). This seems to rely on the higher adaptation level that control subjects reached in the outward condition compared to the inward condition. Differences in adaptation levels between target step directions have been reported before [12, 26, 58, 59].

As a tendency for a change of the visual gain in learning direction was observed in control subjects while remaining completely absent in patients, the cerebellum may play a role in the mediation of the visuospatial target map. However, it seems that a larger sample size is necessary to reliably examine whether cerebellar integrity plays a crucial role in this adaptive process.

Cerebellar pathology reduces short-term learning of the inverse model

The fit of the motor command M to the saccade vectors (Fig 5A-1 an d 5B-1, blue) quantifies how the inverse model, i.e. the motor gain ωm (Fig 5A-3 and 5B-3, blue), adapts when errors are assigned to false computation of the motor command to bring the eyes to the visuospatial target position. Hence, quantification of the motor gain change allows to evaluate inverse model adaptation in isolation, i.e. irrespective of the change of the visual pre-saccadic target representation V1 that is the input to the inverse model. In other words, the quantification of motor gain adaptation allows to differentiate to which cause errors are assigned to, i.e. to a false representation of the motor command M or to a false visuospatial representation of the target V1.

In accordance with previous studies [10, 11, 28, 38, 5962], the motor command M of control subjects decreased by -3.9 ± 1.8° during inward learning (t7 = -6.06, p <.001; Fig 6B) and increased by 3.6 ± 2.0° during outward learning (t7 = 5.10, p = .001; Fig 6D). Accordingly, the inverse model, i.e. the motor gain ωm, adapted significantly in both learning conditions (inward -0.17 ± 0.13, t7 = -3.71, p = .008; outward 0.12 ± 0.13, t7 = 2.50, p = .041). As expected by the model, the motor command changes during inward learning were transposed to saccade kinematics by downregulation of saccadic peak velocity (-94.9 ± 46.2, t7 = -5.81, p = .001; Figs 5A-5 and 3) while saccade duration stayed constant (-2.6 ± 7.3 ms, t7 = -1.01, p = .346; Figs 5A-6 and 4). By contrast, the motor command increase during outward learning was controlled by upregulation of saccade duration (10.4 ± 3.1 ms, t7 = 9.37, p <.001) while saccade peak velocity stayed constant (-11.0 ± 35.3, t7 = -0.89, p = .405). The saccade changes brought control subjects’ eyes significantly closer to the post-saccadic target (V2 change by 3.9 ± 1.8° during inward learning, t7 = 6.06, p <.001; -3.6 ± 2.0° during outward learning, t7 = -5.1, p = .001; Fig 5A-2). Thus, control subjects successfully reduced the postdictive motor error Epost in both learning paradigms (inward 3.8 ± 2.4°, t7 = 4.53, p = .003; outward -3.3 ± 2.2°, t7 = -4.16, p = .004, Fig 5A-4). This corresponds to 64.1 ± 53.0% error reduction in the inward condition (t6 = 3.2, p = .019; Fig 6B) and 50.1 ± 36.3% error reduction in the outward condition (t7 = 3.9, p = .006; Fig 6D).

Compared to control subjects, motor command changes were largely reduced in cerebellar patients. During inward learning, the motor command M decreased by -1.2 ± 2.1° and increased during outward learning by 0.8 ± 2.6° (Fig 5B-1). These changes were not significant (inward t7 = -1.65, p = .143; outward t7 = 0.83, p = .433) and significantly less than in control subjects (inward t14 = -2.75, p = .016; outward t14 = 2.43, p = .029; Fig 6B and 6D). Accordingly, inverse model adaptation was not significant in patients (ωm change inward -0.08 ± 0.14, t7 = -1.73, p = .127, outward 0.07 ± 0.13, t7 = 1.58, p = .158). Saccade peak velocity was significantly decreased in the inward condition (-55.9 ± 55.4, t7 = -2.85, p = .025; Figs 5B-5 and 3) but saccade duration was not significantly increased in the outward condition (4.1 ± 13.7 ms, t7 = 0.84, p = .427; Figs 5B-6 and 4). Hence, reduction of postdictive motor error was less effective than in the control subjects, i.e. only 5.5 ± 78.8% during inward learning (z = 0.84, p = .401; Fig 6B) and 24.5 ± 33.4% during outward learning (t7 = 2.07, p = .077; Fig 6D). Reduction of Epost in the inward condition was significantly less in patients than in controls (t14 = 2.15, p = .049).

In sum, the motor gain adapted significantly in learning direction in control subjects but not in patients. This demonstrates that the cerebellum has a crucial role in the adaptation of the inverse model, i.e. when errors are assigned to a false computation of the motor command such that the saccade is erroneously scaled with respect to the visuospatial target position.

Forward dynamics model is correctly informed about the reduced saccade change

Quantification of CDV in our model (Fig 5A-1 and 5B-1, red) relies on the difference between pre- and post-saccadic localization, hence capturing CDV-based integration between the pre- and the post-saccadic visual scene (see Fig 5A-2 and 5B-2 for changes of the post-saccadic target localization with respect to saccade landing, i.e. ). The CD gain ωcd captures how the forward dynamics model is calibrated when errors are credited to a false internal representation of the saccade. When pre- and post-saccadic visual scene are perfectly integrated, i.e. when the flash is localized at the same position pre- and post-saccadically, ωcd = 1 such that CDV matches M and, hence, correctly reflects the performed saccade.

Fig 7 shows that the motor command M and the CDV signal are highly correlated across control subjects and patients, respectively, in the pre-exposure phase (baseline state) as well as in the post-exposure phase of all three conditions (inward, no step, outward). This shows that the quantification of CDV based on our paradigm is a valid estimate of the internal representation of the saccade.

thumbnail
Fig 7. Correlation of motor command M and CDV signal across subjects at steady state.

At steady state, i.e. in the pre-exposure phase across all three conditions (baseline) and in the post-exposure phase of the inward condition, the no step condition and the outward condition, the motor command M and CDV signal highly correlate across control subjects and patients. This shows that the quantification of the CDV signal by the model based on the trans-saccadic target localization and saccade data is a valid estimate of the internal representation of the saccade.

https://doi.org/10.1371/journal.pcbi.1011322.g007

In control subjects, the CDV signal declined by -4.1 ± 2.1° during inward learning (t7 = -5.54, p <.001; Fig 6B) and increased by 3.9 ± 2.4° during outward learning (z = 2.38, p = .017; Fig 6D). Fig 5A-1 shows that CDV well reflects the actual saccade changes in both learning paradigms. Thus, the initial CD gain (ωcd = 1.03 ± 0.02, averaged as the baseline of all three conditions, Fig 6A) stayed stable throughout learning (inward 0.01 ± 0.08, t7 = 0.39, p = .707, Fig 6B; outward 0.01 ± 0.06, z = -0.42, p = .674, Fig 6D). This means that the forward dynamics model is well informed about the actual saccade changes in control subjects.

Surprisingly, this picture was similar in cerebellar patients. In both the inward and the outward condition, patients showed only a small change of the motor command M in learning direction that did not reach significance. This was well reflected by the CDV signal. During inward learning, CDV of patients was reduced by -1.3 ± 2.9° which did not reach significance (t7 = -1.28, p = .241, Fig 6B), thus well reflecting the actual saccade changes that were less than in control subjects (Fig 5B-1). In the outward condition, CDV of patients stayed roughly constant (change of -0.4 ± 3.5°, t7 = -0.35, p = .737, Fig 6D). Accordingly, the initial CD gain (ωcd = 1.04 ± 0.04, averaged as the baseline of all three conditions, Fig 6A) was not modified during inward learning (-0.01 ± 0.1, t7 = -0.19, p = .856) and only slightly modified during outward learning (-0.07 ± 0.07, t7 = -2.67, p = .032).

In sum, the forward dynamics model seemed correctly informed about the ongoing motor changes, similar in controls and patients (no difference in ωcd change between patients and controls, inward t14 = 0.39, p = .704, outward t14 = 0.39, p = .704). This suggests that the forward dynamics model that is used for learning from post-saccadic errors remains unaffected by cerebellar disease.

Cerebellar pathology is accompanied by overestimation of target eccentricity

Fig 6A presents the baseline state of control subjects and patients, averaged across the three conditions. While control subjects accurately localize the pre-saccadic target on the visual map (V1 = 20.3 ± 0.8°, t-test against P1 = 20°, t7 = 1.18, p = .276), patients showed a significant overestimation of target eccentricity by around 1.5° (V1 = 21.5 ± 1.3°, t-test against P1 = 20°, t7 = 3.13, p = .017). Hence, the visual gain ωv overestimated the actual target distance already when patients entered the laboratory (ωv = 1.07 ± 0.07, t-test against 1, t7 = 3.13, p = .017). This was consistent across all three conditions (Fig 5B-1). Still, patients exhibited a quite accurate baseline saccade (M = 19.5 ± 1.6°, t-test against P1 = 20°, t7 = -0.92, p = .389), i.e. the motor gain ωm was downregulated such that the saccade was not hypermetric as expected from the overestimated target eccentricity (0.91 ± 0.1, t-test against 1, t7 = -2.66, p = .033). Interestingly, target overestimation despite a quite accurate saccade was observed before in a patient with a right posterior ventrolateral and right ventromedial thalamic lesion [63]. In spite of the visuospatial hypermetria, the patients’ visuomotor system seemed to be at a steady state as the postdictive motor error Epost did not significantly differ from zero (t7 = 2.11, p = .073 in patients, t7 = 1.73, p = .127 in control subjects, Fig 6A).

Several patients tend to uncompensated oculomotor fatigue

The no step condition tested subjects’ ability to maintain a stable saccade vector when saccades were repetitively performed to the same, non-stepping target. Fig 5A-5 shows that control subjects underwent oculomotor fatigue, as revealed by the significant decline in saccade peak velocity by -45.7 ± 33.4 (t7 = -3.87, p = .006, γκ = 0.003). However, peak velocity decline was counteracted by a significant upregulation of saccade duration (5.6 ± 5.4 ms, t7 = 2.90, p = .023) such that the saccade vector was maintained (M remained stable with t7 = 0.30, p = .775). The individual main sequences, i.e. the relationship of saccade peak velocity and saccade duration to the saccade vector, are shown in Figs 3 and 4. Fig 3 depicts that the reduction of saccade peak velocity in the no step condition is not accompanied by a reduction of the saccade vector in control subjects. Accordingly, Fig 4 depicts the increase of saccade duration while the saccade vector stays stable. With a duration compensation rate of γλ = 0.964, saccade duration compensated for 96.4% of the within-trial velocity loss. Thus, control subjects successfully kept the visuomotor system calibrated and the error nullified (Fig 5A-4).

The cerebellar patients group experienced a loss of peak velocity by -68.2 ± 102.2 that, however, did not reach significance (t7 = -1.89, p = .101; γκ = 0.008). The high motor variability in conjunction with a small sample size of eight subjects make it difficult to draw clear conclusions from the patients’ data. Though, as the trend in peak velocity loss seemed reminiscent of the cerebellar patients of Golla et al. (2008) [42] (who observed reduced oculomotor fatigue compensation) and Xu-Wilson et al. (2009) [43] (who observed reduced within-saccade compensation for peak velocity fluctuations), we further investigated saccade peak velocity and saccade duration patterns in our sample. In our patients, upregulation of saccade duration could compensate only partially (7.7 ± 14.9 ms, t7 = 1.46, p = .187) such that the saccade vector declined by -0.8 ± 1.5° (t7 = -1.52, p = .172). Please note that also here, the patients’ data were highly variable and both effects did not reach significance at group level. Fig 6C shows that the saccade fell short in five of eight patients (saccade change from the pre- to the post-exposure phase of patient P1 -1.8°, P3 -1.4°, P6 -1.5°, P7 -2.7° and P8 -1.9°). This tendency is spared in patients P2, P4 and P5 (saccade change from the pre- to the post-exposure phase of patient P2 1.6°, P4 0.9° and P5 0.2°). Only one subject of the control group showed a comparable decline of -1.2° (subject C8; control subjects C1-C7 showed a saccade change of -0.4°, -0.3°, -0.1°, 0.2°, 0.3°, 1.0° and 1.0°). The individual main sequences of patients can be seen in Figs 3 and 4 including disease type and disease duration of each patient. Two of three patients with cerebellar ataxia (patients P1 and P3 but not patient P5) and one of two patients with Friedreich ataxia (patient P8 but not patient P2) show the impairment. Neither the disease type nor the specific anatomical involvements of the disease (Table 1) seem to be related to the impairment.

The fit of the saccade duration compensation parameter γλ to the patients’ data suggests that patients could compensate only 45.7% of the within-trial peak velocity loss (γλ = 0.457), resulting in uncompensated oculomotor fatigue (Fig 5B-1). Still, due to the lack of significance at group level, this result describes a tendency that needs further examination. Anyhow, the suggested deficit in saccade duration control during oculomotor fatigue appears consistent with the greater impairment during outward learning in which changes of the inverse model are usually transposed to an increase of saccade duration (Fig 5B-1 outward condition).

Overestimation of target eccentricity may be a result of error reduction to counteract uncompensated oculomotor fatigue

We wondered how the patients’ tendency for uncompensated oculomotor fatigue affected their natural visuomotor behavior in the long run. Saccade kinematics underlie continuous fluctuations and if saccade duration cannot sufficiently compensate for regular declines of saccade peak velocity, the visuomotor system will inevitably tend to be hypometric. Saccade hypometry has been observed as a consequence of vermal cerebellar lesion in monkeys [39, 40, 64, 65] and in humans [66]. However, Takagi et al. (1998) [39] and Barash et al. (1999) [40] reported that saccade vectors recovered within weeks to months. Notably, baseline saccades of cerebellar patients in our experiments were not hypometric (19.5 ± 1.6°, t7 = -0.92, p = .389) and not different from those of control subjects (t14 = 0.30, p = .767). Thus, if hypothetically, saccade hypometria occurred in the early phase of the disease, the visuomotor system would necessarily have experienced significant motor errors due to the saccade vector decline. Even if short-term error reduction is reduced in cerebellar patients, a certain ability for long-term learning is assumed to be preserved [40, 43, 56, 67]. Hence, we simulated how the visuomotor system will behave on a long time course, starting with the visuomotor gains ωv, ωm, ωcd, the learning rates αv, αm, αcd, the peak velocity decay rate γκ and the duration compensation rate γλ of healthy control subjects. Through the course of the simulation, the learning rates αv, αm, αcd, the peak velocity decay rate γκ and the duration compensation rate γλ successively changed towards the respective values obtained from the patients’ model fit.

Fig 8A shows a simulation for 23,3 Mio saccades directed towards a 10° goal. We chose the initial visuomotor gains ωv, ωm and ωcd measured for similarly sized saccades in 35 healthy subjects by Masselink et al. (2021) [38] because the larger sample size of Masselink et al. (2021) provides a more faithful picture than the initial gains of the healthy group of Cheviet et al. (2022) [48]. The simulation shows that a loss in saccade peak velocity and a reduced saccade duration compensation in the very early phase of the disease leads to an immediate decline of the motor command M, hence shortening the saccade vector from 9.8° to 8.5°. In turn, the visuomotor system starts to counteract the growing motor error Epost, inducing outward learning at perceptual (ωv) and motor level (ωm). Consequently, learning to counteract oculomotor fatigue leads to an outward shift of the visual pre-saccadic target localization and simultaneous recovery from saccade hypometry. The steady state at which the visuomotor system arrives is consistent with the baseline gains that we measured when patients came to the laboratory (residual standard error RSEω = 0.05). The simulation depicts a very fast progression of the peak velocity decay rate (scaled by νκ = 0.012) and the duration compensation rate (scaled by νλ = 0.040) towards patients’ values but the learning rates change very slowly across 23,3 Mio saccades (scaled by να = 1.7*10−7). Please note that most signals are shown across shorter timescales because they converge early.

thumbnail
Fig 8. Simulation of long-term visuomotor behavior after disease onset.

(A) Oculomotor fatigue in the very early phase of the disease causes a loss of saccade peak velocity κ that can only be partially compensated by upregulation of saccade duration λ. Consequently, a decline of the saccadic motor command M leads to saccade hypometry, that, in turn, results in a rise of motor errors Epost. Preserved long-term learning at perceptual (αv) and motor level (αm)—that, yet, slowly decays throughout disease progression—counteracts to keep saccadic behavior calibrated. Hence, the visuomotor system stabilizes with an overestimation of the pre-saccadic target eccentricity (V1) and recovery from saccade hypometria (M), matching the baseline visuomotor gains measured in cerebellar patients (ωv and ωm; ωcd deviates; residual standard error RSEω = 0.05). Simulations started with the visuomotor gains ωv, ωm and ωcd, the learning rates αv, αm and αcd, the peak velocity decay rate γκ and the duration compensation rate γλ of healthy subjects; αv, αm, αcd, γκ and γλ change until they match the values of cerebellar patients. Progression rates were νκ = 0.012, νλ = 0.040 and να = 1.7*10−7. Please note that the different subplots show different timescales. (B) Residual standard error RSEω depending on how fast learning rates (scaled by να), peak velocity decay rate (scaled by νκ) and duration compensation rate (scaled by νλ) change towards patients’ values. Simulations with RSEω ≤ 0.05 are marked with a white dot.

https://doi.org/10.1371/journal.pcbi.1011322.g008

We examined how fast the parameters can change towards patients’ values such that the simulation converges on the patients’ baseline visuomotor gains including outward localization of visual targets. Fig 8B depicts the residual standard error of the visuomotor gains at the end of the simulation (RSEω) depending on how fast learning rates, peak velocity decay rate and duration compensation rate change towards patients’ values (scaled by the parameters να, νκ and νλ). A white dot marks the simulations with RSEω ≤ 0.05 (reflecting a mean gain deviation below 4% at the end of the simulation). The simulation can converge on the patients’ baseline visuomotor gains, and, hence, reflect the visual outward localization that we found in the patients, if oculomotor fatigue compensation impairment occurs in the very early phase of the disease while learning is almost still intact, i.e. the learning rates decay slowly towards patients’ values. Only then, preserved learning can counteract the saccade vector decline.

The depicted timescale in Fig 8A would be much longer in a natural real life situation as saccades are usually performed to goals of different directions and eccentricities. Hence, our simulation is a simplification of the changes in oculomotor control during cerebellar disease progression and rather suggests relative timescale differences between the changes in peak velocity decay, duration compensation and learning rates. Moreover, the simulation remains speculative as long as we cannot compare the simulated disease progression to long-term data of cerebellar patients. However, we suggest that intertemporal oculomotor fatigue of multiple occurrence could lead to gradual changes of the visuomotor system similar to our simulation in the long run. Hence, the overestimation of target eccentricity in cerebellar patients may be a consequence of long-term perceptual learning to reduce motor errors.

Discussion

To dissociate cerebellar contributions to three stages of error-based learning, we modeled the visuomotor learning data of eight patients with a neurodegenerative cerebellar disease and eight healthy control subjects. The model, based on Masselink et al. (2021) [38], captures recalibration (1) of the visuospatial target map, (2) of an inverse model that transforms the visuospatial movement goal into a motor command, and (3) of a forward dynamics model that estimates the saccade size from corollary discharge. The model learns from postdictive motor error, i.e. motor error with respect to a postdictive update of the pre-saccadic target position. In patients, short-term learning was reduced in the inverse model and abolished in the visuospatial target map while the forward dynamics model seemed correctly informed about the reduced saccade change. Moreover, some types of cerebellar pathology impaired the ability to upregulate saccade duration in response to oculomotor fatigue. Modeling suggests that long-term error reduction at perceptual level may have counteracted, explaining the overestimation of target eccentricity in the patients’ baseline data. We conclude that the cerebellum recalibrates the inverse model on a short time-scale, particularly via control of saccade duration, while the forward dynamics model does not depend on cerebellar integrity. Long-term perceptual learning may partially compensate for cerebellar motor deficits.

Cerebellum implements short-term learning of the inverse model

Our results confirmed that the inverse model is adapted within the cerebellar microcircuit. Motor gain adaptation was largely restrained in cerebellar patients for both inward and outward errors. This means that the cerebellum is specifically involved in the part of adaptation that relies on credit assignment to a false computation of the motor command, i.e. irrespective of the adaptation of the visual pre-saccadic target representation. Similar impairments have been shown in lesioned monkeys [39, 40], in cerebellar patients [4143] and in healthy subjects with tDCS or magnetic stimulation of the cerebellum [45, 46, 68]. In the patients of Golla et al. (2008) [42], learning from outward target steps was completely suspended while in our study, a small saccade change was preserved. This difference may arise because the patients’ lesion of Golla et al. (2008) [42] was specifically located in the oculomotor vermis, i.e. the particular region assumed to process inverse model adaptation for oculomotor commands [56, 69, 70].

Cerebellar role in short-term learning of the visuospatial target map requires further examination

After visuomotor learning has initially been seen through a pure motor lens, numerous studies have revealed changes in the visuospatial perception of the movement goal [24, 25, 27]. Hence, motor learning is not simply a change of the motor command directed to the same goal location on every trial. Instead, learning underlies recalibration of both the visuospatial goal representation and the motor command required to bring the eyes to that location. This makes sense since motor errors cannot only be due to erroneous computation of the motor command but can also stem from a deficient goal representation. Control subjects showed a learning trend in the visual target representation, yet, in contrast to previous studies [2427, 38, 57] this change did not reach significance, probably due to a small sample size. In patients, a learning trend was not observable. Hence, evidence of differences between patients and control subjects is still lacking and further examination is required to evaluate whether the cerebellum is a neural substrate for visuospatial recalibration during motor learning. Neurophysiologically, the cerebellum is an ideal candidate to modulate visuospatial plasticity because of its inner circular organization [70, 71], its homogenous output organization [17, 18, 72], and its feedback projections to cortical retinotopic areas [1921].

The forward dynamics model is not affected by cerebellar pathology

Consistent with previous work [24, 25, 28, 38], in our study, healthy subjects’ localizations shifted pre- and post-saccadically while both changes remained absent in patients. Thus, CDV seemed to correctly reflect the saccade changes in healthy controls and in patients. This result has two major implications. First, it proposes that the efference copy that is routed into the forward dynamics model reflects the actual adapted motor command and not some kind of initial motor command prior to adaptive changes. Second, it suggests that also in patients, the forward dynamics model accurately captures the actual changes of the inverse model such that the saccade vector is neither over- nor underestimated. This means that the forward dynamics model involved in trans-saccadic spatial integration and feedback motor learning remains unaffected by cerebellar degeneration. Indeed, the cerebellum computes forward models for object motion [7375] and feedforward motor control (movement adjustments during execution [53, 62, 76]). However, the forward dynamics model measured in our task is one involved in trans-saccadic perception and feedback motor learning, i.e. from errors after movement completion. That this forward dynamics model may rely on areas upstream of the cerebellum is supported by studies showing that a disruption of the CD pathway from SC via MD thalamus to the FEF provokes a bias in visual space updating [7780]. The forward dynamics model could also be computed along other CD pathways, e.g. from SC via the thalamic pulvinar to parietal and occipital cortex [81, 82]. The role of parietal cortex in CD processing was recently demonstrated by Cheviet et al. (2022) [48]. At least, saccadic learning involves a distributed network across cerebral cortex [8385].

Cerebellar role in oculomotor fatigue compensation

Our results motivate further examination of whether cerebellar integrity may be crucial for counteracting oculomotor fatigue by saccade duration adjustment. While healthy subjects compensated peak velocity loss by 96% within a trial, patients achieved only 46%. As peak velocity loss and saccade duration compensation did not reach significance at the patients’ group level, the interpretation of these results need to be treated with caution. Still, the observed tendencies are in line with previous data of patients whose disease affected the oculomotor vermis. Firstly, Golla et al. (2008) [42] observed reduced oculomotor fatigue compensation, and, secondly, Xu-Wilson et al. (2009) [43] observed reduced within-saccade compensation for peak velocity fluctuations. Taken together, it seems worth to further examine whether saccade lengthening, within-saccade (due to oculomotor fatigue) or after saccade (due to post-saccadic motor error), is operated by cerebellar-dependent saccade duration control. That the cerebellum accommodates a velocity-duration tradeoff has been shown for hand pointing movements by Markanday et al. (2018) [86]. While movement duration was fine-tuned with respect to movement velocity in healthy controls, the velocity-duration tradeoff was completely absent in patients with different types of global cerebellar degeneration, i.e. comparable to our patient sample. It has to be noted that in our sample, the saccade fell short throughout the no step condition in five of eight patients with different types of cerebellar disease while only one control subject showed a comparable decline. Further examination of the role of the disease type could help to understand which cerebellar structures are critical to adjust saccade duration during oculomotor fatigue.

The cerebellum adjusts movement duration in response to velocity variations also for smooth pursuit [87]. Neurophysiologically, it remains to be examined how duration control is implemented. For manual movements, there is first evidence that the simple spikes of Purkinje cells also code movement velocity and duration [8890].

Long-term perceptual learning may compensate for cerebellar deficiencies in saccade duration control

In the long run, insufficient velocity compensation in cerebellar patients will lead to hypometric saccades. Indeed, saccade hypometry and reduced saccade velocity were observed in cerebellar patients [42, 66] and lesioned monkeys [39, 40, 64, 65]. Even if the affected part of the cerebellum in these studies was restricted to the oculomotor vermis, we also found a hypometric tendency in the course of repetitive saccades in some patients. Our model simulation for long-term behavior of a cerebellar-deficient visuomotor system suggests how saccade hypometry may occur and, with the gradual increase of motor errors, will be counteracted by learning at perceptual and motor level. The simulation can explain the significant overestimation of target eccentricity in the patients’ baseline data on the condition that the impairment in oculomotor fatigue compensation occurs in the very early phase of the disease while learning impairments proceed very slowly, and, hence, become apparent in the later phase of the disease.

Moreover, the simulation depicts a recovery from saccade hypometry as found in monkeys with cerebellar lesion [39, 40]. Yet, unlike the progressive nature of the cerebellar diseases in our patient sample, the cerebellar lesions in these studies occurred abruptly. In addition, recovery concerned saccade accuracy but not saccade precision, like in our patients who exhibited non-hypometric saccades despite high endpoint variability in the baseline. Hence, a tendency for saccade dysmetria due to deficiencies in saccade duration control may induce long-term perceptual and motor learning, probably upstream of the cerebellum. This may enable the visuomotor system to stay calibrated despite cerebellar short-term adaptation deficits.

Neurophysiologically, the visual overestimation of target eccentricity is interesting when put together with Zimmermann et al. (2015) [63] who found a similar result in a patient with a right posterior ventrolateral and ventromedial thalamic lesion. If, in both cases, this visual peculiarity is a consequence of motor adaptation deficits, these results underline the critical role of the cerebellar projection pathway via VL thalamus to frontal cortex for intact motor learning [47, 91].

Implementation of the postdictive motor error in the cerebellum

To model visuomotor behavior we used a postdictive error signal, i.e. the postdictve motor error, that drives learning of the pre-saccadic visual target position, the motor command and the CDV signal. Previous models on oculomotor learning a) have solely been focused on adaptation of the saccade vector and b) were usually driven by visual prediction error, thereby assuming that the saccade is always predicted to land with a stable target undershot as measured in the baseline [53, 92]. In Masselink et al. (2021) [38] we developed a new paradigm that allowed us to estimate the internal saccade size from pre- and trans-saccadic target localization with respect to the actually performed saccade, showing that this assumption cannot hold. While visual prediction error could not explain learning, learning was well explained by postdictive motor error. According to this, the visuomotor system computes a postdictive update of the pre-saccadic target position based on post-saccadic visual input and derives the error of the executed motor command with respect to this position. Supporting this idea, we recently found that adaptive changes to saccade targeting can be induced by post-saccadic visual feedback alone, i.e. when a pre-saccadic target is never shown [93]. It remains to be examined how such a postdictive error signal may be implemented in the cerebellum. The complex spikes of cerebellar Purkinje cells code visual and motor events, including the eccentricity of the post-saccadic target that is measurable 50 to 250 ms after saccade landing [7]. Post-saccadic target eccentricity could be merged with CDV information carried by the parallel fibers [4, 5, 94] to form the postdictive update of the target that is used to evaluate the motor command.

Limitations

Our model provides insights into how learning of visuomotor gains and oculomotor fatigue compensation is impaired in cerebellar patients. However, our model does not explain the increased within-subject and cross-subject saccade variability of the patients compared to the healthy control subjects. Increased within-subject variability of saccade endpoints is known from previous studies on cerebellar pathology [39, 40, 42, 43, 56]. The source of this phenomenon is not resolved yet and its effect on saccade motor learning remains to be examined. However, increased endpoint variability seems consistent with the role of the cerebellum in the within-saccade adjustment of movement duration in response to velocity changes that we found in the no step condition (see also Xu-Wilson et al., 2009 [43]). If the cerebellum cannot sufficiently adjust movement duration, i.e. on the basis of feedforward saccade control [43, 53, 62, 76], saccade endpoints will not only tend to be hypometric (in case of overall velocity decline across saccades) but will additionally be more variable as a result of trial-to-trial velocity fluctuations.

The increased cross-subject saccade variability of cerebellar patients may be due to the different types of cerebellar disease. However, we could not determine a relationship between disease type and saccadic behavior in our patient sample, e.g. by detecting systematic similarities of saccadic behavior between the three patients with spinocerebellar ataxia or between the two patients with Friedreich ataxia. As patients are difficult to recruit for study participation, the sample size was limited to eight patients and eight healthy control subjects. A larger sample size and a systematic analysis of similarities and differences in saccadic behavior is required to further examine how impairments of visuomotor learning and oculomotor fatigue compensation differ depending on the specific disease type and affected cerebellar areas.

In our model, we fit the pre-saccadic visual target position to visual localizations performed with a mouse cursor. The pre-saccadic visual target position defines the spatial saccade goal and is used as the input to the inverse model. The inverse model computes the muscle dynamics in order to bring the eyes to the target. Fitting the pre-saccadic visual target position to the visual localization data implies that the explicitly reported perception of the target is equivalent to the implicit visual representation that is used as a spatial saccade goal by the oculomotor system. The quantification of inverse model adaptation (that is derived from saccade changes relative to pre-saccadic localization changes) is sensitive to this assumption. However, consistent with this assumption, it seems plausible that pre-saccadic localization changes during learning [2327, 38, 57] reflect a change of the spatial saccade goal and not a sole change of visual perception without any effect on saccade execution. The reason lies within the nature of the task itself that serves a motor goal, i.e. to perform accurate movements, not to perform accurate perception. Learning is an optimization problem, and hence, should serve the optimization of the task goal. A sole change of target perception without any effect on saccade execution would not optimize the task goal, i.e. it would not optimize motor function, and, thus, would be useless. Hence, the presence of pre-saccadic localization changes during learning argue for a use by the motor system, i.e. as a spatial movement goal.

Current models of motor learning often include a decay rate that pulls behavior back to a baseline level, thereby introducing a trace of motor memory [92, 9598]. This allows to explain possible short-term forgetting of the learned behavior, e.g. during set breaks, that our model does not provide. In this approach, forgetting relies on the assumption that learned motor behavior is represented as a deviation from a baseline state of visuomotor transformations that is ingrained from the millions of saccades made in natural conditions and stays stable over the time of a typical adaptation experiment [99]. Our model differs from these models in that it is not based on a deviation from a steady state. In our model, the visuomotor system stays in a steady state that is reached, and can flexibly transit to a next steady state when another perturbation is applied without a pull-back process to a previous steady state. This still allows the model to return to baseline in case a perturbation is canceled because then, errors indicate that motor performance needs to be optimized by adaptation in the opposing direction. It would be interesting for future work to see how long-term learning and short-term adaptation can be combined in our model, perhaps by including multiple parallel learning processes at different time scales.

Conclusion

We suggest an expanded view on the cerebellar circuit function for oculomotor learning and fatigue compensation. The cerebellum does adapt the inverse model but the forward dynamics model for trans-saccadic space perception and feedback motor learning may be processed upstream of the cerebellum. Future research needs to examine the cerebellar role in visuospatial recalibration, and where the forward dynamics model is computed, e.g. in the thalamus, in the frontal eye fields or in other cerebral areas.

Methods

Subjects

Our study is based on the dataset of Cheviet et al. (2022) [48] which includes eight healthy control subjects (49.25 ± 7.67 years, four male, four female) and eight patients with a neurodegenerative disorder of the cerebellum (55.5 ± 9.61 years, four male, four female). Please see Table 1 listing the individual patients’ characteristics.

Each patient was clinically examined at the Neurological Hospital Pierre Wertheimer (Bron, France) to verify a chronic progressive cerebellar ataxia, normal or corrected-to-normal vision as well as the ability to concentrate and to remain seated for at least 30 min. Another inclusion criterion was that patients were able to maintain a stable hand position for at least 30 s. This ensured that visual target localizations performed with a mouse cursor were not affected by tremor. Patients with another neurological disease, unstable medical condition, psychotropic medication intake, pronounced nystagmus or ocular instability were excluded from study participation.

Ethics statement

All subjects gave written informed consent to participation. The experiment was performed in accordance with the ethical standards of the 1964 Helsinki Declaration and was approved by the ethics committee (Comité de Protection des Personnes Est-III; ID-RCB: 2017–00942-51; 17.05.09).

Setup

Subjects were seated in a completely dark room in front of a CRT monitor (19 inches, 1280 × 1024 pixels, 85 Hz, 57 cm viewing distance) that was covered with a neutral density filter to prevent visibility of monitor background light and contrast to the surrounding (ND4, 25% transmittance, 0.03 stimuli luminance, 0.01 room luminance).

Movements of the right eye were recorded by an Eyelink 1000 (remote configuration; SR Research, Ontario, Canada) at 1000 Hz with a five-point calibration before each session and drift correction after each session break. Online fixation control was based on a 5° position threshold (horizontally and vertically). Saccade onset was detected online by a 2.5° position threshold, a 22 velocity threshold and a 4000 acceleration threshold. With respect to saccade initiation (0 ms), in the cerebellar group, the saccade was detected online at 16.1 ± 2.3 ms, the target was stepped (in saccade trials) or disappeared (in post-saccadic localization trials) at 28.9 ± 2.3 ms while the saccade ended at 72.7 ± 23.2 ms. In the healthy group, from saccade initiation (0 ms), the saccade was detected online at 15.5 ± 1.6 ms, the target was stepped (or disappeared) at 28.3 ± 1.5 ms while the saccade ended at 64.1 ± 3.8 ms. The procedure was controlled with Experiment Builder (SR Research, Ontario, Canada).

Experimental design

Each subject participated in three experimental conditions recorded in randomized order with at least one week in between. Each condition measured reactive saccades to a 20° rightward target as well as visuospatial localizations of a flashed bar that varied horizontally around the target position by -4, -2, 0, +2 and +4° (Fig 1A). Localization was performed during fixation (pre-saccadic localization trials) and after saccade landing (post-saccadic localization trials).

The pre-exposure phase measured subjects’ baseline state, comprising four repetitions of five pre-saccadic localization trials and five post-saccadic localization trials, followed by a block of 20 saccade trials with peri-saccadic target offset (4*(5+5) + 20 = 60 trials; for an overview of the trial sequence see Fig 2A-B). The exposure phase measured 120 saccade trials with either a 6° peri-saccadic target step opposite to saccade direction (inward condition), or a 6° peri-saccadic target step in saccade direction (outward condition), or without any peri-saccadic target step (no step condition). The relatively large target step of 6° was chosen because cerebellar dysfunction is usually accompanied by a strong variability of saccade amplitudes, thereby leading to a strong variability of the post-saccadic visual target position [39, 40, 42, 43, 56]. In order to keep this disturbance as small as possible, we chose a rather large target step, and, consequently, a rather large pre-saccadic target eccentricity of 20° such the target step was 30% of the pre-saccadic target eccentricity. Control subjects and patients reported that they did not notice the target step. In contrast to the two learning conditions, the no step condition required subjects to maintain the saccade vector, which is usually achieved by upregulating saccade duration to counteract peak velocity decline in the course of repetitive, stereotyped saccades, known as oculomotor fatigue. A break was performed after the 50th and the 100th trial.

The post-exposure phase measured localizations after exposure (four repetitions of five pre-saccadic localization trials and five post-saccadic localization trials), interleaved with three blocks of each 20 saccade refresh trials (including the respective target step) to maintain the achieved saccade vector of the exposure phase (4*(5+5) + 3*20 = 100 trials; Fig 2A-B). Each session began with a 10 min dark adaptation phase to ensure that subjects were able to discriminate the stimuli with low luminance already at the start of the session. During session breaks (before each change of trial type and after the 50th and 100th trial of the exposure phase), the room was dimly illuminated to keep the level of dark adaptation constant. The Eyelink illuminator was not visible to the participants, even after dark adaptation because the 940-nm model was used.

All stimuli (target, flash and pointer) were presented with respect to the average horizontal gaze position of a 50 ms fixation reference period at trial start (fixation reference, see trial descriptions below). This procedure ensured consistent stimulus eccentricity despite reduced fixation accuracy of cerebellar patients. A trial was repeated if fixation criteria were violated. Inter-trial interval was a 800 ms black screen.

Saccade trials.

A red fixation circle of 1° diameter appeared for a random time interval between 800 and 1400 ms (5° to the right of the left screen border, vertically randomized between -1.5, 0 and 1.5° to the horizontal meridian, Fig 1A). If fixation criteria were met during this period, a red target circle was displayed 20° to the right with respect to the fixation reference determined during the last 50 ms of the fixation period. As soon as saccade onset was detected, the target switched off in the pre-exposure phase, and, in the exposure and post-exposure phase, stepped either 6° inward (inward condition), 6° outward (outward condition) or stayed at its position (no step condition). In the exposure and post-exposure phase, the target stayed visible for 100 ms.

Post-saccadic localization trials.

A red fixation circle of 1° diameter was presented 5° to the right of the left screen border on the horizontal meridian (Fig 1A). Subjects were instructed to fixate the circle and to indicate when ready to begin the trial via mouse press. If, after a 450 ms black screen, fixation criteria were met for 50 ms, a blue vertical line was flashed with counterbalanced horizontal distance of 16, 18, 20, 22 or 24° with respect to the fixation reference determined during the 50 ms fixation check period (12 ms, 0.2° width, 29° height). After a 400 ms black screen, a red target circle of 1° diameter appeared 20° to the right of the fixation reference. As soon as saccade onset was detected, the target was extinguished, followed by a 100 ms black screen. A blue line pointer was displayed on the right screen bottom (0.2° width, 10° height, 10° below the horizontal meridian, randomly 17, 19, 21 or 23° to the right of the fixation reference). As soon as subjects started moving the mouse, the pointer was blocked to the horizontal meridian. The subject had to move the pointer to the perceived horizontal flash position, confirming the judgement with a mouse press.

Please note that in case of any dysmetria of arm movements, patients had enough time to adjust the cursor position until it visually matched the position where they had perceived the flashed vertical line. As all patients were able to hold a stable position of the hand for at least 30 s, the mouse click could not be affected by tremor and should solely rely on the patients’ perceived position of the flashed line.

Pre-saccadic localization trials.

A red prohibited direction sign served as a fixation stimulus (with the same spatial parameters as in the post-saccadic localization trials, Fig 1A). The subject had to continuously fixate it from the first mouse press (to indicate she/he is ready to start the trial) to the second mouse press (to confirm the localization judgement). The flash was presented with the same latency and spatial parameters as in the post-saccadic localization trials. The pointer appeared 750 ms after flash offset to roughly match the respective interval of the post-saccadic localization trials, and with the same latency and spatial parameters as in the post-saccadic localization trials. As soon as the subject started moving the mouse, the pointer was blocked to the horizontal meridian and the subject had to move the pointer to the perceived horizontal flash position and press the mouse.

Training phase.

Before each session, subjects practiced a minimum of 5 pre- and 5 post-saccadic localization trials. In case of difficulties, online gaze position and a fixation box around the fixation circle were displayed to help subjects to familiarize with the task requirements.

Model

To quantify the adaptive changes in the visuomotor circuitry of patients and control subjects, we used the model of Masselink et al. (2021) [38] and expanded it with (1) how changes in the motor command are transposed to saccade peak velocity and saccade duration and (2) how the visuomotor system compensates for a fatigue-induced decline in saccade peak velocity (Fig 1B).

Model conception.

The model describes the saccadic circuitry based on three sensorimotor transformations, i.e. synaptic gains [14, 15]. These plastic gains are, first, a visual gain ωv to transform retinal input into target position on a visuospatial map, second, a motor gain ωm to transform spatial target position into a motor command (inverse model), and, third, a CD gain ωcd to transform the corollary discharge CDM of the motor command into CDV, i.e. the computed displacement of visual space due to the saccade (forward dynamics model).

The three gains learn to reduce postdictive motor error, an error of the motor command with respect to a postdictive update of the pre-saccadic target position. We use the postdictive motor error for learning based on the results of Masselink et al. (2021) [38]. Previous models on oculomotor learning have solely been focused on adaptation of the saccade vector, learning either from post-saccadic visual error or from visual prediction error [53, 62, 76, 100]. What argues against the post-saccadic visual error is that the visuomotor system accepts a remaining amount of post-saccadic visual error when learning reaches a steady state, i.e. when it reaches saturation [101103]. This means that post-saccadic visual error is not nullified at steady state. Moreover, the use of visual prediction error is usually based on the assumption that the saccade is predicted to land with a stable target undershoot as measured in the baseline. In Masselink et al. (2021) [38], we quantified the prediction, i.e. the internal visuospatial representation of saccade size CDV from pre- and trans-saccadic target localization with respect to the actually performed saccade. This revealed that visual prediction error is not nullified during learning. If the visuomotor system nullified visual prediction error, this would mean that after adaptation asymptote has been reached, the predicted visual error matches the actual visual error after saccade landing. Hence, subjects would predict the target to appear after the saccade (i.e. localize the pre-saccadic target after saccade landing in the dark) where the target actually appeared in saccade trials when a post-saccadic target is shown. This was obviously not the case. However, learning was well explained by nullification of postdictive motor error. This error relies on a postdictive update of the pre-saccadic target position based on post-saccadic visual input and computes the saccade error with respect to this position. Supporting this idea, in Heins et al. (2023) [93] we have recently shown that adaptive changes to oculomotor behavior can be induced by post-saccadic visual information alone, i.e. without pre-saccadic visual information.

Pre-saccadic calculations.

The pre-saccadic target position V1(n) on the visuospatial map is scaled by the visual gain ωv(n): (1) with the trial number n and physical target eccentricity P1 (Fig 1B-1). If ωv(n) = 1, the target is accurately localized on the visuospatial map. If a postdictive motor error is assigned to an internal failure of visuospatial target representation, ωv(n) is adapted to reduce the error for future movements (see Eq 14).

The inverse model (Fig 1B-2) transforms the pre-saccadic target position into a motor command with the motor gain ωm(n): (2)

Thus, the inverse model transforms the movement goal from visuospatial to motor coordinates based on the assumed state of muscle dynamics. If ωm(n) = 1, the saccade lands at the spatial position V1(n). If a postdictive motor error is assigned to a change in muscle dynamics, ωm(n) is adapted to reduce the error for future movements (see Eq 14). The quantification of changes to the visual gain ωv, on the one hand, and the motor gain ωm, on the other hand, dissociates the contribution of the two processes to behavioral output, i.e. saccade vector change.

The motor command M(n) is copied into CDM(n): (3) and directed into the CD pathway.

Before saccade start, the forward dynamics model (Fig 1B-3) transforms the corollary discharge of the motor command CDM(n) into the CDV(n) signal, i.e. the computed displacement of visual space, to inform the visual system about the upcoming saccade: (4)

If ωcd(n) = 1, CDV(n) matches the actual saccade vector. If a postdictive motor error is assigned to an internal failure of CDV representation, ωcd(n) is adapted to reduce the error for future movements (see Eq 14).

Based on CDV(n), the forward outcome model (Fig 1B-4) shifts the coordinates of the visual pre-saccadic target position V1(n) to predict the visual post-saccadic target position: (5)

If CDV(n) is correct, i.e. if it matches the actual saccade vector, pre- and post-saccadic visual scene are correctly integrated such that the post-saccadic target is predicted to appear () at the same spatiotopic position where the pre-saccadic target was perceived (V1). In the experiment, this is the case when the flash is localized at the same position in the pre- and post-saccadic localization trials.

Saccade kinematics.

Saccadic motor learning has been shown to be differently transposed to saccade kinematics, i.e. to saccade peak velocity and saccade duration, depending on the error direction [42, 56]. Moreover, there is evidence that cerebellar patients are more impaired in learning from outward than from inward errors [42]. Thus, capturing the learning effect on the patients’ saccade kinematics may provide a deeper understanding of how the cerebellum operates the recalibration of motor commands. Hence, we expanded the model of Masselink et al. (2021) [38] with how the motor command is transposed to saccade peak velocity and saccade duration. Accordingly, the saccade is executed with the saccade peak velocity κ(n) and the saccade duration λ(n) that are specific for the motor command (Fig 1B-5, see Eqs 16, 18 and 20 for how the motor command is transposed to saccade peak velocity and saccade duration depending on the paradigm). This specificity relies on a linear relationship between a saccade’s amplitude, its peak velocity and its duration, known as the saccadic main sequence [51, 104, 105]. Hence, we specify the relationship between the motor command M(n), the saccade peak velocity κ(n) and the saccade duration λ(n) by the following equation: (6)

For the modeling, we derive β0, βκ and βλ from fitting the equation to the saccade vectors, saccade peak velocities and saccade durations of patients and healthy controls. Please note that while the saccadic main sequence is usually expressed by two functions describing the relationship of the saccade amplitude with saccade peak velocity and saccade duration separately [51, 103], we unite this relationship in one single function. This allows us to model oculomotor fatigue compensation, i.e. when a decrease of saccade peak velocity does not inevitably result in a decrease of the saccade amplitude but can be compensated by a longer saccade duration such that the saccade amplitude stays stable. Two separate functions expressing the relationship of the saccade amplitude, first, to peak velocity and, second, to duration would not allow to model oculomotor fatigue compensation. We used a linear approximation in sort of a 3D plane as described by Eq 6 that appeared well justified over the rather small range of saccade amplitudes that we measured (Figs 3 and 4).

Saccade execution.

The execution of the motor command produces the physical saccade vector (Fig 1B-6): (7) with random motor noise ϵm(n).

In the inward and outward learning paradigms, the target is shifted during saccade execution. The peri-saccadic target displacement due to the target shift Ps and the motor execution noise ϵm is described by: (8)

Post-saccadic calculations.

After saccade offset, the post-saccadic target is localized on the visuospatial map (Fig 1B-7): (9)

Please note that the visual gain ωv(n) is not applied to V2(n) because visuospatial changes during saccadic motor learning underlie an adaptation field [24, 25]. Hence, a change in the visual gain ωv(n) is a local effect around the pre-saccadic target position and should not be applied to the representation of the post-saccadic target after saccade landing.

The backward outcome model integrates the visual post-saccadic target position V2(n) with CDV(n) to from a postdictive update of the target position in a pre-saccadic reference frame (Fig 1B-8): (10)

We refer to this process as the backward outcome model as it performs a postdictive transformation from post- to pre-saccadic space, analogously to the forward outcome model that performs a predictive transformation, i.e. from pre- to post-saccadic space.

Learning of visuomotor gains and its transposition to saccade kinematics.

In the inward and outward learning paradigms, the visuomotor system experiences a significant motor error due to the peri-saccadic target step. This error, termed the postdictive motor error Epost(n), is calculated as the error of the motor command with respect to the postdicted pre-saccadic target position (Fig 1A-9): (11)

To reduce Epost(n) on the next trial, the visuomotor system adapts its gains following an internal estimate of the Epost(n) gradient [15, 94, 106]. We express this principle by the delta rule [107, 108]: (12) (13) (14)

The learning rates 0 specify the speed of learning.

The changes in the visual gain ωv and in the motor gain ωm result in a changed motor command M(n+1) on the next trial. The transposition of the changed motor command to saccade peak velocity and saccade duration needs to take into account whether the motor command was decreased or increased, i.e. whether changes occurred during inward or during outward learning [42, 56]. Saccade shortening during inward learning is usually characterized by a decline of peak velocity while saccade duration stays constant. Accordingly, during inward learning in our model (Ps < 0), saccade duration λ stays constant while saccade peak velocity κ is downregulated on the next trial according to the changed motor command (Fig 1B-6): (15) (16)

By contrast, saccade lengthening during outward learning is usually characterized by an increase of saccade duration while saccade peak velocity stays constant [42, 56]. Accordingly, during outward learning in our model (Ps > 0), saccade peak velocity κ stays constant while saccade duration λ is upregulated on the next trial according to the changed motor command (Fig 1B-6): (17) (18)

The visuomotor system reaches a steady state, i.e. learning comes to an end, if the saccade lands at the postdicted pre-saccadic target position such that Epost(n) = 0.

Oculomotor fatigue and its compensation by saccade duration.

For a deeper understanding of which calibration processes are perturbed due to cerebellar disease, we also modeled the experimental condition without any target step, i.e. without learning. Although target localizations and saccade vectors should stay roughly constant in healthy subjects, the repetitive execution of saccades to the same target position is usually accompanied by a decline in saccade peak velocity, a sign of oculomotor fatigue. However, this compensating mechanism has been shown to be impaired in patients with vermal cerebellar lesion [42]. Accordingly, in the no step condition of our model (Ps = 0), saccade peak velocity κ declines while saccade duration λ compensates by a certain percentage on the next trial (Fig 1B-6): (19) (20)

Here, γκ ≥ 0 describes the decay rate of peak velocity, κend limits peak velocity loss and γλ describes the percentage by which saccade duration compensates for the peak velocity loss (0 ≤ γλ ≤ 1). According to the changes of peak velocity κ and duration λ in the no step condition (Ps = 0), the motor gain ωm(n+1) of the next trial is: (21)

Hence, during expected oculomotor fatigue in the no step condition, we model changes in saccade peak velocity and saccade duration that are reflected in the motor gain ωm while ωv and ωcd are assumed to stay stable.

Analysis

Data processing.

Data analysis was performed in Matlab R2017a (Mathworks, Natick, MA). Saccades were detected offline with a self-made routine that searched for a saccade after target onset with an onset and offset velocity threshold of 80 deg/s. Saccade traces as well as onset and offset detection were visualized trial by trial in order to correct for detection errors. Based on the verified saccade onset and offset, saccade latency, peak velocity and duration were extracted. The horizontal eye position signal was smoothed and the saccade start position was determined based on the position 50 ms prior to saccade onset detection. The saccade end position was extracted based on the position 50 ms after saccade offset detection to prevent post-saccadic oscillation effects.

Saccades with a latency of less than 100 ms, an amplitude below 10° or above 30° as well as localizations with more than ±7° localization error were excluded from the analysis. Moreover, as our modeling approach explicitly relies on the localization judgement of target eccentricity, we chose only the pre- and post-saccadic localization trials in which the flash position matched the saccade target position, i.e. 20° eccentricity. Thus, we ensured that our modeling approach quantifies perceptual changes exactly at the spatial position of the target and is not disturbed by the process of learning transfer to nearby locations, known as the saccadic adaptation field [12, 24, 25, 52, 60, 69, 109, 110].

For each condition of each subject, we calculated the mean saccade vector, mean saccade peak velocity and mean saccade duration of the 20 pre-exposure saccade trials and the last 20 exposure saccade trials, respectively. The removal of invalid trials resulted in 39.3 ± 2.7 used saccade trials in control subjects and 33.2 ± 3.8 used saccade trials in patients (summed across the pre-exposure and late exposure phase averaged across the three conditions; difference between control subjects and patients t14 = 1.97, p = .069).

For the pre- and the post-saccadic flash localizations we extracted the mean of the pre-exposure and the post-exposure phase. These were eight trials as we did only include those in which the flash matched the distance of the target, i.e. 20° eccentricity. The removal of invalid trials resulted in 7.8 ± 0.5 used pre-saccadic localization trials in control subjects and 7.5 ± 0.8 used pre-saccadic localization trials in patients (difference between control subjects and patients W = 75.5, p = .446). From the post-saccadic localization trials, 7.5 ± 0.2 trials were used in the control group and 7.5 ± 0.5 trials were used in the patients group (difference between control subjects and patients W = 64.5, p = .668).

In addition, we quantified a cross-condition baseline state for each subject, averaging the mean saccade vectors, saccade peak velocities, saccade durations and pre- and post-saccadic localizations across the three conditions. Outliers with more than three standard deviations from the mean were excluded. We averaged the saccade vectors, pre- and post-saccadic localizations at group level, i.e. separately for patients and control subjects and condition. Moreover, we calculated the mean trial-by-trial saccade vectors, saccade peak velocities and saccade durations of the 120 exposure trials, separately for patients and control subjects. After invalid trials were removed, 117.3 ± 11.1 saccade trials could be used in the control group and 96.9 ± 11.1 saccade trials could be used in the patient group (difference between control subjects and patients, Wilcoxon rank sum test, W = 86.0, p = .065).

Model-based analysis.

We calculated the state of the CDV signal, the visuomotor gains ω and the Epost error of each subject’s pre-exposure and post-exposure phase. Based on the pre-saccadic target localization (V1), the saccade vector (M) and the post-saccadic target localization with respect to the saccade landing position (), we derive: (22) (23) (24) (25) (26) with P1 = 20° and Pd = 0° for the pre-exposure phase and the post-exposure phase of the no step condition, Pd = -6° for the post-exposure phase of the inward condition and Pd = 6° for the post-exposure phase of the outward condition. Additionally, we derived the change of each variable from the pre- to the post-exposure phase.

Before fitting the model to the trial-by-trial data, we performed a regression of the saccade vectors on saccade peak velocities and saccade durations of the exposure phase, separately for each group (controls, patients) and each condition (inward, no step, outward). See Table 2 for the fitted regression weights.

Starting from the pre-exposure mean of pre-saccadic localization (V1(1)), post-saccadic localization (with respect to the saccade landing point, ), saccade vector (M(1)), saccade peak velocity (κ(1)) and saccade duration (λ(1)) taken as values of trial n = 1, we fitted the model to the exposure phase, separately for each group (control subjects, patients) and each condition (inward, no step, outward). For the inward and the outward condition, we fitted the learning rates of the visuomotor gains (ϕfit = α), and for the no step condition, we fitted the peak velocity decay rate and the duration compensation rate (ϕfit = (γκ, γλ)). The fitting procedure minimized the weighted sum of squared errors (SSE): with the weights = , ηM = 0.015, ηκ = 0.012 and ηλ = 0.003 to account for the unequal number of data points the signals were fitted to (e.g. only two data points for V1 and but 120 data points for M, κ and λ) and for the unequal units (e.g. a saccade vector error of 3° should be weighted more than a peak velocity error of 3). In the inward and in the outward condition patients showed a small V1 change in the opposite direction of learning. This effect was small, not significant and, thus, rather of random nature (see Fig 5B). In these two cases, we fitted a stable V1 as the mean of the pre-saccadic target localizations between the pre- and the post-exposure phase. For all model fits, we set ϵm(n) = 0.

For goodness of fit, residual standard errors RSE were calculated, collectively for V1, M and (in °) and separately for κ (in ) and λ (in ms): (27) with q as the number of data points used for the respective RSE calculation.

Simulation of long-term visuomotor behavior after disease onset.

We simulated how the visuomotor system would behave on a long timescale after the onset of cerebellar disease, i.e. when a fatigue-induced loss of saccade peak velocity κ can be only partially compensated by upregulation of saccade duration λ while visuomotor learning is progressively reduced. Simulations were performed for a target of P1 = 10° eccentricity, ω(1) = (0.958; 1.023; 0.980) (which are the baseline gains measured in healthy subjects of Masselink et al. (2021) [38]), κ(1) = 300, λ(1) = 43.9 ms, γκ(1) = 0.003, γλ(1) = 0.964 (derived from the control subjects’ model fit to the no step condition), κend = 200, αv(1) = 2.35*10−6, αm(1) = 5.53*10−6, αcd(1) = 5.31*10−13 (derived from control subjects’ model fit to the outward condition), β0 = -11.67, βκ = 0.03, βκ = 0.27 (derived from the saccade vectors, saccade peak velocities and saccade durations of control subjects across all three conditions to capture a minimum range of saccade vectors). The progressive impairment of oculomotor control parameters was reflected by a successive change of the peak velocity decay rate γκ, the duration compensation rate γλ and the learning rates α towards patients’ values: (28) (29) (30)

The progressive change proceeded until (1) the oculomotor fatigue parameters matched those obtained from the patients’ model fit to the no step condition (γκ,end and γλ,end), and (2) until the learning rates matched those obtained from the patients’ model fit to the outward condition (αend). The progression rates νκ, νλ and να determine how fast the respective parameters change from healthy controls’ to patients’ values. We examined which range of progression rates can mirror the course of the disease such that the oculomotor system converges at the patients’ baseline gains ω, hence, leading to an outward localization of visual targets. Thereby, we tested which parameter combinations lead to a residual standard error RSEω ≤ 0.05 at the end of the simulation with respect to the patients’ baseline gains ω. This corresponds to a mean gain deviation below 4% from the patients’ baseline gains at the end of the simulation. For each simulation, we computed the progression of the disease across 23,328 Mio saccades. With 3 saccades per second and 12 awake hours per day this roughly corresponds to the number of saccades that are performed across 6 months. Though, please note that 6 months in our simulation would correspond to a multiple of 6 months in real life as in real life we perform saccades of different amplitudes and different directions. Hence, our simulation is a strong simplification of the changes in oculomotor control during cerebellar disease progression and rather examines relative timescale differences between the changes in peak velocity decay, duration compensation and learning rates.

Statistical analysis.

To compare the data of one group against zero, we used two-sided one-sample t-tests or a Wilcoxon signed rank test in case normality distribution was violated. To compare data between groups, i.e. between patients and control subjects, two-sided two-sample t-tests, or, alternatively, Wilcoxon rank sum tests in case normality distribution was violated, were calculated. All tests were conducted with a significance level of 0.05.

Acknowledgments

We thank Eric Koun and Roméo Salemme for their assistance with the technical implementation of the protocol and their advice about eye movements analysis tools.

References

  1. 1. Marr D. A theory of cerebellar cortex. J Physiol. 1969;202:437–470. pmid:5784296
  2. 2. Albus JS. A theory of cerebellar function. Math Biosci. 1971;10(1-2):25–61.
  3. 3. Herzfeld DJ, Kojima Y, Soetedjo R, Shadmehr R. Encoding of action by the Purkinje cells of the cerebellum. Nature. 2015;526(7573):439–442. pmid:26469054
  4. 4. Bays PM, Wolpert DM. Computational principles of sensorimotor control that minimize uncertainty and variability. J Physiol. 2007;578(2):387–396. pmid:17008369
  5. 5. Herzfeld DJ, Kojima Y, Soetedjo R, Shadmehr R. Encoding of error and learning to correct that error by the Purkinje cells of the cerebellum. Nat Neurosci. 2018;21(5):736–743. pmid:29662213
  6. 6. Raymond JL, Medina JF. Computational principles of supervised learning in the cerebellum. Annu Rev Neurosci. 2018;41:233–253. pmid:29986160
  7. 7. Markanday A, Inoue J, Dicke PW, Thier P. Cerebellar complex spikes multiplex complementary behavioral information. PLoS Biol. 2021;19(9):e3001400. pmid:34529650
  8. 8. Kawato M. Internal models for motor control and trajectory planning. Curr Opin Neurobiol. 1999;9(6):718–727. pmid:10607637
  9. 9. Franklin DW, Wolpert DM. Computational mechanisms of sensorimotor control. Neuron. 2011;72(3):425–442. pmid:22078503
  10. 10. McLaughlin S. Parametric adjustment in saccadic eye movements. Percept Psychophys. 1967;2(8):359–362.
  11. 11. Havermann K, Lappe M. The influence of the consistency of postsaccadic visual errors on saccadic adaptation. J Neurophysiol. 2010;103(6):3302–3310. pmid:20393067
  12. 12. Pelisson D, Alahyane N, Panouilleres M, Tilikete C. Sensorimotor adaptation of saccadic eye movements. Neurosci Biobehav Rev. 2010;34(8):1103–1120. pmid:20026351
  13. 13. Srimal R, Diedrichsen J, Ryklin EB, Curtis CE. Obligatory adaptation of saccade gains. J Neurophysiol. 2008;99(3):1554–1558. pmid:18234985
  14. 14. Shadmehr R, Smith MA, Krakauer JW. Error correction, sensory prediction, and adaptation in motor control. Annu Rev Neurosci. 2010;33:89–108. pmid:20367317
  15. 15. Wolpert DM, Diedrichsen J, Flanagan JR. Principles of sensorimotor learning. Nature Rev Neurosci. 2011;12(12):739–751. pmid:22033537
  16. 16. Middleton FA, Strick PL. Cerebellar projections to the prefrontal cortex of the primate. J Neurosci. 2001;21(2):700–712. pmid:11160449
  17. 17. Hoshi E, Tremblay L, Féger J, Carras PL, Strick PL. The cerebellum communicates with the basal ganglia. Nat Neurosci. 2005;8(11):1491–1493. pmid:16205719
  18. 18. Bostan AC, Dum RP, Strick PL. The basal ganglia communicate with the cerebellum. PNAS. 2010;107(18):8452–8456. pmid:20404184
  19. 19. Buckner RL, Krienen FM, Castellanos A, Diaz JC, Yeo BT. The organization of the human cerebellum estimated by intrinsic functional connectivity. J Neurophysiol. 2011;106(5):2322–2345. pmid:21795627
  20. 20. Brissenden JA, Levin EJ, Osher DE, Halko MA, Somers DC. Functional evidence for a cerebellar node of the dorsal attention network. J Neurosci. 2016;36(22):6083–6096. pmid:27251628
  21. 21. Marek S, Siegel JS, Gordon EM, Raut RV, Gratton C, Newbold DJ, et al. Spatial and temporal organization of the individual human cerebellum. Neuron. 2018;100(4):977–993. pmid:30473014
  22. 22. van Es DM, van der Zwaag W, Knapen T. Topographic maps of visual space in the human cerebellum. Curr Biol. 2019;29(10):1689–1694. pmid:31080082
  23. 23. Moidell BG, Bedell HE. Changes in oculocentric visual direction induced by the recalibration of saccades. Vision Res. 1988;28(2):329–336. pmid:3414020
  24. 24. Collins T, Dore-Mazars K, Lappe M. Motor space structures perceptual space: Evidence from human saccadic adaptation. Brain Res. 2007;1172:32–39. pmid:17803970
  25. 25. Schnier F, Zimmermann E, Lappe M. Adaptation and mislocalization fields for saccadic outward adaptation in humans. J Eye Mov Res. 2010;3(4):1–18.
  26. 26. Zimmermann E, Lappe M. Motor signals in visual localization. J Vis. 2010;10(6):1–11. pmid:20884551
  27. 27. Gremmler S, Bosco A, Fattori P, Lappe M. Saccadic adaptation shapes visual space in macaques. J Neurophysiol. 2014;111(9):1846–1851. pmid:24523523
  28. 28. Bahcall DO, Kowler E. Illusory shifts in visual direction accompany adaptation of saccadic eye movements. Nature. 1999;400(6747):864–866. pmid:10476963
  29. 29. Zimmermann E, Lappe M. Mislocalization of flashed and stationary visual stimuli after adaptation of reactive and scanning saccades. J Neurosci. 2009;29(35):11055–11064. pmid:19726664
  30. 30. Duhamel JR, Colby CL, Goldberg ME. The updating of the representation of visual space in parietal cortex by intended eye movements. Science. 1992;255:90–92. pmid:1553535
  31. 31. Umeno MM, Goldberg ME. Spatial processing in the monkey frontal eye field. I. Predictive visual responses. J Neurophysiol. 1997;78(3):1373–1383. pmid:9310428
  32. 32. Crapse TB, Sommer MA. Frontal eye field neurons with spatial representations predicted by their subcortical input. J Neurosci. 2009;29(16):5308–5318. pmid:19386927
  33. 33. Sommer MA, Wurtz RH. Brain circuits for the internal monitoring of movements. Annu Rev Neurosci. 2008;31:317–338. pmid:18558858
  34. 34. Crapse TB, Sommer MA. Corollary discharge circuits in the primate brain. Curr Opin Neurobiol. 2008;18(6):552–557. pmid:18848626
  35. 35. Sperry RW. Neural basis of the spontaneous optokinetic response produced by visual inversion. J Comp Physiol Psychol. 1950;43(6):482–489. pmid:14794830
  36. 36. von Holst E, Mittelstaedt H. Das Reafferenzprinzip (Wechselwirkung zwischen Zentralnervensystem und Peripherie). Naturwissenschaften. 1950;37:464–476.
  37. 37. Wurtz RH. Neuronal mechanisms of visual stability. Vision Res. 2008;48(20):2070–2089. pmid:18513781
  38. 38. Masselink J, Lappe M. Visuomotor learning from postdictive motor error. eLife. 2021;10:e64278. pmid:33687328
  39. 39. Takagi M, Zee DS, Tamargo RJ. Effects of lesions of the oculomotor vermis on eye movements in primate: saccades. J Neurophysiol. 1998;80(4):1911–1931. pmid:9772249
  40. 40. Barash S, Melikyan A, Sivakov A, Zhang M, Glickstein M, Thier P. Saccadic dysmetria and adaptation after lesions of the cerebellar cortex. J Neurosci. 1999;19(24):10931–10939. pmid:10594074
  41. 41. Alahyane N, Fonteille V, Urquizar C, Salemme R, Nighoghossian N, Pelisson D, et al. Separate neural substrates in the human cerebellum for sensory-motor adaptation of reactive and of scanning voluntary saccades. The Cerebellum. 2008;7(4):595–601. pmid:19009327
  42. 42. Golla H, Tziridis K, Haarmeier T, Catz N, Barash S, Thier P. Reduced saccadic resilience and impaired saccadic adaptation due to cerebellar disease. Eur J Neurosci. 2008;27(1):132–144. pmid:18184318
  43. 43. Xu-Wilson M, Chen-Harris H, Zee DS, Shadmehr R. Cerebellar contributions to adaptive control of saccades in humans. J Neurosci. 2009;29(41):12930–12939. pmid:19828807
  44. 44. Panouilleres M, Neggers SF, Gutteling TP, Salemme R, Stigchel Svd, van der Geest JN, et al. Transcranial magnetic stimulation and motor plasticity in human lateral cerebellum: dual effect on saccadic adaptation. Hum Brain Mapp. 2012;33(7):1512–1525. pmid:21692144
  45. 45. Avila E, van der Geest JN, Kengne Kamga S, Verhage MC, Donchin O, Frens MA. Cerebellar transcranial direct current stimulation effects on saccade adaptation. Neural Plast. 2015;2015. pmid:25821604
  46. 46. Jenkinson N, Miall RC. Disruption of saccadic adaptation with repetitive transcranial magnetic stimulation of the posterior cerebellum in humans. The Cerebellum. 2010;9(4):548–555. pmid:20665254
  47. 47. Middleton FA, Strick PL. Basal ganglia and cerebellar loops: motor and cognitive circuits. Brain Res Rev. 2000;31(2-3):236–250. pmid:10719151
  48. 48. Cheviet A, Masselink J, Koun E, Salemme R, Lappe M, Froment-Tilikete C, et al. Cerebellar Signals Drive Motor Adjustments and Visual Perceptual Changes during Forward and Backward Adaptation of Reactive Saccades. Cereb Cortex. 2022;.
  49. 49. Schmidt D, Abel L, DellOsso L, Daroff R. Saccadic velocity characteristics- Intrinsic variability and fatigue. ASEM. 1979;50(4):393–395. pmid:464963
  50. 50. Fuchs AF, Binder MD. Fatigue resistance of human extraocular muscles. J Neurophysiol. 1983;49(1):28–34. pmid:6827300
  51. 51. Bahill AT, Clark MR, Stark L. The main sequence, a tool for studying human eye movements. Math Biosci. 1975;24(3-4):191–204.
  52. 52. Straube A, Robinson FR, Fuchs AF. Decrease in saccadic performance after many visually guided saccadic eye movements in monkeys. Invest Ophthalmol Vis Sci. 1997;38(13):2810–2816. pmid:9418734
  53. 53. Chen-Harris H, Joiner WM, Ethier V, Zee DS, Shadmehr R. Adaptive control of saccades via internal feedback. J Neurosci. 2008;28(11):2804–2813. pmid:18337410
  54. 54. Xu-Wilson M, Zee DS, Shadmehr R. The intrinsic value of visual information affects saccade velocities. Exp Brain Res. 2009;196(4):475–481. pmid:19526358
  55. 55. Prsa M, Dicke PW, Thier P. The absence of eye muscle fatigue indicates that the nervous system compensates for non-motor disturbances of oculomotor function. J Neurosci. 2010;30(47):15834–15842. pmid:21106822
  56. 56. Thier P, Markanday A. Role of the vermal cerebellum in visually guided eye movements and visual motion perception. Annu Rev Vis. 2019;5:247–268. pmid:31299168
  57. 57. Hernandez TD, Levitan CA, Banks MS, Schor CM. How does saccade adaptation affect visual perception? J Vis. 2008;8(3):1–16. pmid:18831626
  58. 58. Kojima Y, Iwamoto Y, Yoshida K. Memory of learning facilitates saccadic adaptation in the monkey. J Neurosci. 2004;24(34):7531–7539. pmid:15329400
  59. 59. Panouilleres M, Weiss T, Urquizar C, Salemme R, Munoz DP, Pelisson D. Behavioural evidence of separate adaptation mechanisms controlling saccade amplitude lengthening and shortening. J Neurophysiol. 2008;101(3):1550–1559. pmid:19091922
  60. 60. Miller JM, Anstis T, Templeton WB. Saccadic plasticity: Parametric adaptive control by retinal feedback. J Exp Psychol Hum Percept Perform. 1981;7(2):356–366. pmid:6453929
  61. 61. Wallman J, Fuchs AF. Saccadic gain modification: visual error drives motor adaptation. J Neurophysiol. 1998;80(5):2405–2416. pmid:9819252
  62. 62. Ethier V, Zee DS, Shadmehr R. Changes in control of saccades during gain adaptation. J Neurosci. 2008;28(51):13929–13937. pmid:19091981
  63. 63. Zimmermann E, Ostendorf F, Ploner CJ, Lappe M. Impairment of saccade adaptation in a patient with a focal thalamic lesion. J Neurophysiol. 2015;113(7):2351–2359. pmid:25652924
  64. 64. Ignashchenkova A, Dash S, Dicke PW, Haarmeier T, Glickstein M, Thier P. Normal spatial attention but impaired saccades and visual motion perception after lesions of the monkey cerebellum. J Neurophysiol. 2009;102(6):3156–3168. pmid:19759322
  65. 65. Ohki M, Kitazawa H, Hiramatsu T, Kaga K, Kitamura T, Yamada J, et al. Role of primate cerebellar hemisphere in voluntary eye movement control revealed by lesion effects. J Neurophysiol. 2009;101(2):934–947. pmid:19196922
  66. 66. Vahedi K, Rivaud S, Amarenco P, Pierrot-Deseilligny C. Horizontal eye movement disorders after posterior vermis infarctions. J Neurol Neurosurg Psychiatry. 1995;58(1):91–94. pmid:7823077
  67. 67. Tseng Yw, Diedrichsen J, Krakauer JW, Shadmehr R, Bastian AJ. Sensory prediction errors drive cerebellum-dependent adaptation of reaching. J Neurophysiol. 2007;98(1):54–62. pmid:17507504
  68. 68. Panouillères MT, Miall RC, Jenkinson N. The role of the posterior cerebellum in saccadic adaptation: a transcranial direct current stimulation study. J Neurosci. 2015;35(14):5471–5479. pmid:25855165
  69. 69. Hopp JJ, Fuchs AF. The characteristics and neuronal substrate of saccadic eye movement plasticity. Prog Neurobiol. 2004;72(1):27–53. pmid:15019175
  70. 70. Dean P, Porrill J, Ekerot CF, Jörntell H. The cerebellar microcircuit as an adaptive filter: experimental and computational evidence. Nat Rev Neurosci. 2010;11(1):30–43. pmid:19997115
  71. 71. Tanaka H, Ishikawa T, Lee J, Kakei S. The cerebro-cerebellum as a locus of forward model: a review. Front Syst Neurosci. 2020;14:19. pmid:32327978
  72. 72. Sultan F, Mock M, Thier P. Functional architecture of the cerebellar system. Neurol Dis Therapy. 2000;50:1–52.
  73. 73. Nowak DA, Hermsdörfer J, Rost K, Timmann D, Topka H. Predictive and reactive finger force control during catching in cerebellar degeneration. The Cerebellum. 2004;3(4):227–235. pmid:15686101
  74. 74. Miall RC, Christensen LOD, Cain O, Stanley J. Disruption of state estimation in the human lateral cerebellum. PLoS Biol. 2007;5(11):e316. pmid:18044990
  75. 75. Roth MJ, Synofzik M, Lindner A. The cerebellum optimizes perceptual predictions about external sensory events. Curr Biol. 2013;23(10):930–935. pmid:23664970
  76. 76. Albert ST, Shadmehr R. Estimating properties of the fast and slow adaptive processes during sensorimotor adaptation. J Neurophysiol. 2018;119(4):1367–1393. pmid:29187548
  77. 77. Ostendorf F, Liebermann D, Ploner CJ. Human thalamus contributes to perceptual stability across eye movements. PNAS. 2010;107(3):1229–1234. pmid:20080657
  78. 78. Prime SL, Vesia M, Crawford JD. TMS over human frontal eye fields disrupts trans-saccadic memory of multiple objects. Cereb Cortex. 2010;20(4):759–772. pmid:19641017
  79. 79. Ostendorf F, Kilias J, Ploner CJ. Theta-burst stimulation over human frontal cortex distorts perceptual stability across eye movements. Cereb Cortex. 2012;22(4):800–810. pmid:21705393
  80. 80. Cavanaugh J, Berman RA, Joiner WM, Wurtz RH. Saccadic corollary discharge underlies stable visual perception. J Neurosci. 2016;36(1):31–42. pmid:26740647
  81. 81. Wurtz RH, Joiner WM, Berman RA. Neuronal mechanisms for visual stability: progress and problems. Philos Trans R Soc B: Biol Sci. 2011;366(1564):492–503. pmid:21242138
  82. 82. Berman RA, Cavanaugh J, McAlonan K, Wurtz RH. A circuit for saccadic suppression in the primate brain. J Neurophysiol. 2017;117(4):1720–1735. pmid:28003409
  83. 83. Blurton SP, Raabe M, Greenlee MW. Differential cortical activation during saccadic adaptation. J Neurophysiol. 2012;107(6):1738–1747. pmid:22170969
  84. 84. Gerardin P, Miquée A, Urquizar C, Pélisson D. Functional activation of the cerebral cortex related to sensorimotor adaptation of reactive and voluntary saccades. Neuroimage. 2012;61(4):1100–1112. pmid:22465298
  85. 85. Guillaume A, Fuller JR, Srimal R, Curtis CE. Cortico-cerebellar network involved in saccade adaptation. J Neurophysiol. 2018;120(5):2583–2594. pmid:30207858
  86. 86. Markanday A, Messner J, Thier P. A loss of a velocity-duration trade-off impairs movement precision in patients with cerebellar degeneration. Eur J Neurosci. 2018;48(4):1976–1989. pmid:29972715
  87. 87. Dash S, Thier P. Smooth pursuit adaptation (SPA) exhibits features useful to compensate changes in the properties of the smooth pursuit eye movement system due to usage. Front Syst Neurosci. 2013;7:67. pmid:24146639
  88. 88. Roitman AV, Pasalar S, Johnson MT, Ebner TJ. Position, direction of movement, and speed tuning of cerebellar Purkinje cells during circular manual tracking in monkey. J Neurosci. 2005;25(40):9244–9257. pmid:16207884
  89. 89. Pasalar S, Roitman A, Durfee W, Ebner T. Force field effects on cerebellar Purkinje cell discharge with implications for internal models. Nat Neurosci. 2006;9(11):1404–1411. pmid:17028585
  90. 90. Ebner TJ, Pasalar S. Cerebellum predicts the future motor state. The Cerebellum. 2008;7(4):583–588. pmid:18850258
  91. 91. Gaymard B, Rivaud-Péchoux S, Yelnik J, Pidoux B, Ploner C. Involvement of the cerebellar thalamus in human saccade adaptation. Eur J Neurosci. 2001;14(3):554–560. pmid:11553305
  92. 92. Ethier V, Zee DS, Shadmehr R. Spontaneous recovery of motor memory during saccade adaptation. J Neurophysiol. 2008;99(5):2577–2583. pmid:18353917
  93. 93. Heins F, Masselink J, Scherer JN, Lappe M. Adaptive changes to saccade amplitude and target localization do not require pre-saccadic target visibility. Sci Rep. 2023;13(1):8315. pmid:37221275
  94. 94. Taylor JA, Ivry RB. Cerebellar and prefrontal cortex contributions to adaptation, strategies, and reinforcement learning. Prog Brain Res. 2014;210:217–253. pmid:24916295
  95. 95. Smith MA, Ghazizadeh A, Shadmehr R. Interacting adaptive processes with different timescales underlie short-term motor learning. PLoS Biol. 2006;4(6):e179. pmid:16700627
  96. 96. Kording KP, Tenenbaum JB, Shadmehr R. The dynamics of memory as a consequence of optimal adaptation to a changing body. Nat Neurosci. 2007;10(6):779–786. pmid:17496891
  97. 97. Albert ST, Jang J, Sheahan HR, Teunissen L, Vandevoorde K, Herzfeld DJ, et al. An implicit memory of errors limits human sensorimotor adaptation. Nat Hum Behav. 2021;5(7):920–934. pmid:33542527
  98. 98. Orozco SP, Albert ST, Shadmehr R. Adaptive control of movement deceleration during saccades. PLoS Comp Biol. 2021;17(7):e1009176. pmid:34228710
  99. 99. Heald JB, Lengyel M, Wolpert DM. Contextual inference underlies the learning of sensorimotor repertoires. Nature. 2021;600(7889):489–493. pmid:34819674
  100. 100. Cassanello CR, Ostendorf F, Rolfs M. A generative learning model for saccade adaptation. PLoS Comput Biol. 2019;15(8):e1006695. pmid:31398185
  101. 101. Robinson DA. Models of the saccadic eye movement control system. Kybernetik. 1973;14(2):71–83. pmid:4206845
  102. 102. Henson D. Investigation into corrective saccadic eye movements for refixation amplitudes of 10 degrees and below. Vision Res. 1979;19(1):57–61. pmid:419701
  103. 103. Becker W. The neurobiology of saccadic eye movements. Metrics. Rev Oculomot Res. 1989;3:13–67.
  104. 104. Harris CM, Wolpert DM. The main sequence of saccades optimizes speed-accuracy trade-off. Biol Cybern. 2006;95(1):21–29. pmid:16555070
  105. 105. Collins T, Semroud A, Orriols E, Doré-Mazars K. Saccade dynamics before, during, and after saccadic adaptation in humans. IOVS. 2008;49(2):604–612. pmid:18235005
  106. 106. Doya K. What are the computations of the cerebellum, the basal ganglia and the cerebral cortex? Neural Netw. 1999;12(7-8):961–974. pmid:12662639
  107. 107. Widrow B, Hoff M. Adaptive switching circuits, IRE WESCON Conv. Rec, pt. 1960;4:96.
  108. 108. Widrow B, Stearns SD. Adaptive Signal Processing Prentice-Hall. Englewood Cliffs, NJ. 1985.
  109. 109. Deubel H. Adaptivity of gain and direction in oblique saccades. In: Eye movements from physiology to cognition. Elsevier; 1987. p. 181–190.
  110. 110. Frens M, Van Opstal A. Transfer of short-term adaptation in human saccadic eye movements. Exp Brain Res. 1994;100(2):293–306. pmid:7813666