Abstract
Neurons are continually exposed to background synaptic activity in vivo. This is thought to influence neural information processing, but background levels of excitation and inhibition remain controversial. Here we show, using whole-cell recordings in anesthetized rats, that spontaneous depolarizations (“Up states”) in neocortical pyramidal neurons are driven by sparse, mostly excitatory synaptic activity (less than five inputs per millisecond; ∼10% inhibitory). The mean synaptic conductance change is small (<10 nS at the soma) and opposed by anomalous rectification, resulting in a net increase in input resistance during Up states. These conditions enhance the effectiveness of each synapse at depolarized potentials. Hence, neocortical networks are relatively quiet at rest, and the effect of synaptic background is weaker than previously thought.
Introduction
The principal task of every neuron is to interpret the patterns of synaptic input arriving in the dendritic tree and generate appropriate action potential (AP) output, a process termed synaptic integration (Magee, 2000). In the CNS, unitary EPSPs are small, so many inputs must summate to reach AP threshold (Häusser et al., 2000; Magee, 2000). Spontaneous activity within surrounding neural networks makes synaptic integration in vivo even more difficult by producing a noisy background of ongoing synaptic input (Bernander et al., 1991; Destexhe and Paré, 1999; Steriade, 2001; Destexhe et al., 2003). Background activity causes membrane depolarization, often in the form of periodically occurring “Up states” (Steriade et al., 1993, 2001; Wilson and Kawaguchi, 1996), reducing the driving force for EPSCs. In addition, synaptic input can shunt the membrane, reducing input resistance and EPSP amplitudes (Bernander et al., 1991), and this has led to the suggestion that Up states are “high-conductance states” (Destexhe et al., 2003)
Neurons have developed properties to diminish or overcome this shunting problem. These properties include the spatial distribution of synapses throughout the dendritic tree, which compartmentalizes conductance changes (Polsky et al., 2004; Williams, 2004), and the expression of voltage-gated ion channels, which can counteract synaptically induced conductance changes by increasing input resistance with depolarization (anomalous rectification) (Katz, 1949; Wilson, 1992).
Although shunting, synaptic compartmentalization and rectification are well known phenomena, how they interact to determine net input resistance and EPSP summation in the intact brain is unclear. In part, this results from a lack of information about the numbers of background synapses active in vivo. Here, we examine spontaneous background synaptic activity in the neocortex of the anesthetized rat. We show that somatic input resistance increases during Up states and demonstrate that this increase is attributable to anomalous rectification in these neurons. We then estimate the mean numbers of active excitatory and inhibitory synapses. We find that far fewer synapses are spontaneously active in cortical networks than previously thought and that the majority of these active synapses are excitatory.
Materials and Methods
In vivo electrophysiology
All experimental procedures were performed in accordance with the animal welfare guidelines of the Max Planck Society. Wistar rats [postnatal day 27 (P27) to P32] were anesthetized with urethane (1–2 g/kg). Depth of anesthesia was sufficient to eliminate pinch withdrawal, corneal reflex, and vibrissal movements. A small (2 × 2 mm) craniotomy was opened over barrel cortex, and the dura was removed (coordinates for center of craniotomy: 2 mm anterior to bregma, 5–6 mm lateral). The craniotomy was covered with agar (1–1.5%, type III-A; Sigma, St. Louis, MO) in the following solution (in mm): 135 NaCl, 5.4 KCl, 1 MgCl2, 1.8 CaCl2, and 5 HEPES. A glass coverslip was positioned over the agar to reduce brain pulsation.
Whole-cell patch-clamp recordings were obtained using a “blind” technique as described previously (Margrie et al., 2002). Recording pipettes (4–6 MΩ) were filled with an intracellular solution containing the following: 135 mm K-gluconate, 4 mm KCl, 10 mm HEPES, 10 mm Na2-phosphocreatine, 4 mm Mg-ATP, 0.3 mm Na-GTP, 200 μm Oregon Green BAPTA 488-1, 20 μm Alexa 594, 0.2% (w/v) biocytin, pH 7.2, 291–293 mOsm. Positive pressure (200–400 mbar) was applied to the pipette as it was inserted through the agar and the pial surface of the cortex. The positive pressure was reduced to 25–30 mbar when the tip was at approximately the upper limit of L2/3. The pipette was then advanced in 2 μm steps. Voltage pulses were applied to the pipette (10–20 mV; 30 ms; 10 Hz), and the current response monitored. Positive pressure was relieved when the series resistance of the electrode abruptly increased after a 2 μm step, indicating that the tip of the pipette may have been pushed against a neuronal plasma membrane. Gentle suction (up to 100 mbar) was applied where necessary to obtain a gigaohm seal. Initial access resistance in the whole-cell configuration was typically 30–60 MΩ. All recordings were obtained using an Axoclamp 2B amplifier (Molecular Devices, Union City, CA). L2/3 pyramidal neurons were identified as in our previous studies using two-photon microscopy and histological reconstructions (Waters et al., 2003; Waters and Helmchen, 2004). Somatic depths were 150–500 μm below the pia.
Sharp microelectrode recordings were obtained using 100–120 MΩ microelectrodes filled with 2.5 m potassium acetate. Microelectrodes were advanced in 2 μm steps. Typically, their resistance decreased during passage through the tissue to ∼70–90 MΩ before a recording was obtained. After penetration, the neurons were found to be depolarized (to around −40 mV) and current was injected to hold the membrane potential at around −70 mV. In the subsequent 5–10 min, the membrane potential stabilized and the holding current could be reduced (to zero). The neuron was allowed to recover its membrane potential before data collection; no holding current was used during data collection.
Input resistance measurements and data analysis
Brief pulse protocol.
Brief current pulses (−300 pA for 80 ms) were injected through the recording pipette at 200 ms intervals, producing hyperpolarizing voltage deflections of ∼5–10 mV. This allowed sufficient time for the membrane potential to relax to a steady-state value before measurement, both during and after the pulse. Access resistance was monitored every 2–10 s, rebalancing the bridge if necessary. Data were discarded if access resistance changed by >10 MΩ over a 1 min period or >20 MΩ over the entire data collection period (typically 5 min).
We sorted the resulting voltage responses according to whether they occurred during Up or Down states, manually defining Up and Down state thresholds ∼3–5 and 10–15 mV more depolarized than the Down state membrane potential, respectively (see Fig. 1C). We then calculated the mean membrane potential during two 10 ms time windows; one 1–11 ms before the start of the pulse, and the other 30–40 ms after termination of the pulse. Pulses were defined as occurring during Down states when the membrane potential was below the Down state threshold in both time windows (before and after the pulse). Similarly, pulses were defined as occurring during Up states if the membrane potential was above the Up state threshold in both windows. Pulses not meeting either of these criteria, such as those during which a state transition occurred, were sorted into a third group, termed “Discards.” Input resistance was derived for each group according to Ohm’s law, by dividing the mean voltage deflection during the second half of the current pulse by the current amplitude (−300 pA). Decay time constants were derived from the mean voltage responses by fitting a single-exponential curve to the decay of the voltage trace, starting 1 ms after the end of the current pulse.
Constant current injection protocol.
Constant holding currents of up to ±1.6 nA amplitude were injected at 200 pA intervals. The order in which we imposed different currents was semirandomized to avoid systematic artifacts resulting from changes in cellular properties during prolonged depolarization or hyperpolarization. We recorded the membrane potential for 1 min during the injection of each current, starting 10–30 s after the start of the current injection to ensure steady-state conditions. For each period, a membrane potential histogram was calculated as a frequency distribution of the membrane potential, using a bin width of 0.1 mV. The absolute frequency was divided by the total number of points to give the probability for each membrane potential bin. Cumulative frequency histograms were calculated as the integrals of the membrane potential histograms.
All data are presented as mean ± SEM. Except where otherwise stated, all statistical analyses were performed using the two-tailed Wilcoxon matched-pairs signed-rank test in GraphPad InStat, version 3.05 (GraphPad, San Diego, CA).
Reversal potential measurements
The potential at which Up states reversed was measured using the following intracellular solution: 135 mm Cs-gluconate, 4 mm CsCl, 10 mm HEPES, 10 mm Na2-phosphocreatine, 0.3 mm Na-GTP, 0.2 mm Oregon Green BAPTA 488-1, 0.02 mm Alexa 594, 0.2% (w/v) biocytin, 10 mm (2,6-dimethylphenyl)carbamoylmethyl-triethyl-azanium (QX-314). Recordings at various membrane potentials were made using discontinuous current clamp (sampling frequency, 3–4.5 kHz). Distinguishing between Up and Down states in the whole-cell recording became difficult where the Up state amplitude was small, around the Up state reversal potential. We therefore simultaneously monitored Up and Down states in the electrocorticogram (ECoG) as described previously (Waters and Helmchen, 2004) and used the ECoG signal to determine when transitions between Up and Down states occurred. We defined a transition between states as occurring each time the smoothed ECoG signal crossed a threshold equal to the mean of the ECoG signal. This objective approach accurately identified many (although not all) Up and Down state transitions in the whole-cell recordings. The Up and Down state membrane potentials were then calculated as the mean voltages during each state. The membrane potential difference between Up and Down states (ΔV) was then plotted against the mean Down state membrane potential and fit with a straight line. The reversal potential of the Up state was determined as the intersection of the fit with the abscissa.
The ratio of mean inhibition to mean excitation (〈gi〉/〈ge〉) relates to the reversal potential (Erev) as follows: 〈gi〉/〈ge〉 = (Erev − Ee)/(Ei − Erev), where Ee and Ei are the reversal potentials of excitatory and inhibitory synapses, respectively. Excitatory synaptic currents reverse at ∼0 mV. From our intracellular and extracellular solutions, we estimate (from the Nernst equation) Ei to be approximately −95 mV. Hence, a 〈gi〉/〈ge〉 ratio of 0.1 would produce an Up state reversal potential of −9 mV.
Recordings were obtained with a Cs+-based intracellular solution and QX-314 to maximize voltage control and suppress spiking. Nonetheless, broad spikes were sometimes observed near threshold. We discarded traces containing spikes, leading to the absence of V–I data points between approximately −30 and −5 mV.
Mathematical description of anomalous rectification
Anomalous rectification in cortical pyramidal neurons has been reported by many authors (Connors et al., 1982; Stafstrom et al., 1982, 1985; Spain et al., 1987; Sutor and Zieglgänsberger, 1987). We used a novel mathematical description to quantify rectification in pyramidal neurons. The V–I relationship of a passive membrane is a simple linear function ΔV = RN,0ΔI with a voltage-independent input resistance RN,0. Here, we included rectification by introducing an additional quadratic term with coefficient cAR: This simple phenomenological description fit our experimentally determined V–I curves well. To incorporate this rectification model into the NEURON simulation environment (see below), we needed to derive an equation that describes the leak current as a function of voltage deviation from the resting membrane potential. This equation is the inverse function of Equation 1: with depolarizing voltage deflections producing positive (outward) currents. An analytical expression for the voltage dependence of input resistance can then be obtained by taking the derivative of Equation 1, which by definition equals the input resistance: Insertion of Equation 2 into Equation 3 results in a simple expression for the voltage dependence of input resistance, which allows us to calculate the expected input resistance at any membrane potential (e.g., the Up state level) as follows: Input resistance measurements using current pulses yield values slightly different from the value at the initial membrane potential, depending on the direction and amplitude of the resulting voltage deflections (often termed “apparent” input resistances). When comparing input resistances derived from current pulses and from the slopes of V–I relationships, we accounted for the current-induced voltage change by extracting expected input resistances (from V–I relationships) at the midpoint potential (i.e., the potential halfway between the resting potential and that obtained during the current pulse).
Numerical models
The goal of modeling was to obtain an estimate for the number of active synapses during spontaneous activity in vivo. To this end, we reconstructed Up states using simple numerical models constrained by our experimental results. Rather than explicitly including active conductances, we simply simulated a membrane with anomalous rectification as described above. Numerical simulations were performed in the NEURON simulation environment (NEURON 5.7 for Windows) (Hines and Carnevale, 1997). The standard passive leak current was replaced by an anomalously rectifying distributed membrane mechanism including the additional parameter for anomalous rectification cAR.
First, we constructed a point conductance model of an isopotential neuron, which was equipped with an anomalously rectifying membrane and which received both excitatory and inhibitory synaptic input. In accordance with our experimental results, input resistance was adjusted to 30 MΩ, membrane time constant to 9 ms, and the anomalous rectification parameter to cAR = 18 MΩ/nA. Excitatory synaptic conductance changes were simulated using a reversal potential of 0 mV, a peak amplitude of 1 nS, and a double-exponential time course with a rise time constant of 0.2 ms and a decay time constant of 1.7 ms (Häusser and Roth, 1997). Inhibitory inputs were simulated using a reversal potential of −75 mV, a 0.5 nS peak conductance change, and rise and decay time constants of 1 and 10 ms, respectively (Ali et al., 2001). The mean total synaptic conductance change was adjusted to 8 nS (±3%). Synaptic input was modeled using Poisson-distributed presynaptic spike trains. Note that for the point conductance model, many excitatory (or inhibitory) inputs are combined into one input. Because spike intervals were exponentially distributed, the mean spike rate of the merged input is equivalent to n-times the mean rate for n synapses, each activated at the same rate. To compare our results with high-conductance models of Up states, we also constructed Up states under the assumption of an 80% drop in input resistance without anomalous rectification. In this case, a >10-fold higher total mean synaptic conductance change of 138 nS (±3%) had to be used. Under these conditions, to keep the membrane potential depolarization within the experimental bounds, the ratio of inhibition to excitation had to be 3–4.
For a more realistic simulation, we used a compartmental model of a L2/3 pyramidal neuron based on a detailed morphological reconstruction kindly provided by D. Feldmeyer and A. Roth (Lübke et al., 2003). A specific membrane resistance of 4350 Ω cm2 and an axial resistance of 150 Ω cm were used. Specific membrane capacitance was 1 μF/cm2 in the soma and 2 μF/cm2 in dendrites (to account for dendritic spines). Anomalous rectification was adjusted to yield cAR = 18 MΩ/nA. A total of 10,000 excitatory and 2000 inhibitory synapses (DeFelipe and Fariñas, 1992) (with properties as in the point conductance model) was uniformly distributed throughout the dendrites. To determine input resistances in this model, we applied the same brief pulse protocol as in our experiments. In some experiments, additional synaptic input mimicking the afferent synapses from L4 spiny stellate cells to L2/3 neurons was added in the basal dendritic tree (25 synchronized synaptic inputs at a mean distance 68 μm from the soma) (Feldmeyer et al., 2002). The summation of synchronous input was also evaluated during both Down and Up states by increasing the number of synapses on the basal dendrites (maximum number of 125, with a mean distance from the soma of 77.6 μm).
Results
Input resistance increases during Up states in L2/3 neurons
Whole-cell recordings were obtained from L2/3 neocortical pyramidal neurons in young adult rats under urethane anesthesia. Somatic voltage recordings revealed a stereotypical pattern of spontaneous 10–20 mV membrane potential (Vm) fluctuations between a hyperpolarized Down state and a depolarized Up state (Fig. 1A). Up states mostly remained subthreshold with spontaneous APs occurring infrequently (0.089 ± 0.024 APs s−1; range, 0–0.4 APs s−1; n = 29 cells), consistent with previous reports (Margrie et al., 2002; Brecht et al., 2003).
To compare somatic input resistance (RN) during Up and Down states, we used three approaches. First, brief hyperpolarizing current pulses were delivered at regular intervals (Fig. 1B). During Down states, each current pulse produced a similar hyperpolarizing response. In contrast, current injections during Up states resulted in voltage responses with a wide range of amplitudes, some pulses producing larger and others smaller hyperpolarizations compared with the Down states (Fig. 1B). For each neuron, several hundred voltage responses were collected and then sorted into three groups, according to whether they occurred during Down states, Up states, or during state transitions (the latter termed the “Discards” group) (Fig. 1C). Surprisingly, the mean input resistance was higher during Up states than during Down states in every neuron examined (Fig. 1C,D) (29.2 ± 2.0, 31.6 ± 2.2, and 36.6 ± 2.3 MΩ for Down, Discards, and Up groups, respectively; n = 16 neurons; p < 0.0001 for Up vs Down).
Current-evoked voltage responses could be subject to errors caused by changes in access resistance or incorrect bridge balance. We therefore also used a second approach, analyzing the decay of the hyperpolarizing voltage deflections after the end of the current pulse. The time constant for the tail of this decay is insensitive to access resistance, but directly relates to input resistance (Rall, 1969; Koch, 1999). Single exponential curves were fit to the decay of the mean voltage deflections for all three groups (Fig. 1C). The time constant was longer during Up states than during Down states in all neurons examined (Fig. 1E) (mean values, 9.0 ± 0.6, 10.6 ± 0.6, and 12.8 ± 1.1 ms for Down, Discards, and Up groups, respectively; n = 16 neurons; values for Up and Down states were significantly different, p < 0.0001). These results are consistent with an increased mean input resistance during Up states.
A third method for measuring input resistance changes during spontaneous activity is to inject constant current and analyze the resulting steady-state membrane potential changes during Up and Down states. If input resistance differs between Up and Down states, the steady-state membrane potentials of the two states will change by different amounts (Paré et al., 1998). In L2/3 neurons, we observed a significant RN increase during Up states using this approach (Fig. 2) (37.9 ± 7.0 and 42.0 ± 8.2 MΩ for Down and Up states, respectively; p = 0.02). Hence, using three complementary techniques, we consistently found that RN increases by 10–30% during Up states. Similar results were obtained in three additional dendritic recordings and in two recordings of fast-spiking neurons. Data are summarized in Table 1.
Our finding that RN increases during Up states contrasts with several previous studies reporting a decrease of up to 80% in cortical pyramidal neurons (Paré et al., 1998; Steriade et al., 2001; Destexhe et al., 2003; Léger et al., 2005). To examine whether different techniques (whole-cell recording vs sharp microelectrode impalements) might have produced such different results, we repeated measurements with sharp microelectrodes. In four microelectrode recordings, we found no significant difference between Up and Down state input resistances using all three approaches described above (Table 1; supplemental Fig. 1, available at www.jneurosci.org as supplemental material). We conclude that recording technique alone cannot account for the difference between our results and those of previous studies.
Voltage dependence of Up and Down state durations
Up states are driven by volleys of synaptic activity (Cowan and Wilson, 1994; Petersen et al., 2003). How can RN increase during such active periods, which are associated with the opening of synaptic conductances?
A clue came from our constant current injection experiments, in which the proportion of time spent in Up and Down states changed with membrane potential. This voltage dependence was evident in the raw data (Fig. 3a) and in cumulative membrane potential frequency histograms, in which the point of inflection shifted toward lower values with depolarization (Fig. 3b). Hence, when depolarized, neurons spent more time in the Up state and, when hyperpolarized, more time in the Down state. This effect was observed in all five neurons examined. To quantify the voltage dependence, we plotted the ratio of time spent in Up and Down states as a function of the membrane potential shift caused by constant current injection. The slope of a linear fit to the pooled data indicated that the fraction of time spent in the Up state changed by 1.7 ± 0.9% per millivolt of depolarization or hyperpolarization (Fig. 3c).
A similar voltage dependence of Up and Down state durations has been reported for cortical neurons (Cowan and Wilson, 1994) and for neostriatal spiny projection neurons, in which it was attributed to the opening of intrinsic potassium conductances with hyperpolarization, which decreases input resistance. It has been suggested that this increases the number of excitatory synapses that must be activated to initiate an Up state and thereby reduces Up state duration (Wilson and Kawaguchi, 1996). Intrinsic conductances operating in the subthreshold voltage regimen may, likewise, influence Up and Down state durations and input resistance in L2/3 pyramidal neurons. We therefore more closely investigated the intrinsic subthreshold membrane properties in L2/3 pyramidal neurons.
Anomalous rectification in L2/3 pyramidal neurons
Input resistance is defined as the slope of the voltage–current (V–I) relationship. Hence any voltage-dependent changes will appear as rectification in the V–I relationship (Koch, 1999). To examine the role of voltage-dependent currents, we therefore obtained steady-state voltage-current (V–I) relationships from 10 neurons, all of which displayed anomalous (inward) rectification at subthreshold membrane potentials (i.e., increasing RN with depolarization) (Fig. 4A,B).
Unsurprisingly, input resistance dropped dramatically in all neurons when current injections were suprathreshold, causing APs (19.0 ± 4.9 MΩ; n = 10) (Fig. 4). This drop is attributable to the opening of voltage-gated channels during AP generation (Häusser et al., 2001) and probably explains the difference between our results and previous measurements, because the latter were often from neurons firing APs at high frequencies during Up states (5–30 Hz) (Steriade et al., 1993, 2001; Destexhe and Paré, 1999).
Subthreshold V–I curves were fit assuming a quadratic relationship between injected current and the resulting voltage deflection, ΔV = RN,0ΔI + cARΔI2, where RN,0 denotes input resistance at rest and cAR is the rectification coefficient (Fig. 4B). Mean values for RN,0 and cAR were 27.5 ± 2.1 MΩ and 18.7 ± 3.7 MΩ/nA, respectively (n = 10; cAR range, 4–40 MΩ/nA). From this description, we derived an analytical function for the voltage dependence of RN (see Materials and Methods), allowing us to calculate the input resistance expected at the Up state membrane potential based on the V–I relationship of the neuron (RUpVI). We then compared RUpVI to RN measured during Up states (Fig. 4C,D). For the 10 neurons analyzed in this way, RUpVI was 44.1 ± 4.0 MΩ and the Up state input resistance was 38.8 ± 2.7 MΩ. Thus, Up state input resistance was reduced by 12% compared with the input resistance without synaptic input at the same membrane potential. By eliminating the effect of membrane potential, this calculation reveals the true effect of synaptic activity on input conductance.
We also eliminated the effect of membrane potential by using constant current injection to depolarize Down states to the Up state membrane potential. This approach revealed an 8% reduction of RN during Up states (Fig. 5A,B) (35.4 ± 4.7 MΩ at −57.5 ± 2.6 mV during Up states vs 38.1 ± 5.1 MΩ at −57.6 ± 2.0 mV during “isopotential” Down states; n = 7). Finally, we blocked intrinsic voltage-dependent currents by adding 10 mm QX-314 to the pipette solution. This blocked anomalous rectification in the voltage range of Up and Down states, linearizing the V–I relationship (Fig. 5C) (RN,0 = 68.1 ± 4.8 MΩ; cAR = 0.6 ± 6.4, range, −17–27 MΩ/nA; cAR significantly less than without QX-314, p = 0.023, Mann–Whitney test). With QX-314, input resistance was 15% lower during Up than during Down states (Fig. 5D) (52.7 ± 4.8 MΩ during Down state; 44.6 ± 2.2 MΩ during Up state; n = 4). We conclude that prominent anomalous rectification in L2/3 neurons explains the increased RN during Up states. If the voltage dependence of RN is accounted for or removed (by comparing isopotential RN values or using QX-314), RN during Up states is 10–15% lower than during Down states.
Excitatory and inhibitory synaptic conductance changes during Up states
To quantify the mean synaptic conductance change visible at the soma (gsyn), we compared somatic input conductance (the reciprocal of RN) in the absence and presence of synaptic input using the same three approaches used above. Comparison of measured conductances changes (during Up states) with the expected conductance change in the absence of synaptic activity revealed a synaptic conductance change gsyn of 8.2 ± 2.0 nS (range, 0–17 nS; n = 10). Calculating gsyn from the difference of conductances in Up states and isopotential Down states and from the conductance difference with anomalous rectification blocked by QX-314 yielded comparable values of 2.2 ± 1.4 (n = 7) and 3.0 ± 1.0 nS (n = 4), respectively. Thus, only a small conductance change in the range of 2–10 nS is visible at the soma during cortical Up states.
We used this conductance change to estimate the number of active synapses and presynaptic firing rates by reconstructing cortical Up states in numerical simulations. We started with a point conductance model, in which resting membrane properties and anomalous rectification were matched to experimental results (Fig. 6A). Excitatory and inhibitory synaptic inputs were modeled as Poisson-distributed presynaptic spike trains evoking AMPA- and GABAA-type postsynaptic currents, respectively (see Materials and Methods). Assuming a total conductance change of 8 nS we first explored the consequences of changing the mean inhibition/excitation ratio (〈gi〉/〈ge〉). The model reproduced Up states with mean depolarizations and input resistance increases consistent with our experimental results if 〈gi〉/〈ge〉 was <0.3 (Fig. 6B,C). A ratio of 0.1 gave the closest match, resulting in 15 mV mean depolarization and 12.5% net input resistance increase during Up states (Fig. 6B,C).
In contrast, if we excluded anomalous rectification from our model and assumed a large decrease in input resistance during Up states, a 〈gi〉/〈ge〉 ratio of ∼3 was required to produce 15 mV Up states (supplemental Fig. 2, available at www.jneurosci.org as supplemental material). This value is comparable with previous estimates, which have suggested that inhibition dominates excitation (Destexhe et al., 2003).
The ratio of inhibition to excitation also determines the reversal potential of Up states (Wilson and Kawaguchi, 1996; Shu et al., 2003). Under our conditions, a ratio of 0.1 would produce an Up state reversal potential of −9 mV (see Materials and Methods). We directly measured the reversal potential of Up states in vivo by depolarizing QX-314-filled neurons to different membrane potentials and measuring the mean voltage difference between Up and Down states. Up states were defined using a simultaneously measured electrocorticogram (Waters and Helmchen, 2004). The mean reversal potential was −10.4 mV (Fig. 6D) (n = 5 neurons), which corresponds to a mean conductance ratio (〈gi〉/〈ge〉) of 0.12. Thus, both simulations and experiments suggest that inhibition contributes ∼10% to the total conductance change.
Estimate of numbers of activated synapses from a compartmental model
Next, we used a compartmental model of a morphologically reconstructed L2/3 pyramidal neuron to convert the conductance change into an estimate of the number of active synapses (Fig. 7A,B). The neuron was equipped with anomalous rectification so that it reproduced the experimentally measured V–I curves. A total of 10,000 excitatory and 2000 inhibitory synapses was uniformly distributed throughout the dendritic tree (DeFelipe and Farinas, 1992). As in the point conductance model, Up states with a mean depolarization of 15 mV and increased somatic input resistance were generated with low rates of synaptic input (0.33 and 0.067 Hz per input for excitatory and inhibitory synapses, respectively). This corresponds to ∼3 excitatory and 1.3 inhibitory synapses activated per millisecond.
Interestingly, our simple model captured several experimentally observed features of L2/3 neurons. First, although voltage fluctuations during Up states increased with distance from the soma, the mean depolarization was nearly constant throughout the apical dendritic tree (Fig. 7C), similar to our in vivo whole-cell recordings from soma, apical dendrite, and distal tuft branches (Waters and Helmchen, 2004). Second, several studies reported that sensory-evoked EPSPs show reduced amplitudes during Up states (Petersen et al., 2003; Sachdev et al., 2004). Simulating this situation in our model (see Materials and Methods), we found that the compound EPSP amplitude during Up states was reduced to 86% of the Down state amplitude despite a mean input resistance increase of 27% (Fig. 7D). In our model, this attenuation of the compound EPSP is caused exclusively by the reduction in driving force. The larger attenuation observed for sensory-evoked responses (Petersen et al., 2003; Sachdev et al., 2004) presumably results from additional factors such as synaptic depression and feedforward and feedback inhibition.
Finally, we used our compartmental model to estimate how many synapses are required to depolarize the neuron to different membrane potentials. We considered two situations: asynchronous, distributed background activity and synchronous synaptic inputs onto the basal arbors (such as occurs during activation of a specific afferent pathway). The excitation/inhibition ratio was kept constant at 0.1. With asynchronous synaptic input and without anomalous rectification, the mean depolarization per synapse decreased as the frequency of inputs was increased (Fig. 8A). This decrease is attributable to the diminished driving force for excitatory currents at depolarized membrane potentials. In the presence of anomalous rectification, synaptic inputs summated more effectively but a small decrease in mean depolarization per synapse persisted despite the increase in somatic input resistance. Thus, anomalous rectification did not fully linearize the relationship between net depolarization and frequency of synaptic inputs. Nevertheless, only ∼14 active synapses per millisecond were required to reach AP threshold. In the case of synchronous excitatory inputs to the basal dendrites, more than one hundred synapses were required to bring the neuron to AP threshold during Up states (Fig. 8B). We conclude that, during background network activity, spiking is typically driven by tens to little more than a hundred synapses, depending on their rate of activation and their synchrony.
Discussion
We have shown that the input resistance of L2/3 neurons increases during cortical Up states and that this is because anomalous rectification dominates the synaptic conductance increase. We conclude that, first, Up states are not high-conductance states in L2/3 pyramidal neurons but are associated with stable or reduced input conductance. Second, anomalous rectification balances the synaptic conductance increase in vivo, thereby reducing the shunting problem. Finally, few synapses are active during Up states (<5 per millisecond) and the majority (∼90%) are excitatory.
Input resistance during Up states
Using three complementary methods, we found a 10–30% mean increase in input resistance during Up states, equivalent to a 10–25% drop in input conductance. This resembles in vivo results from neostriatal spiny projection neurons (Wilson and Groves, 1981; Wilson and Kawaguchi, 1996) but contrasts with previous studies of cortical neurons that reported twofold to fivefold decreases in input resistance during Up states (Paré et al., 1998; Steriade et al., 2001; Léger et al., 2005) (for review, see Destexhe et al., 2003) (but see Cowan and Wilson, 1994). This discrepancy is not related to the recording technique, because we obtained similar results with whole-cell patch pipettes and sharp microelectrodes. Age and species differences are also unlikely to be responsible, because the conductance change necessary to produce a 50–80% decrease in input resistance would be approximately an order of magnitude larger than we measured here [compare Fig. 6 and supplemental Fig. 2 (available at www.jneurosci.org as supplemental material)], and it seems unlikely that total synaptic input would differ sufficiently between species or between young adult and adult animals.
The most likely explanation is that previously reported decreases in input resistance resulted principally from the opening of voltage-gated channels during APs, not from synaptic conductances. Consistent with this explanation, AP firing caused an apparent input resistance decrease of 50% or more in our recordings (Fig. 4). In agreement with recent in vivo studies (Margrie et al., 2002; Brecht et al., 2003), we observed few APs during Up states. In contrast, other authors have reported high-frequency (tens of hertz) firing during Up states (Steriade et al., 2001). Why different authors report such disparate firing rates during Up states is unclear. Leak around the pipette may be a contributory factor in sharp microelectrode penetrations, although in our microelectrode recordings firing rates were similar to whole-cell recordings (∼0.1 Hz). Low firing rates have been reported in cell-attached recordings in vivo (Margrie et al., 2002), demonstrating that dialysis with patch pipettes is not a critical factor. Furthermore, neocortical AP activity is sparse in unanesthetized animals (Margrie et al., 2002; Petersen et al., 2003). Hence input resistances may be determined principally by subthreshold intrinsic currents and synaptic activity in awake animals.
Anomalous rectification
Anomalous rectification occurs in many cell types in the CNS, including spinal cord motoneurons (Nelson and Frank, 1967), hippocampal pyramidal (Hotson et al., 1979), cerebellar Purkinje (Llinás and Sugimori, 1980), locus ceruleus (Osmanovic and Shefner, 1987), neostriatal spiny (Wilson, 1992), and cortical pyramidal neurons (Connors et al., 1982; Stafstrom et al., 1982; Sutor and Zieglgänsberger, 1987; Cowan and Wilson, 1994). Hence, anomalous rectification counterbalancing synaptic conductance changes is likely to be widespread in the CNS. Although anomalously rectifying currents have been identified as a key factor influencing Up states in neostriatal spiny neurons (Wilson 1992) and despite descriptions of anomalous rectification in the cortical literature, the role of anomalous rectification during Up states in cortical pyramidal neurons has been mostly overlooked.
Several voltage-dependent ion channels may contribute to anomalous rectification. Inwardly rectifying K+ currents are typically responsible at hyperpolarized voltages (Osmanovic and Shefner, 1987; Sutor and Zieglgänsberger, 1987; Nisenbaum and Wilson, 1995). In the depolarizing direction, inward sodium and calcium currents have also been implicated (Stafstrom et al., 1985; Sutor and Zieglgänsberger, 1987). Here, we have shown that QX-314 blocked anomalous rectification in the voltage range of Up and Down states, consistent with previous studies of cortical pyramidal neurons (Stafstrom et al., 1982, 1985; Sutor and Zieglgänsberger, 1987).
Synaptic conductance change during Up states
During Up states, the mean synaptic conductance change visible at the soma is 8 nS. This is >10-fold smaller than estimates from several in vivo studies (Destexhe et al., 2003), but comparable with recent in vitro and in vivo estimates from ferret cortex (4–32 nS) (McCormick et al., 2003; Shu et al., 2003; Haider et al., 2006). In contrast to these ferret studies, which suggest balanced inhibition and excitation, our simulations and reversal potential measurements point to a lower inhibition/excitation ratio of ∼0.1, indicating that inhibition is sparse during cortical Up states. Reversal potential measurements are, however, prone to errors, including junction potentials and misestimation of the reversal potential of IPSPs (Ei). In 3- to 5-week-old rats, Ei has been measured at around −70 mV (Luhmann and Prince, 1991; Gulledge and Stuart, 2003). To estimate the possible error in our calculations, we recalculated 〈gi〉/〈ge〉 with Ei = −70 mV (instead of −95 mV) and assuming a 10 mV junction potential, yielding 〈gi〉/〈ge〉 = 0.41, approximately fourfold higher than our estimate of 〈gi〉/〈ge〉 = 0.1. Our conclusion, that the majority of synaptic drive during Up states is excitatory, not inhibitory, is therefore robust even against sizable possible errors in our calculations.
Our simulations indicate that postsynaptic activity is sparse during Up states (∼0.3 Hz for excitatory synapses). To convert this number into average firing rates of presynaptic neurons, one needs to divide it by the presynaptic release probability, which is unknown for many cortical connections. Assuming a release probability of 0.8, as estimated for connections between L5 pyramidal neurons and from L4 spiny stellates to L2/3 pyramidal neurons (Markram et al., 1997; Silver et al., 2003), would give a presynaptic firing rate of ∼0.4 Hz during Up states. Neurons spent approximately equal amounts of time in Up and in Down states (Table 1), so presynaptic firing rates averaged over time (including both states) would be 0.2 Hz. This is consistent with the low spontaneous firing rates observed in our whole-cell recordings and in previous studies (Margrie et al., 2002; Brecht et al., 2003; Kerr et al., 2005).
Our recordings were obtained from somata, far from most active synapses. Presumably, conductance changes in some locations in the dendrites will be greater. Distance-dependent attenuation of conductance along dendrites is severe (Koch et al., 1990; Williams, 2004). Hence, a drop in input resistance might occur locally in a dendrite, causing local shunting but little effect at the soma (Rall, 1967; Koch et al., 1990). In addition, activation of different combinations of synapses during different Up states (either close to the soma or more distally) might produce different conductance changes at the soma. This could explain the variability that we observed in the amplitudes of voltage deflections between individual Up states (Fig. 1B).
Our simple model reproduced the experimental results assuming a uniform distribution of synaptic input throughout the dendrites. This assumption is probably reasonable for excitatory but not for inhibitory synapses (DeFelipe and Fariñas, 1992; Markram et al., 2004; Waters and Helmchen, 2004). Future models might consider the effect of the location of inhibitory synapses in the dendritic tree. Furthermore, our model needs to be extended to the suprathreshold regime, because local dendritic spikes generated by synchronous, clustered inputs (Gasparini and Magee, 2006) might play an important role in synaptic integration.
Functional implications
Sparse background activity in cortical networks may enable the cortex to use computationally and energetically efficient sparse coding schemes despite ongoing network activity (Olshausen and Field, 2004). In particular, the influence of synaptic background on single-cell computation may be less than previously thought. With few active synapses and anomalously rectifying currents, dendrites become electronically more compact during synaptic background activity (the opposite is expected for a linear membrane) (Koch, 1999). This may explain why AP backpropagation is conserved in vivo, even during synaptic activity (Buzsáki and Kandel, 1998; Charpak et al., 2001; Quirk et al., 2001; Waters and Helmchen, 2004; Waters et al., 2005). Likewise, both forward propagation of dendritically initiated spikes (Golding and Spruston, 1998; Larkum et al., 2001) and amplification of EPSPs by voltage-gated channels (Stuart and Sakmann, 1995) should be preserved during network activity in vivo.
Interestingly, in L5 pyramidal neurons, anomalous rectification becomes stronger during postnatal development, compensating for the decrease in input resistance that occurs in parallel (Kasper et al., 1994). Thus, anomalous rectification could tune network activity during development. Furthermore, neuromodulatory regulation of voltage-dependent currents operating in the subthreshold range could affect the emergence of different network states.
In summary, we have shown that surprisingly few synapses drive spontaneous depolarizations in L2/3 pyramidal neurons in the anesthetized rat. Their contribution to membrane conductance is opposed by anomalous rectification, which enhances the effectiveness of each synapse during volleys of synaptic activity and decreases the influence of synaptic activity on dendritic properties. Against such a low synaptic background, a few precisely timed synaptic inputs (on the order of 100) should be sufficient to produce and maintain specific activation patterns in cortical circuits.
Footnotes
-
We thank Bert Sakmann for his encouragement and generous support. We are especially grateful to Arnd Roth for discussions and help with simulations. The L2/3 pyramidal cell morphology was kindly provided by Dirk Feldmeyer. We thank Jason Kerr, Arnd Roth, and Bert Sakmann for critically reading this manuscript and Marlies Kaiser for expert technical assistance.
-
↵F. Helmchen’s present address: Department of Neurophysiology, Brain Research Institute, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland.
- Correspondence should be addressed to Jack Waters, Department of Physiology, Feinberg School of Medicine, Northwestern University, 303 East Chicago Avenue, Chicago, IL 60611. Email: jackwaters{at}northwestern.edu