HYPOTHESIS AND THEORY article

Front. Neurol., 04 August 2017

Sec. Neuro-Ophthalmology

Volume 8 - 2017 | https://doi.org/10.3389/fneur.2017.00372

A GABAergic Dysfunction in the Olivary–Cerebellar–Brainstem Network May Cause Eye Oscillations and Body Tremor. II. Model Simulations of Saccadic Eye Oscillations

  • Laboratory of Sensorimotor Research, IRP, National Eye Institute, National Institutes of Health, Bethesda, MD, United States

Abstract

Eye and body oscillations are shared features of several neurological diseases, yet their pathophysiology remains unclear. Recently, we published a report on two tennis players with a novel presentation of eye and body oscillations following self-administration of performance-enhancing substances. Opsoclonus/flutter and limb tremor were diagnosed in both patients. Common causes of opsoclonus/flutter were excluded. High-resolution eye movement recordings from one patient showed novel spindle-shaped, asymmetric saccadic oscillations (at ~3.6 Hz) and ocular tremor (~40–60 Hz). Based on these findings, we proposed that the oscillations are the result of increased GABAA receptor sensitivity in a circuit involving the cerebellum (vermis and fastigial nuclei), the inferior olives, and the brainstem saccade premotor neurons (excitatory and inhibitory burst neurons, and omnipause neurons). We present a mathematical model of the saccadic system, showing that the proposed dysfunction in the network can reproduce the types of saccadic oscillations seen in these patients.

Introduction

Oscillations of the head, body, limbs, or eyes characterize several neurological conditions. Nevertheless, their underlying mechanisms are not understood well enough to guide therapy (1). Oculo- and somatomotor systems are similarly organized, and some disorders involve both eye and body oscillations. The anatomy and physiology of the oculomotor system have been studied more than that of other systems. If we could understand ocular oscillations, it might provide insights into the pathophysiology of oscillatory dysfunctions of somatomotor circuits. Here, we present a mathematical model of the saccadic system that has sufficient anatomical and physiological detail to simulate saccadic eye oscillations in patients with opsoclonus. A model of ocular and limb tremor will be presented in a companion paper.

Previously, we studied two patients with a novel presentation of eye and body oscillations following self-administration of performance-enhancing substances (2). Opsoclonus consists of large saccadic oscillations lasting a few or many cycles, without intersaccadic intervals, around all three axes (36). If oscillations occur only in the horizontal plane it is called flutter. Opsoclonus can be triggered by both saccadic and non-saccadic eye movements, by eye closure, and can persist in the dark (7). How neural circuits generate these oscillations is not clear. Note that in our patient, there is a diverse range of waveforms, including opsoclonus (quasi-sinusoidal movements), square-wave oscillations, and square-pulse oscillations (Figure 1). Such a diversity of waveforms has been seen before, and it has been proposed that a common mechanism accounts for all the waveforms (8), but no mechanism has yet been proposed that can do so.

Figure 1

Opsoclonus arises in various diseases [including paraneoplastic, parainfectious, toxic-metabolic, and idiopathic causes (9)], thus one mechanism may not explain all forms of opsoclonus. Nonetheless, opsoclonus is often associated with cerebellar disease (10, 11). Saccades are generated by high-gain burst neurons that are gated by dominant, inhibitory neurons in the brainstem, which pause during saccades [omnipause neurons (OPN)]. Thus, Zee and Robinson (12) proposed that any saccadic oscillations without an intervening interval would require that the OPN be shut off.

Autopsy of one of Cogan’s patients with opsoclonus found encephalitis with lymphocytic infiltration chiefly in the hypothalamus, midbrain, and pons (3). However, studies have shown no consistent pathology in the raphe interpositus (site of the OPN) in patients with opsoclonus (13). Although structural imaging of the brain of opsoclonus patients has not shown any consistent abnormalities, functional imaging has been linked to increased activity in the fastigial nuclei (14). Oguro et al. (15) recorded SPECT images in two patients with opsoclonus and found hyperperfusion in the midline cerebellum (CB) in one, and hypoperfusion in the other patient.

At least three hypotheses for the pathomechanism of opsoclonus/flutter have been proposed, based on different clinical and experimental observations. One line of evidence would suggest a dysfunction of the cerebellar Purkinje cells (PuC). Mutant mice with a modified glutamate receptor on PuC show clustered PuC action potentials, likely induced by climbing fiber activation, and opsoclonus-like eye movements (16). Jen et al. (17) found antibodies to PuC in patients with opsoclonus. They proposed that those antibodies blocked the parallel fiber (PF) input to PuC, allowing spontaneous oscillations generated in the inferior olives (IO) to be passed to the oculomotor vermis through the flocculus.

Alternatively, disinhibition of caudal fastigial nuclei (cFN) might induce unwanted saccades through excitatory projections to brainstem burst neurons (6, 14). Wong et al. (6) modified a lumped model of the saccadic system by adding a negative feedback path from a high pass filtered efference copy of eye position, through the ipsilateral cFN, to the motor error comparator. When the delay and the gain of this feedback pathway were increased (to 20 ms and ~8.5, respectively) to simulate disinhibition of the cFN, ~15 Hz sinusoidal oscillations occurred. However, this mechanism cannot explain how the brain generates non-sinusoidal saccadic oscillations, such as square waves, square pulses, and spindles. This mechanism also does not suggest how a disease process could change the loop’s delay by 20 ms, or raise its gain from normal (~0.7) to ~8, or how it could hold off the OPN, which prevent saccades. Furthermore, lesions of the vermis that disinhibit the cFN cause hypometric saccades, not oscillations (18, 19).

A third idea is that reduction of glycinergic inhibition generates oscillations in the positive feedback loop between left and right saccadic brainstem inhibitory burst neurons (20). This gives rise to small, sinusoidal oscillations, but alone would not be enough to simulate all the waveforms seen in opsoclonus. Thus, none of these theories have been fully supported by lesion studies in animals, clinical findings, or model simulations (21).

Unfortunately, high temporal resolution recordings of eye movements from patients with opsoclonus/flutter, analysis of which could clarify the underlying mechanisms, are extremely rare, because of the difficulty in recording them (given the severity of the clinical symptoms and/or the inability to calibrate the recordings). In our previous study (2), a detailed analysis of one patient’s movements (Figure 1) suggested that ocular oscillations might be generated by a dysfunction of the cerebellar–olivary–brainstem network.

That patient had taken anabolic–androgenic steroids (AAS), which are allosteric modulators of GABAA receptors (GABAAR), with both acute and chronic effects. These effects can enhance or diminish the sensitivity of the chloride channel to γ-aminobutyric acid (GABA), depending on the receptor’s subunit composition and the GABA concentration at the synaptic level (22, 23). We shall assume that in our patient, the AAS increased the gain of the chloride current in the GABAAR (24). Indeed, the key point of our model is that opsoclonus requires increased inhibition of oculomotor vermis (OMV), cFN, and OPN. Increased inhibition of cFN is not necessarily incompatible with the finding of increased BOLD or fMRI signal of cFN during opsoclonus (14). Increasing cFN inhibition does not necessarily mean that they fire less. As we will see below, in our model the main effect of increasing cFN inhibition is to delay its burst onset, which could be tested in animal models. Thus, even if cFN excitability were significantly reduced, the intense saccadic activity during opsoclonus oscillations (averaged over time) would increase the blood flow and oxygen consumption of the cFN. Our hypothesis of inhibitory receptor dysfunction is very different from prior models, wherein the hypothesis was either damage to the OPN, or damage to the CB causing hyperactivation of the cFN. This is also consistent with the lack of evidence of damage to the CB and pons in opsoclonus. Here, we present model simulations of the cerebellar–olivary–brainstem network that support our hypothesis.

Subjects and Methods

Two patients described in our previous article (2) developed opsoclonus/flutter after self-administration of performance-enhancing substances. This study was carried out in accordance with the recommendations of the ethics committee of the University of Siena, Italy. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The ethics committee of the University of Siena, Italy, approved this study. Briefly, neuro-ophthalmological examinations revealed horizontal saccadic intrusions and intermittent ocular flutter. Both patients’ brain MRI were normal. CSF showed few oligoclonal bands. Common infectious, toxic, paraneoplastic, and metabolic causes of opsoclonus/flutter were excluded by negative blood and CSF exams. No tumors were found.

Both patients reported a few months of abuse of substances to improve performance. Patients stopped using the drugs after the symptoms began. Patient 1 provided a sample of the compound for testing, which led to the identification of the AAS nandrolone, stanozolol, and testosterone propionate. Treatment with intravenous IgG and benzodiazepine led to recovery in 3–4 weeks in both patients.

Eye Movement Recording

Eye movement recording during fixation, horizontal (10 and 18° amplitudes), and vertical (8°) saccades was possible in patient 2 because he had inter-oscillatory intervals of steady fixation that allowed a satisfactory calibration. Eye position was recorded at 240 Hz with an ASL 504 eye-tracker device (Applied Science Laboratories, Bedford, MA, USA). All analyses were performed with custom Matlab (The Mathworks, Natick, MA, USA) scripts.

Model

A new model was created by extending the saccadic system model used to simulate eye movement disorders in a cerebellar patient (25). The basic architecture is an adaptive, velocity feedback, integral controller, but without a motor error comparator (Figure 2). For simplicity, only the horizontal direction of saccades is modeled. Both the superior colliculi (SCs) and the caudal fastigial nuclei (cFN) drive the brainstem premotor burst neurons [PBN; consisting of the short-lead excitatory burst neurons (EBN), inhibitory burst neurons (IBN), and long-lead inhibitory burst neurons (LIBN)]. An efference copy of eye velocity (from the EBN) is fed back to the OMV, which projects to the fastigial oculomotor region in the cFN, closing the loop. A pause in activity starts at a locus in the contraversive (with respect to the saccade direction, e.g., left side for a rightward saccade) OMV corresponding to the amount of drive that the cFN should provide for an individual saccade, based on the context of the movement. If the saccade does not get on target, the CB learns to initiate the pause at a different locus under the same context the next time (2628). At saccade start, the OMV activity at the initial locus pauses, releasing the contraversive cFN from inhibition. The cFN then fires, exciting the ipsiversive EBN and IBN. During the saccade, a wave of inhibition (driven by feedback of a velocity efference copy from the EBN) spreads across the OMV. Thus, the OMV acts as a spatial integrator of eye velocity. When the inhibition spreads to the ipsiversive OMV, it disinhibits the ipsiversive cFN, which activates the contraversive IBN. The contraversive IBN then inhibits the ipsiversive EBN and IBN. This stops the movement. At this point, the excitation from the ipsiversive cFN also reactivates the OPN. The EBN and IBN have a very high gain, thus, without OPN reactivation the saccade would be followed by an oscillation of back-to-back saccades with no intersaccadic interval. The caption to Figure 2 gives a color-coded explanation of how the model makes a saccade. Details of the model are given in the Appendix in Supplementary Material.

Figure 2

Although our main hypothesis is that opsoclonus results from an increased sensitivity of GABAAR, here we simulate a lumped model. This is necessary for two reasons. First, we think that the opsoclonus is caused by abnormal levels of activity in a large network, encompassing the CB, IO, and brainstem. Second, to understand the effects of the GABAAR dysfunction at the biophysical level, we would need to know the types of subunits that make up the receptor in the diverse types of neurons in the circuit, in particular, their gains and time constants. Unfortunately, these values are unknown. Thus, we implement the suspected changes in GABAAR function simply by changing various gains in a lumped model. Details of which gains are changed are given with each of the simulations in the Section “Results.” Model parameters are given in Table 1, and parameters for the different neuron types are given in Table 2.

Table 1

ParameterValueFunction (units)ParameterValueFunction (units)
SynDel0.0008Synaptic delay (s)NI G20Neural integrator gain
VisDelay0.05Visual delay (s)NI Tc20Neural integrator time constant (s)
Refrac0.2Refractory disable (s)Plant Te0.008Small plant time constant (s)
ErrThr0.5Retinal error threshold (°)SC Gdisfa2Disfacilitation gain from cFN
MaxDis800Maximum discharge rate (sp/s)SC Tdisfa0.025cFN disfacilitation time constant (s)
BrstDel0.05SC burst delay (s)SC Gx1cSC cross-inhibitory gain
IO Gl1IO low pass gainSC Tx0.01cSC cross-inhibitory time constant (s)
IO Tc0.1IO low pass time constant (s)rSC Delay0.05Delay time before burst (s)
IO Iu5IO output upper limitOMV Gp1Cerebellar plant model gain
IO Il0.1IO output lower limitOMV Tp0.008Cerebellar plant model time constant (s)
IO PW0 or 1IO noise powerOMV Tf0.25OMV fatigue time constant (s)
IO GW1IO noise gainOMV Gl0.5OMV fatigue low pass gain
IO TW0.002IO noise time constant (s)OMV Tl1.5OMV fatigue low pass time constant (s)
IO Fio30IO sine frequency (Hz)OMV Cl0.5OMV fatigue clip low
IO Gio0.35IO sine gainOMV Cu1OMV fatigue clip high
OMV Gf0.34OMV fatigue gainOMV Glk0.5OMV leaky integrator gain
LIBN Grsc6Inhibitory gain from rSC to LIBNOMV Tlk0.5OMV leaky integrator time constant (s)
OPN Grsc3Excitatory gain from rSC to OPNIBN2EBN1.2Gain from IBN to EBN
OPN Tone200Bias on OPNIBN2IBN0.3Gain from IBN to contra IBN
cFN2EBN0.5Gain from cFN to contra EBNIBN2LIBN200Gain from IBN to LIBN
cFN2IBN50Gain from cFN to contra IBNEBN2IBN0.01Gain from EBN to ipsi IBN

Important model parameters.

LIBN, long-lead inhibitory burst neurons; EBN, excitatory burst neurons; IBN, inhibitory burst neurons; OPN, omnipause neurons; IO, inferior olive; SC, superior colliculus; OMV, oculomotor vermis; cSC, caudal SC; rSC, rostral SC.

Table 2

ParameterEBNIBNLIBNOPNcFNrSCcSC
Gain1.0001.0001.0001.0001.0001.0001.000
Time constant (s)0.0020.0020.0030.0020.0030.0100.010
Adaptation gain0.0500.0500.0000.0100.0001.0001.000
Adaptation TC (s)0.1000.1000.0100.0060.0030.0060.006
Inhibitory gain4.00012.0000.0151.0002.0000.10015.000
Excitatory gain1.0000.8001.0001.0001.0001.0003.100
OPN gain10.000100.000.0000.0000.0001.0001.000

Model parameters for different neuron types.

EBN, excitatory burst neurons; IBN, inhibitory burst neurons; LIBN, long-lead inhibitory burst neurons; OPN, omnipause neurons; cSC, caudal SC; rSC, rostral SC.

Results

Normal Saccades

This model assumes that microsaccades (<2°) and macrosaccades (≥2°) are made by the same circuit (2932). It can simulate amplitudes from 0.5 to 50° in both leftward and rightward directions (Figure S1 in Supplementary Material). For amplitudes above 10°, the cortical circuit assumes that the actual goal of the saccade is only 90% of the target jump. Thus, for large saccades, the model undershoots the target but automatically makes corrective saccades to get on target.

A simulated eye movement and the activity in major model neurons are shown in Figure 3. In this highly simplified, one-dimensional, model of the CB, we can see the latency differences between contra- and ipsiversive OMV and cFN bursts (3335). Indeed, it is our central hypothesis of saccadic system function that the role of the oculomotor vermis is to create this timing difference (2628). To simulate all three dimensions of opsoclonus (horizontal, vertical, and torsional), the model presented here would have to be duplicated, once for each axis.

Figure 3

Opsoclonus/Flutter

Opsoclonus/flutter oscillations can exhibit many waveforms. For example, Figure 1 shows eye movements recorded in one session while patient 2 was fixating (abscissa shows time in seconds from the beginning of a single record). Figures 1A–C show examples of the quasi-sinusoidal movements that define classic opsoclonus. Examples of largely horizontal square-wave macrosaccadic oscillations (SWMSO, half-cycles in blue, marked S), a pulse saccadic oscillation without an intersaccadic interval (red, P), and combined half-cycles of SWMSO followed without an intersaccadic interval by a return saccade, which we call a square-pulse oscillation (SP), are shown in Figure 1D. These different waveforms occurred over an interval of just 3 s. We also see an example where the square-pulse changes direction from rightward to leftward (marked −SP). Note that there seems to be a continuum of intersaccadic interval durations (e.g., saccade marked S at about 69.55 s), thus the difference between opsoclonus, square-wave and square-pulse oscillations may only be due to a small change in a few parameters. The vertical component of the eye movement (green) is small but is phase-locked to the horizontal movement. This phase locking is characteristic of opsoclonus. Square-pulse oscillations have been recorded before in patients with opsoclonus, although not commented upon (6, 14).

Figure 1E shows a 2.5 s example of a square-pulse oscillation that grows in amplitude for about 1 s, and then decreases in amplitude [numbers at each extremum indicate gain of the movement, assuming that the eye is trying to get back to the central fixation target (gray line)]. As in Figure 1D, the vertical eye movement is phase locked to the horizontal eye movement. Because of the shape of the movement envelope, these types of oscillations are called spindles.

Figure 1 also reveals an asymmetry in the oscillations. In Figures 1D,E, we see that the saccades to the left usually have higher gains than those to the right. We infer from this that the left cFN projections to the right IBN and OPN are weaker than the right cFN projections to the left IBN and OPN, causing a delay in stopping the leftward saccade.

Tests of Two Prior Hypotheses

Here, we use the new model to test the hypothesis that fastigial disinhibition by Purkinje cell malfunction results in saccadic oscillations, as in the Wong et al. (6) model. In their model, the role of the CB is like that of most models, in that it accelerates the saccade and stops it on target (26, 27, 36). In their simulations, macrosaccadic flutter (back-to-back saccades with no intersaccadic interval) occurs when the loss of OMV inhibition of the cFN increases the gain and the delay in the feedback loop around the brainstem and through the cFN. Their model cannot make other types of oscillations, such as SWMSO, tremor, or square-pulse oscillations. It also faces problems pointed out with the similar Zee and Robinson model (12), in that to make oscillations requires a change in the loop delay that must be set according to the patient’s oscillation frequency, which can cover a very wide range (37).

In the present model, the CB helps accelerate and steer the saccade, and then stops it on target, but it is an adaptive, velocity feedback, integral controller, and not a motor error controller (i.e., its goal is not to reduce the motor error to 0). However, oscillations are controlled by membrane properties of neurons, and not loop delays (37, 38). Our model behaves differently from that of Wong et al. because its OMV and cFN essentially act as a switching network to control the timing of when the ipsiversive and contraversive EBN and IBN turn on and off. Intuition can help us understand a linear feedback controller, but intuition is not helpful in a switching/timing model. Thus, we need simulations to understand the effects of GABAAR dysfunction in these models.

Before a saccade starts, the locus of initial inhibition in the OMV must be determined (Figure 2, green circle). As the saccade progresses, a wave of inhibition must spread across the OMV (cyan wavy line) until it reaches the ipsiversive side (Figure 2, red diamond), disinhibiting the ipsiversive cFN. Figure 4 shows the effect of making the OMV less active than normal. This will cause the ipsiversive cFN to restart too soon, and the saccade (Figure 4B) will be smaller than normal (Figure 4A). This result contradicts the hypothesis of Wong et al. but is consistent with neurophysiological results, which found that cerebellar vermis lesions result in hypometria, not opsoclonus or sinusoidal oscillations (18, 19).

Figure 4

Another hypothesis of how opsoclonus might be generated suggests that antibodies block PF to PuC synapses, thus reducing inhibition on the flocculus. This allows spontaneous oscillatory activity in the IOs to be passed to the ocular motor nuclei (17). This hypothesis has already been tested in the context of oculopalatal tremor (OPT) (39, 40). OPT has waveforms that look like large, random oscillations, but the movements around each axis are independent. These studies showed that if oscillatory activity from a normal IO projected through the flocculus to the brainstem, the resulting eye movements were very small and had pulsatile waveforms. The development of OPT required a pulsatile oscillator (caused by abnormally tight electrotonic coupling in the IO), and a learned response from the cerebellar cortex to enhance the movement’s gain. Thus, these results suggest that IO oscillations would not, alone, be enough to cause saccadic oscillations in opsoclonus.

Square-Wave Macrosaccadic Oscillations

Figure 1D shows an example of SWMSO in our patient (marked S). SWMSO and spindles also occur in cerebellar disorders (41). Thus, we look for the effects of an increased gain of GABAAR in the CB as the mechanism for causing SWMSO. In the Wong et al. (6) model, the cFN, but not the OMV, is inside the feedback loop, so loss of OMV inhibition to the cFN causes the loop gain to increase. They also had to significantly increase the loop delay (no mechanism for which was proposed). In our model, the behavior is very different, because both the cFN and the OMV are inside the feedback loop and the GABAAR dysfunction results in an increased inhibition in both structures.

Purkinje cells are GABAergic, and we assume that increasing the gain of GABAAR in the OMV would mitigate the loss of GABA from the PuC when they become inhibited at saccade start. This would have the effect of slowing the spread of inhibition during the saccade. Thus, to simulate SWMSO, we reduced the gain of the feedback integration within the OMV (Figure 5A). Here, the saccade gains are less than twice normal size, so they decay in size as the oscillation progresses. Also shown is the gain of the fatigue circuit in the CB, which reduces the gain if neuronal activity is too high for too long. In this case, the fatigue gain declines only slightly during the oscillation (green line).

Figure 5

Figure 5B shows the effect of lowering the OMV feedback gain further, which makes the saccade gain higher, and the circuit unstable. The saccade gain is more than twice normal size, and the oscillations grow for a few saccades. However, the high rate of firing and high frequency of saccades result in marked fatigue of the cerebellar activity (green line), and the saccades begin to shrink. The envelope of the oscillation looks like a spindle.

Although the cFN activity is delayed relative to normal during each saccade, the cFN activity over the whole oscillation is very high. This is consistent with findings from imaging studies that the midline CB (Figures 5C,D) and deep cerebellar nuclei (Figures 5E,F) are strongly activated during opsoclonus (14, 15).

It is important to note the difference between the results in Figure 4 (hypometria) and Figure 5 (hypermetria). In Figure 4, we increased the gain of the OMV’s forward model of the plant (from 1 to 300), which caused a rapid spread of inhibition from contraversive to ipsiversive. This caused the ipsiversive OMV to turn off too soon (removal of OMV inhibition is Wong et al.’s hypothesis), allowing the ipsiversive cFN to turn on too soon, resulting in hypometria. In Figure 5, we decreased the forward model’s gain to represent increased GABAAR activity in OMV, thus slowing the spread of inhibition across the OMV and resulting in hypermetria.

Square-Pulse Macrosaccadic Oscillations and Ocular Flutter

One of the unusual waveforms found in our patient is the square-pulse oscillation (SP in Figure 1D). To simulate this waveform requires two sets of changes. First, the saccades must be hypermetric (as in Figure 5). In addition, it is necessary to delay the onset of the OPN in one direction (here, after leftward saccades). Thus, a rightward movement results in a hypermetric saccade that is followed by a hypermetric leftward saccade. However, after the leftward saccade, the OPN reactivation is delayed, resulting in a return movement with no intersaccadic interval, driven by post-inhibitory rebound (PIR) in the brainstem EBN. Importantly, PIR of the EBN must last at least as long as the return pulse, making the adaptation time constant, Ta, an important parameter (Figure S2 in Supplementary Material). A short delay of OPN reactivation (magenta) allows for a half-cycle (tick marks and red part of trace) of ocular flutter (Figure 6A). If the reactivation of the OPN is further delayed, three half-cycles of flutter can be obtained (Figure 6B). In addition, if the CB does not shut down after the saccade to the left, the OMV and cFN can participate in the oscillation, increasing its amplitude (Figure 6C and compare Figures 6E,F). Thus, random fluctuations in the delay until reactivation of the OPN, and whether or not there is cerebellar involvement, can account for the varying size and number of saccadic pulse waveforms in opsoclonus. Experimental studies will be needed to determine whether EBN PIR and oscillations in the CB contribute to ocular flutter.

Figure 6

What might cause the delay of the OPN reactivation? Under our hypothesis, the opsoclonus is caused by an abnormally high GABAAR gain. If the cFN were abnormally inhibited, saccades would be hypermetric, because the ipsiversive cFN would not turn on in time to stop the saccade on target. They might also be too weak to turn the OPN back on at the end of the saccade. OPN receive both GABA and Gly inhibitory transmitters (42). If their GABAAR currents were also enhanced, it might take more excitation to reactivate the OPN, causing them to turn on late. This is consistent with the inference from Figure 1 that the left cFN is weaker than the right cFN.

If the OPN reactivation is delayed even more, as may happen during blinks or large off-vertical saccades (37), ocular flutter (back-to-back saccades with no intersaccadic interval) occurs (as in Figure 6B). Pathological flutter has been associated with lesions in the region of the OPN, and with lesions of the projection from cFN to the brainstem (43, 44). Here, we see that the functional effect that causes flutter is the delay in reactivation of the OPN.

Figure 7 shows a combination of the effect of lowering cerebellar feedback gain and fatigue (as in Figure 5) and slowing OPN reactivation (as in Figure 6). This results in a square-pulse oscillation spindle (Figure 7A, like the one in Figure 1E). Here, the OPN are only delayed after leftward saccades. In Figure 7B, we simulate a classic opsoclonus oscillation by delaying OPN after both leftward and rightward saccades. This simulation looks very much like the patient’s eye movement in Figure 1B. The patient’s quasi-sinusoidal oscillations were usually short (0.5–2.5 cycles), with amplitude about 5–10°. There were very few of these movements, compared to square-pulse oscillations. Importantly, these waveforms are not pure sinusoids. They appear to be back-to-back saccades, which is how they were simulated.

Figure 7

Discussion

The purpose of this model is to show the possible interactions between brain circuits that determine eye movement waveforms. We found that it is not the absolute level of activity, but the relative timing of different areas in the brain, which is important. Here, although the activities in the lumped neurons have highly simplified waveforms, they have the correct timing.

The main predictions of this model are that opsoclonus saccadic waveforms can result from increased GABA inhibition of neurons in the vermis, fastigial nuclei, and brainstem (omnipause and burst neurons), due to increased sensitivity of GABAAR.

GABAAR Mechanisms

The GABAAR consists of five subunits with many subtypes (α1–6, β1–3, γ1–3, δ, ε, π, θ, and ρ1–3) providing diverse receptor functions (45, 46). AAS are allosteric modulators of the GABAAR. Their modulatory activity depends upon which steroid is administered and the GABAAR subunit composition, being greater for the α2 than the α1 subunit, but also acting through the δ and ε subunits (22, 24, 47). Subunit expression is different in different brain areas but is poorly defined because of the lack of specificity of markers and the inattention to cell types important for eye movements, e.g., OPN. IO dendrites contain the α2 subunit, but their somas contain α3 subunits (23, 48, 49). The α2 subunit is more prominent than the α1 subunit in the cerebellar granule cell and molecular layers (50). Purkinje cells and brainstem reticular formation express the α1 subunit, and deep cerebellar nuclei express both subunits (45, 51). We assume that increasing GABAAR modulation and decreasing excitation from cFN to OPN would delay OPN reactivation. Despite these speculations, exactly how AAS affected the GABAAR in our patients cannot be known. A more detailed biophysical model of opsoclonus awaits further experiments on the effects of GABAAR dysfunction in identified cell types.

We have shown how our hypothesis that GABAAR are dysfunctional in a cerebellar–olivary–brainstem network (2) can be implemented by varying parameters in the model. Here, we focused on dysfunction that resulted in a higher gain of GABAAR. Opsoclonus has many causes, and so other mechanisms besides GABAAR modulation may also cause the diverse types of opsoclonus seen in other patients. Whatever the underlying biophysical mechanism, our model allows us to hypothesize how different regions in the brain must be affected to obtain the various waveforms observed in opsoclonus.

Opsoclonus and Oscillatory Eye Movements

Given the large differences in waveforms seen in patients with opsoclonus, including quasi-sinusoidal, square-wave and square-pulse waveforms, it is not surprising that they have been regarded as different types of movements. However, we have shown that a single model can simulate all of these types of movements, simply by changing a few parameters. Thus, we agree with the hypothesis of Ellenberger et al. (8), who emphasized that these waveforms occur together in the same patients, and thus might be unified as dyskinesias of the saccadic system. Furthermore, although the classic waveform for opsoclonus is a large, sinusoidal oscillation, these may, in fact, simply be back-to-back saccades with no intervening interval, i.e., quasi-sinusoidal oscillations. With low bandwidth recordings (e.g., from electroculograms) these waveforms would be low pass filtered and thus would look sinusoidal. However, higher quality recordings from video or eye coil systems reveal their quasi-sinusoidal nature.

These waveforms can all be obtained from a model of the saccadic system by making appropriate parameter changes in both cerebellar and brainstem circuits, but not by making changes in either alone. We infer from our model that the mixture of opsoclonus (quasi-sinusoidal), square-wave and square-pulse oscillations resulting from cerebellar/brainstem dysfunction may commonly co-occur in patients with opsoclonus. Thus, as earlier studies have argued, opsoclonus is not caused by a cerebellar deficit alone. Nor is it caused by lesions of the OPN region. Instead, we hypothesize that opsoclonus occurs when neuronal activity in the CB and brainstem are mistimed.

Statements

Ethics statement

This study was carried out in accordance with the recommendations of the ethics committee of the University of Siena, Italy. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The ethics committee of the University of Siena, Italy, approved this study.

Author contributions

EP proposed the underlying theory. LO implemented the model and ran the simulations. Both authors wrote the manuscript.

Funding

This research was supported by the Intramural Research Program of the National Eye Institute, NIH.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at http://journal.frontiersin.org/article/10.3389/fneur.2017.00372/full#supplementary-material.

Figure S1

Family of saccades simulated with normal parameter values. (A) Saccades from −50 to +50°. Note that saccades larger than 10° undershoot the target by 10%. The model automatically makes a corrective saccade that gets on target (gray lines). (B) Inset magnifies saccades to left and right 0.5 and 2° targets. (C) Velocity traces for saccades in panels (A,B). Peak speeds were approximately matched to those found by Clark and Stark (52) in normal human subjects.

Figure S2

Simulation of post-inhibitory rebound (PIR) in a model neuron (Eq. 1). PIR depends heavily on three parameters: the gain and time constant of the adaptation element, and the amount of hyperpolarization before the rebound. Here, the excitatory input was set to 0. The neuron has two inhibitory inputs, one for γ-aminobutyric acid (GABA) and one for glycine [NB: the omnipause neurons (OPN) and the IBN are glycinergic, and some LIBN are GABAergic]. The total inhibition to the neuron is thus the sum of the GABA and Gly inputs. For convenience, the sum is set to 1. GABAergic inhibitory input changed briefly from 0.8 to 0. Glycinergic OPN inhibition changed briefly from 0.2 to 0. Thus, the total inhibition to the cell changed from 1.0 to 0. (A) Effect of changing the adaptation gain (Ga) on the peak of the rebound activity (adaptive time constant was set to 6 ms). (B) Effect of changing the adaptive time constant (Ta, with adaptive gain set to 1.0). Other parameters for the neuron were: Gl = 1, Tl = 2 ms, Go = 1.0, δ = 0.8 ms. As the time constant increases from 2 to 21 ms, the amplitude and width of the rebound activity increase. Above 21 ms, however, the amplitude begins to decrease. The duration of the rebound is truncated when the OPN resume firing. (C) Effect of reducing γ-aminobutyric acid (GABA) inhibition on rebound amplitude. The maximum amount of GABA ranged from 0.0 to 0.8 (accounting for from 0 to 80% of the inhibition on the neuron). As the proportion of GABA was reduced, the hyperpolarization decreased, and the rebound amplitude decreased.

Abbreviations

BN, SC burst neurons; BUN, SC buildup neurons; CB, cerebellum; cFN, caudal fastigial nucleus; cSC, caudal SC; EBN, excitatory burst neurons; GABA, γ-aminobutyric acid; PuC, Purkinje cells (vermis); Glu, glutamate; Gly, glycine; IBN, inhibitory burst neurons; IO, inferior olive; LIBN, long-lead inhibitory burst neurons; NRTP, nucleus reticularis tegmenti pontis; NRTPc, NRTP → contra OMV; NRTPi, NRTP → ipsi OMV; PF, parallel fiber; OMV, oculomotor vermis; OPN, omnipause neurons; PIR, post-inhibitory rebound; rSC, rostral SC; SC, superior colliculus; SWMSO, square-wave macrosaccadic oscillations; FOR, fastigial oculomotor region.

References

  • 1

    HallettM. Tremor: pathophysiology. Parkinsonism Relat Disord (2014) 20(Suppl 1):S11822.10.1016/S1353-8020(13)70029-4

  • 2

    PretegianiERosiniFRocchiRGinanneschiFVinciguerraCOpticanLMet alGABAAergic dysfunction in the olivary-cerebellar-brainstem network may cause eye oscillations and body tremor. Clin Neurophysiol (2017) 128(3):40810.10.1016/j.clinph.2016.12.014

  • 3

    CoganDG. Ocular dysmetria; flutter-like oscillations of the eyes, and opsoclonus. AMA Arch Ophthalmol (1954) 51(3):31835.10.1001/archopht.1954.00920040324006

  • 4

    CoganDG. Opsoclonus, body tremulousness, and benign encephalitis. Arch Ophthalmol (1968) 79(5):54551.10.1001/archopht.1968.03850040547009

  • 5

    Averbuch-HellerLRemlerB. Opsoclonus. Semin Neurol (1996) 16(1):216.10.1055/s-2008-1040955

  • 6

    WongAMMusallamSTomlinsonRDShannonPSharpeJA. Opsoclonus in three dimensions: oculographic, neuropathologic and modelling correlates. J Neurol Sci (2001) 189(1–2):7181.10.1016/S0022-510X(01)00564-0

  • 7

    ZangemeisterWHMuller-JensenAZschockeS. Benign encephalitis: electro-oculographic analysis of opsoclonus. J Neurol (1979) 222(2):95108.10.1007/BF00313003

  • 8

    EllenbergerCJrKeltnerJLStroudMH. Ocular dyskinesia in cerebellar disease. Evidence for the similarity of opsoclonus, ocular dysmetria and flutter-like oscillations. Brain (1972) 95(4):68592.10.1093/brain/95.4.685

  • 9

    WongA. An update on opsoclonus. Curr Opin Neurol (2007) 20(1):2531.10.1097/WCO.0b013e3280126b51

  • 10

    RossATZemanW. Opsoclonus, occult carcinoma, and chemical pathology in dentate nuclei. Arch Neurol (1967) 17(5):54651.10.1001/archneur.1967.00470290100013

  • 11

    EllenbergerCJrCampaJFNetskyMG. Opsoclonus and parenchymatous degeneration of the cerebellum. The cerebellar origin of an abnormal ocular movement. Neurology (1968) 18(11):10416.10.1212/WNL.18.11.1041

  • 12

    ZeeDSRobinsonDA. A hypothetical explanation of saccadic oscillations. Ann Neurol (1979) 5(5):40514.10.1002/ana.410050502

  • 13

    RidleyAKennardCScholtzCLButtner-EnneverJASummersBTurnbullA. Omnipause neurons in two cases of opsoclonus associated with oat cell carcinoma of the lung. Brain (1987) 110(Pt 6):1699709.10.1093/brain/110.6.1699

  • 14

    HelmchenCRamboldHSprengerAErdmannCBinkofskiFfMRI study. Cerebellar activation in opsoclonus: an fMRI study. Neurology (2003) 61(3):4125.10.1212/01.WNL.0000073271.66866.51

  • 15

    OguroKKobayashiJAibaHHojoH. Opsoclonus-myoclonus syndrome with abnormal single photon emission computed tomography imaging. Pediatr Neurol (1997) 16(4):3346.10.1016/S0887-8994(97)00031-3

  • 16

    YoshidaTKatohAOhtsukiGMishinaMHiranoT. Oscillating Purkinje neuron activity causing involuntary eye movement in a mutant mouse deficient in the glutamate receptor delta2 subunit. J Neurosci (2004) 24(10):24408.10.1523/JNEUROSCI.0783-03.2004

  • 17

    JenJCLopezIBalohRW. Opsoclonus: clinical and immunological features. J Neurol Sci (2012) 320(1–2):615.10.1016/j.jns.2012.06.017

  • 18

    TakagiMZeeDSTamargoRJ. Effects of lesions of the oculomotor vermis on eye movements in primate: saccades. J Neurophysiol (1998) 80(4):191131.

  • 19

    BarashSMelikyanASivakovAZhangMGlicksteinMThierP. Saccadic dysmetria and adaptation after lesions of the cerebellar cortex. J Neurosci (1999) 19(24):109319.

  • 20

    ShaikhAGMiuraKOpticanLMRamatSLeighRJZeeDS. A new familial disease of saccadic oscillations and limb tremor provides clues to mechanisms of common tremor disorders. Brain (2007) 130(Pt 11):302031.10.1093/brain/awm240

  • 21

    LemosJEggenbergerE. Saccadic intrusions: review and update. Curr Opin Neurol (2013) 26(1):5966.10.1097/WCO.0b013e32835c5e1d

  • 22

    YangPJonesBLHendersonLP. Role of the alpha subunit in the modulation of GABA(A) receptors by anabolic androgenic steroids. Neuropharmacology (2005) 49(3):30016.10.1016/j.neuropharm.2005.03.017

  • 23

    OberlanderJGPenattiCAPorterDMHendersonLP. The Buzz about anabolic androgenic steroids: electrophysiological effects in excitable tissues. Neuroendocrinology (2012) 96(2):14151.10.1159/000339123

  • 24

    BitranDHilversRJFryeCAErskineMS. Chronic anabolic-androgenic steroid treatment affects brain GABA(A) receptor-gated chloride ion transport. Life Sci (1996) 58(7):57383.10.1016/0024-3205(95)02326-7

  • 25

    DayePMOpticanLMRozeEGaymardBPougetP. Neuromimetic model of saccades for localizing deficits in an atypical eye-movement pathology. J Transl Med (2013) 11:125.10.1186/1479-5876-11-125

  • 26

    LefevrePQuaiaCOpticanLM. Distributed model of control of saccades by superior colliculus and cerebellum. Neural Netw (1998) 11(7–8):117590.10.1016/S0893-6080(98)00071-9

  • 27

    QuaiaCLefevrePOpticanLM. Model of the control of saccades by superior colliculus and cerebellum. J Neurophysiol (1999) 82(2):9991018.

  • 28

    OpticanLMQuaiaC. Distributed model of collicular and cerebellar function during saccades. Ann N Y Acad Sci (2002) 956:16477.10.1111/j.1749-6632.2002.tb02817.x

  • 29

    Otero-MillanJTroncosoXGMacknikSLSerrano-PedrazaIMartinez-CondeS. Saccades and microsaccades during visual fixation, exploration, and search: foundations for a common saccadic generator. J Vis (2008) 8(14):118.10.1167/8.14.21

  • 30

    HafedZMGoffartLKrauzlisRJ. A neural mechanism for microsaccade generation in the primate superior colliculus. Science (2009) 323(5916):9403.10.1126/science.1166112

  • 31

    Otero-MillanJMacknikSLSerraALeighRJMartinez-CondeS. Triggering mechanisms in microsaccade and saccade generation: a novel proposal. Ann N Y Acad Sci (2011) 1233:10716.10.1111/j.1749-6632.2011.06177.x

  • 32

    SunZJunkerMDickePWThierP. Individual neurons in the caudal fastigial oculomotor region convey information on both macro- and microsaccades. Eur J Neurosci (2016) 44(8):253142.10.1111/ejn.13289

  • 33

    OhtsukaKNodaH. Saccadic burst neurons in the oculomotor region of the fastigial nucleus of macaque monkeys. J Neurophysiol (1991) 65(6):142234.

  • 34

    FuchsAFRobinsonFRStraubeA. Role of the caudal fastigial nucleus in saccade generation. I. Neuronal discharge pattern. J Neurophysiol (1993) 70(5):172340.

  • 35

    OhtsukaKNodaH. Discharge properties of Purkinje cells in the oculomotor vermis during visually guided saccades in the macaque monkey. J Neurophysiol (1995) 74(5):182840.

  • 36

    DeanP. Modelling the role of the cerebellar fastigial nuclei in producing accurate saccades: the importance of burst timing. Neuroscience (1995) 68(4):105977.10.1016/0306-4522(95)00239-F

  • 37

    RamatSLeighRJZeeDSOpticanLM. Ocular oscillations generated by coupling of brainstem excitatory and inhibitory saccadic burst neurons. Exp Brain Res (2005) 160(1):89106.10.1007/s00221-004-1989-8

  • 38

    RamatSLeighRJZeeDSShaikhAGOpticanLM. Applying saccade models to account for oscillations. Prog Brain Res (2008) 171:12330.10.1016/S0079-6123(08)00616-X

  • 39

    HongSLeighRJZeeDSOpticanLM. Inferior olive hypertrophy and cerebellar learning are both needed to explain ocular oscillations in oculopalatal tremor. Prog Brain Res (2008) 171:21926.10.1016/S0079-6123(08)00631-6

  • 40

    ShaikhAGHongSLiaoKTianJSolomonDZeeDSet alOculopalatal tremor explained by a model of inferior olivary hypertrophy and cerebellar plasticity. Brain (2010) 133(Pt 3):92340.10.1093/brain/awp323

  • 41

    SelhorstJBStarkLOchsALHoytWF. Disorders in cerebellar ocular motor control. II. Macrosaccadic oscillation. An oculographic, control system and clinico-anatomical analysis. Brain (1976) 99(3):50922.10.1093/brain/99.3.509

  • 42

    HornAKButtner-EnneverJAWahlePReichenbergerI. Neurotransmitter profile of saccadic omnipause neurons in nucleus raphe interpositus. J Neurosci (1994) 14(4):203246.

  • 43

    AsheJHainTCZeeDSSchatzNJ. Microsaccadic flutter. Brain (1991) 114(Pt 1B):46172.10.1093/brain/114.1.461

  • 44

    TsutsumiTMurakamiMKawaishiJChidaWWatanabeK. Ocular flutter associated with a lesion of the right upper pontine tegmentum. Auris Nasus Larynx (2009) 36(6):6957.10.1016/j.anl.2009.02.007

  • 45

    FritschyJMBenkeDMertensSOertelWHBachiTMohlerH. Five subtypes of type A gamma-aminobutyric acid receptors identified in neurons by double and triple immunofluorescence staining with subunit-specific antibodies. Proc Natl Acad Sci U S A (1992) 89(15):672630.10.1073/pnas.89.15.6726

  • 46

    SieghartW. Structure, pharmacology, and function of GABAA receptor subtypes. Adv Pharmacol (2006) 54:23163.10.1016/S1054-3589(06)54010-4

  • 47

    JonesBLWhitingPJHendersonLP. Mechanisms of anabolic androgenic steroid inhibition of mammalian epsilon-subunit-containing GABAA receptors. J Physiol (2006) 573(Pt 3):57193.10.1113/jphysiol.2006.106534

  • 48

    LaurieDJSeeburgPHWisdenW. The distribution of 13 GABAA receptor subunit mRNAs in the rat brain. II. Olfactory bulb and cerebellum. J Neurosci (1992) 12(3):106376.

  • 49

    DevorAFritschyJMYaromY. Spatial distribution and subunit composition of GABA(A) receptors in the inferior olivary nucleus. J Neurophysiol (2001) 85(4):168696.

  • 50

    SieghartWSperkG. Subunit composition, distribution and function of GABA(A) receptor subtypes. Curr Top Med Chem (2002) 2(8):795816.10.2174/1568026023393507

  • 51

    PirkerSSchwarzerCWieselthalerASieghartWSperkG. GABA(A) receptors: immunocytochemical distribution of 13 subunits in the adult rat brain. Neuroscience (2000) 101(4):81550.10.1016/S0306-4522(00)00442-5

  • 52

    ClarkMRStarkL. Control of human eye movements: III. Dynamic characteristics of the eye tracking mechanism. Math Biosci (1974) 20(3–4):23965.10.1016/0025-5564(74)90003-0

Summary

Keywords

saccade, vermis, fastigial nuclei, inferior olive, omnipause neurons, eye movement, flutter, opsoclonus

Citation

Optican LM and Pretegiani E (2017) A GABAergic Dysfunction in the Olivary–Cerebellar–Brainstem Network May Cause Eye Oscillations and Body Tremor. II. Model Simulations of Saccadic Eye Oscillations. Front. Neurol. 8:372. doi: 10.3389/fneur.2017.00372

Received

14 April 2017

Accepted

17 July 2017

Published

04 August 2017

Volume

8 - 2017

Edited by

Janine Leah Johnston, University of Manitoba, Canada

Reviewed by

Stefano Ramat, University of Pavia, Italy; Kenichiro Miura, Kyoto University, Japan

Updates

Copyright

*Correspondence: Lance M. Optican,

Specialty section: This article was submitted to Neuro-Ophthalmology, a section of the journal Frontiers in Neurology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics