Biophysical Properties of Somatic and Axonal Voltage-Gated Sodium Channels in Midbrain Dopaminergic Neurons

Spiking activities of midbrain dopaminergic neurons are critical for key brain functions including motor control and affective behaviors. Voltage-gated Na+ channels determine neuronal excitability and action potential (AP) generation. Previous studies on dopaminergic neuron excitability mainly focused on Na+ channels at the somatodendritic compartments. Properties of axonal Na+ channels, however, remain largely unknown. Using patch-clamp recording from somatic nucleated patches and isolated axonal blebs from the axon initial segment (AIS) of dopaminergic neurons in mouse midbrain slices, we found that AIS channel density is approximately 4–9 fold higher than that at the soma. Similar voltage dependence of channel activation and inactivation was observed between somatic and axonal channels in both SNc and VTA cells, except that SNc somatic channels inactivate at more hyperpolarized membrane potentials (Vm). In both SNc and VTA, axonal channels take longer time to inactivate at a subthreshold depolarization Vm level, but are faster to recover from inactivation than somatic channels. Moreover, we found that immunosignals of Nav1.2 accumulate at the AIS of dopaminergic neurons. In contrast, Nav1.1 and Nav1.6 immunosignals are not detectible. Together, our results reveal a high density of Na+ channels at the AIS and their molecular identity. In general, somatic and axonal channels of both SNc and VTA dopaminergic neurons share similar biophysical properties. The relatively delayed inactivation onset and faster recovery from inactivation of axonal Na+ channels may ensure AP initiation at high frequencies and faithful signal conduction along the axon.

Spiking activities of midbrain dopaminergic neurons are critical for key brain functions including motor control and affective behaviors. Voltage-gated Na + channels determine neuronal excitability and action potential (AP) generation. Previous studies on dopaminergic neuron excitability mainly focused on Na + channels at the somatodendritic compartments. Properties of axonal Na + channels, however, remain largely unknown. Using patch-clamp recording from somatic nucleated patches and isolated axonal blebs from the axon initial segment (AIS) of dopaminergic neurons in mouse midbrain slices, we found that AIS channel density is approximately 4-9 fold higher than that at the soma. Similar voltage dependence of channel activation and inactivation was observed between somatic and axonal channels in both SNc and VTA cells, except that SNc somatic channels inactivate at more hyperpolarized membrane potentials (V m ). In both SNc and VTA, axonal channels take longer time to inactivate at a subthreshold depolarization V m level, but are faster to recover from inactivation than somatic channels. Moreover, we found that immunosignals of Nav1.2 accumulate at the AIS of dopaminergic neurons. In contrast, Nav1.1 and Nav1.6 immunosignals are not detectible. Together, our results reveal a high density of Na + channels at the AIS and their molecular identity. In general, somatic and axonal channels of both SNc and VTA dopaminergic neurons share similar biophysical properties. The relatively delayed inactivation onset and faster recovery from inactivation of axonal Na + channels may ensure AP initiation at high frequencies and faithful signal conduction along the axon.

INTRODUCTION
Dopaminergic neurons in the midbrain play critical roles in many key brain functions including motor control, motivation, reward, and cognition (Nieoullon, 2002;Wise, 2004;Fields, 2007;Tobler et al., 2007). Dysfunction of dopaminergic systems will lead to major brain disorders, such as Parkinson's disease, depression, schizophrenia, and addiction (Koob, 1998;Lewis and Lieberman, 2000;Dauer and Przedborski, 2003). There are two distinct groups of dopaminergic neurons located in substantia nigra pars compacta (SNc) and the ventral tegmental area (VTA). SNc neurons are essential circuit elements for controlling movement and learning new motor skills, while VTA neurons play a central role in emotive and motivational brain functions. In comparison with those in VTA, dopaminergic neurons in SNc are more vulnerable to degeneration in Parkinson's disease.
Although they locate at different midbrain regions, dopaminergic neurons share similar electrophysiological properties. They discharge action potentials (APs) with broad voltage waveform in a pacemaking manner at a low frequency of 1-5 Hz and are able to generate bursts of APs at higher frequencies (10-30 Hz) in vivo (Grace et al., 2007;Tsai et al., 2009;Cao et al., 2010) and in vitro (Grace and Bunney, 1984;Hounsgaard et al., 1992;Johnson et al., 1992;Neuhoff et al., 2002;Korotkova et al., 2003;Kuwahara et al., 2006;Margolis et al., 2006;Khaliq and Bean, 2010). It has been proposed that tonic firing may determine the background dopamine level, while phasic (burst) firing is responsible for a large amount of dopamine release in response to reward signals (Grace et al., 2007;Tsai et al., 2009;Bermudez and Schultz, 2010). Optogenetic manipulation of the spiking activity in dopaminergic neurons would produce substantial behavior changes (Chaudhury et al., 2013;Tye et al., 2013;Gunaydin et al., 2014). Therefore, it is of interest to know what determines the excitability (i.e., AP generation) of dopaminergic neurons. Previous studies revealed differential mechanisms for pacemaking in these neurons. Spontaneous firing in VTA neurons is determined by a background Na + current (Khaliq and Bean, 2010), whereas that in SNc neurons is driven by a subthreshold Ca 2+ current (Chan et al., 2007(Chan et al., , 2009). These findings suggest that ionic mechanisms governing AP generation in dopaminergic neurons could be different.
Voltage-gated Na + channels (Nav) are critical for AP generation in CNS neurons. Biophysical properties of Nav in dopaminergic neurons differ from those of midbrain GABAergic neurons (Seutin and Engel, 2010;Ding et al., 2011). The predominant pore-forming α subunits in the brain are Nav1.1, Nav1.2, and Nav1.6. Subunit Nav1.3 is also expressed but mainly during earlier development (Shah et al., 2001). Different channel subunits may have distinct biophysical properties. For example, Nav1.6 activates at much lower membrane potential (V m ) levels than Nav1.2 (Rush et al., 2005). In neocortical pyramidal cells, Nav1.6 accumulates at distal regions of the axon initial segment (AIS) and contributes to AP initiation, whereas Nav1.2 concentrates at peri-somatic regions and promotes AP backpropagation to the soma and dendrites (Hu et al., 2009). Previous findings from single-cell RT-PCR suggest expression of Nav1.1, Nav1.2, and Nav1.6 in substantia nigral dopaminergic neurons (Ding et al., 2011). It remains unknown whether midbrain dopaminergic neurons express these channel subunits at the protein level and whether there is any difference in the biophysical properties between SNc and VTA cells. Since AIS is the AP initiation site, it is of interest to examine whether AIS Nav channels differ from those at the soma.
To address these questions, we performed electrophysiological recording and channel subunit antibody staining experiments. We found that, in general, Na + channels of dopaminergic neurons share similar biophysical properties in both SNc and VTA. In comparison with those at soma, axonal channels take longer time to inactivate at a depolarizing V m level but show faster recovery from inactivation. As compared to those at the soma, AIS channels have a much higher density. These features ensure AP initiation at high frequencies at the AIS and faithful conduction. In contrast, properties of somatic channels may help gate AP backpropagation to the somatodendritic compartments. Interestingly, different from previous findings, our immunostaining experiments reveal a high-level expression of Nav1.2 at the AIS of both SNc and VTA dopaminergic neurons, while Nav1.6 and Nav1.1 immunosignals are undetectable.

AP Initiation and Propagation in Midbrain Dopaminergic Neurons
In acute midbrain slices from TH-GFP mice maintained at room temperature (24-26 • C), we performed whole-cell recordings from the soma and/or the axon of GFP-positive cells. Most of the recorded neurons emitted axons from one of their proximal thick dendrites ( Figure 1A), known as the axon-bearing dendrite (ABD). Axon recordings were obtained from the blebs formed at AIS, less than 120 µm away from the ABD branching site. In current-clamp mode, 56% of the recorded SNc cells (n = 10/18) showed spontaneous firing with an average frequency of 1.5 ± 0.6 Hz, 43% of VTA cells discharged spontaneously (n = 6/14) at a frequency of 1.35 ± 1.2 Hz. As shown in Figure 1B, the spontaneous APs of a SNc cell could be detected at both the somatic and axonal recording sites (n = 10 cells). Hyperpolarizing the V m by injecting a small constant negative current (−10 to −90 pA) to the soma would prevent spontaneous firing. Delivery of 500-ms positive current pulses to the axon ( Figure 1C) or the soma ( Figure 1D) could evoke trains of APs at both recording sites.
Close examination of the AP waveforms revealed that the AIS APs preceded the somatic APs no matter whether they are spontaneous and where the current pulses were applied ( Figure 1E), indicating that AP initiated first at the axon. Previous studies (McCormick et al., 2007;Hu et al., 2009) indicated that, if APs initiate near the recording site, the V m will show smooth transition to the AP upstroke because of accumulative activation of voltage-gated Na + channels; however, if APs initiate at a remote site, the V m will first depolarize due to AP propagation and then generate a rapid upstroke due to the recruitment of local Na + channels. Indeed, the first derivative (dV/dt) of individual somatic APs showed two components in the AP rising phase, corresponding to the arrival of AIS potential and the generation of somatodendritic (SD) potential, respectively (blue, Figure 1E). In sharp contrast, there was no bump at the rising phase of AIS APs (red, Figure 1E). These observations agree well with an initiation zone at the AIS (McCormick et al., 2007).
Previous studies in other types of neurons showed that the voltage-gated Na + channels at the AIS possess higher density (Kole et al., 2008;Hu et al., 2009Hu et al., , 2014Li et al., 2014) and unique biophysical properties (Hu et al., 2014), which are critical Dual whole-cell recording from the soma (blue) and the axon (red, 10 µm away from the soma) showing pace making activity at both recording sites. (C) Negative DC current injection to the soma prevents spontaneous firing. Positive current pulses (500 ms in duration) at the AIS cause repetitive firing at the axon and the soma. Note the absence and presence of the initial AP (arrows) in response to weak and strong current pulses (100 vs. 150 pA). (D) Injection of current pulses to the soma evoke repetitive firing at both recording sites. Also note the absence of presence of the initial AP. (E) Top, expansion of the APs indicated in B-D. Axonal APs precede the somatic APs in all cases. Bottom, the first derivative (dV/dt) of the corresponding APs. The two components at the rising phase of the somatic dV/dt trajectory (blue) correspond to the AIS and SD potentials, respectively. Note the difference of the rising phase between somatic and axonal dV/dt traces.
for AP initiation at the AIS. We, therefore, sought to investigate differences in the density and the voltage dependence between somatic and axonal Na + channels.
Together, these results reveal a high density of Na + channels at the AIS in both SNc and VTA cells. Interestingly, the channel density in SNc AIS is significantly higher than that of VTA AIS (p = 0.03), whereas the density of somatic channels in SNc is similar to that of VTA (p = 0.15).
To examine the kinetics of channel activation and inactivation, we fitted the rising and the falling phases of transient currents with single exponential functions (Figures 4A,B). We observed no significant difference in the rise time between somatic and axonal currents evoked by steps to a range of voltage levels (−25 to +20 mV) in both SNc and VTA cells (SNc: n = 8 nucleated patches and 14 isolated blebs; VTA: n = 8 and 13, respectively, Figures 4C,D). However, the decay time constants of axonal FIGURE 2 | Na + channel density at the soma and the AIS of midbrain TH neurons. (A) Fluorescence and DIC images of the recorded somatic nucleated patch and isolated AIS bleb. The recording pipettes contained Alexa Fluo-594, which could label the recorded nucleated patch and axonal bleb.
(B,C) Example Na + currents obtained from SNc (B) and VTA cells (C) by stepping the holding voltage from −120 mV (35 ms) to 0 mV (25 ms). (D) Averaged peak amplitudes of Na + currents. (E) Averaged channel density at the soma and the AIS of SNc and VTA cells. * p < 0.05, comparison between SNc and VTA axonal channel density.
currents were significantly smaller than those of somatic currents for all the tested voltage levels ( Figures 4C,D), indicating that axonal channels inactivate faster than somatic ones. These results indicate that, although somatic and axonal Na + channels share similar voltage dependence, they do differ in some aspects including the channel inactivation kinetics. We next measured the time constant of channel deactivation. Na + channels were activated by stepping the V m from −120 to a 0.2-ms test pulse (0 mV), which was followed by a series of voltage steps (from −100 to −30 mV, Figures 5A,B). We fitted the tail currents with single exponential functions to obtain the time constants at different V m levels ( Figure 5C) and found that axonal channels deactivate more rapidly than somatic channels. The time constant of axonal channels at −30 mV was significantly shorter than that of somatic channels in SNc cells (0.17 ± 0.02 vs. 0.23 ± 0.01 ms, n = 5 nucleated patches and 14 isolated blebs, p = 0.04, Figure 5D). Similar results were observed in VTA cells (soma: 0.27 ± 0.01 ms, n = 5; AIS: 0.17 ± 0.01 ms, n = 14, p = 0.002, Figure 5D).

Steady-State Inactivation and Recovery
Next, we examined the steady-state inactivation (channel availability) at different V m levels. Transient Na + currents were evoked by stepping the V m from a series of conditioning pulses (from −120 to −20 mV, 35 ms) to the 20-ms test pulse at 0 mV (Figures 6A,B). The half-inactivation voltage of SNc somatic channels were −77.5 ± 2.0 mV (n = 7), significantly lower than that of axonal channels (−70.4 ± 2.0 mV, p = 0.03); no significant difference in the slope was observed (9.0 ± 0.8 mV for soma vs. 10.4 ± 1.2 mV for axon, p = 0.4, n = 8, Figures 6A,C,D). For VTA cells, the V 1/2 of inactivation and slope of somatic channels were −65.3 ± 3.1 mV and 8.7 ± 0.3 mV (n = 6), and those of axonal channels were −64.4 ± 2.5 mV (p = 0.9) and 8.4 ± 0.4 mV (n = 12, p = 0.4, Figures 6B,D), respectively. As shown in Figure 6C, the half-inactivation voltage of SNc somatic channels was also significantly hyperpolarized than that of VTA somatic ones (p = 0.01). These results suggest that, at subthreshold V m levels (from −80 to −55 mV), somatic channels of SNc cells show much lower availability for AP generation than their axonal channels; in addition, they are also less available than channels at both the soma and the axon of VTA cells.
To investigate the time course of the onset of channel inactivation, we applied a prepulse to −120 mV for 15 ms and then stepped the V m to −55 mV, followed by a test pulse to 0 mV after a delay ranging from 0 to 38 ms (Figures 7A,B). Single exponential fits of the normalized currents yielded time constants of 8.04 ± 0.52 and 6.40 ± 0.41 ms for SNc (n = 8 cells) and VTA somatic channels (n = 7 cells), respectively, significantly faster than those of axonal channels (SNc: 11.7 ± 1.2 ms, n = 11, p = 0.02; VTA: 10.6 ± 0.8 ms, n = 15, p = 4.0 × 10 −4 , Figures 7C-E).
We then examined the time course of recovery from channel inactivation. After the prepulse (−120 mV), paired pulses with different intervals (from 1 to 110 ms at −120 mV) were applied to the patches (Figure 8). Single exponential fits revealed a tendency of axonal Nav channels to be faster to recover from inactivation than somatic channels (soma: 1.2 ± 0.1 ms, n = 10, vs. axon: 0.9 ± 0.1 ms, n = 14, p = 0.11, Figures 8A,C). Similarly, somatic channels of VTA cells showed a slower recovery from inactivation than axonal channels (soma: 2.2 ± 0.2 ms vs. axon: 0.9 ± 0.1 ms, n = 10 nucleated patches and 19 isolated blebs, p = 1.0 × 10 −5 , Figures 8B,D,E) as well as SNc somatic channels ( Figure 8E). These results suggest that axonal channels recover from inactivation faster than somatic channels, ensuring highfrequency AP generation and conduction.

Molecular Identity of Axonal Na + Channels
We performed triple immunostaining using antibodies of TH, Ankyrin G (AnkG), and different Na + channel α subunits, including Nav1.2, Nav1.1, and Nav1.6. TH and AnkG were used to identify the dopaminergic neurites and AIS, respectively. We performed Nav1.2 antibody staining in the neocortex as a positive control (Hu et al., 2009;Li et al., 2014; Figure 9A). In both SNc and VTA (Figures 9B,C), most of the neurites positive to TH and AnkG were also positive to Nav1.2. Immunosignals of Nav1.2 were found in 93.9% and 96.4% of the neurites positive to both TH and AnkG (n = 46 in SNc and 54 in VTA from 6 mice, Figure 9D). However, none of the TH/AnkG-positive  neurites in both SNc and VTA showed Nav1.1 signals (n = 50 from 3 mice, Figures 10A-C), although positive signals were found in the neocortex ( Figure 10A) and in some of the neurites positive to AnkG but negative to TH (Figures 10B,C). These results indicate that the AIS of dopaminergic neurons express the α subunits Nav1.2 but not Nav1.1. Similar to that in the neocortex (Li et al., 2014) the AIS of TH-negative GABAergic neurons within SNc and VTA may express the subtype Nav1.1. Surprisingly, the most abundant α subunit Nav1.6 in the brain was not found in the AIS of TH-positive neurons (n = 50 from 3 mice, Figures 10D-F). Similar to the pattern of Nav1.1, Nav1.6 signals were found in TH-negative but AnkG-positive neurites, suggesting an expression of this channel subtype in other types of neurons in the midbrain. We found no detectible signals at the soma using channel antibodies, possibly due to the low abundance of channels. The molecular identity of somatic channels remains to be further examined. Together, these results indicate that Nav1.2 is the predominant α subunit in the AIS of midbrain dopaminergic neurons.

DISCUSSION
In this study, we reveal a high density of voltage-gated Na + channels at the AIS of putative midbrain dopaminergic neurons, which may determine the AIS as the AP initiation zone.
Although the activation of voltage-gated Na + channels at the soma and the AIS in midbrain dopaminergic neurons share similar voltage dependence, channel inactivation differs in many aspects. First, the inactivation of transient Na + currents at the axon is much faster than that of somatic channels. Second, at a subthreshold depolarizing V m level, inactivation onset of axonal channels is more delayed than that of somatic channels. Finally, axonal channels show faster recovery from inactivation. These inactivation properties allow higher channel availability, faster kinetics, and more rapid recovery from inactivation, to ensure the generation of APs at higher frequencies and depolarized V m levels. Interestingly, our immunostaining experiments indicate that Nav1.2 is the predominant channel subtype at the AIS that governs AP initiation.

AP Initiation and Channel Density at the AIS
Similar to other types of neurons in the brain including neocortical pyramidal cells (Stuart et al., 1997;Chan et al., 2009) and cerebellar Purkinje cells (Clark et al., 2005), the AP initiation site of midbrain dopaminergic neurons locates at the AIS. Our dual whole-cell recordings at the soma and the AIS in both SNc and VTA neurons indicate that AIS APs always precede the somatic APs no matter whether the APs are generated spontaneously or evoked by current pulses (Figure 1). The phase plots of somatic APs have two components during the upstrokes, corresponding to the AIS potential and SD potential, respectively. However, the upstrokes of AIS APs show smooth initial depolarization and absence of the double components (McCormick et al., 2007), indicating that APs are initiated at the AIS. In agreement with this, simultaneous recordings from the soma and dendrites of dopaminergic neuron show that APs in ABD occur earlier than those at the soma and non-ABD, suggesting that APs originate at the axon (Hausser et al., 1995;Gentet and Williams, 2007). Previous studies suggest roles of specific channel subtypes and high channel density in determining the AP initiation site (Stuart et al., 1997;Li et al., 2014). The Na + channel density at the AIS is about ∼34 fold higher than that at the soma in neocortical pyramidal cells (Hu et al., 2009). The AIS channel density is also higher in neocortical (∼60 fold) (Hu et al., 2009) and hippocampal parvalbumin-expressing neurons (Hu et al., 2009) ensuring fast signal conduction along the axon. In comparison with these types of cells, dopaminergic neurons express relatively less Na + channels at their AIS, the channel density is only 4-9 fold that of somatic ones. Distinct from our estimation of channel density at the soma, previous studies revealed a higher density in rat midbrain cells (Seutin and Engel, 2010;Moubarak et al., 2019). The discrepancy could be attributable to the identification of putative dopaminergic neuron, recording configuration, and animal species. Since dopaminergic neurons are spontaneously active, changes in mean frequency of their tonic firing may control and regulate the background dopamine level in the target brain regions (Sulzer et al., 2016), leading to relatively slow neuromodulation and alteration of brain states. Therefore, precision of individual APs could be not essential for dopaminergic neuron signaling. Activation of AIS channels not only initiates APs but also provides currents to charge the somatodendritic membrane and determine whether or not an AP can successfully backpropagate to the soma and dendrites. Considering that AP backpropagation play critical roles in mediating somatodendritic dopamine release (Geffen et al., 1976;Rice and Patel, 2015) and is subject to modulation by D2 receptors (Cragg and Greenfield, 1997;Gentet and Williams, 2007), we speculate that the lower channel density at the AIS of dopaminergic neurons, as compared with that of pyramidal cells, may ensure a higher sensitivity of AP backpropagation to modulation by dopamine and other neuromodulators.

Voltage Dependence of Channel Activation and Inactivation
We found surprisingly that somatic and AIS Na + channels share similar voltage dependence of channel activation in both SNc and VTA cells (Figure 3). The half-activation voltage shows no significant difference (ranging from −23 to −27 mV). Activation curves of somatic and axonal channels differ slightly in their slopes; small but significant difference was only found in VTA cells (soma: 6.6 mV vs. AIS: 5.4 mV). In addition, activation kinetics of somatic channels are similar to those of AIS channels, as reflected by their similar rising time (Figure 4). The voltage dependence of channel inactivation, however, is more complicated. In the experiments examining the steady-state channel inactivation (availability), we found that SNc somatic channels have the most hyperpolarized half-inactivation voltage than AIS channels in SNc and both somatic and AIS channels in VTA, indicating that at similar V m levels SNc soma has the lowest channel availability (Figure 6). Different from SNc, VTA somatic channels share similar voltage dependence with AIS channels. In general, AIS channels of both SNc and VTA cells show slower onset of inactivation and faster recovery from inactivation as compared with somatic channels (Figures 7, 8).
The difference in the inactivation curves may reflect differences in channel subunits or modulation of channel proteins. Since the expression level of Na + channels at the soma is extremely low and the antibody staining is not sensitive enough, the molecular identity of somatic channels remains unclear. However, our immunostaining results (Figures 9, 10) together with findings from previous in situ hybridization experiments (Gonzalez-Cabrera et al., 2017) indicate that Nav1.2 subunits preferentially express in dopaminergic neurons. In contrast, Nav1.1 and Nav1.6 may only express in other types of neurons, possibly GABAergic inhibitory neurons. Therefore, it is unlikely that the channel subunit difference is attributable to the difference in the inactivation curves of SNc somatic and axonal channels. In comparison to Nav1.6 channels that are refractory to modulation by PKA or PKC (Maurice et al., 2001;Chen et al., 2009), Nav1.2 channels are much more vulnerable to modulation by neuromodulators and protein kinases (Maurice et al., 2001). PKC causes phosphorylation of serine 1506 at the inactivation gate (West et al., 1991) while PKA induces phosphorylation in the I-II loop but not the III-IV loop that is responsible for channel inactivation (Cantrell and Catterall, 2001). Phosphorylation of Nav1.2 induced by PKA and PKC would increase the slow channel inactivation (Carr et al., 2003;Chen et al., 2006Chen et al., , 2008. Therefore, we speculate that the differences in channel inactivation curves as well as those of inactivation and deactivation time courses could be due to differential phosphorylation of somatic and axonal channels. However, it remains unknown what determines the phosphorylation level and whether somatodendritic release of dopamine plays a role in the differential regulation. Previous studies in neocortical pyramidal cells and inhibitory interneurons reveal distinct voltage-dependent properties of Na + channels in the soma and the axon (Hu et al., 2009(Hu et al., , 2014FIGURE 7 | Time course of Na + channel inactivation onset in TH neurons. (A) Example currents evoked by the voltage commands (inset) in SNc neurons. The V m was stepped from a prepulse of -120 mV (20 ms) to -55 mV with varying intervals followed by a test pulse to 0 mV (5 ms). (B) Similar as in A but from VTA neurons. (C,D) Plots of the currents as a function of intervals. Currents were normalized to the Na + currents without time interval. The plots were fitted with single exponential functions. Inset, expanded plots showing differences between the two groups. (E) Time constants for the development of Na + channel inactivation at -55 mV. * p < 0.05; * * * p < 0.001. Li et al., 2014;Ye et al., 2018). The half-activation voltages of channels distributed in distal AIS regions are more hyperpolarized than those in the soma. In midbrain dopaminergic neurons, the similarities in voltage dependence of channel activation at the soma and AIS suggest a more important role of the high channel density and the thin AIS structure in AP initiation (Stuart et al., 1997). In comparison with that of AIS channels of pyramidal cells, the relatively depolarized half-activation voltage of AIS channels would increase the AP threshold and slow down the conduction velocity along the axon. The slower development of channel inactivation and faster recovery ensure higher channel availability in prolonged depolarization levels and with a high firing rate, such as during phasic firing. It remains to be further examined whether biophysical properties of somatic and axonal channels are subject to modulation by neuromodulator receptors.

Molecular Identity of AIS Na + Channels
Previous findings (Ding et al., 2011) from single-cell RT-PCR suggest that dopaminergic neurons express all four Na + channel subtypes in the central nervous system, including Nav1.1, Nav1.2, Nav1.3, and Nav1.6. In addition, quantitative RT-PCR reveal relatively higher expression levels of Nav1.2 in dopaminergic neurons but Nav1.1 and Nav1.6 in GABAergic neurons of the substantial nigra. In situ hybridization experiments also revealed that most of the SNc neurons express mRNA of Nav1.2, only a small population contains mRNA of Nav1.6 (Gonzalez-Cabrera et al., 2017). In agreement with these findings, our immunostaining results reveal that, at the protein level, Nav1.2 is the predominant channel subtype at the AIS of TH-expressing neurons in both SNc and VTA (Figure 9). In sharp contrast, Nav1.1 and Nav1.6 immunosignals are not detectible in THpositive cells (Figure 10). Immunosignals of Nav1.1 and Nav1.6 are found in some TH-negative neurites, most likely from local GABAergic neurons. Considering that Nav1.3 is mainly expressed during earlier developmental stages (Shah et al., 2001) and has a much lower expression level in juvenile animals as compared with other subunits in single-cell RT-PCR experiments (Ding et al., 2011), we did not examine the expression of this channel subtype in the immunostaining experiments. The exclusive expression of Nav1.2 at the AIS of dopaminergic neurons suggest a critical role of Nav1.2 in determining the excitability of these neurons. Dysfunction of channels containing the α subunit Nav1.2 may lead to alteration of brain functions, such as motor control, motivation and reward. Indeed, mutations of Nav1.2 have been associated with autism (Sanders et al., 2012;Ben-Shalom et al., 2017) and epilepsy (Kamiya et al., 2004;Lauxmann et al., 2013;Berecki et al., 2018). Thus, this channel subunit could be a potential drug target for the treatment of related brain disorders.

Slice Preparation
The use and care of animals in this study were approved by the Animal Advisory Committee at the State Key Laboratory of Cognitive Neuroscience and Learning, Beijing Normal University. The number of animals used in the experiments was minimized to comply the rules of ethics committee. We used TH-GFP mice (P16-25) for electrophysiological experiments. Mice were anesthetized by sodium pentobarbital (100 mg/kg, i.p.) before decapitation. The brain was dissected out and immersed in an ice-cold slicing solution (in mM: Sucrose 213, KCl 2.5, NaH 2 PO 4 1.25, NaHCO 3 26, Dextrose 10, MgSO 4 2, CaCl 2 2), which was bubbled with mixed gas (95% O 2 , 5% CO 2 ). In this sucrose-based solution, brain slices (200 µm in thickness) containing the midbrain were cut on a vibratome (VT1200S, Leica). Slices were then transferred to an incubator containing the artificial cerebrospinal fluid (ACSF) (in mM: NaCl 126, KCl 2.5, NaH 2 PO 4 1.25, NaHCO 3 26, Dextrose 25, MgSO 4 2, CaCl 2 2; 315-325 mOsm, pH = 7.2-7.3) and maintained at 34 • C for ∼30-60 min. After incubation, slices were kept in the same solution at the room temperature until use.

Electrophysiological Recordings
To obtain electrophysiological recordings, we transferred slices to a recording chamber perfused with aerated ACSF at 24-26 • C. The GFP-positive cells were considered as putative dopaminergic neurons. They were visualized under an infrared differential interference contrast (IR-DIC) fluorescence microscope (BX51WI, Olympus). Patch pipettes had an impedance of ∼5 M for somatic recording and ∼10 M for axonal bleb recording. Blebs formed at the AIS (60 µm away from the branching point of the ABD) were chosen for recording in this study. The TH-positive axons together with their blebs were identified by the expression of GFP. For current clamp recording, pipettes were filled with a K + -based internal solution (in mM: Kgluconate 140, KCl 3, MgCl2 2, HEPES 10, EGTA 0.2, Na2ATP 2; 295-305 mOsm, pH = 7.2). For voltage clamp recordings, pipettes were filled with a Cs + -based solution (in mM: CsCl 140, MgCl 2 2, HEPES 10, EGTA 10, Na 2 ATP 2; 295-305 mOsm, pH = 7.2). To obtain local Na + currents at the soma (nucleated patch) and the axon (isolated axonal bleb), we added 4-AP (3 mM), TEA-Cl (96 mM) and CdCl 2 (100 µM) to the ACSF but reduced the concentration of NaCl to 50 mM. Isolation of blebs from the axon was achieved as described previously (Chaudhury et al., 2013). In brief, we swept a patch pipette just below the selected GFP-positive bleb and disconnect it from the main axon. A MultiClamp 700B amplifier (Molecular Devices) was used for whole-cell recording and giant patch (somatic nucleated patch and isolated axonal bleb) recording. Voltage and current signals were acquired by pClamp10 software at a sampling rate of 50 kHz. Na + current traces were averaged and then filtered with a Gaussian filter. The liquid junction potentials were not corrected for Vm values shown in the text and figures.
The holding potential of somatic nucleated patches and isolated axonal blebs was −60 mV. For the measurement of Na + channel density and the investigation of voltagedependent channel activation, the V m was stepped from a prepulse (−120 mV, 30 ms) to a series of testing pulses ranging from −70 to +30 mV with increment of 5 mV. The maximal peak amplitude among the evoked families of currents and the membrane area (based on the measured diameter) were used for the estimation of Na + channel conductance density. The peak amplitudes at different voltage levels were normalized to the maximal amplitude and plotted as a function of the membrane potential (activation curves). We fitted the rise and decay of the evoked currents with single exponential functions to examine the kinetics of channel activation and inactivation. For channel deactivation, 0.2 ms after stepping the membrane potential from −120 to 0 mV, a series of 25-ms voltage pulses from −100 to −30 mV were applied to the patches. The deactivation time constants could be obtained from the tail currents at different voltage levels.
To examine the steady-state inactivation of Na + channels, we evoked transient Na + currents by stepping the V m from a series of prepulses (from −120 to −60 mV, 35 ms) to the 20ms test pulse at 0 mV. The evoked currents were normalized to the maximal peak amplitude to generate the inactivation curves, i.e., channel availability at different voltages. To examine the time course of channel inactivation onset, we stepped the (C) Similar as in B but for VTA. (D-F) Similar to A-C but for Nav1.6 staining in the neocortex (D), SNc (E) and VTA (F). Note the absence of Nav1.6 immunosignals from the TH-positive neurites (arrowheads). Arrows indicate neurites positive to AnkG and Nav1.6 but negative to TH. V m from a prepulse at −120 to −55 mV and stepped the V m further to a test pulse of 0 mV after a delay ranging from 0 to 38 ms. Because the recordings from isolated blebs were normally stable for much longer time than that of nucleated patches, we chose to apply relatively longer time intervals in bleb recordings. Channels were supposed to inactivate at −55 mV but in a time-dependent manner. For recovery from inactivation, we applied two paired pulses from −120 to 0 mV (20 ms) with time intervals ranging from 1 to 110 ms. The first pulse in each paired steps lasted for 30 ms. Currents were normalized to the peak amplitude of the current evoked by the first pulse.
We chose to image non-successive sections (about 2 mm (2 mm) for each experiment so that no cell on the surface would be over-counted. Images were taken on a laser scanning confocal microscope (A1plus, Nikon) with 10×, 20×, and 60× objectives. The acquisition parameters were carefully adjusted to make the fluorescence signals linearly displayed and fall into the maximum dynamic range of the detectors. Z-stack images were collected with a voxel interval of 1 µm.

Statistics
All values were presented as mean ± SEM. Statistical significance of difference was examined using non-parametric Wilcoxon rank-sum test.

DATA AVAILABILITY
All datasets generated for this study are included in the manuscript and/or the supplementary files.

ETHICS STATEMENT
This study was carried out in accordance with the recommendations of the Animal Advisory Committee at the State Key Laboratory of Cognitive Neuroscience and Learning, Beijing Normal University. The protocol was also approved by this committee.

AUTHOR CONTRIBUTIONS
YS initiated and designed the experiments. JY performed all the patch-clamp recordings and analyzed the data. YX helped in dual soma-axon recordings. LL and ML conducted the immunostaining experiments. QH helped in data analysis. JY and YS wrote the manuscript.