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

A multidimensional investigation of sleep and biopsychosocial profiles with associated neural signatures

Sleep is essential for optimal functioning and health. Interconnected to multiple biological, psychological and socio-environmental factors (i.e., biopsychosocial factors), the multidimensional nature of sleep is rarely capitalized on in research. Here, we deployed a data-driven approach to identify sleep-biopsychosocial profiles that linked self-reported sleep patterns to inter-individual variability in health, cognition, and lifestyle factors in 770 healthy young adults. We uncovered five profiles, including two profiles reflecting general psychopathology associated with either reports of general poor sleep or an absence of sleep complaints (i.e., sleep resilience) respectively. The three other profiles were driven by the use of sleep aids and social satisfaction, sleep duration and cognitive performance, and sleep disturbance linked to cognition and mental health. Furthermore, identified sleep-biopsychosocial profiles displayed unique patterns of brain network organization. In particular, somatomotor network connectivity alterations were involved in the relationships between sleep and biopsychosocial factors. These profiles can potentially untangle the interplay between individuals' variability in sleep, health, cognition and lifestyle - equipping research and clinical settings to better support individual's well-being.

Two Signaling Modes Are Better than One: Flux-Independent Signaling by Ionotropic Glutamate Receptors Is Coming of Age

Glutamate is the major excitatory neurotransmitter in the central nervous system. Glutamatergic transmission can be mediated by ionotropic glutamate receptors (iGluRs), which mediate rapid synaptic depolarization that can be associated with Ca2+ entry and activity-dependent change in the strength of synaptic transmission, as well as by metabotropic glutamate receptors (mGluRs), which mediate slower postsynaptic responses through the recruitment of second messenger systems. A wealth of evidence reported over the last three decades has shown that this dogmatic subdivision between iGluRs and mGluRs may not reflect the actual physiological signaling mode of the iGluRs, i.e., α-amino-3-hydroxy-5-methyl-4-isoxasolepropionic acid (AMPA) receptors (AMPAR), kainate receptors (KARs), and N-methyl-D-aspartate (NMDA) receptors (NMDARs). Herein, we review the evidence available supporting the notion that the canonical iGluRs can recruit flux-independent signaling pathways not only in neurons, but also in brain astrocytes and cerebrovascular endothelial cells. Understanding the signaling versatility of iGluRs can exert a profound impact on our understanding of glutamatergic synapses. Furthermore, it may shed light on novel neuroprotective strategies against brain disorders.

Data-driven multiscale computational models of cortical and subcortical regions

Data-driven computational models of neurons, synapses, microcircuits, and mesocircuits have become essential tools in modern brain research. The goal of these multiscale models is to integrate and synthesize information from different levels of brain organization, from cellular properties, dendritic excitability, and synaptic dynamics to microcircuits, mesocircuits, and ultimately behavior. This article surveys recent advances in the genesis of data-driven computational models of mammalian neural networks in cortical and subcortical areas. I discuss the challenges and opportunities in developing data-driven multiscale models, including the need for interdisciplinary collaborations, the importance of model validation and comparison, and the potential impact on basic and translational neuroscience research. Finally, I highlight future directions and emerging technologies that will enable more comprehensive and predictive data-driven models of brain function and dysfunction.

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.

Orthogonal IMiD-Degron Pairs Induce Selective Protein Degradation in Cells

Immunomodulatory imide drugs (IMiDs) including thalidomide, lenalidomide, and pomalidomide, can be used to induce degradation of a protein of interest that is fused to a short zinc finger (ZF) degron motif. These IMiDs, however, also induce degradation of endogenous neosubstrates, including IKZF1 and IKZF3. To improve degradation selectivity, we took a bump-and-hole approach to design and screen bumped IMiD analogs against 8380 ZF mutants. This yielded a bumped IMiD analog that induces efficient degradation of a mutant ZF degron, while not affecting other cellular proteins, including IKZF1 and IKZF3. In proof-of-concept studies, this system was applied to induce efficient degradation of TRIM28, a disease-relevant protein with no known small molecule binders. We anticipate that this system will make a valuable addition to the current arsenal of degron systems for use in target validation.

Cell-cycle and Age-Related Modulations in Mouse Chromosome Stiffness

The intricate structure of chromosomes is complex, and many aspects of chromosome configuration/organization remain to be fully understood. Measuring chromosome stiffness can provide valuable insights into their structure. However, the nature of chromosome stiffness, whether static or dynamic, remains elusive. In this study, we analyzed chromosome stiffness in MI and MII oocytes. We revealed that MI oocytes had a ten-fold increase in stiffness compared to mitotic chromosomes, whereas chromosome stiffness in MII oocytes was relatively low chromosome. We then investigated the contribution of meiosis-specific cohesin complexes to chromosome stiffness in MI and MII oocytes. Surprisingly, the Young's modulus of chromosomes from the three meiosis-specific cohesin mutants did not exhibit significant differences compared to the wild type, indicating that these proteins may not play a substantial role in determining chromosome stiffness. Additionally, our findings revealed an age-related increase in chromosome stiffness in MI oocytes. Age correlates with elevated DNA damage levels, so we investigated the impact of etoposide-induced DNA damage on chromosome stiffness, discovering a reduction in stiffness in response to such damage in MI oocytes. Overall, our study underscores the dynamic nature of chromosome stiffness, subject to changes influenced by the cell cycle and age.

Human Purkinje cells outperform mouse Purkinje cells in dendritic complexity and computational capacity

Purkinje cells in the cerebellum are among the largest neurons in the brain and have been extensively investigated in rodents. However, their morphological and physiological properties remain poorly understood in humans. In this study, we utilized high-resolution morphological reconstructions and unique electrophysiological recordings of human Purkinje cells ex vivo to generate computational models and estimate computational capacity. An inter-species comparison showed that human Purkinje cell had similar fractal structures but were larger than those of mouse Purkinje cells. Consequently, given a similar spine density (2/μm), human Purkinje cell hosted approximately 7.5 times more dendritic spines than those of mice. Moreover, human Purkinje cells had a higher dendritic complexity than mouse Purkinje cells and usually emitted 2-3 main dendritic trunks instead of one. Intrinsic electro-responsiveness was similar between the two species, but model simulations revealed that the dendrites could process ~6.5 times (n = 51 vs. n = 8) more input patterns in human Purkinje cells than in mouse Purkinje cells. Thus, while human Purkinje cells maintained spike discharge properties similar to those of rodents during evolution, they developed more complex dendrites, enhancing computational capacity.

Methods for decoding cortical gradients of functional connectivity

Macroscale gradients have emerged as a central principle for understanding functional brain organization. Previous studies have demonstrated that a principal gradient of connectivity in the human brain exists, with unimodal primary sensorimotor regions situated at one end and transmodal regions associated with the default mode network and representative of abstract functioning at the other. The functional significance and interpretation of macroscale gradients remains a central topic of discussion in the neuroimaging community, with some studies demonstrating that gradients may be described using meta-analytic functional decoding techniques. However, additional methodological development is necessary to fully leverage available meta-analytic methods and resources and quantitatively evaluate their relative performance. Here, we conducted a comprehensive series of analyses to investigate and improve the framework of data-driven, meta-analytic methods, thereby establishing a principled approach for gradient segmentation and functional decoding. We found that a two-segment solution determined by a k-means segmentation approach and an LDA-based meta-analysis combined with the NeuroQuery database was the optimal combination of methods for decoding functional connectivity gradients. Finally, we proposed a method for decoding additional components of the gradient decomposition. The current work aims to provide recommendations on best practices and flexible methods for gradient-based functional decoding of fMRI data.

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.

Noncoding genetic variation in ISPD distinguishes gamecocks from nongame chickens

Chickens were domesticated >4,000 years ago, probably first for fighting them and only later as a source of food. Fighting chickens, commonly known as gamecocks, continue to be bred throughout the world, but the genetic relationships among geographically diverse gamecocks and with nongame chickens are not known. Here, we sequenced the genomes of 44 geographically diverse gamecocks and of 62 nongame chickens representing a variety of breeds. We combined these sequences with published genomes to generate the most diverse chicken genomes dataset yet assembled, at 307 samples. We found that gamecocks do not form a homogeneous group, yet they share genetic similarities that distinguish them from nongame chickens. Such similarities are likely the result of a common origin before their local diversification into, or mixing with, nongame chickens. Particularly noteworthy is a variant in an intron of ISPD, an extreme outlier present at a frequency of 90% in gamecocks but only 4% in nongame chickens. The ISPD locus has the strongest signal of selection in gamecocks, suggesting it is important for fighting performance. Because ISPD variants that are highly prevalent in gamecocks are still segregating in nongame chickens, selective breeding may help reduce its frequency in farm conditions in which aggression is not a desired trait. Altogether, our work provides genomic resources for agricultural genetics, uncovers a common origin for gamecocks from around the world and what distinguishes them genetically from chickens bred for purposes other than fighting, and points to ISPD as the most important locus related to fighting performance.

pH-responsive upconversion mesoporous silica nanospheres for combined multimodal diagnostic imaging and targeted photodynamic and photothermal cancer therapy

Photodynamic therapy (PDT) and photothermal therapy (PTT) have garnered considerable interest as non-invasive cancer treatment modalities. However, these approaches remain limited by low solubility, poor stability and inefficient targeting of many common photosensitizers (PSs) and photothermal agents (PTAs). To overcome these limitations, we have designed biocompatible and biodegradable tumor-targeted upconversion nanospheres with imaging capabilities. The multifunctional nanospheres consist of a sodium yttrium fluoride core doped with lanthanides (ytterbium, erbium and gadolinium) and bismuth selenide (NaYF 4 :Yb/Er/Gd,Bi 2 Se 3 ) within a mesoporous silica shell that encapsulates a PS, Chlorin e6 (Ce6), in its pores. NaYF 4 :Yb/Er converts deeply penetrating near-infrared (NIR) light to visible light, which excites the Ce6 to generate cytotoxic reactive oxygen species (ROS), while the PTA Bi 2 Se 3 efficiently converts absorbed NIR light to heat. Additionally, Gd enables magnetic resonance imaging (MRI) of the nanospheres. The mesoporous silica shell is coated with lipid/polyethylene glycol (DPPC/cholesterol/DSPE-PEG) to ensure retention of the encapsulated Ce6 and minimize interactions with serum proteins and macrophages that impede tumor targeting. Finally, the coat is functionalized with the acidity-triggered rational membrane (ATRAM) peptide, which promotes specific and efficient internalization into cancer cells within the mildly acidic tumor microenvironment. Following uptake by cancer cells in vitro , NIR laser irradiation of the nanospheres caused substantial cytotoxicity due to ROS production and hyperthermia. The nanospheres facilitated tumor MRI and thermal imaging, and exhibited potent NIR laser light-induced antitumor effects in vivo via combined PDT and PTT, with no observable toxicity to healthy tissue, thereby substantially prolonging survival. Our results demonstrate that the ATRAM-functionalized, lipid/PEG-coated upconversion mesoporous silica nanospheres (ALUMSNs) offer multimodal diagnostic imaging and targeted combinatorial cancer therapy.

Systematically characterizing the roles of E3-ligase family members in inflammatory responses with massively parallel Perturb-seq

E3 ligases regulate key processes, but many of their roles remain unknown. Using Perturb-seq, we interrogated the function of 1,130 E3 ligases, partners and substrates in the inflammatory response in primary dendritic cells (DCs). Dozens impacted the balance of DC1, DC2, migratory DC and macrophage states and a gradient of DC maturation. Family members grouped into co-functional modules that were enriched for physical interactions and impacted specific programs through substrate transcription factors. E3s and their adaptors co-regulated the same processes, but partnered with different substrate recognition adaptors to impact distinct aspects of the DC life cycle. Genetic interactions were more prevalent within than between modules, and a deep learning model, comβVAE, predicts the outcome of new combinations by leveraging modularity. The E3 regulatory network was associated with heritable variation and aberrant gene expression in immune cells in human inflammatory diseases. Our study provides a general approach to dissect gene function.

Model simulations unveil the structure-function-dynamics relationship of the cerebellar cortical microcircuit

The cerebellar network is renowned for its regular architecture that has inspired foundational computational theories. However, the relationship between circuit structure, function and dynamics remains elusive. To tackle the issue, we developed an advanced computational modeling framework that allows us to reconstruct and simulate the structure and function of the mouse cerebellar cortex using morphologically realistic multi-compartmental neuron models. The cerebellar connectome is generated through appropriate connection rules, unifying a collection of scattered experimental data into a coherent construct and providing a new model-based ground-truth about circuit organization. Naturalistic background and sensory-burst stimulation are used for functional validation against recordings in vivo, monitoring the impact of cellular mechanisms on signal propagation, inhibitory control, and long-term synaptic plasticity. Our simulations show how mossy fibers entrain the local neuronal microcircuit, boosting the formation of columns of activity travelling from the granular to the molecular layer providing a new resource for the investigation of local microcircuit computation and of the neural correlates of behavior.

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.

Black-box and surrogate optimization for tuning spiking neural models of striatum plasticity

The basal ganglia (BG) is a brain structure that has long been proposed to play an essential role in action selection, and theoretical models of spiking neurons have tried to explain how the BG solves this problem. A recently proposed functional and biologically inspired network model of the striatum (an important nucleus of the BG) is based on spike-timing-dependent eligibility (STDE) and captured important experimental features of this nucleus. The model can recognize complex input patterns and consistently choose rewarded actions to respond to such sensory inputs. However, model tuning is challenging due to two main reasons. The first is the expert knowledge required, resulting in tedious and potentially biased trial-and-error procedures. The second is the computational cost of assessing model configurations (approximately 1.78 h per evaluation). This study addresses the model tuning problem through numerical optimization. Considering the cost of assessing solutions, the selected methods stand out due to their low requirements for solution evaluations and compatibility with high-performance computing. They are the SurrogateOpt solver of Matlab and the RBFOpt library, both based on radial basis function approximations, and DIRECT-GL, an enhanced version of the widespread black-box optimizer DIRECT. Besides, a parallel random search serves as a baseline reference of the outcome of opting for sophisticated methods. SurrogateOpt turns out to be the best option for tuning this kind of model. It outperforms, on average, the quality of the configuration found by an expert and works significantly faster and autonomously. RBFOpt and the random search share the second position, but their average results are below the option found by hand. Finally, DIRECT-GL follows this line becoming the worst-performing method.

Endoplasmic Reticulum in Metaplasticity: From Information Processing to Synaptic Proteostasis

The ER (endoplasmic reticulum) is a Ca²⁺ reservoir and the unique protein-synthesizing machinery which is distributed throughout the neuron and composed of multiple different structural domains. One such domain is called EMC (endoplasmic reticulum membrane protein complex), pleiotropic nature in cellular functions. The ER/EMC position inside the neurons unmasks its contribution to synaptic plasticity via regulating various cellular processes from protein synthesis to Ca²⁺ signaling. Since presynaptic Ca²⁺ channels and postsynaptic ionotropic receptors are organized into the nanodomains, thus ER can be a crucial player in establishing TMNCs (transsynaptic molecular nanocolumns) to shape efficient neural communications. This review hypothesized that ER is not only involved in stress-mediated neurodegeneration but also axon regrowth, remyelination, neurotransmitter switching, information processing, and regulation of pre- and post-synaptic functions. Thus ER might not only be a protein-synthesizing and quality control machinery but also orchestrates plasticity of plasticity (metaplasticity) within the neuron to execute higher-order brain functions and neural repair.

Pharmacological Modulation and (Patho)Physiological Roles of TRPM4 Channel-Part 2: TRPM4 in Health and Disease

Transient receptor potential melastatin 4 (TRPM4) is a unique member of the TRPM protein family and, similarly to TRPM5, is Ca²⁺ sensitive and permeable for monovalent but not divalent cations. It is widely expressed in many organs and is involved in several functions; it regulates membrane potential and Ca²⁺ homeostasis in both excitable and non-excitable cells. This part of the review discusses the currently available knowledge about the physiological and pathophysiological roles of TRPM4 in various tissues. These include the physiological functions of TRPM4 in the cells of the Langerhans islets of the pancreas, in various immune functions, in the regulation of vascular tone, in respiratory and other neuronal activities, in chemosensation, and in renal and cardiac physiology. TRPM4 contributes to pathological conditions such as overactive bladder, endothelial dysfunction, various types of malignant diseases and central nervous system conditions including stroke and injuries as well as in cardiac conditions such as arrhythmias, hypertrophy, and ischemia-reperfusion injuries. TRPM4 claims more and more attention and is likely to be the topic of research in the future.

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.

Predicting Synaptic Connectivity for Large-Scale Microcircuit Simulations Using Snudda

Simulation of large-scale networks of neurons is an important approach to understanding and interpreting experimental data from healthy and diseased brains. Owing to the rapid development of simulation software and the accumulation of quantitative data of different neuronal types, it is possible to predict both computational and dynamical properties of local microcircuits in a ‘bottom-up’ manner. Simulated data from these models can be compared with experiments and ‘top-down’ modelling approaches, successively bridging the scales. Here we describe an open source pipeline, using the software Snudda, for predicting microcircuit connectivity and for setting up simulations using the NEURON simulation environment in a reproducible way. We also illustrate how to further ‘curate’ data on single neuron morphologies acquired from public databases. This model building pipeline was used to set up a first version of a full-scale cellular level model of mouse dorsal striatum. Model components from that work are here used to illustrate the different steps that are needed when modelling subcortical nuclei, such as the basal ganglia.

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.

Projection-dependent heterogeneity of cerebellar granule cell calcium responses

Cerebellar granule cells (GCs) relay mossy fiber (MF) inputs to Purkinje cell dendrites via their axons, the parallel fibers (PFs), which are individually located at a given sublayer of the molecular layer (ML). Although a certain degree of heterogeneity among GCs has been recently reported, variability of GC responses to MF inputs has never been associated with their most notable structural variability, location of their projecting PFs in the ML. Here, we utilize an adeno-associated virus (AAV)-mediated labeling technique that enables us to categorize GCs according to the location of their PFs, and compare the Ca²⁺ responses to MF stimulations between three groups of GCs, consisting of either GCs having PFs at the deep (D-GCs), middle (M-GCs), or superficial (S-GCs) sublayer. Our structural analysis revealed that there was no correlation between position of GC soma in the GC layer and location of its PF in the ML, confirming that our AAV-mediated labeling was important to test the projection-dependent variability of the Ca²⁺ responses in GCs. We then found that the Ca²⁺ responses of D-GCs differed from those of M-GCs. Pharmacological experiments implied that the different Ca²⁺ responses were mainly attributable to varied distributions of GABA_A receptors (GABA_ARs) at the synaptic and extrasynaptic regions of GC dendrites. In addition to GABA_AR distributions, amounts of extrasynaptic NMDA receptors appear to be also varied, because Ca²⁺ responses were different between D-GCs and M-GCs when glutamate spillover was enhanced. Whereas the Ca²⁺ responses of S-GCs were mostly equivalent to those of D-GCs and M-GCs, the blockade of GABA uptake resulted in larger Ca²⁺ responses in S-GCs compared with D-GCs and M-GCs, implying existence of mechanisms leading to more excitability in S-GCs with increased GABA release. Thus, this study reveals MF stimulation-mediated non-uniform Ca²⁺ responses in the cerebellar GCs associated with the location of their PFs in the ML, and raises a possibility that combination of inherent functional variability of GCs and their specific axonal projection contributes to the information processing through the GCs.

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.

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.

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.