Inferior olivary activity causes both short-term and long-term changes in cerebellar output underlying motor performance and motor learning. Many of its neurons engage in coherent subthreshold oscillations and are extensively coupled via gap junctions. Studies in reduced preparations suggest that these properties promote rhythmic, synchronized output. However, the interaction of these properties with torrential synaptic inputs in awake behaving animals is not well understood. Here we combine electrophysiological recordings in awake mice with a realistic tissue-scale computational model of the inferior olive to study the relative impact of intrinsic and extrinsic mechanisms governing its activity. Our data and model suggest that if subthreshold oscillations are present in the awake state, the period of these oscillations will be transient and variable. Accordingly, by using different temporal patterns of sensory stimulation, we found that complex spike rhythmicity was readily evoked but limited to short intervals of no more than a few hundred milliseconds and that the periodicity of this rhythmic activity was not fixed but dynamically related to the synaptic input to the inferior olive as well as to motor output. In contrast, in the long-term, the average olivary spiking activity was not affected by the strength and duration of the sensory stimulation, while the level of gap junctional coupling determined the stiffness of the rhythmic activity in the olivary network during its dynamic response to sensory modulation. Thus, interactions between intrinsic properties and extrinsic inputs can explain the variations of spiking activity of olivary neurons, providing a temporal framework for the creation of both the short-term and long-term changes in cerebellar output.
A multitude of behavioral studies leave little doubt that the olivo-cerebellar system organizes appropriate timing in motor behavior [1–3], perceptual function [4–6] and motor learning [7–10]. Furthermore, the role of the inferior olive in motor function is evinced in (permanent and transient) clinical manifestations, such as tremors, resulting from olivary lesions and deficits [11–16]. Although the consequences of olivary dysfunctions are rather clear, the network dynamics producing functional behavior are controversial. At the core of the controversy is the question whether inferior olive cells are oscillating during the awake state and whether these oscillations affect the timing of the inferior olivary output [17–19]. The inferior olive is the sole source of the climbing fibers, the activity of which dictates complex spike firing by cerebellar Purkinje cells (for review, see [20]). Climbing fiber activity is essential for motor coordination, as it contributes to both initiation and learning of movements [8,10,21–26], and it may also be involved in sensory processing and regulating more cognitive tasks [27–30]. Understanding the systemic consequences of inferior olivary spiking is therefore of great importance. The dendritic spines of inferior olivary neurons are grouped in glomeruli, in which they are coupled by numerous gap junctions [10,31–33], which broadcast the activity state of olivary neurons. Due to their specific set of conductances [34–40], the neurons of the inferior olive can produce subthreshold oscillations (STOs) [41–43]. The occurrence of STOs does not require gap junctions per se [44], but the gap junctions appear to affect the amplitude of STOs and engage larger networks in synchronous oscillation [10,16,42]. Both experimental and theoretical studies have demonstrated that STOs may mediate phase-dependent gating where the phase of the STO helps to determine whether excitatory input can or cannot evoke a spike [45,46]. Indeed, whole cell recordings of olivary neurons in the anesthetized preparation indicate that their STOs can contribute to the firing rhythm [42,43] and extracellular recordings of Purkinje cells in the cerebellar cortex under anesthesia often show periods of complex spike firing around the typical olivary rhythm of 10 Hz [17,47–49]. However, several attempts to capture clues to these putative oscillations in the absence of anesthesia have, so far, returned empty handed [19,50]. It has been shown that in the anesthetized state both the amplitude and phase of the STOs can be altered by synaptic inputs [10]. Inhibitory inputs to the inferior olive originate in the cerebellar nuclei and have broadly distributed terminals onto compact sets of olivary cells [51–53]. Excitatory terminals predominantly originate in the spinal cord and lower brainstem, mainly carrying sensory information, and in the nuclei of the meso-diencephalic junction in the higher brainstem, carrying higher-order input from the cerebral cortex (Fig 1A) [15,54,55]. In addition, the inferior olive receives modulating, depolarizing, level-setting inputs from areas like the raphe nuclei [55]. Unlike most other brain regions, the inferior olive is virtually devoid of interneurons [56,57]. Thus, the long-range projections to the inferior olive in conjunction with presumed STOs and gap junctions jointly determine the activity pattern of the complex spikes in Purkinje cells. How these factors contribute to functional dynamics of the inferior olive in awake mammals remains to be elucidated. Here, we combine recordings in awake mice–in the presence and absence of gap junctions–with network simulations using a novel inferior olivary model to study the functional relevance of STOs in terms of resonant spikes. We are led to propose a view of inferior olivary function that is more consistent with the interplay between STOs, gap junctions and inputs to the inferior olive. Rather than acting as a strictly periodic metronome, the inferior olive appears more adequately described as a quasiperiodic ratchet, where cycles with variable short-lasting periods erase long-term phase dependencies.
To study the conditions for, and consequences of, rhythmic activity of the inferior olive, we made single-unit recordings of cerebellar Purkinje cells in lobules crus 1 and crus 2 (n = 52 cells in 16 awake mice) in the presence and absence of short-duration (30 ms) whisker air puff stimulation (Fig 1B and 1C). In the absence of sensory stimulation, the complex spikes of 35% of the Purkinje cells (18 out of 52) showed rhythmic activity (Fig 2A–2C; S1 Fig) with a median frequency of 8. 5 Hz (inter-quartile range (IQR): 4. 7–11. 9 Hz). Upon sensory stimulation, 46 out of the 52 cells (88%) showed statistically significant complex spike responses. Of these, 31 (67%) had sensory-induced rhythmicity (Fig 2D–2F), which was a significantly larger proportion than during spontaneous behavior (p = 0. 002; Fisher’s exact test). The median frequency of the oscillatory activity following stimulation was 9. 1 Hz (IQR: 7. 9–13. 3 Hz). Hence, the preferred frequencies in the presence and absence of sensory stimulation were similar (p = 0. 22; Wilcoxon rank sum test) (Figs 2C, 2F and 3). The duration of the enhanced rhythmicity following stimulation was relatively short in that it lasted no more than 250 ms. With our stringent Z-score criterion (>3), only a single neuron showed 3 consecutive significant peaks in the peri-stimulus time histogram (PSTH). The minimum inter-complex spike interval (ICSI) across cells was around 50 ms, putatively representing the refractory period. We conclude that complex spikes also display rhythmicity in awake behaving mice, and that sensory stimulation can amplify these resonances in periods of a few hundred milliseconds, even though stimulation is not required for the occurrence of rhythmicity per se. The pattern of rhythmic complex spike responses that was apparent for a couple of hundred milliseconds after a particular air puff stimulus repeated itself in a stable manner across the 1,000 trials (applied at 0. 25 Hz) during which we recorded (Fig 4A and 4B). For example, the level of rhythmicity of the first 100 trials was not significantly different from that during the last 100 trials (comparing spike counts in first PSTH peak, p = 0. 824; χ2 test, or latency to first spike, p = 0. 727, t test). This strongly indicates that there is–in a substantial fraction of the Purkinje cells–persistent oscillatory gating of the probability for complex spikes after a sensory stimulus resulting in time intervals (“windows of opportunity”) during which complex spikes preferentially occur (Fig 4B–4E). These windows of opportunity become even more apparent when sorting the trials on the basis of response latency: the first complex spikes with a long latency following the stimulus align with the second spikes of the short latency responses. Similarly, there are trials during which complex spikes appear only at the third cycle (Fig 4C, seen as a steeper rise around trial no. 650). The occurrence of spikes during later cycles, not predicated on prior spikes, argues against refractory periods or rebound spiking as the sole explanations for such rhythmic firing [58] and highlight the putative existence of network-wide coherent oscillations. Since sensory stimulation of the whiskers can trigger a reflexive whisker protraction [59–61] and complex spike firing is known to correlate with the amplitude of this protraction [61], we examined the relation between periodic complex spike firing and whisker protraction. To this end, we further analyzed an existing dataset of simultaneously recorded Purkinje cells and whisker movements during 0. 5 Hz air puff stimulation of the ipsilateral whisker pad. In line with our previous findings [61], trials during which a single complex spike occurred within 100 ms of whisker pad stimulation had on average a slightly, but significantly stronger protraction (from 6. 1 ± 5. 4° to 6. 8 ± 5. 3° (medians ± IQR), n = 35 Purkinje cells, p = 0. 033, Wilcoxon-matched pairs test after Benjamini-Hochberg correction; S2A–S2C Fig). Our new analysis revealed that also the occurrence of a second complex spike was correlated with a stronger whisker protraction. This could be observed as a second period of increased protraction during trials with two complex spikes. When compared to the increase in trials with a single complex spike, this second protraction was highly significant (p < 0. 001, Wilcoxon-matched pairs test after Benjamini-Hochberg correction; S2D Fig). The second complex spike was unlikely a mere reflection of stronger protraction following the first complex spike, as there was no difference in whisker protraction between the trials that had a complex spike during the first, but not the second 100 ms after stimulus onset, and the trials with the opposite pattern (a complex spike during the second, but not the first 100 ms; p = 0. 980, Wilcoxon-matched pairs test after Benjamini-Hochberg correction). The rhythmic firing pattern of complex spikes was thus reflected in the behavioral output of mice. The existence of windows of opportunity for complex spike activity is compatible with the assumption of an underlying STO, and cannot solely be explained by rebound activity without invoking circuit-wide extrinsic mechanisms. To test the implications of assuming olivary STOs, we proceeded to reproduce a detailed network with a tissue-scale computer model of the inferior olive neuropil. The model is constituted by 200 biophysically plausible model cells [40,46,62] embedded in a topographically arranged 3D-grid (Fig 5A–5C). It has the scale of a sheet of olivary neurons of about 10% of the murine principal olivary nucleus (cf. [63]). The model was designed to test hypotheses about the interaction between intrinsic parameters of olivary neurons, such as STOs and gap junctional coupling, and extrinsic parameters including synaptic inputs during the generation of complex spike patterns. Each neuron in the model is composed of a somatic, an axonal and a dendritic compartment, each endowed with a particular set of conductances, including a somatic low threshold Ca2+ channel (Cav3. 1; T-type), a dendritic high threshold Ca2+ channel (Cav2. 1; P/Q-type) and a dendritic Ca2+-activated K+ channel, chiefly regulating STO amplitudes, while a somatic HCN channel partially determines the STO period (Fig 5B; see also Methods). The dendrites of each neuron are connected to the dendrites of, on average, eight nearby neighbors (within a radius of three nodes in the grid, representing a patch of about 400 μm x 400 μm of the murine inferior olive), simulating anisotropic and local gap junctional coupling (Fig 5C). As the inferior olive itself, our model has boundaries which have impact on local connectivity characteristics, such as the clustering coefficient, though these did not have significant impact on the average firing rate between edge and center cells (p = 0. 812, comparing edge and center cells, Kolmogorov-Smirnov test; S3 Fig). The coupling coefficient between model cells varied between 0 and 10%, as reported for experimental data [45,64,65]. Sensory input was implemented as excitatory synaptic input, simulating the whisker signals originating from the sensory trigeminal nuclei that were synchronously delivered to a subset of model neurons. Additionally, a “contextual input” was implemented as a combination of inhibitory feedback from the cerebellar nuclei and a level setting modulating input (Figs 1A and 5A). This contextual input is modeled after an Ohrstein-Uhlenbeck process, essentially a random exploration with a decay parameter that imposes a well-defined mean yet with controllable temporal correlations (see Methods). The amplitude of the contextual input drives the firing rate of the model neurons, which we set around 1 Hz (S4 Fig), corresponding to what has been observed in vivo [28,43,66]. Thus, our model network recapitulates at least part of the neural behavior observed in vivo due to biophysically plausible settings of intrinsic conductances, gap junctional coupling and synaptic inputs. Whether a model neuron at rest displays STOs or not is largely determined by its channel conductances. Activation of somatic T-type Ca2+ channels can trigger dendritic Ca2+-dependent K+ channels that can induce Ih, which in turn can again activate T-type Ca2+ channels, and so forth. This cyclic pattern can cause STOs that could occasionally produce spikes (Fig 5B). In our model, the conductance parameters were randomized (within limits, see Methods) so as to obtain an approximate 1: 3 ratio of oscillating to non-oscillating cells (S5 Fig) guided by proportions observed in vivo [43]. Sensitivity analysis with smaller ratios (down to 1: 5) did not qualitatively alter the results (data and analyses scripts are available online in https: //osf. io/9hmpy/). In the absence of contextual input, model neurons were relatively silent, but when triggered by sensory input, as occurred in our behavioral data (Figs 2 and 4), STOs synchronized by gap junctions would occur for two or three cycles (Fig 6A and 6B). Our network model confirms that gap junctional coupling can broaden the distribution of STO frequencies and that even non-oscillating cells may, when coupled, collectively act as oscillators (S6 Fig) [67]. Adding contextual input to the model network can lead to more spontaneous spiking in between two sensory stimuli. Compared to the situation in the absence of contextual input, the STOs are much less prominent and the post-spike reverberation is even shorter (Fig 6C). Accordingly, despite the significant levels of correlation in the contextual input (10%), the periods between oscillations are more variable due to the interaction of the noisy current and the phase response properties of the network. In addition, in the presence of contextual input our model could readily reproduce the appearance of preferred time windows for spiking upon sensory stimulation as observed in vivo (Fig 6D–6F, cf. Fig 4). This was particularly true for the model cells that directly received sensory inputs (Fig 6G–6H). Moreover, the observed rhythmicity in model cells as observed in their STO activity was in tune with that of the auto-correlogram (Fig 6D–6I) in that the timing of the STOs and that of the spiking were closely correlated (cf. Fig 5B). It should be noted though that model cells adjacent to cells directly receiving sensory input showed only a minor effect of stimulation. Thus, even though the gap junction currents in the model were chosen as the ceiling physiological value for the coupling coefficient (≤10%) [45,67], these currents alone were not enough to trigger spikes in neighboring cells. Both directly stimulated model cells and those receiving only contextual input exhibited phase preferences, seen in the spike-triggered membrane potential average as well as in the spike-triggered average of the input currents (Fig 6G–6I). Spike-triggered averages of membrane potentials for any cell showed depolarization followed by hyperpolarization. In contrast, trials in which no spike was generated showed a depolarization just before the occurrence of the input. Similarly, the average of the input showed a long-lived phase preference, not only for a hyperpolarization before the spike, but also a preference for a depolarization in the previous peak of the STO, more than 100 ms earlier. These results are in line in vitro experiments under dynamic clamp and noisy input [68,69]. Likewise, the model indicates that for short durations STOs can induce clear phase dependencies for spiking, which fades under the variation of period durations dependent on the trial-specific contextual input (as seen in our data). Depolarizing sensory input delivered onto a subset of the model cells can reset the STO phase in oscillatory cells and create resonant transients in others (Fig 7A and 7B; see also S7 Fig on the appearance of rebound firing). If a second stimulus is delivered during this short-lasting transient, the response probability is increased. As in most cells with resonant short-lasting dynamics, inputs delivered during different phases can cause phase advances or delays. Hyperpolarization advanced the phase between 0−π and delayed the phase between π−2π, whereas depolarization had roughly the opposite effect, in addition to phase advancements with spikes in later cycles between π−2π (Fig 7C–7E). Thus, there is a mutual influence of synaptic inputs and STOs on periodicity. While STOs can lead to phase-dependent gating, synaptic input can either modulate or reset the phase of the STOs, generating variable periods that range between 40–160 ms for the chosen amplitude of the contextual input (Fig 8; S6 Fig). The only means of settling the question about the prevalence of STOs in awake and behaving mice would be intracellular recordings of inferior olivary neurons, which remains a daunting experimental challenge. We therefore looked for a less invasive method that could read out, from indirect and infrequent complex spikes, the presence or absence of STOs. We have developed one such paradigm inspired by auditory studies [70,71] using a rhythmic gallop stimulus that we first applied to the network model (Fig 7F). In the gallop paradigm, stimuli are applied in quick succession with alternating intervals, comparable to the putative period of the underlying oscillation. Enough stimuli should be applied such that after multiple presentations the stimuli sample a uniform distribution of phases. In the context of auditory stimuli, the standard gallop experiment involves different tones and is used to test perceptual separation of auditory streams. Such rhythmic stimuli can help indicate resonances or physical limitations of the system, and distinguish across possible models for this separation (such as in neural resonance theory [71]). One possible mechanism of auditory stream separation is an underlying oscillatory process which resets in certain phases and is less responsive in others. According to the in vitro inferior olive literature [41–43] this behavior is to be expected, and hence, such a stimulus can help distinguish underlying processes. If spikes are modulated by an oscillatory process, the presence of spikes on a short interval should be able to predict, in the next interval, the absence or presence of spikes. Indeed, if the underlying process producing spikes has oscillatory components and a relatively stable period, the probability of spikes in each interval is systematically different, which would appear as asymmetric ratios of response in the different intervals (Fig 7F). This can be inspected as the length of the empty and filled vertical bars representing ratio of probabilities of spiking for long or short stimulation intervals (Fig 7G). Thus, if the period of the STO rhythm would be regular and cause phase-dependent gating, complex spike responses following each stimulus interval are expected to show preferences for the short or long window of stimulation; these preferences were indeed observed (Fig 7G, left). However, this clear phase dependency only appears in the noiseless model scenario. After adding a moderate amount of contextual input, this dependency washes off, rendering the responses in the two windows more symmetric (Fig 7G, right), with only a few cells (5/200) displaying significant ratio differences (tested against bootstrap with shuffled spikes). In line with the experimental in vivo data (e. g. , Figs 2 and 4), the olivary spike rhythmicity in the network model was steadily present over longer periods, and for a wide range of contextual input parameters (S5 Fig). In addition, it also comprised, as in the experimental data, variations in frequency and amplitude during shorter epochs (Figs 8 and S6). Analysis of the network parameters indicates that these latter variations in oscillatory behavior can be readily understood by their sensitivity to both the amplitude (parameter' sigma' ) and kinetics (parameter' tau' , temporal decay) of the contextual input. Indeed, because of the underlying Ornstein-Uhlenbeck process, the generation of contextual input converges to a specified mean and standard deviation, but in short intervals the statistics including the average network STO frequency can drift considerably (Figs 6C and 8). Since relatively small differences in oscillation parameters such as frequency can accumulate, they can swiftly overrule longer-term dependencies created by periodically resetting stimuli, as an analysis of phase distributions shows (Fig 9). Thus, based upon the similar outcomes of the network model and in vivo experiments, we are led to propose that (1) the STOs in the inferior olive may well contribute to the continuous generation of short-lasting patterns of complex spikes in awake behaving animals, and that (2) the synaptic input to the inferior olive may modify the main parameters of these STOs. Note that in the absence of input, periodic rhythmic behavior should be the default behavior of oscillating cells. Thus, in all likelihood, even if the inferior olive oscillates endemically, sustained but variable input should induce highly contextual spike responses to variable periods and render the olivary responses quasiperiodic, rather than regularly periodic as observed in reduced preparations. In line with in vivo whole cell recordings made under anesthesia [10,43] our awake data support the possibility that the moment of spiking may be related to the phase of olivary STOs, especially during the period of several hundred ms following stimulation (Figs 2 and 4). As discussed above, a gallop stimulus would expose such an oscillatory process underlying the response probabilities. Four idealized scenarios about the expected results can be constructed, as follows: first, one can start with a complete absence of STOs, which would result in a response probability unrelated to stimulus intervals; second, it could be that there were STOs, but no phase-dependent firing (to be expected if the STO amplitude is small), which would also lead to complex spike firing irrespective of stimulus intervals; third, there could be STOs, but each stimulus would evoke a phase-reset, which again, would not lead to interval dependencies; and fourth, there could be STOs in combination with phase-dependent gating, which would result in a clear dependency of complex spike firing on the previous interval length (Fig 7G, left panel). It should be noted that the large majority of studies on inferior olivary physiology, especially in reduced preparations, found evidence for the fourth situation (STOs + phase-dependent gating) [41,43,67]. To study whether phase-dependent gating in conjunction with an underlying oscillatory process could shape complex spike response timing in vivo we applied both a 250 vs. 400 ms and a 250 vs. 300 ms gallop stimulation using air puffs to the whiskers. Using only trials with a CS in the previous trial to calculate the ratio of responses (' conditional firing' ) a slight bias could be observed in the 250 vs. 400 ms paradigm (Fig 10A) and to a lesser extent in the 250 vs. 300 ms (Fig 10B). Analysis including all trials (' non-conditional' ) is included in S8 Fig and shows no significant bias for any of the cells tested. Hence, our in vivo data are in line with the results from the network model subjected to synaptic noise, and show that the timing of complex spike responses to sensory stimulation is biased but not strongly determined by STOs. Our experimental data provided evidence for phase-dependent complex spike firing during brief intervals, but gallop stimulation did not expose a strong impact of STOs on complex spike response probabilities. Therefore, we sought an alternative approach to study the impact–if any–of STOs on complex spike firing in vivo. We reasoned that, if a sensory stimulus triggers a complex spike response with a certain probability, higher stimulation frequency should result in a proportional increase in complex spike firing. In particular, stimulus frequencies that would be in phase with the underlying STO would be expected to show signs of resonance and result in disproportionally increased complex spike firing. However, over periods of tens of seconds the complex spike frequency was resilient to varying the stimulus frequency between 1 and 4 Hz (linear regression = -0. 02; R2 = 0. 1) (Fig 10C) and did not show signs of resonance with any of the stimulus frequencies, as there were no frequencies at which the complex spike firing was substantially increased. Only a very high rate of sensory stimulation (10 Hz), commensurable with the average duration of windows of opportunity, could induce a mild increase in complex spike firing frequency, albeit at the cost of a highly reduced response probability (average increase: 71 ± 64% corresponding to an average increase from 1. 12 Hz to 1. 92 Hz; n = 5; p < 0. 05; paired t test). This examination indicates that the average complex spike frequency is robust and stiff to modulation over longer time periods, imposing a hard limit on the frequency with which complex spikes can respond to sensory stimuli, confirming recent reports on complex spike homeostasis [72]. As stimulus triggered resonances were not observed at any of the stimulation frequencies, we turned to a more sensitive measure for the detection of oscillatory components in complex spike firing. We developed a statistical model that extrapolated from frequencies inferred through inter-complex spike intervals and stimulus triggered histograms (Fig 11A and 11B). We reasoned that phase-dependent gating would imply that the interval between the last complex spike before and the first one after sensory stimulation aligns to the preferred frequency. In contrast, if sensory stimulation would typically evoke a phase reset, as suggested by our network model (Fig 7), no such relation would be found. The method was applied only to Purkinje cells with highly rhythmic complex spike firing. For each of those, we calculated their preferred frequency in the absence (Fig 11C) or presence of sensory stimulation (Fig 11D). We used that frequency to construct statistical models representing idealized extremes of phase-dependent (oscillatory) and -independent (uniform) responses. For the oscillatory component we employed an oscillatory gating model, where the timing of the first complex spike after stimulus onset would be in-phase with the ongoing oscillation. This model was contrasted to a linear response model in which sensory stimulus could evoke a complex spike independent of the moment of the last complex spike before that stimulus, apart from a refractory period. For each Purkinje cell, we compared the distribution of the intervals between the last complex spike before and the first complex spike after stimulus onset with the predicted distributions based on the linear model, the oscillatory model and nine intermediate models, mixing linear and oscillatory components with different relative weights (Fig 11A–11E). For the two extreme models as well as for the nine intermediate models we calculated a goodness-of-fit per Purkinje cell. Overall, when using these relatively long periods (300 ms), the linear model was superior to the oscillatory model, although a contribution of the oscillatory model could often improve the goodness-of-fit (Fig 11F–11H). Despite the apparent failure of the oscillatory model to fit the data, the data did show an oscillatory profile for many of the cells (Fig 11E). This lends support to our observations that short-lived, but reliable, oscillations are apparent in complex spike timing, although they have little impact on the timing or probability of sensory triggered CS responses. Apart from the STOs, extensive gap junctional coupling between dendrites is a second defining feature of the cyto-architecture of the inferior olive [31,33,73]. Absence of these gap junctions leads to relatively mild, but present deficits in reflex-like behavior and learning thereof [10,74]. We analyzed the inter-complex spike interval times in Purkinje cells of mutant mice that lack the Gjd2 (Cx36) protein and are hence unable to form functional gap junctions in their inferior olive [69]. In line with the predictions made by our network model (Figs 5H and S6), the absence of gap junctions did not quench rhythmic complex spike firing during spontaneous activity (Fig 12A). In fact, the fraction of Purkinje cells showing significant rhythmicity in the Gjd2 KO mice was larger than that in the wild-type littermates (Gjd2 KO: 38 out of 65 Purkinje cells (58%) vs. WT: 15 out of 46 Purkinje cells (33%); p = 0. 0118; Fisher’s exact test), with their average rhythmicity being significantly stronger (p = 0. 003; Kolmogorov-Smirnov test), measured by Z-scores of side peaks (Fig 12B). Indeed, the variation in oscillatory frequencies across Purkinje cells of the mutants was significantly less than that in their wild-type littermates in that the latency to peak times per Purkinje cell were less variable (p = 0. 0431; Mann-Whitney test; Fig 12C). This latter finding is at first sight contradictory to our findings in the network model, where we show that gap junctions promote more uniform firing rates through increased synchrony between neurons (Fig 5H). These simulations were run in the absence of synaptic input, though. Addition of contextual input also creates more variability in the wild type cells (S6B Fig). As the lack of gap junctions increases cell excitability [10,69], it is likely that synaptic input has a larger impact in the absence of gap junctions, leaving less room for inter-cell heterogeneity. Overall, removal of gap junctions affected the temporal and spatial dynamics by increasing the stereotypical rhythmicity of complex spike firing. We made paired recordings of Purkinje cells in awake mice to study the temporal relations of their complex spikes during spontaneous activity. The cell pairs were recorded with two electrodes randomly placed in a grid of 8 x 4, with 300 μm between electrode centers. For each pair of simultaneously recorded Purkinje cells, we made a cross-correlogram. The median number of complex spikes in the reference cell used for these cross-correlograms was 827 (range: 74–2174). Cell pairs showed coherent activity in that they could show a central peak and/or a side peak in their cross-correlogram (Fig 13A–13C). The side peaks could appear at different latencies, similar to the range observed in auto-correlograms of single Purkinje cells (cf. Fig 2B). Moreover, Purkinje cell pairs that did not produce signs of synchronous spiking in the center peak could still produce an “echo” in the side peak after 50–150 ms. Counter-intuitively, cross-correlograms of Purkinje cell pairs of the wild type mice showed less often a significant center peak than those of Gjd2 KOs (WT: 51 out of 96 pairs (53%; N = 4 mice); Gjd2 KO: 44 out of 61 pairs (72%; N = 7 mice); p = 0. 0305; Fisher’s exact test). In line with the more stereotypic firing observed in single cells in the absence of gap junctions (Fig 12), the strength of the center peak was on average enhanced in the mutants (Z-scores of significant center peaks (median ± IQR): WT: 3. 47 ± 1. 82; Gjd2 KO: 5. 75 ± 5. 58; p = 0. 0002; Mann-Whitney test) (Fig 13D–13F). Instead, the side peak of Gjd2 KO Purkinje cell pairs was not stronger than that of WTs (Z-scores of significant side peaks (median ± IQR): WT: 3. 01 ± 0. 89; Gjd2 KO: 3. 04 ± 1. 52; p = 0. 194; Mann-Whitney test), leading to a lower ratio between center and side peak (mean ± SEM: WT: 90. 70 ± 5. 17%; Gjd2 KO: 72. 86 ± 6. 52%; p = 0. 036; t = 2. 143; df = 67; t test) (Fig 13E and 13F). Interestingly, the occurrence of side peaks in Purkinje cell pairs was unidirectional in approximately half the cell pairs (WT: 47 out of 82 pairs with at least one side peak (57%); Gjd2 KO: 25 out of 47 pairs (53%); p = 0. 714; Fisher’s exact test), which means that one of the neurons of a pair was leading the other, but not vice versa. As this was consistent in the Gjd2 KO as well as the WT Purkinje cells, these data could reflect traveling waves across the inferior olive, which, however, must have extrinsic sources [44,75]. Thus, the paired recordings are compatible with the findings highlighted above in that the presence of coupling can affect the coherence of STOs for short periods up to a few hundred milliseconds, while leaving the window for later correlated events open.
During spontaneous activity, Purkinje cells generally fire a complex spike roughly once a second, but this frequency can be increased to about 10 Hz by systemically applying drugs, like harmaline, which directly affect conductances mediating STOs in the inferior olive [41,81]. Since these drugs also induce tremorgenic movements beating at similar frequencies, it has been proposed that the inferior olive may serve as a temporal framework for motor coordination [11,82]. This oscillatory firing behavior of the olivary neurons may mirror limb resonant properties and act as an inverse controller, for example by dampening the dynamics of the muscles involved [83]. In line with previous recordings [22–24,28,61,72,76,78,84], the current data indicate that only a small fraction of Purkinje neurons respond to sensory stimulation with a complex spike response probability larger than 50%. This probability falls substantially with increasing frequency of stimulation, as the overall spike frequency only marginally increases to high frequency stimulation. Even after applying different temporal patterns of sensory stimulation for longer epochs, we observed no substantial deviation from the stereotypic 1 Hz firing rate. Moreover, it should be noted that even if the frequency of underlying oscillations has bearing on the pattern of responses of the gallop stimuli, conditional dependencies should be expected for most STO frequencies, unless the ratio of the interval of gallop and the STO period has no remainder. Given the seemingly consistent frequencies predicted by PSTH' s and autocorrelograms of single cells (Figs 2 and 3, but also seen in cross correlograms, as in Fig 13), we chose gallop intervals with periods commensurate with a representative frequency of 8 Hz, each of which should sample different phases in the oscillation. If at all present, we should have observed conditional dependencies on at least a few cells. In our study, complex spikes remain as unpredictable as ever. Thus, regulatory mechanisms keep the complex spike rate relatively stable over longer time periods [72]. No resonance is exhibited, irrespective of an enduring powerful sensory stimulus in a variety of frequencies. Save few exceptions, the presence of a complex spike in an interval is compensated by the absence in another. Thus, it looks as if the complex spikes rearrange themselves in time in order to keep close to its proverbial 1 Hz frequency. It remains to be shown to what extent the mechanisms involved are intrinsic (cell-dependent) and/or extrinsic (network-dependent). A possible candidate for setting the overall level of excitability through intrinsic mechanisms is given by Ca2+-activated Cl− channels, which are prominently expressed in olivary neurons along with Ca2+-dependent BK and SK K+ channels [85,86]. In addition, the olivo-cerebellar module itself could partly impose this regulation [86–88]. Indeed, the long-term dynamics within the closed olivo-cortico-nuclear loop may well exert homeostatic control, given that increases in complex spikes lead to enhanced inhibition of the inferior olive via the cerebellar nuclei [20]. The impact of such a network mechanism may even be more prominent when changes in synchrony are taken into account [89]. We propose that a closed-loop experiment conducted while imaging from a wide field, producing stimulation as a function of the degree of complex spike synchrony, could tease out conditional complex spike probabilities. Increasing our capability of predicting complex spikes is instrumental to elucidate the control of inferior olivary firing. The existence of temporal windows of opportunity for complex spike responses following sensory stimulation highlighted a potential impact of STOs on conditional complex spike gating [10,41,43,90]. Indeed, autocorrelogram peaks correlated well with interspike intervals following stimulation, arguing for an underlying rhythm. Complex spikes could appear in a particular window even when they were not preceded by a complex spike in a previous window during a single trial, arguing against a prominent role of refractory periods in creating rhythmic complex spike responses. Comparing actual firing patterns with statistical models mixing linear or oscillatory interval distributions indicated a potential impact of oscillations. The mild impact of the oscillatory component on explaining the data may in part depend on the assumption that cells have a well-defined frequency. In other words, a variable rebound time could offset the phase response by a couple of milliseconds, reducing the contribution of the oscillatory model, though phase preferences due to prior spikes may still occur (i. e. , Fig 11E). Our biophysical model suggests that fluctuating inputs, such as those mediating inhibition from the cerebellar nuclei or those relaying depolarizing modulation from the raphe nuclei [91,92], may induce variations in the oscillation period on a cycle-by-cycle basis (Figs 8 and 9). As these contextual inputs are absent or suppressed in decerebrate or anesthetized preparations, as well as in vitro, they may also explain why many earlier studies systematically encountered cells with well-defined STO frequencies [10,41,43,45,47,79,93,94]. In the network model, in which we mimicked the contextual input as an Ornstein-Uhlenbeck process with local variations but no long-term drifts of the mean [95], the results agree well with the experimental observations in terms of synchronous firing, phase shifts, cross-correlogram peaks and side peaks, as well as overall firing frequency. Indeed, the absence of resonant responses over longer time windows and the inconsistency of individual olivary cells to fire on every trial or cycle indicate that the STOs are not regularly periodic, but rather quasiperiodic, while still being synchronous. Even though several lines of evidence suggest a role for STOs (see above), we did not observe an unequivocal, significant conditional dependence of complex spikes in the gallop paradigm, as expected by a noiseless model. How can a system with rhythmic responses at least partially fail to be phase modulated by such periodic stimuli? An attractive alternative explanation for rhythmicity might be the occurrence of high-threshold Cav2. 1 P/Q-type Ca2+ channel-dependent rebound spikes (S7 Fig) [12,62]. If impulse-like input to the olive can evoke a spike, and if this spike produces a rebound spike some tens of milliseconds later, this could explain the alignment between the PSTHs and cross-correlograms. However, this argument cannot explain stimulus triggered spikes at the second or third window of opportunity, without an earlier spike as observed in Fig 4. As the occurrence of the rebound spike is predicated on a prior spike, a spike in the second or third window without a prior spike cannot be explained by the rebound spiking phenomenon, at least not within the same cell. In other words, the spikes happening exclusively in the second (or third) window of opportunity cannot be the result of a previous spike in the same cell, unless there is a shared rhythm in the network. It is also conceivable that strong hyperpolarization that is synchronized with the complex spike rhythm could promote reverberating firing, but this is an extrinsic mechanism, discussed below. As they stand, our findings do not support the idea that the post-spike hyperpolarization is a prerequisite for the complex spike pattern observed. Multiple windows of opportunity could, according to our model, be enhanced by transient oscillations induced by resets relayed by gap junctions to the local olivary circuit. Apart from the almost complete absence of interneurons, the presence of STOs and the exclusive projection to the cerebellum, the abundance of dendro-dendritic gap junctions is another defining feature of the inferior olive. The absence of these gap junctions does not lead to gross motor deficits, but prevents proper acquisition and execution of more challenging tasks [10,16,74], which is in line with the relatively minor impact found on complex spike activity in Gjd2 KO mice. At first sight, the effects of deleting gap junctions seem counterintuitive. Synchronous and rhythmic patterns are exacerbated, rather than diminished by the loss of gap junctions. However, the side peak of the auto-correlogram is significantly squashed, indicating that the gap junctions have a role in the increased coherence of the upcoming oscillation. Gap junctions do not only facilitate synchronization of coupled neurons, they also lower their excitability by increasing the membrane resistance [69]. Together, this results in less direct coupling, observed as reduced synchrony of direct neighbors [16,96], and increased responsiveness to synaptic input. This leads to more long-range coherence and as a consequence gap junction networks may act as a “noise filter”, promoting short-range quorum-voting on phase (a term coined by Winfree [97]). This effect is visible in our model as spikes are most likely to occur when excitation follows inhibition (Fig 6H). This is in line with the finding that complex spikes of nearby Purkinje cells have a preference to fire together [72,98,99]. This concept also agrees with the possibility that coupled olivary neurons may control movements by dampening the dynamics of the muscles involved at an appropriate level [83,100], as both the resonances and movement oscillations increase shortly after sensory stimulation in Gjd2 KO mice [16]. Network resonances are a pervasive feature of brain circuits and they can be induced by subthreshold oscillations of particular cell types [101,102]. In addition to the autochthonous dynamics of the inferior olive, reverberating loops through the circuit could help explain some features of complex spike firing, including the occurrence of complex spike doublets and side peaks in cross-correlograms. Such phenomena could be explained by" network echoes" , where complex spikes in one cycle would induce complex spikes in the next cycle [87,103,104]. The most obvious candidate loop to produce is that via the cerebellum and the nuclei of the meso-diencephalic junction [55,105]. The output of the inferior olive is mainly directed via exceptionally strong synapses to the Purkinje cells [106]. These Purkinje cells in turn inhibit neurons of the cerebellar nuclei that can show rebound firing after a period of inhibition [88,107]. This rebound activity can excite the inferior olive again via a disynaptic connection via the nuclei of the meso-diencephalic junction. While an isolated complex spike is unlikely to evoke such a rebound activity, a larger group of Purkinje cells could be successful in doing so [20,107,108]. The travel time for this loop (around 50–100 ms) has been indirectly assessed in the awake preparation [8,10,26], and corresponds to the latency of the rebound firing in the cerebellar nuclei under anesthesia [87,88,107,109]. This implies that the travel times for the entire loop would be in the same order as found for the preferred frequencies of complex spike firing. Other, more elaborate loops involving for instance the forebrain may also exist [110] and could play an additional role in shaping complex spike patterns. A putative impact of reverberating loops on rebound activity could be a network phenomenon, as the impact of an isolated complex spike may not be sufficient to trigger this loop. This is in line with the reduced “echo” in the cross-correlograms of the Gjd2 KO mice and enhanced doublets following lesions of the nucleo-olivary tract as occurs in olivary hypertrophy [111]. Taken together, rebound spiking, STOs and reverberating loops all seem to promote in a cooperative manner complex spike rhythmicity at a time scale of about 200 ms. Through modeling, we found that not only the state of the inferior olivary oscillations determines which inputs are transmitted, but that these inputs also determine the state of the network. Thus, inputs from both the cerebellum and the cerebrum determine the probability of complex spike responses on a cycle-by-cycle basis providing a quasiperiodic framework to align synchronous groups. This sharply contrasts with a view in which the inferior olive is a clock with regular periodicity. A circuit-wide understanding of cerebellar resonances on the basis of such a mechanism could open a novel pathway to explore the cerebellar gating by other brain regions. The combination of delayed gap junctions and delayed inhibition, as found in the olivo-cerebellar loop [104,112], can affect oscillatory behavior [113]. The interplay between STOs and delayed inhibition is therefore also relevant for other neural circuitries, for instance for creating filter settings for the perception of sounds with specific oscillatory properties [114–116] or orchestrating rhythmic movements as shown in the present study (see also [117]). Well-coordinated movement sequences are not timed rigidly; they must be enacted flexibly and contextually. In order to catch a ball, or a prey, or to perform any other appropriately timed movement, it is essential to fine-tune the duration and onsets of multiple coordinated output systems. An inferior olive that responds contextually to time varying input by advancing and delaying cycles does not act as a rigid clock or metronome, but more contextually, as a ratchet-pole system, with the frequency of' clicks' of the ratchet reflecting the recent history of applied torque. The properties we have encountered in this study are consistent with a' ratchet-like' dynamics for the inferior olive, which integrates time-varying stimulus in a phase-dependent manner. According to this view, the inferior olive responds to all inputs (sensory and otherwise), by producing phase changes that are informative about the recent history of input, and dictate the appearance of coherent complex spike waves arriving at the cerebellar cortex.
All experimental procedures were approved a priori by an independent animal ethical committee (DEC-Consult, Soest, The Netherlands) as required under Dutch law. Experiments were performed on 16 adult (9 males and 7 females of 25 ± 14 weeks old) homozygous Gjd2tm1Kwi (Gjd2 KO, formerly known as Cx36 KO mice [10]) mice which were compared to 15 wild-type littermates (8 males and 7 females of 26 ± 13 weeks old; means ± sd). The generation of these mice has been described previously [118]. The data described in S2 Fig originated from previously published recordings in 35 wild-type mice [61]. All mice had a C57BL6/J background. The mice received a magnetic pedestal that was attached to the skull above bregma using Optibond adhesive (Kerr Corporation, Orange, CA) and a craniotomy of the occipital bone above lobules crus 1 and crus 2. The surgery was performed under isoflurane anesthesia (2–4% V/V in O2). Post-surgical pain was treated with 5 mg/kg carprofen (“Rimadyl”, Pfizer, New York, NY) and 1 μg lidocaine (Braun, Meisingen, Germany). Mice were habituated during 2 daily sessions of 30–60 min. Extracellular recordings of Purkinje cells were made in the cerebellar lobules crus 1 and 2 of awake mice as described previously [28]. Briefly, an 8 x 4 matrix of quartz-platinum electrodes (2–4 MΩ; Thomas Recording, Giessen, Germany) was used to make recordings that were amplified and digitized at 24 kHz using an RZ2 BioAmp processor (Tucker-Davis Technologies, Alachua, FL). The signals were analyzed offline with SpikeTrain (Neurasmus, Rotterdam, The Netherlands) using a digital band-pass filter (30–6,000 Hz). Complex spikes were recognized based on their waveform consisting of an initial spike followed by one or more spikelets. A recording was accepted as that of a single Purkinje cell when a discernible pause of at least 8 ms in simple spike firing followed the complex spikes and when the complex spikes were of similar shape and amplitude throughout the recording. Sensory stimulation was applied as air puffs of 20 psi and 25 ms duration directed at the whisker pad ipsilateral to the side of recording. The stimuli were given in trains of 100 or 360 pulses either at regular or alternating intervals. During a recording, trains with different stimulus intervals were played in a random sequence. Whisker videos were made from above using a bright LED panel as backlight (λ = 640 nm) at a frame rate of 1,000 Hz (480 x 500 pixels using an A504k camera from Basler Vision Technologies, Ahrensburg, Germany). The whiskers were not trimmed or cut. Whisker movements were tracked offline as described previously [119] using a method based on the BIOTACT Whisker Tracking Tool [120]. We used the average angle of all trackable large facial whiskers for further quantification of whisker behavior. Of each Purkinje cell we computed the probability density function (PDF) of both its complex spike autocorrelogram and its distribution of intervals between consecutive complex spikes (inter-complex spike intervals (ICSIs) ). PDFs were calculated with an Epanechnikov kernel (with finite support) with a width of 10 ms. In order to exclude stimulus-induced alterations in complex spike firing, complex spikes detected between 20 and 200 ms after a stimulus were omitted from this phase of the analysis. PDFs were calculated from 0 up till 500 ms. The peak in the ICSI PDF was considered as the “preferred ICSI interval” and its strength was expressed as the Z-score by dividing the peak value by the standard deviation of the PDF. To understand the impact of Purkinje cell with little or no preference for specific ICSI intervals on the analysis, we chose to look both at the Purkinje cells with high and low Z-scores. Thus, we grouped Purkinje cells into high and low level Z-scores, using a threshold of 3. Air puff stimulations frequently triggered double complex spike response peaks, suggestive of an underlying inferior olivary oscillation. For further analysis of the conditional responses, an estimate of the putative inferior olivary frequency was derived from the interval between these two response peaks. First, it was established for each Purkinje cell whether two peaks were present in the PSTH. To this end, we set a threshold for each of these two peaks. For the first peak, this was calculated by reshuffling the ICSIs over the recording followed by calculating a stimulus-triggered pseudo-PSTH, repeating this procedure 10,000 times and selecting the 99% upper-bound. We considered the first response peak to be significant if it crossed the upper-bound uninterruptedly for at least 10 ms. Since the second response peak typically is much smaller than the first one, we calculated a new threshold for the second peak by excluding the time-window for the first response peak. This window was set from the time of the stimulus until where the response probability drops to the average response frequency, the response frequency as expected if stimuli do not trigger complex spikes, following the significant ‘first’ responsive peak. In 5 out of 98 Purkinje cells, the PDF of the response rate between clear peaks remained above the average response rate, in which case we used the time point where the amplitude drop in the PDF was more than twice the difference between upper bound and average response probability. The rest of the bootstrap method was identical to that for the first response peak. Only peaks up to 0. 5 s after the stimulus were included in the population analysis. In order to test whether the phase of the inferior olivary oscillations affected the complex spike response probability, we compared the complex spike intervals over an Air puff for each stimulus that triggered a complex spike. To this end, we analyzed the recordings of 25 Purkinje cells (10 WT and 15 Gjd2 KO) that were measured previously in crus 1 and crus 2 of awake, adult mice. We included only Purkinje cells that displayed clear oscillatory complex spike firing indicated by the display of a secondary complex spike response peak, as evaluated according to the bootstrap method described above, and/or significant peaks in the ICSI histogram. Only stable recordings covering at least 500 stimuli at frequencies below 1 Hz were considered for this analysis. For each recording we compared two idealized statistical models of the observed ICSI distributions: an oscillatory model showing phase-dependent spiking and stable olivary oscillation frequencies and a uniform model lacking phase-dependencies. For the oscillatory model, we created complex spike probability functions for the pre-stimulus interval (-300 to 0 ms) based on the oscillatory period established either for the ICSI distribution or from the interval between the two complex spike response peaks. We fitted a sine wave with the observed frequency, having its peak at the moment of the first complex spike in the stimulus response window (20–200 ms after the air puff) and derived spike probability levels during the pre-stimulus interval from these fits, with the trough representing zero probability. Frequency and amplitude of every cycle were kept constant for the whole recording. In the uniform model, we calculated the pre-stimulus spiking probability with a uniform distribution based on the complex spike frequency of each Purkinje cell. We did include a refractory period, being the shortest ICSI observed for each recording, to reflect the inability of consecutive complex spikes to occur with a very short time interval. Refractory periods were comparable between mutants and wild types; 49 ± 15 ms for Gjd2 KO cells and 50 ± 20 ms for WT cells. Subsequently, we constructed compound fits consisting of linear summations of the two models. One extreme was the oscillatory model and the other the uniform model and we considered nine intermediate combinations (e. g. , 0. 3 x the oscillatory model + 0. 7 x the uniform model). Every compound fit was run for 10,000 times. The goodness of fit was computed as the absolute differences of every single run of the model with the actual ICSI distribution. The model networks used here are comprised of a topographical grid of 200 coupled cells, in a 10x10x2 lattice arrangement, which may resemble an area of about 400 μm x 400 μm of the inferior olive, for instance, the rostral portion of the dorsal lamella of the principle olive of the mouse. It is available online at https: //github. com/MRIO/OliveTree, branch' Warnaar' . For instructions on how to run the model and reproduce analysis, check README_Warnaar. txt' . Each cell within these networks was modelled according to the single cell model described in [46], which is an elaboration of a previous model [62] with an added axon and modified fast sodium channel. Equations are provided in the appendix of that publication at (https: //doi. org/10. 1371/journal. pcbi. 1002814. s002), and can be checked in the MATLAB functions IOcell and createDefaultNeurons in the codebase. The model includes three compartments (soma, dendrite, axon hillock) with 12 conductances. In addition to the ionic mechanisms, the dendrite of the model cell has a Ca2+ concentration state variable, which is related to the intrusion through the Cav2. 1 channels. The main ionic conductances responsible for the oscillation are the somatic T-type Ca2+ and the Ca2+-activated K+ (SK) channels present in the dendritic compartment. The crucial parameters governing the emergence of subthreshold oscillations are randomized, reflecting the experimental facts that about one third of the cells oscillate endemically (in vivo) with intrinsic variations in oscillatory frequencies [43]. The behavior of the STO of the model cells in our network as a function of their parameters, for models with and without gap junctions, are included in S5 and S6 Figs. Cell parameters are found in S1 Table. Connectivity is created with the function' createW. m' in the MATLAB codebase. Briefly, cells within a specified radius of each other were connected according to a probability function such as to ensure the specified mean degree in the network (n = 8), chosen to resemble the observed connectivity distributions reported in the literature. The connectivity parameters (distance and average connection probability) were chosen to match experimental values (radius ≤ 120 μm) and average connection probability (~8 neighbors). The procedure to obtain connectivity is as follows. First, pairwise distances between all cells are calculated. Then, a binary adjacency matrix is created by thresholding those distances within a specified radius. Thereafter, we assign a random number between 0 and 1 to each link from a uniform distribution. Finally, this matrix is made binary by comparing each entry with a probability so that the average number of connections approximates a given mean connectivity. This binary adjacency matrix is then multiplied by the mean gap junction conductance parameter. Finally, gap junction conductance values are then randomized by a uniform jittering of the conductance by 10% of their original value. The conductance of gap junctions was normalized with a saturating factor by difference of potential between the neighboring cells, according to [62] based on findings from [121] with the following function: gc¯ (∆V) =gc (0. 8e (−∆V2/100) +0. 2) FORMULA 1 Where ΔV is the voltage difference between the connected cells, gc is the nominal coupling and gc¯ is the effective coupling. Two inputs are given to the model, one emulating the sensory input from whisker pad stimulation and the other representing a stimulus-independent background reflecting diverse excitatory and inhibitory inputs to the inferior olive. The latter consisted of a continuous stochastic process with known mean and standard deviation with a relaxation parameter following the Ornstein-Uhlenbeck process [95], succinctly described underneath. Only one subset of cells in the center of the network (40% of the cells in a mask spanning a radius of 3 cells from the center of the network) representing efferent arborization, receives the “sensory input”, with “sensory” currents being delivered to the soma of modelled cells (gAMPA = 0. 15 mS/cm2). “Sensory input” was modeled according to O’Donnel et al. [122]. The mask is represented in Fig 5A. The cells of the inferior olivary network most likely share input sources due to overlapping arborizations of efferent projections [123]. To represent both shared and independent input, we have modeled the current source in each cell as having an independent process and a shared process, with a mix parameter (alpha) of 10% input correlation shared by all the cells in the network. This level of correlation leads to a coherent background oscillation in the cells of the network, which is exacerbated in the presence of gap junctions (S5 Fig). Ornstein-Uhlenbeck (OU) is a noise process that ensures that the mean current delivered is well behaved and that the integral of delivered current over time converges to a constant value [95]. The OU current is a good approximation for synaptic inputs originating in a large number of uncorrelated sources, where synaptic events are generated randomly and each event decays with a given rate (τ). We use a recursive implementation according to the following recursive formula: ηi (t+1) =ηi (n, t) *exp (−δ/τ) + (1/τ) (μ−ηi (t) ) +σ*√δ*ξi FORMULA 2 Where ηi (t) is the noise amplitude of neuron i at time t. The noise process is parameterized by τ, σ, μ where τ represents the synaptic decay time constant, δ is the integration step time for our forward Euler integrator, σ is the standard deviation of the noise process and μ is its mean. The random draw from a Gaussian distribution at every time step is represented by ξi. Neurons in the inferior olive receive broad arborizations, leading to input correlations across nearby neurons. In our model this is represented via a mixture of an independent process for each neuron nindependent and a shared process, nall, common to all the neurons in the network, parameterized by a mixing parameter α, called' noise correlation' : ni (t) =α*nindependent+ (1−α) *nall (t) FORMULA 3 Simulation results throughout the article come from simulations with noise use an α where neurons share 10% of their noise input. Reported results are qualitatively robust to changes in this value (S5 Fig). To examine the dependence of network dynamics on the characteristics of the incoming input, we computed the 200 neuron network sweeping a grid of the main input parameters (τ, σ, μ, α) of the Ornstein-Uhlenbeck noise process. The network response in terms of STO frequency, population firing rates, proportion of firing neurons was analyzed with respect to a grid of input parameters. For comparability of statistics and reproducibility of results, all results displayed in this article were obtained from a single random seed. We have tested the network with multiple seeds and the results are qualitatively indistinguishable. The parameters of the Ornstein-Uhlenbeck process were tuned such that the network emulating the wildtype network (with gap junctions, WT) produced an average frequency of 1 Hz and more than 95% of the model cells fire at least once every 5 seconds (the parameter space for the network responses including STO, population firing rate and proportion of cells that fire within 3s is found in S5B Fig). The parameters to achieve these criteria are dependent on the total leak through the gap junctions. There are multiple methods to compensate the absent leak in the gapless network. In the present case, the network without gap junctions has been tuned to produce the same firing frequency as the network with gaps by increasing the membrane leak currents from 0. 010 to 0. 013 mS/cm2. This results in a similar excitability but slightly lower STO frequency in the “mutant”. The average firing rate behavior of the network shows a linear relationship with the standard deviation of the OU process (S4 Fig). For the present network with balanced connectivity and a single gap conductance of 0. 04 mS/cm2, the Ornstein-Uhlenbeck parameters are (τ, σ, μ, α), μ = −0. 6 pA/cm2, σ = 0. 6 pA/cm2 and τ = 20 ms. τ is a decay parameter that represents the synaptic decay times expected for olivary inputs, in this case chosen to emulate dendritic GABA according to Devor and Yarom [124]. Both synchrony and instantaneous frequency were estimated on the basis of a novel phase transformation of the membrane potential, which is more robust than the standard Hilbert transform, and can produce a linear phase response to the non-linear shape of the subthreshold oscillations [125]. This transformation improves the estimation of the momentary phase and compensates for the fact that ionic mechanisms induce different rates of membrane potential change at different phases of the oscillation. This phases analysis was conducted with the DAMOCO toolbox [126]. From the instantaneous phase, the instantaneous frequency is simply the inverse of the first order finite difference of phase. Synchrony across cells is estimated with the Kuramoto order parameter (K): K (t) =|1N∑ei (Ψ (t) −ϕn (t) ) | FORMULA 4 Where ϕn is the phase of each neuron, N is the number of neurons and Ψ is the phase average of all oscillators. To estimate a phase response curve of the stimulated neurons, first a “sensory stimulus” is delivered at a phase known to produce an action potential (and resetting). The location of the first peak after stimulation is recorded. Subsequently, eight more simulations receive another stimulus, with same parameters as the resetting stimulus, but at different phases (at incremental intervals of 2π/8). The effect of that stimulation (delay or advance) on the next peak is recorded as a phase delta. Results are plotted in Fig 7A–7E.
Activity of the inferior olive, transmitted via climbing fibers to the cerebellum, regulates initiation and amplitude of movements, signals unexpected sensory feedback, and directs cerebellar learning. It is characterized by widespread subthreshold oscillations and synchronization promoted by strong electrotonic coupling. In brain slices, subthreshold oscillations gate which inputs can be transmitted by inferior olivary neurons and which will not—dependent on the phase of the oscillation. We tested whether the subthreshold oscillations had a measurable impact on temporal patterning of climbing fiber activity in intact, awake mice. We did so by recording neural activity of the postsynaptic Purkinje cells, in which complex spike firing faithfully represents climbing fiber activity. For short intervals (<300 ms) many Purkinje cells showed spontaneously rhythmic complex spike activity. However, our experiments designed to evoke conditional responses indicated that complex spikes are not predominantly predicated on stimulus history. Our realistic network model of the inferior olive explains the experimental observations via continuous phase modulations of the subthreshold oscillations under the influence of synaptic fluctuations. We conclude that complex spike activity emerges from a quasiperiodic rhythm that is stabilized by electrotonic coupling between its dendrites, yet dynamically influenced by the status of their synaptic inputs.
olives
cell physiology
action potentials
medicine and health sciences
neural networks
nervous system
membrane potential
junctional complexes
electrophysiology
neuroscience
gap junctions
network analysis
plants
computer and information sciences
animal cells
fruits
cellular neuroscience
eukaryota
cell biology
anatomy
synapses
physiology
neurons
biology and life sciences
cellular types
neurophysiology
purkinje cells
organisms
Quasiperiodic rhythms of the inferior olive