Response Properties and Synchronization of Rhythmically Firing Dendritic Neurons

Optogenetic Interrogation of Electrophysiological Dendritic Properties and Their Effect on Pacemaking Neurons from Acute Rodent Brain Slices

Understanding dendritic excitability is essential for a complete and precise characterization of neurons' input-output relationships. Theoretical and experimental work demonstrates that the electrotonic and nonlinear properties of dendrites can alter the amplitude (e.g., through amplification) and latency of synaptic inputs as viewed in the axosomatic region where spike timing is determined. The gold-standard technique to study dendritic excitability is using dual-patch recordings with a high-resistance electrode used to patch a piece of distal dendrite in addition to a somatic patch electrode. However, this approach is often impractical when distal dendrites are too fine to patch. Therefore, we developed a technique that utilizes the expression of Channelrhodopsin-2 (ChR2) to study dendritic excitability in acute brain slices through the combination of a somatic patch electrode and optogenetic activation. The protocol describes how to prepare acute slices from mice that express ChR2 in specific cell types, and how to use two modes of light stimulation: proximal (which activates the soma and proximal dendrites in a ~100 µm diameter surrounding the soma) with the use of a high-magnification objective and full-field stimulation through a low-magnification objective (which activates the entire somato-dendritic field of the neuron). We use this technique in conjunction with various stimulation protocols to estimate model-based spectral components of dendritic filtering and the impact of dendrites on phase response curves, peri-stimulus time histograms, and entrainment of pacemaking neurons. This technique provides a novel use of optogenetics to study intrinsic dendritic excitability through the use of standard patch-clamp slice physiology. Key features • A method for studying the effects of electrotonic and nonlinear dendritic properties on the sub- and suprathreshold responses of pacemaking neurons. • Combines somatic patch clamp or perforated patch recordings with optogenetic activation in acute brain slices to investigate dendritic linear transformation without patching the dendrite. • Oscillatory illumination at various frequencies estimates spectral properties of the dendrite using subthreshold voltage-clamp recordings and studies entrainment of pacemakers in current clamp recordings. • This protocol uses Poisson white noise illumination to estimate dendritic phase response curves and peri-stimulus time histograms.

How neuronal morphology impacts the synchronisation state of neuronal networks

The biophysical properties of neurons not only affect how information is processed within cells, they can also impact the dynamical states of the network. Specifically, the cellular dynamics of action-potential generation have shown relevance for setting the (de)synchronisation state of the network. The dynamics of tonically spiking neurons typically fall into one of three qualitatively distinct types that arise from distinct mathematical bifurcations of voltage dynamics at the onset of spiking. Accordingly, changes in ion channel composition or even external factors, like temperature, have been demonstrated to switch network behaviour via changes in the spike onset bifurcation and hence its associated dynamical type. A thus far less addressed modulator of neuronal dynamics is cellular morphology. Based on simplified and anatomically realistic mathematical neuron models, we show here that the extent of dendritic arborisation has an influence on the neuronal dynamical spiking type and therefore on the (de)synchronisation state of the network. Specifically, larger dendritic trees prime neuronal dynamics for in-phase-synchronised or splayed-out activity in weakly coupled networks, in contrast to cells with otherwise identical properties yet smaller dendrites. Our biophysical insights hold for generic multicompartmental classes of spiking neuron models (from ball-and-stick-type to anatomically reconstructed models) and establish a connection between neuronal morphology and the susceptibility of neural tissue to synchronisation in health and disease.

Non-uniform distribution of dendritic nonlinearities differentially engages thalamostriatal and corticostriatal inputs onto cholinergic interneurons

The tonic activity of striatal cholinergic interneurons (CINs) is modified differentially by their afferent inputs. Although their unitary synaptic currents are identical, in most CINs cortical inputs onto distal dendrites only weakly entrain them, whereas proximal thalamic inputs trigger abrupt pauses in discharge in response to salient external stimuli. To test whether the dendritic expression of the active conductances that drive autonomous discharge contribute to the CINs’ capacity to dissociate cortical from thalamic inputs, we used an optogenetics-based method to quantify dendritic excitability in mouse CINs. We found that the persistent sodium (NaP) current gave rise to dendritic boosting, and that the hyperpolarization-activated cyclic nucleotide-gated (HCN) current gave rise to a subhertz membrane resonance. This resonance may underlie our novel finding of an association between CIN pauses and internally-generated slow wave events in sleeping non-human primates. Moreover, our method indicated that dendritic NaP and HCN currents were preferentially expressed in proximal dendrites. We validated the non-uniform distribution of NaP currents: pharmacologically; with two-photon imaging of dendritic back-propagating action potentials; and by demonstrating boosting of thalamic, but not cortical, inputs by NaP currents. Thus, the localization of active dendritic conductances in CIN dendrites mirrors the spatial distribution of afferent terminals and may promote their differential responses to thalamic vs. cortical inputs.

How to correctly quantify neuronal phase-response curves from noisy recordings

At the level of individual neurons, various coding properties can be inferred from the input-output relationship of a cell. For small inputs, this relation is captured by the phase-response curve (PRC), which measures the effect of a small perturbation on the timing of the subsequent spike. Experimentally, however, an accurate experimental estimation of PRCs is challenging. Despite elaborate measurement efforts, experimental PRC estimates often cannot be related to those from modeling studies. In particular, experimental PRCs rarely resemble the characteristic theoretical PRC expected close to spike initiation, which is indicative of the underlying spike-onset bifurcation. Here, we show for conductance-based model neurons that the correspondence between theoretical and measured phase-response curve is lost when the stimuli used for the estimation are too large. In this case, the derived phase-response curve is distorted beyond recognition and takes on a generic shape that reflects the measurement protocol and masks the spike-onset bifurcation. We discuss how to identify appropriate stimulus strengths for perturbation and noise-stimulation methods, which permit to estimate PRCs that reliably reflect the spike-onset bifurcation - a task that is particularly difficult if a lower bound for the stimulus amplitude is dictated by prominent intrinsic neuronal noise.

Synchronization of oscillators via active media

In this paper, we study pairs of oscillators that are indirectly coupled via active (excitable) cells. We introduce a scalar phase model for coupled oscillators and excitable cells. We first show that one excitable and one oscillatory cell will exhibit phase locking at a variety of m:n patterns. We next introduce a second oscillatory cell and show that the only attractor is synchrony between the oscillators. We will also study the robustness to heterogeneity when the excitable cell fires or is quiescent. We next examine the dynamics when the oscillators are coupled via two excitable cells. In this case, the dynamics are very complicated with many forms of bistability and, in some cases, chaotic behavior. We also apply weak-coupling analysis to this case and explain some of the degeneracies observed in the bifurcation diagram. Further, we look at pairs of oscillators coupled via long chains of excitable cells and show that small differences in the frequency of the oscillators makes their locking more robust. Finally, we demonstrate many of the same phenomena seen in the phase model for a gap-junction coupled system of Morris-Lecar neurons.

Population dynamics and entrainment of basal ganglia pacemakers are shaped by their dendritic arbors

The theory of phase oscillators is an essential tool for understanding population dynamics of pacemaking neurons. GABAergic pacemakers in the substantia nigra pars reticulata (SNr), a main basal ganglia (BG) output nucleus, receive inputs from the direct and indirect pathways at distal and proximal regions of their dendritic arbors, respectively. We combine theory, optogenetic stimulation and electrophysiological experiments in acute brain slices to ask how dendritic properties impact the propensity of the various inputs, arriving at different locations along the dendrite, to recruit or entrain SNr pacemakers. By combining cable theory with sinusoidally-modulated optogenetic activation of either proximal somatodendritic regions or the entire somatodendritic arbor of SNr neurons, we construct an analytical model that accurately fits the empirically measured somatic current response to inputs arising from illuminating the soma and various portions of the dendritic field. We show that the extent of the dendritic tree that is illuminated generates measurable and systematic differences in the pacemaker’s phase response curve (PRC), causing a shift in its peak. Finally, we show that the divergent PRCs correctly predict differences in two major features of the collective dynamics of SNr neurons: the fidelity of population responses to sudden step-like changes in inputs; and the phase latency at which SNr neurons are entrained by rhythmic stimulation, which can occur in the BG under both physiological and pathophysiological conditions. Our novel method generates measurable and physiologically meaningful spatial effects, and provides the first empirical demonstration of how the collective responses of SNr pacemakers are determined by the transmission properties of their dendrites. SNr dendrites may serve to delay distal striatal inputs so that they impinge on the spike initiation zone simultaneously with pallidal and subthalamic inputs in order to guarantee a fair competition between the influence of the monosynaptic direct- and polysynaptic indirect pathways.

The substantia nigra pars reticulata (SNr) is a main output nucleus of the basal ganglia (BG), where inputs from the competing direct and indirect pathways converge onto the same neurons. Interestingly, these inputs are differentially distributed with direct and indirect pathway projections arriving at distal and proximal regions of the dendritic arbor, respectively. We employ a novel method combining theory with electrophysiological experiments and optogenetics to study the distinct effects of inputs arriving at different locations along the dendrite. Our approach represents a useful compromise between complexity and reduction in modelling. Our work addresses the question of high fidelity encoding of inputs by networks of neurons in the new context of pacemaking neurons, which are driven to fire by their intrinsic dynamics rather than by a network state. We provide the first empirical demonstration that dendritic delays can introduce latencies in the responses of a population of neurons that are commensurate with synaptic delays, suggesting a new role for SNr dendrites with implications for BG function.

The role of negative conductances in neuronal subthreshold properties and synaptic integration

Based on passive cable theory, an increase in membrane conductance produces a decrease in the membrane time constant and input resistance. Unlike the classical leak currents, voltage-dependent currents have a nonlinear behavior which can create regions of negative conductance, despite the increase in membrane conductance (permeability). This negative conductance opposes the effects of the passive membrane conductance on the membrane input resistance and time constant, increasing their values and thereby substantially affecting the amplitude and time course of postsynaptic potentials at the voltage range of the negative conductance. This paradoxical effect has been described for three types of voltage-dependent inward currents: persistent sodium currents, L- and T-type calcium currents and ligand-gated glutamatergic N-methyl-D-aspartate currents. In this review, we describe the impact of the creation of a negative conductance region by these currents on neuronal membrane properties and synaptic integration. We also discuss recent contributions of the quasi-active cable approximation, an extension of the passive cable theory that includes voltage-dependent currents, and its effects on neuronal subthreshold properties.

Neurons as oscillators

Regularly spiking neurons can be described as oscillators. In this article we review some of the insights gained from this conceptualization and their relevance for systems neuroscience. First, we explain how a regularly spiking neuron can be viewed as an oscillator and how the phase-response curve (PRC) describes the response of the neuron's spike times to small perturbations. We then discuss the meaning of the PRC for a single neuron's spiking behavior and review the PRCs measured from a variety of neurons in a range of spiking regimes. Next, we show how the PRC can be related to a number of common measures used to quantify neuronal firing, such as the spike-triggered average and the peristimulus histogram. We further show that the response of a neuron to correlated inputs depends on the shape of the PRC. We then explain how the PRC of single neurons can be used to predict neural network behavior. Given the PRC, conduction delays, and the waveform and time course of the synaptic potentials, it is possible to predict neural population behavior such as synchronization. The PRC also allows us to quantify the robustness of the synchronization to heterogeneity and noise. We finally ask how to combine the measured PRCs and the predictions based on PRC to further the understanding of systems neuroscience. As an example, we discuss how the change of the PRC by the neuromodulator acetylcholine could lead to a destabilization of cortical network dynamics. Although all of these studies are grounded in mathematical abstractions that do not strictly hold in biology, they provide good estimates for the emergence of the brain's network activity from the properties of individual neurons. The study of neurons as oscillators can provide testable hypotheses and mechanistic explanations for systems neuroscience.

Active subthreshold dendritic conductances shape the local field potential

Key points

The local field potential (LFP), the low‐frequency part of extracellular potentials recorded in neural tissue, is often used for probing neural circuit activity. Interpreting the LFP signal is difficult, however.

While the cortical LFP is thought mainly to reflect synaptic inputs onto pyramidal neurons, little is known about the role of the various subthreshold active conductances in shaping the LFP.

By means of biophysical modelling we obtain a comprehensive qualitative understanding of how the LFP generated by a single pyramidal neuron depends on the type and spatial distribution of active subthreshold currents.

For pyramidal neurons, the h‐type channels probably play a key role and can cause a distinct resonance in the LFP power spectrum.

Our results show that the LFP signal can give information about the active properties of neurons and imply that preferred frequencies in the LFP can result from those cellular properties instead of, for example, network dynamics.

Abstract

The main contribution to the local field potential (LFP) is thought to stem from synaptic input to neurons and the ensuing subthreshold dendritic processing. The role of active dendritic conductances in shaping the LFP has received little attention, even though such ion channels are known to affect the subthreshold neuron dynamics. Here we used a modelling approach to investigate the effects of subthreshold dendritic conductances on the LFP. Using a biophysically detailed, experimentally constrained model of a cortical pyramidal neuron, we identified conditions under which subthreshold active conductances are a major factor in shaping the LFP. We found that, in particular, the hyperpolarization‐activated inward current, I_h, can have a sizable effect and cause a resonance in the LFP power spectral density. To get a general, qualitative understanding of how any subthreshold active dendritic conductance and its cellular distribution can affect the LFP, we next performed a systematic study with a simplified model. We found that the effect on the LFP is most pronounced when (1) the synaptic drive to the cell is asymmetrically distributed (i.e. either basal or apical), (2) the active conductances are distributed non‐uniformly with the highest channel densities near the synaptic input and (3) when the LFP is measured at the opposite pole of the cell relative to the synaptic input. In summary, we show that subthreshold active conductances can be strongly reflected in LFP signals, opening up the possibility that the LFP can be used to characterize the properties and cellular distributions of active conductances.

The local field potential (LFP), the low‐frequency part of extracellular potentials recorded in neural tissue, is often used for probing neural circuit activity. Interpreting the LFP signal is difficult, however.

While the cortical LFP is thought mainly to reflect synaptic inputs onto pyramidal neurons, little is known about the role of the various subthreshold active conductances in shaping the LFP.

By means of biophysical modelling we obtain a comprehensive qualitative understanding of how the LFP generated by a single pyramidal neuron depends on the type and spatial distribution of active subthreshold currents.

For pyramidal neurons, the h‐type channels probably play a key role and can cause a distinct resonance in the LFP power spectrum.

Our results show that the LFP signal can give information about the active properties of neurons and imply that preferred frequencies in the LFP can result from those cellular properties instead of, for example, network dynamics.

Predicting the responses of repetitively firing neurons to current noise

We used phase resetting methods to predict firing patterns of rat subthalamic nucleus (STN) neurons when their rhythmic firing was densely perturbed by noise. We applied sequences of contiguous brief (0.5–2 ms) current pulses with amplitudes drawn from a Gaussian distribution (10–100 pA standard deviation) to autonomously firing STN neurons in slices. Current noise sequences increased the variability of spike times with little or no effect on the average firing rate. We measured the infinitesimal phase resetting curve (PRC) for each neuron using a noise-based method. A phase model consisting of only a firing rate and PRC was very accurate at predicting spike timing, accounting for more than 80% of spike time variance and reliably reproducing the spike-to-spike pattern of irregular firing. An approximation for the evolution of phase was used to predict the effect of firing rate and noise parameters on spike timing variability. It quantitatively predicted changes in variability of interspike intervals with variation in noise amplitude, pulse duration and firing rate over the normal range of STN spontaneous rates. When constant current was used to drive the cells to higher rates, the PRC was altered in size and shape and accurate predictions of the effects of noise relied on incorporating these changes into the prediction. Application of rate-neutral changes in conductance showed that changes in PRC shape arise from conductance changes known to accompany rate increases in STN neurons, rather than the rate increases themselves. Our results show that firing patterns of densely perturbed oscillators cannot readily be distinguished from those of neurons randomly excited to fire from the rest state. The spike timing of repetitively firing neurons may be quantitatively predicted from the input and their PRCs, even when they are so densely perturbed that they no longer fire rhythmically.

Most neurons receive thousands of synaptic inputs per second. Each of these may be individually weak but collectively they shape the temporal pattern of firing by the postsynaptic neuron. If the postsynaptic neuron fires repetitively, its synaptic inputs need not directly trigger action potentials, but may instead control the timing of action potentials that would occur anyway. The phase resetting curve encapsulates the influence of an input on the timing of the next action potential, depending on its time of arrival. We measured the phase resetting curves of neurons in the subthalamic nucleus and used them to accurately predict the timing of action potentials in a phase model subjected to complex input patterns. A simple approximation to the phase model accurately predicted the changes in firing pattern evoked by dense patterns of noise pulses varying in amplitude and pulse duration, and by changes in firing rate. We also showed that the phase resetting curve changes systematically with changes in total neuron conductance, and doing so predicts corresponding changes in firing pattern. Our results indicate that the phase model may accurately represent the temporal integration of complex patterns of input to repetitively firing neurons.

Spectral reconstruction of phase response curves reveals the synchronization properties of mouse globus pallidus neurons

The propensity of a neuron to synchronize is captured by its infinitesimal phase response curve (iPRC). Determining whether an iPRC is biphasic, meaning that small depolarizing perturbations can actually delay the next spike, if delivered at appropriate phases, is a daunting experimental task because negative lobes in the iPRC (unlike positive ones) tend to be small and may be occluded by the normal discharge variability of a neuron. To circumvent this problem, iPRCs are commonly derived from numerical models of neurons. Here, we propose a novel and natural method to estimate the iPRC by direct estimation of its spectral modes. First, we show analytically that the spectral modes of the iPRC of an arbitrary oscillator are readily measured by applying weak harmonic perturbations. Next, applying this methodology to biophysical neuronal models, we show that a low-dimensional spectral reconstruction is sufficient to capture the structure of the iPRC. This structure was preserved reasonably well even with added physiological scale jitter in the neuronal models. To validate the methodology empirically, we applied it first to a low-noise electronic oscillator with a known design and then to cortical pyramidal neurons, recorded in whole cell configuration, that are known to possess a monophasic iPRC. Finally, using the methodology in conjunction with perforated-patch recordings from pallidal neurons, we show, in contrast to recent modeling studies, that these neurons have biphasic somatic iPRCs. Biphasic iPRCs would cause lateral somatically targeted pallidal inhibition to desynchronize pallidal neurons, providing a plausible explanation for their lack of synchrony in vivo.

Neural dynamics and information representation in microcircuits of motor cortex

The brain has to analyze and respond to external events that can change rapidly from time to time, suggesting that information processing by the brain may be essentially dynamic rather than static. The dynamical features of neural computation are of significant importance in motor cortex that governs the process of movement generation and learning. In this paper, we discuss these features based primarily on our recent findings on neural dynamics and information coding in the microcircuit of rat motor cortex. In fact, cortical neurons show a variety of dynamical behavior from rhythmic activity in various frequency bands to highly irregular spike firing. Of particular interest are the similarity and dissimilarity of the neuronal response properties in different layers of motor cortex. By conducting electrophysiological recordings in slice preparation, we report the phase response curves (PRCs) of neurons in different cortical layers to demonstrate their layer-dependent synchronization properties. We then study how motor cortex recruits task-related neurons in different layers for voluntary arm movements by simultaneous juxtacellular and multiunit recordings from behaving rats. The results suggest an interesting difference in the spectrum of functional activity between the superficial and deep layers. Furthermore, the task-related activities recorded from various layers exhibited power law distributions of inter-spike intervals (ISIs), in contrast to a general belief that ISIs obey Poisson or Gamma distributions in cortical neurons. We present a theoretical argument that this power law of in vivo neurons may represent the maximization of the entropy of firing rate with limited energy consumption of spike generation. Though further studies are required to fully clarify the functional implications of this coding principle, it may shed new light on information representations by neurons and circuits in motor cortex.

Parameterized phase response curves for characterizing neuronal behaviors under transient conditions

Phase response curves (PRCs) are a simple model of how a neuron's spike time is affected by synaptic inputs. PRCs are useful in predicting how networks of neurons behave when connected. One challenge in estimating a neuron's PRCs experimentally is that many neurons do not have stationary firing rates. In this article we introduce a new method to estimate PRCs as a function of firing rate of the neuron. We call the resulting model a parameterized PRC (pPRC). Experimentally, we perturb the neuron applying a current with two parts: 1) a current held constant between spikes but changed at the onset of a spike, used to make the neuron fire at different rates, and 2) a pulse to emulate a synaptic input. A model of the applied constant current and the history is made to predict the interspike interval (ISI). A second model is then made to fit the modulation of the spike time from the expected ISI by the pulsatile stimulus. A polynomial with two independent variables, the stimulus phase and the expected ISI, is used to model the pPRC. The pPRC is validated in a computational model and applied to pyramidal neurons from the CA1 region of the hippocampal slices from rat. The pPRC can be used to model the effect of changing firing rates on network synchrony. It can also be used to characterize the effects of neuromodulators and genetic mutations (among other manipulations) on network synchrony. It can also easily be extended to account for more variables.

Phase response curves of subthalamic neurons measured with synaptic input and current injection

Infinitesimal phase response curves (iPRCs) provide a simple description of the response of repetitively firing neurons and may be used to predict responses to any pattern of synaptic input. Their simplicity makes them useful for understanding the dynamics of neurons when certain conditions are met. For example, the sizes of evoked phase shifts should scale linearly with stimulus strength, and the form of the iPRC should remain relatively constant as firing rate varies. We measured the PRCs of rat subthalamic neurons in brain slices using corticosubthalamic excitatory postsynaptic potentials (EPSPs; mediated by both AMPA- and NMDA-type receptors) and injected current pulses and used them to calculate the iPRC. These were relatively insensitive to both the size of the stimulus and the cell's firing rate, suggesting that the iPRC can predict the response of subthalamic nucleus cells to extrinsic inputs. However, the iPRC calculated using EPSPs differed from that obtained using current pulses. EPSPs (normalized for charge) were much more effective at altering the phase of subthalamic neurons than current pulses. The difference was not attributable to the extended time course of NMDA receptor-mediated currents, being unaffected by blockade of NMDA receptors. The iPRC provides a good description of subthalamic neurons' response to input, but iPRCs are best estimated using synaptic inputs rather than somatic current injection.

Robustness, variability, phase dependence, and longevity of individual synaptic input effects on spike timing during fluctuating synaptic backgrounds: a modeling study of globus pallidus neuron phase response properties

A neuron's phase response curve (PRC) shows how inputs arriving at different times during the spike cycle differentially affect the timing of subsequent spikes. Using a full morphological model of a globus pallidus (GP) neuron, we previously demonstrated that dendritic conductances shape the PRC in a spike frequency-dependent manner, suggesting different functional roles of perisomatic and distal dendritic synapses in the control of patterned network activity. In the present study we extend this analysis to examine the impact of physiologically realistic high conductance states on somatic and dendritic PRCs and the time course of spike train perturbations. First, we found that average somatic and dendritic PRCs preserved their shapes and spike frequency dependence when the model was driven by spatially-distributed, stochastic conductance inputs rather than tonic somatic current. However, responses to inputs during specific synaptic backgrounds often deviated substantially from the average PRC. Therefore, we analyzed the interactions of PRC stimuli with transient fluctuations in the synaptic background on a trial-by-trial basis. We found that the variability in responses to PRC stimuli and the incidence of stimulus-evoked added or skipped spikes were stimulus-phase-dependent and reflected the profile of the average PRC, suggesting commonality in the underlying mechanisms. Clear differences in the relation between the phase of input and variability of spike response between dendritic and somatic inputs indicate that these regions generally represent distinct dynamical subsystems of synaptic integration with respect to influencing the stability of spike time attractors generated by the overall synaptic conductance.

Engineering the synchronization of neuron action potentials using global time-delayed feedback stimulation

We experimentally demonstrate the use of continuous, time-delayed, feedback stimulation for controlling the synchronization of neuron action potentials. Phase-based models were experimentally constructed from a single synaptically isolated cultured hippocampal neuron. These models were used to determine the stimulation parameters necessary to produce the desired synchronization behavior in the action potentials of a pair of neurons coupled through a global time-delayed interaction. Measurements made using a dynamic clamp system confirm the generation of the synchronized states predicted by the experimentally constructed phase model. This model was then utilized to extrapolate the feedback stimulation parameters necessary to disrupt the action potential synchronization of a large population of globally interacting neurons.

Responses of a bursting pacemaker to excitation reveal spatial segregation between bursting and spiking mechanisms

Central pattern generators (CPGs) frequently include bursting neurons that serve as pacemakers for rhythm generation. Phase resetting curves (PRCs) can provide insight into mechanisms underlying phase locking in such circuits. PRCs were constructed for a pacemaker bursting complex in the pyloric circuit in the stomatogastric ganglion of the lobster and crab. This complex is comprised of the Anterior Burster (AB) neuron and two Pyloric Dilator (PD) neurons that are all electrically coupled. Artificial excitatory synaptic conductance pulses of different strengths and durations were injected into one of the AB or PD somata using the Dynamic Clamp. Previously, we characterized the inhibitory PRCs by assuming a single slow process that enabled synaptic inputs to trigger switches between an up state in which spiking occurs and a down state in which it does not. Excitation produced five different PRC shapes, which could not be explained with such a simple model. A separate dendritic compartment was required to separate the mechanism that generates the up and down phases of the bursting envelope (1) from synaptic inputs applied at the soma, (2) from axonal spike generation and (3) from a slow process with a slower time scale than burst generation. This study reveals that due to the nonlinear properties and compartmentalization of ionic channels, the response to excitation is more complex than inhibition.

Functional connectivity and integrative properties of globus pallidus neurons

The globus pallidus consists of the external (GPe) and the internal (GPi) segments. The GPe and GPi have different functional roles. The GPe is located centrally within multiple basal ganglia feedforward and feedback connections. The GPi is an output nucleus of the basal ganglia. A complex interplay between intrinsic pacemaking conductances and the balance of glutamatergic and GABAergic input largely determines the rate and pattern of firing of pallidal neurons. The initial part of this article introduces recent findings made in vivo that are related to the roles of glutamatergic and GABAergic inputs in the control of pallidal activity. The latter part describes the roles of intrinsic mechanisms of GPe neurons in the integration of the synaptic inputs. The presence of dendritic voltage-gated sodium channels may allow the initiation of dendritic spikes, giving distal inputs on the long and thin GPe dendrites an opportunity to strongly shape spiking activity. Basal ganglia disorders including Parkinson's disease, hemiballismus, and dystonias are accompanied by increased irregularity and synchronized bursts of pallidal activity. These changes may be in part due to changes in the GABA release in the GPe and GPi, but also involve intrinsic cellular changes in pallidal neurons.

Cellularly-driven differences in network synchronization propensity are differentially modulated by firing frequency

Spatiotemporal pattern formation in neuronal networks depends on the interplay between cellular and network synchronization properties. The neuronal phase response curve (PRC) is an experimentally obtainable measure that characterizes the cellular response to small perturbations, and can serve as an indicator of cellular propensity for synchronization. Two broad classes of PRCs have been identified for neurons: Type I, in which small excitatory perturbations induce only advances in firing, and Type II, in which small excitatory perturbations can induce both advances and delays in firing. Interestingly, neuronal PRCs are usually attenuated with increased spiking frequency, and Type II PRCs typically exhibit a greater attenuation of the phase delay region than of the phase advance region. We found that this phenomenon arises from an interplay between the time constants of active ionic currents and the interspike interval. As a result, excitatory networks consisting of neurons with Type I PRCs responded very differently to frequency modulation compared to excitatory networks composed of neurons with Type II PRCs. Specifically, increased frequency induced a sharp decrease in synchrony of networks of Type II neurons, while frequency increases only minimally affected synchrony in networks of Type I neurons. These results are demonstrated in networks in which both types of neurons were modeled generically with the Morris-Lecar model, as well as in networks consisting of Hodgkin-Huxley-based model cortical pyramidal cells in which simulated effects of acetylcholine changed PRC type. These results are robust to different network structures, synaptic strengths and modes of driving neuronal activity, and they indicate that Type I and Type II excitatory networks may display two distinct modes of processing information.

Synchronization of the firing of neurons in the brain is related to many cognitive functions, such as recognizing faces, discriminating odors, and coordinating movement. It is therefore important to understand what properties of neuronal networks promote synchrony of neural firing. One measure that is often used to determine the contribution of individual neurons to network synchrony is called the phase response curve (PRC). PRCs describe how the timing of neuronal firing changes depending on when input, such as a synaptic signal, is received by the neuron. A characteristic of PRCs that has previously not been well understood is that they change dramatically as the neuron's firing frequency is modulated. This effect carries potential significance, since cognitive functions are often associated with specific frequencies of network activity in the brain. We showed computationally that the frequency dependence of PRCs can be explained by the relative timing of ionic membrane currents with respect to the time between spike firings. Our simulations also showed that the frequency dependence of neuronal PRCs leads to frequency-dependent changes in network synchronization that can be different for different neuron types. These results further our understanding of how synchronization is generated in the brain to support various cognitive functions.

Effects of dendritic load on the firing frequency of oscillating neurons

We study the effects of passive dendritic properties on the dynamics of neuronal oscillators. We find that the addition of a passive dendrite can sometimes have counterintuitive effects on firing frequency. Specifically, the addition of a hyperpolarized passive dendritic load can either increase, decrease, or have negligible effects on firing frequency. We use the theory of weak coupling to derive phase equations for "ball-and-stick" model neurons and two-compartment model neurons. We then develop a framework for understanding how the addition of passive dendrites modulates the frequency of neuronal oscillators. We show that the average value of the neuronal oscillator's phase response curves measures the sensitivity of the neuron's firing rate to the dendritic load, including whether the addition of the dendrite causes an increase or decrease in firing frequency. We interpret this finding in terms of to the slope of the neuronal oscillator's frequency-applied current curve. We also show that equivalent results exist for constant and noisy point-source input to the dendrite. We note that the results are not specific to neurons but are applicable to any oscillator subject to a passive load.

Role of active dendritic conductances in subthreshold input integration

Dendrites of many types of neurons contain voltage-dependent conductances that are active at subthreshold membrane potentials. To understand the computations neurons perform it is key to understand the role of active dendrites in the subthreshold processing of synaptic inputs. We examine systematically how active dendritic conductances affect the time course of postsynaptic potentials propagating along dendrites, and how they affect the interaction between such signals. Voltage-dependent currents can be classified into two types that have qualitatively different effects on subthreshold input responses: regenerative dendritic currents boost and broaden EPSPs, while restorative currents attenuate and narrow EPSPs. Importantly, the effects of active dendritic currents on EPSP shape increase as the EPSP travels along the dendrite. The effectiveness of active currents in modulating the EPSP shape is determined by their activation time constant: the faster it is, the stronger the effect on EPSP amplitude, while the largest effects on EPSP width occur when it is comparable to the membrane time constant. We finally demonstrate that the two current types can differentially improve precision and robustness of neural computations: restorative currents enhance coincidence detection of dendritic inputs, whereas direction selectivity to sequences of dendritic inputs is enhanced by regenerative dendritic currents.

Phase-response curves and synchronized neural networks

We review the principal assumptions underlying the application of phase-response curves (PRCs) to synchronization in neuronal networks. The PRC measures how much a given synaptic input perturbs spike timing in a neural oscillator. Among other applications, PRCs make explicit predictions about whether a given network of interconnected neurons will synchronize, as is often observed in cortical structures. Regarding the assumptions of the PRC theory, we conclude: (i) The assumption of noise-tolerant cellular oscillations at or near the network frequency holds in some but not all cases. (ii) Reduced models for PRC-based analysis can be formally related to more realistic models. (iii) Spike-rate adaptation limits PRC-based analysis but does not invalidate it. (iv) The dependence of PRCs on synaptic location emphasizes the importance of improving methods of synaptic stimulation. (v) New methods can distinguish between oscillations that derive from mutual connections and those arising from common drive. (vi) It is helpful to assume linear summation of effects of synaptic inputs; experiments with trains of inputs call this assumption into question. (vii) Relatively subtle changes in network structure can invalidate PRC-based predictions. (viii) Heterogeneity in the preferred frequencies of component neurons does not invalidate PRC analysis, but can annihilate synchronous activity.

Neurophysiological and computational principles of cortical rhythms in cognition

Synchronous rhythms represent a core mechanism for sculpting temporal coordination of neural activity in the brain-wide network. This review focuses on oscillations in the cerebral cortex that occur during cognition, in alert behaving conditions. Over the last two decades, experimental and modeling work has made great strides in elucidating the detailed cellular and circuit basis of these rhythms, particularly gamma and theta rhythms. The underlying physiological mechanisms are diverse (ranging from resonance and pacemaker properties of single cells to multiple scenarios for population synchronization and wave propagation), but also exhibit unifying principles. A major conceptual advance was the realization that synaptic inhibition plays a fundamental role in rhythmogenesis, either in an interneuronal network or in a reciprocal excitatory-inhibitory loop. Computational functions of synchronous oscillations in cognition are still a matter of debate among systems neuroscientists, in part because the notion of regular oscillation seems to contradict the common observation that spiking discharges of individual neurons in the cortex are highly stochastic and far from being clocklike. However, recent findings have led to a framework that goes beyond the conventional theory of coupled oscillators and reconciles the apparent dichotomy between irregular single neuron activity and field potential oscillations. From this perspective, a plethora of studies will be reviewed on the involvement of long-distance neuronal coherence in cognitive functions such as multisensory integration, working memory, and selective attention. Finally, implications of abnormal neural synchronization are discussed as they relate to mental disorders like schizophrenia and autism.

Phase-resetting curve determines how BK currents affect neuronal firing

BK channels are large conductance potassium channels gated by calcium and voltage. Paradoxically, blocking these channels has been shown experimentally to increase or decrease the firing rate of neurons, depending on the neural subtype and brain region. The mechanism for how this current can alter the firing rates of different neurons remains poorly understood. Using phase-resetting curve (PRC) theory, we determine when BK channels increase or decrease the firing rates in neural models. The addition of BK currents always decreases the firing rate when the PRC has only a positive region. When the PRC has a negative region (type II), BK currents can increase the firing rate. The influence of BK channels on firing rate in the presence of other conductances, such as I(m) and I(h), as well as with different amplitudes of depolarizing input, were also investigated. These results provide a formal explanation for the apparently contradictory effects of BK channel antagonists on firing rates.

Democracy-independence trade-off in oscillating dendrites and its implications for grid cells

Dendritic democracy and independence have been characterized for near-instantaneous processing of synaptic inputs. However, a wide class of neuronal computations requires input integration on long timescales. As a paradigmatic example, entorhinal grid fields have been thought to be generated by the democratic summation of independent dendritic oscillations performing direction-selective path integration. We analyzed how multiple dendritic oscillators embedded in the same neuron integrate inputs separately and determine somatic membrane voltage jointly. We found that the interaction of dendritic oscillations leads to phase locking, which sets an upper limit on the timescale for independent input integration. Factors that increase this timescale also decrease the influence that the dendritic oscillations exert on somatic voltage. In entorhinal stellate cells, interdendritic coupling dominates and causes these cells to act as single oscillators. Our results suggest a fundamental trade-off between local and global processing in dendritic trees integrating ongoing signals.

Bayesian estimation of phase response curves

Phase response curve (PRC) of an oscillatory neuron describes the response of the neuron to external perturbation. The PRC is useful to predict synchronized dynamics of neurons; hence, its measurement from experimental data attracts increasing interest in neural science. This paper introduces a Bayesian method for estimating PRCs from data, which allows for the correlation of errors in explanatory and response variables of the PRC. The method is implemented with a replica exchange Monte Carlo technique; this avoids local minima and enables efficient calculation of posterior averages. A test with artificial data generated by the noisy Morris-Lecar equation shows that the proposed method outperforms conventional regression that ignores errors in the explanatory variable. Experimental data from the pyramidal cells in the rat motor cortex is also analyzed with the method; a case is found where the result with the proposed method is considerably different from that obtained by conventional regression.

Phase response curve analysis of a full morphological globus pallidus neuron model reveals distinct perisomatic and dendritic modes of synaptic integration

Synchronization of globus pallidus (GP) neurons and cortically entrained oscillations between GP and other basal ganglia nuclei are key features of the pathophysiology of Parkinson's disease. Phase response curves (PRCs), which tabulate the effects of phasic inputs within a neuron's spike cycle on output spike timing, are efficient tools for predicting the emergence of synchronization in neuronal networks and entrainment to periodic input. In this study we apply physiologically realistic synaptic conductance inputs to a full morphological GP neuron model to determine the phase response properties of the soma and different regions of the dendritic tree. We find that perisomatic excitatory inputs delivered throughout the interspike interval advance the phase of the spontaneous spike cycle yielding a type I PRC. In contrast, we demonstrate that distal dendritic excitatory inputs can either delay or advance the next spike depending on whether they occur early or late in the spike cycle. We find this latter pattern of responses, summarized by a biphasic (type II) PRC, was a consequence of dendritic activation of the small conductance calcium-activated potassium current, SK. We also evaluate the spike-frequency dependence of somatic and dendritic PRC shapes, and we demonstrate the robustness of our results to variations of conductance densities, distributions, and kinetic parameters. We conclude that the distal dendrite of GP neurons embodies a distinct dynamical subsystem that could promote synchronization of pallidal networks to excitatory inputs. These results highlight the need to consider different effects of perisomatic and dendritic inputs in the control of network behavior.

A bushy cell network in the rat ventral cochlear nucleus

Geometry of the dendritic tree and synaptic organization of afferent inputs are essential factors in determining how synaptic input is integrated by neurons. This information remains elusive for one of the first brainstem neurons involved in processing of the primary auditory signal from the ear, the bushy cells (BCs) of the ventral cochlear nucleus (VCN). Here, we labeled the BC dendritic trees with retrograde tracing techniques to analyze their geometry and synaptic organization after immunofluorescence for excitatory and inhibitory synaptic markers, electron microscopy, morphometry, double tract-tracing methods, and 3D reconstructions. Our study revealed that BC dendrites provide space for a large number of compartmentalized excitatory and inhibitory synaptic interactions. The dendritic inputs on BCs are of cochlear and noncochlear origin, and their proportion and distribution are dependent on the branching pattern and orientation of the dendritic tree in the VCN. Three-dimensional reconstructions showed that BC dendrites branch and cluster with those of other BCs in the core of the VCN. Within the cluster, incoming synaptic inputs establish divergent multiple-contact synapses (dyads and triads) between BCs. Furthermore, neuron-neuron connections including puncta adherentia, sarcoplasmic junctions, and gap junctions are common between BCs, which suggests that these neurons are electrically coupled. Overall, our study demonstrates the existence of a BC network in the rat VCN. This network may establish the neuroanatomical basis for acoustic information processing by individual BCs as well as for enhanced synchronization of the output signal of the VCN.

The role of ongoing dendritic oscillations in single-neuron dynamics

The dendritic tree contributes significantly to the elementary computations a neuron performs while converting its synaptic inputs into action potential output. Traditionally, these computations have been characterized as both temporally and spatially localized. Under this localist account, neurons compute near-instantaneous mappings from their current input to their current output, brought about by somatic summation of dendritic contributions that are generated in functionally segregated compartments. However, recent evidence about the presence of oscillations in dendrites suggests a qualitatively different mode of operation: the instantaneous phase of such oscillations can depend on a long history of inputs, and under appropriate conditions, even dendritic oscillators that are remote may interact through synchronization. Here, we develop a mathematical framework to analyze the interactions of local dendritic oscillations and the way these interactions influence single cell computations. Combining weakly coupled oscillator methods with cable theoretic arguments, we derive phase-locking states for multiple oscillating dendritic compartments. We characterize how the phase-locking properties depend on key parameters of the oscillating dendrite: the electrotonic properties of the (active) dendritic segment, and the intrinsic properties of the dendritic oscillators. As a direct consequence, we show how input to the dendrites can modulate phase-locking behavior and hence global dendritic coherence. In turn, dendritic coherence is able to gate the integration and propagation of synaptic signals to the soma, ultimately leading to an effective control of somatic spike generation. Our results suggest that dendritic oscillations enable the dendritic tree to operate on more global temporal and spatial scales than previously thought; notably that local dendritic activity may be a mechanism for generating on-going whole-cell voltage oscillations.

A central issue in biology is how local processes yield global consequences. This is especially relevant for neurons since these spatially extended cells process local synaptic inputs to generate global action potential output. The dendritic tree of a neuron, which receives most of the inputs, expresses ion channels that can generate nonlinear dynamics. A prominent phenomenon resulting from such ion channels are voltage oscillations. The distribution of the active membrane channels throughout the cell is often highly non-uniform. This can turn the dendritic tree into a network of sparsely spaced local oscillators. Here we analyze whether local dendritic oscillators can produce cell-wide voltage oscillations. Our mathematical theory shows that indeed even when the dendritic oscillators are weakly coupled, they lock their phases and give global oscillations. We show how the biophysical properties of the dendrites determine the global locking and how it can be controlled by synaptic inputs. As a consequence of global locking, even individual synaptic inputs can affect the timing of action potentials. In fact, dendrites locking in synchrony can lead to sustained firing of the cell. We show that dendritic trees can be bistable, with dendrites locking in either synchrony or asynchrony, which may provide a novel mechanism for single cell-based memory.

Predicting synchronous and asynchronous network groupings of hippocampal interneurons coupled with dendritic gap junctions

Direct electrical communication between central nervous system (CNS) neurons including those in the hippocampus is well-established. This form of communication is mediated by gap junctions and it is known that this coupling is important for brain rhythms such as gamma (20-80 Hz) which occur during active behavioural states. It is also known that gap junctions are present at several locations along the dendrites of hippocampal interneurons including parvalbumin-positive basket cell types. Weakly coupled oscillator theory, which uses phase response curves (PRCs), has been used to understand and predict the dynamics of electrically coupled networks. Here we use compartmental models of hippocampal basket cells with different levels of basal and apical spike attenuation together with the theory to show that network output can be broken down into three groupings: synchronous, asynchronous and antiphase-like patterns. Moreover, quantified PRCs can be used as a rule of thumb to determine the occurrence of a particular grouping under weak coupling conditions, which in turn implies that spike delays are critical factors in determining network output. In moving beyond weak coupling to encompass the full physiological regime of coupling strengths with network simulations, we note that it is important to be able to differentiate between these different groupings as it affects how the network responds with modulation. Specifically, an asynchronous grouping provides more dynamic richness as a larger range of phase-locked states can be expressed with strength changes. From a functional viewpoint it may be that modulation of electrically coupled networks are key to controlling cell assemblies that contribute to information coding brain substrates.

Spectral analysis of time series of events: effect of respiration on heart rate in neonates

Certain types of biomedical processes such as the heart rate generator can be considered as signals that are sampled by the occurring events, i.e. QRS complexes. This sampling property generates problems for the evaluation of spectral parameters of such signals. First, the irregular occurrence of heart beats creates an unevenly sampled data set which must either be pre-processed (e.g. by using trace binning or interpolation) prior to spectral analysis, or analyzed with specialized methods (e.g. Lomb's algorithm). Second, the average occurrence of events determines the Nyquist limit for the sampled time series. Here we evaluate different types of spectral analysis of recordings of neonatal heart rate. Coupling between respiration and heart rate and the detection of heart rate itself are emphasized. We examine both standard and data adaptive frequency bands of heart rate signals generated by models of coupled oscillators and recorded data sets from neonates. We find that an important spectral artifact occurs due to a mirror effect around the Nyquist limit of half the average heart rate. Further we conclude that the presence of respiratory coupling can only be detected under low noise conditions and if a data-adaptive respiratory band is used.

Modeling thalamocortical cell: impact of ca channel distribution and cell geometry on firing pattern

The influence of calcium channel distribution and geometry of the thalamocortical cell upon its tonic firing and the low threshold spike (LTS) generation was studied in a 3-compartment model, which represents soma, proximal and distal dendrites as well as in multi-compartment model using the morphology of a real reconstructed neuron. Using an uniform distribution of Ca²⁺ channels, we determined the minimal number of low threshold voltage-activated calcium channels and their permeability required for the onset of LTS in response to a hyperpolarizing current pulse. In the 3-compartment model, we found that the channel distribution influences the firing pattern only in the range of 3% below the threshold value of total T-channel density. In the multi-compartmental model, the LTS could be generated by only 64% of unequally distributed T-channels compared to the minimal number of equally distributed T-channels. For a given channel density and injected current, the tonic firing frequency was found to be inversely proportional to the size of the cell. However, when the Ca²⁺ channel density was elevated in soma or proximal dendrites, then the amplitude of LTS response and burst spike frequencies were determined by the ratio of total to threshold number of T-channels in the cell for a specific geometry.

The effects of cholinergic neuromodulation on neuronal phase-response curves of modeled cortical neurons

The response of an oscillator to perturbations is described by its phase-response curve (PRC), which is related to the type of bifurcation leading from rest to tonic spiking. In a recent experimental study, we have shown that the type of PRC in cortical pyramidal neurons can be switched by cholinergic neuromodulation from type II (biphasic) to type I (monophasic). We explored how intrinsic mechanisms affected by acetylcholine influence the PRC using three different types of neuronal models: a theta neuron, single-compartment neurons and a multi-compartment neuron. In all of these models a decrease in the amount of a spike-frequency adaptation current was a necessary and sufficient condition for the shape of the PRC to change from biphasic (type II) to purely positive (type I).

Interplay between a phase response curve and spike-timing-dependent plasticity leading to wireless clustering

A phase response curve (PRC) characterizes the signal transduction between oscillators such as neurons on a fixed network in a minimal manner, while spike-timing-dependent plasiticity (STDP) characterizes the way of rewiring networks in an activity-dependent manner. This paper demonstrates that these two key properties both related to the interaction times of oscillators work synergetically to carve functionally useful circuits. STDP working on neurons that prefer asynchrony converts the initial asynchronous firing to clustered firing with synchrony within a cluster. They get synchronized within a cluster despite their preference to asynchrony because STDP selectively disrupts intracluster connections, which we call wireless clustering. Our PRC analysis reveals a triad mechanism: the network structure affects how the PRC is read out to determine the synchrony tendency, the synchrony tendency affects how the STDP works, and STDP affects the network structure, closing the loop.

Class-II neurons display a higher degree of stochastic synchronization than class-I neurons

We describe the relationship between the shape of the phase-resetting curve (PRC) and the degree of stochastic synchronization observed between a pair of uncoupled general oscillators receiving partially correlated Poisson inputs in addition to inputs from independent sources. We use perturbation methods to derive an expression relating the shape of the PRC to the probability density function (PDF) of the phase difference between the oscillators. We compute various measures of the degree of synchrony and cross correlation from the PDF's and use the same to compare and contrast differently shaped PRCs, with respect to their ability to undergo stochastic synchronization. Since the shape of the PRC depends on underlying dynamical details of the oscillator system, we utilize the results obtained from the analysis of general oscillator systems to study specific models of neuronal oscillators. It is shown that the degree of stochastic synchronization is controlled both by the firing rate of the neuron and the membership of the PRC (type I or type II). It is also shown that the circular variance for the integrate and fire neuron and the generalized order parameter for a hippocampal interneuron model have a nonlinear relationship to the input correlation.

Layer and frequency dependencies of phase response properties of pyramidal neurons in rat motor cortex

It is postulated that synchronous firing of cortical neurons plays an active role in cognitive functions of the brain. An important issue is whether pyramidal neurons in different cortical layers exhibit similar tendencies to synchronise. To address this issue, we performed intracellular and whole-cell recordings of regular-spiking pyramidal neurons in slice preparations of the rat motor cortex (18-45 days old) and analysed the phase response curves of these pyramidal neurons in layers 2/3 and 5. The phase response curve represents how an external stimulus affects the timing of spikes immediately after the stimulus in repetitively firing neurons. The phase response curve can be classified into two categories, type 1 (the spike is always advanced) and type 2 (the spike is advanced or delayed depending on the stimulus phase), and are important determinants of whether or not rhythmic synchronization of neuron pairs occurs. We found that pyramidal neurons in layer 2/3 tend to display type-2 phase response curves whereas those in layer 5 tend to exhibit type-1 phase response curves. The differences were prominent particularly in the gamma-frequency range (20-45 Hz). Our results imply that the layer-2/3 pyramidal neurons, when coupled mutually through fast excitatory synapses, may exhibit a much stronger tendency for rhythmic synchronization than layer-5 neurons in the gamma-frequency range.

Synchronous and asynchronous bursting states: role of intrinsic neural dynamics

Brain signals such as local field potentials often display gamma-band oscillations (30-70 Hz) in a variety of cognitive tasks. These oscillatory activities possibly reflect synchronization of cell assemblies that are engaged in a cognitive function. A type of pyramidal neurons, i.e., chattering neurons, show fast rhythmic bursting (FRB) in the gamma frequency range, and may play an active role in generating the gamma-band oscillations in the cerebral cortex. Our previous phase response analyses have revealed that the synchronization between the coupled bursting neurons significantly depends on the bursting mode that is defined as the number of spikes in each burst. Namely, a network of neurons bursting through a Ca(2+)-dependent mechanism exhibited sharp transitions between synchronous and asynchronous firing states when the neurons exchanged the bursting mode between singlet, doublet and so on. However, whether a broad class of bursting neuron models commonly show such a network behavior remains unclear. Here, we analyze the mechanism underlying this network behavior using a mathematically tractable neuron model. Then we extend our results to a multi-compartment version of the NaP current-based neuron model and prove a similar tight relationship between the bursting mode changes and the network state changes in this model. Thus, the synchronization behavior couples tightly to the bursting mode in a wide class of networks of bursting neurons.