QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

HOME  |   SEARCH  |   ARCHIVE  |   SUBSCRIBE  |   CONTACT  |   HELP

This Article Abstract Full Text (PDF) Submit an eLetter Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (20) Citing Articles via Google Scholar Google Scholar Articles by Leznik, E. Articles by Llinás, R. Search for Related Content PubMed PubMed Citation Articles by Leznik, E. Articles by Llinás, R.

 Previous Article  |  Next Article 

The Journal of Neuroscience, April 1, 2002, 22(7):2804-2815

Electrotonically Mediated Oscillatory Patterns in Neuronal Ensembles: An In Vitro Voltage-Dependent Dye-Imaging Study in the Inferior Olive

Elena Leznik, Vladimir Makarenko, and Rodolfo Llinás

Department of Physiology and Neuroscience, New York University School of Medicine, New York, New York 10016

	ABSTRACT

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

Spatiotemporal profiles of ensemble subthreshold neuronal oscillation were studied in brainstem slices using high-speed voltage-sensitive dye imaging. After local electrical stimuli, the overall voltage profile demonstrated coherent oscillatory waves that spread over the inferior olive (IO). These oscillations were also observed in concurrently obtained intracellular recordings from IO neurons. Over the first few seconds after the stimuli, the optically recorded oscillations clustered into coherent groups comprising hundreds of neurons. Statistical analysis of the spatial profiles of these clusters revealed size fluctuation around stable core regions that were surrounded by a rim the diameter of which varied in time during the oscillation period.

The neuronal ensemble oscillations were calcium derived and had an average frequency range of 1-7 Hz. This rhythmic response demonstrated a different spatiotemporal distribution in the presence of picrotoxin, which induced the merging of neuronal clusters into larger areas of coherent activity. The possibility that such clustering is a consequence of intrinsic oscillations in ensembles of coupled neurons was tested using mathematical modeling.

Key words: inferior olive; oscillations; voltage-sensitive dyes; imaging; picrotoxin; mathematical modeling

	INTRODUCTION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

The inferior olive (IO) is a bilateral symmetrical nucleus located in the ventral aspect of the caudal bulbar region of the brainstem. IO neurons project to the contralateral cerebellum and form one of the two major cerebellar afferents known as climbing fibers (Szentagothai and Rajkovits, 1959). Climbing fibers synapse onto the cerebellar Purkinje cells (PCs) in the cerebellar cortex (Ramón y Cajal, 1888) and form one of the largest synaptic junctions in the nervous system (Eccles et al., 1966). Functionally, the IO nucleus is known to be essential in motor coordination (Wilson and Magoun, 1945; Murphy and O'Leary, 1971; Kennedy et al., 1982); its lesions abolish coordinated movements and result in motor disorders similar to those that follow cerebellar damage (Llinás et al., 1975; Soechting et al., 1976). However, the manner in which the olivocerebellar system contributes to cerebellar function is still a subject of extensive debate (for review, see Bloedel and Bracha, 1998; De Zeeuw et al., 1998).

Controversy regarding the function of the olivocerebellar system is attributable partly to its wide responsiveness to peripheral stimuli and its low (1-10 Hz) single-cell firing rate (Eccles et al., 1966). On the basis of its low firing rate, it has been argued that the olivocerebellar system cannot affect cerebellar output directly and its only function is to modify the strength of parallel fiber-Purkinje cell synapses (Keating and Thach, 1995). An alternative view suggests that the olivocerebellar system modifies cerebellar output directly by temporally binding widespread regions of the cerebellar-nuclear system and thereby by specifying combinations of cereberellar outputs in time (Llinás et al., 1975). According to the latter hypothesis, IO activity functions as a timing signal for movement by generating synchronous activation of Purkinje cell outputs (Llinás and Sasaki, 1989; Welsh and Llinás, 1997). Experiments with multiple electrode recordings have further supported this hypothesis. Indeed, rhythmic inhibitory potentials can be recorded at the cerebellar nuclei after olivocerebellar activity (Llinás and Muhlethaler, 1988).

The ability of the olivocerebellar system to generate ensemble synchronous activity has been attributed to the intrinsic properties of the IO neurons (Llinás and Yarom, 1981a,b, 1986; Benardo and Foster, 1986; Bal and McCormick, 1997) and their electrotonic coupling (Llinás et al., 1974; Sotelo et al., 1974; Llinás and Yarom, 1981a; Makarenko and Llinás, 1998). In particular, the interplay between several types of voltage-dependent calcium and potassium conductances enables the IO cells to fire rhythmically at 1-10 Hz. The combination of these intrinsic oscillatory properties with the electrotonic coupling between olivary cells results in the subthreshold membrane potential oscillations (Lampl and Yarom, 1997; Schweighofer et al., 1999). The coherence of such oscillations enables a group of IO neurons to fire in synchrony and simultaneously activate several regions of the cerebellum during movement coordination (Welsh et al., 1995). Therefore, knowledge about the spatiotemporal profiles of the oscillations in the IO nucleus provides additional information regarding the mechanisms responsible for the temporal binding of motoneuron activation and thus of motor coordination.

This question is addressed here by using a voltage-sensitive dye-imaging technique. This technique is presently the methodology of choice in addressing the geometrical distribution of activity in a large neuronal ensemble (Cohen et al., 1978). We report here that ensemble oscillations in the IO emanate from clusters of coherent activity, where each cluster is composed of hundreds of cells. Given the distribution of complex spike (CS) activity in the cerebellar cortex, the clusters are the most likely origin for synchronous activation of Purkinje cells observed in previous in vivo multiple recording experiments. In addition, we compared our experimental results with those obtained by computational modeling of IO neuronal ensembles endowed with oscillatory electrical properties and electrotonic coupling. These modeling results indicate that neuronal oscillatory clustering is a direct consequence of the combined electrotonic/intrinsic properties of coupled IO neurons.

	MATERIALS AND METHODS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

Inferior olivary slices

Brainstem slices were prepared from postnatal day 13-20 Sprague Dawley rats following protocols from previous in vitro studies with some modifications (Llinás and Yarom, 1981a; Bleasel and Pettigrew, 1992). Animals were anesthetized with 15-20 mg of ketamine and decapitated. The brainstem was isolated and placed in cold, oxygenated Krebs'-Ringer's solution containing (in mM): 126 NaCl, 5 KCl, 1.25 NaH₂PO₄, 2 MgSO₄, 2 CaCl₂, 10 glucose, and 26 NaHCO₃. Parasagittal slices (300 µm in thickness) were sectioned using a vibratome. The major portion of the IO was contained in one or two slices. In an effort to increase the viability of slices, they were sectioned in sucrose-substituted Krebs'-Ringer's solution in which NaCl was replaced by sucrose while the osmolarity of the solution was kept constant (Aghajanian and Rasmussen, 1989). Slices were then transferred to a chamber containing normal oxygenated Krebs'-Ringer's solution and incubated at room temperature for at least 1 hr. Optical recording was implemented on slices stained with voltage-sensitive dye RH414 (Molecular Probes, Eugene, OR) dissolved in HEPES-based Ringer's solution containing (in mM): 130 NaCl, 6.25 KCl, 2 MgCl₂, 2 CaCl₂, and 25 HEPES, to a final concentration of 2.5 mg/ml. Slices were incubated for 10 min in the dye solution before being transferred to the interface recording chamber. Experiments were conducted at a bath temperature of 33 ± 1°C.

Recordings of optical signals

A schematic diagram of the recording setup is shown in Figure 1. The recording chamber, attached to an upright microscope (Olympus BX50WI), was illuminated with a halogen lamp (12 W) driven by a stable power supply (Kepco, Inc.). The microscope was mounted onto an air table (Warner Instruments). The excitation light passed to the preparation through a bandpass filter (515 ± 35 nm) via a dichroic mirror. The emitted light returned through a long-pass filter (>590 nm). Optical signals were monitored with a fast CCD camera (HR Deltaron 1700, Fujix) with a 128 × 128 pixel spatial resolution. Images were sampled every 0.6-4.8 msec depending on the experiment. Most of the data presented were taken at 4.8 msec sampling rate. The total area imaged was 2.7 × 2.7 mm, and each pixel collected light from a surface of ~21 × 21 µm. Optical recordings were averaged over four or eight trials. For each trial, the base fluorescence level (F_o) was initially calculated by averaging 64 frames. Changes in membrane potentials were evaluated as F/F (Grinvald et al., 1982). Bleaching and irregularities of staining and illumination were corrected off-line by using Matlab-based software (The Mathworks, Inc.). In brief, the optical signals were first detrended to compensate for bleaching of the dye and slow responses from glia cells (Lev-Ram and Grinvald, 1986; Konnerth et al., 1987). The signals were then filtered with a three-dimensional moving average (3 × 3 × 3) and with the Gaussian low-pass filter. Finally, the optical signals were displayed using the RGB (red-green-blue) 256 color scale in such a way that their maximum amplitude equaled the maximum red color intensity of the RGB scale. The recordings were analyzed using power-spectrum density and autocorrelation methods written in Matlab.

View larger version (69K): [in this window] [in a new window]   Figure 1.   Schematic diagram of the experimental setup. Light from a 12 V halogen source was passed through an excitation filter (515 ± 35 nm), dichoic mirror, and microscope objective (5×) before reaching the slice stained with the voltage-sensitive dye RH-414. Emitted fluorescent light was projected onto a CCD camera after passing through the objective, dichroic mirror, and cutoff filter (>590 nm). The CCD camera (Fujix, HR Deltaron 1700) consisted of 128 × 128 pixels, and each pixel collected light from a surface of ~21 × 21 µm. Images were sampled every 4.8 msec. The optical data were analyzed off-line using Matlab-based software. An example of analyzed image is shown. Its color calibration bar corresponds to 0.001% change in F/F.

One or two electrical stimuli (0.1-2 mA, 0.2 msec) were applied to the dorsal border of the IO nucleus via bipolar concentric electrodes. In five experiments, intracellular recording and optical recordings were simultaneously acquired to calibrate the optical signals with respect to directly recorded membrane potential. Intracellular recordings were obtained using glass micropipettes filled with 3 M KAcetate (60-100 M). Recordings were amplified with an Axoclamp-2A amplifier (Axon Instruments) and acquired at 10 kHz with an Instrunet board (Omega Engineering, Inc.).

Mathematical modeling of inferior olivary oscillations

The motivation behind modeling of neuronal activity clusters was that of testing the assumption that such clustering, and cluster size, are electrotonic coupling dependent. Indeed, although this assumption may be considered obvious at first glance, it has been demonstrated that activity profiles other than clusters (such as spirals, waves, concentric rings) can occur in coupled neuronal networks (Makarenko, 1994).

The experimental data provided parameters such as membrane potential oscillation (from intracellular recordings) and cluster size and shape and their modulation by pharmacological intervention (from voltage-dependent dye imaging). These parameters stipulated the type of model to be implemented. Picrotoxin effect was modeled as change in electrical coupling, although the drug does not affect coupling between IO neurons directly but rather indirectly by shunting electrotonic current flow at the IO glomerulus (Lang et al., 1996). The single-cell model was developed to mimic the quasi-periodic subthreshold membrane potential oscillation observed electrophysiologically (Velarde et al., 2002). Network properties were modeled by interconnecting a set of these modeled elements into a two-dimensional matrix with isotropic coupling properties. Such a model was more easily implemented and scaled than models that incorporate single-cell cable and ionic conductance properties (Manor et al., 1997; Schweighofer et al., 1999).

Single element model. We represent single neurons as distinct elements (points) in a rectangular connectivity lattice connected to four other elements.

Each element is endowed with quasi-periodic oscillatory properties defined by a set of equations having a periodic and a noise term as follows:

(1)

where the variable z characterizes neuron dynamics (with x = Re (z), i.e., z = x + iy; prime implies differentiation), i is an imaginary unit, x and y are real numbers,  is a damping constant, w₀ is angular oscillation frequency in the absence of noise and damping (w₀ = 2 * 10 Hz), D is a parameter determining noise intensity scaling, and (t) is a noise term that has zero mean with a time correlation function given by:

(2)

where brackets imply average, and  denotes a function that is zero when its argument has a non-zero value.

An example of oscillatory behavior of a single element and its power spectrum profile is given in Figure 2. Relatively regular oscillations have variable amplitude with a sharp frequency peak in the region of 10 Hz, which is in good agreement with in vivo experimental data (Sasaki et al., 1989).

View larger version (18K): [in this window] [in a new window]   Figure 2.   Time series and power spectrum of a model neuron (A). Subthreshold oscillations are almost periodic and have a sharp frequency peak at ~10 Hz (B). The width of the power spectrum peak and some irregularity in the oscillations are caused by introduction of noise with zero mean.

The network model. The network consisted of a two-dimensional (n × n) lattice of single-modeled neurons, as described above, with their electrotonic coupling properties represented by periodic boundary conditions given by:

(3)

where the pair (jk) denotes site in the lattice and d^lm_jk accounts for the electrotonic coupling coefficient (jk) and (lm). The d^lm_jk is assumed to be the same over the network, i.e., isotropic. The actual value of the coefficient determines the dynamics of the network. The sum in the right part of Equation 3 is taken over neighboring neurons:

(4)

where R accounts for the radius of neuronal interaction. For nearest-neighbor coupling r = 1. The lattice (Eq. 3) is able to produce oscillations with a well defined frequency band peaked around w₀ with relatively slowly varying amplitudes.

Markov random field as an objective pattern descriptor. The introduction of voltage-dependent dye imaging to determine the geometry of the neuronal activity in vitro raised the necessity to qualify ensemble activity patterns in an objective manner. Indeed, in most imaging studies the interpretation of imaging pattern is based on subjective criteria that are open to ambiguous interpretation. In an attempt to alleviate this problem, we used a statistical approach known as the Markov random field (MRF) methodology, originally developed to characterize levels of spatial organization in two-dimensional matrices. This procedure, used often in structural biology (Grundy et al., 1997) and solid-state physics (Besag, 1974), has also proved useful in the analysis of brain imaging (Held et al., 1997; Wang et al., 2001) and multidimensional neuronal ensemble activity (Makarenko et al., 1997).

MRF offers the possibility of quantifying spatial configuration of clustering and stability of a particular spatial configuration with respect to their fluctuations. Thus, if pixel amplitudes are treated as values for the two-dimensional matrix elements, the MRF estimations can be determined unambiguously. The application of MRF estimates two main parameters:  and . However, in this study only  is relevant, because it presents a quantitative estimation of the level of clustering in the image, whereas  relates to the probability that a particular value may occur in a given element of the matrix. When the matrix is presented as a rectangular lattice with values of element x_i,j placed at the i,j node of the lattice, then:

(5)

where y _i,j = x_i1,j + x_i+1,j + x_i,j1 + x_i,j+1.

is one of two sublattices that constitutes the original lattice and each of which is distributed as the black and white spaces of a chessboard.

The eventual estimation of the  is taken as an arithmetic mean over the two subsets. The detailed description of the method with an illustration of how  correlates to the spatial organization in an image has been addressed previously (Makarenko et al., 1997). Here we review parameter : the more the absolute value of  differs from zero, the higher is the level of spatial organization in the image. The "spatial organization" is codependent on the values of the neighboring pixels. The degree of the codependence averaged over the whole matrix reflects the level of spatial organization that exists at a given moment in the total system. Calculating  for a set of consecutive image sequences as a function of time in the course of an experiment we can see how it characterizes the temporal evolution of the spatial organization, bringing together the spatial and temporal features of the network activity. The example illustrated in Figure 7 compares three configurations consisting of zeros and ones with different degree of clustering from absolutely random distribution (no clusters) and how such clusters reflected in the value of . Values of  for each spatial configuration in Figure 7A are given in the figure legend. It must be remembered that the actual values of  are model dependent and only meaningful within a given model, i.e., within a matrix of the same dimension and value ranges. Thus, comparison of modeled and experimental findings must be based on statistics. In the present case, the comparison is based on the variance of .

	RESULTS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

Characteristics of the optical signals

The spatiotemporal characteristics of the ensemble neuronal oscillations in inferior olive slices were studied using optical recordings of voltage-sensitive dye signals. Several steps were taken to increase the signal-to-noise ratio of our imaging responses. First, optical recordings were averaged over four or eight successive trials. Second, the recording sequence was structured in such a way that the electrical stimuli to the IO nucleus were delivered a short time (50 msec) after the beginning of image acquisition. This allowed the electrical stimuli to reset subthreshold oscillation phase and thereby to entrain a large proportion of neurons to in-phase oscillations (Llinás and Yarom, 1986). Because each camera pixel integrates voltage signals over 21 × 21 µm and includes several neurons, only in-phase synchronous activity can be detected readily. Indeed, synchronization of oscillatory activity over the IO network increased the amplitude of the optical signal to the level that can be easily detected.

Because the present experiments were designed to record stimulus-evoked oscillations averaged over several trials, an issue then arises concerning the relation of spontaneous oscillations with stimulus-evoked, averaged oscillations in the IO. This question was addressed directly with intracellular recordings (Fig. 3).

View larger version (28K): [in this window] [in a new window]   Figure 3.   Comparison of spontaneous and stimulus-evoked oscillations in the inferior olive. A, Intracellular recording of spontaneous oscillations at 2 Hz interrupted by an extracellular stimulus. After extracellular stimulation (marked with an arrowhead), the oscillations disappeared for 750 msec (boxed area) and then resumed. The membrane potential was 60 mV. B, Intracellular recordings of spontaneous (dashed black line) and stimulus-evoked (solid black line) oscillations from the same cell are superimposed. Their corresponding power spectra are shown below. Note that extracellular stimulation only modified the phase of spontaneous oscillations without affecting their amplitude and frequency. C, Six individual intracellular traces of stimulus-evoked oscillations from the same cell are superimposed on the left. Each trace is shown in a different color. Their corresponding power spectra are displayed below. In every recording, the frequency of stimulation-evoked oscillation was the same (2.0 Hz). Note that in each trace the stimulation- induced shift in the oscillatory rhythm of the cell is remarkably similar. Oscillations are clearly seen after the stimulus-induced reset but can be barely detected before the stimulation. D, The average of six traces of stimulus-evoked oscillations (red line) and the recording of spontaneous oscillations (dashed black line) are superimposed. The stimulus-evoked oscillations in the average trace have the same frequency and amplitude as the spontaneous oscillations and differ only in the phase shift. Calibration bar: 1 mV; A, 1 sec; B, D, 500 msec; C, 415 msec.

In accordance with previous results (Llinás and Yarom, 1986), an extracellular stimulus delivered at the dorsal border of the IO nucleus generated a full action potential followed by a membrane hyperpolarization in nearby neurons (Fig. 3A, boxed area). The findings also demonstrate that if the cell was oscillating at the time of the stimulus, its rhythmicity was momentarily stopped but then resumed ~750 msec after the stimulation with a different phase (n = 7). It is important to note that the extracellular stimulation only reset the phase of oscillations, without affecting their frequency or amplitude (Fig. 3B). This electrical behavior could be obtained repeatedly for any given cell, and the stimulation-induced shift in the oscillation phase was remarkably similar (Fig. 3D) (n = 6). Moreover, for a given cell, the average of six individual stimulus-evoked oscillations had the same frequency as that of the spontaneous oscillations (Fig. 3D). When the cell was stimulated by a train of stimuli, the results were similar to those shown for a single stimulus, but the reset time was prolonged (data not shown). Thus, the stimulus-evoked IO oscillations averaged over several trials have the same frequency and amplitude as spontaneous oscillations and differ only in a phase shift. Therefore, characteristics of these oscillations can be used to further understand spontaneous oscillations in the IO.

All the slices used in this imaging study showed oscillatory activity (n = 42, 30 rats). As stated above, the evoked optical response comprised an initial change in background fluorescence near the stimulating electrode (up to 0.027%) followed by several oscillatory cycles (Fig. 4A, Control) (n = 14).

View larger version (49K): [in this window] [in a new window]   Figure 4.   Spatiotemporal patterns of optically recorded inferior olivary oscillations. A, A complete time course of optical responses for four pixels under control conditions (left) and in the presence of 100 µM CdCl2 (right). Locations of the pixels are marked with circles of different color in B. Under control conditions, the optical response consisted of an initial response to a train of extracellular stimuli (2 stimuli at 10 Hz) followed by three oscillatory cycles. The frequency of the optically recorded oscillations was 4 Hz. Addition of CdCl2 blocked optically recorded oscillations; only direct depolarization caused by the extracellular stimulation was detected. B, Several multiframe displays for each of the regions marked in A. Only the area of the IO is shown. Position of the stimulating electrode is indicated. The color bar gives the color scale with dark purple corresponding to lowest fluorescence and red corresponding to highest fluorescence. All images were filtered by a low-pass (30 Hz) Gaussion filter. 1, The response of the IO nucleus to the first stimulus. 2, The response of the IO to the second stimulus. Note that the time course and location of the response to the second stimulus are identical to that of the first, but with a larger change in the background fluorescence. 3, 4, The location and spatial spread of two cycles of the ensemble oscillations. The oscillations emanate from several clusters of closely spaced cells throughout the IO (see the change in amplitude with time). These clusters are highly repeatable from cycle to cycle over an oscillation sequence. 5, An average of four oscillatory cycles further underlines the stability of fluorescent clusters. 6, Same as in 4, but in the presence of 100 µM CdCl2.

The initial change in fluorescence indicated the location, amplitude, and time course of the response of the IO neurons to the stimulus (Fig. 4A, Control, B1). When a train of two stimuli was applied to the IO, the time course and location of the response to the second stimulus were identical to that to the first stimulus, but the second stimulus induced a larger change in the background fluorescence (Fig. 4, compare B1, B2). These initial responses were followed by several ensemble oscillations that spread throughout the IO nucleus with little spatial or timing variance (Fig. 4A,B3,B4). For the most part, the oscillatory sequence repeated itself with a small degree of variability after each stimulus train. The fluorescence clusters emanated from several independent regions of the slice and consisted of closely spaced synchronous neuronal ensemble activity. All regions had the same oscillation frequency and were phase coherent. This phase coherence was most clearly observed after the phase reset triggered by the electrical stimulus. The oscillations were highly repeatable in time and space: the average of three oscillatory cycles had the same spatiotemporal profile as that of the individual oscillations (Fig. 4, compare B3, B4, B5). The average oscillation had the same fluorescent clusters at the same regions of the slice as the individual oscillations.

Within our temporal resolution (4.8 msec between frames), the optical oscillations started synchronously everywhere in the slice (Fig. 4A, Control, 1, 2, 3, 4), suggesting that the conduction velocity of the oscillation spread was higher than the temporal resolution of the recording camera. In addition, because of temporal resolution, we were low-filtering the neuronal activity and could see only slow subthreshold events.

Ionic mechanisms for cluster oscillation

To determine the ionic mechanism supporting the oscillatory activity of these neuronal clusters, three experiments were conducted in the presence of CdCl₂ (100 µM) in the bathing solution. CdCl₂ is known to block Ca²⁺ ionic conductance subserved by calcium channels and thereby to stop subthreshold oscillations in the IO (Llinás and Yarom, 1986). Addition of cadmium to the bathing solution completely blocked optically recorded oscillations (Fig. 4A, right panel, B6). Under this condition, only short-lasting focal activation of the area around the stimulating electrode was detected (Fig. 4A, right panel). This result indicates that the optically recorded ensemble oscillations are Ca²⁺ dependent and thus are generated by the same conductances responsible for single-cell subthreshold oscillations recorded intracellularly in the IO neurons (Benardo and Foster, 1986; Llinás and Yarom, 1986; Bleasel and Pettigrew, 1992).

Comparison of optical signals and intracellular recordings

To correlate the time course of our optical signals with intracellular recordings, we compared these two kinds of measurements in five experiments. The subthreshold oscillations were first recorded intracellularly in the presence of the dye (Fig. 5A). Location of the recording electrode is marked with an asterisk in Figure 5B. The slice was then electrically stimulated with a bipolar electrode at the border of the IO nucleus, and the evoked ensemble IO oscillations were imaged (Fig. 5B). Note that the optically observed oscillations consisted of several clusters that changed its dimensions depending on the oscillation phase. They embraced the largest area during the upward phase of the oscillations (Fig. 5B, second and fourth image from left).

View larger version (37K): [in this window] [in a new window]   Figure 5.   Comparison between optically recorded and intracellularly recorded IO oscillations. A, Subthreshold oscillations recorded intracellularly from an IO cell in the presence of dye (RH-414). Position of the recording electrode is marked with an asterisk in B. B, Frames from optically recorded oscillations from the same slice as in A. Two cycles of oscillations are shown. The position of the stimulating electrode is indicated. C, Autocorrelograms and power spectra of the optically recorded ensemble oscillations (dashed black line) and intracellular recorded oscillations (solid black line). Note that the clusters, seen as spots of fluorescence in the image panel, have oscillatory voltage profiles at frequencies similar to those observed intracellularly.

The intracellularly and optically recorded signals were analyzed using power-spectrum density and autocorrelation methods written in Matlab. As shown in Figure 5C, the frequencies of the intracellular recorded subthreshold oscillations (dashed black line) and the optically recorded neuronal ensemble oscillations (solid black line) (2.0 and 1.8 Hz, respectively) were closely matched. A slightly slower frequency of the cluster oscillations with respect to the intracellular oscillations was observed. This small difference in oscillatory frequency is probably attributable to the fact that the borders of the clusters oscillated at a slower frequency than the individual cells. Because movement of the cluster borders involves sequential activation of neighboring cells and such activation can occur only during the upward phase of the oscillations, the fluctuation of the borders would be slower than oscillations of the individual cells.

The overall similarity of intracellular and cluster oscillations indicates that the subthreshold IO oscillations observed intracellularly can be studied using voltage-sensitive dyes. Moreover, given the temporal and spatial characteristics of such rhythmic activity and its phase reset properties (Makarenko and Llinás, 1998), such optical measurements most probably reflect the ensemble electrophysiological properties observed previously with multiple simultaneous recordings at the Purkinje cells level (Llinás and Sasaki, 1989; De Zeeuw et al., 1998; Lang et al., 1999; Fukuda et al., 2001; Yamamoto et al., 2001).

The spatial and temporal analysis of the clusters

IO cells are interconnected by gap junctions and can exhibit spontaneous oscillatory activity at a frequency that ranges from 1 to 10 Hz (Llinás et al., 1974; Sotelo et al., 1974, Benardo and Foster, 1986; Llinás and Yarom, 1986; Bleasel and Pettigrew, 1992; De Zeeuw et al., 1996). It has been suggested that the IO nucleus consists of several functionally coupled oscillating clusters of cells (Llinás and Yarom, 1986; Lampl and Yarom, 1993; Devor and Yarom, 2000). In support of this view, in all experiments (n = 14) the recorded optical responses produced oscillating clusters that were time coherent (Figs. 4A,B, 5B). Given the possible functional implications of such dynamic structures, the spatial and temporal parameters of the clusters were calculated. In particular, the frequency and average size of a cluster for each experiment were analyzed in detail. Statistical analysis of the properties of a set of clusters from 14 different experiments determined mean size, variances, SE, and oscillatory frequency (Table 1). The oscillation frequency was determined by using power spectrum density analysis (see Materials and Methods). The range of frequencies observed was from 1.6 to 6.5 Hz. These results demonstrate similarities between the optically recorded oscillatory frequency range and those reported for intracellular subthreshold oscillations (Benardo and Foster, 1986; Llinás and Yarom, 1986; Bleasel and Pettigrew, 1992).

View this table: [in this window] [in a new window]   Table 1.   Description of dimensions of the optically recorded inferior olivary clusters

The mean cluster size was determined by choosing a threshold fluorescence level that was 2 SDs above the mean background noise. In every case the chosen threshold level distinguished very adequately areas of optical activity from background noise. A cluster was then defined as an area that showed oscillatory activity with pixel values above the threshold. Each cluster consisted of a core region and the adjoining area (Figs. 4B, 5B). The core region demonstrated a close to constant size, but the extent of the adjoining area was found to be phase dependent. The core area and maximum area (i.e., the core region plus the adjoining area at its uppermost extent) were calculated for several representative clusters in each experiment (Table 1). The mean core area was 11 ± 6 × 10³ mm² , and the mean maximum area was 20 ± 10 × 10³ mm². Because the size of the cluster is thought to be directly related to the level of functional coupling among IO neurons, the difference in the number of intact gap junctions between slices can account for high SDs in the size of observed clusters.

Because most of the IO neurons are coupled via dendrodendritic gap junctions, which are not located within the same plane (Sotelo et al., 1986; De Zeeuw et al., 1990), we assumed that the optically recorded IO clusters were three-dimensional structures. We then estimated the number of cells in each cluster by multiplying the area of the cluster by the thickness of the response region by the density of the IO cells (Table 1). Both the core area and the maximum area of a cluster consisted of hundreds of cells. On average, there were 260 ± 140 cells in the core region and 490 ± 250 cells in the maximum area of the cluster (Table 1). Although this measurement is only an approximation, it does show statistical consistence across all the experiments.

Thus, our optical data indicate that at the network level, the IO nucleus is organized in functionally coupled activity clusters. Each cluster is composed of several hundreds of cells that may act in unison to simultaneously activate groups of cerebellar Purkinje cells.

Factors affecting cluster dimensions

In the final set of experiments. we attempted to define the mechanisms that determine cluster dimensions. Several possibilities were considered. The first possibility was that the level of functional coupling among the IO neurons defined cluster boarders (Welsh and Llinás, 1997). If this were the case, addition of drugs that modulate electrical coupling among olivary neurons, such as picrotoxin (Lang et al., 1996), would change the area of optically recorded clusters. The second possibility was that a relatively small number of all the cells survived in the slice, and thus, the recorded responses were the product of surviving clusters of neurons.

To distinguish between these two possibilities, several experiments were conducted in the presence of picrotoxin (10 µM). At the single-cell level, addition of picrotoxin reduced the frequency and amplitude of spontaneous subthreshold oscillations without affecting the membrane potential of the cell (Fig. 6A,B). The frequency of oscillations was reduced by 15 ± 4%, and the amplitude was decreased by 51 ± 8% (n = 12). At the network level (n = 4), application of picrotoxin increased the size of the clusters by partly merging them into larger areas of activity and by partly activating additional neuronal regions (Fig. 6B-D). In our experiments, application of picrotoxin increased the cluster size by 207 ± 17%. These results demonstrate that clustering was not an artifact resulting from sparcity of neuronal survival, because neurons previously occupying quiescent sites became active after the drug administration.

View larger version (32K): [in this window] [in a new window]   Figure 6.   Effects of picrotoxin on intracellularly and optically recorded oscillations in the IO. A, Intracellular recording from an IO neuron showing spontaneous subthreshold oscillations before (left) and after (right) bath application of 20 µM picrotoxin. B, Power spectra of control (solid black line) and picrotoxin-modulated oscillations (dashed black line) are superimposed. Addition of picrotoxin reduced amplitude and frequency of subthreshold oscillations without changing the resting membrane potential of the cell. C, A frame of imaged ensemble neuronal oscillating clusters in control conditions (left) and after addition of 20 µm of picrotoxin (right). Several representative clusters in the bottom right part of the IO are outlined in black. The clusters were defined as areas with pixel values above the selected threshold level (in this case, 0.007%). After picrotoxin, an additional threshold level was chosen to delineate the highest areas of activity. D, Contours of cluster within the boxed area in A marked with white discontinuous lines (1 × 1 mm) are enlarged. E, Areas shown in B are superimposed. Control clusters are shown in blue; picrotoxin clusters of lower threshold are shown in yellow, and those of higher threshold are shaded in red. Overlapping regions of control clusters and picrotoxin clusters of lower threshold are shown in green. Note that picrotoxin significantly increased the size of clusters by merging several small areas into larger activation areas.

Our experiments with picrotoxin showed that the drug that modulated the level of the effective electrotonic coupling among IO cells changed the dimensions of inferior olivary clusters. In accordance with data from previous in vivo studies in which local application of picrotoxin (Llinás and Sasaki, 1989; Lang et al., 1996) or the destruction of cerebellar nuclei resulted in increments of coherent Purkinje cell activity, in our experiments application of picrotoxin significantly increased the size of fluorescent clusters by partly merging several small clusters into larger areas of activation. Thus, our results demonstrate that there is a basal level of inhibition in the IO slices supported by the presence of the presynaptic inhibitory terminals of brainstem and cerebellar origin, which are known to synapse at the IO glomeruli (Sotelo et al., 1986; De Zeeuw et al., 1996, 1998). Indeed, it is now recognized that about half of the cerebellar nuclei cells are GABAergic, and they project directly to the IO. Blocking of such inhibition within the IO glomeruli, where the gap junctions are located, expands the cluster size and increases the electrical load of the network and within the cluster. This increase is the most likely origin for the reduction in amplitude and frequency of spontaneous subthreshold oscillations observed with intracellular recordings.

In conclusion, our results indicate that the cluster size is probably determined by the IO electrical coupling coefficient and thus by the magnitude and distribution of the return inhibition from the cerebellar nuclear feedback as demonstrated in previous in vivo experiments (Nelson et al., 1989; Ruigrok and Voogd, 1995; Lang et al., 1996).

Model results: IO cluster size is coupling-coefficient dependent

To test the hypothesis that the modulation of electrical coupling is responsible for cluster size modeling, results were obtained by running the simulation protocol described in Materials and Methods. First, as illustrated in Figure 7, variations of the coupling d^lm_jk in the neuronal ensemble model (Eq. 3) were simulated. The results from direct observation of the phase coherence and the calculated value of such statistics as variance for MRF parameter  indicate that clustering of activity among the modeled matrix elements is coupling-coefficient dependent.

View larger version (83K): [in this window] [in a new window]   Figure 7.   Model results illustrate that  is a measure of clustering and that electrotonic coupling determines the value of . The different values for  in A and B reflect model differences (A, binary; B, nonbinary). A, Clustering property dependence on the value of the MRF parameter . Each square takes a binary value of 0 (white square) or 1 (black square). (1), Randomly distributed configuration with no clusters; corresponding 1 ~ 0. (2), Small clusters, 2 = 2. (3), Big clusters, 3 = 3.5 > 2 > 0. B, Clustering in the modeled network as function coupling dlmjk degree. (1), In the presence of reduced coupling (inhibition) in the system,  = 1.0. (2), With increased coupling (absence of inhibition),  = 1.5. Values of  differ for A and B because they relate to different models addressing two different issues. A addresses general model properties and B addresses specific IO modeling (for further explanation, see Materials and Methods).

The modeling results support the conclusion that modulation of electrotonic coupling in the IO by feedback inhibition from the cerebellar nuclei (Lang et al., 1996) regulates IO cluster distribution and size. That is, increased inhibition results in cluster size reduction and in an increase in the number of independent clusters (Fig. 7A). Conversely, decrease of inhibition results in increased cluster size and a simultaneous decrease in the number of independent clusters (Fig. 7B). Moreover, MRF results indicate that the variance (²) of the parameter  increases with reduction of coupling coefficient so that inequality _{weak inh} ~ 0.4 > _{strong inh} ~ 0.25 takes place.

These modeling results are in agreement with the imaging data. We calculated time series for the Markov parameter  and its variance from sets of imaging frames from seven experiments. All frames from complete sequences of optical images were analyzed during oscillatory sequences under normal conditions and after picrotoxin administration. The variance of  reflected cluster size fluctuation probability. Indeed, in the presence of background inhibition, the average  value (n = 7) was 1.7 × 10¹ and its variance (²) was ~1 × 10³, whereas in the presence of picrotoxin (no inhibition) (n = 3) the average  value was 1.4 × 10¹ and its variance (²) was twice the control value (~2 × 10³).

Thus in both modeling and experimental results, blocking of inhibition significantly increased the variance of , which indicated that although the number of individual clusters was decreased, the probability of cluster size fluctuations was increased.

	DISCUSSION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

Inferior olivary neurons are characterized by intrinsic rhythmic activity that results from the interplay between the low-threshold calcium spikes, the calcium-dependent potassium, and hyper-polarization-activated cationic conductances (Llinás et al., 1974; Sotelo et al., 1974; Llinás and Yarom, 1981a,b, 1986; Benardo and Foster, 1986; Bleasel and Pettigrew, 1992; Bal and McCormick, 1997). The combination of these intrinsic resonance properties and electrotonic coupling between the olivary cells via gap junctions results in subthreshold membrane potential oscillations (Llinás and Yarom, 1981b, 1986; Benardo and Foster, 1986; Bleasel and Pettigrew, 1992; Bal and McCormick, 1997; Lampl and Yarom, 1997). The olivary gap junctions are mainly located within the structures known as glomeruli (Llinás et al., 1974; Sotelo et al., 1974; King et al., 1975). A high density of gap junctions (De Zeeuw et al., 1996) and their specific location provide an anatomical substrate for generating simultaneous electrical discharges among neighboring IO neurons. Synchronous firing of a group of IO cells, in turn, can result in synchronized climbing fiber activation of Purkinje cells (complex spikes) throughout the cerebellar cortex. The potential for such global synchronization is shown by the fact that concurrent complex spike activity can be detected between the Purkinje cells located in widely separated regions of the cerebellar cortex (De Zeeuw et al., 1996; Yamamoto et al., 2001).

The ability of the olivocerebellar system to generate rhythmic synchronous discharges has been proposed to play a central role in motor coordination for several decades (Llinás et al., 1975; Llinás, 1991). However, although the IO oscillations are a neuronal ensemble event, they have been studied mainly at a single-cell level. In an attempt to understand these ensemble properties, high-speed voltage-sensitive dye imaging was implemented to define the spatiotemporal properties of the IO oscillations at the network level. Our results demonstrate that IO ensemble oscillations result from closely spaced, synchronously oscillating clusters comprising hundreds of neurons. Each cluster is a dissipative functional entity that acts in unison with other clusters and through its connectivity results in a synchronous afferent input volley to the cerebellum.

Description of the IO clusters

On the basis of the anatomical organization of the IO nucleus and widespread synchronicity of the complex spike activity in the cerebellar cortex, it has been proposed that on the functional level, the IO nucleus consists of clusters of electrotonically coupled neurons that can spontaneously generate rhythmic and synchronous outputs to the cerebellum (Welsh and Llinás, 1997). In agreement with that hypothesis, the present study visualizes these clusters and describes their spatial distribution and dynamic behavior in vitro. Extracellular electrical stimuli applied directly to the brainstem slice at the level of the IO were shown to reset the phase of spontaneous subthreshold oscillations and thereby entrain the IO neurons toward coherent in-phase oscillations (Llinás and Yarom, 1986). We showed with intracellular recordings that the IO stimulus-evoked oscillations had the same amplitude and frequency as spontaneous oscillations and differed only in phase shift.

IO oscillations were recently imaged in vitro for the first time (Leznik et al., 1999; Manor et al., 2000). In the study by Manor et al. (2000), the size of their imaging area allowed them to observe only a small number of IO cells. By increasing the number of IO neurons that can be observed simultaneously, we were able to increase the amplitude of the optical signals, improve the signal-to-noise ratio, and observe the entire extent of the IO nucleus.

Our optical responses consisted of initial activation of the area around the stimulating electrode followed by stimulus reset ensemble neuronal oscillations. There was always a close correspondence between the frequency of intracellularly recorded subthreshold oscillations (Benardo and Foster, 1986; Llinás and Yarom, 1986; Bleasel and Pettigrew, 1992) and that of our optical ensemble neuronal oscillations. The optically recorded oscillations emanated from several spatially separated fluorescent clusters. Every cluster had the same oscillation frequency and was phase coherent.

The observed clusters had a core region and the adjoining area that most likely consisted of loosely connected cells. The total area of a cluster, which included the core area and the adjoining area at its uppermost extent, was defined as a maximum area of the cluster. The core of the cluster was constant in size, but its amplitude was modified depending on the phase of the ensemble oscillations. In contrast, the adjoining area of loosely grouped cells changed its dimensions depending on the phase of the oscillations, so that its size was maximal at the peak oscillation amplitude. This phenomenon can be explained by the fact that with relatively few gap junctions that link cells in the adjoining area of the cluster, the amount of current flow needed to entrain these cells in oscillatory activity might be present only during the upward phase of the oscillations (Lang, 2001). As for the tightly coupled cells within the core of the cluster, the high density of the gap junctions may allow them to function in an all-or -one manner throughout the whole oscillatory cycle. Our data are in general agreement with those of Manor et al. (2000), who used a similar optical imaging technique and observed spontaneous oscillatory activity in the IO slices over the area of 80 × 10³ mm². These authors, however, did not detect multiple clusters in the IO nucleus, partly because their imaging system had a relatively low spatial resolution and covered only a small area of 600 × 600 µm versus the imaging area of 2700 × 2700 µm used in our study, which was about 20 times larger.

Implications of modeling results

The modeling results support our two main working hypotheses. First, the level of inhibition (i.e., coupling) is a main factor in regulating the average size of a cluster of synchronous activity. Second, with the higher level of inhibition, the number of clusters increases and the variance () of the MRF parameter  decreases. Because the variance of  is proportional to the fluctuations of the average cluster size, our results imply that inhibition decreases and stabilizes the average size of the clusters.

A significant implication of the first finding is that as the degrees of freedom increase with an increased number of clusters, the probability of large fluctuations in cluster size decreases. In fact, although counterintuitive, the system becomes "less liquid" with a larger number of smaller clusters. That is, despite the large number of clusters, the system robustly preserves a particular level of spatial organization. More fundamentally, the results indicate the existence of a maximum number of possible functional states, implying that increased inhibitory feedback does not result in an ever-increasing control ability. Equally interesting, even at high levels of clustering, the IO network demonstrates a short-phase resetting time (~70 msec) without losing oscillatory robustness, thereby optimizing reset agility without sacrificing control (Makarenko and Llinás, 1998).

Functional significance of IO cluster activity

Taken together, our results demonstrate that the voltage-sensitive dye-imaging technique can be used successfully to study the geometrical distribution of activity in the IO nucleus. We report here that on a functional level, the IO is organized in several coupled clusters of cells with borders that are most likely determined by the level of effective electrotonic coupling and therefore by the state of the IO network.

In an effort to determine the functional significance of the observed clusters, we calculated the maximum number of cells that contributed to each cluster. Because the clusters must be three-dimensional structures, we estimated the volume of each cluster by multiplying the area of the cluster by the thickness of the slice. This can be done because fluorescence is detectable up to 700 µm in depth (Grinvald et al., 1994), and our slices were only 300 µm thick. When multiplying such volume by the density of the IO neurons in rat (80,000 cells/mm³) (Cunningham et al., 1999), the product corresponds to the approximate number of cells in each cluster. We then made the parsimonious assumption that 50 µm of the exterior surfaces of the slice was probably damaged, and the response region was at maximum 200 µm thick and contained at least one cluster. On the basis of the above assumptions, we calculated that the IO clusters consisted of several hundred cells.

Questions then arise concerning the correlation between the size of the IO clusters and the extent of spontaneous synchronous climbing fiber activity in the cerebellum. Results from a number of previous in vivo experiments have shown that synchronized complex spike activity is primarily observed among Purkinje cells located within parasagittally oriented strips of the cerebellar cortex (Bell and Kawasaki, 1972; Llinás and Sasaki, 1989; Sugihara et al., 1995; Lang et al., 1997, 1999; Hanson et al., 2000; Yamamoto et al., 2001). These parasagittal bands are typically narrow structures with a width of only 0.25-0.5 mm. However, in the anterior-posterior direction, they can extend across several successive folia and are at least 4-5 mm long (Yamamoto et al., 2001). Therefore, each parasagittal strip can occupy an area from 1 to 2.5 mm². Because there are ~1200 PCs per square millimeter (Armstrong and Schild, 1970), the maximum number of PCs per parasagittal band would be 3000. Given that in rat there is an 8:1 ratio of the Purkinje cells to the IO cells, the IO clusters that define the parasagittal bands consist of ~400 cells. This number is of the same order of magnitude as the calculated number of cells in the core of optically registered IO clusters (260 ± 140 cells). Thus, the optical measurements might reflect the ensemble physiological properties observed previously with multiple recordings at the PC level. We propose that each IO cluster controls the activity of a single parasagittal band in the cerebellar cortex. Simultaneous activation of several IO clusters would result in synchronous CS activity in different parasagittal bands, as observed during coordinated movements.

	FOOTNOTES

Received Oct. 3, 2001; revised Jan. 18, 2002; accepted Jan. 23, 2002.

This work was supported by National Institutes of Health/National Institute of Neurological Disorders and Stroke Grant NS13742 and Department of Defense/Office of Naval Research Grant N00149911081. We thank Dr. Diego Contreras for the use of his image analysis programs and Dr. Francisco J. Urbano for critical reading of this manuscript.

Correspondence should be addressed to Rodolfo Llinás, Department of Physiology and Neuroscience, New York University Medical Center, 550 First Avenue, New York, NY 10016. E-mail: llinar01{at}popmail.med.nyu.edu.

	REFERENCES

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

• Aghajanian GK, Rasmussen K (1989) Intracellular studies in the facial nucleus illustrating a simple new method for obtaining viable motoneurons in adult rat brain slices. Synapse 3:331-338[ISI][Medline].
• Armstrong DM, Schild RF (1970) A quantitative study of the Purkinje cells in the cerebellum of the albino rat. J Comp Neurol 139:449-456[Medline].
• Bal T, McCormick DA (1997) Synchronized oscillations in the inferior olive are controlled by the hyperpolarization-activated cation current I(h). J Neurophysiol 77:3145-3156[Abstract/Free Full Text].
• Bell CC, Kawasaki T (1972) Relations among climbing fiber responses of nearby Purkinje Cells. J Neurophysiol 35:155-169[Free Full Text].
• Benardo LS, Foster RE (1986) Oscillatory behavior in inferior olive neurons: mechanism, modulation, cell aggregates. Brain Res Bull 17:773-784[ISI][Medline].
• Besag JE (1974) Spatial interaction and the statistical analysis of lattice systems. J R Stat Soc [Ser B 60] 36:172-236.
• Bleasel AF, Pettigrew AG (1992) Development and properties of spontaneous oscillations of the membrane potential in inferior olivary neurons in the rat. Brain Res Dev Brain Res 65:43-50[Medline].
• Bloedel JR, Bracha V (1998) Current concepts of climbing fiber function. Anat Rec 253:118-126[Medline].
• Cohen LB, Salzberg BM, Grinvald A (1978) Optical methods for monitoring neuron activity. Annu Rev Neurosci 1:171-182[ISI][Medline].
• Cunningham JJ, Sherrard RM, Bedi KS, Renshaw GM, Bower AJ (1999) Changes in the numbers of neurons and astrocytes during the postnatal development of the rat inferior olive. J Comp Neurol 406:375-383[Medline].
• Devor A, Yarom Y (2000) GABAergic modulation of olivary oscillations. Prog Brain Res 124:213-220[Medline].
• De Zeeuw CI, Holstege JC, Ruigrok TJ, Voogd J (1990) Mesodiencephalic and cerebellar terminals terminate upon the same dendritic spines in the glomeruli of the cat and rat inferior olive: an ultrastructural study using a combination of [3H]leucine and wheat germ agglutinin coupled horseradish peroxidase anterograde tracing. Neuroscience 34:645-655[ISI][Medline].
• De Zeeuw CI, Lang EJ, Sugihara I, Ruigrok TJ, Eisenman LM, Mugnaini E, Llinás R (1996) Morphological correlates of bilateral synchrony in the rat cerebellar cortex. J Neurosci 16:3412-3426[Abstract/Free Full Text].
• De Zeeuw CI, Simpson JI, Hoogenraad CC, Galjart N, Koekkoek SK, Ruigrok TJ (1998) Microcircuitry and function of the inferior olive. Trends Neurosci 21:391-400[ISI][Medline].
• Eccles JC, Llinás R, Sasaki K (1966) The excitatory synaptic action of climbing fibers on the Purkinje cells of the cerebellum. J Physiol (Lond) 182:268-296[Abstract/Free Full Text].
• Fukuda M, Yamamoto T, Llinás R (2001) The isochronic band hypothesis and climbing fiber regulation of motricity: an experimental study. Eur J Neurosci 13:315-326[ISI][Medline].
• Grinvald A, Hildesheim R, Farber IC, Anglister L (1982) Improved fluorescent probes for the measurement of rapid changes in membrane potential. Biophys J 39:301-308[Abstract/Free Full Text].
• Grinvald A, Lieke EE, Frostig RD, Hildesheim R (1994) Cortical point-spread function and long-range lateral interactions revealed by real-time optical imaging of macaque monkey primary visual cortex. J Neurosci 14:2545-2568[Abstract].
• Grundy WN, Bailey TL, Elkan CP, Baker ME (1997) Meta-MEME: motif-based hidden Markov models of protein families. Comput Appl Biosci 13:397-406[Abstract/Free Full Text].
• Hanson CL, Chen G, Ebner TJ (2000) Role of climbing fibers in determining the spatial patterns of activation in the cerebellar cortex to peripheral stimulation: an optical imaging study. Neuroscience 96:317-331[ISI][Medline].
• Held K, Kops ER, Krause BJ, Wells III WM, Kikinis R, Muller-Gartner HW (1997) Markov random field segmentation of brain MR images. IEEE Trans Med Imaging 16:878-886[ISI][Medline].
• Keating JG, Thach WT (1995) Nonclock behavior of inferior olive neurons: interspike interval of Purkinje cell complex spike discharge in the awake behaving monkey is random. J Neurophysiol 73:1329-1340[Abstract/Free Full Text].
• Kennedy PR, Ross HG, Brooks VB (1982) Participation of the principal olivary nucleus in neocerebellar motor control. Exp Brain Res 47:95-104[Medline].
• King JS, Martin GF, Bowman MH (1975) The direct spinal area of the inferior olivary nucleus: an electron microscopic study. Exp Brain Res 22:13-24