Neuronal morphologies built for reliable physiology in a rhythmic motor circuit
Abstract
It is often assumed that highly-branched neuronal structures perform compartmentalized computations. However, previously we showed that the Gastric Mill (GM) neuron in the crustacean stomatogastric ganglion (STG) operates like a single electrotonic compartment, despite having thousands of branch points and total cable length >10 mm (Otopalik et al., 2017a; 2017b). Here we show that compact electrotonic architecture is generalizable to other STG neuron types, and that these neurons present direction-insensitive, linear voltage integration, suggesting they pool synaptic inputs across their neuronal structures. We also show, using simulations of 720 cable models spanning a broad range of geometries and passive properties, that compact electrotonus, linear integration, and directional insensitivity in STG neurons arise from their neurite geometries (diameters tapering from 10-20 µm to < 2 µm at their terminal tips). A broad parameter search reveals multiple morphological and biophysical solutions for achieving different degrees of passive electrotonic decrement and computational strategies in the absence of active properties.
https://doi.org/10.7554/eLife.41728.001Introduction
Neurons often present complex and highly-branched morphologies. How synaptic voltage events propagate within and across neurite branches is determined by the structure’s geometrical and biophysical properties. Passive voltage propagation is influenced by the neurite’s diameter, membrane resistance, and axial resistance (Rall, 1959; Rall, 1960; Rall, 1969; Jack et al., 1975; Holmes, 1989). Rall and colleages were the first to apply passive cable theory and develop the equivalent cylinder model for the study of electrotonus in single dendrites and branched dendritic trees (Rall, 1959; Rall, 1960; Rall, 1969). Schierwagen then devised a broader mathematical description of membrane voltage distributions in complex, highly-branched neurite trees with non-uniform boundary conditions and geometries (Schierwagen, 1989). These seminal theoretical studies have provided an invaluable framework for understanding passive neuronal physiology. However, in a functioning neuron, synaptic inputs, receptors, and ion channels may shunt or amplify propagating voltage signals (London and Häusser, 2005). Thus, the transformation from neuronal morphology to electrophysiological activity patterns is often unpredictable in the absence of direct experimental assessment. To date, measuring voltage attenuation across the many neurite paths presented in complex neuronal structures using electrophysiological techniques has proven difficult. Thus, electrotonus has been experimentally assessed in only a handful of neuron types (for example: Spruston and Johnston, 1992; Spruston et al., 1994; Rapp et al., 1994; Carnevale et al., 1997; Stuart and Spruston, 1998; Chitwood et al., 1999; Jaffe and Carnevale, 1999; Otopalik et al., 2017b; Medan et al., 2018), and this greatly restricts our understanding of the breadth of biophysical organizations utilized in different neuron types and circuit contexts.
In two recent studies, we characterized the morphology (Otopalik et al., 2017a) and passive electrophysiology (Otopalik et al., 2017b) of the identified neurons of the crustacean stomatogastric ganglion (STG), a small central pattern-generating circuit mediating the rhythmic contractions of the animal’s foregut. The 14 identified neurons of the STG present distinct, cell-type-specific electrophysiological waveforms, firing patterns and circuit functions (Harris-Warrick et al., 1992). We quantified numerous morphological features pertaining to the macroscopic branching patterns and fine cable properties of four neuron types (Otopalik et al., 2017a). Interestingly, the four neuron types did not adhere to optimal wiring principles (Cuntz et al., 2010) or Rall’s 3/2 rule (Rall, 1959) and exhibited expansive neurite trees that sum to >10 mm of total cable length, tortuous and long individual branches (ranging between 100 µm and 1 mm in length) and thousands of branch points with complex geometries. There was quantifiable inter-animal variability in many features within neuron types, and no single metric or combination of metrics distinguished the four neuron types.
In a second study, we then asked: How do STG neurons produce reliable firing patterns across animals, given their apparently inefficient and highly variable structures? As a first examination of how neuronal morphology maps to physiology in the STG, we characterized electrotonus, or passive voltage signal propagation, in one STG neuron type, the Gastric Mill (GM) neuron. We were surprised to find that, despite their expansive and complex neuronal structures, GM neurons are relatively electrotonically compact and operate much like single electrical compartments (Otopalik et al., 2017b). We suggested that compact electrotonic structures may effectively counteract the potential physiological consequences of morphological variability observed in GM neurons across animals (Otopalik et al., 2017a; Otopalik et al., 2017b).
We first motivate the present study with an empirical numerical model that measures electrotonus in a library of cable models with varying geometrical and passive properties. In doing so, we recapitulate passive electrotonic decrement described in the aforementioned seminal theoretical studies (Rall, 1959; Rall, 1960; Rall, 1969). Yet, we also show that a subset of cable models with geometries consistent with those observed in multiple STG neuron types (Otopalik et al., 2017a), are relatively resilient to electrotonic decrement. Given this prediction, we asked whether compact electrotonus is a generalizable feature in the STG, and therefore common among multiple, distinct neuron types. Using glutamate photo-uncaging in tandem with intracellular electrophysiology (as in Otopalik et al., 2017b), we measure electrotonus in four STG neuron types. We then investigate how neurite geometry shapes voltage integration in these complex neuronal structures. By complementing these experiments with validating computational simulations, we demonstrate that STG neurite geometries and passive properties are sufficient to account for the compact electrotonus and voltage integration observed in this pattern-generating circuit. Furthermore, our numerical simulations suggest that different neurite geometries may be suitable for circuits that subserve different functions.
Results
Simulating electrotonus in diverse neurites
In previous work, we were surprised to find that STG neurons exhibit neurite lengths as long as 1 mm and neurite diameters between 10–20 µm at their primary neurite junctions, that decrement to <1 µm at their terminating tips (Otopalik et al., 2017a), suggesting geometries with a 1–2% taper. Passive cable theory suggests that voltage signals passively propagating such long distances are likely to undergo a great deal of attenuation (Rall, 1960; Rall, 1964; Rall, 1969; Jack et al., 1975; Schierwagen, 1989). Experimental studies to date have validated predictions presented in these theoretical studies in dendrites with uniform and/or smaller diameters (typically <3 µm; Holmes, 1989; Jaffe and Carnevale, 1999; Stuart and Spruston, 1998). These studies suggested that diameter is a critical parameter in determining electrotonic decrement in the absence of amplifying or shunting mechanisms. Yet, none of these studies examined electrotonus in neurites with the wide diameters exhibited by STG neurites. Thus, we first asked whether there may be a boundary at which a neurite’s diameter is wide enough to overcome distance-dependent attenuation.
We conducted a proof-of-concept computational characterization of electrotonus in a library of 720 passive cable models with a much broader range of geometries (Figure 1A) and passive properties (specific Ra values between 50–300 Ω x cm and specific Rm values between 1,000–20,000 Ω x cm2; see Materials and methods) than explored in earlier studies. These geometries range from narrow cables with uniform diameters (Figure 1A, gray) to broad cables with diameters reminiscent of STG neurites (Figure 1A, blue), thereby spanning the broad range of neurite geometries existing in diverse nervous systems. Electrotonus was characterized by measurement of an effective electrotonic length constant (λeffective in µm; see Materials and methods), which is equivalent to the distance at which a voltage signal decrements to 37% of the maximal voltage amplitude (at the activation site). Thus, greater λeffective values are suggestive of less electrotonic decrement. Figure 1B illustrates the measurement of λeffective in a classic cable model with a uniform diameter of 0.5 µm, Ra = 100 Ω x cm, Rm = 10,000 Ω x cm2. In this case, the voltage event attenuates greatly with distance (λeffective is approximately 300 µm). This voltage attenuation is robust and λeffective is less than 1 mm for a broad range of passive properties (Figure 1C). Altering the morphology of this cable model to reflect the geometry of an STG neurite (Figure 1D) results in a smaller response amplitude at the site of activation (compare maximum amplitudes in Figure 1B and D). However, this is accompanied by a robust increase in λeffective and similar signal amplitudes across the cable. Larger λeffective values are observed for a broad range of passive parameters (Figure 1E). Ra values between 50–150 Ω x cm and Rmvalues > 10,000 Ω x cm2 yield λeffective values greater than 1 mm. Examining the entire morphological and biophysical parameter space (Figure 1—figure supplement 1) demonstrates that a range of λeffective values can be achieved across neurites varying not only in their passive properties, but also in their geometries. Neurites with wide proximal diameters present relatively long λeffective values even in the presence of a large load at the proximal end of the cable (d0), made to mimic a putative shunt imposed by the rest of the neurite tree (Figure 1—figure supplement 2). Importantly, these simulations demonstrate that long λeffective values can be achieved as a consequence of neurite geometry alone, in the absence of voltage-gated ion channels or other amplifying mechanisms that have been implicated in boosting distally-evoked events in other neuron types (Magee and Cook, 2000; Andrasfalvy and Magee, 2001; Smith et al., 2003; Gulledge et al., 2005; Lavzin et al., 2012).
Linking electrotonus and neurite geometry in STG neurons
Characterization of electrotonus in our cable model library suggests that the wide diameters of STG neurites may effectively equalize passive voltage signal propagation. Consistent with this prediction, we recently demonstrated that the Gastric Mill neuron exhibits a relatively uniform electrotonic structure (Otopalik et al., 2017b). To determine if this is a generalizable feature across different STG neuron types, we characterized electrotonus in multiple neurites of three additional STG neuron types: Lateral Pyloric (LP), Ventricular Dilator (VD), and Pyloric Dilator (PD) neurons, while also corroborating our previous findings in GM neurons (N = 5–6 neurons of each type). These four neuron types present equally complex and expansive morphologies, but distinct voltage waveforms and circuit functions (Figure 2). PD, LP, VD, and GM neurons were unambiguously identified by their innervation patterns (Figure 2A) and by matching their intracellular spiking patterns with concurrent extracellular recordings of nerves known to contain their axons (Figure 2B). PD and LP neurons innervate two muscles in the pylorus of the foregut (Figure 2A) and participate in the ongoing, triphasic pyloric rhythm (Figure 2B). PD and LP can both be identified by matching their intracellular firing patterns with spiking units on the lateral ventricular nerve (lvn; Figure 2B). The VD neuron innervates the cv1 muscle of the pylorus and can be identified on the medial ventricular nerve (mvn; Figure 2B). The GM neuron participates in the episodic gastric mill rhythm, innervates gm1, 2, and 3 muscles (Figure 2A), and can be identified on the dorsal gastric nerve (dgn; Figure 2B). When filled with fluorescent dye, each neuron presents highly branched and expansive neurite trees (Figure 2C). The morphological features of these neuron types have been described quantitatively and in detail in previous studies (Wilensky et al., 2003; Bucher et al., 2007; Thuma et al., 2009; Otopalik et al., 2017a). To examine neurite geometries in these four neuron types, we completed volumetric reconstructions and continuous measurement of neurite diameters (from soma to terminating tip) for 23 neurite paths, for which we had confocal image stacks with sufficient resolution. Recapitulating previous work, we observed neurite paths that taper from 10 to 20 µm at the soma-primary neurite junction to <2 µm at terminating tips (Figure 2D–E, Figure 2—figure supplement 1–2 to Figure 2; Otopalik et al., 2017a). Geometrical taper from proximal to distal end is thought to increase the electrotonic length constant for long neurite paths (Holmes and Rall, 1992). Yet, we found that, while some STG neurites tapered gradually, others presented diameters that decrease in an abrupt step. Previous study has demonstrated that geometrical taper influences action potential shape and velocity (Goldstein and Rall, 1974). But, it is unclear whether this geometrical feature influences the electrotonic decrement of slower inhibitory potentials. Thus, we simulated and measured λeffective (as in Figure 1 and Figure 1—figure supplement 1–2 to Figure 1) in cables with varying degrees of taper from proximal to distal end (Figure 2—figure supplement 3 to Figure 2, Part A). We ran this simulation in cables with varying proximal diameters (ranging between 1 and 10 µm) and maintained the same 80% reduction in diameter from proximal to distal end (Figure 2—figure supplement 3 to Figure 2, Part B). In brief, we found that cables with a gradual taper presented longer λeffective than those with abrupt step reductions (consistent with Holmes and Rall, 1992). Yet, cables with wide diameters consistent with those measured in STG neurons (Figure 2E), presented λeffective> 1 mm for even abrupt step-reductions in diameter for a range of passive properties (Ravalues < 100 Ω x cm and Rmvalues > 10,000 Ω x cm2). This contrasts with finer cables (tapering from 1 to 2 µm to sub-µm diameters), wherein any step reductions along the path of propagation resulted in λeffective < 0.5 mm.
Measuring electrotonus in four STG neuron types
We experimentally assessed passive voltage signal propagation by evoking inhibitory potentials at numerous sites on the neurite tree with focal photo-uncaging of MNI-glutamate and two-electrode current clamp recordings at the soma (as in Otopalik et al., 2017b). Figure 3A shows example traces of evoked inhibitory potentials at six sites on an individual PD neurite. When the somatic membrane potential is at rest (approximately −50 mV) evoked events at the six sites are inhibitory but vary in magnitude. Two-electrode current clamp was used to manipulate the somatic membrane potential (between −40 and −100 mV) and the apparent reversal potentials (Erevs) for events evoked at each site were determined by plotting response amplitude (measured at the soma) as a function of somatic membrane potential (Figure 3B). The x-intercept of the linear fit of these data serves as our measure of the apparent Erev for each site. Across the six sites evaluated in Figure 3A and B, apparent Erevs ranged between −59 and −70 mV.
Maximal response amplitudes (as measured at a somatic membrane potential of −50 mV) and apparent Erevs were measured for numerous sites across the neurite trees of each neuron type (for 10–30 distinct sites across the neurite trees of five LP, VD, GM neurons and six PD neurons; Tables 1 and 2). These data are summarized in Figure 3C and D, where maximal response amplitudes and apparent Erevs for individual sites are plotted as a function of distance from the somatic recording site for each neuron. For comparison across multiple neurons, maximal response magnitudes (which varied between 0–4 mV) were normalized to the minimum (always 0 mV) and maximum response amplitudes within each neuron (Figure 3—figure supplement 1–4 to Figure 3 show raw maximum response amplitudes and apparent Erevs for individual PD, LP, VD, and GM neurons, respectively). Figure 3C shows that, across all neuron types, maximal response amplitudes show no quantitative trend as a function of distance, nor is there any evidence of normalization of response amplitude with distance (such that response amplitudes are uniform across sites). This is reflected in poorly-fit and insignificant linear regressions of these data (Table 1).
Figure 3D shows normalized apparent Erevs as a function of distance from the somatic recordings site. For comparison of apparent Erevs evoked at sites on multiple neurons, apparent Erevs were normalized to and plotted as a percent of the mean apparent Erev across sites within each neuron. Horizontal lines denote 0.05, or 5%, above and below the mean Erev. For PD, LP, and GM neurons, there appears to be no substantial hyperpolarization in apparent Erevs with distance from the somatic recording site (this is validated with linear regressions shown in Table 1). Figure 3B shows that GM neurons present exceptionally invariant apparent Erevs across sites on their neurite trees (mean coefficient of variance (CV) within individual GM neurons was 0.04; Table 2). This result is consistent with previous findings (Otopalik et al., 2017b). Likewise, LP, and PD neurons exhibit relatively invariant apparent Erevs (mean CVs were: 0.06, 0.05, and 0.08, respectively; Table 2). It should be noted that a subset of VD neurons shows higher standard deviations and CVs, suggesting heterogeneity of voltage signal propagation across the neuronal structure. Yet, statistical comparison of CVs (ANOVA, [F(3, 17)=1.2 p=0.341]) and mean apparent Erevs (ANOVA, [F(3, 17)=2.29, p=0.1154]) revealed no statistically significant differences across neuron types. This suggests that the four neuron types are similarly electrotonically compact neuronal structures, wherein apparent Erevs typically varied by <10% of the mean within individual neurons. This translates to a range of mean apparent Erev ± 6–8 mV within each individual neuron.
Directional sensitivity and voltage integration in diverse neurites
Next, we examined (i) the directional sensitivity of voltage signal propagation and (ii) how multiple voltage events are integrated in single neurites varying in their geometry and passive properties. We first characterized these two properties in our library of cable models with varying geometries and passive properties. Figure 4 illustrates the measurement of directional sensitivity and summation arithmetic in two cable models: a classic cable model with a narrow diameter (0.5 µm) and a cable model reminiscent of an STG neurite (with an axial diameter tapering from 20 µm to 0.5 µm; as in Figure 1).
Inhibitory events of equivalent conductance magnitude and kinetics (see Materials and methods) were evoked at 500, 600, 700, 800, and 900 µm from the recording site individually and in sequence (5 Hz) in the inward and outward directions (Figure 4A–D).
Directional bias was calculated as the integral of the inward response minus the integral of the outward response (Figure 4E and F, top). Thus, positive directional biases are indicative of an inward bias in voltage signal propagation (toward the recording electrode) and negative directional biases are indicative of an outward bias in voltage signal propagation. Figure 4E and F illustrate a notable difference in the directional sensitivity of passive voltage propagation in these two cable models; the classic cable model shows a robust inward bias and the STG neurite shows no directional sensitivity. These biases persist across a broad range of passive properties. However, the inward bias in the classic cable model does attenuate with increasing axial resistivity (Ra) and decreasing membrane resistivity (Rm).
The arithmetic of voltage integration was characterized by comparing the combined response to the arithmetic sum of the individual events (example traces depicted in Figure 4C and D, bottom). Integration is described as sublinear, linear, or supralinear if the combined response is less than, equal to, or larger than the predicted arithmetic sum. We calculated the arithmetic, or linearity, of the evoked responses as the integral of the inward response minus the integral of the combined individual responses (Figure 4E and F, bottom). The classic cable presents sublinear integration of inhibitory voltage events. Yet, as the Rm decreases and Ra increases, voltage summation transitions from sublinear to linear. In contrast, the STG neurite presents relatively linear integration across a broad range of passive properties. Examination of directional sensitivity (Figure 4—figure supplement 1) and summation arithmetic in the entire cable library (Figure 4—figure supplement 2) demonstrates that a broad range of computations can arise as a consequence of neurite morphology. Taken together, these simulations predict that STG neurites will present directional insensitivity and near-linear voltage integration if they operate predominantly by passive propagation and in the absence of active properties serving to shunt or amplify voltage signals.
Direction insensitivity in STG neurites
To test for directional bias in voltage signal propagation, sequential voltage events were evoked at multiple sites within the same secondary branches (from tip to primary neurite junction; Figure 5Ai). The integrals of the summed responses for inward and outward activation were measured at the soma (Figure 5Aii). Directional preference for each branch was assessed by plotting the response integrals for the inward and outward directions against each other and comparison with the identity line, which is indicative of direction insensitivity wherein the inward and outward response integrals are equal (Figure 5Aiii). Any branches with points left of the identity line present an inward, or centripetal bias, whereas any points right of the identity line are suggestive of a centrifugal, or outward, bias. Interestingly, all four neuron types show little directional selectivity (this is supported by root-mean-square error values (RMSE) <0.5 mV*s, a measure of goodness-of-fit to the identity line). Example traces and direction selectivity plots for each cell type can be found in Figure 5—figure supplement 1–4 to Figure 5. Taken together, these results suggest that neurites in each of these four STG neuron types do not exhibit the directional selectivity that has been described in other neuron types (Barlow and Levick, 1965; Euler et al., 2002; London and Häusser, 2005; Branco et al., 2010).
Arithmetic of voltage signal integration
If voltage signal propagation in STG neurons is predominantly shaped by passive properties and tapered neurite geometries, and less so by active biophysical properties, which may shunt or amplify propagating voltage signals, we would expect to observe linear voltage summation as predicted by our simulation (Figure 4 and Figure 4—figure supplement 2). To assess the arithmetic of voltage summation, we calculated the arithmetic sum of responses evoked at individual sites across single neurites (Figure 6Ai–iii; offset by 200 ms to mimic a 5 Hz sequential activation rate). The integrals of the measured response and expected arithmetic sum for a given branch were plotted against each other and compared with the identity line, which is indicative of linear summation, wherein the measured and expected response integrals are equal (Figure 6Aiii). Any branches with points left of the identity line present sublinear summation, whereas any points right of the identity line are suggestive of supralinear summation. Figure 6B shows these plots for more than 20 branches for each neuron type. Across all neuron types, the majority of branches showed linear summation. RMSE values were less than 1.5 mV*s; thus, the measured response integrals were within 1.5 mV*s of the integral expected of linear summation. GM neurons (Figure 6Biv) exhibited particularly uniform linear summation across all branches evaluated and this is reflected in a small RMSE value of 0.32 mV*s. LP branches and PD branches show slightly higher RMSE values (greater than 1 mV*s), perhaps providing some evidence of variable shunting or amplifying mechanisms. Even so, RMSE values less than 2 mV*s across all cell types suggests relatively linear summation across most branches evaluated. Voltage summation arithmetic was similarly linear in the centrifugal direction (consult Figure 5—figure supplement 1–4 to Figure 5). It should be noted that these experiments were performed in the absence of TTX-sensitive sodium channels, which have been abolished with 10−7 M TTX in the bath.
Discussion
Compact computing in STG neurons
In the present study, we find that multiple STG neuron types present electrotonically compact structures and within-neurite voltage summation that is relatively linear and directionally insensitive. Taken together, these findings suggest that STG neurons are built to linearly sum and unify synaptic inputs distributed across their expansive and complex structures, rather than perform distributed computations on a branch-by-branch, or subtree-by-subtree basis. Our computational simulations suggest that this biophysical architecture may be achieved passively and as a simple consequence of neurite geometry in these neurons, which present neurites with wide diameters that taper from 10 to 20 µm near the soma, to sub-micron diameters at their terminating tips.
This computing strategy stands in contrast to that which has been observed in other neuron types with similarly complex neuronal structures. When voltage signals propagate long distances in the absence of amplifying mechanisms, electrotonic decrement is thought to result in sublinear voltage integration and directional bias within single branches (Rall, 1964; Gulledge et al., 2005; London and Häusser, 2005). This has been demonstrated experimentally in a variety of vertebrate neuron types in different circuit contexts, from cortex to retina (Barlow and Levick, 1965; Cash and Yuste, 1998; Euler et al., 2002; Poirazi et al., 2003; Polsky et al., 2004; Branco et al., 2010). Linear and supralinear integration are thought to require amplifying mechanisms, such as voltage-gated ion channels or distance-dependent scaling of receptors (Magee and Cook, 2000; Andrasfalvy and Magee, 2001; Smith et al., 2003; Gulledge et al., 2005; Lavzin et al., 2012). Collectively, these studies support the notion that highly-branched and expansive neuronal structures, are likely to present some degree of passive voltage attenuation (in the absence of amplifying mechanisms) and perform compartmentalized computations as a consequence. Yet, this framework for dendritic computation does not rely on the study of, or account for, the diverse dendrite and neurite morphologies observed in nature.
To make sense of our findings in the STG, we revisited basic cable theory and simulated electrotonus and passive voltage integration in a library of neurites exhibiting diverse geometries and passive biophysical properties. We observed a surprising range of electrotonic decrement, directional bias, and voltage summation arithmetic across the surveyed morphological and biophysical parameter space. Thus, this simulation revealed multiple morphological and biophysical solutions for achieving varying degrees of passive electrotonic decrement and different computational strategies. We also show that linear summation and directional insensitivity, as observed in four STG neuron types, can be achieved in the absence of active properties altogether and in the face of potential heterogeneity of passive biophysical properties. Altogether, this proof-of-concept simulation both recapitulates the principles describe in early theoretical studies and also demonstrates that different biophysical and morphological strategies are likely to be utilized by different neuron types, to suit their unique physiological function. Thus, the widely-accept principles derived from the study of so-called canonical neuron types may be less general than previously thought.
A general solution for reliable pacemaking physiology
Our findings are interesting in light of the synaptic organization of the STG. Synaptic sites are sparsely distributed throughout the neuropil and pre- and post-synaptic sites are closely apposed on the same neurites (King, 1976a; King, 1976b; Kilman and Marder, 1996). Thus, sufficiently large synaptic potentials may originate anywhere on the neurite tree, propagate in any direction, and achieve consistent neuronal output.
The motor rhythm generated by the STG relies on slow oscillations and graded inhibitory transmission (Eisen and Marder, 1982; Marder and Eisen, 1984; Maynard and Walton, 1975; Graubard et al., 1980; Manor et al., 1997; Manor et al., 1999; Bose et al., 2014; Golowasch et al., 2017). The biophysical features we have described in STG neurons may allow for the averaging of very large and slow synaptic conductances in space and time, thereby sustaining pattern generation in this motor circuit. Rhythmic networks that compute with slow oscillations and/or graded transmission (Walsh et al., 1972; Wilson and Wachtel, 1974; Pearson and Fourtner, 1975; Robertson and Pearson, 1985; Angstadt and Calabrese, 1991; Dicaprio, 1989; Dicaprio et al., 1997; Dale, 1995; DiCaprio, 2003; Smarandache-Wellmann et al., 2013) may benefit from a similar biophysical and morphological architecture. Moreover, we suggest that wide neurites may aid STG neurons in generating consistent physiological output, in the face of variable conductance magnitudes and subcellular distributions of intrinsic and synaptic properties across animals (Prinz et al., 2004; Schulz et al., 2006; Marder and Goaillard, 2006; Goaillard et al., 2009; Marder, 2011).
In the present study, we assessed the integration of inhibitory glutamate events in the presence of TTX, which blocks voltage-gated sodium currents and silences circuit activity. Thus, we can only speculate how action potentials, arising from such TTX-sensitive currents, may influence the integration of graded inhibitory glutamatergic synaptic events in the intact circuit. We would not expect TTX-sensitive currents to substantially alter inhibitory voltage signal propagation or integration at the range of membrane potentials probed here (−100 to −40 mV for apparent reversal potentials measurements and −50 mV for the voltage summation and directional sensitivity measurements). Moreover, spike initiations zones, where such TTX-sensitive channels are most likely to reside, are located just outside the neuropil, where the axons exit the ganglion (Raper, 1979; Miller, 1980). PD neurons exhibit a second, dopamine-sensitive axonal spike initiation zone between the upper dvn and its split into the bilateral lvns (refer to Figure 2A; Bucher et al., 2003). Given the peripheral locations of the spike initiation zones in these neurons, 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). Spatial separation of spike initiation zones from synaptic integration and slow wave generation in the neuropil may reduce shunting of synaptic currents. A thorough investigation of this possibility requires experiments in varying pharmacological and modulatory conditions, which may reveal how various currents and ongoing circuit activity may influence these voltage propagation and integration in these otherwise compact passive neuronal structures.
Taken together with previous work (Otopalik et al., 2017a; Otopalik et al., 2017b), the present study suggests that, given the relatively compact electrotonic architecture of STG neurons arising from their neurite geometries, other features of their morphologies and the exact spatial organization of synaptic contacts may not be critical determinants of the distinct physiological waveforms and firing patterns exhibited by different STG neuron types. Thus, one is left with the expectation that their different activity patterns are predominantly determined by their cell-type-specific ion channel and receptor expression profiles. And, indeed, there is much evidence that the 14 different neuron types express different palettes of receptors (Swensen et al., 2000; Swensen and Marder, 2000; Swensen and Marder, 2001) and cell-type-specific ratios of ion channels (Schulz et al., 2006; Schulz et al., 2007).
Of course, one is left wondering why these neurons present such expansive and complex morphologies, if they instead act as single electrotonic compartments. One possibility is that this allows STG neurons to make appropriate synaptic contacts wherever partner neurons ramifying throughout the neuropil find each other. Additionally, it is likely that STG neurons grow to fill distinct spatial fields in the neuropil (Otopalik et al., 2017a) for the purpose of maximizing surface area for the reception of the many diffuse neuromodulatory substances that are released in the hemolymph and by descending modulatory inputs of the stomatogastric nerve (Marder and Bucher, 2007; Blitz and Nusbaum, 2011).
Materials and methods
Animals and dissections
Request a detailed protocolWild-caught adult male Jonah Crabs (Cancer borealis) were acquired and maintained by the Marine Resources Center at the Marine Biological Laboratories in Woods Hole, MA. Animals were maintained on a 12 hr dark/12 hr light cycle without food and in chilled natural seawater (10–13 deg C) in a 2000-liter tank at a density of no more than 30 crabs per tank. STG dissections were executed as in Otopalik et al. (2017b) and 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.4–7.6). The intact stomatogastric nervous system, including: two bilateral commissural ganglia, esophageal ganglion, and stomatogastric ganglion (STG), as well as the lvn, mvn, dgn were dissected from the animal’s foregut and pinned down in a Sylgard-coated petri dish (10 mL). The preparation was continuously superfused with chilled saline (11–13 degrees C) for the duration of the experiment using a bipolar temperature control system (Harvard Apparatus, CL-100).
Electrophysiology and Dye-fills
Request a detailed protocolAll electrophysiology and dye-fill methods are consistent with those utilized in Otopalik et al. (2017b). The STG was desheathed for access to somata for intracellular recordings. These recordings were executed with glass micropipettes (20–30 MΩ) filled with internal solution: 10 mM MgCl2, 400 mM potassium gluconate, 10 mM HEPES buffer, 15 mM NaSO4, 20 mM NaCl (Hooper et al., 2015). Intracellular recordings signals were amplified with an Axoclamp 900A amplifier (Molecular Devices, as described in Otopalik et al. (2017b). For extracellular nerve recordings, Vaseline wells were built around the lvn, mvn, and dgn nerves and stainless-steel pin electrodes were used to monitor extracellular nerve activity (as indicated in Figure 1A). Extracellular nerve recordings were amplified using a Model 3500 extracellular amplifier (A-M Systems). All recordings were acquired with a Digidata 1550 (Molecular Devices) digitizer and visualized with pClamp data acquisition software (Axon Instruments, version 10.7). Neuron types were identified by matching concurrent intracellular spiking patterns with units on nerves known to contain their axons (as in Figure 1B) and verified with positive and negative current injections. After identification, a single neuron was filled with dilute alexa488 dye (2 mM Alexa Fluor 488-hyrazide sodium salt (ThermoFisher Scientific, catalog no. A-10436, dissolved in internal solution)) with negative current pulses (−4 nA, 500 ms at 0.5 Hz) for 15–25 min. Following the dye-fill, input resistance was measured at the soma in two-electrode current clamp (neurons with input resistances < 5 MΩ were discarded). For two-electrode current clamp, the electrode containing dilute alexa488 was used for recording and amplified on a 0.1xHS headstage. The electrode used for cell identification was used for current injection and amplified with a 1xHS headstage. Input resistance was measured throughout the experiment and neurons with input resistances < 5 MΩ were discarded. Reversal potentials for glutamate-evoked responses were determined by evoking responses at >eight membrane potentials between −100 and −40 mV. In some experiments, neurons were filled with 2% Lucifer Yellow CH dipotassium salt (LY; Sigma, catalog no. L0144; diluted in filtered water) for post-hoc imaging. LY was injected with a low-resistance (10–15 MΩ) glass micropipette for 20–50 min with negative current pulses (−six to −8 nA, 500 ms at 0.5 Hz).
Focal glutamate Photo-uncaging
Request a detailed protocolFocal glutamate photo-uncaging methodology was consistent with the methods used in Otopalik et al. (2017b), although different instrumentation was used. For photo-uncaging experiments, preparations were superfused with a re-circulating peristaltic pump to maintain a stable bath volume. 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. Alexa Fluor 488-filled neurons were visualized and focal photo-uncaging was achieved with a Laser Applied Stimulation and Uncaging (LASU) system (Scientifica). In brief: this system was composed of an epifluorescence microscope (SliceScope, Scientifica) equipped with a 4x magnification air and 40x magnification water-immersion objective lenses (Olympus; PLN 4X and LUMPLFLN 40XW, respectively). A 780 IR-LED was used to visualize the stomatogastric ganglion and locate neurons of interest. A white fluorescence illumination system (CoolLED) and FITC/Alexa Fluor 488/Fluo3/Oregon Green filter set (Chroma) were used to excite and visualize fluorescent emission from neuronal dye-fills. Images were captured with a monochrome CCD camera (Scientifica, SciCam Pro; 1360 × 1025 array and 6.54 µm2 pixel size). Focal photo-activation of MNI-glutamate was achieved with a 405 nm laser (1 ms pulses, 35 mW, spot size <1.5 µm with the 40x objective). The preparation platform and micromanipulators were mounted on a motorized movable base plate, allowing for smooth re-positioning (in the X-Y plane) of the objective over different neurites. For photo-activation at multiple sites within the field of view, the laser spot was re-positioned in quick sequence (5 Hz) using a set of X-Y galvanometers (Cambridge Technology, 6251H). Photo-uncaging sites within the field of view, laser pulse duration, and pulse rate were selected with the assistance of the LASU system software (Scientifica).
Electrophysiology analysis
Request a detailed protocolElectrophysiological responses to focal glutamate photo-uncaging were visualized and analyzed offline, as previously described (Otopalik et al., 2017b), using a set of custom MATLAB (Mathworks, version 2017b) scripts that will be made available on the Marder Lab GitHub (https://github.com/marderlab) (Otopalik, 2019; copy archived at https://github.com/elifesciences-publications/Otopalik-Pipkin-Marder-2019). Maximal response amplitudes, directional bias, and voltage integration arithmetic were calculated as an average of 3 trials of photo-stimulation at individual sites or of individual branches, accordingly. Apparent reversal potentials (Erevs) were calculated for individual sites by plotting response amplitudes as a function of somatic membrane potential and fitting these scatter plots with a linear regression. The apparent Erev for each site was calculated as the x-intercept of this fit. These methods are described in detail in Otopalik et al. (2017b). The For the linear regressions for all sites depicted in Figure 3D and supplements (plots in Part C), R-values were >0.85.
Post-hoc imaging and morphological analysis
Request a detailed protocolFor each neuron, all photo-uncaging sites were re-located in fluorescence images of the Alexa Fluor 488 and/or Lucifer Yellow dye-fills acquired with the LASU microscope (at 40x and 20x magnification; described above). The distance between each site and the somatic recording site was measured utilizing a combination of 2-D and 3-D image stacks with the assistance of Simple Neurite Tracer on ImageJ/FIJI (Longair et al., 2011). Although most experiments were exclusively conducted using the LASU microscope, a subset of neurons were imaged at 40x magnification (Zeiss C-apochromat40x/1.2 W) on a VIVO microscope system equipped with a Yokogawa (CSU-X1, Japan) spinning disk confocal scan head mounted on a Zeiss Examiner microscope. Fluorescent dye-fills were visualized with standard GFP filters and images were captured with a Prime CMOS camera (Photometrics, 95B). Multiple image stacks in the z-dimension, spanning the STG, were stitched together with the assistance of the Stitching tool in ImageJ/FIJI (Preibisch et al., 2009). Maximum projections of these stacks are shown in Figure 2C.
Volumetric neurite reconstructions
Request a detailed protocolFor high-resolution imaging and volumetric neurite reconstructions, a subset of preparations were fixed in 4% paraformaldehyde following photo-uncaging experiments and kept at 4° C until immuno-staining. For immuno-staining, each preparation was incubated with a rabbit anti-Lucifer Yellow IgG antibody (Invitrogen A5750) overnight at room temperature in phosphate-buffered saline (PBS) containing 0.1% Tris, then washed and subsequently incubated with a goat anti-rabbit IgG secondary antibody conjugated to Alexa-488 (Invitrogen A11034) in PBS at room temperature. Preparations were mounted on slides using VectaShield (Vector Laboratories H-1000). Seven preparations survived this full sequence, from electrophysiology and glutamate photo-uncaging, through immuno-staining, mounting, and imaging: one PD, three LPs, two VDs, and one GM.
High-resolution confocal images of each preparation were taken on a Leica SP5 system. We used a 63x objective to collect a montage of Z-stacks spanning the entire neuropilar arbor of the labeled cell at a voxel resolution of 0.114 μm x 0.114 μm in the X-Y plane and 0.797 μm in Z dimension. We found this resolution was sufficient to relocate and trace the branches that were subjected to glutamate photo-uncaging. Image montages were aligned and reconstructions acquired within the TrakEM2 (fiji.sc/TrakEM2, RRID: SCR_008954; Cardona et al., 2012) environment of ImageJ/FIJI. We fully traced each neuron from the origin of the primary neurite near the soma to the tips of each branch. As our analyses concerned only the branches visited for photo-uncaging, our reconstructions therefore do not recapitulate the entirety of the neuronal arbor. We traced each preparation in two ways. First, we generated a skeleton wherein the neuron is represented by a tree of connected nodes. Second, we traced the full volume of each branch containing the skeleton by manually coloring in a region in each z-layer that contained that branch.
Using these two reconstructions for each neurite, we measured the path distance of each node in the skeleton to the soma along with the cross-sectional area of the branch at each node. From these cross-sectional areas, we calculate an approximated diameter of the branch assuming the branch were cylindrical; in most cases, of course, the branches are not perfectly cylindrical. To generate cross-sectional area, we used the built-in coding environment of FIJI to write a custom Python script. Briefly, the algorithm generates a list of all points on the edge of branches in the arbor, then finds those which lie within a small threshold (0.5 μm) of a plane orthogonal to the arbor of the branch at each node. The number of these points are further reduced by excluding those which lie outside a certain distance of the node (which corresponds to the maximum expected radius of the branch at that node). For most branches, we used a distance threshold of 3 μm, while for thicker branches we used 11 μm. Finally, the projection of these points onto the orthogonal plane through the branch at each node defines a polygon of which the area is calculated using the shoelace formula.
Passive cable models
Request a detailed protocolThe library of cable models utilized in Otopalik et al. (2017b) was adapted to explore electrotonus, directional bias, and voltage integration in neurites with varying geometries (depicted in Figure 6) and passive properties: six specific axial resistances (Ra): 10, 50, 100, 150, 200, 300 Ω*cm; and six specific membrane resistances (Rm): 2 × 104, 1.6 × 104, 1 × 104, 5 × 103, 1 × 103, 1 × 102 Ω*cm2. All cables were 1000 µm in length and had a membrane capacitance of 1 µF*cm−2. Electrotonus was assessed as described in Figure 1 and voltage summation was assessed by simulating our experimental glutamate photo-uncaging procedure using the simulation platform NEURON (Hines and Carnevale, 2001). All possible combinations of neurite geometries, membrane resistances, and axial resistances were assessed in Figure 1—figure supplement 1 to Figure 1, Figure 4, and Figure 4—figure supplement 1and 2 to Figure 4, whereas a representative cross-section of cable models are presented in Figure 1—figure supplement 2 to Figure 1 and Figure 2—figure supplement 3 to Figure 2. To assess voltage integration directional bias and arithmetic, voltage was recorded 100 µm from the proximal end (d0) and inhibitory potentials (Erev = −75 mV, 𝜏 = 70 ms, gmax = 10 nS) were evoked at five sites with increasing distance from the recording site (as depicted in Figure 4A,B). Sites were activated individually, and then at 5 Hz in the inward (toward the recording site) or outward (away from the recording site) directions. Using MATLAB (Mathworks, version 2018b), integrals were calculated for the inward and outward summed responses and the expected linear sum of the individual events in either direction. Directional bias was calculated as the inward integral minus the outward integral. Linearity was calculated as the expected arithmetic sum minus the recorded voltages sum for the inward direction. Custom scripts written to generate cable models and simulate experimental procedures in NEURON were composed in Sublime Text (Sublime HQ Pty Ltd, Sydney) and can be found on the Marder Lab GitHub (https://github.com/marderlab), along with all simulation output.
Data availability
All computational scripts used for: analysis, visualization, and model simulations are available on the Marder Lab GitHub site (https://github.com/marderlab/Otopalik-Pipkin-Marder-2019), where they are freely available to the public. All source data (electrophysiology recordings and raw images) are publicly available on Dryad (https://dx.doi.org/10.5061/dryad.48pt6jd).
-
Dryad Digital RepositoryData from: Neuronal morphologies built for reliable physiology in a rhythmic motor circuit.https://doi.org/10.5061/dryad.48pt6jd
References
-
Calcium currents and graded synaptic transmission between heart interneurons of the leechThe Journal of Neuroscience 11:746–759.https://doi.org/10.1523/JNEUROSCI.11-03-00746.1991
-
The mechanism of directionally selective units in rabbit's retinaThe Journal of Physiology 178:477–504.https://doi.org/10.1113/jphysiol.1965.sp007638
-
Neural circuit flexibility in a small sensorimotor systemCurrent Opinion in Neurobiology 21:544–552.https://doi.org/10.1016/j.conb.2011.05.019
-
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
-
Axonal dopamine receptors activate peripheral spike initiation in a stomatogastric motor neuronThe Journal of Neuroscience 23:6866–6875.https://doi.org/10.1523/JNEUROSCI.23-17-06866.2003
-
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
-
Electrophysiological characterization of remote chemical synapsesJournal of Neurophysiology 47:606–621.https://doi.org/10.1152/jn.1982.47.4.606
-
Comparative electrotonic analysis of three classes of rat hippocampal neuronsJournal of Neurophysiology 78:703–720.https://doi.org/10.1152/jn.1997.78.2.703
-
Input summation by cultured pyramidal neurons is linear and position-independentThe Journal of Neuroscience 18:10–15.https://doi.org/10.1523/JNEUROSCI.18-01-00010.1998
-
Passive electrotonic properties of rat hippocampal CA3 interneuronesThe Journal of Physiology 515:743–756.https://doi.org/10.1111/j.1469-7793.1999.743ab.x
-
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
-
Experimentally derived model for the locomotor pattern generator in the xenopus embryoThe Journal of Physiology 489:489–510.https://doi.org/10.1113/jphysiol.1995.sp021067
-
Nonspiking interneurons in the ventilatory central pattern generator of the shore crab, carcinus maenasThe Journal of Comparative Neurology 285:83–106.https://doi.org/10.1002/cne.902850108
-
Maintenance of motor pattern phase relationships in the ventilatory system of the crabThe Journal of Experimental Biology 200:963–974.
-
Nonspiking and spiking proprioceptors in the crab: nonlinear analysis of nonspiking TCMRO afferentsJournal of Neurophysiology 89:1826–1836.https://doi.org/10.1152/jn.00978.2002
-
Functional consequences of animal-to-animal variation in circuit parametersNature Neuroscience 12:1424–1430.https://doi.org/10.1038/nn.2404
-
A balance of outward and linear inward ionic currents is required for generation of slow-wave oscillationsJournal of Neurophysiology 118:1092–1104.https://doi.org/10.1152/jn.00240.2017
-
Synaptic integration in dendritic treesJournal of Neurobiology 64:75–90.https://doi.org/10.1002/neu.20144
-
Cancer borealis stomatogastric nervous system dissectionJournal of visualized experiments : JoVE.https://doi.org/10.3791/1207
-
BookDynamical Biological Networks: The Stomatogastric Nervous System (CambridgeMA: MIT Press).
-
NEURON: a tool for neuroscientistsThe Neuroscientist 7:123–135.https://doi.org/10.1177/107385840100700207
-
Electrotonic length estimates in neurons with dendritic tapering or somatic shuntJournal of Neurophysiology 68:1421–1437.https://doi.org/10.1152/jn.1992.68.4.1421
-
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.https://doi.org/10.1152/jn.1999.82.6.3268
-
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-#
-
Dendritic computationAnnual Review of Neuroscience 28:503–532.https://doi.org/10.1146/annurev.neuro.28.061604.135703
-
Temporal dynamics of graded synaptic transmission in the lobster stomatogastric ganglionThe Journal of Neuroscience 17:5610–5621.https://doi.org/10.1523/JNEUROSCI.17-14-05610.1997
-
Network oscillations generated by balancing graded asymmetric reciprocal inhibition in passive neuronsThe Journal of Neuroscience 19:2765–2779.https://doi.org/10.1523/JNEUROSCI.19-07-02765.1999
-
Understanding circuit dynamics using the stomatogastric nervous system of lobsters and crabsAnnual Review of Physiology 69:291–316.https://doi.org/10.1146/annurev.physiol.69.031905.161516
-
Transmitter identification of pyloric neurons: electrically coupled neurons use different transmittersJournal of Neurophysiology 51:1345–1361.https://doi.org/10.1152/jn.1984.51.6.1345
-
Variability, compensation and homeostasis in neuron and network functionNature Reviews Neuroscience 7:563–574.https://doi.org/10.1038/nrn1949
-
Effects of maintained depolarization of presynaptic neurons on inhibitory transmission in lobster neuropilJournal of Comparative Physiology ? A 97:215–243.https://doi.org/10.1007/BF00617544
-
Differential processing in modality-specific mauthner cell dendritesThe Journal of Physiology 596:667–689.https://doi.org/10.1113/JP274861
-
BookMechanisms Underlying Pattern Generation in the Lobster Stomatogastric GanglionSan Diego: University of California.
-
Nonspiking interneurons in walking system of the cockroachJournal of Neurophysiology 38:33–52.https://doi.org/10.1152/jn.1975.38.1.33
-
Computational subunits in thin dendrites of pyramidal cellsNature Neuroscience 7:621–627.https://doi.org/10.1038/nn1253
-
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
-
Membrane potential transients and membrane time constant of motoneuronsExperimental Neurology 2:503–532.https://doi.org/10.1016/0014-4886(60)90029-7
-
Neural Theory and Modeling: Proceedings of the 1962 Ojai SymposiumTheoretical Significance of Dendritic Trees for Neuronal Input-Output Relations., Neural Theory and Modeling: Proceedings of the 1962 Ojai Symposium, Palo Alto: Stanford University Press.
-
Time constants and electrotonic length of membrane cylinders and neuronsBiophysical Journal 9:1483–1508.https://doi.org/10.1016/S0006-3495(69)86467-2
-
Physiology, morphology and detailed passive models of guinea-pig cerebellar purkinje cellsThe Journal of Physiology 474:101–118.https://doi.org/10.1113/jphysiol.1994.sp020006
-
Neural circuits in the flight system of the locustJournal of Neurophysiology 53:110–128.https://doi.org/10.1152/jn.1985.53.1.110
-
A non-uniform equivalent cable model of membrane voltage changes in a passive dendritic treeJournal of Theoretical Biology 141:159–179.https://doi.org/10.1016/S0022-5193(89)80015-3
-
Five types of nonspiking interneurons in local pattern-generating circuits of the crayfish swimmeret systemJournal of Neurophysiology 110:344–357.https://doi.org/10.1152/jn.00079.2013
-
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
-
Perforated patch-clamp analysis of the passive membrane properties of three classes of hippocampal neuronsJournal of Neurophysiology 67:508–529.https://doi.org/10.1152/jn.1992.67.3.508
-
Determinants of voltage attenuation in neocortical pyramidal neuron dendritesThe Journal of Neuroscience 18:3501–3510.https://doi.org/10.1523/JNEUROSCI.18-10-03501.1998
-
GABA and responses to GABA in the stomatogastric ganglion of the crab cancer borealisThe Journal of Experimental Biology 203:2075–2092.
-
Modulators with convergent cellular actions elicit distinct circuit outputsThe Journal of Neuroscience 21:4050–4058.https://doi.org/10.1523/JNEUROSCI.21-11-04050.2001
-
Pyloric neuron morphology in the stomatogastric ganglion of the lobster, panulirus interruptusBrain, Behavior and Evolution 73:26–42.https://doi.org/10.1159/000202988
-
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 (R35NS097343)
- Eve Marder
Grass Foundation
- Adriane G Otopalik
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Jennifer Bestman for assistance in spinning disk and confocal microscopy; the Marine Resources Center at the Marine Biological Laboratories for acquiring and maintaining animals; Louie Kerr at the Central Microscopy Facility; Dana Mock-Munoz de Luna for administrative support; Kamran Kodhakhah, Heather Rhodes, and the 2017 Grass Fellows for their support and feedback; and lastly, Edward Dougherty at the Brandeis University Confocal Imaging Lab for support and microscope maintenance. This study was funded by the Grass Foundation and NINDS awards to F31NS092126 to AO and R35NS097343 to EM.
Copyright
© 2019, 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
-
- 1,657
- views
-
- 236
- downloads
-
- 21
- 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
Long-range axonal projections of diverse classes of neocortical excitatory neurons likely contribute to brain-wide interactions processing sensory, cognitive and motor signals. Here, we performed light-sheet imaging of fluorescently labeled axons from genetically defined neurons located in posterior primary somatosensory barrel cortex and supplemental somatosensory cortex. We used convolutional networks to segment axon-containing voxels and quantified their distribution within the Allen Mouse Brain Atlas Common Coordinate Framework. Axonal density was analyzed for different classes of glutamatergic neurons using transgenic mouse lines selectively expressing Cre recombinase in layer 2/3 intratelencephalic projection neurons (Rasgrf2-dCre), layer 4 intratelencephalic projection neurons (Scnn1a-Cre), layer 5 intratelencephalic projection neurons (Tlx3-Cre), layer 5 pyramidal tract projection neurons (Sim1-Cre), layer 5 projection neurons (Rbp4-Cre), and layer 6 corticothalamic neurons (Ntsr1-Cre). We found distinct axonal projections from the different neuronal classes to many downstream brain areas, which were largely similar for primary and supplementary somatosensory cortices. Functional connectivity maps obtained from optogenetic activation of sensory cortex and wide-field imaging revealed topographically organized evoked activity in frontal cortex with neurons located more laterally in somatosensory cortex signaling to more anteriorly located regions in motor cortex, consistent with the anatomical projections. The current methodology therefore appears to quantify brain-wide axonal innervation patterns supporting brain-wide signaling.
-
- Neuroscience
- Structural Biology and Molecular Biophysics
More than 20 recurrent missense gain-of-function (GOF) mutations have been identified in the sodium-activated potassium (KNa) channel gene KCNT1 in patients with severe developmental and epileptic encephalopathies (DEEs), most of which are resistant to current therapies. Defining the neuron types most vulnerable to KCNT1 GOF will advance our understanding of disease mechanisms and provide refined targets for precision therapy efforts. Here, we assessed the effects of heterozygous expression of a Kcnt1 GOF variant (Kcnt1Y777H) on KNa currents and neuronal physiology among cortical glutamatergic and GABAergic neurons in mice, including those expressing vasoactive intestinal polypeptide (VIP), somatostatin (SST), and parvalbumin (PV), to identify and model the pathogenic mechanisms of autosomal dominant KCNT1 GOF variants in DEEs. Although the Kcnt1Y777H variant had no effects on glutamatergic or VIP neuron function, it increased subthreshold KNa currents in both SST and PV neurons but with opposite effects on neuronal output; SST neurons became hypoexcitable with a higher rheobase current and lower action potential (AP) firing frequency, whereas PV neurons became hyperexcitable with a lower rheobase current and higher AP firing frequency. Further neurophysiological and computational modeling experiments showed that the differential effects of the Kcnt1Y777H variant on SST and PV neurons are not likely due to inherent differences in these neuron types, but to an increased persistent sodium current in PV, but not SST, neurons. The Kcnt1Y777H variant also increased excitatory input onto, and chemical and electrical synaptic connectivity between, SST neurons. Together, these data suggest differential pathogenic mechanisms, both direct and compensatory, contribute to disease phenotypes, and provide a salient example of how a pathogenic ion channel variant can cause opposite functional effects in closely related neuron subtypes due to interactions with other ionic conductances.