Theta-frequency bursting and resonance in cerebellar granule cells: experimental evidence and modeling of a slow k+-dependent mechanism

Theta-tACS modulates cerebellar-related motor functions and cerebellar-cortical connectivity

OBJECTIVE

To evaluate the effects of cerebellar transcranial alternating current stimulation (tACS) delivered at cerebellar-resonant frequencies, i.e., theta (θ) and gamma (γ), on upper limb motor performance and cerebellum-primary motor cortex (M1) connectivity, as assessed by cerebellar-brain inhibition (CBI), in healthy subjects.

METHODS

Participants underwent cerebellar-tACS while performing three cerebellar-dependent motor tasks: (i) rhythmic finger-tapping, (ii) arm reaching-to-grasp ('grasping') and (iii) arm reaching-to-point ('pointing') an object. Also, we evaluated possible changes in CBI during cerebellar-tACS.

RESULTS

θ-tACS decreased movement regularity during the tapping task and increased the duration of the pointing task compared to sham- and γ-tACS. Additionally, θ-tACS increased the CBI effectiveness (greater inhibition). The effect of θ-tACS on movement rhythm correlated with CBI changes and less tapping regularity corresponded to greater CBI.

CONCLUSIONS

Cerebellar-tACS delivered at the θ frequency modulates cerebellar-related motor behavior and this effect is, at least in part, mediated by changes in the cerebellar inhibitory output onto M1. The effects of θ-tACS may be due to the modulation of cerebellar neurons that resonate to the θ rhythm.

SIGNIFICANCE

These findings contribute to a better understanding of the physiological mechanisms of motor control and provide new evidence on cerebellar non-invasive brain stimulation.

A system of feed-forward cerebellar circuits that extend and diversify sensory signaling

Sensory signals are processed by the cerebellum to coordinate movements. Numerous cerebellar functions are thought to require the maintenance of a sensory representation that extends beyond the input signal. Granule cells receive sensory input, but they do not prolong the signal and are thus unlikely to maintain a sensory representation for much longer than the inputs themselves. Unipolar brush cells (UBCs) are excitatory interneurons that project to granule cells and transform sensory input into prolonged increases or decreases in firing, depending on their ON or OFF UBC subtype. Further extension and diversification of the input signal could be produced by UBCs that project to one another, but whether this circuitry exists is unclear. Here we test whether UBCs innervate one another and explore how these small networks of UBCs could transform spiking patterns. We characterized two transgenic mouse lines electrophysiologically and immunohistochemically to confirm that they label ON and OFF UBC subtypes and crossed them together, revealing that ON and OFF UBCs innervate one another. A Brainbow reporter was used to label UBCs of the same ON or OFF subtype with different fluorescent proteins, which showed that UBCs innervate their own subtypes as well. Computational models predict that these feed-forward networks of UBCs extend the length of bursts or pauses and introduce delays-transformations that may be necessary for cerebellar functions from modulation of eye movements to adaptive learning across time scales.

A system of feed-forward cerebellar circuits that extend and diversify sensory signaling

Sensory signals are processed by the cerebellum to coordinate movements. Numerous cerebellar functions are thought to require the maintenance of a sensory representation that extends beyond the input signal. Granule cells receive sensory input, but they do not prolong the signal and are thus unlikely to maintain a sensory representation for much longer than the inputs themselves. Unipolar brush cells (UBCs) are excitatory interneurons that project to granule cells and transform sensory input into prolonged increases or decreases in firing, depending on their ON or OFF UBC subtype. Further extension and diversification of the input signal could be produced by UBCs that project to one another, but whether this circuitry exists is unclear. Here we test whether UBCs innervate one another and explore how these small networks of UBCs could transform spiking patterns. We characterized two transgenic mouse lines electrophysiologically and immunohistochemically to confirm that they label ON and OFF UBC subtypes and crossed them together, revealing that ON and OFF UBCs innervate one another. A Brainbow reporter was used to label UBCs of the same ON or OFF subtype with different fluorescent proteins, which showed that UBCs innervate their own subtypes as well. Computational models predict that these feed-forward networks of UBCs extend the length of bursts or pauses and introduce delays-transformations that may be necessary for cerebellar functions from modulation of eye movements to adaptive learning across time scales.

Simulations predict differing phase responses to excitation vs. inhibition in theta-resonant pyramidal neurons

Rhythmic activity is ubiquitous in neural systems, with theta-resonant pyramidal neurons integrating rhythmic inputs in many cortical structures. Impedance analysis has been widely used to examine frequency-dependent responses of neuronal membranes to rhythmic inputs, but it assumes that the neuronal membrane is a linear system, requiring the use of small signals to stay in a near-linear regime. However, postsynaptic potentials are often large and trigger nonlinear mechanisms (voltage-gated ion channels). The goals of this work were to 1) develop an analysis method to evaluate membrane responses in this nonlinear domain and 2) explore phase relationships between rhythmic stimuli and subthreshold and spiking membrane potential (Vmemb) responses in models of theta-resonant pyramidal neurons. Responses in these output regimes were asymmetrical, with different phase shifts during hyperpolarizing and depolarizing half-cycles. Suprathreshold theta-rhythmic stimuli produced nonstationary Vmemb responses. Sinusoidal inputs produced "phase retreat": action potentials occurred progressively later in cycles of the input stimulus, resulting from adaptation. Sinusoidal current with increasing amplitude over cycles produced "phase advance": action potentials occurred progressively earlier. Phase retreat, phase advance, and subthreshold phase shifts were modulated by multiple ion channel conductances. Our results suggest differential responses of cortical neurons depending on the frequency of oscillatory input, which will play a role in neuronal responses to shifts in network state. We hypothesize that intrinsic cellular properties complement network properties and contribute to in vivo phase-shift phenomena such as phase precession, seen in place and grid cells, and phase roll, also observed in hippocampal CA1 neurons.NEW & NOTEWORTHY We augmented electrical impedance analysis to characterize phase shifts between large-amplitude current stimuli and nonlinear, asymmetric membrane potential responses. We predict different frequency-dependent phase shifts in response excitation vs. inhibition, as well as shifts in spike timing over multiple input cycles, in theta-resonant pyramidal neurons. We hypothesize that these effects contribute to navigation-related phenomena such as phase precession and phase roll. Our neuron-level hypothesis complements, rather than falsifies, prior network-level explanations of these phenomena.

A multi-layer mean-field model of the cerebellum embedding microstructure and population-specific dynamics

Mean-field (MF) models are computational formalism used to summarize in a few statistical parameters the salient biophysical properties of an inter-wired neuronal network. Their formalism normally incorporates different types of neurons and synapses along with their topological organization. MFs are crucial to efficiently implement the computational modules of large-scale models of brain function, maintaining the specificity of local cortical microcircuits. While MFs have been generated for the isocortex, they are still missing for other parts of the brain. Here we have designed and simulated a multi-layer MF of the cerebellar microcircuit (including Granule Cells, Golgi Cells, Molecular Layer Interneurons, and Purkinje Cells) and validated it against experimental data and the corresponding spiking neural network (SNN) microcircuit model. The cerebellar MF was built using a system of equations, where properties of neuronal populations and topological parameters are embedded in inter-dependent transfer functions. The model time constant was optimised using local field potentials recorded experimentally from acute mouse cerebellar slices as a template. The MF reproduced the average dynamics of different neuronal populations in response to various input patterns and predicted the modulation of the Purkinje Cells firing depending on cortical plasticity, which drives learning in associative tasks, and the level of feedforward inhibition. The cerebellar MF provides a computationally efficient tool for future investigations of the causal relationship between microscopic neuronal properties and ensemble brain activity in virtual brain models addressing both physiological and pathological conditions.

Long-Term Synaptic Plasticity Tunes the Gain of Information Channels through the Cerebellum Granular Layer

A central hypothesis on brain functioning is that long-term potentiation (LTP) and depression (LTD) regulate the signals transfer function by modifying the efficacy of synaptic transmission. In the cerebellum, granule cells have been shown to control the gain of signals transmitted through the mossy fiber pathway by exploiting synaptic inhibition in the glomeruli. However, the way LTP and LTD control signal transformation at the single-cell level in the space, time and frequency domains remains unclear. Here, the impact of LTP and LTD on incoming activity patterns was analyzed by combining patch-clamp recordings in acute cerebellar slices and mathematical modeling. LTP reduced the delay, increased the gain and broadened the frequency bandwidth of mossy fiber burst transmission, while LTD caused opposite changes. These properties, by exploiting NMDA subthreshold integration, emerged from microscopic changes in spike generation in individual granule cells such that LTP anticipated the emission of spikes and increased their number and precision, while LTD sorted the opposite effects. Thus, akin with the expansion recoding process theoretically attributed to the cerebellum granular layer, LTP and LTD could implement selective filtering lines channeling information toward the molecular and Purkinje cell layers for further processing.

A cerebellum inspired spiking neural network as a multi-model for pattern classification and robotic trajectory prediction

Spiking neural networks were introduced to understand spatiotemporal information processing in neurons and have found their application in pattern encoding, data discrimination, and classification. Bioinspired network architectures are considered for event-driven tasks, and scientists have looked at different theories based on the architecture and functioning. Motor tasks, for example, have networks inspired by cerebellar architecture where the granular layer recodes sparse representations of the mossy fiber (MF) inputs and has more roles in motor learning. Using abstractions from cerebellar connections and learning rules of deep learning network (DLN), patterns were discriminated within datasets, and the same algorithm was used for trajectory optimization. In the current work, a cerebellum-inspired spiking neural network with dynamics of cerebellar neurons and learning mechanisms attributed to the granular layer, Purkinje cell (PC) layer, and cerebellar nuclei interconnected by excitatory and inhibitory synapses was implemented. The model’s pattern discrimination capability was tested for two tasks on standard machine learning (ML) datasets and on following a trajectory of a low-cost sensor-free robotic articulator. Tuned for supervised learning, the pattern classification capability of the cerebellum-inspired network algorithm has produced more generalized models than data-specific precision models on smaller training datasets. The model showed an accuracy of 72%, which was comparable to standard ML algorithms, such as MLP (78%), Dl4jMlpClassifier (64%), RBFNetwork (71.4%), and libSVM-linear (85.7%). The cerebellar model increased the network’s capability and decreased storage, augmenting faster computations. Additionally, the network model could also implicitly reconstruct the trajectory of a 6-degree of freedom (DOF) robotic arm with a low error rate by reconstructing the kinematic parameters. The variability between the actual and predicted trajectory points was noted to be ± 3 cm (while moving to a position in a cuboid space of 25 × 30 × 40 cm). Although a few known learning rules were implemented among known types of plasticity in the cerebellum, the network model showed a generalized processing capability for a range of signals, modulating the data through the interconnected neural populations. In addition to potential use on sensor-free or feed-forward based controllers for robotic arms and as a generalized pattern classification algorithm, this model adds implications to motor learning theory.

Computational models of neurotransmission at cerebellar synapses unveil the impact on network computation

The neuroscientific field benefits from the conjoint evolution of experimental and computational techniques, allowing for the reconstruction and simulation of complex models of neurons and synapses. Chemical synapses are characterized by presynaptic vesicle cycling, neurotransmitter diffusion, and postsynaptic receptor activation, which eventually lead to postsynaptic currents and subsequent membrane potential changes. These mechanisms have been accurately modeled for different synapses and receptor types (AMPA, NMDA, and GABA) of the cerebellar cortical network, allowing simulation of their impact on computation. Of special relevance is short-term synaptic plasticity, which generates spatiotemporal filtering in local microcircuits and controls burst transmission and information flow through the network. Here, we present how data-driven computational models recapitulate the properties of neurotransmission at cerebellar synapses. The simulation of microcircuit models is starting to reveal how diverse synaptic mechanisms shape the spatiotemporal profiles of circuit activity and computation.

Purkinje Cell Activity Resonation Generates Rhythmic Behaviors at the Preferred Frequency of 8 Hz

Neural activity exhibits oscillations, bursts, and resonance, enhancing responsiveness at preferential frequencies. For example, theta-frequency bursting and resonance in granule cells facilitate synaptic transmission and plasticity mechanisms at the input stage of the cerebellar cortex. However, whether theta-frequency bursting of Purkinje cells is involved in generating rhythmic behavior has remained neglected. We recorded and optogenetically modulated the simple and complex spike activity of Purkinje cells while monitoring whisker movements with a high-speed camera of awake, head-fixed mice. During spontaneous whisking, both simple spike activity and whisker movement exhibit peaks within the theta band. Eliciting either simple or complex spikes at frequencies ranging from 0.5 to 28 Hz, we found that 8 Hz is the preferred frequency around which the largest movement is induced. Interestingly, oscillatory whisker movements at 8 Hz were also generated when simple spike bursting was induced at 2 and 4 Hz, but never via climbing fiber stimulation. These results indicate that 8 Hz is the resonant frequency at which the cerebellar-whisker circuitry produces rhythmic whisking.

The quest for multiscale brain modeling

Addressing the multiscale organization of the brain, which is fundamental to the dynamic repertoire of the organ, remains challenging. In principle, it should be possible to model neurons and synapses in detail and then connect them into large neuronal assemblies to explain the relationship between microscopic phenomena, large-scale brain functions, and behavior. It is more difficult to infer neuronal functions from ensemble measurements such as those currently obtained with brain activity recordings. In this article we consider theories and strategies for combining bottom-up models, generated from principles of neuronal biophysics, with top-down models based on ensemble representations of network activity and on functional principles. These integrative approaches are hoped to provide effective multiscale simulations in virtual brains and neurorobots, and pave the way to future applications in medicine and information technologies.

The voltage and spiking responses of subthreshold resonant neurons to structured and fluctuating inputs: persistence and loss of resonance and variability

We systematically investigate the response of neurons to oscillatory currents and synaptic-like inputs and we extend our investigation to non-structured synaptic-like spiking inputs with more realistic distributions of presynaptic spike times. We use two types of chirp-like inputs consisting of (i) a sequence of cycles with discretely increasing frequencies over time, and (ii) a sequence having the same cycles arranged in an arbitrary order. We develop and use a number of frequency-dependent voltage response metrics to capture the different aspects of the voltage response, including the standard impedance (Z) and the peak-to-trough amplitude envelope ([Formula: see text]) profiles. We show that Z-resonant cells (cells that exhibit subthreshold resonance in response to sinusoidal inputs) also show [Formula: see text]-resonance in response to sinusoidal inputs, but generally do not (or do it very mildly) in response to square-wave and synaptic-like inputs. In the latter cases the resonant response using Z is not predictive of the preferred frequencies at which the neurons spike when the input amplitude is increased above subthreshold levels. We also show that responses to conductance-based synaptic-like inputs are attenuated as compared to the response to current-based synaptic-like inputs, thus providing an explanation to previous experimental results. These response patterns were strongly dependent on the intrinsic properties of the participating neurons, in particular whether the unperturbed Z-resonant cells had a stable node or a focus. In addition, we show that variability emerges in response to chirp-like inputs with arbitrarily ordered patterns where all signals (trials) in a given protocol have the same frequency content and the only source of uncertainty is the subset of all possible permutations of cycles chosen for a given protocol. This variability is the result of the multiple different ways in which the autonomous transient dynamics is activated across cycles in each signal (different cycle orderings) and across trials. We extend our results to include high-rate Poisson distributed current- and conductance-based synaptic inputs and compare them with similar results using additive Gaussian white noise. We show that the responses to both Poisson-distributed synaptic inputs are attenuated with respect to the responses to Gaussian white noise. For cells that exhibit oscillatory responses to Gaussian white noise (band-pass filters), the response to conductance-based synaptic inputs are low-pass filters, while the response to current-based synaptic inputs may remain band-pass filters, consistent with experimental findings. Our results shed light on the mechanisms of communication of oscillatory activity among neurons in a network via subthreshold oscillations and resonance and the generation of network resonance.

Coherent theta oscillations in the cerebellum and supplementary motor area mediate visuomotor adaptation

The cerebellum and its interaction with cortical areas play a key role in our ability to flexibly adapt a motor program in response to sensory input. Current knowledge about specific neural mechanisms underlying the process of visuomotor adaptation is however lacking. Using a novel placement of EEG electrodes to record electric activity from the cerebellum, we studied local cerebellar activity, as well as its coupling with neocortical activity to obtain direct neurophysiological markers of visuomotor adaptation in humans. We found increased theta (4-8Hz) power in "cerebellar" as well as cortical electrodes, when subjects first encountered a visual manipulation. Theta power decreased as subjects adapted to the perturbation, and rebounded when the manipulation was suddenly removed. This effect was observed in two distinct locations: a cerebellar cluster and a central cluster, which were localized in left cerebellar crus I (lCB) and right supplementary motor area (rSMA) using linear constrained minimum variance beamforming. Importantly, we found that better adaptation was associated with increased theta power in left cerebellar electrodes and a right sensorimotor cortex electrode. Finally, increased rSMA -> lCB connectivity was significantly decreased with adaptation. These results demonstrate that: (1) cerebellar theta power is markedly modulated over the course of visuomotor adaptation and (2) theta oscillations could serve as a key mechanism for communication within a cortico-cerebellar loop.

Transcranial direct current stimulation of cerebellum alters spiking precision in cerebellar cortex: A modeling study of cellular responses

Transcranial direct current stimulation (tDCS) of the cerebellum has rapidly raised interest but the effects of tDCS on cerebellar neurons remain unclear. Assessing the cellular response to tDCS is challenging because of the uneven, highly stratified cytoarchitecture of the cerebellum, within which cellular morphologies, physiological properties, and function vary largely across several types of neurons. In this study, we combine MRI-based segmentation of the cerebellum and a finite element model of the tDCS-induced electric field (EF) inside the cerebellum to determine the field imposed on the cerebellar neurons throughout the region. We then pair the EF with multicompartment models of the Purkinje cell (PC), deep cerebellar neuron (DCN), and granule cell (GrC) and quantify the acute response of these neurons under various orientations, physiological conditions, and sequences of presynaptic stimuli. We show that cerebellar tDCS significantly modulates the postsynaptic spiking precision of the PC, which is expressed as a change in the spike count and timing in response to presynaptic stimuli. tDCS has modest effects, instead, on the PC tonic firing at rest and on the postsynaptic activity of DCN and GrC. In Purkinje cells, anodal tDCS shortens the repolarization phase following complex spikes (-14.7 ± 6.5% of baseline value, mean ± S.D.; max: -22.7%) and promotes burstiness with longer bursts compared to resting conditions. Cathodal tDCS, instead, promotes irregular spiking by enhancing somatic excitability and significantly prolongs the repolarization after complex spikes compared to baseline (+37.0 ± 28.9%, mean ± S.D.; max: +84.3%). tDCS-induced changes to the repolarization phase and firing pattern exceed 10% of the baseline values in Purkinje cells covering up to 20% of the cerebellar cortex, with the effects being distributed along the EF direction and concentrated in the area under the electrode over the cerebellum. Altogether, the acute effects of tDCS on cerebellum mainly focus on Purkinje cells and modulate the precision of the response to synaptic stimuli, thus having the largest impact when the cerebellar cortex is active. Since the spatiotemporal precision of the PC spiking is critical to learning and coordination, our results suggest cerebellar tDCS as a viable therapeutic option for disorders involving cerebellar hyperactivity such as ataxia.

Transcranial direct current stimulation (tDCS) of the cerebellum is gaining momentum as a neuromodulation tool for the treatment of neurological diseases like movement disorders. Nonetheless, the response of cells in the cerebellum to tDCS is unclear and hardly generalizes from our understanding of tDCS of the cerebral cortex. We use computational models to investigate the response of several types of cerebellar neurons to the electric field induced by tDCS and show that, differently from the cerebral cortex, tDCS has significant acute effects on the cerebellar cortex. These effects (i) primarily alter the way Purkinje cells encode synaptic stimuli from the molecular layer and (ii) can help hyperactive cells regain postsynaptic spiking precision. Since the spatiotemporal precision of the Purkinje cell spiking is critical to learning and coordination, the study shows how tDCS can operate at the cellular level to treat movement disorders like tremor and ataxia.

Altered Cerebellar Oscillations in Parkinson's Disease Patients during Cognitive and Motor Tasks

Structural and functional abnormalities in the cerebellar region have been shown in patients with Parkinson's disease (PD). Since the cerebellar region has been associated with cognitive and lower-limb motor functions, it is imperative to study cerebellar oscillations in PD. Here, we evaluated cerebellar electroencephalography (EEG) during cognitive processing and lower-limb motor performances in PD. Cortical and cerebellar EEG were collected from 74 PD patients and 37 healthy control subjects during a 7-second interval timing task, 26 PD patients and 13 controls during a lower-limb pedaling task, and 23 PD patients during eyes-open/closed resting conditions. Analyses were focused on the mid-cerebellar Cbz electrode and further compared to the mid-occipital Oz and mid-frontal Cz electrodes. Increased alpha-band power was observed during the eyes-closed resting-state condition over Oz, but no change in alpha power was observed over Cbz. PD patients showed higher dispersion when performing the 7-second interval timing cognitive task and executed the pedaling motor task with reduced speed compared to controls. PD patients exhibited attenuated cue-triggered theta-band power over Cbz during both the interval timing and pedaling motor tasks. Connectivity measures between Cbz and Cz showed theta-band differences, but only during the pedaling motor task. Cbz oscillatory activity also differed from Oz across multiple frequency bands in both groups during both tasks. Our cerebellar EEG data along with previous magnetoencephalography and animal model studies clearly show alterations in cerebellar oscillations during cognitive and motor processing in PD.

State-dependent entrainment of cerebellar nuclear neurons to the local field potential during voluntary movements

Understanding the relationship between the local field potential (LFP) and single neurons is essential if we are to understand network dynamics and the entrainment of neuronal activity. Here, we investigated the interaction between the LFP and single neurons recorded in the rat cerebellar nuclei (CN), which are part of the sensorimotor network, in freely moving rats. During movement, the LFP displayed persistent oscillations in the theta band frequency, whereas CN neurons displayed intermittent oscillations in the same frequency band contingent on the instantaneous LFP power; the neurons oscillated primarily when the concurrent LFP power was either high or low. Quantification of the relative instantaneous frequency and phase locking showed that CN neurons exhibited phase locked rhythmic activity at a frequency similar to that of the LFP or at a shifted frequency during high and low LFP power, respectively. We suggest that this nonlinear interaction between cerebellar neurons and the LFP power, which occurs solely during movement, contributes to the shaping of cerebellar output patterns.NEW & NOTEWORTHY We studied the interaction between single neurons and the LFP in the cerebellar nuclei of freely moving rats. We show that during movement, the neurons oscillated in the theta frequency band contingent on the concurrent LFP oscillation power in the same band; the neurons oscillated primarily when the LFP power was either high or low. We are the first to demonstrate a nonlinear, state-dependent entrainment of single neurons to the LFP.

On the Use of a Multimodal Optimizer for Fitting Neuron Models. Application to the Cerebellar Granule Cell

This article extends a recent methodological workflow for creating realistic and computationally efficient neuron models whilst capturing essential aspects of single-neuron dynamics. We overcome the intrinsic limitations of the extant optimization methods by proposing an alternative optimization component based on multimodal algorithms. This approach can natively explore a diverse population of neuron model configurations. In contrast to methods that focus on a single global optimum, the multimodal method allows directly obtaining a set of promising solutions for a single but complex multi-feature objective function. The final sparse population of candidate solutions has to be analyzed and evaluated according to the biological plausibility and their objective to the target features by the expert. In order to illustrate the value of this approach, we base our proposal on the optimization of cerebellar granule cell (GrC) models that replicate the essential properties of the biological cell. Our results show the emerging variability of plausible sets of values that this type of neuron can adopt underlying complex spiking characteristics. Also, the set of selected cerebellar GrC models captured spiking dynamics closer to the reference model than the single model obtained with off-the-shelf parameter optimization algorithms used in our previous article. The method hereby proposed represents a valuable strategy for adjusting a varied population of realistic and simplified neuron models. It can be applied to other kinds of neuron models and biological contexts.

Granular layEr Simulator: Design and Multi-GPU Simulation of the Cerebellar Granular Layer

In modern computational modeling, neuroscientists need to reproduce long-lasting activity of large-scale networks, where neurons are described by highly complex mathematical models. These aspects strongly increase the computational load of the simulations, which can be efficiently performed by exploiting parallel systems to reduce the processing times. Graphics Processing Unit (GPU) devices meet this need providing on desktop High Performance Computing. In this work, authors describe a novel Granular layEr Simulator development implemented on a multi-GPU system capable of reconstructing the cerebellar granular layer in a 3D space and reproducing its neuronal activity. The reconstruction is characterized by a high level of novelty and realism considering axonal/dendritic field geometries, oriented in the 3D space, and following convergence/divergence rates provided in literature. Neurons are modeled using Hodgkin and Huxley representations. The network is validated by reproducing typical behaviors which are well-documented in the literature, such as the center-surround organization. The reconstruction of a network, whose volume is 600 × 150 × 1,200 μm³ with 432,000 granules, 972 Golgi cells, 32,399 glomeruli, and 4,051 mossy fibers, takes 235 s on an Intel i9 processor. The 10 s activity reproduction takes only 4.34 and 3.37 h exploiting a single and multi-GPU desktop system (with one or two NVIDIA RTX 2080 GPU, respectively). Moreover, the code takes only 3.52 and 2.44 h if run on one or two NVIDIA V100 GPU, respectively. The relevant speedups reached (up to ~38× in the single-GPU version, and ~55× in the multi-GPU) clearly demonstrate that the GPU technology is highly suitable for realistic large network simulations.

Stellate cell computational modeling predicts signal filtering in the molecular layer circuit of cerebellum

The functional properties of cerebellar stellate cells and the way they regulate molecular layer activity are still unclear. We have measured stellate cells electroresponsiveness and their activation by parallel fiber bursts. Stellate cells showed intrinsic pacemaking, along with characteristic responses to depolarization and hyperpolarization, and showed a marked short-term facilitation during repetitive parallel fiber transmission. Spikes were emitted after a lag and only at high frequency, making stellate cells to operate as delay-high-pass filters. A detailed computational model summarizing these physiological properties allowed to explore different functional configurations of the parallel fiber-stellate cell-Purkinje cell circuit. Simulations showed that, following parallel fiber stimulation, Purkinje cells almost linearly increased their response with input frequency, but such an increase was inhibited by stellate cells, which leveled the Purkinje cell gain curve to its 4 Hz value. When reciprocal inhibitory connections between stellate cells were activated, the control of stellate cells over Purkinje cell discharge was maintained only at very high frequencies. These simulations thus predict a new role for stellate cells, which could endow the molecular layer with low-pass and band-pass filtering properties regulating Purkinje cell gain and, along with this, also burst delay and the burst-pause responses pattern.

Increased responsiveness at the cerebellar input stage in the PRRT2 knockout model of paroxysmal kinesigenic dyskinesia

PRoline-Rich Transmembrane protein-2 (PRRT2) is a recently described neuron-specific type-2 integral membrane protein with a large cytosolic N-terminal domain that distributes in presynaptic and axonal domains where it interacts with several presynaptic proteins and voltage-gated Na⁺ channels. Several PRRT2 mutations are the main cause of a wide and heterogeneous spectrum of paroxysmal disorders with a loss-of-function pathomechanism. The highest expression levels of PRRT2 in brain occurs in cerebellar granule cells (GCs) and cerebellar dysfunctions participate in the dyskinetic phenotype of PRRT2 knockout (KO) mice. We have investigated the effects of PRRT2 deficiency on the intrinsic excitability of GCs and the input-output relationships at the mossy fiber-GC synapses. We show that PRRT2 KO primary GCs display increased expression of Na⁺ channels, increased amplitude of Na⁺ currents and increased length of the axon initial segment, leading to an overall enhancement of intrinsic excitability. In acute PRRT2 KO cerebellar slices, GCs were more prone to action potential discharge in response to mossy fiber activation and exhibited an enhancement of transient and persistent Na⁺ currents, in the absence of changes at the mossy fiber-GC synapses. The results support a key role of PRRT2 expressed in GCs in the physiological regulation of the excitatory input to the cerebellum and are consistent with a major role of a cerebellar dysfunction in the pathogenesis of the PRRT2-linked paroxysmal pathologies.

Cerebellar Golgi cell models predict dendritic processing and mechanisms of synaptic plasticity

The Golgi cells are the main inhibitory interneurons of the cerebellar granular layer. Although recent works have highlighted the complexity of their dendritic organization and synaptic inputs, the mechanisms through which these neurons integrate complex input patterns remained unknown. Here we have used 8 detailed morphological reconstructions to develop multicompartmental models of Golgi cells, in which Na, Ca, and K channels were distributed along dendrites, soma, axonal initial segment and axon. The models faithfully reproduced a rich pattern of electrophysiological and pharmacological properties and predicted the operating mechanisms of these neurons. Basal dendrites turned out to be more tightly electrically coupled to the axon initial segment than apical dendrites. During synaptic transmission, parallel fibers caused slow Ca-dependent depolarizations in apical dendrites that boosted the axon initial segment encoder and Na-spike backpropagation into basal dendrites, while inhibitory synapses effectively shunted backpropagating currents. This oriented dendritic processing set up a coincidence detector controlling voltage-dependent NMDA receptor unblock in basal dendrites, which, by regulating local calcium influx, may provide the basis for spike-timing dependent plasticity anticipated by theory.

Cerebellar Cortex 4-12 Hz Oscillations and Unit Phase Relation in the Awake Rat

Oscillations in the granule cell layer (GCL) of the cerebellar cortex have been related to behavior and could facilitate communication with the cerebral cortex. These local field potential (LFP) oscillations, strong at 4–12 Hz in the rodent cerebellar cortex during awake immobility, should also be an indicator of an underlying influence on the patterns of the cerebellar cortex neuronal firing during rest. To address this hypothesis, cerebellar cortex LFPs and simultaneous single-neuron activity were collected during LFP oscillatory periods in the GCL of awake resting rats. During these oscillatory episodes, different types of units across the GCL and Purkinje cell layers showed variable phase-relation with the oscillatory cycles. Overall, 74% of the Golgi cell firing and 54% of the Purkinje cell simple spike (SS) firing were phase-locked with the oscillations, displaying a clear phase relationship. Despite this tendency, fewer Golgi cells (50%) and Purkinje cell’s SSs (25%) showed an oscillatory firing pattern. Oscillatory phase-locked spikes for the Golgi and Purkinje cells occurred towards the peak of the LFP cycle. GCL LFP oscillations had a strong capacity to predict the timing of Golgi cell spiking activity, indicating a strong influence of this oscillatory phenomenon over the GCL. Phase-locking was not as prominent for the Purkinje cell SS firing, indicating a weaker influence over the Purkinje cell layer, yet a similar phase relation. Overall, synaptic activity underlying GCL LFP oscillations likely exert an influence on neuronal population firing patterns in the cerebellar cortex in the awake resting state and could have a preparatory neural network shaping capacity serving as a neural baseline for upcoming cerebellar operations.

Functional Abnormalities of Cerebellum and Motor Cortex in Spinal Muscular Atrophy Mice

Spinal muscular atrophy (SMA) is a devastating genetic neuromuscular disease. Diffuse neuropathology has been reported in SMA patients and mouse models, however, functional changes in brain regions have not been studied. In the SMNΔ7 mouse model, we identified three types of differences in neuronal function in the cerebellum and motor cortex from two age groups: P7-9 (P7) and P11-14 (P11). Microelectrode array studies revealed significantly lower spontaneous firing and network activity in the cerebellum of SMA mice in both age groups, but it was more profound in the P11 group. In the motor cortex, however, neural activity was not different in either age group. Whole-cell patch-clamp was used to study the function of output neurons in both brain regions. In cerebellar Purkinje cells (PCs) of SMA mice, the input resistance was larger at P7, while capacitance was smaller at P11. In the motor cortex, no difference was observed in the passive membrane properties of layer V pyramidal neurons (PN5s). The action potential threshold of both types of output neurons was depolarized in the P11 group. We also observed lower spontaneous excitatory and inhibitory synaptic activity in PN5s and PCs respectively from P11 SMA mice. Overall, these differences suggest functional alterations in the neural network in these motor regions that change during development. Our results also suggest that neuronal dysfunction in these brain regions may contribute to the pathology of SMA. Comprehensive treatment strategies may consider motor regions outside of the spinal cord for better outcomes.

Cellular-resolution mapping uncovers spatial adaptive filtering at the rat cerebellum input stage

Long-term synaptic plasticity is thought to provide the substrate for adaptive computation in brain circuits but very little is known about its spatiotemporal organization. Here, we combined multi-spot two-photon laser microscopy in rat cerebellar slices with realistic modeling to map the distribution of plasticity in multi-neuronal units of the cerebellar granular layer. The units, composed by ~300 neurons activated by ~50 mossy fiber glomeruli, showed long-term potentiation concentrated in the core and long-term depression in the periphery. This plasticity was effectively accounted for by an NMDA receptor and calcium-dependent induction rule and was regulated by the inhibitory Golgi cell loops. Long-term synaptic plasticity created effective spatial filters tuning the time-delay and gain of spike retransmission at the cerebellum input stage and provided a plausible basis for the spatiotemporal recoding of input spike patterns anticipated by the motor learning theory.

Casali, Tognolina et al. use two-photon laser microscopy to spatially map long-term synaptic plasticity in rat cerebellar granular cells following stimulation of mossy fibers. Their data allow them to apply realistic modeling to test hypotheses about the synaptic spiking dynamics and reveal the importance of synaptic inhibition to defining these microcircuits.

Optimization of Efficient Neuron Models With Realistic Firing Dynamics. The Case of the Cerebellar Granule Cell

Biologically relevant large-scale computational models currently represent one of the main methods in neuroscience for studying information processing primitives of brain areas. However, biologically realistic neuron models tend to be computationally heavy and thus prevent these models from being part of brain-area models including thousands or even millions of neurons. The cerebellar input layer represents a canonical example of large scale networks. In particular, the cerebellar granule cells, the most numerous cells in the whole mammalian brain, have been proposed as playing a pivotal role in the creation of somato-sensorial information representations. Enhanced burst frequency (spiking resonance) in the granule cells has been proposed as facilitating the input signal transmission at the theta-frequency band (4–12 Hz), but the functional role of this cell feature in the operation of the granular layer remains largely unclear. This study aims to develop a methodological pipeline for creating neuron models that maintain biological realism and computational efficiency whilst capturing essential aspects of single-neuron processing. Therefore, we selected a light computational neuron model template (the adaptive-exponential integrate-and-fire model), whose parameters were progressively refined using an automatic parameter tuning with evolutionary algorithms (EAs). The resulting point-neuron models are suitable for reproducing the main firing properties of a realistic granule cell from electrophysiological measurements, including the spiking resonance at the theta-frequency band, repetitive firing according to a specified intensity-frequency (I-F) curve and delayed firing under current-pulse stimulation. Interestingly, the proposed model also reproduced some other emergent properties (namely, silent at rest, rheobase and negligible adaptation under depolarizing currents) even though these properties were not set in the EA as a target in the fitness function (FF), proving that these features are compatible even in computationally simple models. The proposed methodology represents a valuable tool for adjusting AdEx models according to a FF defined in the spiking regime and based on biological data. These models are appropriate for future research of the functional implication of bursting resonance at the theta band in large-scale granular layer network models.

Parameter tuning differentiates granule cell subtypes enriching transmission properties at the cerebellum input stage

The cerebellar granule cells (GrCs) are classically described as a homogeneous neuronal population discharging regularly without adaptation. We show that GrCs in fact generate diverse response patterns to current injection and synaptic activation, ranging from adaptation to acceleration of firing. Adaptation was predicted by parameter optimization in detailed computational models based on available knowledge on GrC ionic channels. The models also predicted that acceleration required additional mechanisms. We found that yet unrecognized TRPM4 currents specifically accounted for firing acceleration and that adapting GrCs outperformed accelerating GrCs in transmitting high-frequency mossy fiber (MF) bursts over a background discharge. This implied that GrC subtypes identified by their electroresponsiveness corresponded to specific neurotransmitter release probability values. Simulations showed that fine-tuning of pre- and post-synaptic parameters generated effective MF-GrC transmission channels, which could enrich the processing of input spike patterns and enhance spatio-temporal recoding at the cerebellar input stage.

Intrinsic and Synaptic Properties Shaping Diverse Behaviors of Neural Dynamics

The majority of neurons in the neuronal system of the brain have a complex morphological structure, which diversifies the dynamics of neurons. In the granular layer of the cerebellum, there exists a unique cell type, the unipolar brush cell (UBC), that serves as an important relay cell for transferring information from outside mossy fibers to downstream granule cells. The distinguishing feature of the UBC is that it has a simple morphology, with only one short dendritic brush connected to its soma. Based on experimental evidence showing that UBCs exhibit a variety of dynamic behaviors, here we develop two simple models, one with a few detailed ion channels for simulation and the other one as a two-variable dynamical system for theoretical analysis, to characterize the intrinsic dynamics of UBCs. The reasonable values of the key channel parameters of the models can be determined by analysis of the stability of the resting membrane potential and the rebound firing properties of UBCs. Considered together with a large variety of synaptic dynamics installed on UBCs, we show that the simple-structured UBCs, as relay cells, can extend the range of dynamics and information from input mossy fibers to granule cells with low-frequency resonance and transfer stereotyped inputs to diverse amplitudes and phases of the output for downstream granule cells. These results suggest that neuronal computation, embedded within intrinsic ion channels and the diverse synaptic properties of single neurons without sophisticated morphology, can shape a large variety of dynamic behaviors to enhance the computational ability of local neuronal circuits.

Generating Brain Waves, the Power of Astrocytes

Synchronization of neuronal activity in the brain underlies the emergence of neuronal oscillations termed “brain waves”, which serve various physiological functions and correlate with different behavioral states. It has been postulated that at least ten distinct mechanisms are involved in the formulation of these brain waves, including variations in the concentration of extracellular neurotransmitters and ions, as well as changes in cellular excitability. In this mini review we highlight the contribution of astrocytes, a subtype of glia, in the formation and modulation of brain waves mainly due to their close association with synapses that allows their bidirectional interaction with neurons, and their syncytium-like activity via gap junctions that facilitate communication to distal brain regions through Ca²⁺ waves. These capabilities allow astrocytes to regulate neuronal excitability via glutamate uptake, gliotransmission and tight control of the extracellular K⁺ levels via a process termed K⁺ clearance. Spatio-temporal synchrony of activity across neuronal and astrocytic networks, both locally and distributed across cortical regions, underpins brain states and thereby behavioral states, and it is becoming apparent that astrocytes play an important role in the development and maintenance of neural activity underlying these complex behavioral states.

Frequency-dependent responses of neuronal models to oscillatory inputs in current versus voltage clamp

Action potential generation in neurons depends on a membrane potential threshold and therefore on how subthreshold inputs influence this voltage. In oscillatory networks, for example, many neuron types have been shown to produce membrane potential ([Formula: see text]) resonance: a maximum subthreshold response to oscillatory inputs at a nonzero frequency. Resonance is usually measured by recording [Formula: see text] in response to a sinusoidal current ([Formula: see text]), applied at different frequencies (f), an experimental setting known as current clamp (I-clamp). Several recent studies, however, use the voltage clamp (V-clamp) method to control [Formula: see text] with a sinusoidal input at different frequencies [[Formula: see text]] and measure the total membrane current ([Formula: see text]). The two methods obey systems of differential equations of different dimensionality, and while I-clamp provides a measure of electrical impedance [[Formula: see text]], V-clamp measures admittance [[Formula: see text]]. We analyze the relationship between these two measurement techniques. We show that, despite different dimensionality, in linear systems the two measures are equivalent: [Formula: see text]. However, nonlinear model neurons produce different values for Z and [Formula: see text]. In particular, nonlinearities in the voltage equation produce a much larger difference between these two quantities than those in equations of recovery variables that describe activation and inactivation kinetics. Neurons are inherently nonlinear, and notably, with ionic currents that amplify resonance, the voltage clamp technique severely underestimates the current clamp response. We demonstrate this difference experimentally using the PD neurons in the crab stomatogastric ganglion. These findings are instructive for researchers who explore cellular mechanisms of neuronal oscillations.

Anatomical and physiological foundations of cerebello-hippocampal interaction

Multiple lines of evidence suggest that functionally intact cerebello-hippocampal interactions are required for appropriate spatial processing. However, how the cerebellum anatomically and physiologically engages with the hippocampus to sustain such communication remains unknown. Using rabies virus as a retrograde transneuronal tracer in mice, we reveal that the dorsal hippocampus receives input from topographically restricted and disparate regions of the cerebellum. By simultaneously recording local field potential from both the dorsal hippocampus and anatomically connected cerebellar regions, we additionally suggest that the two structures interact, in a behaviorally dynamic manner, through subregion-specific synchronization of neuronal oscillations in the 6-12 Hz frequency range. Together, these results reveal a novel neural network macro-architecture through which we can understand how a brain region classically associated with motor control, the cerebellum, may influence hippocampal neuronal activity and related functions, such as spatial navigation.

Complex Electroresponsive Dynamics in Olivocerebellar Neurons Represented With Extended-Generalized Leaky Integrate and Fire Models

The neurons of the olivocerebellar circuit exhibit complex electroresponsive dynamics, which are thought to play a fundamental role for network entraining, plasticity induction, signal processing, and noise filtering. In order to reproduce these properties in single-point neuron models, we have optimized the Extended-Generalized Leaky Integrate and Fire (E-GLIF) neuron through a multi-objective gradient-based algorithm targeting the desired input–output relationships. In this way, E-GLIF was tuned toward the unique input–output properties of Golgi cells, granule cells, Purkinje cells, molecular layer interneurons, deep cerebellar nuclei cells, and inferior olivary cells. E-GLIF proved able to simulate the complex cell-specific electroresponsive dynamics of the main olivocerebellar neurons including pacemaking, adaptation, bursting, post-inhibitory rebound excitation, subthreshold oscillations, resonance, and phase reset. The integration of these E-GLIF point-neuron models into olivocerebellar Spiking Neural Networks will allow to evaluate the impact of complex electroresponsive dynamics at the higher scales, up to motor behavior, in closed-loop simulations of sensorimotor tasks.

Reconstruction and Simulation of a Scaffold Model of the Cerebellar Network

Reconstructing neuronal microcircuits through computational models is fundamental to simulate local neuronal dynamics. Here a scaffold model of the cerebellum has been developed in order to flexibly place neurons in space, connect them synaptically, and endow neurons and synapses with biologically-grounded mechanisms. The scaffold model can keep neuronal morphology separated from network connectivity, which can in turn be obtained from convergence/divergence ratios and axonal/dendritic field 3D geometries. We first tested the scaffold on the cerebellar microcircuit, which presents a challenging 3D organization, at the same time providing appropriate datasets to validate emerging network behaviors. The scaffold was designed to integrate the cerebellar cortex with deep cerebellar nuclei (DCN), including different neuronal types: Golgi cells, granule cells, Purkinje cells, stellate cells, basket cells, and DCN principal cells. Mossy fiber inputs were conveyed through the glomeruli. An anisotropic volume (0.077 mm³) of mouse cerebellum was reconstructed, in which point-neuron models were tuned toward the specific discharge properties of neurons and were connected by exponentially decaying excitatory and inhibitory synapses. Simulations using both pyNEST and pyNEURON showed the emergence of organized spatio-temporal patterns of neuronal activity similar to those revealed experimentally in response to background noise and burst stimulation of mossy fiber bundles. Different configurations of granular and molecular layer connectivity consistently modified neuronal activation patterns, revealing the importance of structural constraints for cerebellar network functioning. The scaffold provided thus an effective workflow accounting for the complex architecture of the cerebellar network. In principle, the scaffold can incorporate cellular mechanisms at multiple levels of detail and be tuned to test different structural and functional hypotheses. A future implementation using detailed 3D multi-compartment neuron models and dynamic synapses will be needed to investigate the impact of single neuron properties on network computation.

Consensus paper: Decoding the Contributions of the Cerebellum as a Time Machine. From Neurons to Clinical Applications

Time perception is an essential element of conscious and subconscious experience, coordinating our perception and interaction with the surrounding environment. In recent years, major technological advances in the field of neuroscience have helped foster new insights into the processing of temporal information, including extending our knowledge of the role of the cerebellum as one of the key nodes in the brain for this function. This consensus paper provides a state-of-the-art picture from the experts in the field of the cerebellar research on a variety of crucial issues related to temporal processing, drawing on recent anatomical, neurophysiological, behavioral, and clinical research.The cerebellar granular layer appears especially well-suited for timing operations required to confer millisecond precision for cerebellar computations. This may be most evident in the manner the cerebellum controls the duration of the timing of agonist-antagonist EMG bursts associated with fast goal-directed voluntary movements. In concert with adaptive processes, interactions within the cerebellar cortex are sufficient to support sub-second timing. However, supra-second timing seems to require cortical and basal ganglia networks, perhaps operating in concert with cerebellum. Additionally, sensory information such as an unexpected stimulus can be forwarded to the cerebellum via the climbing fiber system, providing a temporally constrained mechanism to adjust ongoing behavior and modify future processing. Patients with cerebellar disorders exhibit impairments on a range of tasks that require precise timing, and recent evidence suggest that timing problems observed in other neurological conditions such as Parkinson's disease, essential tremor, and dystonia may reflect disrupted interactions between the basal ganglia and cerebellum.The complex concepts emerging from this consensus paper should provide a foundation for further discussion, helping identify basic research questions required to understand how the brain represents and utilizes time, as well as delineating ways in which this knowledge can help improve the lives of those with neurological conditions that disrupt this most elemental sense. The panel of experts agrees that timing control in the brain is a complex concept in whom cerebellar circuitry is deeply involved. The concept of a timing machine has now expanded to clinical disorders.

Membrane potential resonance in non-oscillatory neurons interacts with synaptic connectivity to produce network oscillations

Several neuron types have been shown to exhibit (subthreshold) membrane potential resonance (MPR), defined as the occurrence of a peak in their voltage amplitude response to oscillatory input currents at a preferred (resonant) frequency. MPR has been investigated both experimentally and theoretically. However, whether MPR is simply an epiphenomenon or it plays a functional role for the generation of neuronal network oscillations and how the latent time scales present in individual, non-oscillatory cells affect the properties of the oscillatory networks in which they are embedded are open questions. We address these issues by investigating a minimal network model consisting of (i) a non-oscillatory linear resonator (band-pass filter) with 2D dynamics, (ii) a passive cell (low-pass filter) with 1D linear dynamics, and (iii) nonlinear graded synaptic connections (excitatory or inhibitory) with instantaneous dynamics. We demonstrate that (i) the network oscillations crucially depend on the presence of MPR in the resonator, (ii) they are amplified by the network connectivity, (iii) they develop relaxation oscillations for high enough levels of mutual inhibition/excitation, and (iv) the network frequency monotonically depends on the resonators resonant frequency. We explain these phenomena using a reduced adapted version of the classical phase-plane analysis that helps uncovering the type of effective network nonlinearities that contribute to the generation of network oscillations. We extend our results to networks having cells with 2D dynamics. Our results have direct implications for network models of firing rate type and other biological oscillatory networks (e.g, biochemical, genetic).

Spike burst-pause dynamics of Purkinje cells regulate sensorimotor adaptation

Cerebellar Purkinje cells mediate accurate eye movement coordination. However, it remains unclear how oculomotor adaptation depends on the interplay between the characteristic Purkinje cell response patterns, namely tonic, bursting, and spike pauses. Here, a spiking cerebellar model assesses the role of Purkinje cell firing patterns in vestibular ocular reflex (VOR) adaptation. The model captures the cerebellar microcircuit properties and it incorporates spike-based synaptic plasticity at multiple cerebellar sites. A detailed Purkinje cell model reproduces the three spike-firing patterns that are shown to regulate the cerebellar output. Our results suggest that pauses following Purkinje complex spikes (bursts) encode transient disinhibition of target medial vestibular nuclei, critically gating the vestibular signals conveyed by mossy fibres. This gating mechanism accounts for early and coarse VOR acquisition, prior to the late reflex consolidation. In addition, properly timed and sized Purkinje cell bursts allow the ratio between long-term depression and potentiation (LTD/LTP) to be finely shaped at mossy fibre-medial vestibular nuclei synapses, which optimises VOR consolidation. Tonic Purkinje cell firing maintains the consolidated VOR through time. Importantly, pauses are crucial to facilitate VOR phase-reversal learning, by reshaping previously learnt synaptic weight distributions. Altogether, these results predict that Purkinje spike burst-pause dynamics are instrumental to VOR learning and reversal adaptation.

Cerebellar Purkinje cells regulate accurate eye movement coordination. However, it remains unclear how cerebellar-dependent oculomotor adaptation depends on the interplay between Purkinje cell characteristic response patterns: tonic, high frequency bursting, and post-complex spike pauses. We explore the role of Purkinje spike burst-pause dynamics in VOR adaptation. A biophysical model of Purkinje cell is at the core of a spiking network model, which captures the cerebellar microcircuit properties and incorporates spike-based synaptic plasticity mechanisms at different cerebellar sites. We show that Purkinje spike burst-pause dynamics are critical for (1) gating the vestibular-motor response association during VOR acquisition; (2) mediating the LTD/LTP balance for VOR consolidation; (3) reshaping synaptic efficacy distributions for VOR phase-reversal adaptation; (4) explaining the reversal VOR gain discontinuities during sleeping.

Proprioceptive Recognition with Artificial Neural Networks Based on Organizations of Spinocerebellar Tract and Cerebellum

Muscle kinematics and kinetics are nonlinearly encoded by proprioceptors, and the changes in muscle length and velocity are integrated into Ia afferent. Besides, proprioceptive signals from multiple muscles are probably mixed in afferent pathways, which all lead to difficulties in proprioceptive recognition for the cerebellum. In this study, the artificial neural networks, whose organizations are biologically based on the spinocerebellar tract and cerebellum, are utilized to decode the proprioceptive signals. Consistent with the controversy of the proprioceptive division in the dorsal spinocerebellar tract, the spinocerebellar tract networks incorporated two distinct inferences, (1) the centralized networks, which mixed Ia, II, and Ib and processed them together; (2) the decentralized networks, which processed Ia, II, and Ib afferents separately. The cerebellar networks were based on the Marr-Albus model to recognize the kinematic states. The networks were trained by a specific movement, and the trained networks were subsequently required to predict kinematic states of six movements. The results demonstrated that the centralized networks, which were more consistent with the physiological findings in recent years, had better recognition accuracy than the decentralized networks, and the networks were still effective even when proprioceptive afferents from multiple muscles were integrated.

Potentiation of cerebellar Purkinje cells facilitates whisker reflex adaptation through increased simple spike activity

Cerebellar plasticity underlies motor learning. However, how the cerebellum operates to enable learned changes in motor output is largely unknown. We developed a sensory-driven adaptation protocol for reflexive whisker protraction and recorded Purkinje cell activity from crus 1 and 2 of awake mice. Before training, simple spikes of individual Purkinje cells correlated during reflexive protraction with the whisker position without lead or lag. After training, simple spikes and whisker protractions were both enhanced with the spiking activity now leading behavioral responses. Neuronal and behavioral changes did not occur in two cell-specific mouse models with impaired long-term potentiation at their parallel fiber to Purkinje cell synapses. Consistent with cerebellar plasticity rules, increased simple spike activity was prominent in cells with low complex spike response probability. Thus, potentiation at parallel fiber to Purkinje cell synapses may contribute to reflex adaptation and enable expression of cerebellar learning through increases in simple spike activity.

Rodents use their whiskers to explore the world around them. When the whiskers touch an object, it triggers involuntary movements of the whiskers called whisker reflexes. Experiencing the same sensory stimulus multiple times enables rodents to fine-tune these reflexes, e.g., by making their movements larger or smaller. This type of learning is often referred to as motor learning.

A part of the brain called cerebellum controls motor learning. It contains some of the largest neurons in the nervous system, the Purkinje cells. Each Purkinje cell receives input from thousands of extensions of small neurons, known as parallel fibers. It is thought that decreasing the strength of the connections between parallel fibers and Purkinje cells can help mammals learn new movements. This is the case in a type of learning called Pavlovian conditioning. It takes its name from the Russian scientist, Pavlov, who showed that dogs can learn to salivate in response to a bell signaling food. Pavlovian conditioning enables animals to optimize their responses to sensory stimuli.

But Romano et al. now show that increasing the strength of connections between parallel fibers and Purkinje cells can also support learning. To trigger reflexive whisker movements, a machine blew puffs of air onto the whiskers of awake mice. After repeated exposure to the air puffs, the mice increased the size of their whisker reflexes. At the same time, their Purkinje cells became more active and the connections between Purkinje cells and parallel fibers grew stronger. Artificially increasing Purkinje cell activity triggered the same changes in whisker reflexes as the air puffs themselves.

Textbooks still report that only weakening of connections within the cerebellum enables animals to learn and modify movements. The data obtained by Romano al. thus paint a new picture of how the cerebellum works in the context of whisker learning. They show that strengthening these connections can also support movement-related learning.

Spontaneous recovery of conditioned eyeblink responses is associated with transiently decreased cerebellar theta activity in guinea pigs

Behavioral studies have demonstrated that extinguished conditioned eyeblink responses (CR) can spontaneously recover after extinction. However, the neural mechanisms underlying this process are still unclear. We have shown that spontaneous cerebellar theta activity was predictive of subsequent CR extinction. Here, we sought to further evaluate the association between spontaneous recovery and cerebellar theta activity in behaving guinea pigs. It was found that trace conditioning training significantly diminished the degree of spontaneous recovery during extinction sessions as compared to delay training. Moreover, by recording local field potential in the cerebellum of guinea pigs undergoing an eyeblink conditioning extinction task, we found that spontaneous recovery of delay-paradigm CRs was associated with transiently decreased CS-evoked theta activity in the cerebellum. These findings suggest that decreased CS-evoked cerebellar theta activity may contribute to the neural process that is important for the spontaneous recovery of extinguished motor memory. Future studies are needed to clarify the neural mechanism underlying changed cerebellar theta activity during altered behavioral contingencies.

Neural Coding With Bursts-Current State and Future Perspectives

Neuronal action potentials or spikes provide a long-range, noise-resistant means of communication between neurons. As point processes single spikes contain little information in themselves, i.e., outside the context of spikes from other neurons. Moreover, they may fail to cross a synapse. A burst, which consists of a short, high frequency train of spikes, will more reliably cross a synapse, increasing the likelihood of eliciting a postsynaptic spike, depending on the specific short-term plasticity at that synapse. Both the number and the temporal pattern of spikes in a burst provide a coding space that lies within the temporal integration realm of single neurons. Bursts have been observed in many species, including the non-mammalian, and in brain regions that range from subcortical to cortical. Despite their widespread presence and potential relevance, the uncertainties of how to classify bursts seems to have limited the research into the coding possibilities for bursts. The present series of research articles provides new insights into the relevance and interpretation of bursts across different neural circuits, and new methods for their analysis. Here, we provide a succinct introduction to the history of burst coding and an overview of recent work on this topic.

Gap Junction Modulation of Low-Frequency Oscillations in the Cerebellar Granule Cell Layer

Local field potential (LFP) oscillations in the granule cell layer (GCL) of the cerebellar cortex have been identified previously in the awake rat and monkey during immobility. These low-frequency oscillations are thought to be generated through local circuit interactions between Golgi cells and granule cells within the GCL. Golgi cells display rhythmic firing and pacemaking properties, and also are electrically coupled through gap junctions within the GCL. Here, we tested if gap junctions in the rat cerebellar cortex contribute to the generation of LFP oscillations in the GCL. We recorded LFP oscillations under urethane anesthesia, and examined the effects of local infusion of gap junction blockers on 5-15 Hz oscillations. Local infusion of the gap junction blockers carbenoxolone and mefloquine resulted in significant decreases in the power of oscillations over a 30-min period, but the power of oscillations was unchanged in control experiments following vehicle injections. In addition, infusion of gap junction blockers had no significant effect on multi-unit activity, suggesting that the attenuation of low-frequency oscillations was likely due to reductions in electrical coupling rather than a decreased excitability within the granule cell layer. Our results indicate that electrical coupling among the Golgi cell networks in the cerebellar cortex contributes to the local circuit mechanisms that promote the occurrence of GCL LFP slow oscillations in the anesthetized rat.

Physiology of the cerebellum

The cerebellum is a central brain structure deeply integrated into major loops with the cerebral cortex, brainstem, and spinal cord. The cerebellum shows a complex regional organization consisting of modules with sagittal orientation. The cerebellum takes part in motor control and its lesions cause a movement incoordination syndrome called ataxia. Recent observations also imply involvement of the cerebellum in cognition and executive control, with an impact on pathologies like dyslexia and autism. The cerebellum operates as a forward controller learning to predict the precise timing of correlated events. The physiologic mechanisms of cerebellar functioning are still the object of intense research. The signals entering the cerebellum through the mossy fibers are processed in the granular layer and transmitted to Purkinje cells, while a collateral pathway activates the deep cerebellar nuclei (DCN). Purkinje cells in turn inhibit DCN, so that the cerebellar cortex operates as a side loop controlling the DCN. Learning is now known to occur through synaptic plasticity at multiple synapses in the granular layer, molecular layer, and DCN, extending the original concept of the Motor Learning Theory that predicted a single form of plasticity at the synapse between parallel fibers and Purkinje cells under the supervision of climbing fibers deriving from the inferior olive. Coordination derives from the precise regulation of timing and gain in the different cerebellar modules. The investigation of cerebellar dynamics using advanced physiologic recordings and computational models is now providing new clues on how the cerebellar network performs its internal computations.

Spiking resonances in models with the same slow resonant and fast amplifying currents but different subthreshold dynamic properties

The generation of spiking resonances in neurons (preferred spiking responses to oscillatory inputs) requires the interplay of the intrinsic ionic currents that operate at the subthreshold voltage level and the spiking mechanisms. Combinations of the same types of ionic currents in different parameter regimes may give rise to different types of nonlinearities in the voltage equation (e.g., parabolic- and cubic-like), generating subthreshold (membrane potential) oscillations patterns with different properties. These nonlinearities are not apparent in the model equations, but can be uncovered by plotting the voltage nullclines in the phase-plane diagram. We investigate the spiking resonant properties of conductance-based models that are biophysically equivalent at the subthreshold level (same ionic currents), but dynamically different (parabolic- and cubic-like voltage nullclines). As a case study we consider a model having a persistent sodium and a hyperpolarization-activated (h-) currents, which exhibits subthreshold resonance in the theta frequency band. We unfold the concept of spiking resonance into evoked and output spiking resonance. The former focuses on the input frequencies that are able to generate spikes, while the latter focuses on the output spiking frequencies regardless of the input frequency that generated these spikes. A cell can exhibit one or both types of resonances. We also measure spiking phasonance, which is an extension of subthreshold phasonance (zero-phase-shift response to oscillatory inputs) to the spiking regime. The subthreshold resonant properties of both types of models are communicated to the spiking regime for low enough input amplitudes as the voltage response for the subthreshold resonant frequency band raises above threshold. For higher input amplitudes evoked spiking resonance is no longer present in these models, but output spiking resonance is present primarily in the parabolic-like model due to a cycle skipping mechanism (involving mixed-mode oscillations), while the cubic-like model shows a better 1:1 entrainment. We use dynamical systems tools to explain the underlying mechanisms and the mechanistic differences between the resonance types. Our results demonstrate that the effective time scales that operate at the subthreshold regime to generate intrinsic subthreshold oscillations, mixed-mode oscillations and subthreshold resonance do not necessarily determine the existence of a preferred spiking response to oscillatory inputs in the same frequency band. The results discussed in this paper highlight both the complexity of the suprathreshold responses to oscillatory inputs in neurons having resonant and amplifying currents with different time scales and the fact that the identity of the participating ionic currents is not enough to predict the resulting patterns, but additional dynamic information, captured by the geometric properties of the phase-space diagram, is needed.

Single Neuron Optimization as a Basis for Accurate Biophysical Modeling: The Case of Cerebellar Granule Cells

In realistic neuronal modeling, once the ionic channel complement has been defined, the maximum ionic conductance (G_i-max) values need to be tuned in order to match the firing pattern revealed by electrophysiological recordings. Recently, selection/mutation genetic algorithms have been proposed to efficiently and automatically tune these parameters. Nonetheless, since similar firing patterns can be achieved through different combinations of G_i-max values, it is not clear how well these algorithms approximate the corresponding properties of real cells. Here we have evaluated the issue by exploiting a unique opportunity offered by the cerebellar granule cell (GrC), which is electrotonically compact and has therefore allowed the direct experimental measurement of ionic currents. Previous models were constructed using empirical tuning of G_i-max values to match the original data set. Here, by using repetitive discharge patterns as a template, the optimization procedure yielded models that closely approximated the experimental G_i-max values. These models, in addition to repetitive firing, captured additional features, including inward rectification, near-threshold oscillations, and resonance, which were not used as features. Thus, parameter optimization using genetic algorithms provided an efficient modeling strategy for reconstructing the biophysical properties of neurons and for the subsequent reconstruction of large-scale neuronal network models.

Competition between Persistent Na+ and Muscarine-Sensitive K+ Currents Shapes Perithreshold Resonance and Spike Tuning in CA1 Pyramidal Neurons

Neurons from many brain regions display intrinsic subthreshold theta-resonance, responding preferentially to theta-frequency oscillatory stimuli. Resonance may contribute to selective communication among neurons and to orchestrate brain rhythms. CA1 pyramidal neurons receive theta activity, generating place fields. In these neurons the expression of perithreshold frequency preference is controversial, particularly in the spiking regime, with evidence favoring either non-resonant (integrator-like) or resonant behavior. Perithreshold dynamics depends on the persistent Na⁺ current I_NaP developing above −70 mV and the muscarine-sensitive K⁺ current I_M activating above −60 mV. We conducted current and voltage clamp experiments in slices to investigate perithreshold excitability of CA1 neurons under oscillatory stimulation. Around 20% of neurons displayed perithreshold resonance that is expressed in spiking. The remaining neurons (~80%) acted as low-pass filters lacking frequency preference. Paired voltage clamp measurement of I_NaP and I_M showed that perithreshold activation of I_M is in general low while I_NaP is high enough to depolarize neurons toward threshold before resonance expression, explaining the most abundant non-resonant perithreshold behavior. Partial blockade of I_NaP by pharmacological tools or dynamic clamp changed non-resonant to resonant behavior. Furthermore, shifting I_M activation toward hyperpolarized potentials by dynamic clamp also transformed non-resonant neurons into resonant ones. We propose that the relative levels of I_NaP and I_M control perithreshold behavior of CA1 neurons constituting a gating mechanism for theta resonance in the spiking regime. Both currents are regulated by intracellular signaling and neuromodulators which may allow dynamic switching of perithreshold behavior between resonant and non-resonant.

BIOPHYSICAL PROPERTIES OF SUBTHRESHOLD RESONANCE OSCILLATIONS AND SUBTHRESHOLD MEMBRANE OSCILLATIONS IN NEURONS

Subthreshold-level activities in neurons play a crucial role in neuronal oscillations. These small-amplitude oscillations have been suggested to be involved in synaptic plasticity and in determining the frequency of network oscillations. Subthreshold membrane oscillations (STOs) and subthreshold resonance oscillations (SROs) are the main constituents of subthreshold-level activities in neurons. In this study, a general theoretical framework for analyzing the mechanisms underlying STOs and SROs in neurons is presented. Results showed that the resting membrane potential and the hyperpolarization-activated potassium channel (h-channel) affect the subthreshold-level activities in stellate cells. The contribution of h-channel on resonance is attributed to its large time constant, which produces the time lag between I_h and the membrane potential. Conversely, the persistent sodium channels (Nap-channels) only play an amplifying role in these neurons.

Ionic mechanisms maintaining action potential conduction velocity at high firing frequencies in an unmyelinated axon

The descending contralateral movement detector (DCMD) is a high‐performance interneuron in locusts with an axon capable of transmitting action potentials (AP) at more than 500 Hz. We investigated biophysical mechanisms for fidelity of high‐frequency transmission in this axon. We measured conduction velocities (CVs) at room temperature during exposure to 10 mmol/L cadmium, a calcium current antagonist, and found significant reduction in CV with reduction at frequencies >200 Hz of ~10%. Higher temperatures induced greater CV reductions during exposure to cadmium across all frequencies of ~20–30%. Intracellular recordings during 15 min of exposure to cadmium or nickel, also a calcium current antagonist, revealed an increase in the magnitude of the afterhyperpolarization potential (AHP) and the time to recover to baseline after the AHP (Medians for Control: −19.8%; Nickel: 167.2%; Cadmium: 387.2%), that could be due to a T‐type calcium current. However, the removal of extracellular calcium did not mimic divalent cation exposure suggesting calcium currents are not the cause of the AHP increase. Computational modeling showed that the effects of the divalent cations could be modeled with a persistent sodium current which could be blocked by high concentrations of divalent cations. Persistent sodium current shortened the AHP duration in our models and increased CV for high‐frequency APs. We suggest that faithful, high‐frequency axonal conduction in the DCMD is enabled by a mechanism that shortens the AHP duration like a persistent or resurgent sodium current.

Modeling the Cerebellar Microcircuit: New Strategies for a Long-Standing Issue

The cerebellar microcircuit has been the work bench for theoretical and computational modeling since the beginning of neuroscientific research. The regular neural architecture of the cerebellum inspired different solutions to the long-standing issue of how its circuitry could control motor learning and coordination. Originally, the cerebellar network was modeled using a statistical-topological approach that was later extended by considering the geometrical organization of local microcircuits. However, with the advancement in anatomical and physiological investigations, new discoveries have revealed an unexpected richness of connections, neuronal dynamics and plasticity, calling for a change in modeling strategies, so as to include the multitude of elementary aspects of the network into an integrated and easily updatable computational framework. Recently, biophysically accurate “realistic” models using a bottom-up strategy accounted for both detailed connectivity and neuronal non-linear membrane dynamics. In this perspective review, we will consider the state of the art and discuss how these initial efforts could be further improved. Moreover, we will consider how embodied neurorobotic models including spiking cerebellar networks could help explaining the role and interplay of distributed forms of plasticity. We envisage that realistic modeling, combined with closed-loop simulations, will help to capture the essence of cerebellar computations and could eventually be applied to neurological diseases and neurorobotic control systems.

Computational Modeling of Single Neuron Extracellular Electric Potentials and Network Local Field Potentials using LFPsim

Local Field Potentials (LFPs) are population signals generated by complex spatiotemporal interaction of current sources and dipoles. Mathematical computations of LFPs allow the study of circuit functions and dysfunctions via simulations. This paper introduces LFPsim, a NEURON-based tool for computing population LFP activity and single neuron extracellular potentials. LFPsim was developed to be used on existing cable compartmental neuron and network models. Point source, line source, and RC based filter approximations can be used to compute extracellular activity. As a demonstration of efficient implementation, we showcase LFPs from mathematical models of electrotonically compact cerebellum granule neurons and morphologically complex neurons of the neocortical column. LFPsim reproduced neocortical LFP at 8, 32, and 56 Hz via current injection, in vitro post-synaptic N_2a, N_2b waves and in vivo T-C waves in cerebellum granular layer. LFPsim also includes a simulation of multi-electrode array of LFPs in network populations to aid computational inference between biophysical activity in neural networks and corresponding multi-unit activity resulting in extracellular and evoked LFP signals.

Oscillation-Driven Spike-Timing Dependent Plasticity Allows Multiple Overlapping Pattern Recognition in Inhibitory Interneuron Networks

The majority of operations carried out by the brain require learning complex signal patterns for future recognition, retrieval and reuse. Although learning is thought to depend on multiple forms of long-term synaptic plasticity, the way this latter contributes to pattern recognition is still poorly understood. Here, we have used a simple model of afferent excitatory neurons and interneurons with lateral inhibition, reproducing a network topology found in many brain areas from the cerebellum to cortical columns. When endowed with spike-timing dependent plasticity (STDP) at the excitatory input synapses and at the inhibitory interneuron–interneuron synapses, the interneurons rapidly learned complex input patterns. Interestingly, induction of plasticity required that the network be entrained into theta-frequency band oscillations, setting the internal phase-reference required to drive STDP. Inhibitory plasticity effectively distributed multiple patterns among available interneurons, thus allowing the simultaneous detection of multiple overlapping patterns. The addition of plasticity in intrinsic excitability made the system more robust allowing self-adjustment and rescaling in response to a broad range of input patterns. The combination of plasticity in lateral inhibitory connections and homeostatic mechanisms in the inhibitory interneurons optimized mutual information (MI) transfer. The storage of multiple complex patterns in plastic interneuron networks could be critical for the generation of sparse representations of information in excitatory neuron populations falling under their control.

Thiamine Deficiency Increases Ca2+ Current and CaV1.2 L-type Ca2+ Channel Levels in Cerebellum Granular Neurons

Thiamine (vitamin B1) is co-factor for three pivotal enzymes for glycolytic metabolism: pyruvate dehydrogenase, α-ketoglutarate dehydrogenase, and transketolase. Thiamine deficiency leads to neurodegeneration of several brain regions, especially the cerebellum. In addition, several neurodegenerative diseases are associated with impairments of glycolytic metabolism, including Alzheimer's disease. Therefore, understanding the link between dysfunction of the glycolytic pathway and neuronal death will be an important step to comprehend the mechanism and progression of neuronal degeneration as well as the development of new treatment for neurodegenerative states. Here, using an in vitro model to study the effects of thiamine deficiency on cerebellum granule neurons, we show an increase in Ca²⁺ current density and Ca_V1.2 expression. These results indicate a link between alterations in glycolytic metabolism and changes to Ca²⁺ dynamics, two factors that have been implicated in neurodegeneration.