When complex neuronal structures may not matter
Abstract
Much work has explored animal-to-animal variability and compensation in ion channel expression. Yet, little is known regarding the physiological consequences of morphological variability. We quantify animal-to-animal variability in cable lengths (CV = 0.4) and branching patterns in the Gastric Mill (GM) neuron, an identified neuron type with highly-conserved physiological properties in the crustacean stomatogastric ganglion (STG) of Cancer borealis. We examined passive GM electrotonic structure by measuring the amplitudes and apparent reversal potentials (Erevs) of inhibitory responses evoked with focal glutamate photo-uncaging in the presence of TTX. Apparent Erevs were relatively invariant across sites (mean CV ± SD = 0.04 ± 0.01; 7–20 sites in each of 10 neurons), which ranged between 100–800 µm from the somatic recording site. Thus, GM neurons are remarkably electrotonically compact (estimated λ > 1.5 mm). Electrotonically compact structures, in consort with graded transmission, provide an elegant solution to observed morphological variability in the STG.
https://doi.org/10.7554/eLife.23508.001Introduction
Neuronal circuits can generate stable output despite variable underlying parameters across animals (Marder and Goaillard, 2006; Marder, 2011). Work in invertebrate central pattern-generating circuits shows that circuit function can be maintained across animals despite variable synaptic, intrinsic, and modulator-induced currents in identified constituent neurons (Prinz et al., 2004; Schulz et al., 2006; Marder and Goaillard, 2006; Goaillard et al., 2009; Norris et al., 2011; Marder, 2011; Roffman et al., 2012; Williams et al., 2013; Gutierrez et al., 2013; Rodriguez et al., 2013). This previous work has revealed the principle that neurons compensate for variable ionic conductances at the circuit (Grashow et al., 2010) and single-neuron levels (Tobin et al., 2009; Ball et al., 2010; O'Leary et al., 2013, 2014).
The unique physiology of any given neuron is a consequence of its palette of ionic conductances, and its morphology (Mainen and Sejnowski, 1996; Vetter et al., 2001). A neuron’s distributed, geometric cable properties: the length, diameter, taper, and branching of its neurites, shape passive current flow and voltage propagation (Rall, 1959, 1977; Goldstein and Rall, 1974). This resulting electrotonic structure plays a central role in determining whether voltage signals arising at disparate sites across the dendritic tree are integrated or segregated (Rall, 1959, 1967, 1969a, 1969b; Rall and Rinzel, 1973, 1974; Goldstein and Rall, 1974; Agmon-Snir and Segev, 1993; Vetter et al., 2001). In this way, ion channel expression, when superimposed on morphology, gives rise to neuronal physiology, input-output computations, and circuit-level function (for reviews, see Koch and Segev, 2000; London and Häusser, 2005).
In the present study, we investigated the physiological consequences of animal-to-animal variability in neuronal morphology in the crab stomatogastric ganglion (STG), a central pattern-generating circuit composed of 26–27 neurons. The identified neuron types of the STG exhibit highly conserved physiological waveforms and circuit-level functions (Harris-Warrick et al., 1992), despite their complex and variable morphologies across animals (Wilensky et al., 2003; Baldwin and Graubard, 1995; Bucher et al., 2007; Goeritz et al., 2013; Otopalik et al., 2017). STG neurons display large somata (50–100 µm in diameter) and primary neurites that ramify throughout the STG neuropil (Selverston et al., 1976; King, 1976a, 1976b; Baldwin and Graubard, 1995; Kilman and Marder, 1996).
Synaptic transmission between neurons is predominantly graded, inhibitory cholinergic and glutamatergic transmission (Eisen and Marder, 1982; Marder and Eisen, 1984; Maynard and Walton, 1975; Graubard et al., 1980; Manor et al., 1997, 1999). Synaptic sites are sparsely distributed throughout finer processes in the neuropil region (King, 1976a, 1976b). There are no synapses on somata in these neurons (King, 1976a, 1976b) and spike initiation zones are located in the periphery where axons exit the ganglion (Raper, 1979; Miller, 1980). Thus, synaptic integration and release occurs predominantly in the neuropil. In individual STG neurons, pre- and post-synaptic sites are often tightly apposed on the same neurites (King, 1976a, 1976b). This juxtaposition of synaptic input and output suggests that current will flow in all directions across the neurite tree, centripetally and centrifugally, in the intact circuit, and allow for integration of voltage signals arising from disparate loci on the neurite tree, should the neuron be sufficiently electrotonically compact.
Previous investigation of electrotonic structure in STG neurons employed simultaneous electrophysiological recordings at the soma and primary neurite, wherein electrodes were separated by several hundred microns (Miller, 1980; Golowasch and Marder, 1992). These recordings showed low-pass filtering of high-frequency voltage events across the primary neurite, as is typically imposed by the membrane capacitance (Rall, 1977). In contrast, slow voltage oscillations were subject to less electrotonic decrement. This frequency-dependent decay of electrical signals was consistent with the findings of contemporary theoretical and experimental works (Rall, 1977; Johnston and Brown, 1983; Spruston et al., 1993, 1994; Jaffe and Carnevale, 1999). This early characterization of electrotonic structure in the STG may have been a satisfactory description at the time, given that STG circuit function is mediated predominantly by graded transmission (Eisen and Marder, 1982; Marder and Eisen, 1984; Maynard and Walton, 1975; Graubard et al., 1980; Manor et al., 1997, 1999) and slow oscillations (Graubard and Ross, 1985; Ross and Graubard, 1989), and continues to oscillate in the absence of spikes (Graubard, 1978; Raper, 1979; Graubard et al., 1983;Anderson and Barker, 1981). However, these early experiments left the electrotonic properties of more distal sites and higher-order branches to the imagination.
Relevant voltage events must arise at more distal, finer processes, where pre- and post-synaptic connections are located (King, 1976a, 1976b; Kilman and Marder, 1996). Thus, a full experimental characterization of the electrotonic structure of STG neurons requires electrophysiologically sampling numerous sites across the dendritic tree. Of course, recording at many sites, let alone on the tiniest of neurite processes, presents technical challenges that few have overcome in situ. In one STG neuron type, the gastric mill (GM) neuron, we have utilized focal photo-uncaging techniques in tandem with electrophysiology to examine propagation of voltage events evoked at processes that vary in size and distance from the somatic recording site. We present a surprising case wherein geometrical complexity and variability appear not to constrain passive physiology.
Results
The stomatogastric ganglion (STG) of the crab Cancer borealis is composed of 26–27 neurons situated around a dense neuropil region, wherein each of the neurons branches extensively. This central-pattern-generating circuit mediates the coordinated rhythmic contractions of the animal’s foregut. The gastric mill (GM) neuron is one of fourteen identified neuron types in the STG. There are typically four GM neurons in each animal. Figure 1A shows a schematic of one GM neuron and its axonal projections in the intact stomatogastric nervous system, as dissected for in vitro experiments. Filling the neuron with fluorescent dye reveals the complex morphology of the GM neuron in situ (Figure 1B). The GM neuron ramifies throughout the neuropil region and can be distinguished by multiple axons, projecting to and innervating extrinsic gastric muscles 1a, 1b, 2, 3a, and 3b via the anterior lateral (aln) and dorsal ventral (dvn) nerves (Maynard and Dando, 1974; Selverston and Mulloney, 1974; Weimann et al., 1991). GM neurons participate in the episodic gastric mill rhythm (Hartline and Maynard, 1975). When bursting (Figure 1C), GM activity evokes rhythmic contractions of its target muscles, resulting in grinding movements of the gastric mill ossicles and attached teeth for internal chewing of food (Russell, 1985; Heinzel, 1988). GM neurons are unambiguously identified by matching their spiking activity with spiking units on extracellular nerves known to contain GM axonal projections (Figure 1D; Maynard and Dando, 1974).
Variability in GM neuron morphology
We generated three-dimensional confocal stacks of Lucifer Yellow dye-fills of 14 GM neurons situated in the STGs of 14 different animals. Figure 2A shows six examples of maximum z-projection images generated from GM neuronal dye-fills. GM neurons exhibit large somata (with a mean diameter ± SD of 74 ± 13 µm) and expansive neurite branches that span the neuropil region. Non-axonal neurite branches spanned an average ellipsoid volume of 4.98 ± 1.43 × 106 µm3 (mean ± SD).
Using KNOSSOs software (freely available online at knossostool.org, see Materials and methods), we manually traced and generated 3-dimensional skeletal reconstructions of these neuronal dye-fills (skeletal reconstructions of the six neuronal dye-fills in Figure 2A are shown in Figure 2B; Supplementary file 1-Neuronal Structures Hoc files). From these skeletal reconstructions, we were able to measure total cable lengths, number of branch points, and soma-to-tip path lengths and tortuosities (excluding axons) using a suite of custom morphology analytical tools (see Materials and methods). These neurons present expansive structures, with a mean total cable length of 8840 ± 3678 µm (Figure 2C) and numerous branch points (mean number of branch points ± SD was 155 ± 126; Figure 2D). Notably, across animals, total cable lengths and branch point numbers varied by 40% of the mean (CV = 0.4 for both metrics). Note that these numbers are underestimates and smaller than reported in Otopalik et al. (2017). This is due to the lower resolution used for reconstruction of this large set of neurons, resulting in loss of some of the smallest profiles. Otopalik et al. (2017), reported similar animal-to-animal variability in these metrics across a smaller sample size.
We also measured all soma-to-tip paths within each neuron (Figure 2E). These path length distributions varied within and across neurons, with distances ranging between 200–1000 µm and a mean coefficient of variation of 0.81. The mean soma-to-tip path length was 450 ± 80 µm. Thus, if synaptic voltage events arising at the neurite tip are measured electrophysiologically at the soma, such signals would travel this distance, on average.
In Figure 2A and B, it is visually evident that GM neurons display tortuous three-dimensional structures. This is quantified in Figure 2F, which shows the tortuosities of soma-to-tip paths (path length/Euclidean distance from soma to tip location). For a given soma-to-tip path, a tortuosity of one suggests a minimal, Euclidean path from soma to tip, whereas tortuosities greater than one suggest winding paths that deviate from the minimal Euclidian distance (Figure 2F, right). Across GM neurons, tortuosity distributions vary. But, each neuron presents a broad tortuosity distribution and has a mean tortuosity that is modestly greater than one (Figure 2E; the mean across all neurons is 2.1 ± 0.5). Taken together, these data suggest that GM neurons are tortuous and expansive, and somewhat variable in their macroscopic morphology.
The functional consequences of morphological complexity and variability are best considered in light of how voltage signals arising from presynaptic inputs are integrated and transformed into spike patterns. Previous work suggests that pre- and post-synaptic connections are distributed throughout the finer process within the STG neuropil (King, 1976b; Kilman and Marder, 1996) and that spike initiation zones are located in the periphery, where the axons exit the neuropil (Raper, 1979; Miller, 1980). Thus, axon number, location, and branching patterns relative to each neuron’s neurite tree may speak to how voltage signals arising from presynaptic inputs are integrated. GM neurons typically have between 3–5 axons (Figure 3A), with at least one projection to the aln and one projection to the dvn (this is consistent with previous findings; Maynard and Dando, 1974). In 2/14 neurons, only two axons were identified. This is likely due to an incomplete dye-fill or atypical branching of axonal projections beyond the microscope’s field of view. Branch order distributions varied across neurons, as did axonal branch point orders (Figure 3B). As is illustrated in Figure 3C–G, in some GM neurons all axons projected from the same branch point (as in Figure 3D), whereas in other GM neurons each axon projected from a distinct subtree (as in Figure 3G). Taken together, variable axon numbers, locations, and branching patterns suggest that each GM neuron differentially integrates voltage signals arising from varying proportions of lower and higher branch orders.
Variable glutamate responses across the neuronal structure
One might expect that complex and variable neuronal morphologies would yield complex and variable electrotonic structures. To probe the distributed cable properties of GM neurons, we employed a custom-built microscope to focally photo-uncage MNI-glutamate (Tocris) with ultraviolet light at many sites across the GM neuronal structure (in TTX to attenuate spiking activity and overall circuit activity). In each GM neuron (N = 10), we used two-electrode current clamp at the soma to measure voltage responses arising from glutamate receptor activation at the soma, primary neurite, and more distal sites (Figure 4; 7–20 photo-uncaging sites per neuron with 1 ms, 30 mW UV pulses and a bath concentration of 250 µM MNI-glutamate). Experiments probing the spatial resolution of our photo-uncaging system demonstrated that the effective photo-uncaging radius was approximately 15 µm, a resolution high enough to target individual neurites (Figure 4—figure supplement 1A–D).
Focal glutamate photo-uncaging evoked hyperpolarizing voltage deflections across the neuronal structure (Figure 4). These inhibitory potentials are consistent with ionotropic glutamate-gated inhibitory currents described in previous work (Marder and Paupardin-Tritsch, 1978; Eisen and Marder, 1982; Marder and Eisen, 1984; Cleland and Selverston, 1995). Although voltage responses were uniformly inhibitory across the neuronal structure, their amplitudes varied (Figure 4A–C). Figure 4A–C show dye-fills of three individual GM neurons with photo-uncaging sites indicated with colored circles (left) and maximal responses to focal glutamate photo-uncaging that vary in amplitude across these sites (right). These focal glutamate responses were relatively stable, desensitizing only when the inter-pulse period (IPPs) was less than 10 s (Figure 4—figure supplement 1E–G).
Across GM neurons (N = 10), a non-zero response (measured at the soma) to glutamate applied to the soma was observed in only one GM neuron. Response amplitudes varied across positions with a mean maximum amplitude of 2.15 ± 1 mV. The mean coefficient of variation of non-zero responses, across preparations, was 0.4 ± 0.1 (N = 10), suggesting that, on average, non-zero responses varied by approximately 40% of the mean response amplitude within each neuron.
The electrotonic decrement of a voltage signal, from photo-stimulation site to recording site, is dependent on a number of factors: the diameter of the neurite through which it propagates (which influences the axial resistance to current flow), the membrane resistance, the extent of branching and associated loss of conductance at branch points, and the absolute distance the signal travels (Rall, 1959; Goldstein and Rall, 1974). To quantify the dependence of response amplitude on each neuron’s cable properties, we generated ‘lolliplots’ of response amplitude, measured at the soma, to each photo-stimulated site, as a function of their distance from the somatic recording site, branch order, and neurite size (measured in terms of neurite diameter in the x-y plane; Figure 5). The lolliplots show that the greatest response amplitudes can occur from photo-stimulation at sites with any distance from the soma, of any branch order, and of any size. Linear regression analyses of response amplitudes as a function of each of these geometric parameters showed no significant linear relationships, with insignificant p-values>0.1 (Figure 5—figure supplement 1, Table 1). Taken together, these data suggest that there is no apparent dependence of somatic response amplitude on where the responses were evoked.
In each neuron, heterogeneous response amplitudes are likely a consequence of variable receptor densities and cable properties. Thus, simply measuring the maximal response amplitudes arising from photo-stimulation at each of these sites and regressing them against their varying cable properties does not negate the possibility that the measured amplitudes are a consequence of both of these factors.
Probing electrotonic structure with reversal potential measurements
A neuron’s distributed cable properties, or electrotonic structure, can be quantitatively assessed by measuring the reversal potentials (more specifically, the apparent Erevs measured at the soma) of local glutamate responses evoked at sites varying in distance from the somatic recording site (Calvin, 1969; Carnevale and Johnston, 1982). This approach distinguishes the degree to which the passive cable properties cause electrotonic decrement of voltage signals in their path of propagation to the somatic recording site, independent of the maximal conductance or receptor density at the photo-stimulation site. To demonstrate this logic, we built a library of passive cable models in the NEURON simulation platform (Hines and Carnevale, 1997; see Materials and methods) and simulated activation of local inhibitory chloride currents (actual Erev = −70 mV) at varying distances from the recording site. Apparent Erevs were measured by manipulating the membrane potential (Vm) with current injections between −8 and +2 nA at the recording site. First, we demonstrate that the apparent Erev of an inhibitory event evoked 200 µm away from the recording site (as shown in Figure 6A) is independent of its maximal conductance (or receptor density; Figure 6B–C). The apparent Erevs were measured with four maximal conductance (gmax) values: 1, 5, 10, and 50 nS. Although the response amplitudes measured at the recording site vary with gmax, all voltage responses flip their sign at −77 mV (Figure 6B). This is shown graphically in Figure 6Ci, where response amplitude (deltaV) is plotted as a function of Vm (as measured at the recording site at 0 µm). For each gmax value, the apparent Erev was calculated as the x-intercept of the linear regression of the deltaV versus Vm curve (R > 0.9 and p<0.01 in each case). The inset in Figure 6Ci and Figure 6Cii clearly show that the x-intercepts are the same for all four gmax values. The apparent Erev is independent of gmax magnitude, regardless of the activation site’s distance from the recording site (Figure 6—figure supplement 1).
Second, we show that the apparent Erev for a locally activated inhibitory current changes as a function of distance from the recording site (Figure 6D–F). The apparent Erevs were measured for inhibitory currents of the same gmax (5 nS) but with increasing distance from the recording site: 0 to 1000 µm, at increments of 200 µm (Figure 6D). It is evident from the traces in Figure 6E that the apparent Erev, or membrane potential at which the responses change sign, occurs at increasingly hyperpolarized membrane potentials as a function of increasing distance from the recording site. This hyperpolarizing shift in apparent Erev is illustrated in Figure 6Fi, where the x-intercepts for the deltaV versus membrane potential curves (as measured at the recording site at 0 µm) shift leftward with increasing activation site distance (shown clearly in the inset, a magnification of the x-intercepts). The apparent Erev shifts from an accurate measure of the actual Erev of −70 mV when activated at the recording site (0 µm) to −110 mV when activated 1000 µm away (Figure 6Fii). Thus, from considering one cable model with one set of passive properties, we can deduce that the apparent Erev of a locally-activated inhibitory current is independent of its maximal conductance (or receptor density) and is dependent on the site’s distance from the recording site.
Figure 7 demonstrates that the dependence of the apparent Erev on distance is contingent upon the effective electrotonic structure of the voltage signal’s path of propagation, from activation site to recording site. With a similar set-up to the simulation shown in Figure 6C–F, apparent Erevs were measured for responses arising from activation of the inhibitory current at sites with increasing distance from the recording site (Figure 7A). We conducted this simulation in a library of cable models with uniform diameters (5 µm) and lengths (1000 µm), but varying passive properties. Specifically, this set of 20 cable models varied in their combinations of passive leak conductances (gpas): 5, 10, 20, or 50 nS/cm2 and axial resistances (Ra): 1, 5, 10, 30, or 50 (Ω·cm). Thus, each cable model is distinguished by its electrotonic length constant (λ), determined by the expression: where . λ, in µm or mm, is equivalent to the distance at which a propagating voltage signal decrements to 37% of the maximal voltage signal (as would be measured at the site of activation). Figure 7B shows inhibitory voltage events evoked at increasing distances from the recording site (at 0 µm) in three cable models with different λ values (200, 460, and 800 µm). With distance, the apparent Erev undergoes a hyperpolarizing shift (as was the case Figure 6E). However, the rate at which this apparent Erev hyperpolarizes, as a function of distance, is dependent on the electrotonic length constant of the cable (Figure 7C). When λ = 800 µm, all apparent Erevs are −70 mV, regardless of activation site distance. When λ = 460 µm, apparent Erevs shift from −70 mV at the recording site (0 µm) to −88 mV at 1000 µm. When λ = 200 µm, the apparent Erev shifts from −70 mV at 0 µm to well below −100 mV when the activation site is 1000 µm away. It is evident that the apparent Erev undergoes a greater hyperpolarizing shift with distance with decreasing λ values. Figure 7D shows the dependence of apparent Erev as a function of activation site location for all 20 cable models with λ values ranging between 200–5000 µm. Considering a subset of these cable models, with λ values between 300–1700 µm, the apparent Erevs measured at different activation site distances diverge as λ decreases to 300 µm and converge as λ increases beyond 1000 µm (Figure 7E).
Taken together, these proof-of-principle simulations show the unmistakable relationship between electrotonic structure and apparent Erevs measured for activation sites varying in their distance from the recording site. Invariant apparent Erevs suggest a relatively high electrotonic length constant, whereas heterogeneous Erevs suggest a lower electrotonic length constant. As is shown in Figure 6, this approach to characterizing the electrotonic structure, or passive cable properties, of a neuron is independent of differential gmax values, or receptor densities, across the neuronal structure.
In this simulation paradigm, current was injected at the recording site and flowed from the recording site to the stimulation site, changing the membrane potential at the distal site. The difference between the apparent Erev and the actual Erev is indicative of the ease of current flow in this direction and ability to manipulate the membrane potential at this distal site. Thus, there will always be a discrepancy between the voltage at the soma and the voltage at the photo-uncaging site and this discrepancy will depend on the effective electrotonic length constant of the neurite path. Therefore, this assay tests the passive cable properties of the path of propagation as is most relevant to current flow from recording site to stimulation site. Even so, the observations of sizeable voltage events at the recording site and reasonable reversal potentials are suggestive of a level of electrotonic compactness that is relevant to voltage signal propagation in either direction.
Distributed reversal potentials in GM neurons are nearly invariant
Using two-electrode current clamp at the soma, we measured apparent Erevs of local inhibitory responses evoked by focal photo-uncaging of glutamate at positions varying in distance from the somatic recording site. In these experiments, current was injected at the somatic recording site and flowed centrifugally from the recording site to the photo-uncaging site, changing the membrane potential at this distal site. Figure 8 illustrates the apparent Erevs of local inhibitory responses evoked at 7–15 sites, varying in their cable properties, across the same GM neuronal structure described in Figure 4C (a second example is shown in the Figure 8—figure supplement 1). For each photo-uncaging site, response amplitudes measured at the soma were plotted as a function of somatic membrane potential (Figure 8B,C; Figure 8—figure supplement 1B and C). These data were fit with linear regression analyses (R > 0.9 in all cases) and the reversal potentials were determined by calculating the x-intercepts of the linear fits. The linear fits for each position show little variation in the reversal potential across all positions within each preparation, with a within-neuron coefficient of variation of 0.04 ± 0.01 (mean ± SD; Figure 8D; Figure 8—figure supplement 1D). Although apparent Erevs were nearly invariant within each neuron, mean reversal potentials across neurons did vary, with a pooled mean of –78.6 ± 7 mV (Table 2). This could be attributed to real differences in the ionic current, which could be carried by potassium, chloride, or a combination of the two ions (Marder and Paupardin-Tritsch, 1978; Eisen and Marder, 1982). A consequence of variable response amplitudes across positions yielded variable linear fit slopes within neurons (as shown in Figure 8D and Figure 8—figure supplement 1D). This is likely a reflection of variable receptor densities or maximal conductances at the different photo-uncaging sites (consistent with the simulations shown in Figure 6A–C). There was no significant linear relationship between input resistances measured at the soma (Table 2) and the mean within-neuron reversal potentials (p=0.2; data not shown graphically, but available in Table 2).
Lolliplots for each preparation, showing apparent Erevs for each photo-uncaging site as a function of their distance from the somatic recording site, diameter, and branch order (Figure 9), confirm no dependence of the apparent Erev on these cable properties. Linear regression analyses revealed near-zero slopes for all 10 neurons (Figure 9—figure supplement 1; Table 2). This suggests that, even though these sites vary in their absolute distance, diameter, and branch order, they do not vary substantially in their electrotonic distance from the somatic recording site. The minimal hyperpolarizing shift for apparent Erevs measured for activation sites between 100 and 800 µm from the recording site is consistent with a λ > 1.5 mm (referencing Figure 7D and E). Taken together, these results demonstrate that GM neurons are surprisingly electrotonically compact, despite their expansive structures and morphological complexity.
Discussion
Neuronal circuits function reliably despite remarkable animal-to-animal variability in the synaptic and intrinsic conductances of their constituent neurons (Goaillard et al., 2009; Norris et al., 2011; Roffman et al., 2012 ; Sakurai et al., 2014); see Calabrese et al. (2011) and Marder et al. (2015) for reviews). Given that neuronal physiology also depends on the passive cable properties arising from geometry, we examined the physiological consequences of animal-to-animal variability in neuronal morphology. Despite their expansive and complex morphologies, GM neurons have electrotonically compact structures. This effectively compensates for morphological variability and contributes to consistent neuronal and circuit-level function across animals.
Complex yet compact
Here, we present the ostensible conundrum wherein an identifiable neuron type, despite its complex, highly-branched, neurite tree, is surprisingly electrotonically compact. This result differs from studies in a variety of neuron types that attribute specific physiological computations and plasticity rules to compartmentalized electrotonic structures. Early work in insect identified neurons with stereotyped dendritic branching patterns showed that electrotonically distinct dendritic subtrees result in the weighted integration of sensory inputs (Murphey et al., 1984; Bacon and Murphey, 1984; Miller and Jacobs, 1984; Jacobs and Miller, 1985). Studies in hippocampal pyramidal neurons (Spruston and Johnston, 1992; Carnevale et al., 1997; Mainen and Sejnowski, 1996; Jaffe and Carnevale, 1999), have attributed Hebbian plasticity in part to the passive normalization of postsynaptic potentials arising from the electrotonically distant apical and basal dendritic tufts. Work in medium spiny neurons (MacAskill et al., 2012) and thalamocortical neurons (Connelly et al., 2016) have demonstrated that the activation pattern of spatially-distributed, electrotonically distant, synaptic inputs produces different neuronal and circuit-level computations. The present work differs from these other studies of neuronal electrotonus and enriches our framework for understanding how morphology maps (or does not map) to physiological function. Otherwise stated: structural complexity does not necessarily yield compartmentalized computations.
Physiological implications
In GM neurons, synaptic voltage events may propagate tortuous neurite paths that extend beyond half a millimeter in length (Figure 2). Yet, current can be injected at the somatic recording site and effectively alter the membrane potential at such distal and distant sites allowing for reasonable apparent Erevs (approximately −80 mV). Likewise, the amplitude of the voltage response at the photo-uncaging site must be larger than what we observe at the soma. Because distally evoked events, initiated as far as 1 mm away, can be observed at the somatic recording site, the electrotonic decrement of these distal events must be small enough so that apparent reversal potentials can be recorded. Lastly, the invariance in apparent Erevs across sites (a measure independent of receptor density) suggests relatively little variance in the electrotonic decrement of signals coming from disparate sites across the neurite tree. In this sense, these neurons function almost like a single compartment, despite their complex structures.
The geometric and/or passive cable properties provide the most plausible explanation for this electrotonically compact structure. Recent work has quantified the fine anatomical properties of GM neurons and other STG neuron types, showing that the primary neurite can be as large as 15–20 µm in diameter and that the most distal neurite branches may taper to diameters between <1 and 10 µm (Otopalik et al., 2017). GM neurites may have relatively low axial resistances as a consequence of large neurite diameters. It is possible that voltage signal decrement may also be minimized by high membrane resistances across neurite branches. Given that these neurons are electrotonically compact, the input resistance as measured at the soma is likely a reflection of the resistance of the membrane surface area of much of the entire neuron. Thus, it is not surprising that input resistances measured at the soma are relatively low (mean of approximately 10 MΩ (Table 2), consistent with many years of recordings from STG neurons). If the neuron were less electrotonically compact, the input resistance measured at the soma would be higher, as the measurement would be restricted to the surface area of the local, somatic membrane.
Other neuron types compensate for passive attenuation of voltage responses with distance-dependent scaling of synaptic receptor density (Andrasfalvy and Magee, 2001; Magee and Cook, 2000; Smith et al., 2003). Our experiments showed heterogeneous maximal response amplitudes across the neuronal structure. Due to the uniformity of their apparent Erevs, it is likely that these variable response amplitudes arise from local variations in receptor densities. Because the response amplitudes across sites vary in a manner that is independent of distance from the recording site (Figure 5; Figure 5—figure supplement 1; Table 1), it is unlikely that receptor densities are scaling with distance in a systematic way.
These experiments provide an ‘upper bound’ on the effective electrotonic length constant and, therefore, the compactness of GM neurons. These experiments were done in TTX, wherein TTX-sensitive, voltage-gated sodium channels are blocked and circuit activity is silenced. While the overall membrane conductance may be higher in the absence of TTX (although, no change in input resistance as measured at the soma was detected before and after TTX addition), we would not expect these TTX-sensitive currents to substantially alter inhibitory voltage signal propagation at the range of membrane potentials probed here (−120 to −40 mV). Furthermore, spike initiation zones, where TTX-sensitive channels are most likely to reside, are located in the periphery, where the axons exit the neuropil (Raper, 1979; Miller, 1980). Thus, it is unlikely that TTX-sensitive voltage-gated channels would shunt the current arising from these evoked events in the same way as has been seen in other systems (Laurent, 1990). In STG neurons, it is an intriguing possibility that separation of synaptic integration and slow waves from spike initiation zones may be a morphological strategy established to avoid shunting of synaptic currents. Future experiments in varying pharmacological and modulatory conditions could shed light on how different voltage-gated currents, modulatory currents, and ongoing synaptic input during rhythmic activity, may effectively compartmentalize these otherwise compact passive neuronal structures.
Electrotonic structure in circuit context
GM electrotonic structure is best understood in light of STG circuit architecture. The identified neurons of the STG exhibit complex morphologies (Wilensky et al., 2003; Baldwin and Graubard, 1995; Bucher et al., 2007; Goeritz et al., 2013; Otopalik et al., 2017). Like the GM neuron, all STG neurons display large somata (50–150 µm in diameter) and primary neurites that ramify throughout the STG neuropil, wherein synaptic partners form numerous, sparse synapses (King, 1976a, 1976b; Baldwin and Graubard, 1995). If each neuron type were highly electrotonically compartmentalized, yet variable across animals, wiring this circuit would be a puzzling developmental task. The fact that these neurons present electrotonically compact structures simplifies our understanding of the developmental wiring rules that may be required and how such a circuit can be relatively immune to structural differences.
The neurons of the STG rely on graded transmission and slow oscillations, rather than fast spikes, to maintain phase relationships at the circuit level (Graubard, 1978; Raper, 1979; Graubard et al., 1980, 1983; Anderson and Barker, 1981; Manor et al., 1997, 1999; Bose et al., 2014). Given the relatively slow temporal precision of this circuit, an electrotonically compact structure is sufficient for integrating activity of many, yet sparsely distributed, synaptic inputs from each presynaptic neuron (King, 1976a, 1976b), independent of synaptic site locations. In this way, electrotonic compactness both masks the observed heterogeneity in glutamate sensitivity across the neuronal structure (Figures 4 and 5) and diminishes the consequences of presumed variability in synaptic site location arising from observed animal-to-animal variability in GM morphology (also see Otopalik et al., 2017).
In this scenario, presynaptic inputs may influence GM neuron activity with equivalent efficacy regardless of synaptic site location. This synaptic democracy (Häusser, 2001) is achieved by combining an electrotonically compact structure with graded transmission resilient to electrotonic decrement across sparsely distributed, synchronous presynaptic sites. This strategy is in stark contrast to the tight tuning of receptor or ion channel distributions employed by some neuron types. For example, CA1 pyramidal neurons compensate for passive attenuation of voltage responses with distance-dependent scaling of synaptic receptor density (Andrasfalvy and Magee, 2001; Magee and Cook, 2000; Smith et al., 2003). It is feasible that strategies for achieving synaptic democracy vary across circuit contexts. The input-output computations of pyramidal neurons are typically dependent on spikes and fast voltage transients, whereas the neurons of the stomatogastric ganglion rely more heavily on graded transmission and slow oscillations to serve their circuit-level function (Graubard, 1978; Raper, 1979; Graubard et al., 1980, 1983; Anderson and Barker, 1981; Manor et al., 1997, 1999; Bose et al., 2014). Electrotonic structure may reflect the temporal precision of the neuronal and circuit-level computations performed.
Morphologies that are ‘good enough’ rather than optimal
Numerous works have argued that specific neuronal geometries are optimal for precise neuronal computations (Mainen and Sejnowski, 1996; Stiefel and Sejnowski, 2007; Cuntz et al., 2010). Experimentalists and theorists alike have suggested that neurons employ developmental growth rules that fine-tune neuronal geometry for optimal current transfer and wiring costs (Chklovskii, 2000, 2004; Chen et al., 2006; Wen and Chklovskii, 2008; Cuntz et al., 2007, 2010; Kim et al., 2012). Many of these works rely on studies in neuron types with both recognizable morphologies and known computations and/or plasticity rules. In this way, such rules hinge on a somewhat circular premise that specific neuronal functions arise from specialized geometries. As is evident in the present work, not all neuron types exhibit conserved morphologies across animals, yet show stereotyped physiological properties and circuit-level functions. We present a case in which the solution to the morphology-to-physiology transform is many-to-one.
Meaning in morphology
Neuron types can have characteristic, recognizable morphologies. Many studies have explored stereotypy in macroscopic dendritic and axonal arborization patterns in a variety of systems, including cricket (Miller and Jacobs, 1984) and grasshopper (Goodman, 1976, 1978) sensory interneurons, the insect (Cuntz et al., 2008) and mammalian (Bloomfield and Miller, 1986; Hong et al., 2011) retina, and somatosensory (Wang et al., 2002), motor (Ghosh and Porter, 1988), and visual (Martin et al., 1983; Martin and Whitteridge, 1984a, 1984b) cortices. That said, there are remarkably few instances (Cuntz et al., 2008; Wang et al., 2002) in which multiple examples of relatively complete reconstructions have been published in enough detail to judge whether the ranges of neuronal morphological features shown here, pertinent to a neuron’s cable properties, are typical or more pronounced than in other systems.
Here, we argue that the fine structural details of complex morphology may not matter for the neuronal and circuit-level function of an identified neuron type. In the STG, neurons rely predominantly on slow oscillations for circuit function, and electrotonically compact structures elegantly compensate for a high degree of animal-to-animal variability in morphology. The degree of animal-to-animal variability in neuronal morphology, and whether it is compensated for, may depend on the system and the precision of the neuronal and circuit-level computation(s) to be performed.
Materials and methods
Animals and dissections
Request a detailed protocolAdult male Jonah Crabs (Cancer borealis) were purchased from Commercial Lobster (Boston, MA) and maintained in artificial seawater at 10–13°C on a 12 hr light/12 hr dark cycle without food. On average, animals were acclimated at this temperature for one week before use. Prior to dissection, animals were anesthetized for 30 min on ice. Dissections were performed as previously described (Gutierrez and Grashow, 2009) in saline solution (440 mM NaCl, 11 mM KCl, 26 mM MgCl2, 13 mM CaCl2, 11 mM Trizma base, 5 mM maleic acid, pH 7.45). In brief, the stomach was dissected from the animal. The intact stomatogastric nervous system (STNS) was isolated from the stomach, including: the two bilateral commissural ganglia, esophageal ganglion, and stomatogastric ganglion (STG), as well as the lvn, mvn, dgn. The STNS was pinned down in a Sylgard-coated petri dish (10 mL) and continuously superfused with chilled saline.
Electrophysiology and dye-fills
Request a detailed protocolThe STG was desheathed and intracellular recordings from somata were performed with 20–30 MΩ glass microelectrodes filled with internal solution (10 mM MgCl2, 400 mM potassium gluconate, 10 mM HEPES buffer, 15 mM NaSO4, 20 mM NaCl as in Hooper et al., 2015). Intracellular signals were amplified with an Axoclamp 900A amplifier (Molecular Devices). For extracellular nerve recordings, Vaseline wells were built around the lvn, mvn, and dgn and stainless steel pin electrodes were used to monitor extracellular nerve activity (Figure 1A). Extracellular nerve recordings were amplified using model 3500 extracellular amplifiers (A-M Systems). Data were acquired using a Digidata 1440 digitizer (Axon Instruments) and pClamp data acquisition software (Axon Instruments, version 10.5). For GM identification, one of two electrodes was impaled into the soma and spiking activity was matched with GM spike units on the dgn. GM identity was verified with positive and negative current injections (Figure 1D). Following unambiguous identification, the GM soma was impaled with a second electrode containing dilute alexa488 dye (2 mM Alexa Fluor 488-hyrazide sodium salt (ThermoFisher Scientific, catalog no. A-10436, dissolved in internal solution)). The GM neuron was iontophoretically dye-filled with negative current pulses (−4 nA, 500 ms at 0.5 Hz) for 15–25 min. For two-electrode current clamp, the electrode containing alexa488 was typically used for recording and amplified with a 0.1xHS headstage. The electrode used for cell identification was used for current injection and amplified with a 1xHS headstage. Resting membrane potential and input resistance were monitored throughout the experiment to ensure the integrity of the preparation (neurons with input resistances <5 MΩ were discarded). In TTX, the mean input resistance, as measured at the somata and in the linear range of the current-voltage curve, was 9.9 ± 2.9 MΩ across preparations (Table 2). Input resistances did not change significantly before and after addition of TTX to the bath (data not shown). Reversal potentials for the glutamate response were determined by evoking at least three responses at >6 membrane potentials spanning −110 mV to −40 mV. In three preparations, an offset in the membrane voltage recording occurred during the dye-fill. This offset remained unchanged for the remainder of the experiment. This offset was corrected post-hoc during analysis of the recordings.
Focal glutamate uncaging
Request a detailed protocolFor photo-uncaging experiments, preparations were superfused with a multi-channel Ecoline re-circulating pump (Ismatec/Harvard Apparatus, catalog no. PY2 72–6432) to maintain a stable bath volume and superfusion rate. 250 µM MNI-caged-L-glutamate (dissolved in saline; Tocris Bioscience, catalog no. 1490) was bath applied. 10−7 M teterodotoxin (TTX) was also superfused to minimize spike-driven synaptic activity. Alexa488-filled GM neurons were visualized with a custom-built epifluorescence microscope (Figure 10) equipped with a 40x water-immersion UV fluorescence objective (Olympus, LUMPLFLN 40XW) and a 470 nm LED (Thor Labs, M470L2). The emitted fluorescence was imaged with a monochrome CCD camera (Scientifica, SciCam). Focal photo-activation of MNI-glutamate was achieved with a small ultraviolet (UV) spot (~10 µm in diameter; Figure 4—figure supplement 1) projected through this same 40x objective lens. By situating the custom microscope on a micromanipulator (Sutter MPC-200), this UV spot could be lased at different positions on the GM neuronal structure. Three-dimensional coordinates for each photo-uncaging site were tabulated during the experiment. Glutamate responses were evoked with 30 mW (as measured at the back aperture of the objective), 1 ms UV pulses with inter-pulse periods no less than 30 s to minimize photo-damage and desensitization (Figure 4—figure supplement 1). To achieve this small UV spot, the UV laser beam (DPSS Lasers, model no. 35-07–100, 0.96 W, 100 kHz repetition rate) was coupled to a 50 µm diameter fiber optic cable (Thor Labs, M50L02S-A) with a UV lens (Thor Labs) and collimated with a 50 mm focal length plano-convex lens (Thor Labs, LA4148-UV). This collimated UV beam was delivered into the 40x objective lens to produce a focused spot of UV light on the preparation. A neutral density filter wheel was situated in the beam path for manipulation of the beam intensity (Thor Labs, NDM2). To ensure consistent beam intensities across experiments, the intensity was monitored in real-time using a photodiode (Thor Labs, PDA25K) previously calibrated with a power meter (Thor Labs, S302C). For precise temporal control of UV stimuli, a shutter (Thor Labs, SH1) was situated in the beam path. Both the laser Q-Switch and shutter were triggered by a set of coupled model 2100 isolated pulse stimulators (A-M Systems). The effective photo-uncaging radius in the x-y plane (15 µm) was determined by photo-uncaging at peripheral, distal neurites and moving away at 5 µm increments in the x-y plane (Figure 4—figure supplement 1) in a number of neuron types: pyloric dilater (PD), lateral pyloric (LP), and GM neurons.
Dye-fill amplification and immunohistochemistry
Request a detailed protocolFollowing photo-uncaging experiments, GM neurons were secondarily dye-filled with 2% Lucifer Yellow CH dipotassium salt (LY; Sigma, catalog no. L0144) in filtered water using a low-resistance electrode (10–15 MΩ). LY was injected for 20–50 min with negative current pulses (−6 to −8 nA, 500 ms at 0.5 Hz). Once fine neurites of the cell could be visualized with a fluorescent stereomicroscope (Leica MF165 F), a preliminary image was acquired at 11.5x magnification with an attached monochrome digital camera (Leica DFC365 FC). LY-filled preparations were fixed for 40 min at 21°C or overnight at 4°C in 2% paraformaldehyde in phosphate-buffered saline (PBS; 440 mM NaCl, 11 mM KCl, 10 mM Na2HPO4, 2 mM KH2PO4, pH 7.4). Preparations were washed with 0.1 M PBS-T ((0.1–0.3%% Triton X-100 in PBS) and stored in PBS for 0–3 days prior to immunohistochemistry. The LY signal was amplified by 16 hr incubation with a polyclonal rabbit anti-LY antibody (1:500; Molecular Probes). After washing 5 × 15 min in PBS-T at room temperature, preparations were incubated in a secondary Alexa Fluor-488-conjugated goat-anti-rabbit antibody (1:500; Molecular Probes) for 1.5 hr at room temperature. Preparations were washed 5 × 15 min in PBS at room temperature before mounting on pre-cleaned slides (25 × 75 × 1 mm, superfrost, VWR) in Vectashield (Vector Laboratories, Burlingame, CA), with 9 mm diameter, 0.12 mm depth silicone seal spacers (Electron Microscopy Sciences, Hatfield, PA) under #1.5 coverslips (Fisher Scientific). Mounting in Vectashield with a spacer was sufficient to maintain the 3-dimensional structure of the neuron and ganglion (as in Goeritz et al., 2013).
Confocal imaging and 3D reconstructions
Request a detailed protocolConfocal stacks of the LY-filled neurons were acquired with a SP2 Leica Microscope and Leica Application Suite Advanced Fluorescence (LAS AF) software. Image stacks were acquired with a 20x dry objective (Leica HC PL APO CS 20x) at 1024 × 1024 resolution in 0.5 µm steps. Image stacks were visualized in both FIJI (ImageJ) software and KNOSSOS 3D image visualization and annotation software (developed by teams at the Heidelberg University and Karlsruhe Institute of Technology, employed by the Max Planck Institute for Medical Research, and freely distributed at: http://www.knossostool.org/). KNOSSOS software was used to manually trace and generate skeletons of the GM neuronal structures in three dimensions (Supplementary file 1-Neuronal Structures Hoc files; as in Otopalik et al., 2017). It is important to note that the resolution used here allowed us to reconstruct many neurons, but resulted in smaller total cable lengths and branch point numbers than reported in Otopalik et al. (2017), where neuronal dye-fills were imaged at 60x magnification. This higher magnification would have precluded completion of the reconstructions and photo-uncaging experiments reported here.
Morphological analysis
Request a detailed protocolFollowing 3D skeleton generation, morphological analyses were completed across all GM neurons (n = 14), using a suite of custom analysis scripts written in Python using the iPython command line (freely available at: https://python.org and https://ipython.org, respectively) by AS. For each skeleton, branch points, lengths, and orders were measured in reference to the soma and used to generate dendrogram representations with normalized path lengths. Branch, or path, lengths were measured as the most direct neurite path from the soma to each branch tip. Tortuosities were calculated for each path length, as the ratio of the path length over the Euclidean distance from soma to branch tip. Axon locations were identified as the last branch points without terminating branch tips. Photo-uncaging positions were re-located on preliminary fluorescence images of the Lucifer yellow dye-fill (as situated during experiment, at 11.5x magnification) using a custom alignment script written in MATLAB (Mathworks, version 2015b) by AO, and then manually re-located in the 3D skeleton using KNOSSOS (Supplementary file 2-Uncaging Coordinates Hoc files). The branch order, path length, and diameter were determined for each photo-uncaging site, based on the confocal image stack. All quantitative morphology analysis scripts are freely available at the Marder Lab GitHub website (https://github.com/marderlab/Quantifying_Morphology).
Electrophysiology analysis
Request a detailed protocolRecordings acquired using Clampex software (pClamp Suite by Molecular Devices, version 10.5) and were visualized offline using a MATLAB waveform analysis toolbox written by Ted Brookings and analyzed with custom MATLAB scripts written by AO. Briefly, this pipeline of analysis scripts was used to detect and browse evoked glutamate responses, measure voltage response amplitudes and membrane potentials, plot raw recordings and processed data, and perform some statistical analyses. To determine reversal potentials at a given photo-uncaging position, raw response amplitudes were plotted as a function of membrane potential. These data were fit with a linear function (with an R value > 0.9 in all cases). In some cases, response amplitudes saturated and the linear regression analyses were performed only in the linear range of these data. The x-intercept of the resulting linear function indicated the reversal potential at that position. All electrophysiology analysis scripts are available at the Marder lab GitHub (https://github.com/marderlab).
Passive cable models
Request a detailed protocolA library of equivalent cylinder models was constructed in NEURON (freely available at: https://www.neuron.yale.edu/neuron), to simulate apparent reversal potential (Erev) measurements of inhibitory responses generated at positions varying in distance from the recording site. The cables were uniform in their geometric properties: 1000 µm in length and 5 µm in diameter. All cables had a membrane capacitance of 1 µF·cm−2. However, their passive properties were varied combinatorically with four passive conductance (gpas) values (5, 10, 20, and 50 nS/cm2) and five axial resistance (Ra) values (1, 5, 10, 30, 50 Ω·cm). Consequently, each cable can be distinguished by its electrotonic length constant (λ in µm) such that, where . λ values ranged between 200 µm and 5 mm (Figure 11). To measure the apparent Erevs at the 0 µm end of the cable, the membrane potential was manipulated at the recording site with current injections between −8 and +2 nA. An inhibitory current (actual Erev = −70 mV, 𝜏 = 3 ms, gmax = 1, 5, 10, or 50 nS) was activated at sites varying in distance from the recording site (at 0 µm). For a given cable and activation site, the apparent Erev was calculated by plotting the response amplitude (deltaV) as a function of membrane potential at the recording site. For each curve, linear regression analysis was completed (R > 0.9 and p<0.01 in all cases). The apparent Erevs were identified as the x-intercept of the linear fit.
References
-
Signal delay and input synchronization in passive dendritic structuresJournal of Neurophysiology 70:2066–2151.
-
Synaptic mechanisms that generate network oscillations in the absence of discrete postsynaptic potentialsJournal of Experimental Zoology 216:187–191.https://doi.org/10.1002/jez.1402160121
-
Distance-dependent increase in AMPA receptor number in the dendrites of adult hippocampal CA1 pyramidal neuronsThe Journal of Neuroscience 21:9151–9160.
-
Receptive fields of cricket giant interneurones are related to their dendritic structureThe Journal of Physiology 352:601–623.https://doi.org/10.1113/jphysiol.1984.sp015312
-
Distribution of fine neurites of stomatogastric neurons of the crab cancer borealis: evidence for a structured neuropilThe Journal of Comparative Neurology 356:355–367.https://doi.org/10.1002/cne.903560304
-
A functional organization of ON and OFF pathways in the rabbit retinaThe Journal of Neuroscience 6:1–13.
-
The role of linear and voltage-dependent ionic currents in the generation of slow wave oscillationsJournal of Computational Neuroscience 37:229–242.https://doi.org/10.1007/s10827-014-0498-4
-
Neuronal morphology and neuropil structure in the stomatogastric ganglion of the lobster, homarus americanusThe Journal of Comparative Neurology 501:185–205.https://doi.org/10.1002/cne.21169
-
Coping with variability in small neuronal networksIntegrative and Comparative Biology 51:845–855.https://doi.org/10.1093/icb/icr074
-
Electrophysiological characterization of remote chemical synapsesJournal of Neurophysiology 47:606–621.
-
Comparative electrotonic analysis of three classes of rat hippocampal neuronsJournal of Neurophysiology 78:703–723.
-
Optimal sizes of dendritic and axonal arbors in a topographic projectionJournal of Neurophysiology 83:2113–2119.
-
Glutamate-gated inhibitory currents of central pattern generator neurons in the lobster stomatogastric ganglionThe Journal of Neuroscience 15:6631–6639.
-
Optimization principles of dendritic structureTheoretical Biology and Medical Modelling 4:21.https://doi.org/10.1186/1742-4682-4-21
-
One rule to grow them all: a general theory of neuronal branching and its practical applicationPLoS Computational Biology 6:e1000877.https://doi.org/10.1371/journal.pcbi.1000877
-
The morphological identity of insect dendritesPLoS Computational Biology 4:e1000251.https://doi.org/10.1371/journal.pcbi.1000251
-
Mechanisms underlying pattern generation in lobster stomatogastric ganglion as determined by selective inactivation of identified neurons. III. synaptic connections of electrically coupled pyloric neuronsJournal of Neurophysiology 48:1392–1415.
-
Functional consequences of animal-to-animal variation in circuit parametersNature Neuroscience 12:1424–1430.https://doi.org/10.1038/nn.2404
-
Ionic currents of the lateral pyloric neuron of the stomatogastric ganglion of the crabJournal of Neurophysiology 67:318–331.
-
Isogenic grasshoppers: genetic variability in the morphology of identified neuronsThe Journal of Comparative Neurology 182:681–705.https://doi.org/10.1002/cne.901820408
-
Compensation for variable intrinsic neuronal excitability by circuit-synaptic interactionsJournal of Neuroscience 30:9145–9156.https://doi.org/10.1523/JNEUROSCI.0980-10.2010
-
Graded synaptic transmission between identified spiking neuronsJournal of Neurophysiology 50:508–529.
-
Synaptic transmission without action potentials: input-output properties of a nonspiking presynaptic neuronJournal of Neurophysiology 41:1014–1025.
-
Cancer borealis stomatogastric nervous system dissectionJournal of Visualized Experiments 1207.https://doi.org/10.3791/1207
-
BookDynamical Biological Networks: The Stomatogastric Nervous SystemCambridge, MA: MIT Press.
-
Motor patterns in the stomatogastric ganglion of the lobster panulirus argusThe Journal of Experimental Biology 62:405–420.
-
Synaptic function: dendritic democracyCurrent Biology 11:R10–R12.https://doi.org/10.1016/S0960-9822(00)00034-8
-
Gastric mill activity in the lobster. I. spontaneous modes of chewingJournal of Neurophysiology 59:528–550.
-
The NEURON simulation environmentNeural Computation 9:1179–1209.https://doi.org/10.1162/neco.1997.9.6.1179
-
Stereotyped axonal arbors of retinal ganglion cell subsets in the mouse superior colliculusThe Journal of Comparative Neurology 519:1691–1711.https://doi.org/10.1002/cne.22595
-
Cell dialysis by sharp electrodes can cause nonphysiological changes in neuron propertiesJournal of Neurophysiology 114:1255–1271.https://doi.org/10.1152/jn.01010.2014
-
Passive normalization of synaptic integration influenced by dendritic architectureJournal of Neurophysiology 82:3268–3285.
-
Interpretation of voltage-clamp measurements in hippocampal neuronsJournal of Neurophysiology 50:464–486.
-
Ultrastructure of the stomatogastric ganglion neuropil of the crab, cancer borealisThe Journal of Comparative Neurology 374:362–375.https://doi.org/10.1002/(SICI)1096-9861(19961021)374:3<362::AID-CNE5>3.0.CO;2-#
-
Geometric theory predicts bifurcations in minimal wiring cost trees in biology are flatPLoS Computational Biology 8:e1002474.https://doi.org/10.1371/journal.pcbi.1002474
-
The role of single neurons in information processingNature Neuroscience 3 Suppl:1171–1177.https://doi.org/10.1038/81444
-
Voltage-dependent nonlinearities in the membrane of locust nonspiking local interneurons, and their significance for synaptic integrationThe Journal of Neuroscience 10:2268–2280.
-
Dendritic computationAnnual Review of Neuroscience 28:503–532.https://doi.org/10.1146/annurev.neuro.28.061604.135703
-
Subcellular connectivity underlies pathway-specific signaling in the nucleus accumbensNature Neuroscience 15:1624–1626.https://doi.org/10.1038/nn.3254
-
Temporal dynamics of graded synaptic transmission in the lobster stomatogastric ganglionThe Journal of Neuroscience 17:5610–5621.
-
Network oscillations generated by balancing graded asymmetric reciprocal inhibition in passive neuronsThe Journal of Neuroscience 19:2765–2779.
-
Transmitter identification of pyloric neurons: electrically coupled neurons use different transmittersJournal of Neurophysiology 51:1345–1361.
-
Variability, compensation and homeostasis in neuron and network functionNature Reviews Neuroscience 7:563–574.https://doi.org/10.1038/nrn1949
-
Robust circuit rhythms in small circuits arise from variable circuit components and mechanismsCurrent Opinion in Neurobiology 31:156–163.https://doi.org/10.1016/j.conb.2014.10.012
-
Variability, compensation, and modulation in neurons and circuitsPNAS 108 Suppl 3:15542–15548.https://doi.org/10.1073/pnas.1010674108
-
Physiological and morphological properties of identified basket cells in the cat's visual cortexExperimental Brain Research 50:193–200.https://doi.org/10.1007/bf00239183
-
Form, function and intracortical projections of spiny neurones in the striate visual cortex of the catThe Journal of Physiology 353:463–504.https://doi.org/10.1113/jphysiol.1984.sp015347
-
The structure of the stomatogastric neuromuscular system in Callinectes Sapidus, homarus americanus and panulirus argus (Decapoda crustacea)Philosophical Transactions of the Royal Society B: Biological Sciences 268:161–220.https://doi.org/10.1098/rstb.1974.0024
-
Effects of maintained depolarization of presynaptic neurons on inhibitory transmission in lobster neuropilJournal of Comparative Physiology 97:215–243.https://doi.org/10.1007/BF00617544
-
Relationships between neuronal structure and functionThe Journal of Experimental Biology 112:129–145.
-
BookMechanisms Underlying Pattern Generation in the Lobster Stomatogastric GanglionSan Diego: University of California.
-
Neurospecificity in the cricket cercal systemThe Journal of Experimental Biology 112:7–25.
-
Constancy and variability in the output of a central pattern generatorJournal of Neuroscience 31:4663–4674.https://doi.org/10.1523/JNEUROSCI.5072-10.2011
-
Similar network activity from disparate circuit parametersNature Neuroscience 7:1345–1352.https://doi.org/10.1038/nn1352
-
Branching dendritic trees and motoneuron membrane resistivityExperimental Neurology 1:491–527.https://doi.org/10.1016/0014-4886(59)90046-9
-
Distinguishing theoretical synaptic potentials computed for different soma-dendritic distributions of synaptic inputJournal of Neurophysiology 30:1138–1168.
-
Time constants and electrotonic length of membrane cylinders and neuronsBiophysical Journal 9:1483–1508.https://doi.org/10.1016/S0006-3495(69)86467-2
-
Handbook of Physiology: The Nervous System. Cellular Biology of Neurons39–97, Core conductor theory and cable properties of neurons, Handbook of Physiology: The Nervous System. Cellular Biology of Neurons, Bethesda, MA, American Physiological Society.
-
Convergent rhythm generation from divergent cellular mechanismsJournal of Neuroscience 33:18047–18064.https://doi.org/10.1523/JNEUROSCI.3217-13.2013
-
Animal-to-animal variability of connection strength in the leech heartbeat central pattern generatorJournal of Neurophysiology 107:1681–1693.https://doi.org/10.1152/jn.00903.2011
-
Neural basis of teeth coordination during gastric mill rhythms in spiny lobsters, panulirus interruptusThe Journal of Experimental Biology 114:99–119.
-
Organization of the stomatogastric ganglion of the spiny lobsterJournal of Comparative Physiology 91:33–51.https://doi.org/10.1007/BF00696155
-
The stomatogastric nervous system: structure and function of a small neural networkProgress in Neurobiology 7:215–289.https://doi.org/10.1016/0301-0082(76)90008-3
-
Mechanism of the distance-dependent scaling of schaffer collateral synapses in rat CA1 pyramidal neuronsThe Journal of Physiology 548:245–258.https://doi.org/10.1113/jphysiol.2002.036376
-
Voltage- and space-clamp errors associated with the measurement of Electrotonically remote synaptic eventsJournal of Neurophysiology 70:781–802.
-
Perforated patch-clamp analysis of the passive membrane properties of three classes of hippocampal neuronsJournal of Neurophysiology 67:508–537.
-
Mapping function onto neuronal morphologyJournal of Neurophysiology 98:513–526.https://doi.org/10.1152/jn.00865.2006
-
Propagation of action potentials in dendrites depends on dendritic morphologyJournal of Neurophysiology 85:926–937.
-
Neurons that form multiple pattern generators: identification and multiple activity patterns of gastric/pyloric neurons in the crab stomatogastric systemJournal of Neurophysiology 65:111–122.
-
A cost-benefit analysis of neuronal morphologyJournal of Neurophysiology 99:2320–2328.https://doi.org/10.1152/jn.00280.2007
-
Stereotyped neuropil branching of an identified stomatogastric motor neuronThe Journal of Comparative Neurology 466:554–563.https://doi.org/10.1002/cne.10903
Article and author information
Author details
Funding
National Institute of Neurological Disorders and Stroke (F31NS092126)
- Adriane G Otopalik
National Institute of Neurological Disorders and Stroke (R37NS017813)
- Eve Marder
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank: Frank Mello for assistance in constructing mechanical components of rig; Bernardo Sabatini for optics expertise while constructing the custom microscope; Matthew Stenerson and Richard Ho for manual tracing of neuronal dye-fills; Philipp Rosenbaum for completion of several technical experiments; Cosmo Guerini for generating additional analytical tools in Python; Edward Dougherty and the Confocal Imaging Lab at Brandeis University.
Copyright
© 2017, Otopalik et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,696
- views
-
- 596
- downloads
-
- 39
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Neuronal physiology depends on a neuron’s ion channel composition and unique morphology. Variable ion channel compositions can produce similar neuronal physiologies across animals. Less is known regarding the morphological precision required to produce reliable neuronal physiology. Theoretical studies suggest that moraphology is tightly tuned to minimize wiring and conduction delay of synaptic events. We utilize high-resolution confocal microscopy and custom computational tools to characterize the morphologies of four neuron types in the stomatogastric ganglion (STG) of the crab Cancer borealis. Macroscopic branching patterns and fine cable properties are variable within and across neuron types. We compare these neuronal structures to synthetic minimal spanning neurite trees constrained by a wiring cost equation and find that STG neurons do not adhere to prevailing hypotheses regarding wiring optimization principles. In this highly modulated and oscillating circuit, neuronal structures appear to be governed by a space-filling mechanism that outweighs the cost of inefficient wiring.
-
- Developmental Biology
- Neuroscience
Two major ligand-receptor signaling axes – endothelin Edn3 and its receptor Ednrb, and glial-derived neurotrophic factor (GDNF) and its receptor Ret – are required for migration of enteric nervous system (ENS) progenitors to the hindgut. Mutations in either component cause colonic aganglionosis, also called Hirschsprung disease. Here, we have used Wnt1Cre and Pax2Cre in mice to show that these driver lines label distinct ENS lineages during progenitor migration and in their terminal hindgut fates. Both Cre lines result in Hirschsprung disease when combined with conditional Ednrb or conditional Ret alleles. In vitro explant assays and analysis of lineage-labeled mutant embryos show that GDNF but not Edn3 is a migration cue for cells of both lineages. Instead, Edn3-Ednrb function is required in both for GDNF responsiveness albeit in different ways: by expanding the Ret+ population in the Pax2Cre lineage, and by supporting Ret function in Wnt1Cre-derived cells. Our results demonstrate that two distinct lineages of progenitors give rise to the ENS, and that these divergently utilize endothelin signaling to support migration to the hindgut.