Abstract
Timing is a crucial aspect of synaptic integration. For pyramidal neurons that integrate thousands of synaptic inputs spread across hundreds of microns, it is thus a challenge to maintain the timing of incoming inputs at the axo-somatic integration site. Here we show that pyramidal neurons in the rodent hippocampus use a gradient of inductance in the form of hyperpolarization-activated cation-nonselective (HCN) channels as an active mechanism to counteract location-dependent temporal differences of dendritic inputs at the soma. Using simultaneous multi-site whole-cell recordings complemented by computational modeling, we find that this intrinsic biophysical mechanism produces temporal synchrony of rhythmic inputs in the theta and gamma frequency ranges across wide regions of the dendritic tree. While gamma and theta oscillations are known to synchronize activity across space in neuronal networks, our results identify a new mechanism by which this synchrony extends to activity within single pyramidal neurons with complex dendritic arbors.
This is a preview of subscription content, access via your institution
Access options
Subscribe to this journal
Receive 12 print issues and online access
$209.00 per year
only $17.42 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
References
Singer, W. & Gray, C.M. Visual feature integration and the temporal correlation hypothesis. Annu. Rev. Neurosci. 18, 555–586 (1995).
Engel, A.K., Fries, P. & Singer, W. Dynamic predictions: oscillations and synchrony in top-down processing. Nat. Rev. Neurosci. 2, 704–716 (2001).
Buzsáki, G. Rhythms of the Brain (Oxford University Press, 2011).
Colgin, L.L. et al. Frequency of gamma oscillations routes flow of information in the hippocampus. Nature 462, 353–357 (2009).
Bragin, A. et al. Gamma (40–100 Hz) oscillation in the hippocampus of the behaving rat. J. Neurosci. 15, 47–60 (1995).
Harris, K.D., Csicsvari, J., Hirase, H., Dragoi, G. & Buzsáki, G. Organization of cell assemblies in the hippocampus. Nature 424, 552–556 (2003).
Chrobak, J.J., Lorincz, A. & Buzsaki, G. Physiological patterns in the hippocampo-entorhinal cortex system. Hippocampus 10, 457–465 (2000).
O'Keefe, J. & Recce, M.L. Phase relationship between hippocampal place units and the EEG theta rhythm. Hippocampus 3, 317–330 (1993).
Skaggs, W.E., McNaughton, B.L., Wilson, M.A. & Barnes, C.A. Theta phase precession in hippocampal neuronal populations and the compression of temporal sequences. Hippocampus 6, 149–172 (1996).
Megías, M., Emri, Z., Freund, T.F. & Gulyás, A.I. Total number and distribution of inhibitory and excitatory synapses on hippocampal CA1 pyramidal cells. Neuroscience 102, 527–540 (2001).
Dougherty, K.A., Islam, T. & Johnston, D. Intrinsic excitability of CA1 pyramidal neurones from the rat dorsal and ventral hippocampus. J. Physiol. (Lond.) 590, 5707–5722 (2012).
Rall, W. Distinguishing theoretical synaptic potentials computed for different soma-dendritic distributions of synaptic input. J. Neurophysiol. 30, 1138–1168 (1967).
Alexander, C. Fundamentals of Electric Circuits 2nd edn. (McGraw Hill, 2004).
Mauro, A. Anomalous impedance, a phenomenological property of time-variant resistance. An analytical review. Biophys. J. 1, 353–372 (1961).
Cole, K.S. Transverse impedance of the squid giant axon during current flow. J. Gen. Physiol. 24, 535–549 (1941).
Cole, K.S. Membranes, Ions and Impulses (University of California Press, 1968).
Magee, J.C. Dendritic hyperpolarization-activated currents modify the integrative properties of hippocampal CA1 pyramidal neurons. J. Neurosci. 18, 7613–7624 (1998).
Robinson, R.B. & Siegelbaum, S.A. Hyperpolarization-activated cation currents: from molecules to physiological function. Annu. Rev. Physiol. 65, 453–480 (2003).
Hutcheon, B. & Yarom, Y. Resonance, oscillation and the intrinsic frequency preferences of neurons. Trends Neurosci. 23, 216–222 (2000).
Leung, L.S. & Yim, C.Y. Intracellular records of theta rhythm in hippocampal CA1 cells of the rat. Brain Res. 367, 323–327 (1986).
Pike, F.G. et al. Distinct frequency preferences of different types of rat hippocampal neurones in response to oscillatory input currents. J. Physiol. (Lond.) 529, 205–213 (2000).
Hu, H., Vervaeke, K. & Storm, J.F. Two forms of electrical resonance at theta frequencies, generated by M-current, h-current and persistent Na+ current in rat hippocampal pyramidal cells. J. Physiol. (Lond.) 545, 783–805 (2002).
Narayanan, R. & Johnston, D. The h channel mediates location dependence and plasticity of intrinsic phase response in rat hippocampal neurons. J. Neurosci. 28, 5846–5860 (2008).
Narayanan, R. & Johnston, D. Long-term potentiation in rat hippocampal neurons is accompanied by spatially widespread changes in intrinsic oscillatory dynamics and excitability. Neuron 56, 1061–1075 (2007).
Ulrich, D. Dendritic resonance in rat neocortical pyramidal cells. J. Neurophysiol. 87, 2753–2759 (2002).
Cook, E.P., Guest, J., Liang, Y., Masse, N. & Colbert, C. Dendrite-to-soma input/output function of continuous time-varying signals in hippocampal CA pyramidal neurons. J. Neurophysiol. 98, 2943–2955 (2007).
Hu, H., Vervaeke, K., Graham, L.J. & Storm, J.F. Complementary theta resonance filtering by two spatially segregated mechanisms in CA1 hippocampal pyramidal neurons. J. Neurosci. 29, 14472–14483 (2009).
Golding, N.L. Factors mediating powerful voltage attenuation along CA1 pyramidal neuron dendrites. J. Physiol. (Lond.) 568, 69–82 (2005).
Magee, J.C. Dendritic Ih normalizes temporal summation in hippocampal CA1 neurons. Nat. Neurosci. 2, 508–514 (1999).
Williams, S.R. & Stuart, G.J. Site independence of EPSP time course is mediated by dendritic Ih in neocortical pyramidal neurons. J. Neurophysiol. 83, 3177–3182 (2000).
Angelo, K., London, M., Christensen, S.R. & Häusser, M. Local and global effects of I(h) distribution in dendrites of mammalian neurons. J. Neurosci. 27, 8643–8653 (2007).
Colgin, L.L. & Moser, E.I. Gamma oscillations in the hippocampus. Physiology (Bethesda) 25, 319–329 (2010).
Jensen, O. & Colgin, L.L. Cross-frequency coupling between neuronal oscillations. Trends Cogn. Sci. 11, 267–269 (2007).
Kamondi, A., Acsady, L., Wang, X.J. & Buzsaki, G. Theta oscillations in somata and dendrites of hippocampal pyramidal cells in vivo: activity-dependent phase-precession of action potentials. Hippocampus 8, 244–261 (1998).
Harvey, C.D., Collman, F., Dombeck, D.A. & Tank, D.W. Intracellular dynamics of hippocampal place cells during virtual navigation. Nature 461, 941–946 (2009).
Magee, J.C. & Cook, E.P. Somatic EPSP amplitude is independent of synapse location in hippocampal pyramidal neurons. Nat. Neurosci. 3, 895–903 (2000).
Oren, I., Mann, E.O., Paulsen, O. & Hajos, N. Synaptic currents in anatomically identified CA3 neurons during hippocampal gamma oscillations in vitro. J. Neurosci. 26, 9923–9934 (2006).
Gasparini, S. & Magee, J.C. State-dependent dendritic computation in hippocampal CA1 pyramidal neurons. J. Neurosci. 26, 2088–2100 (2006).
Takahashi, H. & Magee, J.C. Pathway interactions and synaptic plasticity in the dendritic tuft regions of CA1 pyramidal neurons. Neuron 62, 102–111 (2009).
Fan, Y. et al. Activity-dependent decrease of excitability in rat hippocampal neurons through increases in I(h). Nat. Neurosci. 8, 1542–1551 (2005).
Brager, D.H. & Johnston, D. Plasticity of intrinsic excitability during long-term depression is mediated through mGluR-dependent changes in Ih in hippocampal CA1 pyramidal neurons. J. Neurosci. 27, 13926–13937 (2007).
Cassenaer, S. & Laurent, G. Hebbian STDP in mushroom bodies facilitates the synchronous flow of olfactory information in locusts. Nature 448, 709–713 (2007).
Wehr, M. & Laurent, G. Odour encoding by temporal sequences of firing in oscillating neural assemblies. Nature 384, 162–166 (1996).
Hahnloser, R.H.R., Kozhevnikov, A.A. & Fee, M.S. An ultra-sparse code underlies the generation of neural sequences in a songbird. Nature 419, 65–70 (2002).
Buzsáki, G. Neural syntax: cell assemblies, synapsembles, and readers. Neuron 68, 362–385 (2010).
Buzsáki, G. Theta rhythm of navigation: link between path integration and landmark navigation, episodic and semantic memory. Hippocampus 15, 827–840 (2005).
Buzsáki, G. & Moser, E.I. Memory, navigation and theta rhythm in the hippocampal-entorhinal system. Nat. Neurosci. 16, 130–138 (2013).
Chevaleyre, V. & Castillo, P.E. Assessing the role of Ih channels in synaptic transmission and mossy fiber LTP. Proc. Natl. Acad. Sci. USA 99, 9538–9543 (2002).
Carnevale, N.T. & Hines, M.L. The NEURON Book (Cambridge Univ. Press, 2009).
Gasparini, S., Migliore, M. & Magee, J.C. On the initiation and propagation of dendritic spikes in CA1 pyramidal neurons. J. Neurosci. 24, 11046–11056 (2004).
Shah, M.M., Migliore, M., Valencia, I., Cooper, E.C. & Brown, D.A. Functional significance of axonal Kv7 channels in hippocampal pyramidal neurons. Proc. Natl. Acad. Sci. USA 105, 7869–7874 (2008).
Routh, B.N., Johnston, D., Harris, K. & Chitwood, R.A. Anatomical and electrophysiological comparison of CA1 pyramidal neurons of the rat and mouse. J. Neurophysiol. 102, 2288–2302 (2009).
Kemenes, I. et al. Dynamic clamp with StdpC software. Nat. Protoc. 6, 405–417 (2011).
Glickfeld, L.L. & Scanziani, M. Distinct timing in the activity of cannabinoid-sensitive and cannabinoid-insensitive basket cells. Nat. Neurosci. 9, 807–815 (2006).
Torrence, C. & Compo, G. A practical guide to wavelet analysis. Bull. Am. Meteorol. Soc. 79, 61–78 (1998).
Liu, Y., San Liang, X. & Weisberg, R. Rectification of the bias in the wavelet power spectrum. J. Atmos. Ocean. Technol. 24, 2093–2102 (2007).
Acknowledgements
We thank R. Chitwood, N. Dembrow, R. Gray and R. Narayanan for helpful discussions during the course of this study. We also thank members of the Johnston laboratory and L.L. Colgin for comments on earlier versions of this manuscript. This work was supported by grant MH 048432 from US National Institutes of Health to D.J.
Author information
Authors and Affiliations
Contributions
S.P.V. & D.J. designed the experiments, interpreted the results and wrote the manuscript. S.P.V. performed the experiments, computer simulations and analysis of data.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Integrated supplementary information
Supplementary Figure 1 Theta frequency oscillatory synchrony is observed only at the soma
(a) The dendritic voltage response (normalized for amplitude) to a sinusoidal current of 7 Hz/100 pA (---) into the dendrite (red) and at the soma (black). Note that both responses are measured in the dendrite and do not show synchrony of response as is observed at the soma. They show distinct location-dependence differences in both control condition and with ZD7288; (b) Summary data for experiments described as in (a) with response latency defined as the difference in time between current injection at input site and voltage response at 300 mm from soma (in ms). [Control: Dendrite 3.59±1.15 and Soma 9.83±1.33 (n=8; ***, p=1.45e–06, power(1–β)=0.995, paired t-test); ZD7288: Dendrite 14.21±0.79 and Soma 21.23±1.11 (n=6; ***,p=6.34e-05, power(1–β)=0.999, paired t-test).] (c) Impedance phase profile (ZPP) for voltage response at 300 mm measured by 2s/100pA sinusoids of varying frequency when injected at 300 mm (red) and at the soma (black). Note the absence of synchrony (profiles do not intersect); (d) same as (c) but with ZD7288. (e) Summary data for experiments in (c,d) but described as difference in phase between the somatic and the dendritic responses. Negative values suggest dendritic input precedes the somatic input at 300 mm and vice-versa. Scale bars: 10 ms. Dendritic recording locations same as in Fig.1.
Supplementary Figure 2 ZPPsoma comparison between CA1 neuron, the ball-and-stick model and morphologically realistic model
(a) depicts the impedance phase profile at soma (ZPPsoma) for dendritic input at 300 μm (red) and for local somatic input (black) measured from a representative dual whole-cell recording in a CA1 neuron with HCN channels blocked with 20 μM ZD7288. (b) depicts the same measurements of ZPPsoma in normal conditions i.e in the presence of active HCN conductance. (c) and (d) correspond to (a) and (b), respectively, but in the simple ball-and-stick model (inset). Similarly, (e) and (f) correspond to (a) and (b), respectively, but in the morphologically realistic model (inset).
Supplementary Figure 3 Gradient of inductance achieves oscillatory synchrony with less voltage attenuation
(a) & (b) depict the difference in the phase between the dendritic input and the somatic input when measured at the soma. Negative values indicate the appearance of the somatic input before the dendritic input and vice versa. All responses are measured at 7 Hz for an increasing conductance density of leak in (a) and HCN in (b). Values on the side denote the maximum conductance density in S/cm2 in a sigmoidal distribution of the conductances. Note that the dashed line represents perfect synchrony where there is no difference in phase at the soma for input at any location along the dendrite. (c) summarizes the data in (a),(b) for an input at 300 mm (vertical dotted lines). (d) represents the transfer impedance amplitude for the same data points represented in (c).
Supplementary Figure 4 The spatial distribution of HCN channels achieves maximum transfer at synchronization frequency
(a),(b)&(c) depict the difference in the phase between the dendritic input and the somatic input when measured at the soma. Negative values indicate the appearance of the somatic input before the dendritic input and vice versa. All responses are measured at 7 Hz for an increasing conductance density of HCN channels for sigmoidal distribution in (a), linear in (b) and uniform in (c). Values on the side denote the maximum conductance density in S/cm2 in the distal dendritic region of the model. Note that the dashed line represents perfect synchrony where there is no difference in phase at the soma for input at any location along the dendrite. (d) shows the total HCN conductance for a sigmoidal, linear and uniform gradient with maximum conductance density of 10–3 S/cm2, which achieves perfect synchrony in (a),(b) and (c). (e) shows that linear and sigmoidal gradients have significantly higher transfer impedance at synchronization frequency than that in a uniform gradient irrespective of the maximum conductance density of the distribution. (f) shows that transfer impedance amplitude is maximum at synchronization frequency in sigmoidal or linear distribution but not in the unifrom distribution. (e)&(f) describe impedance measurements from 300 mm in the dendrite to the soma.
Supplementary Figure 5 Dependence of oscillatory synchrony on the voltage dependence of activation and kinetics of HCN channels
(a) shows the alterations in the voltage dependence of activation (V1/2) of the HCN conductance tested. (d) describes the scaling of the voltage dependence of activation/deactivation time constants. (b) & (e) depict the difference in the phase between the dendritic input and the somatic input when measured at the soma for a sinusoid of 7 Hz. Negative values indicate the appearance of the somatic input before the dendritic input and vice versa. (c) & (f) depict the influence of these alterations on synchronization frequency. Note that variations in V1/2 significantly altered oscillatory synchrony at 7 Hz in (b) and synchronization frequency in (c) while scaling of the activation/deactivation time constants did not (e,f). HCN conductance in this simulation had a sigmoidal distribution with a maximum conductance value of 1 x 10-3 S/cm2.
Supplementary Figure 6 Voltage dependence of SyncFreq
(a) depicts the voltage dependence of synchronization frequency in a model neuron with maximum HCN conductance of 0.25 x 10–3 S/cm2 (sigmoidal distribution). Note that Synchronization Frequency is restricted to sub-threshold potentials. (b) & (c) show the ZPPsoma at two distinct voltages in (a). Note the presence of oscillatory synchrony at –68 mV in (b) and its absence at –57 mV in (c). (d) shows that the voltage dependence of synchronization frequency is also dependent on the HCN conductance in the model neuron. The numbers depict the HCN conductance in S/cm2.
Supplementary Figure 7 Addition of ZD7288 alters the band-pass filtering of SyncFreqs
(a) & (b) describe the synaptic current injected by the dynamic clamp system at 300 mm to simulate a single synaptic event in (a) and a burst of 5 events at 60 Hz in (b), along with the corresponding voltage waveform recorded at soma in 4 different cells after application of 20 μM ZD7288. (c),(d),(e) ⊕ (f) analyze the frequency components in the corresponding dendritic currents or somatic voltage waveforms for traces depicted in (a) and (b). Note that although the frequency components in the synaptic current input are similar to that in control (Fig. 6), the voltage waveform at the soma is now composed of low-frequency components that do not show synchrony at the soma (Fig. 1d).
Supplementary Figure 8 Accuracy of dynamic clamp over current clamp for high-frequency synaptic inputs
A model cell was injected with synaptic current described either by an alpha-EPSC function in current clamp (red) or by calculating conductance using a single exponential decay model in the dynamic clamp system (black). (a) Conductance amplitude and exponential decay was adjusted in the dynamic clamp system so that both injections produced a similar EPSP waveform. (b), (c) & (d) show that with higher burst frequencies, there are significant differences between current injected by the two systems reflected here in the voltage amplitudes. This is because the dynamic clamp system accounts for the driving force during multiple opening events which is unaccounted for in the current clamp thus injecting a more realistic current waveform.
Supplementary Figure 9 Measurement of impedance phase profile (ZPP) by chirp versus sinusoids
ZPP measurements from 5 different cells are depicted where the solid lines indicate ZPP measured by a chirp stimulus (0 to 15 Hz in 15 s) and ⊕ depicts Impedance phase measured by a 2s sinusoid at a given frequency (0 to 15 Hz; increment 1 Hz). Red indicates measurements of ZPP when the stimulus was injected in the dendrite and black represents ZPP when the stimulus injection was at the soma. Note that no stimulus dependent differences were found in the measurement of ZPP.
Supplementary Figure 10 Time-frequency analysis using continuous wavelet transform
(a)–(d) depicts the time domain signals from Fig. 6 and their scaleogram (Time-Frequency Analysis) obtained by the continuous wavelet transform method. Note that the lower frequencies have a higher spectral resolution but bad temporal resolution and higher frequencies have lower spectral resolution but good temporal resolution. Traces on the right depict the time series average of the scaleogram.
Supplementary information
Supplementary Text and Figures
Supplementary Figures 1–10 (PDF 3805 kb)
Rights and permissions
About this article
Cite this article
Vaidya, S., Johnston, D. Temporal synchrony and gamma-to-theta power conversion in the dendrites of CA1 pyramidal neurons. Nat Neurosci 16, 1812–1820 (2013). https://doi.org/10.1038/nn.3562
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nn.3562
This article is cited by
-
Power-efficient neural network with artificial dendrites
Nature Nanotechnology (2020)
-
PRMT7 deficiency causes dysregulation of the HCN channels in the CA1 pyramidal cells and impairment of social behaviors
Experimental & Molecular Medicine (2020)
-
Robust emergence of sharply tuned place-cell responses in hippocampal neurons with structural and biophysical heterogeneities
Brain Structure and Function (2020)
-
Immediate neurophysiological effects of transcranial electrical stimulation
Nature Communications (2018)
-
Transient potassium channels augment degeneracy in hippocampal active dendritic spectral tuning
Scientific Reports (2016)