A spiking network model for passage-of-time representation in the cerebellum

Mesoscale simulations predict the role of synergistic cerebellar plasticity during classical eyeblink conditioning

According to the motor learning theory by Albus and Ito, synaptic depression at the parallel fibre to Purkinje cells synapse (pf-PC) is the main substrate responsible for learning sensorimotor contingencies under climbing fibre control. However, recent experimental evidence challenges this relatively monopolistic view of cerebellar learning. Bidirectional plasticity appears crucial for learning, in which different microzones can undergo opposite changes of synaptic strength (e.g. downbound microzones-more likely depression, upbound microzones-more likely potentiation), and multiple forms of plasticity have been identified, distributed over different cerebellar circuit synapses. Here, we have simulated classical eyeblink conditioning (CEBC) using an advanced spiking cerebellar model embedding downbound and upbound modules that are subject to multiple plasticity rules. Simulations indicate that synaptic plasticity regulates the cascade of precise spiking patterns spreading throughout the cerebellar cortex and cerebellar nuclei. CEBC was supported by plasticity at the pf-PC synapses as well as at the synapses of the molecular layer interneurons (MLIs), but only the combined switch-off of both sites of plasticity compromised learning significantly. By differentially engaging climbing fibre information and related forms of synaptic plasticity, both microzones contributed to generate a well-timed conditioned response, but it was the downbound module that played the major role in this process. The outcomes of our simulations closely align with the behavioural and electrophysiological phenotypes of mutant mice suffering from cell-specific mutations that affect processing of their PC and/or MLI synapses. Our data highlight that a synergy of bidirectional plasticity rules distributed across the cerebellum can facilitate finetuning of adaptive associative behaviours at a high spatiotemporal resolution.

Emergence of brain-inspired small-world spiking neural network through neuroevolution

Studies suggest that the brain's high efficiency and low energy consumption may be closely related to its small-world topology and critical dynamics. However, existing efforts on the performance-oriented structural evolution of spiking neural networks (SNNs) are time-consuming and ignore the core structural properties of the brain. Here, we introduce a multi-objective Evolutionary Liquid State Machine (ELSM), which blends the small-world coefficient and criticality to evolve models and guide the emergence of brain-inspired, efficient structures. Experiments reveal ELSM's consistent and comparable performance, achieving 97.23% on NMNIST and outperforming LSM models on MNIST and Fashion-MNIST with 98.12% and 88.81% accuracies, respectively. Further analysis shows its versatility and spontaneous evolution of topologies such as hub nodes, short paths, long-tailed degree distributions, and numerous communities. This study evolves recurrent spiking neural networks into brain-inspired energy-efficient structures, showcasing versatility in multiple tasks and potential for adaptive general artificial intelligence.

Positive feedback of efferent copy via pontine nucleus facilitates cerebellum-mediated associative learning

The cerebellum is critical for motor coordination and learning. However, the role of feedback circuitry in this brain region has not been fully explored. Here, we characterize a nucleo-ponto-cortical feedback pathway in classical delayed eyeblink conditioning (dEBC) of rats. We find that the efference copy is conveyed from the interposed cerebellar nucleus (Int) to cerebellar cortex through pontine nucleus (PN). Inhibiting or exciting the projection from the Int to the PN can decelerate or speed up acquisition of dEBC, respectively. Importantly, we identify two subpopulations of PN neurons (PN1 and PN2) that convey and integrate the feedback signals with feedforward sensory signals. We also show that the feedforward and feedback pathways via different types of PN neurons contribute to the plastic changes and cooperate synergistically to the learning of dEBC. Our results suggest that this excitatory nucleo-ponto-cortical feedback plays a significant role in modulating associative motor learning in cerebellum.

Synaptic basis of a sub-second representation of time in a neural circuit model

Temporal sequences of neural activity are essential for driving well-timed behaviors, but the underlying cellular and circuit mechanisms remain elusive. We leveraged the well-defined architecture of the cerebellum, a brain region known to support temporally precise actions, to explore theoretically whether the experimentally observed diversity of short-term synaptic plasticity (STP) at the input layer could generate neural dynamics sufficient for sub-second temporal learning. A cerebellar circuit model equipped with dynamic synapses produced a diverse set of transient granule cell firing patterns that provided a temporal basis set for learning precisely timed pauses in Purkinje cell activity during simulated delay eyelid conditioning and Bayesian interval estimation. The learning performance across time intervals was influenced by the temporal bandwidth of the temporal basis, which was determined by the input layer synaptic properties. The ubiquity of STP throughout the brain positions it as a general, tunable cellular mechanism for sculpting neural dynamics and fine-tuning behavior.

CerebelluMorphic: Large-Scale Neuromorphic Model and Architecture for Supervised Motor Learning

The cerebellum plays a vital role in motor learning and control with supervised learning capability, while neuromorphic engineering devises diverse approaches to high-performance computation inspired by biological neural systems. This article presents a large-scale cerebellar network model for supervised learning, as well as a cerebellum-inspired neuromorphic architecture to map the cerebellar anatomical structure into the large-scale model. Our multinucleus model and its underpinning architecture contain approximately 3.5 million neurons, upscaling state-of-the-art neuromorphic designs by over 34 times. Besides, the proposed model and architecture incorporate 3411k granule cells, introducing a 284 times increase compared to a previous study including only 12k cells. This large scaling induces more biologically plausible cerebellar divergence/convergence ratios, which results in better mimicking biology. In order to verify the functionality of our proposed model and demonstrate its strong biomimicry, a reconfigurable neuromorphic system is used, on which our developed architecture is realized to replicate cerebellar dynamics during the optokinetic response. In addition, our neuromorphic architecture is used to analyze the dynamical synchronization within the Purkinje cells, revealing the effects of firing rates of mossy fibers on the resonance dynamics of Purkinje cells. Our experiments show that real-time operation can be realized, with a system throughput of up to 4.70 times larger than previous works with high synaptic event rate. These results suggest that the proposed work provides both a theoretical basis and a neuromorphic engineering perspective for brain-inspired computing and the further exploration of cerebellar learning.

BiCoSS: Toward Large-Scale Cognition Brain With Multigranular Neuromorphic Architecture

The further exploration of the neural mechanisms underlying the biological activities of the human brain depends on the development of large-scale spiking neural networks (SNNs) with different categories at different levels, as well as the corresponding computing platforms. Neuromorphic engineering provides approaches to high-performance biologically plausible computational paradigms inspired by neural systems. In this article, we present a biological-inspired cognitive supercomputing system (BiCoSS) that integrates multiple granules (GRs) of SNNs to realize a hybrid compatible neuromorphic platform. A scalable hierarchical heterogeneous multicore architecture is presented, and a synergistic routing scheme for hybrid neural information is proposed. The BiCoSS system can accommodate different levels of GRs and biological plausibility of SNN models in an efficient and scalable manner. Over four million neurons can be realized on BiCoSS with a power efficiency of 2.8k larger than the GPU platform, and the average latency of BiCoSS is 3.62 and 2.49 times higher than conventional architectures of digital neuromorphic systems. For the verification, BiCoSS is used to replicate various biological cognitive activities, including motor learning, action selection, context-dependent learning, and movement disorders. Comprehensively considering the programmability, biological plausibility, learning capability, computational power, and scalability, BiCoSS is shown to outperform the alternative state-of-the-art works for large-scale SNN, while its real-time computational capability enables a wide range of potential applications.

Cerebellum Involvement in Dystonia During Associative Motor Learning: Insights From a Data-Driven Spiking Network Model

Dystonia is a movement disorder characterized by sustained or intermittent muscle contractions causing abnormal, often repetitive movements, postures, or both. Although dystonia is traditionally associated with basal ganglia dysfunction, recent evidence has been pointing to a role of the cerebellum, a brain area involved in motor control and learning. Cerebellar abnormalities have been correlated with dystonia but their potential causative role remains elusive. Here, we simulated the cerebellar input-output relationship with high-resolution computational modeling. We used a data-driven cerebellar Spiking Neural Network and simulated a cerebellum-driven associative learning task, Eye-Blink Classical Conditioning (EBCC), which is characteristically altered in relation to cerebellar lesions in several pathologies. In control simulations, input stimuli entrained characteristic network dynamics and induced synaptic plasticity along task repetitions, causing a progressive spike suppression in Purkinje cells with consequent facilitation of deep cerebellar nuclei cells. These neuronal processes caused a progressive acquisition of eyelid Conditioned Responses (CRs). Then, we modified structural or functional local neural features in the network reproducing alterations reported in dystonic mice. Either reduced olivocerebellar input or aberrant Purkinje cell burst-firing resulted in abnormal learning curves imitating the dysfunctional EBCC motor responses (in terms of CR amount and timing) of dystonic mice. These behavioral deficits might be due to altered temporal processing of sensorimotor information and uncoordinated control of muscle contractions. Conversely, an imbalance of excitatory and inhibitory synaptic densities on Purkinje cells did not reflect into significant EBCC deficit. The present work suggests that only certain types of alterations, including reduced olivocerebellar input and aberrant PC burst-firing, are compatible with the EBCC changes observed in dystonia, indicating that some cerebellar lesions can have a causative role in the pathogenesis of symptoms.

Influence of various temporal recoding on pavlovian eyeblink conditioning in the cerebellum

We consider the Pavlovian eyeblink conditioning (EBC) via repeated presentation of paired conditioned stimulus (tone) and unconditioned stimulus (US; airpuff). In an effective cerebellar ring network, we change the connection probability p c from Golgi to granule (GR) cells, and make a dynamical classification of various firing patterns of the GR cells. Individual GR cells are thus found to show various well- and ill-matched firing patterns relative to the US timing signal. Then, these variously-recoded signals are fed into the Purkinje cells (PCs) through the parallel-fibers (PFs). Based on such unique dynamical classification of various firing patterns, we make intensive investigations on the influence of various temporal recoding (i.e., firing patterns) of the GR cells on the synaptic plasticity of the PF-PC synapses and the subsequent learning process for the EBC. We first note that the variously-recoded PF signals are effectively depressed by the (error-teaching) instructor climbing-fiber (CF) signals from the inferior olive neuron. In the case of well-matched PF signals, they are strongly depressed through strong long-term depression (LTD) by the instructor CF signals due to good association between the in-phase PF and the instructor CF signals. On the other hand, practically no LTD occurs for the ill-matched PF signals because most of them have no association with the instructor CF signals. This kind of "effective" depression at the PF-PC synapses coordinates firings of PCs effectively, which then makes effective inhibitory coordination on the cerebellar nucleus neuron [which elicits conditioned response (CR; eyeblink)]. When the learning trial passes a threshold, acquisition of CR begins. In this case, the timing degree T d of CR becomes good due to presence of the ill-matched firing group which plays a role of protection barrier for the timing. With further increase in the number of trials, strength S of CR (corresponding to the amplitude of eyelid closure) increases due to strong LTD in the well-matched firing group, while its timing degree T d decreases. In this way, the well- and the ill-matched firing groups play their own roles for the strength and the timing of CR, respectively. Thus, with increasing the number of learning trials, the (overall) learning efficiency degree L e (taking into consideration both timing and strength of CR) for the CR is increased, and eventually it becomes saturated. Finally, we also discuss dependence of the variety degree for firing patterns of the GR cells and the saturated learning efficiency degree L e of the CR on p c and their relations.

A cerebellar-based solution to the nondeterministic time delay problem in robotic control

[Figure: see text].

Real-Time Simulation of a Cerebellar Scaffold Model on Graphics Processing Units

Large-scale simulation of detailed computational models of neuronal microcircuits plays a prominent role in reproducing and predicting the dynamics of the microcircuits. To reconstruct a microcircuit, one must choose neuron and synapse models, placements, connectivity, and numerical simulation methods according to anatomical and physiological constraints. For reconstruction and refinement, it is useful to be able to replace one module easily while leaving the others as they are. One way to achieve this is via a scaffolding approach, in which a simulation code is built on independent modules for placements, connections, and network simulations. Owing to the modularity of functions, this approach enables researchers to improve the performance of the entire simulation by simply replacing a problematic module with an improved one. Casali et al. (2019) developed a spiking network model of the cerebellar microcircuit using this approach, and while it reproduces electrophysiological properties of cerebellar neurons, it takes too much computational time. Here, we followed this scaffolding approach and replaced the simulation module with an accelerated version on graphics processing units (GPUs). Our cerebellar scaffold model ran roughly 100 times faster than the original version. In fact, our model is able to run faster than real time, with good weak and strong scaling properties. To demonstrate an application of real-time simulation, we implemented synaptic plasticity mechanisms at parallel fiber-Purkinje cell synapses, and carried out simulation of behavioral experiments known as gain adaptation of optokinetic response. We confirmed that the computer simulation reproduced experimental findings while being completed in real time. Actually, a computer simulation for 2 s of the biological time completed within 750 ms. These results suggest that the scaffolding approach is a promising concept for gradual development and refactoring of simulation codes for large-scale elaborate microcircuits. Moreover, a real-time version of the cerebellar scaffold model, which is enabled by parallel computing technology owing to GPUs, may be useful for large-scale simulations and engineering applications that require real-time signal processing and motor control.

Human-scale Brain Simulation via Supercomputer: A Case Study on the Cerebellum

Performance of supercomputers has been steadily and exponentially increasing for the past 20 years, and is expected to increase further. This unprecedented computational power enables us to build and simulate large-scale neural network models composed of tens of billions of neurons and tens of trillions of synapses with detailed anatomical connections and realistic physiological parameters. Such "human-scale" brain simulation could be considered a milestone in computational neuroscience and even in general neuroscience. Towards this milestone, it is mandatory to introduce modern high-performance computing technology into neuroscience research. In this article, we provide an introductory landscape about large-scale brain simulation on supercomputers from the viewpoints of computational neuroscience and modern high-performance computing technology for specialists in experimental as well as computational neurosciences. This introduction to modeling and simulation methods is followed by a review of various representative large-scale simulation studies conducted to date. Then, we direct our attention to the cerebellum, with a review of more simulation studies specific to that region. Furthermore, we present recent simulation results of a human-scale cerebellar network model composed of 86 billion neurons on the Japanese flagship supercomputer K (now retired). Finally, we discuss the necessity and importance of human-scale brain simulation, and suggest future directions of such large-scale brain simulation research.

Effect of diverse recoding of granule cells on optokinetic response in a cerebellar ring network with synaptic plasticity

We consider a cerebellar ring network for the optokinetic response (OKR), and investigate the effect of diverse recoding of granule (GR) cells on OKR by varying the connection probability p_c from Golgi to GR cells. For an optimal value of p_c^∗(=0.06), individual GR cells exhibit diverse spiking patterns which are in-phase, anti-phase, or complex out-of-phase with respect to their population-averaged firing activity. Then, these diversely-recoded signals via parallel fibers (PFs) from GR cells are effectively depressed by the error-teaching signals via climbing fibers from the inferior olive which are also in-phase ones. Synaptic weights at in-phase PF-Purkinje cell (PC) synapses of active GR cells are strongly depressed via strong long-term depression (LTD), while those at anti-phase and complex out-of-phase PF-PC synapses are weakly depressed through weak LTD. This kind of "effective" depression (i.e., strong/weak LTD) at the PF-PC synapses causes a big modulation in firings of PCs, which then exert effective inhibitory coordination on the vestibular nucleus (VN) neuron (which evokes OKR). For the firing of the VN neuron, the learning gain degree L_g, corresponding to the modulation gain ratio, increases with increasing the learning cycle, and it saturates at about the 300th cycle. By varying p_c from p_c^∗, we find that a plot of saturated learning gain degree L_g^∗ versus p_c forms a bell-shaped curve with a peak at p_c^∗ (where the diversity degree in spiking patterns of GR cells is also maximum). Consequently, the more diverse in recoding of GR cells, the more effective in motor learning for the OKR adaptation.

A Cerebellar Computational Mechanism for Delay Conditioning at Precise Time Intervals

The cerebellum is known to have an important role in sensing and execution of precise time intervals, but the mechanism by which arbitrary time intervals can be recognized and replicated with high precision is unknown. We propose a computational model in which precise time intervals can be identified from the pattern of individual spike activity in a population of parallel fibers in the cerebellar cortex. The model depends on the presence of repeatable sequences of spikes in response to conditioned stimulus input. We emulate granule cells using a population of Izhikevich neuron approximations driven by random but repeatable mossy fiber input. We emulate long-term depression (LTD) and long-term potentiation (LTP) synaptic plasticity at the parallel fiber to Purkinje cell synapse. We simulate a delay conditioning paradigm with a conditioned stimulus (CS) presented to the mossy fibers and an unconditioned stimulus (US) some time later issued to the Purkinje cells as a teaching signal. We show that Purkinje cells rapidly adapt to decrease firing probability following onset of the CS only at the interval for which the US had occurred. We suggest that detection of replicable spike patterns provides an accurate and easily learned timing structure that could be an important mechanism for behaviors that require identification and production of precise time intervals.

The Rules of Cerebellar Learning: Around the Ito Hypothesis

As a tribute to Masao Ito, we propose a model of cerebellar learning that incorporates and extends his original model. We suggest four principles that align well with conclusions from multiple cerebellar learning systems. (1) Climbing fiber inputs to the cerebellum drive early, fast, poorly-retained learning in the parallel fiber to Purkinje cell synapse. (2) Learned Purkinje cell outputs drive late, slow, well-retained learning in non-Purkinje cell inputs to neurons in the cerebellar nucleus, transferring learning from the cortex to the nucleus. (3) Recurrent feedback from Purkinje cells to the inferior olive, through interneurons in the cerebellar nucleus, limits the magnitude of fast, early learning in the cerebellar cortex. (4) Functionally different inputs are subjected to plasticity in the cerebellar cortex versus the cerebellar nucleus. A computational neural circuit model that is based on these principles mimics a large amount of neural and behavioral data obtained from the smooth pursuit eye movements of monkeys.

Simulation of a Human-Scale Cerebellar Network Model on the K Computer

Computer simulation of the human brain at an individual neuron resolution is an ultimate goal of computational neuroscience. The Japanese flagship supercomputer, K, provides unprecedented computational capability toward this goal. The cerebellum contains 80% of the neurons in the whole brain. Therefore, computer simulation of the human-scale cerebellum will be a challenge for modern supercomputers. In this study, we built a human-scale spiking network model of the cerebellum, composed of 68 billion spiking neurons, on the K computer. As a benchmark, we performed a computer simulation of a cerebellum-dependent eye movement task known as the optokinetic response. We succeeded in reproducing plausible neuronal activity patterns that are observed experimentally in animals. The model was built on dedicated neural network simulation software called MONET (Millefeuille-like Organization NEural neTwork), which calculates layered sheet types of neural networks with parallelization by tile partitioning. To examine the scalability of the MONET simulator, we repeatedly performed simulations while changing the number of compute nodes from 1,024 to 82,944 and measured the computational time. We observed a good weak-scaling property for our cerebellar network model. Using all 82,944 nodes, we succeeded in simulating a human-scale cerebellum for the first time, although the simulation was 578 times slower than the wall clock time. These results suggest that the K computer is already capable of creating a simulation of a human-scale cerebellar model with the aid of the MONET simulator.

Gradients in the mammalian cerebellar cortex enable Fourier-like transformation and improve storing capacity

Cerebellar granule cells (GCs) make up the majority of all neurons in the vertebrate brain, but heterogeneities among GCs and potential functional consequences are poorly understood. Here, we identified unexpected gradients in the biophysical properties of GCs in mice. GCs closer to the white matter (inner-zone GCs) had higher firing thresholds and could sustain firing with larger current inputs than GCs closer to the Purkinje cell layer (outer-zone GCs). Dynamic Clamp experiments showed that inner- and outer-zone GCs preferentially respond to high- and low-frequency mossy fiber inputs, respectively, enabling dispersion of the mossy fiber input into its frequency components as performed by a Fourier transformation. Furthermore, inner-zone GCs have faster axonal conduction velocity and elicit faster synaptic potentials in Purkinje cells. Neuronal network modeling revealed that these gradients improve spike-timing precision of Purkinje cells and decrease the number of GCs required to learn spike-sequences. Thus, our study uncovers biophysical gradients in the cerebellar cortex enabling a Fourier-like transformation of mossy fiber inputs.

Real-time cerebellar neuroprosthetic system based on a spiking neural network model of motor learning

OBJECTIVE

Damage to the brain, as a result of various medical conditions, impacts the everyday life of patients and there is still no complete cure to neurological disorders. Neuroprostheses that can functionally replace the damaged neural circuit have recently emerged as a possible solution to these problems. Here we describe the development of a real-time cerebellar neuroprosthetic system to substitute neural function in cerebellar circuitry for learning delay eyeblink conditioning (DEC).

APPROACH

The system was empowered by a biologically realistic spiking neural network (SNN) model of the cerebellar neural circuit, which considers the neuronal population and anatomical connectivity of the network. The model simulated synaptic plasticity critical for learning DEC. This SNN model was carefully implemented on a field programmable gate array (FPGA) platform for real-time simulation. This hardware system was interfaced in in vivo experiments with anesthetized rats and it used neural spikes recorded online from the animal to learn and trigger conditioned eyeblink in the animal during training.

MAIN RESULTS

This rat-FPGA hybrid system was able to process neuronal spikes in real-time with an embedded cerebellum model of ~10 000 neurons and reproduce learning of DEC with different inter-stimulus intervals. Our results validated that the system performance is physiologically relevant at both the neural (firing pattern) and behavioral (eyeblink pattern) levels.

SIGNIFICANCE

This integrated system provides the sufficient computation power for mimicking the cerebellar circuit in real-time. The system interacts with the biological system naturally at the spike level and can be generalized for including other neural components (neuron types and plasticity) and neural functions for potential neuroprosthetic applications.

Modeling the Encoding of Saccade Kinematic Metrics in the Purkinje Cell Layer of the Cerebellar Vermis

Recent electrophysiological observations related to saccadic eye movements in rhesus monkeys, suggest a prediction of the sensory consequences of movement in the Purkinje cell layer of the cerebellar oculomotor vermis (OMV). A definite encoding of real-time motion of the eye has been observed in simple-spike responses of the combined burst-pause Purkinje cell populations, organized based upon their complex-spike directional tuning. However, the underlying control mechanisms that could lead to such action encoding are still unclear. We propose a saccade control model, with emphasis on the structure of the OMV and its interaction with the extra-cerebellar components. In the simulated bilateral organization of the OMV, each caudal fastigial nucleus is arranged to receive incoming projections from combined burst-pause Purkinje cell populations. The OMV, through the caudal fastigial nuclei, interacts with the brainstem to provide adaptive saccade gain corrections that minimize the visual error in reaching a given target location. The simulation results corroborate the experimental Purkinje cell population activity patterns and their relation with saccade kinematic metrics. The Purkinje layer activity that emerges from the proposed organization, precisely predicted the speed of the eye at different target eccentricities. Simulated granular layer activity suggests no separate dynamics with respect to shaping the bilateral Purkine layer activity. We further examine the validity of the simulated OMV in maintaining the accuracy of saccadic eye movements in the presence of signal dependent variabilities, that can occur in extra-cerebellar pathways.

A Cerebellar Spiking Neural Model for Phase Reversal of Vestibulo-ocular Reflex

Cerebellum possesses very rich motor control and learning capability which is critical for animals. In this study, we proposed a spiking neural network model of cerebellum for gain and phase adaptation in vestibulo-ocular reflex (VOR). VOR is a critical adaptive reflexive eye movement for maintaining a stable visual field. In this model (with neuron number at the order of 104), synaptic plasticity at parallel fiber-Purkinje cell synapses was considered. In particular, we have shown that the inhibitory inputs from molecular layer interneurons on Purkinje cells play a critical role in phase adaptation of VOR. The inhibitory input from interneurons indirectly affects the strength of long-term potentiation (LTP) and long-term depression (LTD), resulting in more drastic phase shift upon learning and hence allowing phase reversal of VOR. The strength of inhibitory input also affects the maximum phase shift that can be achieved. Our result is consistent with experiments in mutant mice with blocked inhibitory inputs.

Conditioned climbing fiber responses in cerebellar cortex and nuclei

The eyeblink conditioning paradigm captures an elementary form of associative learning in a neural circuitry that is understood to an extraordinary degree. Cerebellar cortical Purkinje cell simple spike suppression is widely regarded as the main process underlying conditioned responses (CRs), leading to disinhibition of neurons in the cerebellar nuclei that innervate eyelid muscles downstream. However, recent work highlights the addition of a conditioned Purkinje cell complex spike response, which at the level of the interposed nucleus seems to translate to a transient spike suppression that can be followed by a rapid spike facilitation. Here, we review the characteristics of these responses at the cerebellar cortical and nuclear level, and discuss possible origins and functions.

Dynamic modulation of activity in cerebellar nuclei neurons during pavlovian eyeblink conditioning in mice

While research on the cerebellar cortex is crystallizing our understanding of its function in learning behavior, many questions surrounding its downstream targets remain. Here, we evaluate the dynamics of cerebellar interpositus nucleus (IpN) neurons over the course of Pavlovian eyeblink conditioning. A diverse range of learning-induced neuronal responses was observed, including increases and decreases in activity during the generation of conditioned blinks. Trial-by-trial correlational analysis and optogenetic manipulation demonstrate that facilitation in the IpN drives the eyelid movements. Adaptive facilitatory responses are often preceded by acquired transient inhibition of IpN activity that, based on latency and effect, appear to be driven by complex spikes in cerebellar cortical Purkinje cells. Likewise, during reflexive blinks to periocular stimulation, IpN cells show excitation-suppression patterns that suggest a contribution of climbing fibers and their collaterals. These findings highlight the integrative properties of subcortical neurons at the cerebellar output stage mediating conditioned behavior.

Spatiotemporal network coding of physiological mossy fiber inputs by the cerebellar granular layer

The granular layer, which mainly consists of granule and Golgi cells, is the first stage of the cerebellar cortex and processes spatiotemporal information transmitted by mossy fiber inputs with a wide variety of firing patterns. To study its dynamics at multiple time scales in response to inputs approximating real spatiotemporal patterns, we constructed a large-scale 3D network model of the granular layer. Patterned mossy fiber activity induces rhythmic Golgi cell activity that is synchronized by shared parallel fiber input and by gap junctions. This leads to long distance synchrony of Golgi cells along the transverse axis, powerfully regulating granule cell firing by imposing inhibition during a specific time window. The essential network mechanisms, including tunable Golgi cell oscillations, on-beam inhibition and NMDA receptors causing first winner keeps winning of granule cells, illustrate how fundamental properties of the granule layer operate in tandem to produce (1) well timed and spatially bound output, (2) a wide dynamic range of granule cell firing and (3) transient and coherent gating oscillations. These results substantially enrich our understanding of granule cell layer processing, which seems to promote spatial group selection of granule cell activity as a function of timing of mossy fiber input.

The cerebellum is an organ of peculiar geometrical properties, and has been attributed the function of applying spatiotemporal transforms to sensorimotor data since Eccles. In this work we have analyzed the spatiotemporal response properties of the first part of the cerebellar circuit, the granule layer. On the basis of a biophysically plausible and large-scale model of the cerebellum, constrained by a wealth of anatomical data, we study the network dynamics and firing properties of individual cell populations in response to 'realistic' input patterns. We make specific predictions about the spatiotemporal features of granule layer processing regarding the effects of the gap junction coupled network of Golgi cells on a spatially restricted input, in an effect we denominate first-takes-all. Furthermore, we calculate that the granule cell layer has a wide dynamic range, indicating that this is a system that can transmit large variations of input intensities.

Optogenetics in Silicon: A Neural Processor for Predicting Optically Active Neural Networks

We present a reconfigurable neural processor for real-time simulation and prediction of opto-neural behaviour. We combined a detailed Hodgkin-Huxley CA3 neuron integrated with a four-state Channelrhodopsin-2 (ChR2) model into reconfigurable silicon hardware. Our architecture consists of a Field Programmable Gated Array (FPGA) with a custom-built computing data-path, a separate data management system and a memory approach based router. Advancements over previous work include the incorporation of short and long-term calcium and light-dependent ion channels in reconfigurable hardware. Also, the developed processor is computationally efficient, requiring only 0.03 ms processing time per sub-frame for a single neuron and 9.7 ms for a fully connected network of 500 neurons with a given FPGA frequency of 56.7 MHz. It can therefore be utilized for exploration of closed loop processing and tuning of biologically realistic optogenetic circuitry.

Distributed Circuit Plasticity: New Clues for the Cerebellar Mechanisms of Learning

The cerebellum is involved in learning and memory of sensory motor skills. However, the way this process takes place in local microcircuits is still unclear. The initial proposal, casted into the Motor Learning Theory, suggested that learning had to occur at the parallel fiber-Purkinje cell synapse under supervision of climbing fibers. However, the uniqueness of this mechanism has been questioned, and multiple forms of long-term plasticity have been revealed at various locations in the cerebellar circuit, including synapses and neurons in the granular layer, molecular layer and deep-cerebellar nuclei. At present, more than 15 forms of plasticity have been reported. There has been a long debate on which plasticity is more relevant to specific aspects of learning, but this question turned out to be hard to answer using physiological analysis alone. Recent experiments and models making use of closed-loop robotic simulations are revealing a radically new view: one single form of plasticity is insufficient, while altogether, the different forms of plasticity can explain the multiplicity of properties characterizing cerebellar learning. These include multi-rate acquisition and extinction, reversibility, self-scalability, and generalization. Moreover, when the circuit embeds multiple forms of plasticity, it can easily cope with multiple behaviors endowing therefore the cerebellum with the properties needed to operate as an effective generalized forward controller.

Mechanisms and functional roles of glutamatergic synapse diversity in a cerebellar circuit

Synaptic currents display a large degree of heterogeneity of their temporal characteristics, but the functional role of such heterogeneities remains unknown. We investigated in rat cerebellar slices synaptic currents in Unipolar Brush Cells (UBCs), which generate intrinsic mossy fibers relaying vestibular inputs to the cerebellar cortex. We show that UBCs respond to sinusoidal modulations of their sensory input with heterogeneous amplitudes and phase shifts. Experiments and modeling indicate that this variability results both from the kinetics of synaptic glutamate transients and from the diversity of postsynaptic receptors. While phase inversion is produced by an mGluR2-activated outward conductance in OFF-UBCs, the phase delay of ON UBCs is caused by a late rebound current resulting from AMPAR recovery from desensitization. Granular layer network modeling indicates that phase dispersion of UBC responses generates diverse phase coding in the granule cell population, allowing climbing-fiber-driven Purkinje cell learning at arbitrary phases of the vestibular input.

DOI:http://dx.doi.org/10.7554/eLife.15872.001

Whether walking, riding a bicycle or simply standing still, we continually adjust our posture in small ways to prevent ourselves from falling. Our sense of balance depends on a set of structures inside the inner ear called the vestibular system. These structures detect movements of the head and relay this information to the brain in the form of electrical signals. A brain area called the vestibulo-cerebellum then combines these signals with sensory input from the eyes and muscles, before sending out further signals to trigger any adjustments necessary for balance.

One of the main cell types within the vestibulo-cerebellum is the unipolar brush cell (or UBC for short). UBCs pass on signals to another type of neuron called Purkinje cells, which support the learning of motor skills such as adjusting posture. Zampini, Liu et al. set out to test the idea that UBCs transform inputs from the vestibular system into a format that makes it easier for cerebellar Purkinje cells to drive this kind of learning.

First, recordings from slices of rodent brain revealed that UBCs respond in highly variable ways to vestibular input, with both the size and timing of responses varying between cells. This is because vestibular signals trigger the release of a chemical messenger called glutamate onto UBCs, but UBCs possess a variety of different types of glutamate receptors. Vestibular input therefore activates distinct signaling cascades from one UBC to the next. According to a computer model, this variability in UBC responses ensures that a subset of UBCs will always be active at any point during vestibular input. This in turn means that Purkinje cells can fire at any stage of a movement, which boosts the learning of motor skills.

The next steps will be to test this hypothesis using mutant mice that lack specific receptor subtypes in UBCs or UBCs completely. A further challenge for the future will be to build a computer model of the vestibulo-cerebellar system that includes all of its component cell types.

DOI:http://dx.doi.org/10.7554/eLife.15872.002

Real-Time Simulation of Passage-of-Time Encoding in Cerebellum Using a Scalable FPGA-Based System

The cerebellum plays a critical role for sensorimotor control and learning. However, dysmetria or delays in movements' onsets consequent to damages in cerebellum cannot be cured completely at the moment. Neuroprosthesis is an emerging technology that can potentially substitute such motor control module in the brain. A pre-requisite for this to become practical is the capability to simulate the cerebellum model in real-time, with low timing distortion for proper interfacing with the biological system. In this paper, we present a frame-based network-on-chip (NoC) hardware architecture for implementing a bio-realistic cerebellum model with  ∼ 100 000 neurons, which has been used for studying timing control or passage-of-time (POT) encoding mediated by the cerebellum. The simulation results verify that our implementation reproduces the POT representation by the cerebellum properly. Furthermore, our field-programmable gate array (FPGA)-based system demonstrates excellent computational speed that it can complete 1sec real world activities within 25.6 ms. It is also highly scalable such that it can maintain approximately the same computational speed even if the neuron number increases by one order of magnitude. Our design is shown to outperform three alternative approaches previously used for implementing spiking neural network model. Finally, we show a hardware electronic setup and illustrate how the silicon cerebellum can be adapted as a potential neuroprosthetic platform for future biological or clinical application.

Real-World-Time Simulation of Memory Consolidation in a Large-Scale Cerebellar Model

We report development of a large-scale spiking network model of the cerebellum composed of more than 1 million neurons. The model is implemented on graphics processing units (GPUs), which are dedicated hardware for parallel computing. Using 4 GPUs simultaneously, we achieve realtime simulation, in which computer simulation of cerebellar activity for 1 s completes within 1 s in the real-world time, with temporal resolution of 1 ms. This allows us to carry out a very long-term computer simulation of cerebellar activity in a practical time with millisecond temporal resolution. Using the model, we carry out computer simulation of long-term gain adaptation of optokinetic response (OKR) eye movements for 5 days aimed to study the neural mechanisms of posttraining memory consolidation. The simulation results are consistent with animal experiments and our theory of posttraining memory consolidation. These results suggest that realtime computing provides a useful means to study a very slow neural process such as memory consolidation in the brain.

Distributed Cerebellar Motor Learning: A Spike-Timing-Dependent Plasticity Model

Deep cerebellar nuclei neurons receive both inhibitory (GABAergic) synaptic currents from Purkinje cells (within the cerebellar cortex) and excitatory (glutamatergic) synaptic currents from mossy fibers. Those two deep cerebellar nucleus inputs are thought to be also adaptive, embedding interesting properties in the framework of accurate movements. We show that distributed spike-timing-dependent plasticity mechanisms (STDP) located at different cerebellar sites (parallel fibers to Purkinje cells, mossy fibers to deep cerebellar nucleus cells, and Purkinje cells to deep cerebellar nucleus cells) in close-loop simulations provide an explanation for the complex learning properties of the cerebellum in motor learning. Concretely, we propose a new mechanistic cerebellar spiking model. In this new model, deep cerebellar nuclei embed a dual functionality: deep cerebellar nuclei acting as a gain adaptation mechanism and as a facilitator for the slow memory consolidation at mossy fibers to deep cerebellar nucleus synapses. Equipping the cerebellum with excitatory (e-STDP) and inhibitory (i-STDP) mechanisms at deep cerebellar nuclei afferents allows the accommodation of synaptic memories that were formed at parallel fibers to Purkinje cells synapses and then transferred to mossy fibers to deep cerebellar nucleus synapses. These adaptive mechanisms also contribute to modulate the deep-cerebellar-nucleus-output firing rate (output gain modulation toward optimizing its working range).

At the Edge of Chaos: How Cerebellar Granular Layer Network Dynamics Can Provide the Basis for Temporal Filters

Models of the cerebellar microcircuit often assume that input signals from the mossy-fibers are expanded and recoded to provide a foundation from which the Purkinje cells can synthesize output filters to implement specific input-signal transformations. Details of this process are however unclear. While previous work has shown that recurrent granule cell inhibition could in principle generate a wide variety of random outputs suitable for coding signal onsets, the more general application for temporally varying signals has yet to be demonstrated. Here we show for the first time that using a mechanism very similar to reservoir computing enables random neuronal networks in the granule cell layer to provide the necessary signal separation and extension from which Purkinje cells could construct basis filters of various time-constants. The main requirement for this is that the network operates in a state of criticality close to the edge of random chaotic behavior. We further show that the lack of recurrent excitation in the granular layer as commonly required in traditional reservoir networks can be circumvented by considering other inherent granular layer features such as inverted input signals or mGluR2 inhibition of Golgi cells. Other properties that facilitate filter construction are direct mossy fiber excitation of Golgi cells, variability of synaptic weights or input signals and output-feedback via the nucleocortical pathway. Our findings are well supported by previous experimental and theoretical work and will help to bridge the gap between system-level models and detailed models of the granular layer network.

The 40-year history of modeling active dendrites in cerebellar Purkinje cells: emergence of the first single cell "community model"

The subject of the effects of the active properties of the Purkinje cell dendrite on neuronal function has been an active subject of study for more than 40 years. Somewhat unusually, some of these investigations, from the outset have involved an interacting combination of experimental and model-based techniques. This article recounts that 40-year history, and the view of the functional significance of the active properties of the Purkinje cell dendrite that has emerged. It specifically considers the emergence from these efforts of what is arguably the first single cell "community" model in neuroscience. The article also considers the implications of the development of this model for future studies of the complex properties of neuronal dendrites.

Information processing in the hemisphere of the cerebellar cortex for control of wrist movement

A region of cerebellar lobules V and VI makes strong loop connections with the primary motor (M1) and premotor (PM) cortical areas and is assumed to play essential roles in limb motor control. To examine its functional role, we compared the activities of its input, intermediate, and output elements, i.e., mossy fibers (MFs), Golgi cells (GoCs), and Purkinje cells (PCs), in three monkeys performing wrist movements in two different forearm postures. The results revealed distinct steps of information processing. First, MF activities displayed temporal and directional properties that were remarkably similar to those of M1/PM neurons, suggesting that MFs relay near copies of outputs from these motor areas. Second, all GoCs had a stereotyped pattern of activity independent of movement direction or forearm posture. Instead, GoC activity resembled an average of all MF activities. Therefore, inhibitory GoCs appear to provide a filtering function that passes only prominently modulated MF inputs to granule cells. Third, PCs displayed highly complex spatiotemporal patterns of activity, with coordinate frames distinct from those of MF inputs and directional tuning that changed abruptly before movement onset. The complexity of PC activities may reflect rapidly changing properties of the peripheral motor apparatus during movement. Overall, the cerebellar cortex appears to transform a representation of outputs from M1/PM into different movement representations in a posture-dependent manner and could work as part of a forward model that predicts the state of the peripheral motor apparatus.

Modeling possible effects of atypical cerebellar processing on eyeblink conditioning in autism

Autism is unique among other disorders in that acquisition of conditioned eyeblink responses is enhanced in children, occurring in a fraction of the trials required for control participants. The timing of learned responses is, however, atypical. Two animal models of autism display a similar phenotype. Researchers have hypothesized that these differences in conditioning reflect cerebellar abnormalities. The present study used computer simulations of the cerebellar cortex, including inhibition by the molecular layer interneurons, to more closely examine whether atypical cerebellar processing can account for faster conditioning in individuals with autism. In particular, the effects of inhibitory levels on delay eyeblink conditioning were simulated, as were the effects of learning-related synaptic changes at either parallel fibers or ascending branch synapses from granule cells to Purkinje cells. Results from these simulations predict that whether molecular layer inhibition results in an enhancement or an impairment of acquisition, or changes in timing, may depend on (1) the sources of inhibition, (2) the levels of inhibition, and (3) the locations of learning-related changes (parallel vs. ascending branch synapses). Overall, the simulations predict that a disruption in the balance or an overall increase of inhibition within the cerebellar cortex may contribute to atypical eyeblink conditioning in children with autism and in animal models of autism.

A scalable FPGA-based cerebellum for passage-of-time representation

The cerebellum plays a critical role for sensorimotor control and learning. However dysmertria or delays in movements' onsets consequent to damages in cerebellum cannot be cured completely at the moment. To foster a potential cure based on neuroprosthetic technology, we present a frame-based Network-on-Chip (NoC) hardware architecture for implementing a bio-realistic cerebellum model with 100,000 neurons, which has been used for studying timing control or passage-of-time (POT) encoding mediated by the cerebellum. The results demonstrate that our implementation can reproduce the POT functionality properly. The computational speed can achieve to 25.6 ms for simulating 1 sec real world activities. Furthermore, we show a hardware electronic setup and illustrate how the silicon cerebellum can be adapted as a potential neuroprosthetic platform for future biological or clinical applications.

Fast convergence of learning requires plasticity between inferior olive and deep cerebellar nuclei in a manipulation task: a closed-loop robotic simulation

The cerebellum is known to play a critical role in learning relevant patterns of activity for adaptive motor control, but the underlying network mechanisms are only partly understood. The classical long-term synaptic plasticity between parallel fibers (PFs) and Purkinje cells (PCs), which is driven by the inferior olive (IO), can only account for limited aspects of learning. Recently, the role of additional forms of plasticity in the granular layer, molecular layer and deep cerebellar nuclei (DCN) has been considered. In particular, learning at DCN synapses allows for generalization, but convergence to a stable state requires hundreds of repetitions. In this paper we have explored the putative role of the IO-DCN connection by endowing it with adaptable weights and exploring its implications in a closed-loop robotic manipulation task. Our results show that IO-DCN plasticity accelerates convergence of learning by up to two orders of magnitude without conflicting with the generalization properties conferred by DCN plasticity. Thus, this model suggests that multiple distributed learning mechanisms provide a key for explaining the complex properties of procedural learning and open up new experimental questions for synaptic plasticity in the cerebellar network.

Long-term depression as a model of cerebellar plasticity

Long-term depression (LTD) here concerned is persistent attenuation of transmission efficiency from a bundle of parallel fibers to a Purkinje cell. Uniquely, LTD is induced by conjunctive activation of the parallel fibers and the climbing fiber that innervates that Purkinje cell. Cellular and molecular processes underlying LTD occur postsynaptically. In the 1960s, LTD was conceived as a theoretical possibility and in the 1980s, substantiated experimentally. Through further investigations using various pharmacological or genetic manipulations of LTD, a concept was formed that LTD plays a major role in learning capability of the cerebellum (referred to as "Marr-Albus-Ito hypothesis"). In this chapter, following a historical overview, recent intensive investigations of LTD are reviewed. Complex signal transduction and receptor recycling processes underlying LTD are analyzed, and roles of LTD in reflexes and voluntary movements are defined. The significance of LTD is considered from viewpoints of neural network modeling. Finally, the controversy arising from the recent finding in a few studies that whereas LTD is blocked pharmacologically or genetically, motor learning in awake behaving animals remains seemingly unchanged is examined. We conjecture how this mismatch arises, either from a methodological problem or from a network nature, and how it might be resolved.

Granule cell ascending axon excitatory synapses onto Golgi cells implement a potent feedback circuit in the cerebellar granular layer

The function of inhibitory interneurons within brain microcircuits depends critically on the nature and properties of their excitatory synaptic drive. Golgi cells (GoCs) of the cerebellum inhibit cerebellar granule cells (GrCs) and are driven both by feedforward mossy fiber (mf) and feedback GrC excitation. Here, we have characterized GrC inputs to GoCs in rats and mice. We show that, during sustained mf discharge, synapses from local GrCs contribute equivalent charge to GoCs as mf synapses, arguing for the importance of the feedback inhibition. Previous studies predicted that GrC-GoC synapses occur predominantly between parallel fibers (pfs) and apical GoC dendrites in the molecular layer (ML). By combining EM and Ca(2+) imaging, we now demonstrate the presence of functional synaptic contacts between ascending axons (aa) of GrCs and basolateral dendrites of GoCs in the granular layer (GL). Immunohistochemical quantification estimates these contacts to be ∼400 per GoC. Using Ca(2+) imaging to identify synaptic inputs, we show that EPSCs from aa and mf contacts in basolateral dendrites display similarly fast kinetics, whereas pf inputs in the ML exhibit markedly slower kinetics as they undergo strong filtering by apical dendrites. We estimate that approximately half of the local GrC contacts generate fast EPSCs, indicating their basolateral location in the GL. We conclude that GrCs, through their aa contacts onto proximal GoC dendrites, define a powerful feedback inhibitory circuit in the GL.

Golgi cell activity during eyeblink conditioning in decerebrate ferrets

Golgi cells have a central position in the cerebellar cortical network and are indirectly connected to Purkinje cells, which are important for the acquisition of learned responses in classical conditioning. In order to clarify the role of Golgi cells in classical conditioning, we made extracellular Golgi cell recordings during different stages of conditioning, using four different conditional stimuli. Our results show that forelimb and superior colliculus stimulation, but not mossy fiber stimulation, evokes a short latency increase in Golgi cell firing. These results suggest that Golgi cells are involved in modulating input to the cerebellar cortex. There were however no differences in Golgi cell activity between naïve and trained animals, which suggests that Golgi cells are not intimately involved in the plastic changes that occur during classical conditioning. The absence of long latency effects of the conditional stimulus also questions whether Golgi cells contribute to the generation of a temporal code in the granule cells.

Perceptual decision making "through the eyes" of a large-scale neural model of v1

Sparse coding has been posited as an efficient information processing strategy employed by sensory systems, particularly visual cortex. Substantial theoretical and experimental work has focused on the issue of sparse encoding, namely how the early visual system maps the scene into a sparse representation. In this paper we investigate the complementary issue of sparse decoding, for example given activity generated by a realistic mapping of the visual scene to neuronal spike trains, how do downstream neurons best utilize this representation to generate a “decision.” Specifically we consider both sparse (L1-regularized) and non-sparse (L2 regularized) linear decoding for mapping the neural dynamics of a large-scale spiking neuron model of primary visual cortex (V1) to a two alternative forced choice (2-AFC) perceptual decision. We show that while both sparse and non-sparse linear decoding yield discrimination results quantitatively consistent with human psychophysics, sparse linear decoding is more efficient in terms of the number of selected informative dimension.

Realtime cerebellum: a large-scale spiking network model of the cerebellum that runs in realtime using a graphics processing unit

The cerebellum plays an essential role in adaptive motor control. Once we are able to build a cerebellar model that runs in realtime, which means that a computer simulation of 1 s in the simulated world completes within 1 s in the real world, the cerebellar model could be used as a realtime adaptive neural controller for physical hardware such as humanoid robots. In this paper, we introduce "Realtime Cerebellum (RC)", a new implementation of our large-scale spiking network model of the cerebellum, which was originally built to study cerebellar mechanisms for simultaneous gain and timing control and acted as a general-purpose supervised learning machine of spatiotemporal information known as reservoir computing, on a graphics processing unit (GPU). Owing to the massive parallel computing capability of a GPU, RC runs in realtime, while reproducing qualitatively the same simulation results of the Pavlovian delay eyeblink conditioning with the previous version. RC is adopted as a realtime adaptive controller of a humanoid robot, which is instructed to learn a proper timing to swing a bat to hit a flying ball online. These results suggest that RC provides a means to apply the computational power of the cerebellum as a versatile supervised learning machine towards engineering applications.

Modeling spike-train processing in the cerebellum granular layer and changes in plasticity reveal single neuron effects in neural ensembles

The cerebellum input stage has been known to perform combinatorial operations on input signals. In this paper, two types of mathematical models were used to reproduce the role of feed-forward inhibition and computation in the granular layer microcircuitry to investigate spike train processing. A simple spiking model and a biophysically-detailed model of the network were used to study signal recoding in the granular layer and to test observations like center-surround organization and time-window hypothesis in addition to effects of induced plasticity. Simulations suggest that simple neuron models may be used to abstract timing phenomenon in large networks, however detailed models were needed to reconstruct population coding via evoked local field potentials (LFP) and for simulating changes in synaptic plasticity. Our results also indicated that spatio-temporal code of the granular network is mainly controlled by the feed-forward inhibition from the Golgi cell synapses. Spike amplitude and total number of spikes were modulated by LTP and LTD. Reconstructing granular layer evoked-LFP suggests that granular layer propagates the nonlinearities of individual neurons. Simulations indicate that granular layer network operates a robust population code for a wide range of intervals, controlled by the Golgi cell inhibition and is regulated by the post-synaptic excitability.

NMDA receptors with incomplete Mg²⁺ block enable low-frequency transmission through the cerebellar cortex

The cerebellar cortex coordinates movements and maintains balance by modifying motor commands as a function of sensory-motor context, which is encoded by mossy fiber (MF) activity. MFs exhibit a wide range of activity, from brief precisely timed high-frequency bursts, which encode discrete variables such as whisker stimulation, to low-frequency sustained rate-coded modulation, which encodes continuous variables such as head velocity. While high-frequency MF inputs have been shown to activate granule cells (GCs) effectively, much less is known about sustained low-frequency signaling through the GC layer, which is impeded by a hyperpolarized resting potential and strong GABA(A)-mediated tonic inhibition of GCs. Here we have exploited the intrinsic MF network of unipolar brush cells to activate GCs with sustained low-frequency asynchronous MF inputs in rat cerebellar slices. We find that low-frequency MF input modulates the intrinsic firing of Purkinje cells, and that this signal transmission through the GC layer requires synaptic activation of Mg²⁺-block-resistant NMDA receptors (NMDARs) that are likely to contain the GluN2C subunit. Slow NMDAR conductances sum temporally to contribute approximately half the MF-GC synaptic charge at hyperpolarized potentials. Simulations of synaptic integration in GCs show that the NMDAR and slow spillover-activated AMPA receptor (AMPAR) components depolarize GCs to a similar extent. Moreover, their combined depolarizing effect enables the fast quantal AMPAR component to trigger action potentials at low MF input frequencies. Our results suggest that the weak Mg²⁺ block of GluN2C-containing NMDARs enables transmission of low-frequency MF signals through the input layer of the cerebellar cortex.

A computational mechanism for unified gain and timing control in the cerebellum

Precise gain and timing control is the goal of cerebellar motor learning. Because the basic neural circuitry of the cerebellum is homogeneous throughout the cerebellar cortex, a single computational mechanism may be used for simultaneous gain and timing control. Although many computational models of the cerebellum have been proposed for either gain or timing control, few models have aimed to unify them. In this paper, we hypothesize that gain and timing control can be unified by learning of the complete waveform of the desired movement profile instructed by climbing fiber signals. To justify our hypothesis, we adopted a large-scale spiking network model of the cerebellum, which was originally developed for cerebellar timing mechanisms to explain the experimental data of Pavlovian delay eyeblink conditioning, to the gain adaptation of optokinetic response (OKR) eye movements. By conducting large-scale computer simulations, we could reproduce some features of OKR adaptation, such as the learning-related change of simple spike firing of model Purkinje cells and vestibular nuclear neurons, simulated gain increase, and frequency-dependent gain increase. These results suggest that the cerebellum may use a single computational mechanism to control gain and timing simultaneously.

Learning the optimal control of coordinated eye and head movements

Various optimality principles have been proposed to explain the characteristics of coordinated eye and head movements during visual orienting behavior. At the same time, researchers have suggested several neural models to underly the generation of saccades, but these do not include online learning as a mechanism of optimization. Here, we suggest an open-loop neural controller with a local adaptation mechanism that minimizes a proposed cost function. Simulations show that the characteristics of coordinated eye and head movements generated by this model match the experimental data in many aspects, including the relationship between amplitude, duration and peak velocity in head-restrained and the relative contribution of eye and head to the total gaze shift in head-free conditions. Our model is a first step towards bringing together an optimality principle and an incremental local learning mechanism into a unified control scheme for coordinated eye and head movements.

Human beings and many other species redirect their gaze towards targets of interest through rapid gaze shifts known as saccades. These are made approximately three to four times every second, and larger saccades result from fast and concurrent movement of the animal's eyes and head. Experimental studies have revealed that during saccades, the motor system follows certain principles such as respecting a specific relationship between the relative contribution of eye and head motor systems to total gaze shift. Various researchers have hypothesized that these principles are implications of some optimality criteria in the brain, but it remains unclear how the brain can learn such an optimal behavior. We propose a new model that uses a plausible learning mechanism to satisfy an optimality criterion. We show that after learning, the model is able to reproduce motor behavior with biologically plausible properties. In addition, it predicts the nature of the learning signals. Further experimental research is necessary to test the validity of our model.

Stimulus-dependent state transition between synchronized oscillation and randomly repetitive burst in a model cerebellar granular layer

Information processing of the cerebellar granular layer composed of granule and Golgi cells is regarded as an important first step toward the cerebellar computation. Our previous theoretical studies have shown that granule cells can exhibit random alternation between burst and silent modes, which provides a basis of population representation of the passage-of-time (POT) from the onset of external input stimuli. On the other hand, another computational study has reported that granule cells can exhibit synchronized oscillation of activity, as consistent with observed oscillation in local field potential recorded from the granular layer while animals keep still. Here we have a question of whether an identical network model can explain these distinct dynamics. In the present study, we carried out computer simulations based on a spiking network model of the granular layer varying two parameters: the strength of a current injected to granule cells and the concentration of Mg²⁺ which controls the conductance of NMDA channels assumed on the Golgi cell dendrites. The simulations showed that cells in the granular layer can switch activity states between synchronized oscillation and random burst-silent alternation depending on the two parameters. For higher Mg²⁺ concentration and a weaker injected current, granule and Golgi cells elicited spikes synchronously (synchronized oscillation state). In contrast, for lower Mg²⁺ concentration and a stronger injected current, those cells showed the random burst-silent alternation (POT-representing state). It is suggested that NMDA channels on the Golgi cell dendrites play an important role for determining how the granular layer works in response to external input.

Intensive studies of Pavlovian delay eyelid conditioning suggest that the cerebellum can memorize a passage-of-time (POT) from the onset of an external stimulus. To account for possible mechanisms of such POT representation, some network models have been proposed to show that granule cells (grcs) in the cerebellar granular layer can exhibit random alternation of burst and silent modes under feedback inhibition from Golgi cells, resulting in non-recurrent generation of active granule cells populations. On the other hand, the oscillation of local field potential (LFP) has been observed in the cerebellar granular layer when animals stay at rest. Some network models have shown that grcs can elicit synchronous spikes in an oscillatory manner. These qualitatively different neural dynamics of the granular layer raises a question of how they can be accounted for by an identical network in the granular layer. Here we report that grc activities of a biologically plausible spiking network model undergo the state transition between synchronized oscillation and random burst-silent alternation, depending on the activation of NMDA channels on the Golgi cell dendrites and the strength of a current injected to grcs.

Sensory prediction or motor control? Application of marr-albus type models of cerebellar function to classical conditioning

Marr-Albus adaptive filter models of the cerebellum have been applied successfully to a range of sensory and motor control problems. Here we analyze their properties when applied to classical conditioning of the nictitating membrane response in rabbits. We consider a system-level model of eyeblink conditioning based on the anatomy of the eyeblink circuitry, comprising an adaptive filter model of the cerebellum, a comparator model of the inferior olive and a linear dynamic model of the nictitating membrane plant. To our knowledge, this is the first model that explicitly includes all these principal components, in particular the motor plant that is vital for shaping and timing the behavioral response. Model assumptions and parameters were systematically investigated to disambiguate basic computational capacities of the model from features requiring tuning of properties and parameter values. Without such tuning, the model robustly reproduced a range of behaviors related to sensory prediction, by displaying appropriate trial-level associative learning effects for both single and multiple stimuli, including blocking and conditioned inhibition. In contrast, successful reproduction of the real-time motor behavior depended on appropriate specification of the plant, cerebellum and comparator models. Although some of these properties appear consistent with the system biology, fundamental questions remain about how the biological parameters are chosen if the cerebellar microcircuit applies a common computation to many distinct behavioral tasks. It is possible that the response profiles in classical conditioning of the eyeblink depend upon operant contingencies that have previously prevailed, for example in naturally occurring avoidance movements.

Model-founded explorations of the roles of molecular layer inhibition in regulating purkinje cell responses in cerebellar cortex: more trouble for the beam hypothesis

For most of the last 50 years, the functional interpretation for inhibition in cerebellar cortical circuitry has been dominated by the relatively simple notion that excitatory and inhibitory dendritic inputs sum, and if that sum crosses threshold at the soma the Purkinje cell generates an action potential. Thus, inhibition has traditionally been relegated to a role of sculpting, restricting, or blocking excitation. At the level of networks, this relatively simply notion is manifest in mechanisms like "surround inhibition" which is purported to "shape" or "tune" excitatory neuronal responses. In the cerebellum, where all cell types except one (the granule cell) are inhibitory, these assumptions regarding the role of inhibition continue to dominate. Based on our recent series of modeling and experimental studies, we now suspect that inhibition may play a much more complex, subtle, and central role in the physiological and functional organization of cerebellar cortex. This paper outlines how model-based studies are changing our thinking about the role of feed-forward molecular layer inhibition in the cerebellar cortex. The results not only have important implications for continuing efforts to understand what the cerebellum computes, but might also reveal important features of the evolution of this large and quintessentially vertebrate brain structure.

Computational models of timing mechanisms in the cerebellar granular layer

A long-standing question in neuroscience is how the brain controls movement that requires precisely timed muscle activations. Studies using Pavlovian delay eyeblink conditioning provide good insight into this question. In delay eyeblink conditioning, which is believed to involve the cerebellum, a subject learns an interstimulus interval (ISI) between the onsets of a conditioned stimulus (CS) such as a tone and an unconditioned stimulus such as an airpuff to the eye. After a conditioning phase, the subject’s eyes automatically close or blink when the ISI time has passed after CS onset. This timing information is thought to be represented in some way in the cerebellum. Several computational models of the cerebellum have been proposed to explain the mechanisms of time representation, and they commonly point to the granular layer network. This article will review these computational models and discuss the possible computational power of the cerebellum.

Timing and plasticity in the cerebellum: focus on the granular layer

Two of the most striking properties of the cerebellum are its control in timing of motor operations and its ability to adapt behavior to new sensorimotor associations. Here, we propose a 'time-window matching' hypothesis for granular layer processing. Our hypothesis states that mossy fiber inputs to the granular layer are transformed into well-timed spike bursts by intrinsic granule cell processing, that feedforward Golgi cell inhibition sets a limit to the duration of such bursts and that these activities are spread over particular fields in the granular layer so as to generate ongoing time-windows for proper control of interacting motor domains. The role of synaptic plasticity would be that of fine-tuning pre-wired circuits favoring activation of specific granule cell groups in relation to particular time windows. This concept has wide implications for processing in the olivo-cerebellar system as a whole.