JN Journal of Applied Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 96: 471-477, 2006. First published April 5, 2006; doi:10.1152/jn.00628.2005
0022-3077/06 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
96/1/471    most recent
00628.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via ISI Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Münch, T. A.
Right arrow Articles by Werblin, F. S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Münch, T. A.
Right arrow Articles by Werblin, F. S.

REPORT

Symmetric Interactions Within a Homogeneous Starburst Cell Network Can Lead to Robust Asymmetries in Dendrites of Starburst Amacrine Cells

Thomas A. Münch1 and Frank S. Werblin2

1Helen Wills Neuroscience Institute, 2Molecular and Cell Biology, University of California, Berkeley, California

Submitted 17 June 2005; accepted in final form 1 April 2006


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DESCRIPTION
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Starburst amacrine cells in the mammalian retina respond asymmetrically to movement along their dendrites; centrifugal movement elicits stronger responses in each dendrite than centripetal movement. It has been suggested that the asymmetrical response can be attributed to intrinsic properties of the processes themselves. But starburst cells are known to release and have receptors for both GABA and acetylcholine. We tested whether interactions within the starburst cell network can contribute to their directional response properties. In a computational model of interacting starburst amacrine cells, we simulated the response of individual dendrites to moving light stimuli. By setting the model parameters for "synaptic connection strength" (cs) to positive or negative values, overlapping starburst dendrites could either excite or inhibit each other. For some values of cs, we observed a very robust inward/outward asymmetry of the starburst dendrites consistent with the reported physiological findings. This is the case, for example, if a starburst cell receives inhibition from other starburst cells located in its surround. For other values of cs, individual dendrites can respond best either to inward movement or respond symmetrically. A properly wired network of starburst cells can therefore account for the experimentally observed asymmetry of their response to movement, independent of any internal biophysical or biochemical properties of starburst cell dendrites.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DESCRIPTION
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The starburst amacrine cell, a retinal interneuron, plays a critical role in the circuit of direction-selective (DS) ganglion cells (Fried et al. 2002Go, 2005Go; Yoshida et al. 2001Go). DS cells spike vigorously when a visual stimulus transverses their receptive field in one direction but remain silent when the same stimulus travels in the opposite direction (Barlow and Levick 1965Go). Individual dendrites of starburst cells also show directional responses (Euler et al. 2002Go): The observed calcium signal in the tip of starburst dendrites was strongest when motion was directed "outward, " or centrifugal, i.e., from the cell body to the tip of the process. The response was minimal for opposite movement (inward, centripetal, from the tip to the cell body). This makes starburst cell dendrites some of the earliest components in the DS circuit that express directional response properties.

Several possible explanations of this phenomenon have been given. These include geometrical properties of starburst dendrites that lead to biophysical properties that favor centripetal movement (Borg-Graham and Grzywacz 1992Go; Tukker et al. 2004Go); asymmetric distribution of chloride co-transporters along the dendrites that lead to spatially asymmetric chloride currents and hence to directional asymmetries (Gavrikov et al. 2003Go); or cell-internal biochemical processes, calcium-induced calcium currents, which lead to stronger calcium signals for outward movement (Barlow 1996Go). These mechanisms are consistent with each other and could function synergistically.

These studies all account for the directional behavior of starburst dendrites by invoking intrinsic directional properties of starburst cell dendrites. However, it is not possible to determine if starburst cells are intrinsically directional because we study them deeply embedded within the retinal circuitry. Theoretically, however, we can define how intrinsic directional behavior could be distinguished from extrinsically imposed directional behavior. To do this, we imagine a starburst cell with all its synaptic input sites. We then activate these synapses sequentially, as if a light bar is sweeping across the cell, and compare the activity of a dendrite when the activity sweeps in opposite directions. It is important that the activity of any individual synapse does not depend on the direction of the sweep. We would consider the starburst cell to have intrinsic directional properties if, under these hypothetical conditions, the response of a dendrite is different for sweeps in opposite directions.

In this report we investigate if the retinal network can impose directional behavior on starburst cell dendrites, when the dendrites do not have intrinsic directional properties. To approach this problem, we constructed a computational model of an interacting network of starburst cells. The following two well-established findings about the geometric arrangement of starburst cells within the retina were incorporated in our model: first, starburst dendrites receive synaptic input along the whole length of the dendrite, but they release transmitter only at their distal third (Famiglietti 1991Go). Second, the dendritic trees of starburst cells overlap strongly (Famiglietti 1983Go; Vaney 1984Go). While the dendritic trees have diameters of ~300 µm, the distance between neighboring cell bodies is on the order of only 30 µm. This creates a dense dendritic network with many possible sites of interaction between starburst cells. Moreover, those interactions could be excitatory and inhibitory because starburst cells release both acetylcholine and gamma-aminobutyric acid (GABA) (Brecha et al. 1988Go; O’Malley et al. 1992Go) and also have receptors for both of these neurotransmitters (Feller 2002Go; Zhou and Fain 1995Go).


    MODEL DESCRIPTION
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DESCRIPTION
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We tested the influence of such excitatory and/or inhibitory network interactions on the directional response properties of starburst dendrites. For simplicity, our model network is one dimensional, consisting of starburst cells with only two dendrites: one dendrite pointing to the left, the other dendrite pointing to the right (Fig. 1A). The two dendrites of a cell are modeled to be functionally independent of each other, i.e., there is no "diffusion" of activity between the two dendrites of a cell. In this one-dimensional case, there are two configurations for a pair of interacting dendrites: the two interacting dendrites can point in the same direction (like cells A and B in Fig. 1B) or they point in opposite directions (cells B and C).


Figure 1
View larger version (68K):
[in this window]
[in a new window]
 
FIG. 1. Layout of the model. A: geometric properties of individual starburst cells. Each dendrite receives inputs at its distal 110 µm and makes output synapses with its distal 50 µm. We chose to exclude input sites from the proximal 40 µm because we wanted to focus on interactions between starburst cells, and the proximal region of starburst dendrites does not lie in the "starburst stratum" of the inner plexiform layer (Famiglietti 1991Go). B, top: geometric properties of the starburst network. The amount of overlap between individual dendrites is determined by the geometry. For example, the release region of the left dendrite of cell A (LA) overlaps with the input region of LB by 20 µm, whereas, conversely, the release region of LB overlaps with the input region of LA by 50 µm. Overlapping starburst dendrites can point in the same direction (the dendrites of cells A and B) or in opposite directions (cells B and C). Starburst cells receive input from an array of bipolar cells (gray) spaced 30 µm apart. Bottom: 2 confocal micrographs showing 2 reconstructed starburst cells each. Left: cell pair is located closely together, therefore its dendrites point in the same direction. Right: cell pair is far apart, their overlapping dendrites point in opposite directions. Images courtesy of Professor David Vaney.

 
Our model is a discrete-time model. In each time step, we calculate the state of all starburst dendrites based on the state of all starburst dendrites at the previous time step and on the current "light stimulus", which is a user-given input to the model.

The state of each dendrite is represented by a number which can take on positive values (interpretation: "activated", "depolarized"), zero ("resting state"), or negative values ("suppressed", "hyperpolarized"). There are no boundaries built into the model as to how large or small this number can become. We chose not to impose absolute boundaries because any boundary would have to be set arbitrarily; interfering with the behavior of the model for some sets of parameters but not for others (see following text for a description of the parameters of the model). Compare, for example, the three graphs in Fig. 3B that show the model behavior for three sets of parameters. An arbitrary boundary of 100 would not affect the behavior in the first and third graph, whereas the cells in the second graph reach a maximum of ~130. However, it is worth pointing out that due to our restriction of the parameter space (see RESULTS), the cells in the model do behave bounded, despite the lack of an absolute imposed boundary. Although the maximum and minimum values that are reached depend on the particular set of parameters, it is therefore possible to interpret the state value as a (not necessarily linearly related) measure of membrane potential, or intracellular calcium, or synaptic release.


Figure 3
View larger version (38K):
[in this window]
[in a new window]
 
FIG. 3. Directional behavior of starburst model. A: we measured the response of the right dendrite of cell 28 and the left dendrite of cell 34 in a model with 61 starburst cells. The stimulus is a white bar (90 µm wide) on a gray background that moves from left to right at 30 µm/time step. B: example model behavior for three parameter combinations (A: d = 0.7, css = 0, cso = 0; B: d = 0.7, css = 3, cso = –1; C: d = 0.7, css = –4, cso = 3). During the first 20 time steps, the model dendrites reach steady state (like in the test for boundedness, Fig. 2A, top). For parameter combination A (no connectivity between starburst cells), R28 and L34 behave identical. For combination B, the dendrites prefer outward movement; for combination C they prefer inward movement. C: distribution of directional index of model behavior over the whole parameter space. Yellow-red areas correspond to preferred outward movement, greenish-blue areas to preferred inward movement. D: summary of model behavior. The directional preference of starburst dendrites is largely independent of the model parameter d.

 
Dendrites in the model are connected with rectifying synapses, i.e., a dendrite in the model can influence other dendrites only if its state is larger than zero. In addition, the effect that a particular dendrite can have on other dendrites is proportional to the length of the physical overlap between the release sites of that dendrite and the input sites of the receiving dendrites (Fig. 1B); in other words, we assume an even distribution of synapses along the contact sites. A trivial consequence of this model property is that nonoverlapping dendrites do not contact each other.

There are three free parameters in our model. The "decay factor" d determines how much of the activity of every dendrite decays from time step to time step. The second and third parameters (css and cso) reflect the synaptic connections between starburst cells, the "connection strength" (cs). cs can be set to positive values (reflecting excitatory interactions through acetylcholine) or negative values (reflecting inhibitory interactions through GABA) or it can be zero, indicating that there is no interaction. Relatively large positive or negative values of cs can be interpreted as high density of synapses along the dendrite and/or as high efficacy of the synapses. In our model, we allow the connection of two dendrites to depend on the direction in which they point. If two dendrites point in the same ("s") direction (like cells A and B in Fig. 1B), their synaptic connection strength is determined by the parameter css. If the dendrites point in opposite ("o") directions (like cells B and C), they interact synaptically with strength cso. The total connection strength between any two dendrites is then determined by the model parameter css or cso (depending on the relative direction of the dendrites), multiplied by the length of overlap between the output and receiving sites of the two dendrites. This geometric overlap, and therefore the degree of interaction between a pair of dendrites, is not necessarily symmetric (see Fig. 1B). We tested the directional behavior of the starburst cells in the model for all plausible combinations (see following text) of d, css, and cso.

The behavior of each dendrite in our model is therefore determined by the following equations (given for a left dendrite in the model)

Formula 1(1)

Formula 2(2)
A corresponding equation can be defined for the right dendrites of starburst cells in the model

Formula 3(3)

Formula 4(4)
with Li(t), Ri(t): state of the left an right dendrites of cell i at time t. Li and Ri are dimensionless numbers; N: number of starburst cells in the model; css, cso, d: parameters of the model, see above. The cs-parameters are given in mm–1; and

Formula 5(5)
overlapxy: length of overlap between input sites of dendrite x and output sites of dendrite y given in micrometer. These numbers are determined purely by the geometric layout of the model (Fig. 1B); the overlap between most dendritic pairs is 0 µm. To avoid boundary effects, the model is both large (large N) and circular, with cell 1 coming to lie next to cell N. Note also that the input to a starburst dendrite from other dendrites depends on their activity at the previous time step (before the decay of their activity determined by parameter d). light(t): stimulus that the model receives at time t. The light stimulus is a user-defined array of values between –1 (black) and +1 (white). The standard background intensity is 0 (gray). The stimulus is passed through a balanced spatial Mexican hat filter (difference of Gaussians, center sigma = 90 µm, surround sigma = 240 µm) before being passed on to an array of "bipolar cells". The output of the bipolar cells is a rectified version of the light stimulus, according to the equation output = 1/(1 + e1.54–7input), yielding no output to black (= –1) stimuli, some baseline release under background (= 0) conditions [corresponding to baseline glutamatergic input that starburst cells receive (Peters and Masland 1996Go)], and maximum release to white (= +1) stimuli. The bipolar cells are evenly spaced every 30 µm, and each starburst dendrite sums all activity of those bipolar cells that overlap with its input sites (Fig. 1B).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DESCRIPTION
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We used the following strategy to investigate the influence of network interactions on directional starburst cell behavior: first, we restricted the parameter space to reasonable values (see following text). Second, we scanned the complete remaining parameter space and quantified the directional behavior for each parameter combination. We then discuss the behavior of the model starburst cells and provide some intuition about the model properties that lead to directional behavior of our cells. We also compare our model to the behavior of real starburst cells, to predict actual connectivity within the starburst cell network, and the influence of that connectivity on the directional behavior of starburst cells.

Restriction of parameter space

The decay constant d has the trivial boundaries 0 (no decay between time steps) and 1 (complete decay). The parameters css and cso have no such natural boundaries, and the abstract and general design of the model does not provide us with boundaries dictated by the biology of the system. We found boundaries for parameters css and cso by applying two different criteria to the model behavior: boundedness and robustness.

Because our starburst model is a network of interacting elements, unbounded behavior is easily achieved. For example, if all dendrites in the model were to excite each other strongly, all dendrites would reinforce each other and the values Li(t) and Ri(t) would soon approach infinity. To test for boundedness, we simply let the system run with the light stimulus set to 0 (gray). This will lead to baseline release from the model bipolar cells; the starburst cells will interact and eventually either reach some steady state, or their behavior will be unbounded (Fig. 2A, top). Parameter combinations d–css–cso, which lead to unbounded behavior, were excluded.


Figure 2
View larger version (20K):
[in this window]
[in a new window]
 
FIG. 2. Restriction of the parameter space. A: example model behavior for 3 parameter combinations (A: d = 0.7, css = 3.5, cso = 2; B: d = 0.7, css = 2, cso = 0; C: d = 0.7, css = 3.5, cso = –3). Top: results of the test for boundedness; bottom: test results for robustness (see text for description of the tests). Only parameter combination B passes both tests, whereas A is not bounded and C is not robust. B: parameter values that are both bounded and robust are shaded gray. The location of these "reasonable" parameters within the css-cso plane depends on the value of d. Overall, the region of "reasonable" parameters is limited in size.

 
With a uniform background illumination, all dendrites in the model should in theory behave absolutely identically because, by design, the model is completely symmetric and homogeneous. In some cases, however, we observed that small numerical rounding errors caused very different behavior of left and right dendrites in the model (i.e., for those parameter combinations the system was not robust against rounding errors). To test for robustness, we therefore started with a system where all right dendrites in the model were set to a different initial value (= 0) than the left dendrites (= 1, but the specific value has no influence on the outcome of this test), whereas the light stimulus was still a uniform gray. Parameter combinations d–css–cso for which the difference between left and right dendrites diminished and approached 0 were considered robust. All other parameter combinations were excluded (Fig. 2A, bottom).

Parameter combinations d–css–cso that are both bounded and robust (according to the tests described in the preceding text) represent an upper bound of parameters that can be expected to reasonably describe "biological" behavior in our model. Those parameters lie within a convex region of the dcsscso space (Fig. 2B). We have therefore restricted the reasonable parameter space to a manageable size and are able to scan the complete space to observe the behavior of the starburst cells for each possible parameter combination.

Directional behavior of starburst network

We tested the directional behavior of the starburst cells with a white bar moving across the model network (Fig. 3A). We compared the behavior of the left and the right dendrites that are located in the center of the model (for n = 61 cells, this is R28 and L34). The bar is moving from left to right; therefore R28 sees outward or centrifugal movement, and L34 sees inward or centripetal movement. Figure 3B shows three individual examples for three different sets of parameters. The first pair of traces shows the behavior of R28 and L34 when css = cso = 0. In other words, there is no interaction within the starburst network; the starburst dendrites receive input exclusively from the bipolar cells. R28 and L34 do not behave differently in this case, consistent with the design of the model in which the starburst dendrites have no intrinsic directional behavior. The other two examples illustrate the two possible asymmetries that the cells in the model can express: the starburst dendrites can either respond more strongly to outward movement (R28 > L34, Fig. 3B, middle), or they can prefer inward movement (R28 < L34, Fig. 3B, bottom). We quantified the directional preference by calculating the directional index Formula 5. Values of DI > 1 indicate preferred outward movement, DI = 1 indicates nondirectional behavior, and DI < 1 indicates preferred inward movement. Figure 3C shows the distribution of the directional index over the parameter space. We find that the value of the decay factor d is not critical for the directional behavior of starburst dendrites in the model. As a general rule, whenever css is larger than cso, the dendrites prefer outward movement. When css is smaller than cso the dendrites respond stronger to inward movement (Fig. 3D).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DESCRIPTION
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
It has long been a hot topic of discussion whether directionally asymmetric behavior in the nervous system requires a neuron with intrinsic directional properties or if directional behavior can emerge as a circuit property, even if there are no intrinsically directional neurons in the circuit. The dendrites of starburst amacrine cells in the retina are one of the earliest building blocks in the circuit of DS cells that show directional behavior. Many studies have explained this behavior with intrinsic properties of starburst cells. In this report, we show that the directional behavior can emerge from network interactions, independent of intrinsic directional properties of starburst dendrites.

How does the model generate directional selectivity?

The starburst dendrites in our model do not behave asymmetrically unless there is some directional bias in the synaptic circuitry. For example, let us consider the right dendrite of cell B in Fig. 1B. It may receive inputs from dendrites that point in the same direction (e.g., the right dendrite of cell A) and from dendrites that point in opposite direction (e.g., the left dendrite of cell C). If these synaptic connections have the same strength (i.e., css = cso), then the starburst network will not impose directional behavior on the right dendrite of cell B (Fig. 3, C and D). If there is a bias, however, to favor excitation (or inhibition) for like-wise oriented dendrites, then this bias will result in a directional behavior of the starburst dendrite: it will either favor outward or inward movement, depending on the bias.

At first glance, it may seem surprising that the dendrites in our model can behave directionally because the detection of movement direction inherently requires the "analysis" of the change of position of an object over sequential time points. How can a dendrite in our model do this, although it is modeled as a single compartment, and does therefore not have the ability to distinguish between different spatial positions? The answer is that a dendrite cannot do this—the directional behavior is a network property (reflected in the behavior of the individual dendrites) and not the property of any single dendrite.

It is worth to provide some intuition about the network interactions that underlie the directional behavior. Let’s first consider the two dendrites A and B in Fig. 4A. They are at the same spatial location but point in opposite directions. We now show a brief stimulus (lasting only 1 time step) that lies just to the left of both dendrites (solid outline, Fig. 4A), so that it causes no additional input to either dendrite A or B beyond the steady-state background input that all dendrites receive. The stimulus will, however, excite the three gray dendrites shown above dendrite A. At the next time step, the stimulus disappears. Dendrites A and B still receive the steady-state background input and, in addition, input from the three gray dendrites that can be either positive, negative, or zero depending on the values of css and cso. For simplicity, we set cso = 0 (so that dendrite B will receive no input from the gray dendrites), and consider the cases css > 0, which leads to positive input to dendrite A (and hence we will get A > B) or css < 0, which leads to negative input to dendrite A (and A < B). In each case, dendrites A and B have now different values. Likewise, if we set css = 0, we can consider similar scenarios for positive and negative values of cso, and will again get different activity of dendrites A and B.


Figure 4
View larger version (32K):
[in this window]
[in a new window]
 
FIG. 4. Principle underlying the generation of directionality in the model. A: illustration of principles underlying network interactions. See text for details. B: stimulus used in C. See Fig. 3A for details. C: responses of R28 (red) and L34 (blue) to a bar moving slowly to the right (moving 1 spatial unit every 12 time steps, indicated by gray vertical lines). Parameters: d = 1, css = 5, cso = 0. Left inset: responses to the same stimulus moving at every time step. Right inset: integrated responses over time of R28 and L34 to the slow stimulus.

 
This example shows the underlying basis for the directional behavior of the model. Even a brief stationary stimulus sets up different activity patterns in the left- and right-pointing dendrites surrounding the stimulus location (as long as css != cso), which will lead to directional behavior for a moving stimulus. The crucial property of the model is therefore the ability of starburst cells to influence other cells that are displaced from the stimulus location; they can perform "action at a distance". By this mechanism, the dendrites A and B of the preceding example integrate information about the activity at a different location at the previous time step (mediated by the gray dendrites).

How does this concept (action at a distance with a delay of 1 time step) lead to directionally selective behavior when the stimulus is actually moving? Figure 4, B and C, shows a specific example that illustrates the underlying interactions. For simplicity, this example is modeled with no outer retinal preprocessing of the stimulus (as opposed to the traces shown in Fig. 3) with cso set to 0 and positive css (css = 5). For illustrative purposes, we set d = 1 (complete decay between time steps, so that the state value of a dendrite reflects only its current inputs and is not directly dependent on its own past). The 17 positions of the rightward moving stimulus bar (positions –8 to +8) are shown in Fig. 4B, top (these are the same positions as for the stimulus used in Fig. 3), and C shows the responses of the central dendrites R28 and L34. We use a slowly moving stimulus [stepping 1 spatial unit (= 30 µm) every 12 time steps, indicated by the gray vertical lines] that clearly reveals the interactions underlying the directional behavior. For comparison, Fig. 4C, left inset, shows the responses when the stimulus moves at the same speed as in Fig. 3, i.e., stepping one spatial unit each time step.

Note the transient peaks in the response of R28 during the first half of the movement (up to position 0). These peaks are signatures of coincidence detection that is happening in the network. Consider for example the transient peak as the stimulus steps from position –1 to position 0. This peak is caused by an increase in the bipolar input to R28 because the stimulus has just moved to the right and is covering more of the input region of dendrite R28. At the same time, R28 is still receiving strong input from the likewise oriented dendrites (compare the gray dendrites in Fig. 4A) where the stimulus has been at the previous time step. The decay after the peak is due to decreased input to R28 from those likewise-oriented dendrites after the stimulus has stepped away. Importantly, this decreased network input "reaches" R28 one time step delayed compared with the increased bipolar input. In other words, the transient peak is due the coincidence of strong input from two sources: the increased direct bipolar input and the not-yet decreased input from the likewise oriented neighboring dendrites. Each time the stimulus steps to the right, it causes such a transient peak in the dendrites directly underlying the stimulus, as described in the preceding text. But even dendrites further away from the stimulus will eventually "see" this event because of successive network transmission through the likewise-oriented dendrites. For example, dendrite R28 eventually reports the step of the stimulus from position –8 to position –7, even though this event is happening far to the left of the input sites of R28.

The transient "negative dips" in the response of L34 in the second half of the movement (after position 0) are not due to any negative connections. Instead, those dips correspond to the stimulus moving away from the left-pointing dendrite L34, so that it gets (instantaneously) less bipolar input. Then the response grows again due to inputs from other left-pointing dendrites located further to the right (network effect that is delayed by 1 time step). Again, such an event is transmitted to L34 through likewise oriented dendrites even if it happens further to the right of L34.

A fast-moving stimulus, like shown in the left inset, will emphasize the transient aspects of the responses (like the peaks in the response of R28), so that the directional difference is more pronounced than during the slowly moving stimulus. In fact, the steady-state values of R28 and L34 (after a stimulus has been at a position for a while) are not different, as is most easily seen at position 0, but also at any other pair of corresponding symmetric positions (e.g., the responses of R28 at position –2 and of L34 at position +2: the steady-state responses are the same, even though the stimulus gets to both positions from the left, i.e., outward for R28 and inward for L34).

It can be shown mathematically that there are two general requirements for the discrimination of movement direction (see for example the classical model of directional selectivity of Barlow and Levick 1965Go). One requirement is the comparison of the activity at two spatial locations ({Delta}s) at two different time points ({Delta}t). This requirement is fulfilled in our model by the "action at a distance" of starburst cells ({Delta}s) that they perform with the delay of one time step ({Delta}t) as described in the preceding text. Physiologically, this delay can be interpreted as the delay and/or persistence of synaptic release. The second requirement is at least one nonlinearity in the analysis. The only nonlinearity in our model is the rectification of release from starburst cells, which only kicks in for negative values of the state value. For most parameter combinations in our model, however, the starburst cells stay in the positive range, that is, there is no nonlinearity involved. How can the cells in our model still distinguish the direction of movement? The answer is that there is a hidden nonlinearity in the way we quantify the response: we compare the maximum of the of the cells’ responses, which is equivalent to a thresholding operation. If we were to compare the total activity of the cells, i.e., the "area under the curve" of the responses, we would find no difference between inward and outward movement. Fig. 4C, right inset, illustrates this linearity. Note also that the areas under the curves shown in Fig. 3B are the same for the responses of R28 and L34.

It should also be mentioned that the Mexican hat filter of the "outer retina" is not crucial for the directional behavior of the model. If we run the simulation with no outer retinal preprocessing, as exemplified by Fig. 4, the results are very similar. The same is true when the bar moves at twice the speed as shown in Fig. 3 (data not shown).

What does the model tell us about the biology of starburst cells?

We know from experimental results that starburst dendrites do indeed respond more strongly to outward movement (Euler et al. 2002Go). The strongest prediction of our model is therefore that there should be no connectivity bias that would favor inward movement (i.e., we should not find cso > css) because then, even if starburst cells do have intrinsic properties, these properties would have to be strong enough to overcome such a bias imposed by the network. It seems unlikely to us that network and internal properties compete with each other.

What is known about connectivity of starburst cells in the retina? Zhou and colleagues (Zheng et al. 2004Go) recently reported that in the maturing rabbit retina, cholinergic nicotinic synapses between starburst cells slowly disappear, whereas GABAergic connections remain present. In the terminology of our model, this means that both css and cso are nonpositive in the mature retina. Negative cso corresponds to an inhibitory surround of the starburst cell (note that cell C in Fig. 1B is in the position to supply surround-input to cell B and vice versa). The question then is whether cso is more negative than css, or, in other words, if any given starburst cell receives stronger inhibitory input from the starburst cells with cell bodies lying in its surround ("surround-cells") than from those starburst cells with cell bodies that lie within its dendritic field ("nonsurround cells"). If this is the case, our model predicts that we can attribute at least some of the directional behavior of starburst dendrites to the retinal circuitry. Zhou and colleagues (Lee and Zhou 2005Go) also performed double patch experiments and found inhibitory connections between starburst cells that were close together (non-surround cells, like cells A and B in Fig. 1B) and also between cells that were far apart (surround-cells, like cells B and C). Unfortunately, these experiments do not allow us to strictly answer the question if cso < css because one can only measure the overall synaptic input to the cells and not the strength of the synaptic input to an individual dendrite. It is likely that a cell pair A–B, which is spatially almost completely co-incident, has more synapses in common than a cell pair B–C, which has a much smaller region of overlap. This would result in stronger inhibitory currents measured in a double-patch experiment between cells A and B than between cells B and C (which could be interpreted as cso > css). If we were able to measure the inhibitory input to an individual dendrite, the result may or may not be opposite.

It was shown that the directional behavior of starburst dendrites remains intact in the presence of blockers of GABAA receptors (Euler et al. 2002Go) as well as GABAC receptors (Hausselt et al. 2004Go). Based on these findings and those of Zhou and colleagues (that there are no nicotinic synapses between starburst cells in the adult retina), one might conclude that network interactions cannot underlie the directional behavior of starburst cells. If we allow only direct and electrogenic synapses between starburst cells (using nicotinic acetylcholine receptors, and GABAA and GABAC receptors), this seems indeed to be the case. The conclusion of our model would then be: the directional behavior of starburst cells has to be, at least in part, due to internal properties of the dendrites.

However, we cannot exclude the possibility that other types of connectivity exist in the starburst network; even indirect connections between starburst cells (through an additional interneuron) have been proposed in the literature. There are therefore other possible mechanisms of how one starburst dendrite can have a "positive" or "negative" influence on another: a positive value of cs in our model could be interpreted as cholinergic enhancement of glutamate release from bipolar cells (Yamada et al. 2003Go), as inhibition of an inhibitory amacrine cell (disinhibition), or as a nonelectrogenic muscarinic excitation between starburst dendrites (through receptors other than M2) (see Wasselius et al. 1998Go). None of these indirect or nonelectrogenic connections could be easily detected by common electrophysiological techniques. Similarly, one starburst cell could have "negative" influence on another starburst cell through nonelectrogenic connections through GABAB receptors (Zucker et al. 2005), the reduction of glutamate release from bipolar cells (Linn and Massey 1992Go), or cholinergic excitation of a presynaptic GABAergic or glycinergic amacrine cell (indirect inhibition) (Neal and Cunningham 1995Go). Again, these indirect connections are not easily observed by standard electrophysiological techniques.

In summary, we conclude that interactions within a homogeneous network of starburst cells can impose a directional bias on (intrinsically nondirectional) starburst cell dendrites. If we only consider direct synaptic connections between starburst cells, such a directional bias can in principle be generated by a symmetrical inhibition from surrounding starburst cells with processes that point toward the central starburst cell. Previous experimental results suggest that this could not be the only mechanism that generates directional behavior in starburst dendrites. But the influence of the circuitry does not have to be confined to direct connections between starburst cells. This would allow for even more sophisticated directional network processing.


    GRANTS
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DESCRIPTION
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
This work was supported by National Eye Institute Grant EY-00561.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DESCRIPTION
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank D. Vaney for supplying the confocal micrographs of starburst cells used in Fig. 1B and N. Ghorashian for help with the initial phase of the model design.

Present address of T. A. Münch: Friedrich Miescher Institute for Biomedical Research, Maulbeerstr. 66,4058 Basel, Switzerland.


    FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Address for reprint requests and other correspondence: F. S. Werblin, Molecular and Cell Biology, University of California at Berkeley, 145 Life Sciences Addition, Berkeley, CA 94720 (E-mail:werblin{at}berkeley.edu)


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MODEL DESCRIPTION
 RESULTS
 DISCUSSION
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Barlow H. Intraneuronal information processing, directional selectivity and memory for spatio-temporal sequences. Network Comput Neural Syst 7: 251–259, 1996.[CrossRef]

Barlow HB and Levick WR. Mechanism of directionally selective units in rabbits retina. J Physiol 178: 477–504, 1965.[Free Full Text]

Borg-Graham L and Grzywacz NM. A model of the direction selectivity circuit in retina: transformations by neurons singly and in concert. In: Single Neuron Computation, edited by T. McKenna, J. Davis, and S. F. Zornetzer. New York: Academic, 1992, p. 347–376.

Brecha N, Johnson D, Peichl L, and Wässle H. Cholinergic amacrine cells of the rabbit retina contain glutamate-decarboxylase and gamma-aminobutyrate immunoreactivity. Proc Natl Acad Sci USA 85: 6187–6191, 1988.[Abstract/Free Full Text]

Euler T, Detwiler PB, and Denk W. Directionally selective calcium signals in dendrites of starburst amacrine cells. Nature 418: 845–852, 2002.[CrossRef][Medline]

Famiglietti EV. On and off pathways through amacrine cells in mammalian retina—the synaptic connections of starburst amacrine cells. Vision Res 23: 1265–1279, 1983.[CrossRef][ISI][Medline]

Famiglietti EV. Synaptic organization of starburst amacrine cells in rabbit retina—analysis of serial thin sections by electron-microscopy and graphic reconstruction. J Comp Neurol 309: 40–70, 1991.[CrossRef][ISI][Medline]

Feller MB. The role of nAChR-mediated spontaneous retinal activity in visual system development. J Neurobiol 53: 556–567, 2002.[CrossRef][ISI][Medline]

Fried SI, Münch TA, and Werblin FS. Mechanisms and circuitry underlying directional selectivity in the retina. Nature 420: 411–414, 2002.[CrossRef][Medline]

Fried SI, Münch TA, and Werblin FS. Directional selectivity is formed at multiple levels by laterally offset inhibition in the rabbit retina. Neuron 46: 117–127, 2005.[CrossRef][ISI][Medline]

Gavrikov KE, Dmitriev AV, Keyser KT, and Mangel SC. Cation-chloride cotransporters mediate neural computation in the retina. Proc Natl Acad Sci USA 100: 16047–16052, 2003.[Abstract/Free Full Text]

Hausselt SE, Detwiler PB, Fantana H, Denk W, and Euler T. Computation of directional motion in starburst amacrine cells dendrites—a role for active currents. Soc Neurosci Abstr 299.4, 2004.

Lee S and Zhou ZJ. Directionally asymmetric synaptic inputs to starburst amacrine cell dendrites during light stimulation (E-Abstract). Invest Ophthalmol Vis Sci 46: 2337, 2005.

Linn DM and Massey SC. GABA inhibits ACh release from the rabbit retina: a direct effect or feedback to bipolar cells? Vis Neurosci 8: 97–106, 1992.[ISI][Medline]

Neal MJ and Cunningham JR. Baclofen enhancement of acetylcholine release from amacrine cells in the rabbit retina by reduction of glycinergic inhibition. J Physiol 483: 363–372, 1995.

O’Malley DM, Sandell JH, and Masland RH. Corelease of acetylcholine and GABA by the starburst amacrine cells. J Neurosci 12: 1394–1408, 1992.[Abstract]

Peters BN and Masland RH. Responses to light of starburst amacrine cells. J Neurophysiol 75: 469–480, 1996.[Abstract/Free Full Text]

Tukker JJ, Taylor WR, and Smith RG. Direction selectivity in a model of the starburst amacrine cell. Vis Neurosci 21: 611–625, 2004.[Medline]

Vaney DI. Coronate amacrine cells in the rabbit retina have the starburst dendritic morphology. Proc R Soc Lond B Biol Sci 220: 501–508, 1984.

Wasselius J, Johansson K, Bruun A, Zucker C, and Ehringer B. Correlations between cholinergic neurons and muscarinic m2 receptors in the rat retina. Neuroreport 9: 1799–1802, 1998.[ISI][Medline]

Yamada ES, Dmitrieva N, Keyser KT, Lindstrom JM, Hersh LB, and Marshak DW. Synaptic connections of starburst amacrine cells and localization of acetylcholine receptors in primate retinas. J Comp Neurol 461: 79–90, 2003.

Yoshida K, Watanabe D, Ishikane H, Tachibana M, Pastan I, and Nakanishi S. A key role of starburst amacrine cells in originating retinal directional selectivity and optokinetic eye movement. Neuron 30: 771–780, 2001.[CrossRef][ISI][Medline]

Zheng JJ, Lee S, and Zhou ZJ. A developmental switch in the excitability and function of the starburst network in the mammalian retina. Neuron 44: 851–864, 2004.[CrossRef][ISI][Medline]

Zhou ZJ and Fain GL. Neurotransmitter receptors of starburst amacrine cells in rabbit retinal slices. J Neurosci 15: 5334–5345, 1995.[Abstract]

Zucker CL, Ehinger B, and Grzywacz NM. Compartmental localization of GABAB receptors on starburst amacrine cell dendrites in the rabbit retina (E-Abstract). Invest Ophthalmol Vis Sci 45: 25362, 2005.





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
96/1/471    most recent
00628.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via ISI Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Münch, T. A.
Right arrow Articles by Werblin, F. S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Münch, T. A.
Right arrow Articles by Werblin, F. S.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2006 by the The American Physiological Society.