Cerebellar Golgi cells (GoCs) efficiently control the spiking activity of granule cells through GABAA receptor-mediated tonic and phasic inhibition. Recent experiments provided compelling evidence for the extensive interconnection of GoCs through electrical synapses, but their chemical inhibitory synaptic inputs are debated. Here, we investigated the GABAergic synaptic inputs of GoCs using in vitro electrophysiology and quantitative light microscopy (LM) and electron microscopy (EM). We characterized GABAA receptor-mediated IPSCs in GoCs and Lugaro cells (LuCs), and found that IPSCs in GoCs have lower frequencies, smaller amplitudes, and much slower decay kinetics. Pharmacological and LM immunolocalization experiments revealed that GoCs express α3, whereas LuCs express α1 subunit-containing GABAA receptors. The selective expression and clustered distribution of the α3 subunit in GoCs allowed the quantitative analysis of GABAergic synapses on their dendrites in the molecular layer (ML). EM and LM experiments in rats, and wild-type and GlyT2-GFP transgenic mice revealed that only one third of axon terminals establishing GABAergic synapses on GoC dendrites contain GlyT2, ruling out LuCs, globular cells, and any noncortical glycinergic inputs as major inhibitory sources. We also show that axon terminals of stellate/basket cells very rarely innervate GlyT2-GFP-expressing GoCs, indicating that only a minority of the inhibitory inputs to GoCs in the ML originates from local interneurons, and the majority of their inhibitory inputs exclusively releases GABA.

NeuroML is an XML-based model description language, which provides a powerful common data format for defining and exchanging models of neurons and neuronal networks. In the latest version of NeuroML, the structure and behavior of ion channel, synapse, cell, and network model descriptions are based on underlying definitions provided in LEMS, a domain-independent language for expressing hierarchical mathematical models of physical entities. While declarative approaches for describing models have led to greater exchange of model elements among software tools in computational neuroscience, a frequent criticism of XML-based languages is that they are difficult to work with directly. Here we describe two Application Programming Interfaces (APIs) written in Python (http://www.python.org), which simplify the process of developing and modifying models expressed in NeuroML and LEMS. The libNeuroML API provides a Python object model with a direct mapping to all NeuroML concepts defined by the NeuroML Schema, which facilitates reading and writing the XML equivalents. In addition, it offers a memory-efficient, array-based internal representation, which is useful for handling large-scale connectomics data. The libNeuroML API also includes support for performing common operations that are required when working with NeuroML documents. Access to the LEMS data model is provided by the PyLEMS API, which provides a Python implementation of the LEMS language, including the ability to simulate most models expressed in LEMS. Together, libNeuroML and PyLEMS provide a comprehensive solution for interacting with NeuroML models in a Python environment.

Two-photon microscopy is widely used to study brain function, but conventional microscopes are too slow to capture the timing of neuronal signalling and imaging is restricted to one plane. Recent development of acousto-optic-deflector-based random access functional imaging has improved the temporal resolution, but the utility of these technologies for mapping 3D synaptic activity patterns and their performance at the excitation wavelengths required to image genetically encoded indicators have not been investigated.

Sustained rate-coded signals encode many types of sensory modalities. Some sensory synapses possess specialized ribbon structures, which tether vesicles, to enable high-frequency signaling. However, central synapses lack these structures, yet some can maintain signaling over a wide bandwidth. To analyze the underlying molecular mechanisms, we investigated the function of the active zone core component Bassoon in cerebellar mossy fiber to granule cell synapses. We show that short-term synaptic depression is enhanced in Bassoon knockout mice during sustained high-frequency trains but basal synaptic transmission is unaffected. Fluctuation and quantal analysis as well as quantification with constrained short-term plasticity models revealed that the vesicle reloading rate was halved in the absence of Bassoon. Thus, our data show that the cytomatrix protein Bassoon speeds the reloading of vesicles to release sites at a central excitatory synapse.

Correlated network activity is important in the development of many neural circuits. Purkinje cells are among the first neurons to populate the cerebellar cortex, where they sprout exuberant axon collaterals. We used multiple patch-clamp recordings targeted with two-photon microscopy to characterize monosynaptic connections between the Purkinje cells of juvenile mice. We found that Purkinje cell axon collaterals projected asymmetrically in the sagittal plane, directed away from the lobule apex. On the basis of our anatomical and physiological characterization of this connection, we constructed a network model that robustly generated waves of activity that traveled along chains of connected Purkinje cells. Consistent with the model, we observed traveling waves of activity in Purkinje cells in sagittal slices from young mice that require GABA(A) receptor-mediated transmission and intact Purkinje cell axon collaterals. These traveling waves are absent in adult mice, suggesting they have a developmental role in wiring the cerebellar cortical microcircuit.

Active invasion of the dendritic tree by action potentials (APs) generated in the axon is essential for associative synaptic plasticity and neuronal ensemble formation. In cortical pyramidal cells (PCs), this AP back-propagation is supported by dendritic voltage-gated Na+ (Nav) channels, whose molecular identity is unknown. Using a highly sensitive electron microscopic immunogold technique, we revealed the presence of the Nav1.6 subunit in hippocampal CA1 PC proximal and distal dendrites. Here, the subunit density is lower by a factor of 35 to 80 than that found in axon initial segments. A gradual decrease in Nav1.6 density along the proximodistal axis of the dendritic tree was also detected without any labeling in dendritic spines. Our results reveal the characteristic subcellular distribution of the Nav1.6 subunit, identifying this molecule as a key substrate enabling dendritic excitability.

Conductance-based neuronal network models can help us understand how synaptic and cellular mechanisms underlie brain function. However, these complex models are difficult to develop and are inaccessible to most neuroscientists. Moreover, even the most biologically realistic network models disregard many 3D anatomical features of the brain. Here, we describe a new software application, neuroConstruct, that facilitates the creation, visualization, and analysis of networks of multicompartmental neurons in 3D space. A graphical user interface allows model generation and modification without programming. Models within neuroConstruct are based on new simulator-independent NeuroML standards, allowing automatic generation of code for NEURON or GENESIS simulators. neuroConstruct was tested by reproducing published models and its simulator independence verified by comparing the same model on two simulators. We show how more anatomically realistic network models can be created and their properties compared with experimental measurements by extending a published 1D cerebellar granule cell layer model to 3D.

The nanoscale topographical arrangement of voltage-gated calcium channels (VGCC) and synaptic vesicles (SVs) determines synaptic strength and plasticity, but whether distinct spatial distributions underpin diversity of synaptic function is unknown. We performed single bouton Ca2+ imaging, Ca2+ chelator competition, immunogold electron microscopic (EM) localization of VGCCs and the active zone (AZ) protein Munc13-1, at two cerebellar synapses. Unexpectedly, we found that weak synapses exhibited 3-fold more VGCCs than strong synapses, while the coupling distance was 5-fold longer. Reaction-diffusion modeling could explain both functional and structural data with two strikingly different nanotopographical motifs: strong synapses are composed of SVs that are tightly coupled (∼10 nm) to VGCC clusters, whereas at weak synapses VGCCs were excluded from the vicinity (∼50 nm) of docked vesicles. The distinct VGCC-SV topographical motifs also confer differential sensitivity to neuromodulation. Thus, VGCC-SV arrangements are not canonical, and their diversity could underlie functional heterogeneity across CNS synapses.

Computational models are powerful tools for exploring the properties of complex biological systems. In neuroscience, data-driven models of neural circuits that span multiple scales are increasingly being used to understand brain function in health and disease. But their adoption and reuse has been limited by the specialist knowledge required to evaluate and use them. To address this, we have developed Open Source Brain, a platform for sharing, viewing, analyzing, and simulating standardized models from different brain regions and species. Model structure and parameters can be automatically visualized and their dynamical properties explored through browser-based simulations. Infrastructure and tools for collaborative interaction, development, and testing are also provided. We demonstrate how existing components can be reused by constructing new models of inhibition-stabilized cortical networks that match recent experimental results. These features of Open Source Brain improve the accessibility, transparency, and reproducibility of models and facilitate their reuse by the wider community.

Whole-cell voltage recordings were obtained from 64 synaptically coupled excitatory layer 4 (L4) spiny neurones and L2/3 pyramidal cells in acute slices of the somatosensory cortex ('barrel' cortex) of 17- to 23-days-old rats. Single action potentials (APs) in the L4 spiny neurone evoked single unitary EPSPs in the L2/3 pyramidal cell with a peak amplitude of 0.7 +/- 0.6 mV. The average latency was 2.1 +/- 0.6 ms, the rise time was 0.8 +/- 0.3 ms and the decay time constant was 12.7 +/- 3.5 ms. The percentage of failures of an AP in a L4 spiny neurone to evoke a unitary EPSP in the L2/3 pyramidal cell was 4.9 +/- 8.8 % and the coefficient of variation (c.v.) of the unitary EPSP amplitude was 0.27 +/- 0.13. Both c.v. and percentage of failures decreased with increased average EPSP amplitude. Postsynaptic glutamate receptors (GluRs) in L2/3 pyramidal cells were of the N-methyl-D-aspartate (NMDA) receptor (NMDAR) and the non-NMDAR type. At -60 mV in the presence of extracellular Mg2+ (1 mM), 29 +/- 15 % of the EPSP voltage-time integral was blocked by NMDAR antagonists. In 0 Mg2+, the NMDAR/AMPAR ratio of the EPSC was 0.50 +/- 0.29, about half the value obtained for L4 spiny neurone connections. Burst stimulation of L4 spiny neurones showed that EPSPs in L2/3 pyramidal cells depressed over a wide range of frequencies (1-100 s(-1) ). However, at higher frequencies (30 s(-1)) EPSP summation overcame synaptic depression so that the summed EPSP was larger than the first EPSP amplitude in the train. The number of putative synaptic contacts established by the axonal collaterals of the L4 projection neurone with the target neurone in layer 2/3 varied between 4 and 5, with an average of 4.5 +/- 0.5 (n = 13 pairs). Synapses were established on basal dendrites of the pyramidal cell. Their mean geometric distance from the pyramidal cell soma was 67 +/- 34 microm (range, 16-196 microm). The results suggest that each connected L4 spiny neurone produces a weak but reliable EPSP in the pyramidal cell. Therefore transmission of signals to layer 2/3 is likely to have a high threshold requiring simultaneous activation of many L4 neurons, implying that L4 spiny neurone to L2/3 pyramidal cell synapses act as a gate for the lateral spread of excitation in layer 2/3.

Synapses exhibit different short-term plasticity patterns and this behaviour influences information processing in neuronal networks. We tested how the short-term plasticity of excitatory postsynaptic currents (EPSCs) depends on the postsynaptic cell type, identified by axonal arborizations and molecular markers in the hippocampal CA1 area. Three distinct types of short-term synaptic behaviour (facilitating, depressing and combined facilitating-depressing) were defined by fitting a dynamic neurotransmission model to the data. Approximately 75 % of the oriens-lacunosum-moleculare (O-LM) interneurones received facilitating EPSCs, but in three of 12 O-LM cells EPSCs also showed significant depression. Over 90 % of the O-LM cells were immunopositive for somatostatin and mGluR1alpha and all tested cells were decorated by strongly mGluR7a positive axon terminals. Responses in eight of 12 basket cells were described well with a model involving only depression, but the other cells displayed combined facilitating-depressing EPSCs. No apparent difference was found between the plasticity of EPSCs in cholecystokinin- or parvalbumin-containing basket cells. In oriens-bistratified cells (O-Bi), two of nine cells showed facilitating EPSCs, another two depressing, and the remaining five cells combined facilitating-depressing EPSCs. Seven of 10 cells tested for somatostatin were immunopositive, but mGluR1alpha was detectable only in two of 11 tested cells. Furthermore, most O-Bi cells projected to the CA3 area and the subiculum, as well as outside the hippocampal formation. Postsynaptic responses to action potentials recorded in vivo from a CA1 place cell were modelled, and revealed great differences between and within cell types. Our results demonstrate that the short-term plasticity of EPSCs is cell type dependent, but with significant heterogeneity within all three interneurone populations.

Synaptic efficacy and precision are influenced by the coupling of voltage-gated Ca(2+) channels (VGCCs) to vesicles. But because the topography of VGCCs and their proximity to vesicles is unknown, a quantitative understanding of the determinants of vesicular release at nanometer scale is lacking. To investigate this, we combined freeze-fracture replica immunogold labeling of Cav2.1 channels, local [Ca(2+)] imaging, and patch pipette perfusion of EGTA at the calyx of Held. Between postnatal day 7 and 21, VGCCs formed variable sized clusters and vesicular release became less sensitive to EGTA, whereas fixed Ca(2+) buffer properties remained constant. Experimentally constrained reaction-diffusion simulations suggest that Ca(2+) sensors for vesicular release are located at the perimeter of VGCC clusters (<30 nm) and predict that VGCC number per cluster determines vesicular release probability without altering release time course. This "perimeter release model" provides a unifying framework accounting for developmental changes in both synaptic efficacy and time course.

The release of GABA from cholecystokinin-containing interneurons is modulated by type-1 cannabinoid receptors (CB1). Here we tested the hypothesis that the strength of CB1-mediated modulation of GABA release is related to the CB1 content of axon terminals. Basket cell boutons have on average 78% higher CB1 content than those of dendritic-layer-innervating (DLI) cells, a consequence of larger bouton surface and higher CB1 density. The CB1 antagonist AM251 caused a 54% increase in action potential-evoked [Ca(2+)] in boutons of basket cells, but not in DLI cells. However, the effect of AM251 did not correlate with CB1 immunoreactivity of individual boutons. Moreover, a CB1 agonist decreased [Ca(2+)] in a cell type- and CB1-content-independent manner. Replica immunogold labelling demonstrated the colocalization of CB1 with the Cav2.2 Ca(2+) channel subunit. Our data suggest that only a subpopulation of CB1s, within nanometre distances from their target Cav2.2 channels, are responsible for endocannabinoid-mediated modulation of GABA release.

In many brain regions, hyperpolarization-activated cationic currents (Ih) are involved in the generation of rhythmic activities, but the role of Ih in olfactory oscillations remains unclear. Knowledge of the cellular and subcellular distributions of hyperpolarization-activated and cyclic nucleotide-gated channel (HCN) subunits is necessary for understanding the role of Ih in olfactory network activities. Using light microscopic immunocytochemistry, we demonstrate strong HCN1 labelling of the glomerular layer and moderate staining of granule cell, internal and external plexiform layers of the rat main olfactory bulb. In the glomerular layer, among many unlabelled neurons, two distinct subpopulations of juxtaglomerular cells are labelled. Approximately 10% of the juxtaglomerular cells strongly express HCN1. These small diameter cells are immunoreactive for GABA and comprise a subpopulation of periglomerular cells. An additional subset of juxtaglomerular cells ( approximately 1%) expresses low levels of HCN1. They are large in diameter, GABA immunonegative but immunopositive for vesicular glutamate transporter 2, characterizing them as external tufted cells. Quantitative immunogold localization revealed that the somatic plasma membranes of periglomerular cells contain approximately four times more HCN1 labelling than those of external tufted cells. Unlike in cortical pyramidal cells, immunogold density for HCN1 does not significantly differ in somatic and dendritic plasma membranes of external tufted cells, indicating that post-synaptic potentials arriving at proximal and distal dendrites are modulated by the same density of Ih. Our results demonstrate a cell type-dependent expression of HCN1 in the olfactory bulb and predict a differential contribution of distinct juxtaglomerular cell types to network oscillations.

Significant inroads have been made to understand cerebellar cortical processing but neural coding at the output stage of the cerebellum in the deep cerebellar nuclei (DCN) remains poorly understood. The DCN are unlikely to just present a relay nucleus because Purkinje cell inhibition has to be turned into an excitatory output signal, and DCN neurons exhibit complex intrinsic properties. In particular, DCN neurons exhibit a range of rebound spiking properties following hyperpolarizing current injection, raising the question how this could contribute to signal processing in behaving animals. Computer modeling presents an ideal tool to investigate how intrinsic voltage-gated conductances in DCN neurons could generate the heterogeneous firing behavior observed, and what input conditions could result in rebound responses. To enable such an investigation we built a compartmental DCN neuron model with a full dendritic morphology and appropriate active conductances. We generated a good match of our simulations with DCN current clamp data we recorded in acute slices, including the heterogeneity in the rebound responses. We then examined how inhibitory and excitatory synaptic input interacted with these intrinsic conductances to control DCN firing. We found that the output spiking of the model reflected the ongoing balance of excitatory and inhibitory input rates and that changing the level of inhibition performed an additive operation. Rebound firing following strong Purkinje cell input bursts was also possible, but only if the chloride reversal potential was more negative than -70 mV to allow de-inactivation of rebound currents. Fast rebound bursts due to T-type calcium current and slow rebounds due to persistent sodium current could be differentially regulated by synaptic input, and the pattern of these rebounds was further influenced by HCN current. Our findings suggest that active properties of DCN neurons could play a crucial role for signal processing in the cerebellum.

Recently developed pharmacogenetic and optogenetic approaches, with their own advantages and disadvantages, have become indispensable tools in modern neuroscience. Here, we employed a previously described knock-in mouse line (GABA(A)Rγ2(77I)lox) in which the γ2 subunit of the GABA(A) receptor (GABA(A)R) was mutated to become zolpidem insensitive (γ2(77I)) and used viral vectors to swap γ2(77I) with wild-type, zolpidem-sensitive γ2 subunits (γ2(77F)). The verification of unaltered density and subcellular distribution of the virally introduced γ2 subunits requires their selective labelling. For this we generated six N- and six C-terminal-tagged γ2 subunits, with which cortical cultures of GABA(A)Rγ2(−/−) mice were transduced using lentiviruses. We found that the N-terminal AU1 tag resulted in excellent immunodetection and unimpaired synaptic localization. Unaltered kinetic properties of the AU1-tagged γ2 ((AU1)γ2(77F)) channels were demonstrated with whole-cell patch-clamp recordings of spontaneous IPSCs from cultured cells. Next, we carried out stereotaxic injections of lenti- and adeno-associated viruses containing Cre-recombinase and the (AU1)γ2(77F) subunit (Cre-2A-(AU1)γ2(77F)) into the neocortex of GABA(A)Rγ2(77I)lox mice. Light microscopic immunofluorescence and electron microscopic freeze-fracture replica immunogold labelling demonstrated the efficient immunodetection of the AU1 tag and the normal enrichment of the (AU1)γ2(77F) subunits in perisomatic GABAergic synapses. In line with this,miniature and action potential-evoked IPSCs whole-cell recorded from transduced cells had unaltered amplitudes, kinetics and restored zolpidem sensitivity. Our results obtained with a wide range of structural and functional verification methods reveal unaltered subcellular distributions and functional properties of γ2(77I) and (AU1)γ2(77F) GABA(A)Rs in cortical pyramidal cells. This transgenic–viral pharmacogenetic approach has the advantage that it does not require any extrinsic protein that might endow some unforeseen alterations of the genetically modified cells. In addition, this virus-based approach opens up the possibility of modifying multiple cell types in distinct brain regions and performing alternative recombination-based intersectional genetic manipulations.

Local circuit GABAergic interneurons comprise the most diverse cell populations of neuronal networks. Interneurons have been characterized and categorized based on their axo-somato-dendritic morphologies, neurochemical content, intrinsic electrical properties and their firing in relation to in-vivo population activity. Great advances in our understanding of their roles have been facilitated by their selective identification. Recently, we have described three major subtypes of deep short-axon cells (dSACs) of the main olfactory bulb (MOB) based on their axo-dendritic distributions and synaptic connectivity. Here, we investigated whether dSACs also display pronounced molecular diversity and whether distinct dSAC subtypes selectively express certain molecules. Multiple immunofluorescent labeling revealed that the most commonly used molecular markers of dSACs (e.g. vasoactive intestinal polypeptide, calbindin and nitric oxide synthase) label only very small subpopulations (< 7%). In contrast, voltage-gated potassium channel subunits Kv2.1, Kv3.1b, Kv4.3 and the GABA(A) receptor alpha1 subunit are present in 70-95% of dSACs without showing any dSAC subtype-selective expression. However, metabotropic glutamate receptor type 1alpha mainly labels dSACs that project to the glomerular layer (GL-dSAC subtype) and comprise approximately 20% of the total dSAC population. Analysing these molecular markers with stereological methods, we estimated the total number of dSACs in the entire MOB to be approximately 13,500, which is around a quarter of the number of mitral cells. Our results demonstrate a large molecular heterogeneity of dSACs and reveal a unique neurochemical marker for one dSAC subtype. Based on our results, dSAC subtype-specific genetic modifications will allow us to decipher the role of GL-dSACs in shaping the dynamic activity of the MOB network.

Potassium channels comprise the most diverse family of ion channels and play critical roles in a large variety of physiological and pathological processes. In addition to their molecular diversity, variations in their distributions and densities on the axo-somato-dendritic surface of neurons are key parameters in determining their functional impact. Despite extensive electrophysiological and anatomical investigations, the exact location and densities of most K(+) channels in small subcellular compartments are still unknown. Here we aimed at providing a quantitative surface map of two delayed-rectifier (Kv1.1 and Kv2.1) and one G-protein-gated inwardly rectifying (Kir3.2) K(+) channel subunits on hippocampal CA1 pyramidal cells (PCs). Freeze-fracture replica immunogold labelling was employed to determine the relative densities of these K(+) channel subunits in 18 axo-somato-dendritic compartments. Significant densities of the Kv1.1 subunit were detected on axon initial segments (AISs) and axon terminals, with an approximately eight-fold lower density in the latter compartment. The Kv2.1 subunit was found in somatic, proximal dendritic and AIS plasma membranes at approximately the same densities. This subunit has a non-uniform plasma membrane distribution; Kv2.1 clusters are frequently adjacent to, but never overlap with, GABAergic synapses. A quasi-linear increase in the Kir3.2 subunit density along the dendrites of PCs was detected, showing no significant difference between apical dendritic shafts, oblique dendrites or dendritic spines at the same distance from the soma. Our results demonstrate that each subunit has a unique cell-surface distribution pattern, and predict their differential involvement in synaptic integration and output generation at distinct subcellular compartments.

In vivo, cortical pyramidal cells are bombarded by asynchronous synaptic input arising from ongoing network activity. However, little is known about how such 'background' synaptic input interacts with nonlinear dendritic mechanisms. We have modified an existing model of a layer 5 (L5) pyramidal cell to explore how dendritic integration in the apical dendritic tuft could be altered by the levels of network activity observed in vivo. Here we show that asynchronous background excitatory input increases neuronal gain and extends both temporal and spatial integration of stimulus-evoked synaptic input onto the dendritic tuft. Addition of fast and slow inhibitory synaptic conductances, with properties similar to those from dendritic targeting interneurons, that provided a 'balanced' background configuration, partially counteracted these effects, suggesting that inhibition can tune spatio-temporal integration in the tuft. Excitatory background input lowered the threshold for NMDA receptor-mediated dendritic spikes, extended their duration and increased the probability of additional regenerative events occurring in neighbouring branches. These effects were also observed in a passive model where all the non-synaptic voltage-gated conductances were removed. Our results show that glutamate-bound NMDA receptors arising from ongoing network activity can provide a powerful spatially distributed nonlinear dendritic conductance. This may enable L5 pyramidal cells to change their integrative properties as a function of local network activity, potentially allowing both clustered and spatially distributed synaptic inputs to be integrated over extended timescales.

Many theories of cerebellar function assume that long-term depression (LTD) of parallel fiber (PF) synapses enables Purkinje cells to learn to recognize PF activity patterns. We have studied the LTD-based recognition of PF patterns in a biophysically realistic Purkinje-cell model. With simple-spike firing as observed in vivo, the presentation of a pattern resulted in a burst of spikes followed by a pause. Surprisingly, the best criterion to distinguish learned patterns was the duration of this pause. Moreover, our simulations predicted that learned patterns elicited shorter pauses, thus increasing Purkinje-cell output. We tested this prediction in Purkinje-cell recordings both in vitro and in vivo. In vitro, we found a shortening of pauses when decreasing the number of active PFs or after inducing LTD. In vivo, we observed longer pauses in LTD-deficient mice. Our results suggest a novel form of neural coding in the cerebellar cortex.