Abstract
The lobula giant motion detector (LGMD) is a wide-field bilateral visual interneuron in North American locusts that acts as an angular threshold detector during the approach of a solid square along a trajectory perpendicular to the long axis of the animal (Gabbiani et al., 1999a). We investigated the dependence of this angular threshold computation on several stimulus parameters that alter the spatial and temporal activation patterns of inputs onto the dendritic tree of the LGMD, across three locust species. The same angular threshold computation was implemented by LGMD in all three species. The angular threshold computation was invariant to changes in target shape (from solid squares to solid discs) and to changes in target texture (checkerboard and concentric patterns). Finally, the angular threshold computation did not depend on object approach angle, over at least 135° in the horizontal plane. A two-dimensional model of the responses of the LGMD based on linear summation of motion-related excitatory and size-dependent inhibitory inputs successfully reproduced the experimental results for squares and discs approaching perpendicular to the long axis of the animal. Linear summation, however, was unable to account for invariance to object texture or approach angle. These results indicate that LGMD is a reliable neuron with which to study the biophysical mechanisms underlying the generation of complex but invariant visual responses by dendritic integration. They also suggest that invariance arises in part from non-linear integration of excitatory inputs within the dendritic tree of the LGMD.
Invariance of neuronal responses is a key aspect of sensory processing. It characterizes the extraction of specific features in a stimulus, independent of their context (Cavanagh, 1978). Some invariant responses, such as those underlying the detection of edges in an image regardless of contrast level, rely on gain control mechanisms that are relatively well understood at the cellular and network levels (Galarreta and Hestrin, 1998; Sanchez-Vives et al., 2000a,b). In the inferotemporal visual area (IT) of the monkey, evidence suggests that the responses of many neurons are invariant to translation and scale transformations (Schwartz et al., 1983; Rolls and Baylis, 1986), a possible prerequisite for object identification (Sary et al., 1993). Although little is known about how this type of invariance arises, some authors have recognized the need for mechanistic investigations and proposed models (Olshausen et al., 1993; Salinas and Abbott, 1997;Riesenhuber and Poggio, 1999). The neural basis of response invariance, however, remains poorly understood. It is not known whether it results from network properties, from dendritic processing within single cells, or from combinations of both. We study two monosynaptically connected neurons in the visual system of the locust and demonstrate that they exhibit invariance properties to many attributes of the stimuli that best excite them. Because the neuronal processing leading to these invariant responses is thought to occur within the dendritic tree of one of these neurons, it might prove appropriate to study the biophysical mechanisms underlying the generation of invariant visual responses in single neurons.
The lobula giant motion detector (LGMD) is located in the third visual neuropil of the locust optic lobe (Strausfeld and Nässel, 1981). Its dendritic tree consists of three subfields that arborize in the lobula (O'Shea and Williams, 1974). The largest subfield receives excitatory inputs from afferents that are sensitive to local motion over the whole visual hemifield, whereas the remaining two subfields are thought to receive feedforward inhibitory projections that are size dependent (Palka, 1967; Rowell et al., 1977). LGMD synapses onto the descending contralateral motion detector (DCMD) (Killmann et al., 1999), a neuron the large axon of which projects to thoracic motor centers responsible for the generation of jump and flight steering maneuvers (Pearson et al., 1980; Simmons, 1980; Robertson and Pearson, 1983). The synaptic connection between LGMD and DCMD is so strong as to cause a one-to-one correspondence between presynaptic and postsynaptic spiking activity under visual stimulation (O'Shea and Williams, 1974; Rind, 1984). LGMD and DCMD are vigorously excited by objects looming toward the animal (Schlotterer, 1977; Rind and Simmons, 1992). The peak firing rate during approach of such looming objects signals the moment when the object reaches a fixed angular threshold size (Gabbiani et al., 1999a). This stimulus variable has in turn been related to the generation of escape and collision avoidance behaviors (Robertson and Johnson, 1993; Hatsopoulos et al., 1995).
Invariance of the peak firing time of LGMD and DCMD to changes of body temperature and stimulus luminance/contrast have already been demonstrated (Gabbiani et al., 1999a). In the present work, we investigate the extent to which the peak firing time of LGMD/DCMD—and thus the angular threshold detection computation—changes as we vary parameters of the approaching object that are expected to directly affect the time course of the spatially distributed excitatory and inhibitory inputs impinging onto the dendrites of the LGMD.
MATERIALS AND METHODS
The experimental materials and methods used were similar to those of Gabbiani et al. (1999a). The following brief account mainly emphasizes the differences with that work.
Preparation. Experiments were performed on three different species of North American and African locusts. Male and female specimens of the species Schistocerca americana (American Bird Grasshopper) were taken from the laboratory colony 3–4 weeks after their final molt. Adult male specimens from the African speciesSchistocerca gregaria (Desert Locust) and Locusta migratoria (African Migratory Locust) were imported (USDA–APHIS permits 36933 and 42816, respectively) from the International Center for Insect Physiology and Ecology (ICIPE, Nairobi, Kenya) and the University of Bielefeld (Germany). After the animals were fixed to a plastic holder, the head was aligned under a microscope with reference points marked on a reticular grid inserted in one of the eyepieces. This procedure and the calibration protocol described inGabbiani et al. (1999a) allowed us to reliably align the center of the locust eye with the center of the video monitor used for visual stimulation. Animals were prepared for electrophysiological recordings in one of the following three manners: (1) mounted dorsal side up on the plastic holder, frontal dissection of the head capsule, performed as in Gabbiani et al. (1999a); (2) mounted ventral side up, connectives exposed by an incision of the neck cuticle (Hatsopoulos et al., 1995); or (3) mounted dorsal side up, rectangular incision of the pronotum. The gut was grabbed with a pair of fine forceps, cut as close as possible to the mouth, and removed through an abdominal incision, thus exposing the connectives. In all cases, the preparation was bathed in locust saline (Laurent and Davidowitz, 1994).
Electrophysiology and data acquisition. Locusts were fixed to a clamp with their longitudinal body axis parallel to the stimulation screen, except for the turntable experiments described below. For preparation 1, the connective contralateral to the stimulated eye was placed in a suction electrode (Gabbiani et al., 1999a). Alternatively, in preparations 2 and 3, two hook electrodes made of 50 μm (0.02 inch) stainless steel wire isolated up to the tip with H-Formvar (California Fine Wire Co., Grover City, CA) were placed around the connective. The two electrodes were electrically isolated from each other with Vaseline (Hatsopoulos et al., 1995). Extracellular signals were amplified with a differential AC amplifier (A-M Systems, Everett, WA) and a Brownlee amplifier (model 210A; Brownlee, San Jose, CA). They were acquired together with transistor transistor logic (TTL) pulses synchronizing the visual stimuli with the recordings (Gabbiani et al., 1999a) using a 12 bit A/D board (win30; United Electronic Industries, Watertown, MA) connected to a personal computer running the QNX real-time operating system (QSSL, Kanata, Ontario, Canada). Each recorded spike train was inspected visually; DCMD action potentials (typically the largest in the nerve cord) (Fig.1A) were selected on-line using custom software written in C.
Stimulation. Stimuli were generated using a fast monochrome video monitor refreshed at a rate of 200 Hz, as described inGabbiani et al. (1999a). The sequences of video images simulated circular and square objects with various textures approaching on a collision course with the animal. Each image was computed by central projection from the center of the stimulated eye onto the screen plane (Gabbiani et al., 1999a, their Fig. 2A and Eq. 2). The distance eye-screen was 120 mm. The parameters characterizing each object approach were the radius (half-size), l, of the approaching disk or square (in centimeters), and the simulated constant speed of approach, v (in meters per second) (see Fig. 1A, inset). If we set t = 0 as the time of expected collision and adopt the convention that t < 0 before collision, then the velocity v is also negative when the object is approaching (Gabbiani et al., 1999a, their Eq. 1). Because the objects approached at a constant speed, elementary trigonometry shows that the image sequence is characterized by the ratio of these two parameters l/‖v‖ (in milliseconds) (Gabbiani et al., 1999a, their Fig. 1A). The range of values used for l/‖v‖ was from 5 to 50 msec (see species, shape, texture, and approach angle protocols below), corresponding to approach speeds between 2 and 10 m/sec for half-sizes between 6 and 14 cm (Gabbiani et al., 1999a, their Table 1).
Species protocols. In the first series of experiments, we compared the responses across different locust species by presenting to each animal 10 repetitions of looming squares approaching at 5 values of l/‖v‖ pseudorandomly interleaved (from l/‖v‖ = 10 –50 msec in steps of 10 msec; 50 protocols total). An interval of 40 sec was inserted between successive stimulus presentations to minimize habituation (Gabbiani et al., 1999a). The squares were black (0 cd/m2) on a bright background (95 cd/m2). This protocol was applied to n = 5 specimens of the species S. gregaria, 3 L. migratoria, and 64 S. americana. An additional eight L. migratoria were stimulated 10 times with 10 different values of l/‖v‖ (from 5 to 50 msec in steps of 5 msec) interleaved pseudorandomly (100 protocols total; an experiment lasted ∼1.25 hr). This protocol was identical to the one used in Gabbiani et al. (1999a) and allowed for a more stringent test of the linear relationship between peak firing time and l/‖v‖, because an additional 50 trials at 5 supplementary l/‖v‖ values were gathered to test the linear model.
Shape protocol. In the second series of experiments we compared the responses to black disks and black squares looming toward the animal by presenting 10 times these two objects at 5 values of l/‖v‖ (l/‖v‖ = 10 –50 msec in steps of 10 msec; two shapes × 5 l/‖v‖ values × 10 repetitions = 100 protocols; n = 5 animals,S. americana). Both the shape of the presented object and the value of l/‖v‖ were interleaved pseudorandomly, with a 40 sec interstimulus interval.
Texture protocols. In the third series of experiments, we studied the effect of target texture by comparing the responses to a black looming square and a square textured with a 3 × 3 checkerboard pattern as illustrated in Figure 5A (inset) (l/‖v‖ = 10 –50 msec in steps of 10 msec; 100 pseudorandomly interleaved protocols; 40 sec interstimulus interval; n = 12 animals, 5 S. gregaria, 7S. americana). The background was bright with a luminance IB = B·Imax (B = 100%, Imax = 95 cd/m2), and the checkerboard pattern had five squares of luminance IO = O·Imax, and four squares of luminance IP = P·Imax with O = 0% and P = 30%. In two experiments, the contrast between the squares of the checkerboard pattern was increased by setting P = 70%.
The fourth series of experiments was a modified version of the previous series, with a textured pattern consisting of four concentric squares as illustrated in Figure 6A (inset). The luminance of the background, IB, and object, IO, IP, were identical to those of the previous series (B = 100%, O = 0%, P = 30%; 100 pseudorandomly interleaved protocols; 40 sec interstimulus interval; n = 12 animals, 9 S. americana, 3 L. migratoria).
Approach angle protocols. The last three series of experiments investigated the effect of target approach angle. Animals were fixed to a clamp attached to a turntable that allowed rotation of both the animal and the micromanipulators holding the recording electrodes around a virtual vertical axis passing through the center of the eye. The turntable was equipped with a graduated dial allowing angles to be read with an accuracy of ±1°. We denote by 0° the position for which the longitudinal body axis is perpendicular to the stimulation screen (i.e., the animal faces the screen; see Fig.7A, inset). Positive angles denote counterclockwise rotation, and negative angles denote clockwise rotations. The turntable could be positioned from 135° (the object approaches from the same side as the stimulated eye, almost from the back) to −45° (the object approaches from the side contralateral to the stimulated eye). These experiments were performed using only the third dissection described above (see Preparation).
In the fifth series of experiments, the effect of approaching angle was tested at two positions: 0° and 90°. During the first six experiments, the angles were alternated between 0° and 90°, and the value of l/‖v‖ (l/‖v‖ = 10 –50 msec in steps of 10 msec) were varied pseudorandomly from trial to trial (2 positions × 5 l/‖v‖ values × 10 repetitions = 100 protocols). During the next seven experiments, the angle was fixed at 0° for the first 50 trials (pseudorandomly interleaved values of l/‖v‖) and then at 90° for the remaining 50 trials (n = 13 animals total, S. americana; 40 sec interstimulus interval).
In the sixth series of experiments, the effect of approaching angle was tested at four positions on the same side as the stimulated eye (0°, 45°, 90°, and 135°, from front to back) and 3 values of l/‖v‖ (l/‖v‖ = 10, 30, and 50 msec). At the end of a trial the animal was moved from one position to the next, and a new value of l/‖v‖ was chosen pseudorandomly. Each combination of position and l/‖v‖ value was repeated 10 times (4 positions × 3 l/‖v‖ values × 10 repetitions = 120 protocols; n = 5 animals,S. americana; 40 sec interstimulus interval).
In the seventh series of experiments, five positions close to 0° (animal facing the screen) were investigated (−45°, −22.5°, 0°, 22.5°, and 45°) at 3 values of l/‖v‖ (l/‖v‖ = 10, 30, and 50 msec). At the end of a trial, the animal was moved from one position to the next, and a new value of l/‖v‖ was chosen pseudorandomly. Each combination of position and l/‖v‖ value was repeated eight times (5 positions × 3 l/‖v‖ values × 8 repetitions = 120 protocols; n = 11 animals, S. americana; 40 sec interstimulus interval). In this experiment, recordings were obtained from both connectives by placing one hook electrode around each connective and amplifying the voltage signal of each channel with respect to a far reference electrode implanted in the body.
Data analysis. Data analysis methods followed closely those of Gabbiani et al. (1999a). The extracellular traces obtained in response to the various stimuli described above were thresholded to extract the spike occurrence times of the DCMD (Fig.1A). Each spike train was smoothed with a 20 msec Gaussian window to obtain an estimate of the instantaneous firing rate (Fig.1B) (Gabbiani et al., 1999a). Estimates of the mean time of peak firing rate and its SD were obtained from the repetitions (from 8 to 10) of each trial (Fig. 1B). The timing of the peak consistently shifted closer to collision time as the parameter l/‖v‖ characterizing the stimulus was decreased (Fig. 1C). This dependence was studied by plotting the peak time as a function of l/‖v‖ (see Figs. 3-7). Linear regressions of peak firing time as a function of l/‖v‖ were performed for each experiment separately (see Figs. 3-7) to obtain the slope, α, and the intercept parameter, δ, of the best linear fit together with their SD (Press et al., 1992). The angular threshold angles (see Figs. 3-7) and their SD were computed from α using Equation 6 of Gabbiani et al. (1999a): Equation 1and by error propagation, respectively (i.e., using the formula Δf = ‖df/dx‖·Δx, for f = f(x)) (Bevington and Robinson, 1992, sect. 3.2). The delay between the time when the angular threshold size is reached and the peak firing rate, δ (see Figs. 3-7), is equal to the intercept parameter obtained by the linear fit described above. Both mean threshold angles and mean delay values are illustrated in Figures4-7 with their SDs. The SDs correspond to 68.3% confidence intervals on the mean values of these parameters, whereas two SDs correspond to 95.4% confidence intervals (Press et al., 1992, sect. 15.6). All data analysis was performed using Matlab 5.2 (The MathWorks, Natick, MA).
Two-dimensional modeling of the responses of the LGMD. The time course of LGMD/DCMD firing rate in response to solid squares or discs approaching toward the center of the eye along a trajectory perpendicular to the body axis may be described by the following one-dimensional model: Equation 2In this equation, (g ○ exp)(·) = g (exp(·)) represents the composition of the exponential function and a static nonlinearity g(·) described inGabbiani et al. (1999a; their Eq. 11 and Discussion). The term log ψ(t − δ), where ψ(t) = ½θ˙(t) is the angular velocity of the edges of the objects on the retina, represents motion-dependent excitation, whereas −αθ˙(t − δ) is an inhibitory term proportional to the angular size of the object. Equation 2 describes the approach of a solid square or circle based on the angle that it subtends with respect to the center of the eye and its temporal derivative. Hence, it cannot predict the responses to more complex (e.g., textured) stimuli or to the different approach angles explored in the present study. We therefore generalized Equation 2 to take explicitly into account the two-dimensional projection of the object on the retina. In the following, we will set the delay parameter δ = 0 for simplicity. This does not affect the generality of the arguments presented below because δ only causes a global translation of the firing rate along the time axis and thus a translation along the vertical axis of the linear relationship between the peak firing time relative to collision and l/‖v‖ (see Figs. 3-7). Similarly, the static nonlinearity (g ○ exp)(·) does not influence the prediction of the peak firing rate time by the model because it is a monotonic function. It was therefore not taken further into consideration.
The eye was modeled as a hemisphere. Because the eye radius did not affect the output of the model, it was normalized to R = 1 in the following description (see Appendix for complete equations, including explicitly R). The object was projected onto the eye surface from the center of the hemisphere (Fig.2A, left diagram). The motion-dependent excitatory input was obtained by integrating along each luminosity edge on the hemisphere the logarithm of the optic flow generated by the expanding object. Biophysically plausible mechanisms to compute this quantity have been described in insects (Hildreth and Koch, 1987; Reichardt et al., 1988; Zanker et al., 1999). Let [p0; p1] be a closed interval and γ(t, p), p ∈ [p0; p1] denote a parametrization of the edge. Each value of the parameter p is required to describe the same physical point on the edge, independent of time. Letγ˙t0(p) = ∂/∂t γ(t, p)‖t=t0 denote the instantaneous velocity of expansion of the point γ(t, p) at time t = t0 and letnt0(p) be the outward pointing normal unit vector (Fig. 2A, right diagram). The excitation, Ex (t0, l/‖v‖), caused by the motion of the edge on the retina was obtained by integrating the logarithm of the optic flow weighted by a factor w(γ(t0, p)) along the luminance edge: Equation 3 In this equation dlγ(·) denotes the line element associated with the curve γ(t0, ·) on the sphere (Dubrovin et al., 1991), and 〈 ·; · 〉 is the Euclidian scalar product. In Equation 3, the optic flow term has been renormalized by a constant Ccutoff: Equation 4The constant Ccutoff was chosen so that log+〈γ˙t0(p);nt0(p)〉 was equal to zero when the object subtended a half-angle ½θcutoff = 1° (corresponding to the optical resolution of the eye). During the remainder of the approach trajectory log+〈γ˙t0(p); nt0(p)〉 was positive. The weight factor w(γ(t0, p)) depended on the distance of the local motion stimulus to the central axis of the eye. The weight factor mimicks the experimental decrease in LGMD sensitivity observed as a local motion stimulus is shifted from the eye center toward the periphery. Experimental evidence suggests that this sensitivity decreases faster in the vertical than in the horizontal plane. Iso-contour lines of equal weight were therefore chosen as elongated ellipsoids (Fig. 2B, right diagram). Let (ϑ ϕ) denote the coordinates of γ(t0, p) with respect to the coordinate system illustrated in Figure 2B, (left diagram). Then: Equation 5with a = 2 and ½θcutoff = 1°. The choice a = 2 yields the ellipsoidal contour lines illustrated in Figure 2B, and ½θcutoff = 1° causes the weight factor to saturate beyond the optical resolution of the eye. The dependence of w(ϑ, ϕ) on the angular distance to the eye center in the vertical and horizontal planes is illustrated in Figure 2C.
The size dependent inhibitory term In (t, l/‖v‖) was obtained by converting the surface, S(t), covered by the object on the retina into an equivalent angle, θin(t): Equation 6 Note that in this equation S(t) is normalized by the surface of the eye (=2π = surface of a hemisphere of radius 1). The value of θin as a function of S is plotted in Figure 2D. For comparison with experimental data, the parameter α was set equal to 5, a value that lies in the middle of the typical range observed for a large sample of LGMD neurons (Gabbiani et al., 1999a, their Fig. 6A).
When simulating approaches from angles other than 90°, the optic flow was integrated only across that portion of the object that effectively projected onto the eye. Similarly, only the surface of the object covering the retina was taken into account in Equation 6.
Appendix gives a precise mathematical analysis of this model for an approaching circle. Appendix describes the parametrization used for the edges of approaching squares, as well as the conditions describing the portion of the approaching square projecting onto the retina for approach angles different from 90°.
For the concentric square pattern, the excitatory term was obtained by linear summation of the excitatory terms corresponding to each of the concentric squares composing the pattern. If we denote by Exsquare(t, l/‖v‖) the excitation at time t attributable to a square of half-size l approaching at speed ‖v‖, then: Similarly, the excitation term for the checkerboard pattern, Excheckerboard(t, l/‖v‖), was obtained by taking into account the two additional expanding vertical and horizontal inner edges. If the corresponding term is denoted by Exinner(t, l/3·‖v‖), then: Appendix explains how Exinner(t, l/3·‖v‖) is computed.
The symbolic algebra package Maple (Maple Inc., Waterloo, Ontario, Canada) was used to obtain numerical values for Ex (t, l/‖v‖) and In (t, l/‖v‖) under the tested experimental conditions. Although In (t, l/‖v‖) could be computed in closed form, the indefinite integral corresponding to Equation 3 could not be obtained in closed form for approaching squares. The integrand of Equation 3 was therefore computed symbolically and then integrated numerically using an 8th order Newton-Côtes quadrature method (Forsythe et al., 1977). Numerical calculation of Ex (t, l/‖v‖) and In (t, l/‖v‖) for all condition tested (l/‖v‖ = 10, 20, 30, 40, 50 msec for solid discs, squares, concentric square patterns, and checkerboard patterns approaching at 90°; same values of l/‖v‖ for solid squares approaching at angles from 90° down to −30° in 15° steps), and all combinations of model parameters tested (a = 1, 2, and ½θcutoff = 0°, 1° in Eq. 5) took approximately 1 week on a computer equipped with a 450 MHz Pentium III processor (Intel, Santa Clara, CA). None of the results depended qualitatively on the choice of a or ½θcutoff in Equation 5. Maple simulation files and Matlab functions reproducing Figures 8 and 9 may be obtained on the World Wide Web at http://glab.bcm.tmc.edu by following the “Publications” link and selecting the page corresponding to the title of this article.
RESULTS
The results described below are based on recordings and complete data analysis from 58 different animals. Part of the data used in Fig.3D was collected in 35 additional animals (S. americana) for different purposes.
Is the angular threshold computation species dependent?
The response of LGMD/DCMD in L. migratoria to the simulated approach of a solid square on a collision course with the animal is illustrated in Fig. 1A. The time course of visual stimulation on the retina is characterized by the ratio of the half-size of the square, l, and its approach velocity, ‖v‖, (Fig. 1A, inset) [see also Materials and Methods and Gabbiani et al. (1999a)]. The response starts early during the approach, when the object reaches approximately 10° in visual angle. It gradually increases and reaches its peak value before collision for large values of l/‖v‖ (≥10 msec) (Fig. 1A–C). Each stimulation was repeated 8–10 times to obtain estimates of the mean time of peak firing rate and its SD (Fig. 1B). As may be seen in Figure 1C, the time of peak firing rate relative to collision (★) decreases systematically as the stimulation parameter l/‖v‖ decreases (i.e., as the approach velocity, ‖v‖, of the moving object increases for a fixed object half-size, l). Figure 3A replots the time of peak firing rate, tpeak, obtained from Figure1C as a function of the kinematic parameter l/‖v‖ characterizing the approach [see alsoGabbiani et al. (1999a), their Figs. 3, 4]. To a very good approximation, tpeak follows a linear relationship with l/‖v‖. Similar results have been reported for S. americana (Gabbiani et al., 1999a) and are illustrated in Figure 3B for one specimen of this species. Furthermore, as reported in Gabbiani et al. (1999a), this linear relationship between tpeak and l/‖v‖ implies that the peak in firing rate occurs a fixed delay (5–40 msec on average in the present experiments) after the object has reached a fixed angular threshold size (between 15° and 40° on average) on the locust's retina, independent of the value of l/‖v‖ (Gabbiani et al., 1999a, their Fig. 5). To test further the linear dependence of peak firing time on l/‖v‖, an additional eight animals were presented looming stimuli at 10 values of l/‖v‖ with a protocol identical to the one used inGabbiani et al. (1999a, their Fig. 4). The results of these experiments were also similar to those reported in S. americana. Thus, in L. migratoria as in S. americana, the peak firing rate can be seen as the output of an angular threshold detector and may contribute to trigger visually evoked escape behaviors [Hatsopoulos et al. (1995); see Gabbiani et al. (1999a) for a detailed discussion]. To ascertain whether this is also the case in a second African locust species, we performed similar experiments onS. gregaria. Figure 3C illustrate the results of experiments performed on one animal of this species. In all experiments analyzed, the same linear relationship was found in both the American and African species. Figure 3D plots the mean angular threshold sizes and delays observed for S. gregaria (n = 5, ★) and L. migratoria (n = 11, ●). These values are well into the range observed in a large sample of animals from the American locust S. americana (n = 65, ○).
Angular threshold computation and target shape
Shape differences between solid objects looming toward the locust's eye translate into temporal and spatial differences in the stimulation of ommatidia on the retina. In the case of a disc and square having the same radius (or half-size), l, the temporal sequence of light intensities during approach is the same only along the two horizontal and vertical symmetry axes passing through the center of the objects. For a disc, the angular velocity of expansion on the retina at any given time points out from the center with a magnitude that is constant along the whole disc boundary. Because this vector is always perpendicular to the disc boundary, it coincides with the optic flow vector (i.e., the projection of the expansion velocity vector onto the unit normal to the boundary) that can be measured locally at the retina (Appendix ). By contrast, the angular velocity of expansion of a square is not always perpendicular to the square edges and its magnitude is greatest at its four corners (Fig.4A, inset). After projection on the unit normal to the boundary, the optic flow vector magnitude is—on the contrary—greatest at the middle of an edge rather than at its corners. Both quantities may be computed from the equations given in Appendix . The ratio of the optic flow magnitude at the corner of a square to that of a disc decreases monotonically during approach from a peak value of 1 when both objects are far away to a minimum value of 0.5 at collision time. Thus, the motion stimulus generated by an expanding square is always less strong than the one generated by a disc, whereas the surface that a square of half-size l subtends at the retina is always larger than the one subtended by a disc of radius l.
We investigated whether such shape differences affected the timing of peak firing by presenting to the same animal square (■) and circular (●) targets interleaved pseudo-randomly over the course of an experiment. Figure 4A illustrates the outcome in one preparation: the differences in mean times of peak firing rate measured with both objects could not be distinguished from the intertrial variability (SD) (Fig. 4A, solid bars). Figure4B summarizes the outcome of five such experiments in different preparations by plotting mean angular thresholds versus mean delays obtained for square (■) and circular (●) targets and their SDs (solid lines parallel to the axes). For each mean angle and delay pair, the SDs outline the boundary of the two-dimensional 68.3% confidence region on their value (Gabbiani et al., 1999a, their Fig. 4B; Press et al., 1992, sect. 15.6). As can be seen from the plot, the two-dimensional confidence regions strongly overlap, and therefore mean angular thresholds and delays were not statistically different between the responses to these two stimulus conditions. Figure 4C plots the mean delays for square versus circular targets and their SDs along the delay axis. No significant trend could be observed, and the data always lay within 1 SD of the diagonal characterizing equal delays under both conditions. By contrast, the mean angular thresholds were consistently above the diagonal (Fig. 4D), corresponding to smaller average threshold angles, or equivalently earlier peak times for square targets. Differences in mean threshold angles were not significant in single preparations. In two cases, the SDs of the threshold angles did not intercept the diagonal, but these two points were located within two SDs of the diagonal and were thus not significantly different at the 95.4% level. The difference in mean threshold angles averaged over all preparations was small: 3.2 ± 2.1° (mean ± SD), but statistically significant (paired t test, p < 0.005). Thus, these results point to a small (i.e., not statistically significant in single preparations) but consistent trend across preparations of higher threshold angles and later peak firing times for disc targets.
Angular threshold computation and target texture
Variations in target texture will also result in different temporal and spatial activation of ommatidia on the locust's retina. We investigated whether such changes affected the time of peak firing rate during approach by comparing the responses obtained for solid square targets to two different textured targets: a checkerboard pattern (CBP) (Fig. 5A,inset) and a concentric squares pattern (CSP) (Fig.6A, inset). For the CBP, the inner edges move in opposite directions as the object grows on the retina, whereas for the CSP, inner edges consistently follow the leading edges of the square although at a different instantaneous angular velocity. In both cases, the luminance of the pattern object was set at 30% (dark gray) and 0% (black) of the background luminance (95 cd/m2); similar results were obtained in two animals with CBP relative luminances of 70% (light gray) and 0%. Figures 5A and 6A illustrate the outcome of one experiment for the CBP (Fig.5A, ▵) and the CSP (Fig. 6A, ▵), respectively. In both cases, the time of peak firing rate was statistically indistinguishable from the controls obtained with black square targets (Figs. 5A, 6A, ■). Figures 5B and 6B plot mean angular threshold versus delay obtained with the CBP (▵, n = 10 animals) and the CSP (▵, n = 12 animals), respectively, and compares them with those obtained in the same preparations with black square targets (■). Experiments performed on the same preparation are connected with a dashed line. In both cases, no consistent changes were observed with respect to the control condition (■). This was confirmed by first plotting the mean delays and SDs obtained in the control condition versus those of the two test conditions (CBP and CSP). Figures 5C and 6C illustrate the results for the CBP and CSP, respectively: the mean delays in the test versus control conditions lie within 1 SD of the diagonal and therefore are not significantly different from one another. Figures 5D and6D illustrate that the mean threshold angles were also not statistically different from those obtained in control experiments. Two mean threshold angles of 12 preparations lay >1 SD away from the diagonal, but <2 (Fig. 6D). Therefore the difference was not statistically significant at the 95.4% confidence level.
Angular threshold computation and target approaching angle
We investigated the effect of the angle of approach on LGMD/DCMD firing by rotating the animal around a virtual vertical axis passing through the center of the stimulated eye. This allowed us to present looming targets from a wide range of approaching angles. Let 0° be the position at which the animal's front faces the stimulation screen and positive angles denote counterclockwise rotation. If we consider the right eye, the angular range covered target approach angles from the back (135°) on the same side as the stimulated eye to −45° on the contralateral side (Fig.7A, inset).
Figure 7A illustrates an experiment in which the animal was stimulated at 0° (●, facing the screen) and at 90° (▵, perpendicular to the screen). No statistical differences were observed in the timing of the peak firing rate under these two conditions. The variability in the time of peak firing rate, however, usually increased from 90° to 0°. In Figure 7A for example, the largest SDs in the peak firing time at each value of l/‖v‖ were in all but one case (l/‖v‖ = 20 msec) attributable to stimulation at 0°. Figure 7B summarizes the results of 5 similar experiments of 13 that had the most reliable responses at 0°. The mean angular thresholds and delays under these two experimental conditions are plotted together with their SDs, and experiments performed on the same animal are connected by a dotted line. There were no statistical differences between the mean delay values in all cases. In one experiment (see ▵, top left), the mean threshold angles were different at the 95.4% level, whereas the remaining threshold angles showed no statistically significant differences (one mean threshold angle, indicated by anarrow, lay within 2 SDs of its control experiment and was therefore not statistically different at the 95.4% confidence level). Here also, in all but one case, variability increased from 90° to 0°. This increase in variability was correlated with the number of action potentials elicited per stimulus presentation at various approach angles. Figure 7C plots for one experiment the mean number of action potentials (±SD) elicited at three values of l/‖v‖ for looming black squares with approach angles 135°, 90°, 45°, and 0°. As in Figure 7, A and B, no statistically significant differences were observed in the timing of the peak firing rate in all cases analyzed (n = 5). Furthermore, the number of action potentials elicited per trial was remarkably constant over a large portion of the visual field ranging from 135° to 45°. Between 45° and 0°, however, the response declined abruptly, by >50% in this experiment; similar results were obtained in the remaining four experiments. In the experiment illustrated in Figure 7C, there was a clear decrease in the number of action potentials elicited per trial as l/‖v‖ decreased from 50 to 10 msec. However, this observation could not always be made so clearly on other preparations (Gabbiani et al., 1999a).
This decrease in the number of spikes per trial produced at 0° was compensated for by the contralateral DCMD. Figure 7D illustrates the number of action potentials (±SD) recorded in each of the two DCMDs during an experiment in which the target approaching angle was varied between 5 values in the frontal visual field: 45°, 22.5°, 0°, −22.5°, and −45°. The response recorded in the left connective (★) decreased from 45° to 0°, as in Figure7C, and essentially vanished at −45°. Concomitantly, the number of action potentials recorded from the right connective (●) increased over the same range, resulting in the total number of action potentials elicited from both DCMDs per trial (▵) being constant over the entire range (−45°–45°). Similar results were obtained in n = 11 additional animals.
Two-dimensional linear summation model for motion-dependent excitation and size-dependent inhibition
The experimental results described above were compared with those obtained in a two-dimensional model of the responses of the LGMD to looming stimuli. The motion-dependent excitatory input of the model was obtained by integrating the logarithm of the optic flow along the luminance edges of the expanding object weighted by a factor depending on the distance of the local motion stimulus to the central axis of the eye (Figs. 2A–C and Eqs. 4 and 5). When the approaching object contained additional luminance edges, such as the CBP and the CSP described above, excitatory inputs were summed linearly for all luminance edges. The size-dependent inhibitory input was obtained at each time point by computing the surface covered by the object on the retina and transforming it by using the nonlinear function plotted in Figure 2D. Figure 8, A and B, illustrates the time course of excitation and inhibition for a disc approaching perpendicular to the body axis toward the center of the eye. The excitatory and inhibitory terms were subtracted from each other (Fig. 8C), and the times of peak firing rate were computed and plotted as a function l/‖v‖ (Fig. 8D, ●). As may be seen from this Figure, the prediction of the model is a linear function of l/‖v‖, although the slope of the linear relation is slightly smaller (4.69) than the prediction of the one-dimensional model [dashed line; slope identical to α = 5) ofGabbiani et al. (1999a)]. This, in turn, corresponds to a slightly larger threshold angle (24° vs 22°). The difference, however, is well within the uncertainty of the experimental estimates for these parameters (Fig. 4B, D). As proven in Appendix , the two-dimensional model reduces exactly to the one-dimensional model for a looming disc if a less realistic choice of parameters is made, corresponding to circular iso-contour lines for the weight factor (a = 1 in Eq. 5) and no weight saturation beyond the optical resolution of the eye (½θcutoff = 0° in Eq. 5).
Figure 9 summarizes the results for the remaining stimuli used in our experiments. As may be seen by comparing Figure 9A with Figure 8D, the time of peak firing rate for a square is closer to the prediction of the one-dimensional model, and accordingly the slope is higher (4.86) corresponding to a smaller threshold angle (23°). These results are consistent with those observed experimentally, although the difference in threshold angles in the model (∼1°) is smaller than the mean value derived from Figure 4D (∼3°). In the model, the earlier peak firing time for a square stimulus versus a disc is not surprising: as explained in Results, the optic flow associated with a square is smaller than that of the corresponding disc and its surface is larger. Consequently, inhibition will start to dominate the excitatory term earlier, leading to earlier peak times.
Figure 9, B and C, illustrates the results for the CBP and the CSP patterns: the time of peak firing still varies linearly with l/‖v‖, but as the strength of the excitatory term is increased by the additional moving edges in the stimulus, the timing of the peak moves closer to collision time. The values of the slopes in these two cases (1.56 and 0.92) correspond to predicted threshold angles of 65° and 95°, respectively, that are incompatible with the experimental observations of Figures 5 and 6. Similarly, presentation of square targets approaching at different angles disrupted significantly the predictions of the model for deviations as small as 15° with respect to the 90° approach direction (Fig. 9D) (the slope equals 3.73 and the value of the threshold angle is 30°). Approach angles of 45° and 0°, such as those used experimentally, could not be explained by the linear model (slopes of 1 and 0 corresponding to threshold angles of 90° and 180°, respectively).
DISCUSSION
Earlier results have shown that the peak firing rate of LGMD/DCMD always occurs at a fixed delay after a solid black square looming toward the animal (perpendicular to the body axis) has reached a threshold angle (or size) on the locust's retina, regardless of the actual size of the square or speed of approach (Gabbiani et al., 1999a). Building on this work, we investigated in three locust species whether this angular threshold computation depends on stimulus changes expected to alter the spatial and temporal activation of inputs onto the dendritic tree of the LGMD. The results indicate that this angular threshold computation remains invariant under a wide range of experimental conditions. Furthermore, this invariance is not compatible with linear summation of excitatory and inhibitory inputs within the dendritic tree of LGMD.
Identical computation implemented in three locust species
Our initial observations were based on experiments with S. americana. Here, they were generalized to two African locust species, S. gregaria and L. migratoria. S. gregaria resembles very closely S. americana, but hybridization studies have shown it to be a distinct species (Harvey, 1981; Jago et al., 1979). It is the only Old World species of the Schistocerca genus that has richly radiated to comprise at least 42 species in the New World (Otte and Nasrecki, 1997). The 42 New World species are thought to have evolved from S. gregaria after one or more crossings of the Atlantic Ocean (Ritchie and Pedgley, 1989; Amedegnato, 1993). In contrast, L. migratoria belongs to a different subfamily of Acridids (Otte and Nasrecki, 1997) that evolved independently for >100 million years (Vickery, 1987; Whitington and Bacon, 1999). Therefore, angular threshold detection might be an escape mechanism common to many other orthopteran species (Rowell, 1971). Sun and Frost (1998)recast this hypothesis in a broader context by showing that a class of neurons in the nucleus rotundus of the pigeon have response properties identical to those of LGMD/DCMD neurons. Whether similar responses also arise in mammalian neurons remains to be investigated (Luksch et al., 1998). Finally, our results rule out interspecies differences as an explanation for the results reported by Rind and Simmons (1997) [see Gabbiani et al. (1999a,b) for a complete discussion].
Invariance to target shape
We investigated the effect of target shape on LGMD/DCMD responses by comparing the timing of the peak firing rate during the approach of circular and square targets. The difference between these two stimuli is detectable in principle when the objects cover nine ommatidia or 4° in visual angle on the retina, that is, before LGMD and DCMD begin to respond during approach. However, simulation results using the two-dimensional model of the responses of the LGMD suggested that only a small decrease (∼1°) in the angular threshold size should be observed for a square target with respect to a circular target, because of a decrease of optic flow stimulation and an increase of size-dependent inhibition. This is in agreement with experimental observations that point to a ∼3° decrease in average angular threshold size for square targets. This change lies within the margin of accuracy of our measurements.
Invariance to target texture
In contrast, the two textured targets used in the present experiments contain inner edges that may be expected to result in a substantial additional spread of excitation over the time course of approach. The two-dimensional model of the responses of the LGMD demonstrates that linear summation of the excitation attributable to inner luminance edges in the CBP and CSP significantly alters the peak firing rate time for both targets. Such changes, however, were not observed experimentally. A possibility is that these additional excitatory inputs are filtered out presynaptically to LGMD. This appears unlikely: the contrast of the CBP and CSP were in the range for which LGMD/DCMDs are known to respond vigorously to expanding objects, thus ruling out early gain control mechanisms. In the case of the CSP (and to a lesser extent of the CBP), one might argue that excitatory inputs triggered by the outer object contour/edge might have caused a strong decrement in the response to inner edges, because the same ommatidia were stimulated successively within a brief interval. However, repeated local stimulation with moving targets showed that habituation to moving stimuli occurs much more slowly than could occur with the present stimuli (Krapp and Gabbiani, 2000). Finally, lateral inhibition is unlikely to be effective in the case of the CSP (and presumably played no role for the CBP) because successive edges were separated by more than one to two interommatidial angles over most of the approach trajectory (as soon as the object subtended 16° in visual angle, assuming an interommatidial angle of 2°) (Horridge, 1978).
Invariance to target approach angle
When the angle of target approach was varied, the timing of the peak LGMD/DCMD firing rate remained invariant over 135° on one side of the animal's longitudinal body axis. Both pairs of LGMD/DCMDs thus provide an invariant warning signal to potentially threatening objects looming toward the animal over at least 270° in the horizontal plane. This is remarkable because the temporal and spatial activation sequence of inputs onto the dendritic tree of the LGMD may be expected to differ completely across these conditions. For an object approaching at 90°, the expansion will be symmetric around the center of the eye, whereas at 135° the expansion is maximally asymmetric, with one vertical edge sweeping over most of the receptor array. In accordance with this, the two-dimensional model of LGMD predicted non-invariant responses over the range of approach angles tested experimentally.
In the frontal visual field, responses decreased rapidly as the approach angle shifted toward the contralateral side. This was presumably caused by two factors: (1) the ommatidia sample only a limited portion of the contralateral visual field (∼10°), implying that the vertical edge of a looming square opposite to the eye is positioned outside of the receptive field of the LGMD over most of a frontal approach and cannot excite the neuron; (2) the sensitivity to moving targets is weaker in the frontal visual field than at the center of the eye (Rowell, 1971; Krapp and Gabbiani, 2000). This decrease in DCMD response leads in turn to an increase in peak firing time variability that could conceivably be reduced by averaging over both responses of the DCMDs because they are independent. Such a procedure is plausible, given the relatively constant number of action potentials produced by both DCMDs over the frontal field of view. Alternatively, preparation for escape in natural conditions might involve repositioning of the animal with respect to the target (Hassenstein and Hustert, 1999).
Robertson and Gray (1997) first recorded the responses of both DCMDs to frontal targets in tethered flying locusts. On the basis of behavioral studies (Robertson and Johnson, 1993), they suggested that the relative number of action potentials in both DCMDs might trigger directional steering maneuvers during flight to avoid frontal obstacles. Our results confirm that both DCMDs can encode target approach direction (left vs right; see the SD in Fig. 7D at ±22.5°), although the increased variability of responses to frontal targets suggests that directional steering behaviors would be triggered outside of their optimal operating range. Accurate directional estimates may require additional signals, probably available from other lobula movement-sensitive neurons (Gewecke and Hou, 1993).
Biophysical mechanisms responsible for invariance to target texture and approach angle
One explanation for the observed texture response invariance could be nonlinear saturation of postsynaptic excitatory inputs. In addition, active membrane conductances might contribute to suppress the excitation expected from inner luminance edges. Invariance to approach angles could result from the geometric arrangement of excitatory inputs onto the dendrites of LGMD as well as passive attenuation and active amplification of contributions associated with specific positions in the visual field. To be tested rigorously, these hypotheses will require a precise mapping of the excitatory retinotopic inputs onto the dendritic tree of the LGMD, as done for wide-field tangential neurons in flies (Egelhaaf and Borst, 1995) and an assessment of the passive and active properties of the neuron.
LGMD as a model for invariant feature extraction
A major threat to survival for locusts and other insects in the wild is embodied by predatory birds. Species foraging on insects vary in shape, body texture, and predatory tactics. They will typically engage in chases from various positions with respect to their prey (Jablonski, 1999). It is therefore not surprising that locusts and other insects have evolved escape behaviors to confront these dangers. In the locust species investigated here, one can extract from the LGMD/DCMD peak firing time the angular threshold of an approaching object regardless of its shape, texture, or direction of approach. Because visual stimuli reaching each ommatidia on the locust's retina are thought to be first integrated temporally and spatially in the dendritic tree of the LGMD, our results suggest that this invariance is implemented within the LGMD dendrites. LGMD is uniquely identifiable from animal to animal and accessible to intracellular recordings; it therefore offers a model for the study of angular threshold computation and its biophysical implementation (Gabbiani et al., 1999a). The present study suggests that LGMD is also an ideal model for studying the emergence of invariant receptive field properties at the dendritic and biophysical levels.
Appendix
In this appendix we show that if a = 1 and ½θcutoff = 0° in Equation 5, the two-dimensional model of Equations 3, 5, and 6 reduces to the one-dimensional model of Equation 2 for a circle approaching at an angle of 90°.
The spherical coordinate system of Figure 2B corresponds to where R is the radius of the eye. Let us denote by Sx>02 the hemisphere used to model the eye. The scalar product induced on Sx>02 by the three-dimensional Euclidean scalar product is described in spherical coordinates by the matrix: Equation A1(Dubrovin et al., 1991). The boundary of the circle on the retina is described by the equation: with ϑ(t) = tan−1l/vt. The subscript (·)T denotes matrix transposition. We compute first the excitatory term: Equation A2Because ½θcutoff = 0°, the second term in Equation 5 vanishes and the first term reduces to Equation A3where we have used a = 1 and the trigonometric identity sin2 ϕ + cos2ϕ = 1. The tangent vector to the curve γ(t0, ·) is given by: and its norm by: The line element in Equation EA2 is given by: Equation A4The expansion vector is given by: Equation A5and the outward pointing unit normal vector is easily seen to be equal to: Equation A6because: and We compute only the first term of Equation 4 (the second term is an easily calculated constant). Using Equations EA5 and EA6: Equation A7 Plugging Equations EA3, EA4, and EA7 in A2 we deduce that: which matches the excitatory term of Equation 2, because ½θ˙(t) = ϑ˙(t) for a circle approaching at 90°.
The surface covered by the disc on the eye at time t0 is given by: (Dubrovin et al., 1991). Inverting this latter result yields: and therefore 2ϑ(t) = θin(t) for a sphere of radius 1.
Appendix
In this appendix we give the formulas used to compute Ex (t, l/‖v‖) and In (t, l/‖v‖) for squares approaching from various angles with respect to the eye's central axis. These formulas were used to write the Maple scripts available athttp://glab.bcm.tmc.edu (under the title of this article after following the Publication link).
Top and bottom boundaries
The three-dimensional approach trajectory of the top and bottom boundaries (Fig. 2B) of a square is described by: where each value of b ∈ [−1; 1] corresponds to a point on the square edge. The projection of the edge on the eye is given by: with and where we have set β = l/v.
Left and right boundaries
With a notation similar to the one used above, the right and left square edges are described by (Fig. 2B). Their projection onto the eye surface is given by: For the right boundary, and For the left boundary, ϕ(b) is given by with ϑ(t, b) unchanged.
Approach angles different from 90°
For approach angles different from 90° (not perpendicular to the body axis), the optic flow was computed in a coordinate system adapted to the approach direction (i.e., with its x -axis passing through the center of the square). The parametrization of the square edges given above therefore remained unchanged.
Weight factor for arbitrary approach angles
To compute the weight factor w(ϑ, ϕ), the coordinates of a point on the square edge with respect to the coordinate system whose x -axis passes through the eye center are needed (Fig. 2). Let (x y z)T be the coordinates of the point in the coordinate system adapted to the approach direction (used to compute the optic flow component of excitation; see last paragraph) and let ψ denote the angle of the approach direction with respect to the eye center (i.e., ψ = 0 for an approach angle of 90°). Let (xψyψ zψ)T be the corresponding coordinates in the coordinate system whose x -axis passes through the eye center. The values of (xψ yψzψ)T are obtained by rotating (x y z)T by an angle ψ about the z -axis, We next project the point onto the surface of the eye: with D = = and compute the weighted distance to the eye center, For the bottom and top boundaries, we have at time t for the point on the edge corresponding to the parameter b, and so that A similar calculation for the right and left boundaries yields, respectively, For approach trajectories with angles 0° < ψ < 90° the time tthres at which the left edge of the object reaches the eye boundary is given by tan ψ = vtthres/ l, or, equivalently, tthres = β tan ψ. The portion of the top and bottom boundaries remaining on the eye after tthres is given by b < bthres(t), bthres(t) = t/β tan ψ > 0. For the approach trajectory ψ = 90°, tthres = −∞ and bthres = 0. Finally, for ψ > π/2 rad (90°), tthres = β/tan(ψ − (π/2)), and
Surface for arbitrary approach angles
To compute the surface of the square projecting on the retina for arbitrary approach angles, we use the following spherical coordinate system: with the x -axis passing through the center of the approaching square. The surface element is given by = R2‖ cos ϑ‖ (Dubrovin et al., 1991), and the surface itself is given by: We distinguish the following cases: (1) for ψ ≤ 90° and t ≤ tthres: ϕmax/min = ±tan−1 β/t and ϑbound(t, ϕ) = tan−1 (β/t cos ϕ); (2) for ψ ≤ 90° and t ≥ tthres: ϕmax = π/2 − ψ and ϕmin, ϑbound are given by the above formulas; (3) for ψ ≥ 90° and t ≤ tthres: S(t) = 0; (4) for ψ ≥ 90° and t ≥ tthres: ϕmax = −(ψ − π/2), ϕmin = −tan−1 β/t and ϑbound is as above.
Checkerboard pattern
It is easy to see that the projection of the horizontal and vertical inner edges of the checkerboard pattern on the eye may be parametrized with similar formulas as the outer edges. The top and bottom edges are given by: with A similar formula holds for the left and right vertical edges. The corresponding excitatory term Exinner(t, l/3·‖v‖) is obtained as usual by computing symbolically the integrand of Equation 3 and subsequent numerical integration.
Footnotes
This work was supported by the German American Academic Council, by a National Institute on Deafness and Other Communication Disorders (National Institutes of Health) grant to G.L., by a National Science Foundation grant to C.K., and by the Sloan Center for Theoretical Neuroscience at Caltech. Experiments and data analysis were carried out by F.G. and C.M.; F.G. and G.L. co-wrote the paper. We thank Dr. A. Hassanali (ICIPE) and Dr. H. Krapp (University of Bielefeld) for kindly furnishing L. migratoria and S. gregariaspecimens, as well as Dr. D. Sheinberg for help in setting up the data acquisition system. F.G. and C.M. thank Dr. C. Koch for many discussions and his encouragement over the course of this project.
Correspondence should be addressed to Dr. F. Gabbiani, Division of Neuroscience, One Baylor Plaza, Houston, TX 77030 (E-mail:gabbiani{at}bcm.tmc.edu) or to Dr. G. Laurent, Division of Biology, Caltech, Pasadena, CA 91125 (E-mail:laurentg{at}caltech.edu).