Abstract
Depending on the behavioral state, hippocampal CA1 pyramidal neurons receive very distinct patterns of synaptic input and likewise produce very different output patterns. We have used simultaneous dendritic and somatic recordings and multisite glutamate uncaging to investigate the relationship between synaptic input pattern, the form of dendritic integration, and action potential output in CA1 neurons. We found that when synaptic input arrives asynchronously or highly distributed in space, the dendritic arbor performs a linear integration that allows the action potential rate and timing to vary as a function of the quantity of the input. In contrast, when synaptic input arrives synchronously and spatially clustered, the dendritic compartment receiving the clustered input produces a highly nonlinear integration that leads to an action potential output that is extraordinarily precise and invariant. We also present evidence that both of these forms of information processing may be independently engaged during the two distinct behavioral states of the hippocampus such that individual CA1 pyramidal neurons could perform two different state-dependent computations: input strength encoding during theta states and feature detection during sharp waves.
Introduction
The rodent hippocampus exhibits two major behaviorally related states, each displaying a characteristic neuronal population pattern. One of these states occurs as an animal actively explores the environment or during rapid eye movement sleep and is characterized by theta frequency (5–10 Hz) network oscillations (Vanderwolf, 1969; O’Keefe and Nadel, 1978). The other is a nontheta state that occurs during awake-immobility and slow-wave sleep and is characterized by irregular large amplitude oscillations [sharp waves (SPWs)] and higher frequency (150–200 Hz) synchronous discharges (ripples) (O’Keefe and Nadel, 1978; Buzsaki, 1986).
These separable behavioral states are also associated with distinct input–output patterns in the various hippocampal subregions, with the network appearing to operate at two distinctive timescales. Indeed, synaptic input to CA1 pyramidal neurons during theta is robust but spread out over hundreds of milliseconds, whereas, during SPWs, input synchrony is estimated to increase ∼10-fold (Csicsvari et al., 1998, 2000). The output patterns produced during these two states are also quite different. During theta, the relative timing of action potential (AP) output advances over many tens of milliseconds as the animal moves toward and through the place field of a particular neuron (O’Keefe and Dostrovsky, 1971; O’Keefe and Recce, 1993; Harris et al., 2002; Mehta et al., 2002). In contrast, an extremely precise and invariant form of output timing is found during SPWs, with AP firing found almost exclusively during a 1 ms time window of the ripple trough (Csicsvari et al., 1999).
We were interested in investigating the idea that a fundamental change in dendritic integration of synaptic input could provide a cellular-level mechanism supporting the state-dependent network activity. CA1 and other pyramidal neurons have been reported to be capable of performing either a linear form of integration (Margulis and Tang, 1998; Cash and Yuste, 1999; Nettleton and Spain, 2000; Tamas et al., 2002) or a nonlinear form stemming from the generation of dendritic spikes (Golding and Spruston, 1998; Larkum et al., 2001; Wei et al., 2001; Ariav et al., 2003; Gasparini et al., 2004; Polsky et al., 2004). Both of these integration styles appear to result from complex interactions among the passive and active dendritic membrane properties, as well as the properties of the synapses themselves. Some current hypotheses have compressed these two integration forms into a single sigmoid-shaped function with linear integration occurring simply because input levels are subthreshold for local spike generation (Poirazi et al., 2003; Polsky et al., 2004; Spruston and Kath, 2004). We propose here, however, that linear and nonlinear dendritic integration can be two fundamentally distinct integration forms, with each selectively engaged by input with different temporal structure. This situation could provide the network with a degree of flexibility that might be useful in producing the dramatically different timescales of network firing found during theta and SPWs. We have used multisite two-photon glutamate uncaging and simultaneous whole-cell recordings from the dendrites and soma of CA1 pyramidal neurons in acute hippocampal slices to determine whether different input patterns can produce distinct single-cell computations most complementary for a given behavioral state.
Materials and Methods
Transverse hippocampal slices (400 μm thick) were prepared from 8- to 12-week old Sprague Dawley rats as described previously (Magee, 1998), according to methods approved by the Louisiana State University Health Science Center Institutional Animal Care and Use Committee. Rats were given a lethal dose of ketamine and xylazine, perfused through the ascending aorta with an oxygenated solution just before death, and decapitated.
Dendrites from pyramidal CA1 pyramidal cells were visualized using a Zeiss (Oberkochen, Germany) Axioskop fit with differential interference contrast (DIC) optics under infrared illumination. The external solution used for recordings contained the following (in mm): 125 NaCl, 3 KCl, 25 NaHCO3, 1.25 NaHPO4, 1.3 CaCl2, 1 MgCl2, and 25 glucose, and was saturated with 95% O2 and 5% CO2 at 34–36°C. Simultaneous whole-cell patch-clamp recordings from distal dendrites and somata were performed using Dagan (Minneapolis, MN) BVC-700 amplifiers in the active “bridge” mode, filtered at 1–10 kHz and digitized at 50 kHz. Patch pipettes had a resistance of 4–6 MΩ for dendritic and 2–4 MΩ for somatic recordings when filled with a solution containing the following (in mm): 120 K-methylsulfate, 20 KCl, 10 HEPES, 0.5 EGTA, 4 NaCl, 4 MgATP, 0.3 Tris2GTP, 14 phosphocreatine, pH 7.25. The series resistance was usually <20 MΩ. The distance of the dendritic electrode from the visually identified soma was measured under DIC conditions.
EPSC-shaped waveforms characterized by AMPA current kinetics (τrise, 0.1 ms; τdecay, 2 ms) (Magee and Cook, 2000; Smith et al., 2003) were injected through the whole-cell electrode. In particular, we injected either 10 identical EPSCs of increasing amplitude at different frequencies (100 Hz to 10 kHz) or a random pattern of 12–260 EPSCs (of 80 pA each) dispersed in 1–100 ms. The value of 80 pA was chosen, because it generated a voltage deflection of 0.8–1 mV when injected in the distal dendrites, as reported for the single synaptic event at these locations (Magee and Cook, 2000). In a set of experiments, a real-time conductance injection amplifier (SM-1; Cambridge Conductance, Cambridge, UK) (Harsch and Robinson, 2000) was used to stimulate the neurons by delivering a combination of AMPA and NMDA conductance changes (AMPA/NMDA amplitude ratio was set at 2) (McDermott et al., 2005). In this case, two dendritic electrodes were used, to independently inject current and record the Vm. The amplifier reproduced the voltage dependence of the block by Mg2+ by multiplying the linear conductance equation by a Boltzmann-type nonlinearity (Harsch and Robinson, 2000), where K1 and K2 were set to 0.3 and 0.06/mV, respectively. τrise and τdecay for NMDA currents were set to 1 and 20 ms, respectively, and the reversal potential for AMPA and NMDA currents was set to 0 mV.
The presentation of variable spatiotemporal input patterns was performed using a dual galvanometer-based scanning system to multiphoton photorelease glutamate at multiple dendritic spines and simultaneously image local Ca2+ in the dendrites of CA1 neurons (Prairie Technologies, Middleton, WI). Neurons were filled with Oregon Green BAPTA-1 (OGB-1) (Molecular Probes, Eugene, OR) (standard internal solution with EGTA replaced by 100 μm OGB-1) and imaged using a 60× objective (Olympus, Melville, NY). Ultrafast, pulsed, laser light at 930 nm (Mira 900F; Coherent, Auburn, CA) was used to excite the OGB-1, while 720 nm ultrafast, pulsed, laser light (Chameleon; Coherent, Auburn, CA) was used to photolyze 4-methoxy-7-nitroindolinyl-caged l-glutamate (MNI-Glu) (Tocris Cookson, Ballwin, MO) (10 mm applied via pipette above slice) (Smith et al., 2003). To quantify input relative to unitary synaptic input for this region, MNI-Glu was uncaged at 7–10 identified spines (0.2 ms exposure) with a 100 ms uncaging interval at a given laser intensity. This exposure duration produced glutamate uncaging-induced EPSPs (glu-EPSPs) whose kinetics matched closely those reported for unitary excitatory input to this region. glu-EPSPs within 10 μm of the recording electrode had a 20–80% rise time of 0.5 ± 0.04 ms and a decay time constant of 8.3 ± 0.8 ms (n = 22), comparable to a ∼0.6 ms rise time and 7 ms decay time constant for distal-trunk miniature EPSPs (Magee and Cook, 2000). Locations were chosen to give a relatively uniform glu-EPSP amplitude (2–3 mV at the dendrite; ∼0.5 mV at the soma) representing approximately three unitary synaptic inputs (each of 0.8 mV at dendrite, 0.2 mV for the soma) (Magee and Cook, 2000) for each location, giving a total of 20–30 for all 7–10 locations. Laser power was gradually increased to ultimately provide the equivalent of 60–120 synapses worth of input (where each individual glu-EPSP had an amplitude of 5–10 mV at the dendrite or 1–2 mV at the soma) (see Fig. 2b). Two different temporal patterns (0.3–0.4 or 5.2 ms uncaging intervals) were then given at the respective laser power intensities, and this procedure was repeated for the different spatial distributions (all input spines within a dendritic region of 20 or 150 μm). Our uncaging scanner was measured to make a highly reproducible step across the entire imaging field (∼130 μm) in <300 μs. We generally used three to four steps to cover the full distance of the spatial profile with the remaining four to seven steps used to fill in the intervening space randomly (see Fig. 2a). The extraordinarily high spine density of the distal trunk region (∼7/μm) made finding more than one or two adequately isolated spines impossible (Smith et al., 2003). This means that there was a significant glutamate spillover between some spines, with a larger than normal NMDA receptor-mediated slow EPSP component resulting. The functional impact of this was a larger amount of summation during the asynchronous input patterns (observable as an approximately twofold larger slope in Fig. 2e than in Fig. 1e).
Data are reported as means ± SEM. Statistical comparisons were performed by using Student’s t test. Means were considered to be significantly different when p < 0.05.
Results
Dendritic integration mode is dependent on the temporal input pattern
To investigate the impact of temporal input pattern on dendritic integration, a random pattern of 12–260 synaptic inputs were delivered over varying durations of time (1–100 ms) through a dendritic electrode placed at >250 μm from the soma (Fig. 1a) (unitary conductance or current produced an ∼1 mV dendritic EPSP that is similar in size and shape to the mean unitary EPSP found in this dendritic region) (see Materials and Methods) (Magee and Cook, 2000). Membrane potential from either nearby dendrite, soma, or both was simultaneously recorded with a second or third electrode.
When progressively increasing amounts of asynchronous input were delivered to the dendrite (input spread over time windows ranging from 10 to 100 ms), the level of resulting depolarization increased in an approximately linear manner, with a slope that was dependent on the amount of input synchrony (Fig. 1a,e). Eventually, enough input was delivered to generate AP output. For all cases of asynchronous input patterns (10–100 ms), APs were first initiated in the soma/axon region of the neuron and then backpropagated into the dendrite (bAP) (Fig. 1c). Although the arrival of bAPs makes a precise quantification difficult, it appeared that dendritic depolarizations continue to increase as a linear function of input well beyond the output AP threshold (Fig. 1a,e,f). Thus, as long as the input pattern was spread over tens of milliseconds, the dendrite integrated this input in a remarkably linear manner, producing progressively larger, slow depolarizations that directly reflected the number of inputs. The slow depolarizations propagated well to the soma where they initiated an action potential output at the soma/axon region. These data show that asynchronous input is associated with a linear integration mode that evokes action potential output from the soma/axon region.
A different pattern of synaptic input produced a different form of dendritic depolarization. When increasing amounts of highly synchronized input (all input arriving within 1–3 ms) were delivered to the dendrite, progressively larger, quite rapid, dendritic depolarizations were produced (Fig. 1b,e). Eventually, however, enough synchronous input was delivered to evoke a local dendritic spike that clearly occurred before the initiation of an output spike at the soma/axon (Fig. 1d) (dendritic spike threshold was 52 ± 5 inputs for a temporal window of <3 ms). The occurrence of a dendritic spike produced a disproportionate increase in dendritic depolarization that resulted in a sigmoid-shaped local input–output curve (Fig. 1e,f). We thus define this dendritic depolarization profile as a nonlinear form of integration.
Because the dendritic spike threshold in CA1 pyramidal neurons is heavily dependent on the rate of depolarization (dV/dt), only tightly synchronized input can produce the rapid depolarization required for local dendritic spike initiation (Gasparini et al., 2004). Thus, the highly nonlinear form of integration produced by local spike generation is seen only in response to synchronous input patterns. Indeed, we never observed the initiation of dendritic spikes for the more asynchronous input patterns even when very large amounts of input were delivered to the dendrite (e.g., >220 inputs over 50 ms producing dendritic depolarizations to approximately −25 mV). It should be noted that this level of depolarization is above threshold for both axonal APs and dendritic spikes evoked by synchronous input. Together, then, the data indicate that the form of dendritic integration (either linear or nonlinear) (Fig. 1f) is highly dependent on the timing of the input, with the dendrite performing an approximately linear integration for all input patterns except those that are extraordinarily well synchronized.
Dendritic integration mode is dependent on the spatial input pattern
To investigate how the spatial distribution of the input affects dendritic integration, we used multisite two-photon glutamate uncaging to deliver varying amounts of input in a spatial pattern that was either clustered or distributed. An input was considered to be clustered when a pattern of 7–10 uncaging events were delivered to spines covering ∼20 μm of dendrite, whereas the distributed input was delivered to spines spread over 125–150 μm of dendrite. To mimic a progressive increase in the number of synaptic inputs for each of these spatial patterns, laser intensity was gradually enhanced (see Materials and Methods), and the resulting membrane potential deflections were recorded with a dendritic pipette located within or near the input. The impact of spatial distribution was determined for both synchronous (all input within 3 ms) and asynchronous (all input within 50 ms) input patterns.
For asynchronous input patterns, the effect of spatial distribution was minimal, with the only observable impact being a statistically insignificant decrease in the slope of the local input–output curve for the more dispersed inputs (clustered, 0.33 ± 0.03 mV/input, n = 5; distributed, 0.28 ± 0.03 mV/input, n = 4) (Fig. 2c,e). For synchronous input patterns, the impact of spatial distribution was more dramatic. When all input was delivered to spines covering an ∼20 μm region of the dendrite, a strongly nonlinear form of integration was observed as the highly synchronous and clustered input initiated a powerful local dendritic spike (amplitude, 64 ± 3 mV; maximum spike dV/dt, 57 ± 5 V/s; n = 8) (Fig. 2d,e, blue traces). When the same temporally synchronous input pattern was delivered to spines in a distributed spatial pattern, the degree of nonlinear integration generated was either greatly reduced or removed completely (amplitude, 52 ± 5 mV; maximum spike dV/dt, 36 ± 4 V/s; n = 6) (Fig. 2d,e, red traces). Thus, for synchronous input patterns, the mode of dendritic integration was highly dependent on the degree of spatial clustering of the input, with the dendrite producing a strong nonlinearity only if the input arrived within a spatial distribution of <150 μm.
These data indicate that, depending on the spatiotemporal pattern of the synaptic input, either a linear or nonlinear form of dendritic depolarization occurs before the initiation of an output AP at the soma/axon region. We refer to these different dendritic depolarization forms as two distinct input pattern-dependent modes of integration, which are produced as the different spatiotemporal properties of the synaptic input differentially engage the active properties of CA1 pyramidal dendrites. It seems likely that the differential shaping of the incoming synaptic input will in turn cause unique patterns of AP output to be associated with the linear and nonlinear integration modes.
Input efficacy depends on the dendritic integration mode
To investigate the impact of integration form on AP output, we first determined how the efficacy of the input, as indicated by the number of synaptic inputs required for axonal output, varied as a function of the different input patterns. Different input patterns were delivered as above (Fig. 3b), and we found that the number of synaptic inputs required for the generation of an axonal output was reduced for more synchronous when compared with less synchronous input patterns (53 ± 3, n = 18, for 1–3 ms duration vs 126 ± 6 for 50 ms duration, n = 18) (Fig. 3c–g). Furthermore, we observed that distributing a highly synchronous input over 125–150 μm of dendrite reduced the efficacy of the input toward that of the less synchronous input (clustered, 53 ± 3 inputs, vs distributed, 79 ± 4 inputs; n = 8 and 6, respectively) (Fig. 3f,g). These data show that CA1 pyramidal neurons preferentially respond to synchronous input patterns only if this input is spatially clustered enough to shift the dendrites into the nonlinear integration mode. Presumably, this effect is attributable to the initiation of a dendritic spike by this pattern of input.
We tested the idea that the generation of a strong dendritic spike is required for increasing the efficacy of a synchronous and clustered input by reducing dendritic excitability through bath application of a low concentration of tetrodotoxin (TTX) (50 nm). Such low concentrations of TTX have repeatedly been shown to effectively dampen dendritic excitability without greatly affecting that of the somatic/axonal region (somatic action potential threshold for 1 ms input: control, −54 ± 1 mV; TTX, −49 ± 1 mV; for 50 ms input, control, −52 ± 1 mV; TTX, −47 ± 1 mV) (Magee and Carruth, 1999; Metz et al., 2005). TTX application completely eliminated the ability of the dendrite to produce the supralinear integration mode that is usually associated with synchronous input without changing the dendritic depolarization profile produced by less synchronous input (data not shown).
In this low TTX condition, the number of inputs required for an axonal output dramatically increased for the most synchronous input pattern (132 ± 10 inputs within 1 ms; n = 9) but increased only slightly for the less synchronous patterns (162 ± 16 for the 50 ms duration; n = 12) (Fig. 3g). A similar TTX-dependent decrease in the efficacy of synchronous synaptic input has been observed in vivo (Canals et al., 2005). These data indicate that it is the robust propagation of the dendritic spike, not simply the synchronous input pattern itself, that is required for the generation of an axonal spike at lower input levels. In fact, the initiation and robust propagation of the dendritic spike make highly synchronous input the preferred pattern, despite the severe low-pass filtering properties of dendrites. Furthermore, the strict spatial requirements for dendritic spike generation allow CA1 pyramidal neurons to select for those patterns that are both synchronous and clustered. All other patterns will be less effective in producing an action potential output, because they require significantly greater amounts of input.
We next examined the impact that increasing the amount of input above threshold has on AP output. As expected, the number of APs evoked by asynchronous input patterns increased as the level of input was enhanced beyond threshold (Fig. 3b–d,h). Conversely, the initiation of a dendritic spike by the synchronous input pattern almost always resulted in only a single AP generated at the soma even at input levels that were twice threshold. These data indicate that the maintenance of the linear integration mode above AP output threshold allows the number of action potentials initiated by CA1 pyramidal neurons to steadily change as a function of input levels during asynchronous or spatially distributed input patterns. In contrast, the transition into the nonlinear integration mode greatly reduced the variation in spike rate such that only a single action potential output was evoked in response to the highly synchronous and clustered input.
Output precision is dependent on dendritic integration mode
We found another aspect of action potential output that was dependent on the mode of dendritic integration. Not surprisingly, the latency of somatic action potential output decreases with increased input synchrony (from 41.5 ± 1.3 for the 50 ms interval to 2.8 ± 0.1 ms for the 1 ms duration; n = 12) (Fig. 4a–d). Furthermore, the less synchronous input produced an action potential output whose latency was quite variable from trial to trial, especially compared with that produced by the most synchronous and clustered input pattern. Indeed, the jitter (defined as the SD of spike latency) associated with the synchronous patterns was more than an order of magnitude less than that of the asynchronous or distributed input patterns (2.1 ± 0.4 ms, n = 11, for the 50 ms duration, vs 0.06 ± 0.01 ms, n = 12, for the 1 ms duration) (Fig. 4e). A similarly low jitter has been reported for APs initiated by basal dendritic spikes in CA1 pyramidal neurons once these more weakly propagating spikes are paired with subthreshold depolarizations (Ariav et al., 2003).
As above, inhibiting the ability of the dendrites to produce a nonlinear integration by spatially distributing synchronous input or by bath-applying low TTX greatly affected AP latency and precision (Fig. 4c). TTX (50 nm) more than doubled the latency of action potential output induced by synchronous input (6.1 ± 0.3 ms; n = 6) (Fig. 4e) while only slightly increasing AP latency for the more asynchronous patterns (47.6 ± 0.9 ms; n = 13). Likewise, the extremely low latency variation observed during highly synchronized input patterns was completely eliminated in the presence of 50 nm TTX (the jitter for the 1 ms duration increased to 0.9 ± 0.1 ms; n = 6) (Fig. 4e). The increased precision of action potential output was also abolished with input patterns that were spread over 125–150 μm (jitter increased from 0.06 ± 0.01 to 1.4 ± 0.1 ms; n = 6 and 5, respectively) (Fig. 4c,e). Again, this indicates that it is not simply the synchronized input pattern itself, but the shift into the nonlinear dendritic integration mode that is required for the extremely precise action potential output pattern associated with this input pattern. Indeed, it seems likely that it is the extraordinarily rapid depolarization of the soma/axon region (10’s mV/ms) (Gasparini et al., 2004) by the actively propagating dendritic spike that allows the remarkably precise spike output. Thus, without the strong nonlinear dendritic integration mode, the neuron loses the uniquely precise action potential output normally associated with the arrival of highly synchronous and well clustered input patterns.
Finally, we again observed that only during asynchronous input patterns was the latency to first action potential firing capable of varying significantly as a function of the number of inputs (Fig. 4f). Indeed, the ability of the neuron to change spike timing as a function of input amount was severely reduced if the nonlinear integration mode was engaged by highly synchronous and clustered input. Thus, spike timing in CA1 pyramidal neurons can significantly vary as a function of the amount of input only when the dendrites integrate incoming input in a linear manner. If, in contrast, the input is highly clustered and synchronous such that a dendritic spike is initiated, AP output timing will be remarkably insensitive to variations in input strength.
Thus far, we have seen that, depending on the spatiotemporal input pattern, the dendrites of CA1 pyramidal neurons can perform two different forms of integration and that a unique action potential output pattern is associated with each of these integration modes. The linear integration mode was most frequently engaged (by a wide variety of less synchronous or spatially distributed input patterns), and it was associated with an output pattern in which the number and timing of evoked action potentials can vary as a function of input intensity. The other nonlinear dendritic integration mode was produced by synaptic input patterns that are both highly synchronized and spatially clustered and was associated with an output pattern characterized by a reduced output threshold and short latency, low jitter action potentials that were relatively unresponsive to changes in the number of inputs. We next wanted to examine the functional impact of these two distinct dendritic integration modes.
Dendritic integration mode during different behavioral states
The hippocampus exhibits two distinct behavioral states, and as detailed below, each of these states is associated with a unique input–output pattern in area CA1. We were interested in determining the role that the input-dependent modification of dendritic integration plays in producing the appropriate network patterns observed during each of these unique behavioral states. To this end, we presented CA1 pyramidal neurons with input patterns similar to those reported to be found during either theta or SPW states, determined whether the appropriate output pattern was produced, and examined the role of the dendritic integration mode in producing this output.
Theta states
During theta states, the synaptic input to CA1 pyramidal neurons increases with the movement of the animal toward and through the place field of that neuron (Harris et al., 2002; Mehta et al., 2002). In response, these neurons increase their firing rates and shift the timing of the spikes progressively forward over many tens of milliseconds (phase precession) (O’Keefe and Dostrovsky, 1971; O’Keefe and Recce, 1993). The magnitude and temporal pattern of synaptic input during this state can be estimated from the firing of the presynaptic CA3 pyramidal neuron population. The number of CA3 pyramidal cells that are firing during the 100–200 ms of a theta cycle has been estimated to be 3–5% of the total population (Csicsvari et al., 1998, 2000). Given the anatomical organization of the Schaffer collateral pathway, this works out to be ∼600 potential synaptic inputs (4% of 15,000 apical Schaffer collateral synapses) (Bannister and Larkman, 1995; Megias et al., 2001). Other factors, such as the total number of action potentials fired per cycle (from one to four per cycle) (Csicsvari et al., 2000) and the low probability of release at each synapse (0.1–0.5) (Raastad et al., 1992; Allen and Stevens, 1994; Smith et al., 2003), make it likely that a group of presynaptic CA3 pyramidal neurons provides a given CA1 pyramidal neuron with an input that ranges from <100 to ∼500 inputs per cycle as the animal moves toward and through the place field of that cell. Simultaneous with this pattern of excitation, the proximal portions of CA1 pyramidal neurons (the soma and proximal trunk) receive a rhythmic 5–10 Hz input from the inhibitory basket cell and axo-axonic cell network (Fox, 1989; Cobb et al., 1995, Klausberger et al., 2003), whereas the dendrites receive a relatively sparse (∼3% of excitatory synapse numbers) but nevertheless potentially effective input from the bistratified and O-LM (oriens-lacunosum moleculare) interneuron populations that is ∼180° out of phase with the axo-axonic and basket cell input (Klausberger et al., 2004).
We therefore examined the pattern of action potential output resulting from the combination of a progressively increasing asynchronous excitatory input to the dendrite with an out-of-phase 10 Hz somatic hyperpolarizing sine wave to mimic proximal inhibitory input (amplitude was adjusted to increase approximately twofold the threshold amount of excitatory input) (Fig. 5a,b) (Harris et al., 2002; Mehta et al., 2002). The phasic dendritic input pattern was set up to mimic the progressive increase in excitation from 80 to 400 synaptic inputs, phasically interrupted by a like-sized dendritic inhibitory input (Klausberger et al., 2003, 2004). Generally, by the third cycle, the dendritic depolarization started to overcome the somatic inhibition and action potential output was generated (for 240 synapses activated in the example in Fig. 5a). As the input levels increased, the level of dendritic depolarization also increased accordingly as did the number of output APs per cycle (Fig. 5a,d,e). Together with this progressive increase in spike firing, the mean spike timing smoothly moved forward with respect to the theta phase demonstrating the characteristic phase precession found during theta states (Fig. 5a,d,e). Thus, the injection of asynchronous input patterns similar to those found during theta states left the dendritic arbor in a linear integration mode (only bAPs were observed) (Fig. 5c) and produced an output whose rate and coarse timing smoothly varied as a function of input intensity, just as occurs in the behaving animal (O’Keefe and Dostrovsky, 1971; O’Keefe and Recce, 1993; Harris et al., 2002; Mehta et al., 2002). These data indicate that, when paired with a phasic proximal inhibitory input, the linear dendritic integration mode is uniquely capable of producing an action potential output whose rate and coarse timing can directly reflect the amount of synaptic input received by the neuron.
Interestingly, in a number of cases (67% of neurons and 38% of the trials within those neurons), bAP duration increased significantly during the periods of highest input because of what is presumably an enhanced activation of dendritic voltage-gated Ca2+ channels (Fig. 5c, asterisks). These “Ca2+ shoulders” appeared to propagate to the soma where they produced an afterdepolarizing potential and burst firing (Golding et al., 1999; Magee and Carruth, 1999). Such an event is reminiscent of BAC (backpropagation-activated calcium spike) firing in layer V pyramidal cells (Larkum et al., 2001), except that here it requires more than a single spike (more than three bAPs) and is only likely to occur as the dendritic A-type K+ current is inactivated by prolonged dendritic depolarizations (Hoffman et al., 1997; Magee and Johnston, 1997). The role of this firing pattern in theta activity will be further explored in future experiments.
Sharp-wave states
CA1 input–output features are quite different during sharp wave/ripple activity. Here, the level of input synchrony is elevated almost tenfold, with ∼4% of the CA3 pyramidal cell population providing input to CA1 within a 20 ms time window (Csicsvari et al., 2000). This synchronous input pattern is also associated with an action potential output reported to occur almost exclusively within a 1 ms time window during the troughs of the extracellularly recorded ripple event (∼200 Hz) of SPW states (jitter, ∼0.3 ms) (Csicsvari et al., 1999). In an attempt to reproduce this in vivo situation, we paired a dendritic excitatory input of 60–100 synaptic inputs arriving within 1–3 ms (1.5% of 15,000 synapses with a mean Prelease ∼0.3) with a 200 Hz hyperpolarizing sine wave at the soma. Again the somatic sine wave is used as an approximation of the phasic arrival of proximal inhibition recorded during the ripple events (Fig. 6a1) (Buzsaki, 1986; Ylinen et al., 1995; Klausberger et al., 2003, 2004). This synchronous input pattern evoked a dendritic spike producing a well timed output (all output occurred within a 1 ms time window) that was characterized by a short latency and high precision (jitter, 0.04 ± 0.01 ms; n = 7) (Fig. 6a2–a5). Importantly, this high level of precision was relatively insensitive to modifications that might be expected to occur quite commonly during in vivo conditions. That is, the variability in the timing of the output APs with respect to the 200 Hz sine was still very low even with changes in resting membrane (±5 mV) and the number of inputs (increases to twofold above threshold) (jitter, 0.16 ± 0.04 ms; n = 5) (Fig. 6a). The jitter in AP timing recorded under these conditions approaches that reported for in vivo recordings during SPWs (0.2 vs ∼0.3 ms) (Csicsvari et al., 1999).
The recruitment of the nonlinear integrative properties of the dendrite by a synchronous input was absolutely required to produce the extremely well timed and invariant AP output during the simulated SPW states (Fig. 6a). When an asynchronous input (180 synapses activated within 50 ms) was injected into the dendrites along with the 200 Hz hyperpolarizing somatic input, the threshold for spike generation was reached at the soma (Fig. 6b). As a consequence of the lower precision associated with the linear integration mode, the output spiking was spread over three cycles instead of a single (Fig. 6b3). Even when the spike timing was plotted relative to the phase of an individual cycle (to remove the cycle variations), the degree of temporal dispersion was much greater than that achieved during the synchronous input (compare Fig. 6b4,b5 with a4,a5). These data indicate that a phasic inhibitory input alone, at least when mimicked as a 200 Hz sine wave, is not capable of providing the reproducibly invariant AP timing reported to be present during SPWs. Additionally, a relatively synchronous input pattern that was nevertheless subthreshold for dendritic spike initiation (50–90 inputs over 4 ms) delivered to the dendrite in combination with the 200 Hz hyperpolarizing sine was also incapable of achieving the invariant AP output timing observed during dendritic spike initiation. AP timing during this input pattern varied across 100° as the amount of input was varied from 50 to 90 inputs (Fig. 6c4,c5) and over 65° with fluctuations in resting membrane potential (±5 mV) (data not shown).
Finally, the strict input pattern requirements of the nonlinear integration mode can be observed to provide a degree of functionality to a neuron that might be especially relevant to the precise sequential firing observed in the CA1 network during SPWs. When a neuron is presented with a spatially clustered and synchronous input pattern, together with a 200 Hz phasic proximal hyperpolarization, action potential output probability is high (∼0.8) and a well timed output is observed (Fig. 7b,d, input labeled “clustered”). However, if a similarly sized input is distributed across spines covering ∼150 μm of dendrite, the input efficacy is reduced (Fig. 3) and action potential output probability is likewise significantly inhibited (∼0.15) (Fig. 7d). Furthermore, the precision of the spike timing for these few APs is quite poor and similar to that seen in Figure 6c. These data indicate that CA1 neurons can use the strict spatiotemporal requirements for dendritic spike generation to distinguish between multiple synchronous inputs, responding with a precise AP output only to input that is both synchronous and clustered. This would be useful if a particular CA1 neuron receives a succession of synchronous input from the different groups of CA3 pyramidal cells being sequentially read out during a single SPW/ripple event (Lee and Wilson, 2002). The ability to ignore altogether input that is not properly clustered (or at most produce an infrequent poorly timed output) will enhance the ability of the neuron to fire in the proper temporal order during the rapid network sequence replay.
Together, these data indicate that behaviorally relevant types of input patterns are capable of shifting the dendritic integration mode of CA1 pyramidal neurons and that this shift plays an important role in producing output patterns that resemble quite closely those associated with each of the two major behavioral states seen in the hippocampus. In light of these data, it seems likely that the characteristic network patterns associated with theta and SPW states are in part produced by this capability of CA1 pyramidal dendrites to change the way they integrate incoming input. Without this, it seems unlikely that the neurons embedded within the CA3–CA1 network could operate under such disparate timescales as those observed during theta and SPW states.
Discussion
In summary, we have directly demonstrated that different spatiotemporal patterns of synaptic input differentially engage the active properties of CA1 pyramidal dendrites to produce two distinct forms of integration. Furthermore, each of these dendritic integration modes induces a unique style of action potential output with significant functional relevance. Finally, we have attempted to associate this relationship between input patterns, dendritic integration mode, and output generation with the two main behavioral states found in the hippocampus. In this respect, we have shown that asynchronous input patterns similar to those found during theta states summate linearly in the dendrites and soma, ultimately producing an action potential output whose rate and coarse timing is a direct function of the amount of input. This is an input–output relationship that closely approximates what is observed in the rodent hippocampus during active exploration of space. In contrast, synchronous inputs that are also clustered in space summate supralinearly through the initiation of a dendritic spike that then propagates actively to the soma/axon. This actively propagating dendritic spike evokes a precise axonal output that reflects the spatiotemporal pattern of the input. This integration mode appears to be uniquely capable of producing the precise level of spike timing required by the much-compressed timescales present during sharp wave states. There is in fact some evidence of dendritic spiking during SPWs (Kamondi et al., 1998). Together, these data indicate that CA1 pyramidal neurons can perform two very different single-cell computations, depending on the pattern of activity within the network of presynaptic CA3 cells (Fig. 8).
Dendritic properties involved in linear and nonlinear integration modes
The dendritic and synaptic properties of CA1 pyramidal neurons allow them to integrate synaptic inputs in either an essentially linear or highly nonlinear manner depending on the type of stimuli received. The linear integration mode is remarkable in that it is an “active linearization” achieved as Na+, NMDA, Ih, and K+ channel activations effectively counterbalance each other and the loss of AMPA current driving force over a window of dendritic Vm spanning approximately from −70 to −25 mV. Indeed, blockade of any one of these channel populations dramatically alters the dendritic current–voltage relationship (Cash and Yuste, 1999; Gasparini et al., 2004). In this study, we have used the distal apical trunk as a model dendritic compartment, because its morphology allows for direct electrical recordings and because of the availability of high-quality published data. Preliminary data from our lab support the idea that other dendritic regions, such as the radial oblique and basal branches, possess very similar functional properties (Magee and Losonczy, 2004).
The ability of a particular input pattern to produce the transition from linear to supralinear integration depends on the rate of membrane depolarization produced by the input. Indeed, local dendritic spike threshold becomes more depolarized as dV/dt decreases, eventually reaching such depolarized levels that no dendritic spikes can be initiated (Gasparini et al., 2004). This property, combined with the small amplitude and short duration of unitary EPSPs in the dendrite (<1 mV; <10 ms half-width) (Magee and Cook, 2000; Smith et al., 2003) determine that only highly synchronized and spatially clustered inputs will provide the level of rapid dendritic depolarization required (Turner et al., 1991; Golding and Spruston, 1998). In contrast, as long as the pattern of input is relatively distributed in space or time, the lower and more stable AP threshold at the output region will be exceeded. Thus, very specific input conditions are required to produce the shift from linear to nonlinear integration, and unless these are met, the arbor will usually operate in a linear integration mode.
An additional factor to be considered is the neuromodulatory state of the neuron. It is currently thought that theta states are associated with relatively high levels of ACh, whereas norepinephrine levels have been reported to be high during SPW activity (for review, see Pace-Schott and Hobson, 2002). Because dendritic ion channels are susceptible to these neuromodulators (Tsubokawa and Ross, 1997; Colbert and Johnston, 1998; Hoffman and Johnston, 1998; Gasparini and Magee, 2002; Magee and Johnston, 2005), the exact spatiotemporal patterns of the input required to shift the processing mode may vary from one behavioral state to another. Moreover, activity-dependent processes that modify dendritic active channels (Frick et al., 2004) could produce similar alterations in the pattern dependence of dendritic integration.
Finally, dendritic spike propagation has been reported to be variable, with some spikes propagating more strongly than others (Golding and Spruston, 1998; Ariav et al., 2003; Gasparini et al., 2004). Obviously, weakly propagating local spikes will not be as effective as stronger spikes in shifting the output pattern of the neuron. However, because the forward propagation of dendritic spikes is very sensitive to subthreshold levels of depolarization (Ariav et al., 2003; Gasparini et al., 2004; Jarsky et al., 2005), background excitatory input may provide a mechanism to equalize the impact of most dendritic spikes.
Dendritic computations during different behavioral states
The distinct input–output patterns found during theta and SPW states are indicative of the very different functions thought to occur during the respective behavioral states. It is currently believed that the temporal sequences of a particular episode occurring over seconds of time are encoded and transiently stored within the hippocampal network during theta states (McNaughton and Morris, 1987; Mehta et al., 2000; Dragoi et al., 2003). During subsequent SPW states, the previously stored information seems to be “read out” to the higher cortical areas for consolidation as a tightly compressed replay (now over only tens of milliseconds) of essentially the same neuronal firing sequences activated during previous theta states (Qin et al., 1997; Siapas and Wilson, 1998; Kudrimoti et al., 1999; Nadasdy et al., 1999; Lee and Wilson, 2002). It is our hypothesis that, as the hippocampus moves from one state to another, a switch in the mode of dendritic integration occurs that allows CA1 pyramidal neurons (and perhaps other hippocampal neurons) to perform the computation that is most appropriate for the unique timescales present during each behavioral state.
Our data suggest that, during theta states, the regions of the neuron receiving Schaffer collateral input should operate in a linear integration mode with all input contributing to the output of the cell (Magee, 2000). Linear integration models have been shown to be quite effective in producing output patterns (changes in firing frequency) that accurately reflect changes in the quantity of inputs (termed “power variation”) (Poirazi et al., 2003). The addition of proximal phasic inhibitory input permits CA1 pyramidal neurons to also use spike timing to encode the location of the animal in space (Harris et al., 2002; Mehta et al., 2002; Huxter et al., 2003). Thus, during theta states, the pyramidal cell network may be in an integration mode that allows the cells to calculate power variations that encode the sequence of an animal’s behavior within the rate and course timing of neuronal firing over timescales that are dictated by the actual behavior (many hundreds of milliseconds). This type of cellular activity may be most advantageous for producing the synaptic and dendritic plasticity required for transient information storage within the hippocampal network. The functioning of these neurons during this time should not be confused with the concept of a “point neuron” simply because the dendrites are linearly integrating synaptic input. Indeed the presence of regionally localized inhibition and excitation as well as the complex interplay between backpropagating spikes and synaptic depolarization are obvious interactions among proximal and distal compartments that would not be present within a point neuron.
Processing in the nonlinear mode could functionally separate the CA1 arbor into a large number of independent nonlinear computational units, each sending its own output to the soma (Poirazi and Mel, 2001; Wei et al., 2001). This form of integration would allow neurons to respond not so much to the quantity of input but instead to the specific spatiotemporal input pattern (termed “configuration variability”) (Poirazi et al., 2003). During the period of heightened synchrony found during SPW states, the lack of an actual constraining behavior allows the network to process essentially the same information concerning the sequence of events ∼20 times more rapidly (tens of milliseconds) (Nadasdy et al., 1999; Lee and Wilson, 2002). The challenge to the individual CA1 neuron within the network is to reproducibly fire APs during a brief time window within the correct ripple phase so that the appropriate firing order among the neurons encoding the sequence can be maintained. Furthermore, this precise firing pattern has to occur many times over as the same neuronal sequences are repeatedly replayed (Lee and Wilson, 2002). The features provided by the nonlinear integration mode allow the postsynaptic CA1 pyramidal neuron to properly respond to the demands imposed on it during SPW states. That is, the cell is not likely to fire until it receives the appropriately timed and clustered input pattern (the proper configuration), at which point the neuron will respond with an extraordinarily well timed and invariant output. The rapid replay and exquisite network timing found during SPWs may facilitate the induction of plasticity in downstream cortical structures, perhaps through the induction of dendritic spiking and nonlinear dendritic integration in those cells.
In conclusion, there appears to be a direct link between the specific synaptic input pattern, the mode of dendritic integration, and ultimately the form of action potential output evoked in hippocampal CA1 pyramidal neurons. Furthermore, the features of the output deriving from the different integration modes appear to be reminiscent of those found during the two main hippocampal behavioral states. This link between input, integration, and output may allow the most appropriate single-cell computation to be performed during each behavioral state.
Footnotes
-
This work was supported by National Institute of Health Grant NS39458 (J.C.M.) and National Alliance for Research on Schizophrenia and Depression Grants YI2002 (S.G.) and II2003 (J.C.M.). We thank Gyorgy Buzsaki, Dan Johnston, and Randy Chitwood for many helpful discussions. We also express our gratitude to all of our friends and colleagues who offered us much needed kindness and support following the Hurricane Katrina disaster.
- Correspondence should be addressed to Jeffrey C. Magee, Neuroscience Center, Louisiana State University Health Science Center, 2020 Gravier Street, New Orleans, LA 70112. Email: jmagee{at}lsuhsc.edu