henhaobio
Sign in

bioRxivNeurosciencePreprint

The Cerebellar Engine: Multiscale Digital Brain Co-simulations Reveal How Cerebellar Spiking Architecture Shapes Cortical Coherence

Alice GeminianiiD, Jil Mona MeieriD, Dionysios PerdikisiD, Sabrine Ouertani, Claudia CasellatoiD, Petra RitteriD, Egidio Ugo D'AngeloiD

PDFDOICodeData

Abstract

Cellular activities shape large-scale brain dynamics determining brain functioning and disease, yet the causal mechanisms across scales remain unclear. In particular, the cerebellum has been reported to modulate whole-brain dynamics during sensorimotor integration through unknown circuit interactions. To investigate the underlying mechanisms, we developed a novel multiscale digital brain simulator, in which a spiking neural network of the olivocerebellar microcircuit is embedded in a mean-field virtual mouse brain and wired using an atlas-based long-range connectome. Parameters were systematically tuned to match multiscale experimental data from primary sensory and motor cortices (S1 and M1) and cerebellum. We analyzed the role of cerebellar circuitry on sensorimotor integration by lesioning critical circuit connections in silico. Results suggested that Purkinje cell inhibition enhances the processing efficiency of the 'cerebellar engine' through decorrelation of cerebellar nuclei activity and that the pathway between mossy fibers and cerebellar nuclei is the specific pathway inside the microcircuit driving M1-S1 coherence. These results indicate a mechanistic link between cerebellar microcircuit and cortical sensorimotor processing. This novel framework opens new perspectives for the broader multiscale investigation of brain physiological and pathological states in relation to specific cellular and microcircuit properties.

In-depth summary

Claude Opus 5.5 · xhigh effort · 8 Nov 2026 · from the full text of version 3 (5 Sept 2026) · CC BY 4.0

In a mouse whole-brain co-simulation with a spiking olivocerebellar network, cutting mossy-fiber input to cerebellar nuclei lowers M1-S1 gamma coherence, while removing Purkinje inhibition raises it.

Inactivating the cerebellar nuclei has been reported to reduce gamma coherence between primary motor (M1) and somatosensory (S1) cortex during whisking, but which cerebellar microcircuit elements drive this is hard to isolate in vivo. According to the authors, earlier multiscale co-simulations were unidirectional and focused on hippocampus or basal ganglia; this work instead embeds a spiking cerebellum in a closed loop with a whole-brain neural-mass model, so cross-scale causal effects can be tested.

Key results

  1. SuggestedRemoving mossy fiber excitation of the cerebellar nuclei (CN) lowered input-output coherence (principal trigeminal nucleus, PrV, to CN) and lowered M1-S1 gamma coherence to about the level seen with the whole cerebellum inactivated. Removing Purkinje cell (PC) inhibition of the CN raised both.

    Virtual lesions (MOStoCN, PCtoCN) in TVMB-NEST co-simulations, 10 runs per condition vs CONTROL, two-sided unpaired t-tests with FDR correction (Fig. 4A-B). Both lesions lowered the gamma-to-theta coherence ratio (Fig. 4C-D). Gamma coherence trends held without normalizing cerebellar output, but the ratio effect was then masked (Supp. Fig. S5). · 2.3.2 Cerebellar microcircuit contributions to…; 2.3.3 Cerebellar circuit elements differentially…; 4.6 Microcircuit lesion co-simulations; Figure legends

  2. SuggestedAt the spike level, PC inhibition appears to decorrelate nuclear output from mossy fiber input. Synchrony between input and output fell with MOStoCN and rose with PCtoCN, while phase locking of nuclear spikes to the input gamma rhythm dropped with both lesions.

    The phase-locking value (PLV) of excitatory CN neuron (CNe) spikes to the gamma-filtered mossy fiber rate fell significantly in both lesions (Fig. 6A). Firing-rate cross-correlation and pairwise spike-train correlation between mossy fibers and CNe fell with MOStoCN and rose with PCtoCN (Fig. 6B-C). Gamma-band cross-spectral density magnitude fell with both lesions (p < 0.001, FDR-corrected; Fig. 5B). · 2.3.3 Cerebellar circuit elements differentially…; 2.3.4 Changes in cerebellar spiking activity…; 3. Discussion; Figure legends

  3. SuggestedGamma power and coherence came apart. Both CN-input lesions cut CNe gamma power relative to mossy fibers, yet PCtoCN raised coherence. Lesioning molecular layer interneuron (MLI) inhibition altered PC gamma dynamics without changing M1-S1 coherence.

    Relative CNe gamma power fell significantly with both MOStoCN and PCtoCN (Fig. 3D). Relative PC gamma power fell only with MOStoCN (Fig. 3C), which the authors attribute to recurrent loops through the rest of the brain. MLItoPC lowered coherence only within the cerebellum and to a lesser extent; MLItoMLI left coherence unchanged (Table 1; Supplementary Note S2). · 2.3.1 Cerebellar microcircuit contributions to…; 2.3.2 Cerebellar microcircuit contributions to…; 2.3.3 Cerebellar circuit elements differentially…; 3. Discussion

  4. ShownAfter fitting, the model reproduced the experimental drop in M1-S1 gamma coherence when the cerebellum is inactivated. This held both with neural masses alone and with the spiking cerebellum embedded.

    Gamma-band (25 to 60 Hz) coherence drop, computed as (cerebON - cerebOFF)/cerebON: mean 53.4% for TVMB-only and 54.8% for TVMB-NEST co-simulations vs 51.3% in the experimental data. cerebON vs cerebOFF was significant (two-sided t-test, p < 0.0001; 10 runs). The pathway gains were fitted to this drop using TVMB-only runs. · 2.2 The model captures experimental data across…; 4.7 Analysis and Statistics; Figure legends

  5. ShownAll spiking cerebellar populations fired at physiological rates and carried whisking-related theta oscillations, although only the granule cell rate was fitted. The network was never tuned on oscillations, yet gamma power in both PCs and nuclear neurons was higher than in mossy fibers.

    The interface weight was set by grid search to a granule cell target of 4.2 Hz. Mean rates of the other populations over 10 runs fell within experimental ranges from rodent single-unit recordings during free whisking (Fig. 2F). Relative to mossy fibers, PC gamma and theta power rose; CNe gamma power rose and theta power fell (Fig. 3A-B). · 2.2 The model captures experimental data across…; 4.4 Inter-scale interface models; Figure legends

Methods, data and reproducibility

A 214-region mouse connectome (Allen Mouse Brain Connectivity Atlas) of Wilson-Cowan neural masses in The Virtual Mouse Brain (TVMB), with two whisker nodes, was co-simulated with a NEST spiking olivocerebellar network. Parameters were fitted to M1/S1 power spectra, the coherence drop on cerebellar inactivation and granule cell rates; four spiking-network synapse classes were then lesioned.

  • Allen Mouse Brain Connectivity Atlas connectome; one hemisphere mirrored, then non-whisking regions merged into a 214-region whole-brain network · mouse · n = 298 areas per hemisphere (596 regions), merged to 214
  • Local field potential power spectra from M1 and S1 during rodent whisking, taken from the literature as fitting targets · rodent (species not specified)
  • Experimental M1-S1 gamma coherence drop after cerebellar nuclei inactivation during free whisking (prior study), used to fit pathway gains · rodent (species not specified)
  • Single-unit firing rates of cerebellar granule cells, MLIs, Purkinje cells and nuclear neurons during free whisking, from the literature · rodent (species not specified)
  • Simulation outputs: TVMB-only and TVMB-NEST runs (cerebON/cerebOFF) and lesion co-simulations (CONTROL plus four lesions) · mouse (in silico model) · n = 10 runs per condition, both hemispheres · Zenodo record 20530965
Methods & tools
The Virtual Mouse Brain (TVMB), The Virtual Brain, NEST 3, Allen Mouse Brain Connectivity Atlas, Allen Software Development Kit, cortico-thalamic Wilson-Cowan neural mass model, extended-generalized leaky integrate-and-fire neuron model, simulation-based inference (SBI), Docker
Software & environment
Docker image dionperd/tvb-multiscale-dev, tag rising-net_publication, pinned by sha256 digest; spiking network in NEST3, neural masses in TVMB (an extension of The Virtual Brain).
Assay conditions & compute
Simulations and SBI fitting ran on a high-performance computing cluster at the Berlin Institute of Health (Berlin, Germany).
Reproducibility gaps
No node counts, run times or memory reported for co-simulations or SBI fitting; the experimental target datasets are cited from earlier studies without accessions.
Limitations (6)
  • The connectome routes cerebello-thalamo-cortical loops only to sensorimotor areas, following the Allen atlas. Cerebellar projections to associative regions such as prefrontal cortex, which the authors say dominate in humans, are not represented. · stated by the authors
  • Whisking is driven through a minimal pathway gain. Parallel sensorimotor loops are not modelled, so how distributed network architecture shapes cerebello-cortical synchronization remains open. · stated by the authors
  • Only free whisking was simulated. Reflexive whisking, where climbing fibers also matter, and MLI effects in frequency bands other than theta and gamma were not examined. · stated by the authors
  • Because the scales form a closed loop, parameters at different scales are interdependent, which the authors say made fitting harder than for single-scale models. · stated by the authors
  • Every mechanistic conclusion rests on in silico lesions; none of the MOStoCN or PCtoCN predictions is tested experimentally. The coherence drop used for validation was also the fitting target for the pathway gains. · noted in this summary
  • Significance comes from t-tests over 10 stochastic runs per condition, so p-values reflect simulation noise rather than biological variability. Spike-correlation bins are 10 ms in the Methods but 20 ms in the Fig. 6C legend. · noted in this summary

Read this if you build multiscale brain models that couple neural masses with spiking networks, or study how inputs to the cerebellar nuclei shape cortical gamma coherence during whisking.

Citations
2

Reproducibility ledger

Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.

Sign in to file a reproducibility report.

Similarity Map
SelectedSame fieldCross-field
THISLindeman '20Gruver '24Lorenzi '22Streng '23Hur '23Yan '24
75% semantic matchNeuroscience
Sander Lindeman, Sungho Hong, Lieke Kros+6 more·Proceedings of the National Academy of Sciences·Dec 22, 2020· 72 citations