|
|
||||||||
REPORT
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 |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
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 1992
; Tukker et al. 2004
); asymmetric distribution of chloride co-transporters along the dendrites that lead to spatially asymmetric chloride currents and hence to directional asymmetries (Gavrikov et al. 2003
); or cell-internal biochemical processes, calcium-induced calcium currents, which lead to stronger calcium signals for outward movement (Barlow 1996
). 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 1991
). Second, the dendritic trees of starburst cells overlap strongly (Famiglietti 1983
; Vaney 1984
). 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. 1988
; OMalley et al. 1992
) and also have receptors for both of these neurotransmitters (Feller 2002
; Zhou and Fain 1995
).
| MODEL DESCRIPTION |
|---|
|
|
|---|
|
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.
|
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)
![]() | (1) |
![]() | (2) |
![]() | (3) |
![]() | (4) |
![]() | (5) |
| RESULTS |
|---|
|
|
|---|
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 dcsscso, which lead to unbounded behavior, were excluded.
|
Parameter combinations dcsscso 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
. 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 |
|---|
|
|
|---|
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 thisthe 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. Lets 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.
|
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 1965
). One requirement is the comparison of the activity at two spatial locations (
s) at two different time points (
t). This requirement is fulfilled in our model by the "action at a distance" of starburst cells (
s) that they perform with the delay of one time step (
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. 2002
). 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. 2004
) 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 2005
) 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 AB, which is spatially almost completely co-incident, has more synapses in common than a cell pair BC, 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. 2002
) as well as GABAC receptors (Hausselt et al. 2004
). 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. 2003
), 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. 1998
). 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 1992
), or cholinergic excitation of a presynaptic GABAergic or glycinergic amacrine cell (indirect inhibition) (Neal and Cunningham 1995
). 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 |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
Present address of T. A. Münch: Friedrich Miescher Institute for Biomedical Research, Maulbeerstr. 66,4058 Basel, Switzerland.
| FOOTNOTES |
|---|
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 |
|---|
|
|
|---|
Barlow HB and Levick WR. Mechanism of directionally selective units in rabbits retina. J Physiol 178: 477504, 1965.
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. 347376.
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: 61876191, 1988.
Euler T, Detwiler PB, and Denk W. Directionally selective calcium signals in dendrites of starburst amacrine cells. Nature 418: 845852, 2002.[CrossRef][Medline]
Famiglietti EV. On and off pathways through amacrine cells in mammalian retinathe synaptic connections of starburst amacrine cells. Vision Res 23: 12651279, 1983.[CrossRef][ISI][Medline]
Famiglietti EV. Synaptic organization of starburst amacrine cells in rabbit retinaanalysis of serial thin sections by electron-microscopy and graphic reconstruction. J Comp Neurol 309: 4070, 1991.[CrossRef][ISI][Medline]
Feller MB. The role of nAChR-mediated spontaneous retinal activity in visual system development. J Neurobiol 53: 556567, 2002.[CrossRef][ISI][Medline]
Fried SI, Münch TA, and Werblin FS. Mechanisms and circuitry underlying directional selectivity in the retina. Nature 420: 411414, 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: 117127, 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: 1604716052, 2003.
Hausselt SE, Detwiler PB, Fantana H, Denk W, and Euler T. Computation of directional motion in starburst amacrine cells dendritesa 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: 97106, 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: 363372, 1995.
OMalley DM, Sandell JH, and Masland RH. Corelease of acetylcholine and GABA by the starburst amacrine cells. J Neurosci 12: 13941408, 1992.[Abstract]
Peters BN and Masland RH. Responses to light of starburst amacrine cells. J Neurophysiol 75: 469480, 1996.
Tukker JJ, Taylor WR, and Smith RG. Direction selectivity in a model of the starburst amacrine cell. Vis Neurosci 21: 611625, 2004.[Medline]
Vaney DI. Coronate amacrine cells in the rabbit retina have the starburst dendritic morphology. Proc R Soc Lond B Biol Sci 220: 501508, 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: 17991802, 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: 7990, 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: 771780, 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: 851864, 2004.[CrossRef][ISI][Medline]
Zhou ZJ and Fain GL. Neurotransmitter receptors of starburst amacrine cells in rabbit retinal slices. J Neurosci 15: 53345345, 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.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |