Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT

User menu

Search

  • Advanced search
eNeuro
eNeuro

Advanced Search

 

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT
PreviousNext
Research ArticleNew Research, Neuronal Excitability

Dramatically Amplified Thoracic Sympathetic Postganglionic Excitability and Integrative Capacity Revealed with Whole-Cell Patch-Clamp Recordings

Michael Lee McKinnon, Kun Tian, Yaqing Li, Alan Joel Sokoloff, Meredith Lucy Galvin, Mi Hyun Choi, Astrid Prinz and Shawn Hochman
eNeuro 30 April 2019, 6 (2) ENEURO.0433-18.2019; https://doi.org/10.1523/ENEURO.0433-18.2019
Michael Lee McKinnon
1Department of Physiology, Emory University, Atlanta, GA 30322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Michael Lee McKinnon
Kun Tian
2Department of Biology, Emory University, Atlanta, GA 30322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yaqing Li
1Department of Physiology, Emory University, Atlanta, GA 30322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alan Joel Sokoloff
1Department of Physiology, Emory University, Atlanta, GA 30322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Meredith Lucy Galvin
1Department of Physiology, Emory University, Atlanta, GA 30322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mi Hyun Choi
1Department of Physiology, Emory University, Atlanta, GA 30322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Mi Hyun Choi
Astrid Prinz
2Department of Biology, Emory University, Atlanta, GA 30322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Astrid Prinz
Shawn Hochman
1Department of Physiology, Emory University, Atlanta, GA 30322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Shawn Hochman
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

This article has a correction. Please see:

  • Erratum: McKinnon et al., Dramatically Amplified Thoracic Sympathetic Postganglionic Excitability and Integrative Capacity Revealed with Whole-Cell Patch-Clamp Recordings - July 03, 2019

Abstract

Thoracic paravertebral sympathetic postganglionic neurons (tSPNs) comprise the final integrative output of the distributed sympathetic nervous system controlling vascular and thermoregulatory systems. Considered a non-integrating relay, what little is known of tSPN intrinsic excitability has been determined by sharp microelectrodes with presumed impalement injury. We thus undertook the first electrophysiological characterization of tSPN cellular properties using whole-cell recordings and coupled results with a conductance-based model to explore the principles governing their excitability in adult mice of both sexes. Recorded membrane resistance and time constant values were an order of magnitude greater than values previously obtained, leading to a demonstrable capacity for synaptic integration in driving recruitment. Variation in membrane resistivity was the primary determinant controlling cell excitability with vastly lower currents required for tSPN recruitment. Unlike previous microelectrode recordings in mouse which observed inability to sustain firing, all tSPNs were capable of repetitive firing. Computational modeling demonstrated that observed differences are explained by introduction of a microelectrode impalement injury conductance. Overall, tSPNs largely linearly encoded injected current magnitudes over a broad frequency range. Thus, whole-cell recordings reveal tSPNs have more dramatically amplified excitability than previously thought, with greater intrinsic capacity for synaptic integration and with the ability for maintained firing to support sustained actions on vasomotor tone and thermoregulatory function. Rather than acting as a relay, these studies support a more responsive role and possible intrinsic capacity for tSPNs to drive sympathetic autonomic function.

  • computational model
  • firing properties
  • membrane properties
  • mouse
  • paravertebral ganglia

Significance Statement

Thoracic sympathetic postganglionic neurons (tSPNs) represent the final neural output for control of vasomotor and thermoregulatory function. We used whole-cell recordings and computational modeling to provide broad insight on intrinsic cellular mechanisms controlling excitability and capacity for synaptic integration. Compared to past intracellular recordings using microelectrode impalement, we observed dramatically higher membrane resistivity with primacy in controlling enhanced tSPN excitability and recruitment via synaptic integration. Compared to reported phasic firing, all tSPNs fire repetitively and linearly encode injected current magnitude to firing frequency over a broad range. Modeling studies suggest microelectrode impalement injury accounts for differences in tSPN properties previously observed. Overall, intrinsic tSPN excitability plays a much greater role in the integration and maintenance of sympathetic output than previously thought.

Introduction

Sympathetic postganglionic neurons (SPNs) within paravertebral chain ganglia represent a large fraction of the final output of the sympathetic nervous system. Whereas prevertebral sympathetic ganglia are typically associated with one or more visceral organs in a discrete location (celiac ganglion, superior/inferior mesenteric ganglion), thoracic paravertebral chain ganglia are associated with control of dispersed tissue systems such as vasculature, brown adipose tissue, sweat glands, and piloerector muscles (Jänig, 2006; Bartness et al., 2010). As such, the sympathetic chain can be thought of as a distribution system for sympathetic activity that spans the body. The vast majority of paravertebral postganglionic neurons in mice are adrenergic (Gibbins, 1991; Jobling and Gibbins, 1999) since sweat glands, innervated by cholinergic postganglionic neurons, are largely absent in the mouse (Lu and Fuchs, 2014).

Traditionally, thoracic SPNs (tSPNs) have been envisioned as passive followers of intraspinal preganglionic neuronal activity. By this viewpoint, postganglionic neurons fire if and only if preganglionics fire and serve as 1:1 relays that pass central commands to the periphery (Jänig, 2006). This relationship is explained by the “n + 1” rule, wherein postganglionic neurons receive n small synaptic inputs, and one major, always suprathreshold input which leads to firing with a high safety factor. The n smaller synaptic inputs are typically sub-threshold and infrequent, and are not thought to contribute appreciably to the firing rate (McLachlan et al., 1998; Karila and Horn, 2000; McLachlan, 2003; Wheeler et al., 2004; Rimmer and Horn, 2010). However, recent evidence from rodent sympathetic ganglia has shown that postganglionic neurons play a more active role in shaping sympathetic output (Bratton et al., 2010; Springer et al., 2015). In light of these findings, we must reconsider the role that SPNs play in synaptic integration and signal transmission.

Despite their critical importance as the final output controlling sympathetic neural commands, surprisingly little is known about the SPNs in thoracic segments (tSPNs) of the sympathetic chain. The most likely reason for this is their near inaccessibility to in vivo study, and the relative difficulty for in vitro cellular characterization. Because of this difficulty, electrophysiological properties of sympathetic neurons have been largely inferred from recordings in other mammalian paravertebral sympathetic ganglia, namely the superior cervical ganglion (SCG; Eccles, 1935; Erulkar and Woodward, 1968; Purves and Wigston, 1983; Li and Horn, 2006) and to a lesser extent the stellate and lumbar ganglia (Jänig, 1985; Cassell et al., 1986; Valli et al., 1989; Bratton et al., 2010). Compared to SCG, mouse tSPNs are smaller, have less elaborate dendritic arbors, are likely more excitable, and differ in measures of action potential (AP) shape (Jobling and Gibbins, 1999). Thoracic ganglia also contain a different subset of molecularly distinct SPN subpopulations and project to different end-organs (Jänig, 2006; Furlan et al., 2016). Unfortunately, few studies have directly characterized electrophysiological properties of thoracic ganglia (Blackman and Purves, 1969; Lichtman et al., 1980; Jobling and Gibbins, 1999). These studies used sharp microelectrodes for recordings, which likely introduce a considerable impalement injury conductance compared to whole-cell patch-clamp recordings (Staley et al., 1992; Springer et al., 2015). This injury-induced conductance alters basic membrane properties, such as input resistance and membrane time constant, which reduce recruitment and synaptic integrative actions according to classical cable theory (Rall, 2011; Springer et al., 2015). The impalement conductance introduced by microelectrode recordings can also prevent expression of repetitive firing properties (Cymbalyuk et al., 2002; Springer et al., 2015). Indeed, while it is generally thought that most paravertebral SPNs fire phasically (Cassell et al., 1986; Jobling and Gibbins, 1999; Li and Horn, 2006), whole-cell recordings in SCG support repetitive rather than phasic firing (Springer et al., 2015). Whether repetitive firing properties are predominant in tSPNs remains unknown.

The aim of the present study is to investigate the electrophysiological properties of tSPNs using whole cell recordings to more accurately characterize the cellular mechanisms that drive and modulate excitability of tSPNs. We furthermore matched recordings to a computational model to better understand how synaptic inputs and passive and active membrane properties interact to recruit neurons and generate the firing properties observed.

Materials and Methods

Animals

All animal procedures were performed in accordance with the Emory University Institutional Animal Care and Use Committee’s regulations and conformed to the Guide for the Care and Use of Laboratory Animals. Experiments were performed on adult (P37–P379) C57BL/6 mice (RRID:IMSR_JAX:000664). Mice were anesthetized with inhaled isoflurane and maintained or killed with urethane (intraperitoneal injection, 40 mg/kg for transcardial perfusions, ∼500 mg/kg for in vitro electrophysiology). Complete sedation or death was confirmed by lack of foot pinch and eye blink reflex.

Immunohistochemistry

Neurotransmitter identity

Two ChAT-eGFP mice (RRID:IMSR_JAX:007902), a male and a female (P91 and P101, respectively) were anesthetized and transcardially perfused with heparinized saline (0.9% NaCl, 0.1% NaNO2, 10-units/ml heparin), followed by 4% paraformaldehyde (0.5 M phosphate, 4% paraformaldehyde, NaOH). Tissue was post-fixed overnight, then transferred to a 15% sucrose solution and stored at 4°C. Sympathetic chains were isolated from stellate (T1 and T2) to T12/13. Tissue was embedded (TissueTek optimal cutting temperature compound), sectioned on a cryostat (−21°C, 8-μm slice thickness), and mounted on glass slides. Tissue was washed in 0.1 M PBS for 1 h and permeabilized with PBS containing 0.3% Triton X-100 (PBS-T) overnight. Sections were subsequently incubated for 2–3 d with primary antibodies: sheep anti-tyrosine hydroxylase (TH; Millipore, 1:100, RRID:AB_90755) and chicken anti-green fluorescent protein (Jackson, 1:100). Preparations were then washed in PBS-T (3 × 30 min) and incubated for 1.5 h with secondary antibodies: Cy3 donkey anti-sheep (Abcam, 1:250) and Alexa Fluor 488 donkey anti-chicken (Abcam, 1:250). Slides were washed a final time in PBS-T (20 min), then 50 mM Tris-HCl (2 × 20 min) and allowed to dry before being coverslipped [SlowFade Gold antifade reagent with 4',6-diamidino-2-phenylindole (DAPI)]. Sections were visualized under a fluorescent microscope (Olympus BX51). Cells with visible nuclei were counted and assessed for neurotransmitter identity. Interanimal cell count variability was substantial (6494 vs 19,721 cells).

Cell diameter

Six C57Bl/6J mice (RRID:IMSR_JAX:000664), 5 males and one female (all ∼P60) were transcardially perfused, as above. T5 Sympathetic ganglia were isolated. Unmounted tissue was washed in PBS-T overnight. Slides were subsequently incubated for 5 d with sheep anti-TH (Millipore, 1:100, RRID:AB_90755). Preparations were then washed in PBS-T (3 × 2 h) and incubated for 3 d with Alexa Fluor 488 donkey anti-sheep (Jackson, 1:100). Slides were washed a final time in PBS-T (2 h), then 50 mM Tris-HCl (2 × 1 h). Intact ganglia were mounted on glass slides and coverslipped (SlowFade Gold antifade reagent with DAPI). TH-immunoreactive cells were visualized under a fluorescent microscope (Olympus BX51, 40× objective) using a Microfire digital camera (Optronics), and traced using Neurolucida software (MBF Bioscience, RRID:SCR_001775). Cell diameters were calculated as the arithmetic mean of minimum and maximum Feret. Diameter was only determined for cells with a discernible perimeter (176 ± 131 cells per ganglion) representing a mean 71% of the total TH+ cell population (range of 36–95% neurons/ganglia measured). As diameter distributions were comparable between ganglia, the possibility of sampling bias in estimated cell diameter is unlikely. Results are reported as mean ± SD.

Electrophysiology

Tissue preparation

Mice were killed and the spinal column was quickly dissected out with sympathetic chain and spinal roots attached. Figure 1A provides a simplified schematic of the anatomic organization of intraspinal preganglionic and paravertebral postganglionic neurons. The remaining tissue was incubated in continually oxygenated ACSF containing collagenase (20-mg Type III per 1-ml ACSF, Worthington Biochemical Corporation) for 1.5 h. ACSF used for incubation was buffered with either bicarbonate or HEPES. No difference was observed as a result of incubation buffer. Following incubation, tissue was vortexed to remove adherent fat and washed with ACSF several times to eliminate residual collagenase. The intact sympathetic chain was removed by severing rami, and was then pinned down into a clear Sylgaard recording dish (Fig. 1B), through which recirculating, oxygenated ACSF was continually perfused.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Cell size and composition. A, Simplified schematic depicting the anatomic organization of preganglionic and postganglionic neurons. B, Recording setup. Sympathetic chains are pinned down in a silicone chamber, superfused with oxygenated ACSF. They are then visualized under a microscope and recorded using a glass patch-clamp electrode. C, Confocal slice through whole-mounted tissue showing TH immunolabeling and nuclear labeling with DAPI. Note the numerous smaller and more intensely labeled nuclei that are presumably non-neuronal cells. Scale bar represents 50 µm. D, Histogram showing distribution of TH-IR cell diameters in T5 ganglia of six animals.

Whole-cell recordings

Whole-cell patch recordings were obtained from postganglionic cells at room temperature. Cells were identified using an upright microscope (Olympus, BX51WI) affixed with a low-light camera (Olympus, OLY-150). Patch electrodes were pulled on a vertical puller (Narishige, PP-83) from 1.5-mm outer diameter filamented, borosilicate glass capillaries (World Precision Instruments, stock #TW150F-4) for a target resistance of 5–9 MOhm. Signals were amplified using a MultiClamp 700A and digitized at 10 kHz using a Digidata 1322A and Clampex software (Molecular Devices, RRID:SCR_011323).

We considered for analysis all cells which displayed clearly defined APs on depolarization by square current steps. Of these, cells were excluded if more than 100 pA was required to hyperpolarize a cell to −70 mV (indicative of a significant leak), if APs appeared stunted (indicative of an incomplete breakthrough), or if membrane potential was highly variable (indicative of improper seal formation). All cells which met these criteria (n = 35) had resting membrane potentials more negative than −50 mV and input resistances higher than 200 MΩ. All recordings were made in ACSF containing: 127.99 mM NaCl, 1.90 mM KCl, 1.30 mM MgSO4·7H2O, 2.40 mM CaCl2·2H2O, 1.20 mM KH2PO4, 9.99 mM glucose, and 26.04 mM NaHCO3. ACSF pH was adjusted to 7.4 after saturation with gas (95%O2, 5%CO2) at room temperature. Intracellular patch-clamp solution contained: 140.0 mM K-gluconate, 11.0 mM EGTA, 10 mM HEPES, and 1.32 mM CaCl2; pH was adjusted to 7.3 using KOH. Target osmolarity was <290 mOsm. In most recordings (25/39 cells), support solution was added consisting of 4.0 mM ATP and 1.0 mM GTP.

The ratio of male to female mice was ∼1:1. Recordings were taken from the right thoracic ganglia, with the majority of recordings coming from T5. The number of cells from ganglia T3 through T12 was 2, 6, 17, 3, 2, 1, 1, 0, 1, and 2, respectively. We initially assessed sex and segment related differences in cell properties. Our sample size precludes a more detailed analysis of segment-specific differences, but we were able to pool data from rostral (T3 and T4) and caudal (T5 to T12) segments (Furlan et al., 2016). Statistical tests revealed no differences with respect to sex or rostro-caudal location (t test, two-tailed, unequal variance, Šidák corrected α = 0.0019), so all data were pooled for additional analyses.

Data analysis

All cellular properties were analyzed in Clampfit (Molecular Devices, RRID:SCR_011323) or MATLAB (MathWorks, RRID:SCR_001622). All parameters were estimated from a single set of current steps for each cell. This ensured that parameters for a given cell were estimated at nearly the same point in time. In current clamp mode, membrane voltage response to hyperpolarizing current steps of at least 1.5 s was fit to an exponential of the form Equation 1 using the Levenberg–Marquardt algorithm built in to Clampfit. The value of membrane time constant (τm) was calculated in this manner (Golowasch et al., 2009). Rin was estimated by dividing maximal voltage deflection (ΔV) by the injected current (Iinj; Eq. 2). Membrane capacitance (Cm), a measure of total cell surface area, was estimated by dividing τm by Rin (Eq. 3).

Embedded Image (1)Embedded Image (2)Embedded Image (3)

Measured rheobase current (Irheo) was taken as the smallest long-duration (1.5 s or longer) positive current injection which elicited a single spike. In the case that an incremental increase in current elicited multiple spikes, rheobase was estimated to be the mean of the adjacent subthreshold and suprathreshold steps, e.g., if 30 pA did not elicit any spikes but 40 pA elicited several, the measured rheobase estimate would be 35 pA. to achieve a more finely-grained estimate of rheobase, we also calculated rheobase based on the equation:Embedded Image where Vth is the AP threshold, taken to be the point at which the first derivative of voltage, dVm/dt, begins to increase (Platkiewicz and Brette, 2010). Measured values related to AP and post-spike afterhyperpolarization (AHP) characteristics were taken from traces elicited at minimal suprathreshold current, i.e. the smallest current magnitude used which elicited at least one AP.

The parameters of the fast AHP (fAHP) varied as a function of firing rate, so analysis of fAHP properties was limited to cells which fired a single spike at minimal suprathreshold current intensity. AP amplitude was defined as the difference between the peak voltage and threshold. AP half-width is the width of the spike at half AP amplitude. fAHP amplitude was defined as the difference between peak negative voltage and steady-state voltage at rheobase current injection. fAHP half-decay is the time it takes for the fAHP to decay to half its amplitude. fAHP duration is the time between spike onset and return to baseline (Hochman and McCrea, 1994). Slow AHP (sAHP) amplitude was defined as the difference between peak negative voltage and baseline (holding voltage). sAHP half-decay is the time it takes for the sAHP to decay to half its amplitude.

Instantaneous firing rate (IFR) was taken as the inverse of the interspike interval. Maximal firing rate was the IFR for the first spike pair at the beginning of current onset. Sustained firing rate was the mean IFR for the last three interspike intervals, given they occur during the last half of the depolarizing current step. Frequency-current (ƒ-I) slope is the slope of the linear regression of the ƒ-I curve. The spike rate adaptation (SRA) ratio is defined as the ratio between the maximal and sustained firing rate at a given current injection (Venance and Glowinski, 2003; Miles et al., 2005). To directly compare firing rate across cells with variable Rin and ƒ-I curves, we used the sustained firing rate at twice minimal suprathreshold current injection. Current step duration was at least 1.5 s for all cells, and 3 s for the majority. Liquid junction potential was calculated to be −9.8 mV and empirically measured to be −13 mV. All values of absolute voltage (resting membrane potential, absolute threshold, peak voltage) were adjusted by −10 mV to approximately account for liquid junction potential. For example, a recorded AP peak of 30 mV would be reported as 20 mV and a recorded RMP of −60 mV would be reported as −70 mV.

Computational modeling

Single neuron model

We built a conductance-based neuron model to help understand observed results in relation to their underlying biophysical mechanisms. While tSPNs do possess dendrites, their dendritic arborizations are relatively simple. We therefore assume that ganglionic cells are electrotonically compact, and that a single-compartment model can replicate all essential physiologic properties observed in experiments. All currents included in the model have been observed in rodent sympathetic ganglia (Galvan and Sedlmeir, 1984; Sacchi et al., 1995; Jobling and Gibbins, 1999; Rittenhouse and Zigmond, 1999) and transcript expression in mouse thoracic ganglia has recently been confirmed by a single-cell RNA sequencing study (Furlan et al., 2016).

The model is based on a model of bullfrog paravertebral sympathetic ganglia (Wheeler et al., 2004), which represents the most complete available computational model of a paravertebral neuron. From this model the following conductances were taken: a fast sodium current, INa; a delayed-rectifier potassium current, IKd; a slow and non-inactivating potassium current, IM; and a voltage-independent leak current, Ileak. Additional conductances were added from models derived in other species. These include the following: a fast transient potassium current, IA (Rush and Rinzel, 1995); a hyperpolarization-activated inward current, Ih (Kullmann et al., 2016); and a calcium-dependent potassium current, IKCa (Ermentrout and Terman, 2010). IKCa depends on intracellular calcium concentration, [Ca2+], so a model of persistent calcium current, ICaL (Bhalla and Bower, 1993) and somatic calcium dynamics (Kurian et al., 2011) were added as well. Model parameters were then tuned to fit recorded data from the present study.

The membrane voltage, V, is updated according to the equation:Embedded Image (4)

Membrane capacitance, Cm, was set at 100 pF to approximate the mean in recorded neurons. Each current, Ii, is described by the equation:Embedded Image (5)where Gi is the maximal conductance, Ei is the reversal potential, and m and h are gating variables for activation and inactivation. A standard model neuron was used to replicate the majority of observed phenomena. Maximal conductances of this standard neuron are indicated in Table 1. The standard model was modified as necessary to fit individual recordings, which comprise a heterogeneous population. The reversal potentials for the various membrane currents are indicated in Table 1.

View this table:
  • View inline
  • View popup
Table 1.

Model parameters

The activation and inactivation variables m and h are updated by the equation:Embedded Image (6)

The intracellular calcium concentration is updated by:Embedded Image (7)where λ = 0.01 is the ratio of free to bound [Ca2+], α = 0.002 μM·ms−1·pA−1 is the conversion factor from current to concentration, and kCaS = 0.024 ms−1 is the somatic [Ca2+] removal rate.

Impalement simulation

To replicate impalement injury, an additional leak conductance was added to the model to simulate microelectrode impalement. This conductance, gimp, was modeled as a non-selective ohmic leak channel with Eimp = −15 mV. The impalement reversal potential was estimated by solving the Goldman–Hodgkin–Katz equation (Eq. 8), with equal permeabilities of the three major ionic species. This estimate agrees well with estimates of impalement reversal potential in bullfrog ganglia (Brown, 1988). For analysis, model neurons were subjected to a bias current and held at −70 mV, unless otherwise stated. gimp was normally set at 0 nS, and was only included where indicated for simulation of impalement.Embedded Image (8)

Synapse simulation

Synaptic input was implemented with equation:Embedded Image (9)where Isyn is synaptic current, A is conductance amplitude, and Esyn is the synaptic reversal potential set at 0 mV. Synaptic conductance, gsyn, was calculated from the equation:Embedded Image (10)where τr and τd are the rise and decay time constants, respectively, and s is a scaling factor to normalize the amplitude to 1 nS. Equations were adapted from Springer et al. (2015). Rise and decay time constants were 1 and 15 ms, respectively, as estimated from voltage clamp recordings of spontaneous synaptic activity.

Code accessibility

Source code for all simulation and analysis are available online at https://github.com/pinewave/tSPN and ModelDB (Hines et al., 2004, accession #245926). Simulation and analysis scripts were written in Python 2.7.10 and executed in PyCharm (CE 2017.1.2) on macOS 10.12.3 with a 1.7-GHz processor. Scripts were also translated into MATLAB code and executed on Windows 10 with a 2.4-GHz processor. All differential equations were integrated using an Exponential Euler method with a time step of 0.1 ms (Prinz et al., 2004).

Experimental design and statistical analysis

The present study used a descriptive design. Statistical analyses were performed in Microsoft Excel. Basic properties are presented as mean ± SD in Table 2. Correlations were determined by Pearson’s correlation coefficient, r. A two-tailed t test was used to calculate each p value. To control for 30 multiple comparisons and maintain an experiment-wise α = 0.05, a Šidák corrected α = 0.0017 was used to assign statistical significance. In some cases, parameter pairs with moderate values of r, (|r|>0.4) failed to reach significance as a result of intrinsic variability inherent within this data. Such correlations are reported as moderate, and should be interpreted cautiously. Exact r, R 2, and p values are presented in Table 3.

View this table:
  • View inline
  • View popup
Table 2.

Basic properties of tSPNs

View this table:
  • View inline
  • View popup
Table 3.

Selected correlations between tSPN parameters

Results

Thoracic ganglia composition

Postganglionic neurons have been shown to be either adrenergic or cholinergic (Jobling and Gibbins, 1999; Jänig, 2006). to assess neurotransmitter identity throughout the sympathetic chain, we used a choline-acetyltransferase (ChAT) transgenic mouse which fluorescently labels putative cholinergic postganglionic neurons (ChAT::eGFP), and co-immunolabelled tissue with an antibody to TH to label putative adrenergic neurons. Neurons were counted from stellate (T1 and T2) to T13 ganglia. We found that TH+ neurons comprised >97% of the population and no ganglion contained >6% presumptive cholinergic neurons. This agrees with prior findings in rodent thoracic ganglia which found few cholinergic neurons (Schäfer et al., 1998; Jobling and Gibbins, 1999; Masliukov and Timmermans, 2004; Schütz et al., 2015; see also Furlan et al., 2016). This indicates that cholinergic neurons in thoracic ganglia are rare. We therefore assume that a large majority of recorded cells were adrenergic.

As the majority of electrophysiological recordings focused on T5, diameter and number of TH-IR cells were examined in T5 ganglia in a separate sample from six adult mice (Fig. 1C). The mean number of TH-IR neurons counted in T5 ganglia was 247 ± 127 (ranging from 106 to 418). tSPNs had a mean cell diameter of 23.8 ± 5.4 μm with cell size distribution shown in Figure 1D. These values are smaller than those reported previously (31.0 ± 1.2 μm; cf. Jobling and Gibbins, 1999). Differing methodology may explain this discrepancy. The aforementioned study measured the diameter of dye-filled cells after microelectrode impalement, which would preferentially target larger neurons sustained f-l slope but not maximal f-l slope (Brown, 1981).

Passive membrane properties

Whole-cell patch-clamp recordings were acquired from 35 tSPNs obtained from 30 adult mice. Basic cellular properties are summarized in Table 2. The distribution of resting membrane potential is shown in Figure 2A. Input resistance (Rin) and membrane time constant (τm) were, on average, an order of magnitude higher than values recorded using microelectrode recordings in mouse (Jobling and Gibbins, 1999; Fig. 2B) and guinea pig (Blackman and Purves, 1969) thoracic ganglia. Rin was strongly correlated with τm (Fig. 2B), but not cell capacitance (Cm), an estimate of cell size. This indicates that membrane resistivity, but not cell size, is primarily responsible for the variability seen in resistance measures (Gustafsson and Pinter, 1984). Cm was also moderately correlated with τm. A summary of correlation parameters is provided in Table 3.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Passive membrane properties. A, Histogram showing distribution of resting membrane potential values. B, Input resistance is highly correlated with membrane time constant. Solid line indicates linear least-squares fit. Filled black circle represents population mean. Red filled circle represents population mean from (Jobling and Gibbins, 1999). C, Example of synaptic summation leading to AP recruitment in a particularly active recording. Shown is a raster of epochs of spontaneous synaptic activity. Cell resting membrane potential was −60 mV. In this neuron, a τm of 109 ms led to comparably long EPSP membrane voltage decay τs. Vertical scale bar is 20 mV; horizontal scale bar is 500 ms. D, Top, Model neuron subjected to simulated synaptic input fires in response to synaptic summation. Middle, If an impalement conductance is added, synaptic summation is no longer effective. Bottom, Simulated gsyn used to generate voltage traces. Horizontal scale bar is 500 ms; vertical scale bars are 20 mV and 5 nS, respectively.

One impact of larger τm is longer duration spontaneous EPSPs (sEPSPs) and consequently greater capacity for temporal summation. Spontaneous synaptic activity is often observed in whole-cell recordings, including instances of sEPSP summation that lead to recruitment of APs (Fig. 2C). In this neuron, a τm of 109 ms led to comparably long sEPSP membrane voltage decay τs.

To explore the impact of preserved passive membrane properties on synaptic summation, we implemented a synaptic conductance in the computational model. A template conductance was constructed with Poisson-distributed events whose amplitudes and mean frequency match values from whole-cell voltage clamp recordings. This template conductance was used to stimulate a standard model neuron (Fig. 2D, top trace) and a model neuron with simulated microelectrode impalement injury (middle trace). In the intact cell, synaptic events are larger in amplitude and synaptic summation can lead to cell recruitment. In the model neuron with simulated impalement, AP recruitment was observed only in response to the largest single synaptic event.

Rheobase

The current required to depolarize a cell from its holding potential to firing threshold (rheobase) was examined in 35 cells by injecting long duration (1.5–3 s) pulses through patch electrodes. to control for the possible influence of a variable resting membrane potential on rheobase, tonic bias current was injected to hold cells at approximately −70 mV before rheobase estimation. Fluctuations in membrane voltage made it difficult to precisely set holding potential before injected current steps, and values ranged from −56 to −83 mV. We compared actual holding voltage against rheobase to determine whether this variability altered rheobase estimation. Rheobase was not correlated with holding potential.

Voltage threshold was assessed at minimal suprathreshold current intensity. Assuming cell depolarization is governed by Ohmic or non-rectifying processes, the ratio of relative voltage threshold to input resistance would predict rheobase (Gustafsson and Pinter, 1984). Indeed, measured and calculated rheobase are well correlated (Fig. 3A; Table 3) and approximately equal, indicating that rectifying currents do not play a major role in determining rheobase for the population. However, deviation of calculated rheobase values above and below those predicted by ohmic processes support a role for voltage-dependent conductances (Gustafsson and Pinter 1984). Values of measured and calculated rheobase are presented in Table 2. As calculated rheobase provides a more precise index of excitability, further analysis focuses on this parameter.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Factors affecting rheobase. A, There is a strong correlation between measured rheobase and calculated rheobase and the two values are approximately equal. This suggests that rheobase is predominantly governed by ohmic phenomena. Dashed line is line of identity. B, Calculated rheobase was well correlated with input conductance in recorded neurons, open circles. Gray line represents the rheobase versus input conductance relationship for a single model neuron chosen to fit experimental data. Standard model cell with GM = 20 nS and GA = 15 nS. A, B, Solid line represents least squares regression. C, Strength-duration curves for model neurons. Black, standard model cell with no impalement conductance. Black arrow indicates rheobase. Gray, standard cell with gimp = 7 nS, comparable to a microelectrode recording. Gray arrow indicates much higher rheobase for impaled cell.

Calculated rheobase current was strongly correlated with input conductance, gin = Rin −1 (Fig. 3B), and moderately correlated with the inverse of time constant, τm −1, but was uncorrelated with capacitance. A summary of correlation parameters is provided in Table 3.

We further investigated the relationship between input conductance and rheobase in a model cell. We adjusted gleak to vary input conductance of a model neuron over most of the range observed in recorded neurons (0.5–3 nS). Bias current was adjusted to hold the model cell at −70 mV. The rheobase was then calculated for each value of input conductance by using a binary search algorithm to find the minimal injected current which produces a spike (Fig. 3B, gray line). The results show that there is a deterministic relationship between rheobase and input conductance that can help to explain some of the correlation observed in recorded neurons. However, given the variability of rheobase measures in recorded cells with comparable values for input conductance, it is clear that input conductance alone does not fully explain the range of rheobase values observed in recorded neurons.

Rheobase values were 80–90% lower than values estimated in tSPNs previously with microelectrode recordings in both mouse and guinea pig (Blackman and Purves, 1969; Jobling and Gibbins, 1999). Reduced rheobase values indicate that tSPNs are much more excitable than previously considered. To more fully explore the influence of microelectrode impalement on cell excitability, we constructed strength-duration curves for model cells. In a standard model cell, the strength-duration curve follows a characteristic inverse curve. After implementation of an impalement conductance consistent with a microelectrode recording, rheobase is increased ∼8-fold (Fig. 3C). This is consistent with the discrepancy between our experimental findings and the aforementioned studies using microelectrodes.

Repetitive firing

Increasing current steps were delivered to assess repetitive firing properties from a holding potential of approximately −70 mV. All cells (n = 35) were capable of repetitive firing in response to sustained current injection. This contradicts an earlier report that tSPNs fire phasically in response to depolarization (Jobling and Gibbins, 1999). Figure 4A shows an example of a recorded cell which fires repetitively at progressively higher frequency in response to increasing depolarizing current steps (top). A model neuron that used known voltage-dependent conductances for paravertebral sympathetic neurons and incorporated values for input conductance obtained from our whole-cell recordings was able to replicate repetitive firing (Fig. 4A, bottom).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

tSPNs exhibit repetitive firing. A, Top, Representative trace from a tSPN showing increases in repetitive firing frequency in response to increasing current steps. Bottom, Model neuron also showing repetitive firing. Standard model with GM = 30, GKCa = 70, GA = 80, Gleak = 2 nS. Injected current from left to right in both recorded neuron and model is 30, 50, 70, 90, 110, 130 pA. Scale bar is 1 s. B, ƒ-I relations for recorded and model neurons. i, Maximal IFR is plotted versus injected current for all cells. ii, Same as Ci with sustained firing rate. Red line in i and ii is maximal and sustained ƒ-I curve from model neuron in A. iii, Maximal ƒ-I curve from a model neuron in which gleak was adjusted from 0 to 3 nS. Note: as varying input conductance also changes holding potential, each model neuron was subjected to a different holding current to hold the initial voltage at –70 mV. Also note that gleak is distinct from gimp. iv, Corresponding sustained ƒ-I curves. Model parameters other than gleak are the same as in A. C, Maximal ƒ-I slope is positively correlated with input resistance. D, Maximal ƒ-I slope is negatively correlated with calculated rheobase.

ƒ-I relations were obtained by plotting the maximal (initial) and sustained firing rate versus injected current magnitude. Figure 4Bi,ii shows the maximal and sustained, respectively, ƒ-I curves for all cells. Maximal IFR did not exceed 28 Hz, while sustained firing rate did not exceed 17 Hz for the highest steps given. ƒ-I curves were approximately linear. to determine the role that input conductance plays in determining ƒ-I relations in tSPNs, we selected a model neuron that matches the experimental ƒ-I curves and then systematically changed input conductance by varying gleak from 0 to 3 nS. Of note, varying input conductance also changes holding potential so to remain consistent with experimental protocol, each model neuron was subjected to a different holding current to hold the initial voltage at −70 mV. Figure 4Biii,iv demonstrates that altering gleak can shift the ƒ-I curve, but it does not appear to change the ƒ-I slope. Thus, input conductance cannot fully account for the range of ƒ-I curves observed. To determine whether other model parameters are capable of changing ƒ-I slope, we systematically varied each parameter and observed its influence on maximal and sustained ƒ-I curves. Most notably, Cm appears to influence the slope of the maximal ƒ-I curve while GCaL and GKCa impact the slope of the sustained ƒ-I curve. Other model parameters (GNa, GK, GM, GA, Gleak) are able to shift ƒ-I curves without significantly altering slope. Thus, we are able to match any realistic ƒ-I curve by adjusting model parameters, which implies that a host of intrinsic cellular properties are responsible for the range of ƒ-I curves we observed.

Slope for both maximal and sustained ƒ-I curves was calculated as a measure of excitability (Zimmerman and Hochman, 2010). In short, a cell with a higher ƒ-I slope would respond to an incremental change in current with a higher change in firing frequency. In this way, ƒ-I slope can be thought of as the gain between input and output of a neuron. Values for maximal and sustained ƒ-I slope are given in Table 2.

We assessed the role of variations in Rin and rheobase in cell excitability based on ƒ-I slope measures. Maximal firing rate at 100-pA current injection was significantly correlated with Rin. Maximal ƒ-I slope was moderately correlated with Rin (Fig. 4C), and moderately and negatively correlated with calculated rheobase (Fig. 4D). No such relationship was found for τm or Cm. A summary of correlation parameters is provided in Table 3. Cells with lower rheobase and higher Rin had higher ƒ-I slopes, suggesting that Ohmic properties contribute to the ƒ-I response.

Impalement simulation

The discrepancy between observations of phasic and repetitive firing likely arises as a result of leak introduced by microelectrode impalement (Springer et al., 2015). We undertook additional modeling to test whether an impalement injury can convert repetitive to phasic firing. An additional impalement conductance, gimp, was added to a standard model cell. Reversal potential of gimp was set at −15 mV. We explored the relationship between gimp and firing type over a range of conductance and current injection combinations (Fig. 5). For a given set of gimp and injected current, a model cell was characterized as non-firing (N), phasic firing (P), or repetitively firing (R). Setting gimp to 7 nS results in an input resistance of ∼100 MΩ, the mean value of input resistance reported by Jobling and Gibbins (1999). With this level of microelectrode leak, non-firing was observed in response to subthreshold current injection (Fig. 5Ai), and phasic firing was observed in response to suprathreshold current injection over the range of values tested by Jobling and Gibbins (Fig. 5Aii). However, when gimp was set at 0 nS, analogous to a whole-cell recording, repetitive firing was observed instead (Fig. 5Aiii).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Simulated impalement can alter firing properties. A, Impact of injected current and impalement conductance, gimp, on firing properties. Voltage response of a model cell to subthreshold (i) and suprathreshold (ii) current injection after gimp is set to 7 nS, analogous to a microelectrode recording. Only phasic firing is observed. iii, Repetitive firing is observed when gimp is set to 0 nS, analogous to a whole-cell recording. When gM is removed from the model, the same parameters used in i, ii lead to non-firing (iv) and repetitive firing (v). B, Shaded regions indicate the set of all parameters which lead to non-firing (N, red), phasic firing (P, blue), and repetitive firing (R, black/gray). At gimp = 0 nS, the model neuron transitions rapidly from N to R, and repetitive firing results from any current injection above ∼20 pA. At gimp = 7 nS, the model neuron transitions from N to P at around 200-pA current injection, and repetitive firing is not observed for injected current less than 500 pA. C, same as B with gM set to 0 nS. Removing IM from the model by setting gM = 0 nS eliminates phasic firing altogether, i.e., cells transition directly from N to R regardless of impalement conductance. Open circles in B, C indicate the gimp and injected current values used to generate traces in A. Standard model with gleak = 0.5 nS.

Prior studies have reported that phasic firing sympathetic neurons could instead fire repetitively if IM was blocked (Brown and Adams, 1980; Cassell et al., 1986). To test this, we blocked IM in our model cell by setting gM to 0 nS. This change completely eliminated phasic firing in the model, and only non-firing (Fig. 5Aiv) or repetitive firing was observed (Fig. 5Av).

To more extensively characterize this phenomenon, the boundaries between each of the three firing types were identified using a binary search algorithm. In the case where IM is included in the model (Fig. 5B), there was a rapid transition from repetitive to phasic firing as gimp is increased. When IM is removed (Fig. 5C), the phasic firing region (P) disappears.

SRA

Implicit in the observation that sustained firing rates were lower than initial observed frequencies is that all cells displayed SRA, or a decrease in firing rate over time. We were able to replicate SRA in our model (Fig. 6A). The time course of adaptation consists of a fast and a slow phase (Fig. 6B).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Modeling suggests that SRA in tSPNs depends on both IM and IKCa. Ai, Representative trace showing tSPN response to 50-pA current injection. Note that the interspike interval increases over time, corresponding to a decrease in instantaneous frequency. ii, Trace from a model cell chosen to fit the recording shows similar SRA for 50-pA current injection. Maximal conductances are: GNa = 400 nS, GK = 3000 nS, GCaL = 1.2 nS, GM = 40 nS, GKCa = 60 nS, GA = 80 nS, GH = 1 nS, Gleak = 2 nS. Scale bar in both panels is 1 s. B, Instantaneous frequency versus time for the same recorded cell at 50-, 70-, 90-, 110-, and 130-pA current injection (from bottom to top). The 50-pA curve (red) corresponds to the trace in Ai. Fast and slow components of adaptation are indicated. C, Maximal and sustained ƒ-I curves match well between recorded and modeled cell over a range of injected currents. Red, maximal (top, solid) and sustained (bottom, dashed) ƒ-I curves for the cell in Ai, B. Blue lines are the corresponding ƒ-I curves from the model cell in Aii. D, Instantaneous frequency versus time curves for the model cell in Aii. The recorded 50-pA curve from B is reproduced for comparison to the analogous curve generated in the model cell in Aii (blue). Black curves numbered 1–3 represent effect of removal of two conductances from the model. Removal of gKCa (curve 1) predominantly influences the slow SRA. Removal of gM (curve 2) predominantly influences the fast SRA. Removal of both (curve 3) eliminates SRA. The ordinate axis is shared among B–D.

The difference between the initial firing rate and the sustained firing rate becomes more pronounced as injected current is increased in all cells. This can be illustrated by comparing the maximal ƒ-I curve to the sustained ƒ-I curve in both recorded and model neuron over a range of current injection (Fig. 6C). This relationship between maximal and sustained firing rate is a common feature of adapting neurons (Benda and Herz, 2003).

Several mechanisms have been proposed to underlie SRA in different neuronal populations including Na+ channel inactivation (Miles et al., 2005), fAHP summation (Powers et al., 1999), activation of IKCa (Miles et al., 2005) and activation of IM (Yi et al., 2015). We selectively removed conductances from the model and determined which were primarily responsible for SRA (Fig. 6D). Removal of IKCa preferentially impaired the later phase of adaptation (curve 1) while removal of IM preferentially impaired the early phase of adaptation (curve 2). Removal of both IKCa and IM completely eliminated SRA (curve 3). Our model supports the conclusion that the combination of these conductances is necessary to replicate SRA.

AHP

AHP dynamics play an important role in regulating neuronal firing. Based on decay time, we identified three types of AHP within the thoracic ganglia. These include the fAHP after a single AP, and the sAHP and ultra-slow AHP (usAHP) after multiple APs (Fig. 7A).

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

AHP. A, Side by side comparison of three types of AHPs. i, fAHP present after single spike. ii, sAHP is present in the same cell only after repetitive firing. The half-decay time of fAHPs and sAHPs are indicated by the gray bar beneath each trace. Scale bar 1 s. iii, Depolarizing current steps (10–130 pA in 10-pA increments) in a different cell showing the progressive hyperpolarization characteristic of the usAHP. Scale bar 10 s. iv, Expanded view of voltage traces in Aiii indicated by vertical arrows. Note that the gray trace is hyperpolarized by 20 mV compared to the black trace and has a characteristic notch (arrowhead) on depolarization. Current injection profile is shown below each trace. B, fAHP half-decay was negatively correlated with maximal firing rate at twice the minimal suprathreshold current. C, SRA ratio is positively correlated with sAHP half-decay. Black line is the linear regression.

Fast post-spike after-hyperpolarization (fAHP) amplitude, half-decay time, and duration were measured at rheobase current injection (Fig. 7Ai). Parameters related to fAHP are summarized in Table 2. Half-decay time was very well correlated with duration and more reliably obtained, so further analysis focused on fAHP half-decay time. fAHP half-decay time was compared to passive membrane properties and rheobase. fAHP half-decay time was not correlated with Rin, Cm, or rheobase, but was moderately correlated with τm. Previous studies have reported an inverse relationship between fAHP duration and firing rate in motoneurons (Brownstone et al., 1992; Stauffer et al., 2007). To determine whether this relationship exists in postganglionic neurons, we plotted fAHP half-decay time versus sustained firing rate at two times the minimal suprathreshold current injection. Note that this corresponds to twice the current magnitude used to estimate fAHP half-decay time. We found that there is indeed a moderate negative correlation between fAHP half-decay time and maximal firing rate (Fig. 7B). fAHP half-decay time was also moderately correlated with sustained ƒ-I slope but not ƒ-I slope.

sAHPs were also observed following larger depolarizing steps that elicited higher repetitive firing frequencies (Fig. 7Aii). Only cells displaying obvious sAHP were analyzed (n = 27 of 35). sAHPs were measured at maximal current injection. Parameters related to sAHP are summarized in Table 2. sAHP half-decay time was four-fold longer on average than fAHP half-decay, but the two were not correlated. To examine the relationship between sAHP and SRA, we plotted the sAHP half-decay versus the SRA ratio for 27 cells (Fig. 7C). We found the two parameters were significantly correlated. A summary of correlation parameters for both fAHP and sAHP is provided in Table 3. As with SRA, our computational model showed that IM and IKCa were capable of reproducing sAHP after repetitive firing (not shown).

Prior work in the rabbit SCG identified a long-lasting AHP following sustained depolarization that was due to the ouabain sensitive Na+/K+-ATPase (Lees and Wallis, 1974). In the neonatal mouse spinal cord, it has been shown to be due to activation of α3 Na+/K+-ATPase (Picton et al., 2017).This AHP is unique in its ability to hyperpolarize a cell membrane beyond the reversal potential of K+. We identified an AHP with a similar time course. We injected depolarizing current to cause the cell to fire repetitively. This repetitive firing led to a steadily increasing hyperpolarization (Fig. 7Aiii). In the example shown, induced epochs of repetitive firing led to a 20-mV membrane hyperpolarization (Fig. 7Aiv). This feature was present in two of 14 cells tested with a current step protocol that would allow for its observation. Of note, the usAHP was observed only in relatively high resistance cells when ATP and GTP were included in the electrode solution. This AHP was also able to achieve a membrane potential of −101.7 ± 11.5 mV, which is more negative than the calculated −98 mV K+ reversal potential. The time course of this hyperpolarization is too long to be due to IM or IKCa.

Subthreshold conductances

Subthreshold conductances can play an important role in determining cell excitability and firing properties. We evaluated activation of these conductances with current steps that included assessment at hyperpolarized membrane potentials seen during the usAHP.

In response to depolarizing current steps, membrane voltage first followed an exponential time course with subsequent recruitment of voltage-gated conductances that altered the trajectory. In 24 of 35 cells, membrane trajectory exhibited a negative deflection from the exponential trajectory which preceded activation of voltage gated Na+ conductance. The observed deflection, or “notch,” in membrane led to a delay in the first AP in a train (Fig. 8Ai) and has been described previously in tSPNs (Jobling and Gibbins, 1999). This phenomenon was often observed at a holding potential of −70 mV, and became more pronounced with greater hyperpolarization (−90 mV). This is consistent with activation of the transient, voltage-gated A-type K+ current (IA). To test the contribution of IA to the notch and delayed firing, we held a model neuron at two different holding potentials and found that the change in trajectory was indeed attributable to de-inactivation of IA (Fig. 8Aii; Rush and Rinzel, 1995). Notably, a similar notch was observed in cells displaying usAHP (Fig. 7Aiv), demonstrating that the usAHP leads to a state of membrane hyperpolarization where IA would delay onset of firing.

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8.

Subthreshold conductances. Ai, A cell depolarized from −90 mV exhibits a characteristic notch (arrow) accompanied by a delay in spiking (black trace). The same cell depolarized from −70 mV does not have a notch (gray trace). ii, Model neuron showing comparable results with pre-spike inflection seen only for hyperpolarized trace. Standard model with GM = 10 nS, GKCa = 10 nS, GA = 90 nS, Gleak = 0 nS. Inset, Magnitude of hA at onset of current injection shows that IA is less inactivated (hA is higher) at hyperpolarized voltage, and IA takes longer to fully inactivate. Scale bars represent 500 ms for all panels. Bi, Voltage sag, indicated by arrow, on hyperpolarization beyond −90 mV in a cell held at −70 mV. Note that the effect becomes more pronounced with greater hyperpolarization. ii, Model neuron showing similar sag. Standard model with GA = 5 nS and Gleak = 0.1 nS. Ci, Hyperpolarizing trace from a different cell held at −50 mV showing rebound spiking associated with voltage sag. ii, Model neuron showing rebound spiking at the same holding voltage and current injection. Maximal conductances are: GNa = 200 nS, GK = 2000 nS, GCaL = 1.2 nS, GM = 20 nS, GKCa = 20 nS, GA = 20 nS, GH = 1 nS, Gleak = 2 nS. Removal of gH (curve 1) does not inhibit rebound firing. Removal of gM eliminates firing (curve 2) as does removal of both currents (curve 3).

During hyperpolarizing current injection, a depolarizing voltage “sag” was often observed. When present, a voltage sag was easily detected with membrane hyperpolarization beyond −100 mV (Fig. 8Bi) but was also observed at less negative hyperpolarization (Fig. 8Ci). We found a voltage sag in 17 of 28 cells hyperpolarized to at least −100 mV from a holding potential of −70 mV. This phenomenon has been previously reported in mouse tSPNs (Jobling and Gibbins, 1999) and other mammalian sympathetic neurons (Cassell et al., 1986) where it has been attributed to the anomalous rectifier, or H-current (IH). To support a role for IH, this conductance was implemented in the computational model and was found to reproduce the observed voltage sag (Fig. 8Bii). IH has also been shown to contribute to a more depolarized membrane potential (Pape, 1996; Lamas, 1998), so we compared resting membrane potential in cells with (n = 17) and without (n = 11) evidence of IH but found no significant differences (Student’s t test, two-tailed, p = 0.17).

IH has also been implicated in post-inhibitory rebound firing (Pape, 1996; Ascoli et al., 2010; Engbers et al., 2011; Ferrante et al., 2017). Sag was seen in 12 of 13 tSPNs exhibiting rebound firing, but rebound firing was only observed when cells were held closer to firing threshold (between −60 and −50 mV; Fig. 8Ci) where IM has been shown to be responsible for inducing a voltage sag and rebound firing (Constanti and Galvan, 1983). We used a computational model to understand the relative contributions of IH and IM and determined that sag is due to IH for significant hyperpolarizations, and IM for more moderate hyperpolarization. Rebound firing can occur in the absence of IH but does not occur in the absence of IM following release from moderate (∼10 mV) hyperpolarization, indicating that IH is neither necessary nor sufficient to induce rebound firing in tSPNs (Fig. 8Cii).

Discussion

Reappraisal of physiologic consequence of passive membrane properties

We obtained high-quality recordings of mouse tSPNs and built a computational model to provide mechanistic insight into their function. Whole-cell recordings preserve membrane properties and provide an accurate representation of tSPN function. This is critically important, as the impalement conductance introduced by microelectrode recordings can change passive membrane properties (Staley et al., 1992; Cymbalyuk et al., 2002; Springer et al., 2015), reduce apparent excitability, underestimate the importance of synaptic convergence (Karila and Horn, 2000; Horn and Kullmann, 2007), and prevent repetitive firing (Springer et al., 2015).

Input resistance (Rin) and membrane time constant (τm) were highly correlated and their values, as well as rheobase, occupy an approximately 10-fold range. Values of Rin and τm are an order of magnitude larger than values previously obtained from the same population using traditional microelectrode recordings (Blackman and Purves, 1969; Jobling and Gibbins, 1999), which indicates that the excitability of tSPNs has been substantially underestimated. Measured cell diameters in the T5 ganglion occupied a five-fold range (cf. Jobling and Gibbins, 1999). Capacitance (Cm) values occupied a three-fold range, and were unrelated to cell recruitment. The strong observed relationship between Rin and measures of firing threshold (i.e., rheobase) demonstrate that ohmic processes dominate tSPN recruitment. These observations suggest that membrane resistivity rather than cell size is the primary determinant of recruitment threshold across the population (Gustafsson and Pinter, 1984), although it is unclear whether the observed variability in excitability represents a population recruitment principle.

The preservation of the passive membrane electrical properties Rin and τm leads to synaptic events of greater amplitude and longer duration, which has important consequences for synaptic recruitment. Paravertebral neurons receive nicotinic EPSPs comprising both sub- and suprathreshold events of variable amplitude (Nishi and Koketsu, 1960; Blackman and Purves, 1969; Karila and Horn, 2000; Bratton et al., 2010). An overall increase in EPSP amplitude would convert many subthreshold events into suprathreshold events, thereby increasing tSPN firing rate (Bratton et al., 2010). Traditionally, summation of EPSPs was not thought to contribute to cell recruitment in paravertebral ganglia (North, 1986; McLachlan et al., 1997; Jänig, 2006). However, recent whole-cell recordings from rat SCG demonstrate long-duration sEPSPs with much greater capacity for summation (Springer et al., 2015). We also observed long-duration sEPSPs with decay time constant comparable to τm and examples of sEPSP summation leading to cell recruitment. This provides direct support for the gain hypothesis for amplification of preganglionic activity (Karila and Horn, 2000; Horn and Kullmann, 2007). The observed τm values indicate that tSPNs could act as integrators during states of strong preganglionic sympathetic drive from individual neurons (Jänig, 1985; Ivanov and Purves, 1989) and could widen the temporal window for coincidence detection and summation of convergent synchronous preganglionic inputs (Skok, 1973; König et al., 1996; Ratté et al., 2013). These observations support the concept that tSPNs do not merely relay preganglionic activity, but rather actively integrate and amplify sympathetic output. Metabotropic receptor-mediated changes in intrinsic membrane conductances may further amplify this process (North, 1986; Karila and Horn, 2000).

Additionally, important was the observation that all tSPNs were capable of firing repetitively, which contrasts traditional observations in all paravertebral neurons, including tSPNs, of phasic firing in response to sustained current injection (Jobling and Gibbins, 1999; Jänig, 2006; Springer et al., 2015). Recent whole-cell recordings in rat SCG similarly found paravertebral neurons were capable of repetitive firing, and suggested the discrepancy was a result of impalement conductance (Springer et al., 2015). We were able to replicate these results using our model; by introducing an impalement conductance consistent with microelectrode impalement, we were able to convert repetitively firing model neurons to phasically firing model neurons. Phasic firing after impalement injury appears to be dependent on the presence of IM, as blocking IM can convert sympathetic neurons from phasic to repetitively firing (Brown and Adams, 1980; Brown and Constanti, 1980; Cassell et al., 1986; Luther and Birren, 2009). This observation was reproduced by subtracting IM in our model. IKCa has also been shown to contribute to the interconversion of sympathetic neuron membrane firing properties (Sacchi et al., 1995; Luther and Birren, 2009). Thus, the firing properties of paravertebral sympathetic neurons that exhibit IM and IKCa are particularly sensitive to impalement leak, which underscores the importance of using whole-cell recordings. Blackman and colleagues were able to observe repetitive firing with microelectrodes, a finding that has been consistently overlooked (Blackman and Purves, 1969). A possible explanation could be differences in ion channel expression between the mouse and guinea pig.

The physiologic relevance of repetitive firing in tSPNs

The physiologic relevance of repetitive firing in tSPNs in response to current stimulation might be dismissed if one assumes that postganglionic neurons are only driven by nicotinic preganglionic input. However, paravertebral neurons can exhibit long-lasting depolarization and sustained firing (Blackman and Purves, 1969; Jänig et al., 1982; Kawatani et al., 1987). Activation of metabotropic muscarinic and various other non-cholinergic receptors are implicated (Jänig et al., 1982; North, 1986; Kawatani et al., 1987; Elfvin et al., 1993). These studies support the idea that tSPNs can generate sustained sympathetic drive with limited influence from preganglionics.

Passive membrane properties and various conductances are responsible for sculpting the firing response of tSPNs. Rin is important in determining firing rate over a range of injected current values. Rin also impacts the slope of the ƒ-I curve. tSPNs with steeper slope may be more effective at amplifying postganglionic output gain (Salinas and Thier, 2000; Zimmerman and Hochman, 2010). Given the relatively low steady-state firing rates of preganglionic neurons observed in vivo (Jänig, 2006), the physiologic relevance of variability in response amplification is unclear (Springer et al., 2015).

However, synaptic drive may contribute to response amplification during bouts of metabotropic receptor-mediated sustained activity described above. We observe spontaneous EPSCs with amplitudes ranging from 10 pA to over 100 pA (data not shown). Comparing these amplitudes to values of rheobase (range, 5–70 pA) supports conditions where synaptic actions are capable of transient response amplification.

Relating observed cellular properties to underlying conductances

While the firing rate of tSPNs is strongly determined by the temporal dynamics of the fAHP, a feature carried by IA and IK in rodent SCG (Belluzzi and Sacchi, 1988), the mechanisms underlying SRA have not been studied in paravertebral ganglia including tSPNs. SRA has been well characterized elsewhere (Benda and Herz, 2003; Benda and Tabak, 2013). Contributions from IM and IKCa are among the proposed mechanisms (Sawczuk et al., 1997; Powers et al., 1999; Miles et al., 2005; Yi et al., 2015), and these currents have been previously identified in rodent paravertebral ganglia (Sacchi et al., 1995; Davies et al., 1996; Haley et al., 2000; Locknar et al., 2004; Maingret et al., 2008). Our modeling found that IM and IKCa were required to replicate the fast and slow components of SRA, respectively. IM and IKCa are also known to contribute to the sAHP in rodent SCG and hippocampus (Storm, 1990; Sacchi et al., 1995), and inclusion of IKCa or IM in the model reproduced the sAHP after repetitive firing. That SRA ratio and sAHP half-decay were correlated further supports co-involvement of these conductances.

Other factors contributing to modulation of tSPN excitability

tSPNs are known to express IA, IH and IM (Jobling and Gibbins, 1999). These currents have been shown to modulate EPSP amplitude, synaptic integration, membrane potential, and repetitive firing rate (Connor and Stevens, 1971; Storm, 1990; Rush and Rinzel, 1995; Hoffman et al., 1997; Lamas, 1998; Prescott et al., 2006; George et al., 2009; Kullmann et al., 2016). We found evidence of IA, IH and IM in our recordings by observing phenomena such as notch, sag, and rebound firing, and we replicated their effects using computational modeling. These phenomena typically require hyperpolarization to emerge. While there are no known inhibitory synapses in sympathetic ganglia (McLachlan, 2007), a slow IPSP due to metabotropic activation of K+ conductances has been observed in SCG (Libet and Kobayashi, 1974; North, 1986). Another method of hyperpolarization observed in a small group of tSPNs is the slowly developing usAHP that follows prolonged activity (Zhang and Sillar, 2012). The usAHP has been observed in rabbit SCG (Lees and Wallis, 1974) and reflects Na+-dependent activation of the ouabain-sensitive α3 Na+/K+-ATPase (Picton et al., 2017). These long-lasting hyperpolarizations may provide a physiologic mechanism by which the aforementioned phenomena may emerge.

Extended Data 1: Computational model code

Python and MATLAB code for the computational model of tSPN. Documentation is provided within the code. Download Extended Data 1, ZIP file.

Footnotes

  • The authors declare no competing financial interests.

  • This work was supported by the National Institutes of Health Grant 5R01NS102871 and the Department of Defense Grant SCI-30225. M.L.M. was supported by a fellowship from the National Science Foundation (GRFP), and K.T. was supported by a fellowship from the Emory Center for Mind, Brain, and Culture (CMBC).

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    Ascoli GA, Gasparini S, Medinilla V, Migliore M (2010) Local control of postinhibitory rebound spiking in CA1 pyramidal neuron dendrites. J Neurosci 30:6434–6442. doi:10.1523/JNEUROSCI.4066-09.2010 pmid:20445069
    OpenUrlAbstract/FREE Full Text
  2. ↵
    Bartness TJ, Vaughan CH, Song CK (2010) Sympathetic and sensory innervation of brown adipose tissue. Int J Obes (Lond) 34 [Suppl 1]:S36–S42. doi:10.1038/ijo.2010.182 pmid:20935665
    OpenUrlCrossRefPubMed
  3. ↵
    Belluzzi O, Sacchi O (1988) The interactions between potassium and sodium currents in generating action potentials in the rat sympathetic neurone. J Physiol 397:127–147. pmid:2457694
    OpenUrlCrossRefPubMed
  4. ↵
    Benda J, Herz AV (2003) A universal model for spike-frequency adaptation. Neural Comput 15:2523–2564. doi:10.1162/089976603322385063 pmid:14577853
    OpenUrlCrossRefPubMed
  5. ↵
    Benda J, Tabak J (2013) Spike-frequency adaptation. In: Encyclopedia of computational neuroscience ( Jaeger D, Jung R , eds), pp 1–12. New York: Springer New York.
  6. ↵
    Bhalla US, Bower JM (1993) Exploring parameter space in detailed single neuron models: simulations of the mitral and granule cells of the olfactory bulb. J Neurophysiol 69:1948–1965. doi:10.1152/jn.1993.69.6.1948 pmid:7688798
    OpenUrlCrossRefPubMed
  7. ↵
    Blackman JG, Purves RD (1969) Intracellular recordings from ganglia of the thoracic sympathetic chain of the guinea-pig. J Physiol 203:173–198. pmid:5821870
    OpenUrlPubMed
  8. ↵
    Bratton B, Davies P, Jänig W, McAllen R (2010) Ganglionic transmission in a vasomotor pathway studied in vivo: vasomotor ganglionic transmission in vivo . J Physiol 588:1647–1659. doi:10.1113/jphysiol.2009.185025
    OpenUrlCrossRefPubMed
  9. ↵
    Brown AG (1981) Organization in the spinal cord: the anatomy and physiology of identified neurons. London: Springer.
  10. ↵
    Brown DA (1988) M currents. In: Ion channels ( Narahashi T , ed), pp 55–94. Boston, MA: Springer.
  11. ↵
    Brown DA, Constanti A (1980) Intracellular observations on the effects of muscarinic agonists on rat sympathetic neurones. Br J Pharmacol 70:593–608. pmid:7470731
    OpenUrlPubMed
  12. ↵
    Brown DA, Adams PR (1980) Muscarinic suppression of a novel voltage-sensitive K+ current in a vertebrate neurone. Nature 283:673–676. doi:10.1038/283673a0
    OpenUrlCrossRefPubMed
  13. ↵
    Brownstone RM, Jordan LM, Kriellaars DJ, Noga BR, Shefchyk SJ (1992) On the regulation of repetitive firing in lumbar motoneurones during fictive locomotion in the cat. Exp Brain Res 90:441–455. pmid:1426105
    OpenUrlCrossRefPubMed
  14. ↵
    Cassell JF, Clark AL, McLachlan EM (1986) Characteristics of phasic and tonic sympathetic ganglion cells of the guinea-pig. J Physiol 372:457–483. pmid:2425087
    OpenUrlCrossRefPubMed
  15. ↵
    Connor JA, Stevens CF (1971) Prediction of repetitive firing behaviour from voltage clamp data on an isolated neurone soma. J Physiol 213:31–53. pmid:5575343
    OpenUrlCrossRefPubMed
  16. ↵
    Constanti A, Galvan M (1983) M-current in voltage-clamped olfactory cortex neurones. Neurosci Lett 39:65–70. pmid:6314193
    OpenUrlCrossRefPubMed
  17. ↵
    Cymbalyuk GS, Gaudry Q, Masino MA, Calabrese RL (2002) Bursting in leech heart interneurons: cell-autonomous and network-based mechanisms. J Neurosci 22:10580–10592. pmid:12486150
    OpenUrlAbstract/FREE Full Text
  18. ↵
    Davies PJ, Ireland DR, McLachlan EM (1996) Sources of Ca2+ for different Ca(2+)-activated K+ conductances in neurones of the rat superior cervical ganglion. J Physiol 495:353–366. doi:10.1113/jphysiol.1996.sp021599
    OpenUrlCrossRefPubMed
  19. ↵
    Eccles JC (1935) The action potential of the superior cervical ganglion. J Physiol 85:179–206. pmid:16994707
    OpenUrlPubMed
  20. ↵
    Elfvin L-G, Lindh B, Hökfelt T (1993) The chemical neuroanatomy of sympathetic ganglia. Annu Rev Neurosci 16:471–507. doi:10.1146/annurev.ne.16.030193.002351 pmid:8384808
    OpenUrlCrossRefPubMed
  21. ↵
    Engbers JD, Anderson D, Tadayonnejad R, Mehaffey WH, Molineux ML, Turner RW (2011) Distinct roles for I(T) and I(H) in controlling the frequency and timing of rebound spike responses. J Physiol 589:5391–5413. doi:10.1113/jphysiol.2011.215632 pmid:21969455
    OpenUrlCrossRefPubMed
  22. ↵
    Ermentrout GB, Terman DH (2010) Mathematical foundations of neuroscience. Berlin: Springer.
  23. ↵
    Erulkar SD, Woodward JK (1968) Intracellular recording from mammalian superior cervical ganglion in situ. J Physiol 199:189–203. doi:10.1113/jphysiol.1968.sp008648
    OpenUrlCrossRefPubMed
  24. ↵
    Ferrante M, Shay CF, Tsuno Y, William Chapman G, Hasselmo ME (2017) Post-inhibitory rebound spikes in rat medial entorhinal layer II/III principal cells: in vivo, in vitro, and computational modeling characterization. Cereb Cortex 27:2111–2125. doi:10.1093/cercor/bhw058 pmid:26965902
    OpenUrlCrossRefPubMed
  25. ↵
    Furlan A, La Manno G, Lübke M, Häring M, Abdo H, Hochgerner H, Kupari J, Usoskin D, Airaksinen MS, Oliver G, Linnarsson S, Ernfors P (2016) Visceral motor neuron diversity delineates a cellular basis for nipple- and pilo-erection muscle control. Nat Neurosci 19:1331–1340. doi:10.1038/nn.4376
    OpenUrlCrossRef
  26. ↵
    Galvan M, Sedlmeir C (1984) Outward currents in voltage-clamped rat sympathetic neurones. J Physiol 356:115–133. pmid:6097667
    OpenUrlCrossRefPubMed
  27. ↵
    George MS, Abbott LF, Siegelbaum SA (2009) HCN hyperpolarization-activated cation channels inhibit EPSPs by interactions with M-type K(+) channels. Nat Neurosci 12:577–584. doi:10.1038/nn.2307 pmid:19363490
    OpenUrlCrossRefPubMed
  28. ↵
    Gibbins IL (1991) Vasomotor, pilomotor and secretomotor neurons distinguished by size and neuropeptide content in superior cervical ganglia of mice. J Auton Nerv Syst 34:171–183. pmid:1680891
    OpenUrlCrossRefPubMed
  29. ↵
    Golowasch J, Thomas G, Taylor AL, Patel A, Pineda A, Khalil C, Nadim F (2009) Membrane capacitance measurements revisited: dependence of capacitance value on measurement method in nonisopotential neurons. J Neurophysiol 102:2161–2175. doi:10.1152/jn.00160.2009 pmid:19571202
    OpenUrlCrossRefPubMed
  30. ↵
    Gustafsson B, Pinter MJ (1984) An investigation of threshold properties among cat spinal alpha-motoneurones. J Physiol 357:453. pmid:6512700
    OpenUrlCrossRefPubMed
  31. ↵
    Haley JE, Delmas P, Offermanns S, Abogadie FC, Simon MI, Buckley NJ, Brown DA (2000) Muscarinic inhibition of calcium current and M current in Galpha q-deficient mice. J Neurosci 20:3973–3979. pmid:10818132
    OpenUrlAbstract/FREE Full Text
  32. ↵
    Hines ML, Morse T, Migliore M, Carnevale NT, Shepherd GM (2004) ModelDB: a database to support computational neuroscience. J Comput Neurosci 17:7–11. doi:10.1023/B:JCNS.0000023869.22017.2e pmid:15218350
    OpenUrlCrossRefPubMed
  33. ↵
    Hochman S, McCrea DA (1994) Effects of chronic spinalization on ankle extensor motoneurons. II. Motoneuron electrical properties. J Neurophysiol 71:1468–1479. doi:10.1152/jn.1994.71.4.1468 pmid:8035228
    OpenUrlCrossRefPubMed
  34. ↵
    Hoffman DA, Magee JC, Colbert CM, Johnston D (1997) K+ channel regulation of signal propagation in dendrites of hippocampal pyramidal neurons. Nature 387:869–875. doi:10.1038/43119 pmid:9202119
    OpenUrlCrossRefPubMed
  35. ↵
    Horn JP, Kullmann PHM (2007) Dynamic clamp analysis of synaptic integration in sympathetic ganglia. Neurophysiology 39:423–429. doi:10.1007/s11062-008-9002-y
    OpenUrlCrossRef
  36. ↵
    Ivanov A, Purves D (1989) Ongoing electrical activity of superior cervical ganglion cells in mammals of different size. J Comp Neurol 284:398–404. doi:10.1002/cne.902840307
    OpenUrlCrossRefPubMed
  37. ↵
    Jänig W (1985) Organization of the lumbar sympathetic outflow to skeletal muscle and skin of the cat hindlimb and tail. Rev Physiol Biochem Pharmacol 102:119–213.
    OpenUrlCrossRefPubMed
  38. ↵
    Jänig W, Krauspe R, Wiedersatz G (1982) Transmission of impulses from pre- to postganglionic vasoconstrictor and sudomotor neurons. J Auton Nerv Syst 6:95–106. pmid:6290564
    OpenUrlCrossRefPubMed
  39. ↵
    Jänig W (2006) The integrative action of the autonomic nervous system. New York: Cambridge University Press.
  40. ↵
    Jobling P, Gibbins IL (1999) Electrophysiological and morphological diversity of mouse sympathetic neurons. J Neurophysiol 82:2747–2764. doi:10.1152/jn.1999.82.5.2747 pmid:10561442
    OpenUrlCrossRefPubMed
  41. ↵
    Karila P, Horn JP (2000) Secondary nicotinic synapses on sympathetic B neurons and their putative role in ganglionic amplification of activity. The J Neurosci 20:908–918. pmid:10648695
    OpenUrlAbstract/FREE Full Text
  42. ↵
    Kawatani M, Rutigliano MJ, de Groat WC (1987) Vasoactive intestinal polypeptide facilitates the late component of the 5-hydroxytryptamine-induced discharge in the cat superior cervical ganglion. Neurosci Lett 73:59–64. pmid:2882445
    OpenUrlPubMed
  43. ↵
    König P, Engel AK, Singer W (1996) Integrator or coincidence detector? The role of the cortical neuron revisited. Trends Neurosci 19:130–137. pmid:8658595
    OpenUrlCrossRefPubMed
  44. ↵
    Kullmann PH, Sikora KM, Clark KL, Arduini I, Springer MG, Horn JP (2016) HCN hyperpolarization-activated cation channels strengthen virtual nicotinic EPSPs and thereby elevate synaptic amplification in rat sympathetic neurons. J Neurophysiol 116:438–447. doi:10.1152/jn.00223.2016 pmid:27146984
    OpenUrlCrossRefPubMed
  45. ↵
    Kurian M, Crook SM, Jung R (2011) Motoneuron model of self-sustained firing after spinal cord injury. J Comput Neurosci 31:625–645. doi:10.1007/s10827-011-0324-1 pmid:21526348
    OpenUrlCrossRefPubMed
  46. ↵
    Lamas JA (1998) A hyperpolarization-activated cation current (Ih) contributes to resting membrane potential in rat superior cervical sympathetic neurones. Pflugers Arch 436:429–435. pmid:9644226
    OpenUrlCrossRefPubMed
  47. ↵
    Lees GM, Wallis DI (1974) Hyperpolarization of rabbit superior cervical ganglion cells due to activity of an electrogenic sodium pump. Br J Pharmacol 50:79–93. pmid:4823465
    OpenUrlPubMed
  48. ↵
    Li C, Horn JP (2006) Physiological classification of sympathetic neurons in the rat superior cervical ganglion. J Neurophysiol 95:187–195. doi:10.1152/jn.00779.2005 pmid:16177176
    OpenUrlCrossRefPubMed
  49. ↵
    Libet B, Kobayashi H (1974) Adrenergic mediation of slow inhibitory postsynaptic potential in sympathetic ganglia of the frog. J Neurophysiol 37:805–814. doi:10.1152/jn.1974.37.4.805 pmid:4546163
    OpenUrlCrossRefPubMed
  50. ↵
    Lichtman JW, Purves D, Yip JW (1980) Innervation of sympathetic neurones in the guinea-pig thoracic chain. J Physiol 298:285–299. pmid:7359403
    OpenUrlPubMed
  51. ↵
    Locknar SA, Barstow KL, Tompkins JD, Merriam LA, Parsons RL (2004) Calcium-induced calcium release regulates action potential generation in guinea-pig sympathetic neurones. J Physiol 555:627–635. doi:10.1113/jphysiol.2003.059485 pmid:14724192
    OpenUrlCrossRefPubMed
  52. ↵
    Lu C, Fuchs E (2014) Sweat gland progenitors in development, homeostasis, and wound repair. Cold Spring Harb Perspect Med 4.
  53. ↵
    Luther JA, Birren SJ (2009) p75 and TrkA signaling regulates sympathetic neuronal firing patterns via differential modulation of voltage-gated currents. J Neurosci 29:5411–5424. doi:10.1523/JNEUROSCI.3503-08.2009 pmid:19403809
    OpenUrlAbstract/FREE Full Text
  54. ↵
    Maingret F, Coste B, Hao J, Giamarchi A, Allen D, Crest M, Litchfield DW, Adelman JP, Delmas P (2008) Neurotransmitter modulation of small-conductance Ca2+-activated K+ channels by regulation of Ca2+ gating. Neuron 59:439–449. doi:10.1016/j.neuron.2008.05.026 pmid:18701069
    OpenUrlCrossRefPubMed
  55. ↵
    Masliukov PM, Timmermans JP (2004) Immunocytochemical properties of stellate ganglion neurons during early postnatal development. Histochem Cell Biol 122:201–209. doi:10.1007/s00418-004-0692-y pmid:15338227
    OpenUrlCrossRefPubMed
  56. ↵
    McLachlan EM (2003) Transmission of signals through sympathetic ganglia–modulation, integration or simply distribution? Acta Physiol Scand 177:227–235. doi:10.1046/j.1365-201X.2003.01075.x pmid:12608993
    OpenUrlCrossRefPubMed
  57. ↵
    McLachlan EM (2007) Diversity of sympathetic vasoconstrictor pathways and their plasticity after spinal cord injury. Clin Auton Res 17:6–12. doi:10.1007/s10286-006-0394-8 pmid:17264977
    OpenUrlCrossRefPubMed
  58. ↵
    McLachlan EM, Davies PJ, Häbler HJ, Jamieson J (1997) On-going and reflex synaptic events in rat superior cervical ganglion cells. J Physiol 501:165–181. doi:10.1111/j.1469-7793.1997.165bo.x
    OpenUrlCrossRefPubMed
  59. ↵
    McLachlan EM, Häbler HJ, Jamieson J, Davies PJ (1998) Analysis of the periodicity of synaptic events in neurones in the superior cervical ganglion of anaesthetized rats. J Physiol 511:461–478. doi:10.1111/j.1469-7793.1998.461bh.x
    OpenUrlCrossRefPubMed
  60. ↵
    Miles GB, Dai Y, Brownstone RM (2005) Mechanisms underlying the early phase of spike frequency adaptation in mouse spinal motoneurones. J Physiol 566:519–532. doi:10.1113/jphysiol.2005.086033 pmid:15878947
    OpenUrlCrossRefPubMed
  61. ↵
    Nishi S, Koketsu K (1960) Electrical properties and activities of single sympathetic neurons in frogs. J Cell Comp Physiol 55:15–30. pmid:14427308
    OpenUrlPubMed
  62. ↵
    North RA (1986) Mechanisms of autonomic integration. In: The nervous system, intrinsic regulatory systems of the brain, pp 115–153. New York: Oxford University Press.
  63. ↵
    Pape HC (1996) Queer current and pacemaker: the hyperpolarization-activated cation current in neurons. Annu Rev Physiol 58:299–327. doi:10.1146/annurev.physiol.58.1.299
    OpenUrlCrossRefPubMed
  64. ↵
    Picton LD, Nascimento F, Broadhead MJ, Sillar KT, Miles GB (2017) Sodium pumps mediate activity-dependent changes in mammalian motor networks. J Neurosci 37:906–921. doi:10.1523/JNEUROSCI.2005-16.2016 pmid:28123025
    OpenUrlAbstract/FREE Full Text
  65. ↵
    Platkiewicz J, Brette R (2010) A threshold equation for action potential initiation. PLoS Comput Biol 6:e1000850. doi:10.1371/journal.pcbi.1000850 pmid:20628619
    OpenUrlCrossRefPubMed
  66. ↵
    Powers RK, Sawczuk A, Musick JR, Binder MD (1999) Multiple mechanisms of spike-frequency adaptation in motoneurones. J Physiol 93:101–114. doi:10.1016/S0928-4257(99)80141-7
    OpenUrlCrossRefPubMed
  67. ↵
    Prescott SA, Ratté S, De Koninck Y, Sejnowski TJ (2006) Nonlinear interaction between shunting and adaptation controls a switch between integration and coincidence detection in pyramidal neurons. J Neurosci 26:9084–9097. doi:10.1523/JNEUROSCI.1388-06.2006 pmid:16957065
    OpenUrlAbstract/FREE Full Text
  68. ↵
    Prinz AA, Bucher D, Marder E (2004) Similar network activity from disparate circuit parameters. Nat Neurosci 7:1345–1352. doi:10.1038/nn1352 pmid:15558066
    OpenUrlCrossRefPubMed
  69. ↵
    Purves D, Wigston DJ (1983) Neural units in the superior cervical ganglion of the guinea-pig. J Physiol 334:169–178. pmid:6864556
    OpenUrlPubMed
  70. ↵
    Rall W (2011) Core conductor theory and cable properties of neurons. In: Handbook of physiology, the nervous system, cellular biology of neurons, pp 39–97. Bethesda: American Physiological Society.
  71. ↵
    Ratté S, Hong S, De Schutter E, Prescott SA (2013) Impact of neuronal properties on network coding: roles of spike initiation dynamics and robust synchrony transfer. Neuron 78:758–772. doi:10.1016/j.neuron.2013.05.030 pmid:23764282
    OpenUrlCrossRefPubMed
  72. ↵
    Rimmer K, Horn JP (2010) Weak and straddling secondary nicotinic synapses can drive firing in rat sympathetic neurons and thereby contribute to ganglionic amplification. Front Neurol 1:130. doi:10.3389/fneur.2010.00130 pmid:21173895
    OpenUrlCrossRefPubMed
  73. ↵
    Rittenhouse AR, Zigmond RE (1999) Role of N- and L-type calcium channels in depolarization-induced activation of tyrosine hydroxylase and release of norepinephrine by sympathetic cell bodies and nerve terminals. J Neurobiol 40:137–148. pmid:10413445
    OpenUrlCrossRefPubMed
  74. ↵
    Rush ME, Rinzel J (1995) The potassium A-current, low firing rates and rebound excitation in Hodgkin-Huxley models. Bull Math Biol 57:899–929. pmid:8528161
    OpenUrlCrossRefPubMed
  75. ↵
    Sacchi O, Rossi ML, Canella R (1995) The slow Ca(2+)-activated K+ current, IAHP, in the rat sympathetic neurone. J Physiol 483:15–27. doi:10.1113/jphysiol.1995.sp020564
    OpenUrlCrossRefPubMed
  76. ↵
    Salinas E, Thier P (2000) Gain modulation: a major computational principle of the central nervous system. Neuron 27:15–21. doi:10.1016/S0896-6273(00)00004-0
    OpenUrlCrossRefPubMed
  77. ↵
    Sawczuk A, Powers RK, Binder MD (1997) Contribution of outward currents to spike-frequency adaptation in hypoglossal motoneurons of the rat. J Neurophysiol 78:2246–2253. doi:10.1152/jn.1997.78.5.2246
    OpenUrlCrossRefPubMed
  78. ↵
    Schäfer MK, Eiden LE, Weihe E (1998) Cholinergic neurons and terminal fields revealed by immunohistochemistry for the vesicular acetylcholine transporter. II. The peripheral nervous system. Neuroscience 84:361–376. pmid:9539210
    OpenUrlCrossRefPubMed
  79. ↵
    Schütz B, Schäfer MK, Gördes M, Eiden LE, Weihe E (2015) Satb2-independent acquisition of the cholinergic sudomotor phenotype in rodents. Cell Mol Neurobiol 35:205–216. doi:10.1007/s10571-014-0113-2 pmid:25239161
    OpenUrlCrossRefPubMed
  80. ↵
    Skok VI (1973) Physiology of autonomic ganglia. Tokyo: Igaku Shoin.
  81. ↵
    Springer MG, Kullmann PHM, Horn JP (2015) Virtual leak channels modulate firing dynamics and synaptic integration in rat sympathetic neurons: implications for ganglionic transmission in vivo. J Physiol 593:803–823. doi:10.1113/jphysiol.2014.284125
    OpenUrlCrossRefPubMed
  82. ↵
    Staley KJ, Otis TS, Mody I (1992) Membrane properties of dentate gyrus granule cells: comparison of sharp microelectrode and whole-cell recordings. J Neurophysiol 67:1346–1358. doi:10.1152/jn.1992.67.5.1346 pmid:1597717
    OpenUrlCrossRefPubMed
  83. ↵
    Stauffer EK, McDonagh JC, Hornby TG, Reinking RM, Stuart DG (2007) Historical reflections on the afterhyperpolarization–firing rate relation of vertebrate spinal neurons. J Comp Physiol A 193:145–158. doi:10.1007/s00359-006-0198-2
    OpenUrlCrossRefPubMed
  84. ↵
    Storm JF (1990) Potassium currents in hippocampal pyramidal cells. Prog Brain Res 83:161–187. pmid:2203097
    OpenUrlCrossRefPubMed
  85. ↵
    Valli P, Zucca G, Botta L, Seghezzi R (1989) Lumbar sympathetic ganglia in man: an electrophysiological study in vitro. J Auton Nerv Syst 28:211–217. pmid:2628463
    OpenUrlPubMed
  86. ↵
    Venance L, Glowinski J (2003) Heterogeneity of spike frequency adaptation among medium spiny neurones from the rat striatum. Neuroscience 122:77–92. pmid:14596850
    OpenUrlCrossRefPubMed
  87. ↵
    Wheeler DW, Kullmann PH, Horn JP (2004) Estimating use-dependent synaptic gain in autonomic ganglia by computational simulation and dynamic-clamp analysis. J Neurophysiol 92:2659–2671. doi:10.1152/jn.00470.2004 pmid:15212430
    OpenUrlCrossRefPubMed
  88. ↵
    Yi G, Wang J, Tsang KM, Wei X, Deng B, Han C (2015) Spike-frequency adaptation of a two-compartment neuron modulated by extracellular electric fields. Biol Cybern 109:287–306. doi:10.1007/s00422-014-0642-2 pmid:25652337
    OpenUrlCrossRefPubMed
  89. ↵
    Zhang HY, Sillar KT (2012) Short-term memory of motor network performance via activity-dependent potentiation of Na+/K+ pump function. Curr Biol 22:526–531. doi:10.1016/j.cub.2012.01.058 pmid:22405867
    OpenUrlCrossRefPubMed
  90. ↵
    Zimmerman A, Hochman S (2010) Heterogeneity of membrane properties in sympathetic preganglionic neurons of neonatal mice: evidence of four subpopulations in the intermediolateral nucleus. J Neurophysiol 103:490–498. doi:10.1152/jn.00622.2009 pmid:19923248
    OpenUrlCrossRefPubMed

Synthesis

Reviewing Editor: Upinder S. Bhalla, National Centre for Biological Sciences, Tata Institute of Fundamental Research;

Decisions are customarily a result of the Reviewing Editor and the peer reviewers coming together and discussing their recommendations until a consensus is reached. When revisions are invited, a fact-based synthesis statement explaining their decision and outlining what is needed to prepare a revision will be listed below. The following reviewer(s) agreed to reveal their identity: Rohit Manchanda. Note: If this manuscript was transferred from JNeurosci and a decision was made to accept the manuscript without peer review, a brief statement to this effect will instead be what is listed below.

The reviewers and I discussed the revised version of this MS and we agreed that it was substantially improved, and the authors had addressed almost all points raised by the reviewers. There is one minor but essential revision that the authors should do:

The authors should drop the classification of cells into three types, as this is neither convincing nor essential to the point of the paper.

For the author's benefit, here are the specific reviewer's comments:

“One point that is not addressed satisfactorily is their classification of these cells into three types based on the peak voltage of the initial spike compared to the peak voltage of subsequent spikes. They state that other AP features were different but none of the basic membrane properties correlated with these features. There is no indication of any functional significance of this classification. It seems to be just as good a way to quantitavely classify neurons as any number of other AP parameters. Moreover, in Figures 7 C and D the data points suggest a continuum rather than discrete groups. Doing a PCA and clustering showing clearly separate clusters would have been more convincing. It is not clear what value this section adds to the study.”

Back to top

In this issue

eneuro: 6 (2)
eNeuro
Vol. 6, Issue 2
March/April 2019
  • Table of Contents
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this eNeuro article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Dramatically Amplified Thoracic Sympathetic Postganglionic Excitability and Integrative Capacity Revealed with Whole-Cell Patch-Clamp Recordings
(Your Name) has forwarded a page to you from eNeuro
(Your Name) thought you would be interested in this article in eNeuro.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Dramatically Amplified Thoracic Sympathetic Postganglionic Excitability and Integrative Capacity Revealed with Whole-Cell Patch-Clamp Recordings
Michael Lee McKinnon, Kun Tian, Yaqing Li, Alan Joel Sokoloff, Meredith Lucy Galvin, Mi Hyun Choi, Astrid Prinz, Shawn Hochman
eNeuro 30 April 2019, 6 (2) ENEURO.0433-18.2019; DOI: 10.1523/ENEURO.0433-18.2019

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Share
Dramatically Amplified Thoracic Sympathetic Postganglionic Excitability and Integrative Capacity Revealed with Whole-Cell Patch-Clamp Recordings
Michael Lee McKinnon, Kun Tian, Yaqing Li, Alan Joel Sokoloff, Meredith Lucy Galvin, Mi Hyun Choi, Astrid Prinz, Shawn Hochman
eNeuro 30 April 2019, 6 (2) ENEURO.0433-18.2019; DOI: 10.1523/ENEURO.0433-18.2019
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Significance Statement
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
    • Synthesis
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • computational model
  • firing properties
  • membrane properties
  • mouse
  • paravertebral ganglia

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

New Research

  • A Very Fast Time Scale of Human Motor Adaptation: Within Movement Adjustments of Internal Representations during Reaching
  • Optogenetic Activation of β-Endorphin Terminals in the Medial Preoptic Nucleus Regulates Female Sexual Receptivity
  • Hsc70 Ameliorates the Vesicle Recycling Defects Caused by Excess α-Synuclein at Synapses
Show more New Research

Neuronal Excitability

  • Novel roles for the GPI-anchor cleaving enzyme, GDE2, in hippocampal synaptic morphology and function
  • Tolerance in Thalamic Paraventricular Nucleus Neurons Following Chronic Treatment of Animals with Morphine
  • Investigating Mechanically Activated Currents from Trigeminal Neurons of Nonhuman Primates
Show more Neuronal Excitability

Subjects

  • Neuronal Excitability
  • Home
  • Alerts
  • Follow SFN on BlueSky
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Latest Articles
  • Issue Archive
  • Blog
  • Browse by Topic

Information

  • For Authors
  • For the Media

About

  • About the Journal
  • Editorial Board
  • Privacy Notice
  • Contact
  • Feedback
(eNeuro logo)
(SfN logo)

Copyright © 2025 by the Society for Neuroscience.
eNeuro eISSN: 2373-2822

The ideas and opinions expressed in eNeuro do not necessarily reflect those of SfN or the eNeuro Editorial Board. Publication of an advertisement or other product mention in eNeuro should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in eNeuro.