Determinants of synaptic integration and heterogeneity in rebound firing explored with data-driven models of deep cerebellar nucleus cells

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.     

Inducing Dendritic Growth in Cultured Sympathetic Neurons

The shape of the dendritic arbor determines the total synaptic input a neuron can receive ^1-3, and influences the types and distribution of these inputs ^4-6. Altered patterns of dendritic growth and plasticity are associated with impaired neurobehavioral function in experimental models ⁷, and are thought to contribute to clinical symptoms observed in both neurodevelopmental disorders ^8-10 and neurodegenerative diseases ^11-13. Such observations underscore the functional importance of precisely regulating dendritic morphology, and suggest that identifying mechanisms that control dendritic growth will not only advance understanding of how neuronal connectivity is regulated during normal development, but may also provide insight on novel therapeutic strategies for diverse neurological diseases.Mechanistic studies of dendritic growth would be greatly facilitated by the availability of a model system that allows neurons to be experimentally switched from a state in which they do not extend dendrites to one in which they elaborate a dendritic arbor comparable to that of their in vivo counterparts. Primary cultures of sympathetic neurons dissociated from the superior cervical ganglia (SCG) of perinatal rodents provide such a model. When cultured in defined medium in the absence of serum and ganglionic glial cells, sympathetic neurons extend a single process which is axonal, and this unipolar state persists for weeks to months in culture ^14,15. However, the addition of either bone morphogenetic protein-7 (BMP-7) ^16,17 or Matrigel ¹⁸ to the culture medium triggers these neurons to extend multiple processes that meet the morphologic, biochemical and functional criteria for dendrites. Sympathetic neurons dissociated from the SCG of perinatal rodents and grown under defined conditions are a homogenous population of neurons ¹⁹ that respond uniformly to the dendrite-promoting activity of Matrigel, BMP-7 and other BMPs of the decapentaplegic (dpp) and 60A subfamilies ^17,18,20,21. Importantly, Matrigel- and BMP-induced dendrite formation occurs in the absence of changes in cell survival or axonal growth ^17,18.Here, we describe how to set up dissociated cultures of sympathetic neurons derived from the SCG of perinatal rats so that they are responsive to the selective dendrite-promoting activity of Matrigel or BMPs.     

Glutamate-Bound NMDARs Arising from In Vivo-like Network Activity Extend Spatio-temporal Integration in a L5 Cortical Pyramidal Cell Model

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.     

Using Multi-Compartment Ensemble Modeling As an Investigative Tool of Spatially Distributed Biophysical Balances: Application to Hippocampal Oriens-Lacunosum/Moleculare (O-LM) Cells

Multi-compartmental models of neurons provide insight into the complex, integrative properties of dendrites. Because it is not feasible to experimentally determine the exact density and kinetics of each channel type in every neuronal compartment, an essential goal in developing models is to help characterize these properties. To address biological variability inherent in a given neuronal type, there has been a shift away from using hand-tuned models towards using ensembles or populations of models. In collectively capturing a neuron''s output, ensemble modeling approaches uncover important conductance balances that control neuronal dynamics. However, conductances are never entirely known for a given neuron class in terms of its types, densities, kinetics and distributions. Thus, any multi-compartment model will always be incomplete. In this work, our main goal is to use ensemble modeling as an investigative tool of a neuron''s biophysical balances, where the cycling between experiment and model is a design criterion from the start. We consider oriens-lacunosum/moleculare (O-LM) interneurons, a prominent interneuron subtype that plays an essential gating role of information flow in hippocampus. O-LM cells express the hyperpolarization-activated current (I _h). Although dendritic I _h could have a major influence on the integrative properties of O-LM cells, the compartmental distribution of I _h on O-LM dendrites is not known. Using a high-performance computing cluster, we generated a database of models that included those with or without dendritic I _h. A range of conductance values for nine different conductance types were used, and different morphologies explored. Models were quantified and ranked based on minimal error compared to a dataset of O-LM cell electrophysiological properties. Co-regulatory balances between conductances were revealed, two of which were dependent on the presence of dendritic I _h. These findings inform future experiments that differentiate between somatic and dendritic I _h, thereby continuing a cycle between model and experiment.     

Dendrites determine the contribution of after depolarization potentials (ADPs) to generation of repetitive action potentials in hypothalamic gonadotropin releasing-hormone (GnRH) neurons

The impact of structure in modulating synaptic signals originating in dendrites is widely recognized. In this study, we focused on the impact of dendrite morphology on a local spike generating mechanism which has been implicated in hormone secretion, the after depolarization potential (ADP). Using multi-compartmental models of hypothalamic GnRH neurons, we systematically truncated dendrite length and determined the consequence on ADP amplitude and repetitive firing. Decreasing the length of the dendrite significantly increased the amplitude of the ADP and increased repetitive firing. These effects were observed in dendrites both with and without active conductances suggesting they largely reflect passive characteristics of the dendrite. In order to test the findings of the model, we performed whole-cell recordings in GnRH neurons and elicited ADPs using current injection. During recordings, neurons were filled with biocytin so that we could determine dendritic and total projection (dendrite plus axon) length. Neurons exhibited ADPs and increasing ADP amplitude was associated with decreasing dendrite length, in keeping with the predictions of the models. Thus, despite the relatively simple morphology of the GnRH neuron’s dendrite, it can still exert a substantial impact on the final neuronal output. This work was supported by HD-45436 to KJS and by NCRR P20 RR16481 to Nigel Cooper.     

A multi-compartment model for interneurons in the dorsal lateral geniculate nucleus

GABAergic interneurons (INs) in the dorsal lateral geniculate nucleus (dLGN) shape the information flow from retina to cortex, presumably by controlling the number of visually evoked spikes in geniculate thalamocortical (TC) neurons, and refining their receptive field. The INs exhibit a rich variety of firing patterns: Depolarizing current injections to the soma may induce tonic firing, periodic bursting or an initial burst followed by tonic spiking, sometimes with prominent spike-time adaptation. When released from hyperpolarization, some INs elicit rebound bursts, while others return more passively to the resting potential. A full mechanistic understanding that explains the function of the dLGN on the basis of neuronal morphology, physiology and circuitry is currently lacking. One way to approach such an understanding is by developing a detailed mathematical model of the involved cells and their interactions. Limitations of the previous models for the INs of the dLGN region prevent an accurate representation of the conceptual framework needed to understand the computational properties of this region. We here present a detailed compartmental model of INs using, for the first time, a morphological reconstruction and a set of active dendritic conductances constrained by experimental somatic recordings from INs under several different current-clamp conditions. The model makes a number of experimentally testable predictions about the role of specific mechanisms for the firing properties observed in these neurons. In addition to accounting for the significant features of all experimental traces, it quantitatively reproduces the experimental recordings of the action-potential- firing frequency as a function of injected current. We show how and why relative differences in conductance values, rather than differences in ion channel composition, could account for the distinct differences between the responses observed in two different neurons, suggesting that INs may be individually tuned to optimize network operation under different input conditions.     

Activity Correlation Imaging: Visualizing Function and Structure of Neuronal Populations

For the analysis of neuronal networks it is an important yet unresolved task to relate the neurons' activities to their morphology. Here we introduce activity correlation imaging to simultaneously visualize the activity and morphology of populations of neurons. To this end we first stain the network's neurons using a membrane-permeable [Ca²⁺] indicator (e.g., Fluo-4/AM) and record their activities. We then exploit the recorded temporal activity patterns as a means of intrinsic contrast to visualize individual neurons' dendritic morphology. The result is a high-contrast, multicolor visualization of the neuronal network. Taking the Xenopus olfactory bulb as an example we show the activities of the mitral/tufted cells of the olfactory bulb as well as their projections into the olfactory glomeruli. This method, yielding both functional and structural information of neuronal populations, will open up unprecedented possibilities for the investigation of neuronal networks.     

Dendritic sodium spikes endow neurons with inverse firing rate response to correlated synaptic activity

Many neurons possess dendrites enriched with sodium channels and are capable of generating action potentials. However, the role of dendritic sodium spikes remain unclear. Here, we study computational models of neurons to investigate the functional effects of dendritic spikes. In agreement with previous studies, we found that point neurons or neurons with passive dendrites increase their somatic firing rate in response to the correlation of synaptic bombardment for a wide range of input conditions, i.e. input firing rates, synaptic conductances, or refractory periods. However, neurons with active dendrites show the opposite behavior: for a wide range of conditions the firing rate decreases as a function of correlation. We found this property in three types of models of dendritic excitability: a Hodgkin-Huxley model of dendritic spikes, a model with integrate and fire dendrites, and a discrete-state dendritic model. We conclude that fast dendritic spikes confer much broader computational properties to neurons, sometimes opposite to that of point neurons.     

Neuron tracing and quantitative analyses of dendritic architecture reveal symmetrical three-way-junctions and phenotypes of git-1 in C. elegans

Complex dendritic trees are a distinctive feature of neurons. Alterations to dendritic morphology are associated with developmental, behavioral and neurodegenerative changes. The highly-arborized PVD neuron of C. elegans serves as a model to study dendritic patterning; however, quantitative, objective and automated analyses of PVD morphology are missing. Here, we present a method for neuronal feature extraction, based on deep-learning and fitting algorithms. The extracted neuronal architecture is represented by a database of structural elements for abstracted analysis. We obtain excellent automatic tracing of PVD trees and uncover that dendritic junctions are unevenly distributed. Surprisingly, these junctions are three-way-symmetrical on average, while dendritic processes are arranged orthogonally. We quantify the effect of mutation in git-1, a regulator of dendritic spine formation, on PVD morphology and discover a localized reduction in junctions. Our findings shed new light on PVD architecture, demonstrating the effectiveness of our objective analyses of dendritic morphology and suggest molecular control mechanisms.     

Synaptic integration in dendritic trees

Most neurons have elaborate dendritic trees that receive tens of thousands of synaptic inputs. Because postsynaptic responses to individual synaptic events are usually small and transient, the integration of many synaptic responses is needed to depolarize most neurons to action potential threshold. Over the past decade, advances in electrical and optical recording techniques have led to new insights into how synaptic responses propagate and interact within dendritic trees. In addition to their passive electrical and morphological properties, dendrites express active conductances that shape individual synaptic responses and influence synaptic integration locally within dendrites. Dendritic voltage-gated Na(+) and Ca(2+) channels support action potential backpropagation into the dendritic tree and local initiation of dendritic spikes, whereas K(+) conductances act to dampen dendritic excitability. While all dendrites investigated to date express active conductances, different neuronal types show specific patterns of dendritic channel expression leading to cell-specific differences in the way synaptic responses are integrated within dendritic trees. This review explores the way active and passive dendritic properties shape synaptic responses in the dendrites of central neurons, and emphasizes their role in synaptic integration.     

Firing pattern and synchronization property analysis in a network model of the olfactory bulb

In the olfactory system, both the temporal spike structure and spatial distribution of neuronal activity are important for processing odor information. In this paper, a biophysically-detailed, spiking neuronal model is used to simulate the activity of olfactory bulb. It is shown that by varying some key parameters such as maximal conductances of Ks and Nap the spike train of single neuron can exhibit various firing patterns. Synchronization in coupled neurons is also investigated as the coupling strength varying in the situation of two neurons and network. It is illustrated that the coupled neurons can exhibit different types of pattern when the coupling strength varies. These results may be instructive to understand information transmission in olfactory system.     

Correlated Conductance Parameters in Leech Heart Motor Neurons Contribute to Motor Pattern Formation

Neurons can have widely differing intrinsic membrane properties, in particular the density of specific conductances, but how these contribute to characteristic neuronal activity or pattern formation is not well understood. To explore the relationship between conductances, and in particular how they influence the activity of motor neurons in the well characterized leech heartbeat system, we developed a new multi-compartmental Hodgkin-Huxley style leech heart motor neuron model. To do so, we evolved a population of model instances, which differed in the density of specific conductances, capable of achieving specific output activity targets given an associated input pattern. We then examined the sensitivity of measures of output activity to conductances and how the model instances responded to hyperpolarizing current injections. We found that the strengths of many conductances, including those with differing dynamics, had strong partial correlations and that these relationships appeared to be linked by their influence on heart motor neuron activity. Conductances that had positive correlations opposed one another and had the opposite effects on activity metrics when perturbed whereas conductances that had negative correlations could compensate for one another and had similar effects on activity metrics.     

Information coding in vasopressin neurons—The role of asynchronous bistable burst firing

The task of the vasopressin system is homeostasis, a type of process which is fundamental to the brain's regulation of the body, exists in many different systems, and is vital to health and survival. Many illnesses are related to the dysfunction of homeostatic systems, including high blood pressure, obesity and diabetes. Beyond the vasopressin system's own importance, in regulating osmotic pressure, it presents an accessible model where we can learn how the features of homeostatic systems generally relate to their function, and potentially develop treatments. The vasopressin system is an important model system in neuroscience because it presents an accessible system in which to investigate the function and importance of, for example, dendritic release and burst firing, both of which are found in many systems of the brain. We have only recently begun to understand the contribution of dendritic release to neuronal function and information processing. Burst firing has most commonly been associated with rhythm generation; in this system it clearly plays a different role, still to be understood fully.     

Automated three-dimensional detection and shape classification of dendritic spines from fluorescence microscopy images

A fundamental challenge in understanding how dendritic spine morphology controls learning and memory has been quantifying three-dimensional (3D) spine shapes with sufficient precision to distinguish morphologic types, and sufficient throughput for robust statistical analysis. The necessity to analyze large volumetric data sets accurately, efficiently, and in true 3D has been a major bottleneck in deriving reliable relationships between altered neuronal function and changes in spine morphology. We introduce a novel system for automated detection, shape analysis and classification of dendritic spines from laser scanning microscopy (LSM) images that directly addresses these limitations. The system is more accurate, and at least an order of magnitude faster, than existing technologies. By operating fully in 3D the algorithm resolves spines that are undetectable with standard two-dimensional (2D) tools. Adaptive local thresholding, voxel clustering and Rayburst Sampling generate a profile of diameter estimates used to classify spines into morphologic types, while minimizing optical smear and quantization artifacts. The technique opens new horizons on the objective evaluation of spine changes with synaptic plasticity, normal development and aging, and with neurodegenerative disorders that impair cognitive function.     

Co-existent activity patterns in inhibitory neuronal networks with short-term synaptic depression

A network of two neurons mutually coupled through inhibitory synapses that display short-term synaptic depression is considered. We show that synaptic depression expands the number of possible activity patterns that the network can display and allows for co-existence of different patterns. Specifically, the network supports different types of n-m anti-phase firing patterns, where one neuron fires n spikes followed by the other neuron firing m spikes. When maximal synaptic conductances are identical, n-n anti-phase firing patterns are obtained and there are conductance intervals over which different pairs of these solutions co-exist. The multitude of n-m anti-phase patterns and their co-existence are not found when the synapses are non-depressing. Geometric singular perturbation methods for dynamical systems are applied to the original eight-dimensional model system to derive a set of one-dimensional conditions for the existence and co-existence of different anti-phase solutions. The generality and validity of these conditions are demonstrated through numerical simulations utilizing the Hodgkin-Huxley and Morris-Lecar neuronal models.     

Stochastic Ion Channel Gating in Dendritic Neurons: Morphology Dependence and Probabilistic Synaptic Activation of Dendritic Spikes

Neuronal activity is mediated through changes in the probability of stochastic transitions between open and closed states of ion channels. While differences in morphology define neuronal cell types and may underlie neurological disorders, very little is known about influences of stochastic ion channel gating in neurons with complex morphology. We introduce and validate new computational tools that enable efficient generation and simulation of models containing stochastic ion channels distributed across dendritic and axonal membranes. Comparison of five morphologically distinct neuronal cell types reveals that when all simulated neurons contain identical densities of stochastic ion channels, the amplitude of stochastic membrane potential fluctuations differs between cell types and depends on sub-cellular location. For typical neurons, the amplitude of membrane potential fluctuations depends on channel kinetics as well as open probability. Using a detailed model of a hippocampal CA1 pyramidal neuron, we show that when intrinsic ion channels gate stochastically, the probability of initiation of dendritic or somatic spikes by dendritic synaptic input varies continuously between zero and one, whereas when ion channels gate deterministically, the probability is either zero or one. At physiological firing rates, stochastic gating of dendritic ion channels almost completely accounts for probabilistic somatic and dendritic spikes generated by the fully stochastic model. These results suggest that the consequences of stochastic ion channel gating differ globally between neuronal cell-types and locally between neuronal compartments. Whereas dendritic neurons are often assumed to behave deterministically, our simulations suggest that a direct consequence of stochastic gating of intrinsic ion channels is that spike output may instead be a probabilistic function of patterns of synaptic input to dendrites.     

Exploring how extracellular electric field modulates neuron activity through dynamical analysis of a two-compartment neuron model

To investigate how extracellular electric field modulates neuron activity, a reduced two-compartment neuron model in the presence of electric field is introduced in this study. Depending on neuronal geometric and internal coupling parameters, the behaviors of the model have been studied extensively. The neuron model can exist in quiescent state or repetitive spiking state in response to electric field stimulus. Negative electric field mainly acts as inhibitory stimulus to the neuron, positive weak electric field could modulate spiking frequency and spike timing when the neuron is already active, and positive electric fields with sufficient intensity could directly trigger neuronal spiking in the absence of other stimulations. By bifurcation analysis, it is observed that there is saddle-node on invariant circle bifurcation, supercritical Hopf bifurcation and subcritical Hopf bifurcation appearing in the obtained two parameter bifurcation diagrams. The bifurcation structures and electric field thresholds for triggering neuron firing are determined by neuronal geometric and coupling parameters. The model predicts that the neurons with a nonsymmetric morphology between soma and dendrite, are more sensitive to electric field stimulus than those with the spherical structure. These findings suggest that neuronal geometric features play a crucial role in electric field effects on the polarization of neuronal compartments. Moreover, by determining the electric field threshold of our biophysical model, we could accurately distinguish between suprathreshold and subthreshold electric fields. Our study highlights the effects of extracellular electric field on neuronal activity from the biophysical modeling point of view. These insights into the dynamical mechanism of electric field may contribute to the investigation and development of electromagnetic therapies, and the model in our study could be further extended to a neuronal network in which the effects of electric fields on network activity may be investigated.     

Parameter Estimation Methods for Single Neuron Models

With the advancement in computer technology, it has become possible to fit complex models to neuronal data. In this work, we test how two methods can estimate parameters of simple neuron models (passive soma) to more complex ones (neuron with one dendritic cylinder and two active conductances). The first method uses classical voltage traces resulting from current pulses injection (time domain), while the second uses measures of the neuron's response to sinusoidal stimuli (frequency domain). Both methods estimate correctly the parameters in all cases studied. However, the time-domain method is slower and more prone to estimation errors in the cable parameters than the frequency-domain method. Because with noisy data the goodness of fit does not distinguish between different solutions, we suggest that running the estimation procedure a large number of times might help find a good solution and can provide information about the interactions between parameters. Also, because the formulation used for the model's response in the frequency domain is analytical, one can derive a local sensitivity analysis for each parameter. This analysis indicates how well a parameter is likely to be estimated and helps choose an optimal stimulation protocol. Finally, the tests suggest a strategy for fitting single-cell models using the two methods examined.