Research Article

Optical Imaging of Neuronal Populations During Decision-Making

See allHide authors and affiliations

Science  11 Feb 2005:
Vol. 307, Issue 5711, pp. 896-901
DOI: 10.1126/science.1103736


We investigated decision-making in the leech nervous system by stimulating identical sensory inputs that sometimes elicit crawling and other times swimming. Neuronal populations were monitored with voltage-sensitive dyes after each stimulus. By quantifying the discrimination time of each neuron, we found single neurons that discriminate before the two behaviors are evident. We used principal component analysis and linear discriminant analysis to find populations of neurons that discriminated earlier than any single neuron. The analysis highlighted the neuron cell 208. Hyperpolarizing cell 208 during a stimulus biases the leech to swim; depolarizing it biases the leech to crawl or to delay swimming.

Understanding the mechanisms of behavioral choice would be a major step in bringing together neuroscience, psychology, and ethology (1). Research into decision-making has used several different strategies. One very productive approach is to have a primate make a sensory discrimination between very similar stimuli while the activity of neurons in various parts of the nervous system is recorded (28). A second approach uses choice competition: presenting an animal with two stimuli that produce mutually exclusive behaviors (choices), to see which behavior predominates (9). This has led to the notion that behavioral choices are hierarchical. The neuronal mechanism originally proposed to underlie behavioral hierarchies was inhibitory interactions among the neurons responsible for triggering the different behaviors (10). Later work has found that neurons capable of eliciting one behavior are often activated during other, sometimes conflicting, behaviors (11, 12). Among other things, this observation suggests that individual decision-making neurons can be multiplexed—they contribute to choosing more than one behavior—and that they trigger behaviors by being active with other combinations of neurons.

We used a third approach to study decision-making: choice variability. We presented a nervous system with identical stimuli that repeatedly produce two different, mutually exclusive behaviors with roughly equal probabilities. This approach allowed us to focus on neurons involved in decision-making that are downstream from neurons used to make sensory discriminations.

We used the isolated central nervous system (CNS) of the medicinal leech. Motor neuron activity patterns characteristic of swimming (13) and crawling (14) can be elicited from isolated preparations by electrically stimulating peripheral nerves. Such sensory stimulation activates mechanosensory neurons in patterns that mimic touching the leech's skin (15). Stimulating the same kinds of mechanosensory neurons in different locations on the leech produces characteristic behaviors like swimming or crawling (16, 17). We follow the terminology proposed by Schall (18), referring to the different behavioral outputs as choices and the process leading up to a choice as decision-making.

Previously, recording from neurons intracellularly one at a time, then stimulating them to determine their effect on the initiation of behavior, has successfully uncovered interneurons that activate swimming (19, 20), crawling (21), and whole-body shortening (22). However, to explore how decisions are made by populations of neurons (11), we needed to record from many neurons at once (23). We therefore used voltage-sensitive dyes (24, 25) that allowed us to record simultaneously from many neurons in a midbody segmental ganglion at a resolution better than 5 mV (26).

The leech CNS makes behavioral choices. The isolated leech CNS consists of a nerve cord connecting 21 segmental ganglia plus a head and tail brain (Fig. 1A). This preparation generates motor patterns that are recognizable as behaviors observed in intact leeches, including swimming and crawling. We exposed one ganglion between ganglia 7 (G7) and 10 (G10) for voltage-sensitive dye imaging. We also recorded extracellularly from at least two peripheral nerves [dorsal posterior (DP) nerves] using suction electrodes. These electrodes were used to both stimulate and record. A train of electrical pulses to a DP nerve mimics a touch to the body wall in an intact leech (15) and can elicit both swimming and crawling (Fig. 1B). By stimulating DP nerves between G13 and G16, we evoked swimming in about half of the trials (blue) and crawling in the other half (red).

Fig. 1.

Recording behavioral choices in the isolated leech CNS. (A) Schematic of the sites of recording and stimulation. A midbody (G7 to G10) ganglion was imaged with a voltage-sensitive dye. Suction electrodes were used to record from and stimulate DP nerves (DP13 to DP15). (B) Eight sequential 60-s trials demonstrated intertrial variability. The stimulus (2 to 3 V, 10-ms pulses at 15 Hz) lasted 300 ms (small black bar). A ganglion was imaged during the initial 10 s of each trial (green bar). Trials are color-coded by behavior: swimming (blue, ∼1- to 2-Hz bursts) or crawling (red, ∼0.05- to 0.1-Hz bursts). (C) The population of neurons (143 in this example) from which we recorded on the ventral surface of a ganglion. (D) We averaged the pixels from each neuron to produce a time-varying record of the percentage change in the fluorescence signal (measured as dF/F). The top two traces are the signals from the two dye molecules, coumarin and oxonol. The largest and least noisy signal is the ratio of these two signals (third trace); it was used for all further analysis. The simultaneous intracellular recording (bottom trace) demonstrates the high sensitivity of the dye (15 to 20% per 100 mV). Pot., potential. (E and F) Raw optical data from two trials, one which elicited (E) swimming and the other (F) crawling. The optical signal for each neuron is plotted versus time. Color encodes the percent change in fluorescence (dF/F): Positive changes (red) correspond to relative depolarization, and negative changes (blue) correspond to relative hyperpolarization. The panels below each raster plot are simultaneous DP nerve recordings. Vertical black lines indicate the onset and duration of the stimulus.

Each trial lasted 60 s, with an intertrial interval of three min. Although the nerve activity was recorded for the entire trial, the ganglion was imaged during only the initial 10 s (Fig. 1B, green bar). We imaged the neurons on the ventral surface of a midbody segmental ganglion (Fig. 1C). In each preparation, we were able to resolve 130 to 150 of the ∼160 known neurons (27) on the ventral surface. The fluorescence resonance energy transfer (FRET)–based voltage-sensitive dye we used (24, 25) is very sensitive to small membrane potential fluctuations (Fig. 1D). For comparison, we show the raw data from a swim trial and a crawl trial (Fig. 1, E and F). Almost every neuron was activated immediately after each stimulus. There was a clear difference in the activity of many neurons once the motor pattern was apparent in the nerve recordings, usually after ∼4 s. These neurons are presumably central pattern generating (CPG) interneurons or motor neurons that generate the swimming or crawling rhythm. We were more interested in activity differences before this time, between the stimulus and 4 s, when the decision between the two behaviors was made.

Discrimination by single neurons. The activity patterns of a subpopulation of single neurons were able to discriminate swimming trials from crawling trials. Of the neurons that responded to the stimulus, we observed four classes of responses: nondiscriminating (ND) cells, early discriminating (ED) cells, late discriminating (LD) cells, and transiently discriminating (TD) cells (Fig. 2A). We quantified the earliest discrimination time of each single cell (tSC) by performing a sliding window analysis of variance (ANOVA) for each cell (Fig. 2B).

Fig. 2.

Single-cell discrimination. (A) Examples of the four classes of discrimination responses observed. Each graph plots overlapping raw fluorescence traces from five swimming (blue) and five crawling (red) trials for a single cell. In this and all subsequent plots, the gray shaded region denotes the time and duration of the stimulus. (B) A sliding-window ANOVA was used to quantify the discrimination times for each cell. The earliest discrimination time (tSC) was the time at which the swimming and crawling trajectories significantly diverged (black arrow). (C) A raster plot of the ANOVA results from all of the cells from one experiment. Black indicates nonsignificant times, and white indicates significantly different times (P < 10–6). The discrimination time based on the nerve recording of the behavior is shown (tNERVE, green line). (D) A histogram of the earliest tSC from the raster plot in (C). Seventeen neurons produced significantly different trajectories before tNERVE.

Neurons were ordered by their earliest discrimination times (Fig. 2C). The number of cells that discriminated at some point in time (ED, LD, or TD cells) ranged between 50 and 75% across experiments. We also performed an ANOVA on the nerve recordings (28) to determine the time (tNERVE) at which we could discriminate the behaviors on the basis of motor neuron activity (Fig. 2, C and D, green line). We were most interested in cells that discriminate before tNERVE because these cells are predictive of the behaviors and are candidate decision-making neurons. A histogram of the earliest tSC for each cell shows that 17 cells (mean ± SD = 19 ± 6 cells; n = 6 preparations) discriminated before tNERVE (Fig. 2D).

Discrimination by populations of neurons. We view the nervous system of the leech as a dynamical system (29, 30). Abstractly, the behavioral state of the nervous system at any instant is a point in the phase space of this system (Fig. 3A). Each axis represents a variable that measures the temporal evolution of trajectories in this space (31, 32). For us, these variables are linear combinations of observed neurons. As in the single-cell analysis, we labeled the time at which the trajectories have significantly diverged as the discrimination time. The decision-making period must occur at or before this time. We recognize that the eventual choice may depend on the behavioral state before stimulation, but we did not address this issue in these experiments.

Fig. 3.

Population discrimination. (A) A schematic of our conceptual framework. The state of the leech CNS at any instant is a point in a d-dimensional phase space (d = 3 in the drawing). The temporal evolution of its state defines a trajectory in this space. Before stimulation, the state of the nervous system is in a rest region. Upon stimulation, the trajectories diverge toward either a swimming or a crawling region on the basis of which choice was made. (B) The first PC from a single experiment. Each PC is a linear combination of the observed neurons in the N-dimensional data set (N = 143 in this experiment). (C) The same data as in (B), projected onto the three PC axes. Each trajectory represents a trial (swimming blue, crawling red). The average swimming and crawling trajectories are in bold. Black dots indicate 1-s intervals along the average trajectories. The right panel plots the simultaneous DP nerve recording for each trial. The black arrow illustrates the direction of a linear discriminant estimated for a 500-ms time bin in the PC space. The direction of the linear discriminant is a weighted vector sum of the three PCs with weights x, y, and z. (D) The resulting linear discriminant, a linear combination of the three PCs. This direction is similar to that of PC3 in (B). (E) The linear discriminant weightings for each of the neurons projected onto a map of the ganglion. Red indicates large-magnitude (both positive and negative) contributions to the linear discriminant direction; blue indicates small contributions.

The single-cell discrimination analysis gave us an idea of the number of potential candidate decision-making cells. However, the purpose of simultaneously recording populations of neurons was to look at dynamic interactions among them. Using the multiple recordings, we asked whether a linear combination of neurons could discriminate earlier in time than any single neuron. To address this question, we used two analysis techniques in conjunction: linear discriminant analysis (LDA) (33) and principal component analysis (PCA) (28, 3436).

PCA identifies dimensions that separate behaviors. For the experiment shown in Fig. 3C, we performed PCA on a data set containing 143 neurons (dimensions) measured across 14 trials. We show the first three principal components (PCs) (Fig. 3B). The bar graphs represent the dimensionless contribution of each neuron to each PC. Values close to zero indicate a small contribution, and high positive or negative values indicate a large contribution. The first three PCs typically account for 60 to 80% of the overall variance in a data set. The remaining 140 PCs are ignored here for visualization purposes. Effectively, we have reduced the dimensionality of the data set from 143 (neurons) to 3 (linear combinations of neurons, the PCs).

We plotted the data three-dimensionally, using the first three PCs as the axes (Fig. 3C). The trajectories all start in one region and then diverge toward two different regions of the space. From this plot, we extract two features: (i) The separation between the trajectories is an objective measure of decision-making. (ii) The neurons that contribute most to this separation (i.e., have large positive or negative values) are most likely to be responsible for making the choice.

Occasionally, we observed a trial such as the one shown in green in Fig. 3C. This trial would have normally been classified as a crawling trial on the basis of the nerve recording. However, the trajectory in the phase space initially diverges in the direction of the swimming trajectories. After a delay, the trajectory turns and moves toward the direction of the crawling trajectories. Although trials such as this were rare (2 out of 60 trials from six experiments), we interpret this as evidence that decision-making is a dynamical process: The leech nervous system can start to make a decision and subsequently change to an alternative choice.

LDA finds neuronal populations that discriminate. Having reduced the dimensionality of the data set with PCA, we asked (i) at what time did the swimming and crawling trajectories significantly diverge, and (ii) which neurons are responsible for the separation at this time? We thus divided the temporal trajectories into time bins and estimated the linear discriminant for each bin. We then performed an ANOVA on the data in each respective time bin projected onto each linear discriminant. Therefore, each time bin has an associated P value. The time at which this P value became significant is denoted tLDA. The time bin width and the number of PCs were varied to find the optimal tLDA for each data set (37) (fig. S1).

We show an estimated linear discriminant in Fig. 3D using three PCs. Neurons contributing strongly to the linear discriminant were those that best helped discriminate between swimming and crawling at tLDA. In this case, PC3 had the largest weighting, so the linear discriminant direction is similar to PC3.

We wished to visualize the spatial locations of the neurons with large-magnitude contributions to the linear discriminant. We used color to encode the magnitude of each neuron and projected the colors onto a map of the ganglion (Fig. 3E). This spatial ganglion map was used as a guide for identifying candidate decision-making neurons.

Discrimination of single neurons versus neuronal populations. We compared the discrimination times of the earliest single cell (tESC) and the linear discriminant (tLDA) from six experiments (Fig. 4A). In all experiments, tLDA occurred before tESC. tSC times (black lines) that occurred before the tNERVE times (green lines) are also plotted.

Fig. 4.

Comparing single-cell discrimination to population discrimination. (A) The distribution of discrimination times from six experiments. Discrimination times for the LDA (tLDA, red), single cells (tSC, black), earliest single cells (tESC, cyan), and nerve recordings (tNERVE, green) are shown. The average difference between tLDA and tESC was 290 ± 60 ms (mean ± SEM). Single-cell discrimination times occurring later than the tNERVE for each experiment are omitted. (B) The color-coded linear discriminant weightings from a single experiment (color-coding as in Fig. 3E). The black arrow indicates cell 208. (C) The color-coded single-cell discrimination times (tSC) from the same experiment. Yellow represents the earliest discrimination times, and red represents later discrimination times. (D) The result from one experiment. Cell 208 was intracellularly impaled and alternately depolarized (+1.5 nA) or hyperpolarized (–1.5 nA) for the initial 10 s of each trial (black bar). The membrane resistance of cell 208 at rest was ∼50 MΩ. A nerve shock was delivered for 300 ms as in Fig. 1B. In the four trials in which cell 208 was hyperpolarized (blue), the preparation swam within 10 s. In the four remaining trials in which cell 208 was depolarized (red), the preparation either produced the crawling motor pattern or it delayed swimming until the intracellular stimulation ended. (E) A contingency table summarizing the pooled results from five preparations. The values in parentheses are the expected counts if the observed behavior did not depend on the stimulus condition.

The neuronal populations contributing highly to the linear discriminant were generally different from the neurons with early single-cell discrimination times (Fig. 4, B and C). Single-cell discrimination times (tSC) in Fig. 4C are color-coded, with yellow representing the earliest tSC times and red representing later tSC times. The ganglion maps from all experiments are given in figs. S2 and S3.

Cell 208 biases decisions. The analyses described to this point were performed within 15 min during an ongoing experiment, so the ganglion maps were used to identify candidate decision-making neurons. To test whether these neurons were sufficient to influence decision-making individually, we passed polarizing current into each of them during the nerve shock. We impaled each candidate neuron and injected hyperpolarizing or depolarizing current before and after the nerve shock. None of the neurons with early single-cell discrimination times (Fig. 4C) significantly affected the elicited behaviors (33 neurons from six preparations). When we tested neurons contributing strongly to the linear discriminant (17 neurons from six preparations), we found a neuron, cell 208, that can selectively bias the decision to swim or to crawl (Fig. 4B, arrow, and fig. S3). When depolarized or hyperpolarized alone, this neuron did not initiate any behaviors (28). However, when stimulated during a nerve shock, cell 208 biased the decision toward swimming or crawling (Fig. 4D). With cell 208 hyperpolarized, the nerve shock reliably evoked swimming (blue trials); with it depolarized, the nerve shock evoked crawling or delayed swimming (red trials). In five preparations, the correlation between the level of current injection and the observed behavior (Fig. 4E) was significant (P < 0.01, Fisher's exact test). We conclude that this neuron plays a role in decision-making. The neuron was labeled with an intracellular dye in all experiments and identified as cell 208 on the basis of its morphology and electrophysiological properties (38).

Discussion. One of the central questions we attempted to address is whether decision-making is performed by single neurons or by neuronal populations (39). One extreme proposal is multiple competing circuits in which decision-making neurons for one behavior act by inhibiting the other behaviors, so that only one behavior occurs at a time (9, 10). At the other extreme would be the complete sharing of decision-making neurons by two or more behaviors, with the dynamics of the network determining which behavior is chosen. In the leech, we hypothesize a middle ground in which decision-making neurons are partially shared, where the dynamics of neuronal populations can determine choices, but individual neurons in these populations can profoundly influence decision-making. This view is supported by our results that (i) neuronal populations discriminated between swimming and crawling earlier in time than did single neurons; (ii) none of the single neurons that discriminated early were able to bias the decision; and (iii) although cell 208 was able to bias the decision, it was part of a population defined by the linear discriminant.

The finding that a linear combination of neurons discriminates early in time shows that there is sufficient information in a population of neurons for an experimenter to predict the ultimate choice. However, this result demonstrates only that the activity pattern of this network correlates with the eventual choice. To show that this information is used during decision-making, we needed to manipulate the network. Ideally, we would like to selectively hyperpolarize or depolarize a population of individual cells; however, we were technologically limited to stimulating one or two cells at a time.

Although we have not exhaustively tested all candidate neurons detected, most of the neurons tested were not able to bias the behavioral choice. However, we were able to reliably bias the choice with one of the linear discriminant candidates, cell 208. A relatively large depolarization or hyperpolarization of cell 208 biased the choice toward crawling or swimming, respectively. Although this result was statistically significant, cell 208 did not determine the choice on all trials. This less-than-perfect control by cell 208 can probably be credited to the cell 208 homo-logs in each of the 20 other ganglia (38) that we did not stimulate.

Cell 208 has previously been described as a CPG neuron (38), but a unique one that connects the swim-initiating network to the swim CPG circuit. It has not been shown to trigger behaviors in intact isolated nerve cords. Our results suggest that it is also part of a decision-making circuit, although we do not yet know its role in this circuit or the mechanism by which it biases the system. Although it is possible that cell 208 is driven by a higher-order decision-making neuron, we have shown that it alone is capable of biasing the entire system.

Most recent research about decision-making has focused on value-based choices in which there is always a right and wrong answer (28). We propose that the term decision-making should refer to a spectrum of goal-driven behaviors. At the most complex level are conscious, introspective choices that incorporate expected reward (1). At the simplest level are reproducible, predictable reflexes. There is a large area of involuntary, subconscious decision-making that has been neglected. For example, suppose your goal is to walk down the street. Which foot do you lead with? This is a choice that is made without conscious effort and will vary from time to time. Although this choice is goal-directed, there is no correct or incorrect choice. In our reduced preparation, there is no obvious value associated with swimming versus crawling, although both choices would achieve the goal of escaping from a stimulus.

Why then does the leech, a relatively simple nervous system, not respond in a predictable manner? We hypothesize two possibilities: Either the choice depends on the rest state before each stimulus, or the state is reset upon stimulation and then diverges stochastically because of noise in the system. The leech may provide a system in which to resolve this important question.

Supporting Online Material

Materials and Methods

SOM Text

Figs. S1 to S3

References and Notes

References and Notes

View Abstract

Navigate This Article