|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1Computation and Neural Systems Program, California Institute of Technology, Pasadena, California; 2Research Imaging Center, University of Texas Health Science Center at San Antonio; and 3Cajal Neuroscience Research Center at the University of Texas San Antonio, San Antonio, Texas
Submitted 12 October 2006; accepted in final form 13 October 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The lack of clear pf beamlike effects in vivo from an input that, in mammals, provides approximately 150,000 excitatory synapses per PC (Gundappa-Sulur et al. 1999
; Harvey and Napper 1991
) has posed a physiological conundrum for almost 50 years (for review see Bower 2002
). All published explanations for the lack of beamlike effects under more natural stimulus conditions have proposed that pfs are less powerful than previously believed (Eccles et al. 1972
), either because of the desynchronization of action potentials as they travel along the pfs (Llinás 1982
) or because of the small number of simultaneously activated granule cells (Braitenberg et al. 1997
).
We first used a computer model of the cerebellar cortex to test these previous explanations. We found that even when the maximum likely pf desynchronization was combined with a very small number of activated granule cells, the model produced a full beam of activated PCs. However, when molecular layer inhibition was added to an otherwise identical network model, beamlike pf effects did not occur, suggesting that molecular layer inhibition is responsible for the lack of beamlike activation of PCs under natural stimulus conditions in vivo. To test this most fundamental prediction of the model, we recorded tactile evoked responses in PCs in vivo located along the pfs in the presence and absence of two different
-aminobutyric acid type A (GABAA) receptor blockers (bicuculline and gabazine). As the model predicted, blocking molecular layer inhibition with either drug resulted in the emergence of beamlike patterns of PC activity consistent with stimulation by pf synapses.
After experimentally testing the model prediction, we completed a more detailed analysis of our modeling results, producing several more specific predictions for the inhibitory mechanisms underlying the ability of feedforward inhibition to counterbalance pf excitation. These results will serve as the basis for future experimental and modeling studies and also have important implications for theories of cerebellar function in general and the role of pfs in particular.
| METHODS |
|---|
|
|
|---|
MODEL STRUCTURE.
Purkinje cell.
We used our previously published PC model (De Schutter and Bower 1994a
,b
) with updated synaptic kinetics. The ion channels, their kinetics, and distributions are described in detail in the Supplementary Material Table S1.
Granule cells.
Granule cells were distributed in a 300 x 2,500-µm area, 20 µm below the PC somas. Each cell gave rise to an axon that coursed vertically into the molecular layer (ascending segment), making a synaptic contact with the overlying PC every 1090 µm (Gundappa-Sulur et al. 1999
; Harvey and Napper 1991
). After ascending for varying distances through the molecular layer, the granule cell axon bifurcates into a pf (Fig. 1 A).
|
Molecular layer interneurons.
The biophysical properties of these cells were not represented explicitly, but instead their influence was parameterized by the activation of their synapses on PCs. Molecular interneurons make two types of synapses on PCs: stellate-type connections on dendrites and basket-type connections on the PC soma (Sultan and Bower 1988
). The type and number of these connections are dependent on the depth of the interneuron's soma in the molecular layer: the deeper the soma the more basket-type synapses it produces (Sultan and Bower 1998
). Reflecting these anatomical constraints, the PC's dendritic tree was divided into three horizontal overlapping layers: from 100 to 200, 100 to 300, and 200 to 400 µm, where 0 is the position of the soma (Fig. 1B). Synapses in the top layer made only stellate-type connections, whereas cells in the middle and bottom layers made both stellate- and basket-type connections, with the bottom layer biased toward basket-type synapses. Stellate- and basket-type synaptic kinetics are described in Supplementary Material Table S2.
MODEL ACTIVATION.
Background excitatory and inhibitory activity.
We implemented 1,600 granule cells (and the same number of passive spines on the PC), which corresponds to roughly 1% of the total granule cell input believed to converge onto an average PC (Harvey and Napper 1991
). We compensated for this reduction in the number of inputs by increasing the firing rate of these synapses (De Schutter and Bower 1994c
) with the updated synaptic kinetics (Supplementary Material Figure S1). The number of modeled inhibitory inputs contacting the PC dendrite in these simulations was 1,695, which anatomical data suggest is close to a realistic value (Sultan and Bower 1998
).
Unless otherwise noted the background activation of the granule cells and inhibitory synapses was modeled as a Poisson process with a mean firing rate of 4.0 Hz for granule cells and 1.1 Hz for inhibitory synapses. This background stimulation reproduced the average 40-Hz PC firing rate found in vivo (Bower and Woolston 1983
). Because the random background pf and inhibitory inputs were assumed to arise from spontaneous granule cell layer activity from throughout the cerebellum folium, we did not correlate these inputs in the model. As described later, however, stimulus-evoked activity in the pfs and molecular layer interneuron was correlated through the appropriate network connections.
Focal granule cell layer activation.
Reflecting the patchy organization of mossy fiber tactile inputs (Shambes et al. 1978
) as well as the width of our PC model dendrite, a focal activation of the granule cell layer was modeled by synchronously activating a nearly 300-µm-wide section of the granule cell layer with variable length. The ascending and pf synapses associated with the granule cells in this patch were activated after the temporal delay resulting from the propagation velocity of the axon and the distance from the site of activation. We opted for a synchronous activation of all granule cells instead of a temporally spread function to rigorously test the desynchronization hypothesis. Any temporal spread in activation of the granule cells would increase the level of desynchronization along the pfs, resulting in a relaxing of the conditions for setting up compensatory feedforward inhibitory inputs on PCs (see RESULTS).
An upper value for the maximum number of activated granule cells was estimated based on data obtained from micromapping tactile responses in the cerebellar granule cell layer (Bower and Kassel 1990
; Shambes et al. 1978
). Specifically, a single upper lip tactile stimulus to the ipsilateral rat face consistently activates a 500 x 500-µm cortical area in the center of folium Crus IIa (Fig. 4). Using published values for the average concentration of granule cells (Harvey and Napper 1991
) and taking into account the 300-µm width of the modeled PC dendritic tree, we estimated that an upper lip tactile stimulus results in at most 30,000 active granule cell synapses per PC. For the pfs, this value corresponds to 20% of the total 150,000 pf excitatory inputs per PC.
|
A summary of all network parameters is shown in Supplementary Material Table S2.
DATA ANALYSIS.
Peristimulus time histograms (PSTHs, 1-ms bin) were constructed from PC spike times 10 ms before and 20 ms after simulated stimulus onset. The significance of changes in the shape of the PSTH, poststimulus, was quantified using t-test, dividing the total number of trials (128 or 256) in subsets of 16 (Santamaria and Bower 2005
; Santamaria et al. 2002
) and comparing the firing rates of the PC before and after the stimulus. The pattern of background excitatory and inhibitory synaptic inputs was randomized on each trial and was therefore completely independent of the pattern of stimulus-evoked excitatory and inhibitory input.
Simulations for each set of data were performed using a fixed range for pf propagation velocities and fixed numbers of activated granule cells. Systematic parameter variations were made until the model produced PSTHs that came closest to matching previous experimental results.
To avoid any biases introduced by the particular spatial distribution of inputs, the physical positions of all the synaptic inputs on PC dendrites were randomized from trial to trial. All simulations were implemented in GENESIS 2.1 (Bower and Beeman 1995
) running on supercomputers operated by the San Diego Super Computer Center at the University of California, San Diego.
Experimental procedures
Modeling predictions were tested, in vivo, using fourteen 3- to 6-mo-old female SpragueDawley rats. Animal handling, surgical procedures, and euthanasia techniques were approved by the Caltech Animal Care and Use Committee as well as the Animal Use Committee at the University of Texas Health Science Center in San Antonio as in compliance with the National Institutes of Health guidelines. Details of surgical procedures and recording techniques can be found in the Supplementary Materials. All recordings were obtained in animals anesthetized using ketaminexylazineacepromazine (ketamine 100 mg/kg; xylazine 5 mg/kg; acepromazine 1 mg/kg). This combination of drugs results in granule cell layer responses most similar to those obtained in awake behaving animals (cf. Fig. 2, trace 2; Hartmann and Bower 2001
).
|
for PCs and 2 to 4 M
for recording of the granule cell layer. PCs responses were recorded before and after topical application of either GABAA receptor blocker bicuculline or gabazine (SR95531), all obtained from Sigma. These drugs were applied using a micropositioner to place a 10-µl microsyringe (Hamilton) within several microns of the surface of the cerebellum. Between 1 and 7 µl of either a 5 mM saline solution of bicuculline or 60 µM gabazine were applied to the surface of the cerebellum, forming a small bubble that spread over the surface, covering the position of all recording electrodes. Data collection was resumed 1015 min after drug application. Multiunit recordings of granule cell layer activity were also recorded before and after inhibitory blockers were applied. Data analyses are reported as SE unless otherwise noted. Changes in PC activity recorded experimentally were quantified using the same statistical procedure as that for the modeled data (see above).
| RESULTS |
|---|
|
|
|---|
Our initial objective was to use the computer model to test the hypothesis that either the desynchronization of action potentials traveling along pfs (Llinás 1982
) or an insufficient number of simultaneously activated granule cells (Braitenberg et al. 1997
) is responsible for the lack of PC beams seen in vivo under natural stimulus conditions (Bower and Woolston 1983
; Kolb et al. 1997
).
To determine the largest likely level of pf desynchronization, we first used the model to estimate the spatial and temporal spread of action potentials along pfs after a focal activation of the granule cell layer. Figure 2A shows the spatial distribution of action potentials in the molecular layer calculated using the widest difference in conduction velocities reported in mammals between superficial (0.15 m/s in rats; Vranesic et al. 1994
) and deep pfs (0.5 m/s in cats; Crepel et al. 1981
). Figure 2A shows that with these values, an originally synchronous volley of action potentials would spread over >1.5 mm of the molecular layer by the time the fastest action potentials reached the end of the pfs (2.5 mm; Harvey and Napper 1991
). Note that we chose the widest reported mammalian values to provide the best possible conditions for the desynchronization hypothesis. A narrowing or lack of difference in superficial and deep pf conduction velocities would reduce the level of action potential desynchronization and generate more stringent conditions for this hypothesis. Also, whereas this graph was made assuming depth-related differences in pf conduction velocities, the lag in propagation of pf action potentials in the superficial layers of the molecular layer also occurs as a consequence of action potentials propagating along the full length of the ascending granule cell axon segment and the spatial distribution of the granule cells.
Figure 2B shows the distribution of the predicted propagation times for pf action potentials converging on PCs at different distances from a focal site of granule cell layer activation. The shortest propagation time (lower bound) is a combination of the fastest propagation velocity with travel along the deepest pfs. Similarly, the upper bound is the propagation time along the most superficial pfs. Because the model assumes that pf propagation velocities vary linearly with respect to their laminar position in the molecular layer the difference between the upper and lower bounds indicates the maximum synaptic desynchronization for PCs located at different distances from a focal site of granule cell layer activation. Accordingly, for a PC 1 mm away from the site of stimulation, the maximum synaptic desynchronization is 4 ms whereas for a PC at 2.5 mm, it is 13 ms.
Simulated PC responses to granule cell activity
The top two rows of PSTHs in Fig. 3 compare the typical spatial distribution of PC responses along the path of pfs recorded in vivo after a focal activation of the granule cell layer (Fig. 3, AD) with responses obtained in the model under maximum conditions of pf desynchronization (Fig. 3, EH). The in vivo data shown in Fig. 3, AD were obtained from four PCs recorded simultaneously at the distances shown from the site of focal granule cell layer activation. The PCs were recorded in a mediolateral line down the center of the crown of Crus IIa and therefore along the course of the pfs. In agreement with previous reports (Bower and Woolston 1983
; Kolb et al. 1997
), the only PC that responds with short-latency excitation is located immediately above the activated region of the granule cell layer (at 0 µm). PCs farther away from the site of granule cell layer activation, and along the course of the pfs, respond at a short distance with a reduction in spiking frequencies and at greater distances with no significant change in firing rate.
|
The generation of propagating beams of PC activity in the model was robust to a wide range of parameter values. For example, beams of PCs were found even when as little as 2% of the simulated granule cells were activated. Further, the results were independent of parameters related to the background levels of excitatory and inhibitory inputs (Supplementary Material Figure S1). Thus our modeling results suggest that under a broad range of parameters, even a highly desynchronized parallel fiber volley with a small number of action potentials is more than capable of inducing beamlike activation of PCs. Clearly, however, the simulations do not match the in vivo experimental results, strongly suggesting that some factor other than pf action potential desynchronization or the number of activated granule cells is responsible for the lack of PC beams in vivo.
Simulated effects of feedforward inhibition on PC responses to granule cell activity
The only stimulus-evoked synaptic inputs to PCs in the modeling results presented to this point were provided by excitatory granule cell synapses. However, pfs also synapse directly on molecular layer interneurons, which themselves provide an inhibitory feedforward inhibition on PCs (Cajal 1904
; Eccles et al. 1966a
). The PSTHs shown in Fig. 3, IL were obtained from simulations identical to those used for Fig. 3, EH, after adding feedforward molecular layer inhibition of both the basket- and stellate-types to the model. By tuning values for the delay only between pf excitation to PC inhibition, the total number of activated inhibitory synapses, and their spatial distribution, the model readily reproduced the short-latency excitation directly above the activated granule cell layer (Fig. 3I), the reduction in firing frequency at slightly longer latencies in overlying and nearby PCs (Fig. 3, I and J), and the lack of PC responses at a distance along the pfs (Fig. 3, K and L), all also seen in the in vivo data (Fig. 3, AD). Similar results were obtained over a wide and robust range of model parameters (see Fig. 6), including a much narrower range of pf propagation velocities (0.200.27 m/s; Supplementary Materials Figures S1 and S2). As described in the final section of the RESULTS, all values for molecular layer inhibition were also within likely physiological values. The model therefore predicts that the lack of PC beams in vivo, after focal activation of the granule cell layer, is a consequence of the presence of feedforward molecular layer inhibition.
|
To test the prediction that the lack of pf activated beams of PCs is attributed to the influence of molecular layer inhibitory interneurons we extracellularly recorded granule and PC activity before and after topical application of two different GABAA receptor blockers: bicuculline (Fig. 4) and gabazine (Fig. 5). As in previous studies (Bower and Woolston 1983
; Lu et al. 2005
), we took advantage of the known fractured somatotopy of afferent tactile projections to the cerebellar folium Crus IIa, in which small areas of the granule cell layer receive tactile input from distinctly different perioral regions (Bower and Kassel 1990
).
|
During the same experiment we recorded the activity of a PC located 1,400 µm lateral to the site of maximum recorded granule cell layer activation. Figure 4B shows the calculated PSTHs (240 trials) for this cell during three different experimental conditions. The top histogram shows a recording obtained before the application of bicuculline, in which there is no statistically significant response to CUL stimulation. This lack of response is expected because this cell overlies a region of the granule cell layer that does not itself respond to stimulation of the CUL (Fig. 3, C and D; see Bower and Woolston 1983
). In contrast, the middle histogram obtained 10 min after bicuculline application shows a strong excitatory response with an onset at 17 ms. If this postbicuculline CUL response arises from pfs, we would expect a longer-latency PC response to stimulation of the CUL than the IUL represented in the immediately underlying granule cell layer. The bottom histogram in Fig. 4B confirms this prediction, showing that IUL stimulation results in a peak excitatory response at 12 ms. The 5-ms difference in the delay of onset between the two responses predicts a putative pf propagation velocity of 0.28 m/s, which is within the reported physiological range (Vranesic et al. 1994
).
Figure 5 shows the results of a similar experiment using the more specific GABAA receptor blocker gabazine (60 µM). Figure 5A shows the diagram of Crus IIa for this experiment with overlaid averaged granule cell layer responses to CUL stimulation. The responses were recorded after 5 µl of gabazine was applied. As in Fig. 4A, the response is restricted to the medial regions of the folium even after application of gabazine. Figure 5B shows PSTHs for two PCs recorded over the IUL patch before gabazine application (PC1 and PC2). Because these PCs overlie a region of the granule cell layer activated by the IUL, CUL does not evoke an excitatory response in these cells under control conditions. Figure 5C shows three PSTHs obtained from another PC recorded 1,200 µm away from the site of stimulation (PC3) to CUL stimulation and IUL before the application of gabazine (top and middle histograms) and to CUL stimulation after application of the drug (bottom histogram). As in the case of bicuculline, this PC responds to CUL stimulation only after GABAA inhibition is blocked (bottom histogram). Based on the difference in the onset of this PC response between IUL and postgabazine CUL stimulation, the calculated pf propagation velocity was 0.12 m/s. In those experiments in which putative pf propagation velocities were calculated, the average value with bicuculline was 0.27 ± 0.02 m/s (n = 3), whereas it was 0.21 ± 0.07 m/s with gabazine (n = 3).
Figure 6 shows the results from a different experimental procedure specifically intended to both replicate the effects of blocking GABAA inhibition and also estimate pf conduction velocities. We simultaneously recorded from two PCs separated either by 380 µm (Fig. 6A) or 760 µm (Fig. 6B) along the course of the pfs before and after bicuculline application. The top two PSTHs in Fig. 6A were recorded from a cell overlying a region of the granule cell layer that was activated by tactile stimulation of the furry bucal pad (FBP). This cell responds with excitation both before (left PSTH) and after (right PSTH) topical application of bicuculline. The second PC (bottom) was close to but not overlapping the FBP granule cell patch and therefore produced only a short-latency inhibitory response to the same stimulus under control conditions. However, after bicuculline was applied, this cell responded with short-latency excitation. The delay between the onsets of excitation between these two cells is 1 ms, which results in a calculated pf conduction velocity of 0.38 m/s. For all seven pairs of cells recorded with a 380-µm separation, the average calculated pf conduction velocity was 0.26 ± 0.05 m/s.
Similarly, Fig. 6B shows data from two PCs separated by 760 µm after IUL stimulation. The cell located over the IUL granule cell patch again responded with increased firing both before (left) and after (right) bicuculline application. The second cell (bottom) showed no change in spiking under control conditions, but showed an increase after the application of bicuculline. In agreement with a delay resulting from action potentials propagating along the pfs, the difference in time to peak of the PCs separated by 760 µm was greater than that for the PCs separated by 360 µm. The calculated pf conduction velocity was 0.15 m/s. The average conduction velocity for the four pairs of cells recorded at a 760 µm distance was of 0.25 ± 0.08 m/s. The similarity in calculated pf conduction velocities in all experiments, using either drug, strongly suggests that the excitatory responses that emerge after blocking molecular layer inhibition are a direct result of pf input.
Although conduction along the pfs is the most likely explanation for the emergence of propagating short-latency excitatory PC responses after application of the GABAA-receptorblocking drugs, we also specifically looked for evidence for any changes in PC excitability. Neither bicuculline nor gabazine resulted in a statistically significant change in basal PC firing rates (Fig. 7, A and B). For bicuculline, average PC firing rates actually showed a slight, although nonsignificant decrease from 41.6 ± 4.3 to 33.8 ± 3.2 Hz (n = 30, eight experiments, t-test). In fact, both increases and decreases in firing frequencies were found for different cells in the same experiment. Similarly, gabazine resulted in a slight average reduction from 40.0 ± 7.8 to 39.2 ± 12.0 Hz (n = 5, four experiments, t-test, nonsignificant). Thus we found that blocking GABAA receptors resulted in no statistically significant change in PC excitability.
|
In summary, neither bicuculline nor gabazine resulted in any significant changes in background PC firing rates, or in the latency, amplitude, or spatial activation of granule cell layer responses to tactile stimulation. Thus the excitatory PC responses at a distance from focal activation of the granule cell layer after blocking inhibition cannot be attributed to a generalized increase in neuronal excitability, but instead is most likely the result of unmasked pf excitation.
Network effects of pf excitation and molecular layer inhibition on PCs
Having experimentally tested the most fundamental prediction of the model, we then explored in more detail the network (Fig. 8) and biophysical (Figs. 9 and 10) interactions underlying the model's ability to replicate in vivo results. Figure 8A shows, once again, the range of pf propagation times from the site of origin of granule cell layer activation to a distance of 2 mm. Figure 8B shows PC firing rates within 15 ms after the first arrival of action potentials along the pfs for PCs located at different distances from the site of granule cell layer activation. Based on in vivo results (cf. Fig. 3), inhibitory model parameters were chosen so that PCs on top of the site of stimulation increased firing rates by >20% followed by inhibition, whereas PCs between 200 and 400 µm from the site of granule cell layer activation showed a reduction in firing rate of >20%. At greater distances (>400 µm) firing rates of modeled PCs were required to be not statistically different from background (P < 0.05).
|
|
|
A fundamental prediction of the analysis presented in Fig. 8, which is a direct result of matching the model to experimental results, is the differential spatial distribution of stellate- and basket-type inhibitory synapses onto PCs (Fig. 8, C and E). For PCs farther away than 400 µm, suppression of the PC beam was not strongly dependent on the range of inhibitory synaptic delays and was mostly achieved by stellate-type synapses. However, between 0 and 400 µm, pf excitation occurred too quickly (<2 ms) to be compensated for solely by stellate-type inhibition. In this region, accurate simulation of experimental data required basket-type inhibitory inputs. Even then, the ability of this inhibition to shape PC spiking occurred over a fairly wide range of synaptic delays (Fig. 8D). Thus our model makes the prediction that the strength of basket- and stellate-type synaptic inhibitions onto PCs is spatially differentiated from the site of granule cell layer stimulation, with basket-type synapses more prominent closer than at greater distances along the pfs.
Biophysical effects of pf excitation and molecular layer inhibition on PC dendrites
The PC spiking output involves a complex interaction between synaptic inputs and the cell's intrinsic electrical and chemical properties (Llinás and Sugimori 1980
; Womack and Khodakhah 2002
). Modeling studies demonstrated that intrinsic properties produce considerable complexities in the action potential generation process in PCs (De Schutter and Bower 1994c
; Etzion and Grossman 1998
; Jaeger and Bower 1999
; Miyasho et al. 2001
; Santamaria and Bower 2005
; Watanabe et al. 1998
). In this section we focus on the contribution to this discussion provided by interactions between the large intrinsic dendritic currents found in PCs and the relative timing and strengths of excitatory and inhibitory synaptic currents generated by the network model.
Figure 9A is a cartoon representation of the timing and spatial relationships between the different excitatory and inhibitory network components of the model. Figure 9, BG shows the average synaptic, dendritic, and dendrosomatic currents, as well as the resulting changes in somatic firing for PCs found above and at three different distances from a focal granule cell layer activation as indicated in Fig. 9A (n = 64 simulations). The traces shown in Fig. 9, BE are the total summed dendritic currents for each set of channels (Jaeger et al. 1997
; see METHODS). Figure 9B shows the total granule cell excitatory input in these simulations. This input is provided primarily by ascending segment synapses for the PC at 0 µm and by pf inputs for all other PCs. The steepness of the rise time and the total amplitude of this excitatory input decrease with distance along the pfs as a direct result of the progressive pf desynchronization. Although excitatory synaptic currents have the fastest rise time and largest amplitude for the PC immediately overlying the activated region of the granule cell layer, the induced currents are still large enough even for the most distant PC to generate action potentials in the absence of feedforward inhibition (Fig. 3, EH).
Figure 9, C and D shows total inhibitory currents generated by stellate- and basket-type synaptic inputs, respectively. As discussed above, replicating the experimental data requires an inverse correlation of the strength of basket- and stellate-type synapses as a function of distance from the site of granule cell layer activation, with basket-type synapses stronger at shorter distances. It should be noted that, even though at short distances (0 and 400 µm) there are no stimulus-activated stellate-type synapses, there is still a net outward flow through the associated channels. This flow arises from current shunting by background-activated stellate-type synapses in response to pf excitatory input that, in effect, results in an immediate damping response to the inward current flowing through these excitatory synapses (Jaeger et al. 1997
). The outward flow of current is, of course, bigger at greater distances when stellate-type synapses are actually active.
The large intrinsic voltage-dependent currents in the PC dendrite and soma have a powerful influence on action potential generation (Jaeger et al. 1997
). Furthermore, the interaction between dendritic and somatic voltage-dependent conductances is complex (Jaeger and Bower 1999
) and its full discussion is beyond the scope of the current paper. However, summed dendritic and somatic currents for these simulations are shown in Fig. 9, E and F, respectively. In Fig. 9E, excitatory granule cell synaptic input results in a local dendritic depolarization that in turn activates dendritic currents, including predominantly CaP (red) and Kca (black) (Santamaria et al. 2002
). Ultimately, the influence of these intrinsic currents in the dendritic tree is reflected in the currents flowing between the soma and the dendrite (Fig. 9F). This somatodendritic current is then directly reflected in the spiking response of the soma (Fig. 9G). For the PC on top of the site of stimulation (0 µm) an initial flow of current from dendrite to soma increases the firing rate at the soma. Activation of basket-type somatic inhibition combined with Kca conductances in the dendrite, however, quickly reverse this current to flow from soma to dendrite, truncating the excitatory spiking response of the PC. For PCs within 400 µm, basket-type inhibition directly on the soma dominates dendrosomatic current flows, resulting in a sharp decrease in PC spiking output even though there is a large excitatory synaptic current in the dendrite (Fig. 9B). At greater distances along the pfs, the temporal balance between pf excitatory inputs and stellate-type inhibition at the level of the PC dendritic tree results in a more effective compensation of CaP currents by Kca currents, which in turn results in no net change in the somatodendritic current and therefore no net change in the average firing frequency of the PC. It is important to note, however, that synaptic input in all modeled PCs generates large changes in dendritic currents even when there is no change in average spiking activity in the soma. In fact, the dendritic currents in Fig. 9E show that considerably more current generally flows into the dendrite than is actually measured at the soma. The model suggests that the most important issue for PC spiking thus is not the size of synaptic inputs, but the relative timing of the balances between both synaptic (pf and molecular layer interneuron) and intrinsic (CaP and Kca) currents.
The importance for somatic spiking of the timing of activation of dendritic conductances is specifically shown in Fig. 10. The averaged conductance traces (n = 64 simulations) associated with granule cell (black) and both basket-type (green) and stellate-type (red) synapses in Fig. 10A show that at 0 µm, the short but still finite intrinsic delays in the feedforward inhibitory circuit do not compensate for excitatory input. As activity progresses down the pfs, however, the relative timing of the excitatory and inhibitory (initially basket-type and eventually stellate-type) conductances come more into register, resulting in an almost coincident average onset of excitatory and inhibitory influences even with an intrinsic inhibitory delay (1,200 µm). This somewhat counterintuitive effect is the result of the temporal dispersion of the pf volley and the bottom-to-top progression of the excitatory and thus the inhibitory wave.
As expected, the intrinsic dendritic conductances show a strong dependency on the timing of network effects, which are very sensitive to the presence of inhibition. Figure 10, B and C compares conductance traces obtained in the presence (solid) and absence (dashed) of feedforward inhibition. Figure 10B shows that PCs directly overlying the activated granule cell layer, whether in the presence or absence of feedforward inhibition, generate a substantial CaP conductance, explained by the fact that this conductance is rapidly activated by the excitatory input, rising faster than the time-delayed Kca conductance, and peaking before basket-type inhibitory influences can affect the soma. Figure 10C shows a stronger consequence of inhibition in the initial onset of the Kca conductance with a 50% net reduction in peak amplitude. At greater distances along the pfs, however, the increasing desynchronization of the pf input coupled with the counteracting effects of feedforward inhibition are sufficient to generate a slower response of the CaP, allowing its drive on membrane potential to be better compensated for by Kca. In the absence of feedforward inhibition, however, pf input drives a much faster onset in the CaP, reducing the effectiveness of the lagging Kca conductance (see channel kinetics in Supplementary Material Table S1). This is what allows inputs of pfs to directly influence somatic spiking. It is important to note that all these processes regulating the somatodendritic current take place at the level of the dendrite and bear no resemblance to classical synaptic integrative mechanisms that assume a simple summing to threshold of synaptic conductances at the soma.
| DISCUSSION |
|---|
|
|
|---|
Modeling parameters
Desynchronization of action potentials traveling along pfs can come from multiple sources and not only from a linear distribution of propagation velocities as assumed in our study (Vranesic et al. 1994
), which has been challenged based on anatomical studies (Wyatt et al. 2005
). For example, volleys of action potentials induced by synchronous activation of the granule cell layer would still be desynchronized arising from the fact that deep pfs branch first before the superficial ones in which action potentials travel longer along the ascending segment of the axon. Experimental evidence shows that even direct electrical activation of the pfs results in a rapidly broadening field potential along the pfs, presumably because of variations in pf conduction velocities (Eccles et al. 1966a
).
Because there is very little information about the relative timing of molecular layer inhibition and pf excitation, this parameter was lumped and treated as a free variable around which we used parameter-searching techniques to arrive at values producing the best matches to experimental data. With the exception for PCs found overlying or near a focus of granule cell layer activation, the interaction between pf and molecular layer interneurons was robust to synaptic timing. In this case the model suggests that fast somatic basket-type inhibitory input may specifically serve to block pf activation of the PC soma.
Several features of cortical circuitry could potentially affect the timing of inhibition that were not included in the model or are not known, such as molecular layer interneurons coupled by gap junctions (Mann-Metzer and Yarom 1999
; van der Giessen et al. 2006
) and inhibitory synaptic effects on each other (Pouzat and Marty 1999
). Although gap junctions could contribute to a more rapid spread of inhibitory influence, synaptic interactions could suppress it. We hope that the current study leads to more experimental and modeling interest in the role of molecular layer inhibition in cerebellar processing.
Although we have begun to study the effects of more complex patterns of granule cell layer activation on our models (Santamaria and Bower 2005
), in the present work we focused on establishing the conditions under which a pf-driven PC beam was most likelythat is, the synchronous activation of a small region of the granule cell layer. This kind of input produces the most synchronous activation of pfs, maximizing the potential consequences of the additional synaptic delay in the feedforward inhibitory circuitry. Even under these conditions, inhibition was fully capable of counteracting the effect of activation along the pfs. A more gradual and temporally or spatially spread out activation of the granule cell layer can be expected to further reduce the consequences of the additional synaptic delay in the inhibitory pathway. It is worth noting that, experimentally, the more prolonged in vivo granule cell layer responses evoked by peripheral stimulation do not produce beams of activated PCs (Bower and Woolston 1983
).
Experimental results
Although it is very easy to remove inhibition from a model, it is more difficult to do so experimentally with precision. We used two different drugs to minimize the likelihood that our results could be attributed to the specific properties of either (e.g., Debarbieux et al. 1998
). We showed that normally unresponsive PCs generated an excitatory response in the presence of either bicuculline or gabazine. We also showed that neither drug affects the basal level of PC activity nor the basic timing of either expected or emergent excitatory PC responses. We also demonstrated that these drugs do not affect the spatial organization of tactile inputs to the underlying granule cell layer. However, previous experiments in vivo (Chadderton et al. 2004
) and in vitro (Hamann et al. 2002
) reported that GABAA-receptor blockers can have an effect on granule cell layer responsiveness, presumably resulting from interference with Golgi cell inhibition. In our data, neither drug affects the onset of the granule cell layer response, but both show a small increase in granule cell layer responsiveness at longer latencies. This result could suggest, for example, that Golgi cell inhibition may more selectively influence secondary activation of the granule cell layer by cortico-cerebellar pathways (Morissette and Bower 1994
). It is also possible that the small effects on granule cell layer responsiveness we see arise from the lack of penetration of the drugs from surface application or are a consequence of some other aspect of the granule cell layerGolgi cell circuitry (Tahon et al. 2005
).
Anesthetics are also known to alter the balance between excitation and inhibition (Sonn and Mayevsky 2006
; Whittington et al. 2000
) and therefore could affect the interpretation of our experimental results. However, the lack of PC beams was previously shown in both anesthetized and unanesthetized preparations (Bower and Woolston 1984; Cohen and Yarom 1998
). Furthermore, there is very little difference between the amplitude of granule cell layer responses to tactile stimulation between unanesthetized and ketaminexylazine-anesthetized rats (compare Fig. 2, trace 2 with plots in Figs. 4A, 5A, and 7D; Hartmann and Bower 2001
). Therefore there is no evidence to suggest that the lack of PC beams is dependent on anesthetic effects.
Although none of the PCs recorded in these experiments showed excitatory responses to distantly activated regions of the granule cell layer under control conditions, not all distant PCs revealed excitatory responses after the application of GABAA-receptor blockers. This lack of response cannot be attributed to GABA-blocker changes in patterns of granule cell layer activity because we specifically controlled for that possibility. However, we previously reported that PCs directly on top of the site of granule cell stimulation were also sometimes unresponsive (Lu et al. 2005
). We interpreted this result to suggest that PCs might not be uniform with respect to their responses to granule cell layer activation. This variation in PC responsiveness does not appear to apply to responses to pfs because, again, there is no direct experimental evidence that individual PCs not overlying activated regions of the granule cell layer respond to pf inputs, unless inhibition has been artificially blocked.
Another issue that could be raised with respect to the interpretation of our experimental results is whether some excitatory influence other than the pfs could be responsible for the emergent responses. However, there are numerous reasons to believe that the emergent PC responses arise from pf activity. First, all recorded PCs were located in the center of Crus IIa, on the mediolateral axis of the pfs, and PC dendritic trees traverse most of the width of the crown of Crus IIa. Therefore it is essentially certain that pfs originating in the activated region of the granule cell layer traversed through the postdrug-excited PCs.
Second, there is no other known pathway that would produce this type of excitatory response in PCs. Climbing fibers, for example, are activated by tactile stimuli at much longer latencies (2035 ms; Brown and Bower 2001
, 2002
). Whereas our single-unit recordings also would have revealed the complex spike response of the climbing fiber, this input originates outside the cerebellum and there is thus no reason why topically applied bicuculline or gabazine would affect their responses. The only other possible synaptic mechanism for inducing short-latency excitatory responses would require a change in the spatial pattern of granule cell activity in response to peripheral stimulation after drug application. If such a change occurred, then the uncovered responses could result from ascending granule cell input, and not from the pfs. However, we explicitly controlled for this possibility by showing that the granule cell layer both nearby and beneath the recorded PCs is not activated by the stimulus either before or after GABAA receptors are blocked. This finding is also consistent with previous experiments demonstrating that even peripheral lesions do not result in an unmasking of new receptive fields within these regions of the granule cell layer (Shumway et al. 1999
).
In addition to these exclusory arguments, using either of the GABAA-receptorblocking drugs and two different calculation methods the values of pf conduction velocities are similar and consistent with previous reports in the literature. This is true whether the estimate is based either on the mediolateral distance from the recorded PC to the largest recorded granule cell layer response to the stimulus or by compared response latencies in two simultaneously recorded PCs.
Finally, the timing and shape of the uncovered responses are also consistent with a direct effect of the pfs. Specially, these responses are generally of lower amplitude and longer duration that those found immediately above activated regions of the granule cell layer. First, this difference is consistent with the fact that pf synapses are less synchronized than the ascending segment synapses converging on overlying PCs. Second, it is likely that fewer pf synapses are activated on each PC from a distant focal excitation than are activated by ascending segment synapses on PCs overlying an activated region of the granule cell layer. Ekerot and Jorntell (2001) estimated that perhaps as few as roughly 2% or 3,000 pf synapses are activated on PCs under these conditions, whereas our electron microscopic studies suggest that a tactile stimulus may activate as many as 25,000 ascending synapses per PC (Gundappa-Sulur et al. 2000). It is also important to point out that these experiments were specifically undertaken to test the predictions of a model, in which the uncovered responses are, for certain, a result of the pfs. Thus the model as well stands as evidence that the uncovered responses are attributed to the pfs.
Functional implications
CONTEXT FOR THESE STUDIES.
Most theories of cerebellar cortical function continue to assume that pfs are the primary influence on PC somatic output (Albus 1971
; Anastasio 2001
; Barto et al. 1999
; Bernard and Axelrad 1991
; Braitenberg and Atwood 1958
; Daya and Chauvet 1999
; Hofstotter 2002
; Marr 1969
; Medina and Mauk 2000
; Schweighofer et al. 1998
; Yamamoto 2002
), even though there is considerable experimental evidence that this is not the case. Although numerous experiments demonstrated that direct electrical molecular layer stimulation of the pfs can produce beamlike responses (Coutinho et al. 2004
; Dunbar et al. 2004
; Eccles et al. 1966b
; Gao et al. 2003
; Ito and Kano 1982
), almost all efforts to obtain beams of PCs under more natural stimulus conditions have failed. For example, Bell and Grimm (1969)
reported a lack of pf-induced correlative firing in PCs separated by more than a few tens of microns, a result also recently reported by Jaeger (2003)
. Eccles and colleagues found patchesnot beamsof activated PCs in the cat cerebellum after peripheral tactile stimulation (Eccles et al. 1972
). Our own tactile stimulation mapping experiments found activated PCs only immediately overlying activated regions of the granule cell layer in the rat (Bower and Woolston 1983
; Lu et al. 2005
), a result also subsequently demonstrated in the cat (Kolb et al. 1997
) and the isolated guinea pig cerebellum (Cohen and Yarom 1998
). In fact, in the last 50 years, we are aware of only one paper, purported to demonstrate beamlike effects in vivo using direct mapping techniques (Garwicz and Andersson 1992
), attributed the previous reports of a lack of PC beamlike responses to possible experimentally induced damage to pfs. However, our present results show that tactile stimuli delivered under exactly the same experimental conditions as those in previous studies (Bower and Woolston 1983
) are fully capable of producing excitatory PC beams as long as inhibition is blocked.
Unfortunately, Garwicz and Andersson (1992)
and several more recent in vitro studies (Coutinho et al. 2004
; Diez-Garcia et al. 2005
; Heck 1999
; Vranesic et al. 1994
), questioning previous in vivo reports of a lack of pf-induced PC beams, also made the assumption that the presence of field potential responses or dendritic currents directly implies a pf excitatory effect on PC somatic firing. The membrane current records fr