Contribution of near-threshold currents to intrinsic oscillatory activity in rat medial entorhinal cortex layer II stellate cells

Tau-mediated synaptic dysfunction is coupled with HCN channelopathy

INTRODUCTION

In tauopathies, altered tau processing correlates with impairments in synaptic density and function. Changes in hyperpolarization-activated cyclic nucleotide-gated (HCN) channels contribute to disease-associated abnormalities in multiple neurodegenerative diseases.

METHODS

To investigate the link between tau and HCN channels, we performed histological, biochemical, ultrastructural, and functional analyses of hippocampal tissues from Alzheimer's disease (AD), age-matched controls, Tau35 mice, and/or Tau35 primary hippocampal neurons.

RESULTS

Expression of specific HCN channels is elevated in post mortem AD hippocampus. Tau35 mice develop progressive abnormalities including increased phosphorylated tau, enhanced HCN channel expression, decreased dendritic branching, reduced synapse density, and vesicle clustering defects. Tau35 primary neurons show increased HCN channel expression enhanced hyperpolarization-induced membrane voltage "sag" and changes in the frequency and kinetics of spontaneous excitatory postsynaptic currents.

DISCUSSION

Our findings are consistent with a model in which pathological changes in tauopathies impact HCN channels to drive network-wide structural and functional synaptic deficits.

HIGHLIGHTS

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are functionally linked to the development of tauopathy. Expression of specific HCN channels is elevated in the hippocampus in Alzheimer's disease and the Tau35 mouse model of tauopathy. Increased expression of HCN channels in Tau35 mice is accompanied by hyperpolarization-induced membrane voltage "sag" demonstrating a detrimental effect of tau abnormalities on HCN channel function. Tau35 expression alters synaptic organization, causing a loosened vesicle clustering phenotype in Tau35 mice.

Calcium-Sensitive Subthreshold Oscillations and Electrical Coupling in Principal Cells of Mouse Dorsal Cochlear Nucleus

In higher sensory brain regions, slow oscillations (0.5-5 Hz) associated with quiet wakefulness and attention modulate multisensory integration, predictive coding, and perception. Although often assumed to originate via thalamocortical mechanisms, the extent to which subcortical sensory pathways are independently capable of slow oscillatory activity is unclear. We find that in the first station for auditory processing, the cochlear nucleus, fusiform cells from juvenile mice (of either sex) generate robust 1-2 Hz oscillations in membrane potential and exhibit electrical resonance. Such oscillations were absent prior to the onset of hearing, intrinsically generated by hyperpolarization-activated cyclic nucleotide-gated (HCN) and persistent Na+ conductances (NaP) interacting with passive membrane properties, and reflected the intrinsic resonance properties of fusiform cells. Cx36-containing gap junctions facilitated oscillation strength and promoted pairwise synchrony of oscillations between neighboring neurons. The strength of oscillations were strikingly sensitive to external Ca2+, disappearing at concentrations >1.7 mM, due in part to the shunting effect of small-conductance calcium-activated potassium (SK) channels. This effect explains their apparent absence in previous in vitro studies of cochlear nucleus which routinely employed high-Ca2+ extracellular solution. In contrast, oscillations were amplified in reduced Ca2+ solutions, due to relief of suppression by Ca2+ of Na+ channel gating. Our results thus reveal mechanisms for synchronous oscillatory activity in auditory brainstem, suggesting that slow oscillations, and by extension their perceptual effects, may originate at the earliest stages of sensory processing.

Entorhinal cortical delta oscillations drive memory consolidation

Long-term memories are formed by creating stable memory representations via memory consolidation, which mainly occurs during sleep following the encoding of labile memories in the hippocampus during waking. The entorhinal cortex (EC) has intricate connections with the hippocampus, but its role in memory consolidation is largely unknown. Using cell-type- and input-specific in vivo neural activity recordings, here we show that the temporoammonic pathway neurons in the EC, which directly innervate the output area of the hippocampus, exhibit potent oscillatory activities during anesthesia and sleep. Using in vivo individual and populational neuronal activity recordings, we demonstrate that a subpopulation of the temporoammonic pathway neurons, which we termed sleep cells, generate delta oscillations via hyperpolarization-activated cyclic-nucleotide-gated channels during sleep. The blockade of these oscillations significantly impaired the consolidation of hippocampus-dependent memory. Together, our findings uncover a key driver of delta oscillations and memory consolidation that are found in the EC.

Heterogeneous stochastic bifurcations explain intrinsic oscillatory patterns in entorhinal cortical stellate cells

Stellate cells (SC) in the medial entorhinal cortex manifest intrinsic membrane potential oscillatory patterns. Although different theoretical frameworks have been proposed to explain these patterns, a robust unifying framework that jointly accounts for intrinsic heterogeneities and stochasticity is missing. Here, we first performed in vitro patch-clamp electrophysiological recordings from rat SCs and found pronounced cell-to-cell variability in their characteristic physiological properties, including peri-threshold oscillatory patterns. We demonstrate that noise introduced into two independent populations (endowed with deterministic or stochastic ion-channel gating kinetics) of heterogeneous biophysical models yielded activity patterns that were qualitatively similar to electrophysiological peri-threshold oscillatory activity in SCs. We developed spectrogram-based quantitative metrics for the identification of valid oscillations and confirmed that these metrics reliably captured the variable-amplitude and arhythmic oscillatory patterns observed in electrophysiological recordings. Using these quantitative metrics, we validated activity patterns from both heterogeneous populations of SC models, with each model assessed with multiple trials of different levels of noise at distinct membrane depolarizations. Our analyses unveiled the manifestation of stochastic resonance (detection of the highest number of valid oscillatory traces at an optimal level of noise) in both heterogeneous populations of SC models. Finally, we show that a generalized network motif comprised of a slow negative feedback loop amplified by a fast positive feedback loop manifested stochastic bifurcations and stochastic resonance in the emergence of oscillations. Together, through a unique convergence of the degeneracy and stochastic resonance frameworks, our unifying framework centered on heterogeneous stochastic bifurcations argues for state-dependent emergence of SC oscillations.

Cell-Type Specific Inhibition Controls the High-Frequency Oscillations in the Medial Entorhinal Cortex

The medial entorhinal cortex (mEC) plays a critical role for spatial navigation and memory. While many studies have investigated the principal neurons within the entorhinal cortex, much less is known about the inhibitory circuitries within this structure. Here, we describe for the first time in the mEC a subset of parvalbumin-positive (PV+) interneurons (INs)-stuttering cells (STUT)-with morphological, intrinsic electrophysiological, and synaptic properties distinct from fast-spiking PV+ INs. In contrast to the fast-spiking PV+ INs, the axon of the STUT INs also terminated in layer 3 and showed subthreshold membrane oscillations at gamma frequencies. Whereas the synaptic output of the STUT INs was only weakly reduced by a μ-opioid agonist, their inhibitory inputs were strongly suppressed. Given these properties, STUT are ideally suited to entrain gamma activity in the pyramidal cell population of the mEC. We propose that activation of the μ-opioid receptors decreases the GABA release from the PV+ INs onto the STUT, resulting in disinhibition of the STUT cell population and the consequent increase in network gamma power. We therefore suggest that the opioid system plays a critical role, mediated by STUT INs, in the neural signaling and oscillatory network activity within the mEC.

The voltage and spiking responses of subthreshold resonant neurons to structured and fluctuating inputs: persistence and loss of resonance and variability

We systematically investigate the response of neurons to oscillatory currents and synaptic-like inputs and we extend our investigation to non-structured synaptic-like spiking inputs with more realistic distributions of presynaptic spike times. We use two types of chirp-like inputs consisting of (i) a sequence of cycles with discretely increasing frequencies over time, and (ii) a sequence having the same cycles arranged in an arbitrary order. We develop and use a number of frequency-dependent voltage response metrics to capture the different aspects of the voltage response, including the standard impedance (Z) and the peak-to-trough amplitude envelope ([Formula: see text]) profiles. We show that Z-resonant cells (cells that exhibit subthreshold resonance in response to sinusoidal inputs) also show [Formula: see text]-resonance in response to sinusoidal inputs, but generally do not (or do it very mildly) in response to square-wave and synaptic-like inputs. In the latter cases the resonant response using Z is not predictive of the preferred frequencies at which the neurons spike when the input amplitude is increased above subthreshold levels. We also show that responses to conductance-based synaptic-like inputs are attenuated as compared to the response to current-based synaptic-like inputs, thus providing an explanation to previous experimental results. These response patterns were strongly dependent on the intrinsic properties of the participating neurons, in particular whether the unperturbed Z-resonant cells had a stable node or a focus. In addition, we show that variability emerges in response to chirp-like inputs with arbitrarily ordered patterns where all signals (trials) in a given protocol have the same frequency content and the only source of uncertainty is the subset of all possible permutations of cycles chosen for a given protocol. This variability is the result of the multiple different ways in which the autonomous transient dynamics is activated across cycles in each signal (different cycle orderings) and across trials. We extend our results to include high-rate Poisson distributed current- and conductance-based synaptic inputs and compare them with similar results using additive Gaussian white noise. We show that the responses to both Poisson-distributed synaptic inputs are attenuated with respect to the responses to Gaussian white noise. For cells that exhibit oscillatory responses to Gaussian white noise (band-pass filters), the response to conductance-based synaptic inputs are low-pass filters, while the response to current-based synaptic inputs may remain band-pass filters, consistent with experimental findings. Our results shed light on the mechanisms of communication of oscillatory activity among neurons in a network via subthreshold oscillations and resonance and the generation of network resonance.

Membrane Resonance in Pyramidal and GABAergic Neurons of the Mouse Perirhinal Cortex

The perirhinal cortex (PRC) is a polymodal associative region of the temporal lobe that works as a gateway between cortical areas and hippocampus. In recent years, an increasing interest arose in the role played by the PRC in learning and memory processes, such as object recognition memory, in contrast with certain forms of hippocampus-dependent spatial and episodic memory. The integrative properties of the PRC should provide all necessary resources to select and enhance the information to be propagated to and from the hippocampus. Among these properties, we explore in this paper the ability of the PRC neurons to amplify the output voltage to current input at selected frequencies, known as membrane resonance. Within cerebral circuits the resonance of a neuron operates as a filter toward inputs signals at certain frequencies to coordinate network activity in the brain by affecting the rate of neuronal firing and the precision of spike timing. Furthermore, the ability of the PRC neurons to resonate could have a fundamental role in generating subthreshold oscillations and in the selection of cortical inputs directed to the hippocampus. Here, performing whole-cell patch-clamp recordings from perirhinal pyramidal neurons and GABAergic interneurons of GAD67-GFP⁺ mice, we found, for the first time, that the majority of PRC neurons are resonant at their resting potential, with a resonance frequency of 0.5–1.5 Hz at 23°C and of 1.5–2.8 Hz at 36°C. In the presence of ZD7288 (blocker of HCN channels) resonance was abolished in both pyramidal neurons and interneurons, suggesting that I_h current is critically involved in resonance generation. Otherwise, application of TTx (voltage-dependent Na⁺ channel blocker) attenuates the resonance in pyramidal neurons but not in interneurons, suggesting that only in pyramidal neurons the persistent sodium current has an amplifying effect. These experimental results have also been confirmed by a computational model. From a functional point of view, the resonance in the PRC would affect the reverberating activity between neocortex and hippocampus, especially during slow wave sleep, and could be involved in the redistribution and strengthening of memory representation in cortical regions.

Biomimetic FPGA-based spatial navigation model with grid cells and place cells

The mammalian spatial navigation system is characterized by an initial divergence of internal representations, with disparate classes of neurons responding to distinct features including location, speed, borders and head direction; an ensuing convergence finally enables navigation and path integration. Here, we report the algorithmic and hardware implementation of biomimetic neural structures encompassing a feed-forward trimodular, multi-layer architecture representing grid-cell, place-cell and decoding modules for navigation. The grid-cell module comprised of neurons that fired in a grid-like pattern, and was built of distinct layers that constituted the dorsoventral span of the medial entorhinal cortex. Each layer was built as an independent continuous attractor network with distinct grid-field spatial scales. The place-cell module comprised of neurons that fired at one or few spatial locations, organized into different clusters based on convergent modular inputs from different grid-cell layers, replicating the gradient in place-field size along the hippocampal dorso-ventral axis. The decoding module, a two-layer neural network that constitutes the convergence of the divergent representations in preceding modules, received inputs from the place-cell module and provided specific coordinates of the navigating object. After vital design optimizations involving all modules, we implemented the tri-modular structure on Zynq Ultrascale+ field-programmable gate array silicon chip, and demonstrated its capacity in precisely estimating the navigational trajectory with minimal overall resource consumption involving a mere 2.92% Look Up Table utilization. Our implementation of a biomimetic, digital spatial navigation system is stable, reliable, reconfigurable, real-time with execution time of about 32 s for 100k input samples (in contrast to 40 minutes on Intel Core i7-7700 CPU with 8 cores clocking at 3.60 GHz) and thus can be deployed for autonomous-robotic navigation without requiring additional sensors.

Kv11 (ether-à-go-go-related gene) voltage-dependent K+ channels promote resonance and oscillation of subthreshold membrane potentials

Key points

Some ion channels are known to behave as inductors and make up the parallel resonant circuit in the plasma membrane of neurons, which enables neurons to respond to current inputs with a specific frequency (so‐called ‘resonant properties’).

Here, we report that heterologous expression of mouse Kv11 voltage‐dependent K⁺ channels generate resonance and oscillation at depolarized membrane potentials in HEK293 cells; expressions of individual Kv11 subtypes generate resonance and oscillation with different frequency properties.

Kv11.3‐expressing HEK293 cells exhibited transient conductance changes that opposed the current changes induced by voltage steps; this probably enables Kv11 channels to behave like an inductor.

The resonance and oscillation of inferior olivary neurons were impaired at the resting membrane potential in Kv11.3 knockout mice.

This study helps to elucidate basic ion channel properties that are crucial for the frequency responses of neurons.

Abstract

The plasma membranes of some neurons preferentially respond to current inputs with a specific frequency, and output as large voltage changes. This property is called resonance, and is thought to be mediated by ion channels that show inductor‐like behaviour. However, details of the candidate ion channels remain unclear. In this study, we mainly focused on the functional roles of Kv11 potassium (K⁺) channels, encoded by ether‐á‐go‐go‐related genes, in resonance in mouse inferior olivary (IO) neurons. We transfected HEK293 cells with long or short splice variants of Kv11.1 (Merg1a and Merg1b) or Kv11.3, and examined membrane properties using whole‐cell recording. Transfection with Kv11 channels reproduced resonance at membrane potentials depolarized from the resting state. Frequency ranges of Kv11.3‐, Kv11.1(Merg1b)‐ and Kv11.1(Merg1a)‐expressing cells were 2–6 Hz, 2–4 Hz, and 0.6–0.8 Hz, respectively. Responses of Kv11.3 currents to step voltage changes were essentially similar to those of inductor currents in the resistor–inductor–capacitor circuit. Furthermore, Kv11 transfections generated membrane potential oscillations. We also confirmed the contribution of HCN1 channels as a major mediator of resonance at more hyperpolarized potentials by transfection into HEK293 cells. The Kv11 current kinetics and properties of Kv11‐dependent resonance suggested that Kv11.3 mediated resonance in IO neurons. This finding was confirmed by the impairment of resonance and oscillation at –30 to –60 mV in Kcnh7 (Kv11.3) knockout mice. These results suggest that Kv11 channels have important roles in inducing frequency‐dependent responses in a subtype‐dependent manner from resting to depolarized membrane potentials.

Some ion channels are known to behave as inductors and make up the parallel resonant circuit in the plasma membrane of neurons, which enables neurons to respond to current inputs with a specific frequency (so‐called ‘resonant properties’).

Here, we report that heterologous expression of mouse Kv11 voltage‐dependent K⁺ channels generate resonance and oscillation at depolarized membrane potentials in HEK293 cells; expressions of individual Kv11 subtypes generate resonance and oscillation with different frequency properties.

Kv11.3‐expressing HEK293 cells exhibited transient conductance changes that opposed the current changes induced by voltage steps; this probably enables Kv11 channels to behave like an inductor.

The resonance and oscillation of inferior olivary neurons were impaired at the resting membrane potential in Kv11.3 knockout mice.

This study helps to elucidate basic ion channel properties that are crucial for the frequency responses of neurons.

Frequency-dependent responses of neuronal models to oscillatory inputs in current versus voltage clamp

Action potential generation in neurons depends on a membrane potential threshold and therefore on how subthreshold inputs influence this voltage. In oscillatory networks, for example, many neuron types have been shown to produce membrane potential ([Formula: see text]) resonance: a maximum subthreshold response to oscillatory inputs at a nonzero frequency. Resonance is usually measured by recording [Formula: see text] in response to a sinusoidal current ([Formula: see text]), applied at different frequencies (f), an experimental setting known as current clamp (I-clamp). Several recent studies, however, use the voltage clamp (V-clamp) method to control [Formula: see text] with a sinusoidal input at different frequencies [[Formula: see text]] and measure the total membrane current ([Formula: see text]). The two methods obey systems of differential equations of different dimensionality, and while I-clamp provides a measure of electrical impedance [[Formula: see text]], V-clamp measures admittance [[Formula: see text]]. We analyze the relationship between these two measurement techniques. We show that, despite different dimensionality, in linear systems the two measures are equivalent: [Formula: see text]. However, nonlinear model neurons produce different values for Z and [Formula: see text]. In particular, nonlinearities in the voltage equation produce a much larger difference between these two quantities than those in equations of recovery variables that describe activation and inactivation kinetics. Neurons are inherently nonlinear, and notably, with ionic currents that amplify resonance, the voltage clamp technique severely underestimates the current clamp response. We demonstrate this difference experimentally using the PD neurons in the crab stomatogastric ganglion. These findings are instructive for researchers who explore cellular mechanisms of neuronal oscillations.

Membrane potential resonance in non-oscillatory neurons interacts with synaptic connectivity to produce network oscillations

Several neuron types have been shown to exhibit (subthreshold) membrane potential resonance (MPR), defined as the occurrence of a peak in their voltage amplitude response to oscillatory input currents at a preferred (resonant) frequency. MPR has been investigated both experimentally and theoretically. However, whether MPR is simply an epiphenomenon or it plays a functional role for the generation of neuronal network oscillations and how the latent time scales present in individual, non-oscillatory cells affect the properties of the oscillatory networks in which they are embedded are open questions. We address these issues by investigating a minimal network model consisting of (i) a non-oscillatory linear resonator (band-pass filter) with 2D dynamics, (ii) a passive cell (low-pass filter) with 1D linear dynamics, and (iii) nonlinear graded synaptic connections (excitatory or inhibitory) with instantaneous dynamics. We demonstrate that (i) the network oscillations crucially depend on the presence of MPR in the resonator, (ii) they are amplified by the network connectivity, (iii) they develop relaxation oscillations for high enough levels of mutual inhibition/excitation, and (iv) the network frequency monotonically depends on the resonators resonant frequency. We explain these phenomena using a reduced adapted version of the classical phase-plane analysis that helps uncovering the type of effective network nonlinearities that contribute to the generation of network oscillations. We extend our results to networks having cells with 2D dynamics. Our results have direct implications for network models of firing rate type and other biological oscillatory networks (e.g, biochemical, genetic).

Degeneracy in the robust expression of spectral selectivity, subthreshold oscillations, and intrinsic excitability of entorhinal stellate cells

Biological heterogeneities are ubiquitous and play critical roles in the emergence of physiology at multiple scales. Although neurons in layer II (LII) of the medial entorhinal cortex (MEC) express heterogeneities in channel properties, the impact of such heterogeneities on the robustness of their cellular-scale physiology has not been assessed. Here, we performed a 55-parameter stochastic search spanning nine voltage- or calcium-activated channels to assess the impact of channel heterogeneities on the concomitant emergence of 10 in vitro electrophysiological characteristics of LII stellate cells (SCs). We generated 150,000 models and found a heterogeneous subpopulation of 449 valid models to robustly match all electrophysiological signatures. We employed this heterogeneous population to demonstrate the emergence of cellular-scale degeneracy in SCs, whereby disparate parametric combinations expressing weak pairwise correlations resulted in similar models. We then assessed the impact of virtually knocking out each channel from all valid models and demonstrate that the mapping between channels and measurements was many-to-many, a critical requirement for the expression of degeneracy. Finally, we quantitatively predict that the spike-triggered average of SCs should be endowed with theta-frequency spectral selectivity and coincidence detection capabilities in the fast gamma-band. We postulate this fast gamma-band coincidence detection as an instance of cellular-scale-efficient coding, whereby SC response characteristics match the dominant oscillatory signals in LII MEC. The heterogeneous population of valid SC models built here unveils the robust emergence of cellular-scale physiology despite significant channel heterogeneities, and forms an efficacious substrate for evaluating the impact of biological heterogeneities on entorhinal network function.

NEW & NOTEWORTHY We assessed the impact of heterogeneities in channel properties on the robustness of cellular-scale physiology of medial entorhinal cortical stellate neurons. We demonstrate that neuronal models with disparate channel combinations were endowed with similar physiological characteristics, as a consequence of the many-to-many mapping between channel properties and the physiological characteristics that they modulate. We predict that the spike-triggered average of stellate cells should be endowed with theta-frequency spectral selectivity and fast gamma-band coincidence detection capabilities.

Subcellular Localization and Activity of TRPM4 in Medial Prefrontal Cortex Layer 2/3

TRPM4 is a Ca²⁺-activated non-selective cationic channel that conducts monovalent cations. TRPM4 has been proposed to contribute to burst firing and sustained activity in several brain regions, however, the cellular and subcellular pattern of TRPM4 expression in medial prefrontal cortex (mPFC) during postnatal development has not been elucidated. Here, we use multiplex immunofluorescence labeling of brain sections to characterize the postnatal developmental expression of TRPM4 in the mouse mPFC. We also performed electrophysiological recordings to correlate the expression of TRPM4 immunoreactivity with the presence of TRPM4-like currents. We found that TRPM4 is expressed from the first postnatal day, with expression increasing up to postnatal day 35. Additionally, in perforated patch clamp experiments, we found that TRPM4-like currents were active at resting membrane potentials at all postnatal ages studied. Moreover, TRPM4 is expressed in both pyramidal neurons and interneurons. TRPM4 expression is localized in the soma and proximal dendrites, but not in the axon initial segment of pyramidal neurons. This subcellular localization is consistent with a reduction in the basal current only when we locally perfused 9-Phenanthrol in the soma, but not upon perfusion in the medial or distal dendrites. Our results show a specific localization of TRPM4 expression in neurons in the mPFC and that a 9-Phenanthrol sensitive current is active at resting membrane potential, suggesting specific functional roles in mPFC neurons during postnatal development and in adulthood.

Resonance Properties in Auditory Brainstem Neurons

Auditory signals carry relevant information on a large range of time scales from below milliseconds to several seconds. Different stages in the auditory brainstem are specialized to extract information in specific frequency domains. One biophysical mechanism to facilitate frequency specific processing are membrane potential resonances. Here, we provide data from three different brainstem nuclei that all exhibit high-frequency subthreshold membrane resonances that are all most likely based on low-threshold potassium currents. Fitting a linear model, we argue that, as long as neurons possess active subthreshold channels, the main determinant for their resonance behavior is the steady state membrane time constant. Tuning this leak conductance can shift membrane resonance frequencies over more than a magnitude and therefore provide a flexible mechanism to tune frequency-specific auditory processing.

Spiking resonances in models with the same slow resonant and fast amplifying currents but different subthreshold dynamic properties

The generation of spiking resonances in neurons (preferred spiking responses to oscillatory inputs) requires the interplay of the intrinsic ionic currents that operate at the subthreshold voltage level and the spiking mechanisms. Combinations of the same types of ionic currents in different parameter regimes may give rise to different types of nonlinearities in the voltage equation (e.g., parabolic- and cubic-like), generating subthreshold (membrane potential) oscillations patterns with different properties. These nonlinearities are not apparent in the model equations, but can be uncovered by plotting the voltage nullclines in the phase-plane diagram. We investigate the spiking resonant properties of conductance-based models that are biophysically equivalent at the subthreshold level (same ionic currents), but dynamically different (parabolic- and cubic-like voltage nullclines). As a case study we consider a model having a persistent sodium and a hyperpolarization-activated (h-) currents, which exhibits subthreshold resonance in the theta frequency band. We unfold the concept of spiking resonance into evoked and output spiking resonance. The former focuses on the input frequencies that are able to generate spikes, while the latter focuses on the output spiking frequencies regardless of the input frequency that generated these spikes. A cell can exhibit one or both types of resonances. We also measure spiking phasonance, which is an extension of subthreshold phasonance (zero-phase-shift response to oscillatory inputs) to the spiking regime. The subthreshold resonant properties of both types of models are communicated to the spiking regime for low enough input amplitudes as the voltage response for the subthreshold resonant frequency band raises above threshold. For higher input amplitudes evoked spiking resonance is no longer present in these models, but output spiking resonance is present primarily in the parabolic-like model due to a cycle skipping mechanism (involving mixed-mode oscillations), while the cubic-like model shows a better 1:1 entrainment. We use dynamical systems tools to explain the underlying mechanisms and the mechanistic differences between the resonance types. Our results demonstrate that the effective time scales that operate at the subthreshold regime to generate intrinsic subthreshold oscillations, mixed-mode oscillations and subthreshold resonance do not necessarily determine the existence of a preferred spiking response to oscillatory inputs in the same frequency band. The results discussed in this paper highlight both the complexity of the suprathreshold responses to oscillatory inputs in neurons having resonant and amplifying currents with different time scales and the fact that the identity of the participating ionic currents is not enough to predict the resulting patterns, but additional dynamic information, captured by the geometric properties of the phase-space diagram, is needed.

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.

Resonance modulation, annihilation and generation of anti-resonance and anti-phasonance in 3D neuronal systems: interplay of resonant and amplifying currents with slow dynamics

Subthreshold (membrane potential) resonance and phasonance (preferred amplitude and zero-phase responses to oscillatory inputs) in single neurons arise from the interaction between positive and negative feedback effects provided by relatively fast amplifying currents and slower resonant currents. In 2D neuronal systems, amplifying currents are required to be slave to voltage (instantaneously fast) for these phenomena to occur. In higher dimensional systems, additional currents operating at various effective time scales may modulate and annihilate existing resonances and generate antiresonance (minimum amplitude response) and antiphasonance (zero-phase response with phase monotonic properties opposite to phasonance). We use mathematical modeling, numerical simulations and dynamical systems tools to investigate the mechanisms underlying these phenomena in 3D linear models, which are obtained as the linearization of biophysical (conductance-based) models. We characterize the parameter regimes for which the system exhibits the various types of behavior mentioned above in the rather general case in which the underlying 2D system exhibits resonance. We consider two cases: (i) the interplay of two resonant gating variables, and (ii) the interplay of one resonant and one amplifying gating variables. Increasing levels of an amplifying current cause (i) a response amplification if the amplifying current is faster than the resonant current, (ii) resonance and phasonance attenuation and annihilation if the amplifying and resonant currents have identical dynamics, and (iii) antiresonance and antiphasonance if the amplifying current is slower than the resonant current. We investigate the underlying mechanisms by extending the envelope-plane diagram approach developed in previous work (for 2D systems) to three dimensions to include the additional gating variable, and constructing the corresponding envelope curves in these envelope-space diagrams. We find that antiresonance and antiphasonance emerge as the result of an asymptotic boundary layer problem in the frequency domain created by the different balances between the intrinsic time constants of the cell and the input frequency f as it changes. For large enough values of f the envelope curves are quasi-2D and the impedance profile decreases with the input frequency. In contrast, for f ≪ 1 the dynamics are quasi-1D and the impedance profile increases above the limiting value in the other regime. Antiresonance is created because the continuity of the solution requires the impedance profile to connect the portions belonging to the two regimes. If in doing so the phase profile crosses the zero value, then antiphasonance is also generated.

Competition between Persistent Na+ and Muscarine-Sensitive K+ Currents Shapes Perithreshold Resonance and Spike Tuning in CA1 Pyramidal Neurons

Neurons from many brain regions display intrinsic subthreshold theta-resonance, responding preferentially to theta-frequency oscillatory stimuli. Resonance may contribute to selective communication among neurons and to orchestrate brain rhythms. CA1 pyramidal neurons receive theta activity, generating place fields. In these neurons the expression of perithreshold frequency preference is controversial, particularly in the spiking regime, with evidence favoring either non-resonant (integrator-like) or resonant behavior. Perithreshold dynamics depends on the persistent Na⁺ current I_NaP developing above −70 mV and the muscarine-sensitive K⁺ current I_M activating above −60 mV. We conducted current and voltage clamp experiments in slices to investigate perithreshold excitability of CA1 neurons under oscillatory stimulation. Around 20% of neurons displayed perithreshold resonance that is expressed in spiking. The remaining neurons (~80%) acted as low-pass filters lacking frequency preference. Paired voltage clamp measurement of I_NaP and I_M showed that perithreshold activation of I_M is in general low while I_NaP is high enough to depolarize neurons toward threshold before resonance expression, explaining the most abundant non-resonant perithreshold behavior. Partial blockade of I_NaP by pharmacological tools or dynamic clamp changed non-resonant to resonant behavior. Furthermore, shifting I_M activation toward hyperpolarized potentials by dynamic clamp also transformed non-resonant neurons into resonant ones. We propose that the relative levels of I_NaP and I_M control perithreshold behavior of CA1 neurons constituting a gating mechanism for theta resonance in the spiking regime. Both currents are regulated by intracellular signaling and neuromodulators which may allow dynamic switching of perithreshold behavior between resonant and non-resonant.

Hierarchical spike clustering analysis for investigation of interneuron heterogeneity

Action potentials represent the output of a neuron. Especially interneurons display a variety of discharge patterns ranging from regular action potential firing to prominent spike clustering or stuttering. The mechanisms underlying this heterogeneity remain incompletely understood. We established hierarchical cluster analysis of spike trains as a measure of spike clustering. A clustering index was calculated from action potential trains recorded in the whole-cell patch clamp configuration from hippocampal (CA1, stratum radiatum) and entorhinal (medial entorhinal cortex, layer 2) interneurons in acute slices and simulated data. Prominent, region-dependent, but also variable spike clustering was detected using this measure. Further analysis revealed a strong positive correlation between spike clustering and membrane potentials oscillations but an inverse correlation with neuronal resonance. Furthermore, clustering was more pronounced when the balance between fast-activating K(+) currents, assessed by the spike repolarisation time, and hyperpolarization-activated currents, gauged by the size of the sag potential, was shifted in favour of fast K(+) currents. Simulations of spike clustering confirmed that variable ratios of fast K(+) and hyperpolarization-activated currents could underlie different degrees of spike clustering and could thus be crucial for temporally structuring interneuron spike output.

Interplay between Subthreshold Oscillations and Depressing Synapses in Single Neurons

In this paper we analyze the interplay between the subthreshold oscillations of a single neuron conductance-based model and the short-term plasticity of a dynamic synapse with a depressing mechanism. In previous research, the computational properties of subthreshold oscillations and dynamic synapses have been studied separately. Our results show that dynamic synapses can influence different aspects of the dynamics of neuronal subthreshold oscillations. Factors such as maximum hyperpolarization level, oscillation amplitude and frequency or the resulting firing threshold are modulated by synaptic depression, which can even make subthreshold oscillations disappear. This influence reshapes the postsynaptic neuron’s resonant properties arising from subthreshold oscillations and leads to specific input/output relations. We also study the neuron’s response to another simultaneous input in the context of this modulation, and show a distinct contextual processing as a function of the depression, in particular for detection of signals through weak synapses. Intrinsic oscillations dynamics can be combined with the characteristic time scale of the modulatory input received by a dynamic synapse to build cost-effective cell/channel-specific information discrimination mechanisms, beyond simple resonances. In this regard, we discuss the functional implications of synaptic depression modulation on intrinsic subthreshold dynamics.

Active dendrites mediate stratified gamma-range coincidence detection in hippocampal model neurons

KEY POINTS

Quantitative metrics for the temporal window of integration/coincidence detection, based on the spike-triggered average, were employed to assess the emergence and dependence of gamma-range coincidence detection in hippocampal pyramidal neurons on various ion channel combinations. The presence of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels decreased the coincidence detection window (CDW) of the neuronal compartment to the gamma frequency range. Interaction of HCN channels with T-type calcium channels and persistent sodium channels further reduced the CDW, whereas interaction with A-type potassium channels broadened the CDW. When multiple channel gradients were co-expressed, the high density of resonating conductances in the distal dendrites led to a slow gamma CDW in the proximal dendrites and a fast-gamma CDW in the distal dendrites. The presence of resonating and spike-generating conductances serve as a mechanism underlying the emergence of stratified gamma-range coincidence detection in the dendrites of CA1 pyramidal neurons, enabling them to perform behaviour- and state-dependent gamma frequency multiplexing.

ABSTRACT

Hippocampal pyramidal neurons exhibit gamma-phase preference in their spikes, selectively route inputs through gamma frequency multiplexing and are considered part of gamma-bound cell assemblies. How do these neurons exhibit gamma-frequency coincidence detection capabilities, a feature that is essential for the expression of these physiological observations, despite their slow membrane time constant? In this conductance-based modelling study, we developed quantitative metrics for the temporal window of integration/coincidence detection based on the spike-triggered average (STA) of the neuronal compartment. We employed these metrics in conjunction with quantitative measures for spike initiation dynamics to assess the emergence and dependence of coincidence detection and STA spectral selectivity on various ion channel combinations. We found that the presence of resonating conductances (hyperpolarization-activated cyclic nucleotide-gated or T-type calcium), either independently or synergistically when expressed together, led to the emergence of spectral selectivity in the spike initiation dynamics and a significant reduction in the coincidence detection window (CDW). The presence of A-type potassium channels, along with resonating conductances, reduced the STA characteristic frequency and broadened the CDW, but persistent sodium channels sharpened the CDW by strengthening the spectral selectivity in the STA. Finally, in a morphologically precise model endowed with experimentally constrained channel gradients, we found that somatodendritic compartments expressed functional maps of strong theta-frequency selectivity in spike initiation dynamics and gamma-range CDW. Our results reveal the heavy expression of resonating and spike-generating conductances as the mechanism underlying the robust emergence of stratified gamma-range coincidence detection in the dendrites of hippocampal and cortical pyramidal neurons.

Subthreshold amplitude and phase resonance in models of quadratic type: nonlinear effects generated by the interplay of resonant and amplifying currents

We investigate the biophysical and dynamic mechanisms of generation of subthreshold amplitude and phase resonance in response to sinusoidal input currents in two-dimensional models of quadratic type. These models feature a parabolic voltage nullcline and a linear nullcline for the recovery gating variable, capturing the interplay of the so-called resonant currents (e.g., hyperpolarization-activated mixed-cation inward and slow potassium) and amplifying currents (e.g., persistent sodium) in biophysically realistic parameter regimes. These currents underlie the generation of resonance in medial entorhinal cortex layer II stellate cells and CA1 pyramidal cells. We show that quadratic models exhibit nonlinear amplifications of the voltage response to sinusoidal inputs in the resonant frequency band. These are expressed as an increase in the impedance profile as the input amplitude increases. They are stronger for values positive than negative to resting potential and are accompanied by a shift in the phase profile, a decrease in the resonant and phase-resonant frequencies, and an increase in the sharpness of the voltage response. These effects are more prominent for smaller values of ∊ (larger levels of the time scale separation between the voltage and the resonant gating variable) and for values of the resting potential closer to threshold for spike generation. All other parameter fixed, as ∊ increases the voltage response becomes "more linear"; i.e., the nonlinearities are present, but "ignored". In addition, the nonlinear effects are strongly modulated by the curvature of the parabolic voltage nullcline (partially reflecting the effects of the amplifying current) and the slope of the resonant current activation curve. Following the effects of changes in the biophysical conductances of realistic conductance-based models through the parameters of the quadratic model, we characterize the qualitatively different effects that resonant and amplifying currents have on the nonlinear properties of the voltage response. We identify different classes of resonant currents, represented by h- and slow potassium, according to whether they enhance (h-) or attenuate (slow potassium) the nonlinear effects. Finally, we use dynamical systems tools to investigate the dynamic mechanisms of generation of resonance and phase-resonance. We show that the nonlinear effects on the voltage response (e.g., amplification of the voltage response in the resonant frequency band and shifts in the resonant and phase-resonant frequencies) result from the ability of limit cycle trajectories to follow the unstable (right) branch of the voltage nullcline for a significant amount of time. This is a canard-related mechanism that has been shown to underlie the generation of intrinsic subthreshold oscillations in quadratic type models such as medial entorhinal cortex stellate cells. Overall, our results highlight the complexity of the voltage response to oscillatory inputs in nonlinear models and the roles that resonant and amplifying currents have in shaping these responses.

Membrane potential-dependent integration of synaptic inputs in entorhinal stellate neurons

Stellate cells (SCs) of the medial entorhinal cortex exhibit robust spontaneous membrane-potential oscillations (MPOs) in the theta (4-12 Hz) frequency band as well as theta-frequency resonance in their membrane impedance spectra. Past experimental and modeling work suggests that these features may contribute to the phase-locking of SCs to the entorhinal theta rhythm and may be important for forming the hexagonally tiled grid cell place fields exhibited by these neurons in vivo. Among the major biophysical mechanisms contributing to MPOs is a population of persistent (non-inactivating or slowly inactivating) sodium channels. The resulting persistent sodium conductance (GNaP ) gives rise to an apparent increase in input resistance as the cell approaches threshold. In this study, we used dynamic clamp to test the hypothesis that this increased input resistance gives rise to voltage-dependent, and thus MPO phase-dependent, changes in the amplitude of excitatory and inhibitory post-synaptic potential (PSP) amplitudes. We find that PSP amplitude depends on membrane potential, exhibiting a 5-10% increase in amplitude per mV depolarization. The effect is larger than-and sums quasi-linearly with-the effect of the synaptic driving force, V - Esyn . Given that input-driven MPOs 10 mV in amplitude are commonly observed in MEC stellate cells in vivo, this voltage- and phase-dependent synaptic gain is large enough to modulate PSP amplitude by over 50% during theta-frequency MPOs. Phase-dependent synaptic gain may therefore impact the phase locking and phase precession of grid cells in vivo to ongoing network oscillations. © 2014 Wiley Periodicals, Inc.

Expression and functional roles of Kv7/KCNQ/M-channels in rat medial entorhinal cortex layer II stellate cells

The medial entorhinal cortex (MEC) is important for spatial navigation and memory. Stellate cells (SCs) of MEC layer II provide major input to the hippocampus, and are thought to be the neuronal correlate of the grid cells. Their electrophysiological properties have been used to explain grid field formation. However, little is known about the functional roles of potassium channels in SCs. M-current is a slowly activating potassium current, active at subthreshold potentials. Although some studies have suggested that Kv7/M-channels may affect subthreshold resonance in SCs, others have found no Kv7/M-current in these cells, so the expression and roles of Kv7/M-channels in SCs are still debated. Using whole-cell voltage-clamp, we have identified a typical M-current with pharmacological properties characteristic of Kv7/M-channels in rat MEC SCs. Current-clamp experiments showed that the specific Kv7/M-channel blocker XE991 increased SCs excitability, and reduced spike frequency adaptation. Our results demonstrate that Kv7/M-channels are expressed in SCs and contribute substantially to regulation of excitability in these cells.

Electrical resonance in the θ frequency range in olfactory amygdala neurons

The cortical amygdala receives direct olfactory inputs and is thought to participate in processing and learning of biologically relevant olfactory cues. As for other brain structures implicated in learning, the principal neurons of the anterior cortical nucleus (ACo) exhibit intrinsic subthreshold membrane potential oscillations in the θ-frequency range. Here we show that nearly 50% of ACo layer II neurons also display electrical resonance, consisting of selective responsiveness to stimuli of a preferential frequency (2–6 Hz). Their impedance profile resembles an electrical band-pass filter with a peak at the preferred frequency, in contrast to the low-pass filter properties of other neurons. Most ACo resonant neurons displayed frequency preference along the whole subthreshold voltage range. We used pharmacological tools to identify the voltage-dependent conductances implicated in resonance. A hyperpolarization-activated cationic current depending on HCN channels underlies resonance at resting and hyperpolarized potentials; notably, this current also participates in resonance at depolarized subthreshold voltages. KV7/KCNQ K⁺ channels also contribute to resonant behavior at depolarized potentials, but not in all resonant cells. Moreover, resonance was strongly attenuated after blockade of voltage-dependent persistent Na⁺ channels, suggesting an amplifying role. Remarkably, resonant neurons presented a higher firing probability for stimuli of the preferred frequency. To fully understand the mechanisms underlying resonance in these neurons, we developed a comprehensive conductance-based model including the aforementioned and leak conductances, as well as Hodgkin and Huxley-type channels. The model reproduces the resonant impedance profile and our pharmacological results, allowing a quantitative evaluation of the contribution of each conductance to resonance. It also replicates selective spiking at the resonant frequency and allows a prediction of the temperature-dependent shift in resonance frequency. Our results provide a complete characterization of the resonant behavior of olfactory amygdala neurons and shed light on a putative mechanism for network activity coordination in the intact brain.