Two-Photon Na+ Imaging Reports Somatically Evoked Action Potentials in Rat Olfactory Bulb Mitral and Granule Cell Neurites

Dendrodendritic synaptic interactions are a hallmark of neuronal processing in the vertebrate olfactory bulb. Many classes of olfactory bulb neurons including the principal mitral cells (MCs) and the axonless granule cells (GCs) dispose of highly efficient propagation of action potentials (AP) within their dendrites, from where they can release transmitter onto each other. So far, backpropagation in GC dendrites has been investigated indirectly via Ca2+ imaging. Here, we used two-photon Na+ imaging to directly report opening of voltage-gated sodium channels due to AP propagation in both cell types. To this end, neurons in acute slices from juvenile rat bulbs were filled with 1 mM SBFI via whole-cell patch-clamp. Calibration of SBFI signals revealed that a change in fluorescence ΔF/F by 10% corresponded to a Δ[Na+]i of ∼22 mM. We then imaged proximal axon segments of MCs during somatically evoked APs (sAP). While single sAPs were detectable in ∼50% of axons, trains of 20 sAPs at 50 Hz always resulted in substantial ΔF/F of ∼15% (∼33 mM Δ[Na+]i). ΔF/F was significantly larger for 80 Hz vs. 50 Hz trains, and decayed with half-durations τ1/2 ∼0.6 s for both frequencies. In MC lateral dendrites, AP trains yielded small ΔF/F of ∼3% (∼7 mM Δ[Na+]i). In GC apical dendrites and adjacent spines, single sAPs were not detectable. Trains resulted in an average dendritic ΔF/F of 7% (16 mM Δ[Na+]i) with τ1/2 ∼1 s, similar for 50 and 80 Hz. Na+ transients were indistinguishable between large GC spines and their adjacent dendrites. Cell-wise analysis revealed two classes of GCs with the first showing a decrease in ΔF/F along the dendrite with distance from the soma and the second an increase. These classes clustered with morphological parameters. Simulations of Δ[Na+]i replicated these behaviors via negative and positive gradients in Na+ current density, assuming faithful AP backpropagation. Such specializations of dendritic excitability might confer specific temporal processing capabilities to bulbar principal cell-GC subnetworks. In conclusion, we show that Na+ imaging provides a valuable tool for characterizing AP invasion of MC axons and GC dendrites and spines.

Dendrodendritic synaptic interactions are a hallmark of neuronal processing in the vertebrate olfactory bulb. Many classes of olfactory bulb neurons including the principal mitral cells (MCs) and the axonless granule cells (GCs) dispose of highly efficient propagation of action potentials (AP) within their dendrites, from where they can release transmitter onto each other. So far, backpropagation in GC dendrites has been investigated indirectly via Ca 2+ imaging. Here, we used two-photon Na + imaging to directly report opening of voltage-gated sodium channels due to AP propagation in both cell types. To this end, neurons in acute slices from juvenile rat bulbs were filled with 1 mM SBFI via whole-cell patch-clamp. Calibration of SBFI signals revealed that a change in fluorescence F/F by 10% corresponded to a [Na + ] i of ∼22 mM. We then imaged proximal axon segments of MCs during somatically evoked APs (sAP). While single sAPs were detectable in ∼50% of axons, trains of 20 sAPs at 50 Hz always resulted in substantial F/F of ∼15% (∼33 mM [Na + ] i ). F/F was significantly larger for 80 Hz vs. 50 Hz trains, and decayed with half-durations τ 1/2 ∼0.6 s for both frequencies. In MC lateral dendrites, AP trains yielded small F/F of ∼3% (∼7 mM [Na + ] i ). In GC apical dendrites and adjacent spines, single sAPs were not detectable. Trains resulted in an average dendritic F/F of 7% (16 mM [Na + ] i ) with τ 1/2 ∼1 s, similar for 50 and 80 Hz. Na + transients were indistinguishable between large GC spines and their adjacent dendrites. Cell-wise analysis revealed two classes of GCs with the first showing a decrease in F/F along the dendrite with distance from the soma and the second an increase. These classes clustered with morphological parameters. Simulations of [Na + ] i replicated these behaviors via negative and positive gradients in Na + current density, assuming faithful AP backpropagation. Such specializations of dendritic excitability might confer specific temporal processing capabilities to bulbar principal cell-GC subnetworks. In conclusion, we show that Na + imaging provides a valuable tool for characterizing AP invasion of MC axons and GC dendrites and spines.

INTRODUCTION
Many neuronal interactions in the vertebrate olfactory bulb are mediated by dendrodendritic synapses at both stages of the olfactory network, the glomerular input layer and the EPL. Accordingly, the dendrites of many classes of olfactory bulb neurons dispose of several mechanisms that support dendrodendritic processing, including release of transmitter (Sheperd et al., 1990;Wachowiak and Shipley, 2006). In particular, bulbar dendrites often are endowed with high densities of active conductances since they feature an efficient propagation of action potentials.
The principal MCs and the axonless GCs are especially interesting because they interact via dendrodendritic reciprocal synapses and thus both rely on dendritic transmitter release to exchange information. The specific functional consequences of this local interaction are not fully resolved yet. MC dendrites are accessible to whole-cell patch-clamp; in combination with Ca 2+ and voltage-sensitive dye imaging it has been shown that (1) action potentials can be initiated in their apical dendritic tuft and then propagate to the soma and that (2) under certain conditions action potentials will travel far out into their lateral dendrites (Chen et al., 1997;Margrie et al., 2001;Xiong and Chen, 2002;Christie and Westbrook, 2003;Debarbieux et al., 2003;Djurisic et al., 2004).
So far, propagation in GC dendrites has been investigated only indirectly via Ca 2+ imaging (Egger et al., 2003;Zelles et al., 2006), since the small sizes of their dendrites and soma have precluded dendritic patching and voltage-sensitive dye imaging. Therefore other techniques such as Na + imaging might yield useful additional information on their active properties. This approach has been introduced by Jaffe et al. (1992) and Ross et al. (1993), who used the fluorescent Na + -indicator SBFI for detection of Na + influx through voltage-gated sodium channels (Na v ) and thereby provided first direct evidence for active backpropagation of action potentials into dendrites of hippocampal CA1 pyramidal neurons. Later on, two-photon imaging with SBFI demonstrated the existence of Na + transients in response to backpropagating action potentials in dendritic spines of CA1 neurons (Rose et al., 1999). Na + imaging also allows detection of Na + influx into active axons (Kole et al., 2008;Bender and Trussell, 2009;Fleidervish et al., 2010;Scott et al., 2014;Miyazaki and Ross, 2015).
Olfactory bulb GCs represent a special case since they do not bear a classical axon. Rather, their only output occurs onto MC and tufted cell (TC) lateral dendrites from their apical dendrite (Price and Powell, 1970). While this neurite bears spines, it also exhibits many properties of axons for which Na v play a particularly important functional role. GC dendrites are known to dispose of Na v -dependent backpropagation that can be blocked by TTX (Egger et al., 2003), and show dendritic Na + spikelets both in vitro and in vivo (Mori and Takagi, 1978;Wellis and Scott, 1990;Luo and Katz, 2001;Pinato and Midtgaard, 2005;Zelles et al., 2006;Labarrera et al., 2013;Burton and Urban, 2015). There is also accumulating evidence for an essential role of Na v in reciprocal processing at the GC spine, where their local activation triggered by AMPAR-mediated EPSPs can further depolarize the spine to activate the classical voltage-gated presynaptic N/P/Q Ca 2+ channels (Bywalez et al., 2015). Thus it seems likely that GC APs evoke robust Na + signals detectable by Na + imaging in both apical GC dendrites and spines.
Based on this assumption, we chose to probe dendritic Na + signaling in MCs and GCs using both single APs and short AP trains at frequencies of 50 and 80 Hz, which correspond to the maximal frequencies of network oscillations in the olfactory bulb that dendrodendritic MC-GC synapses are involved in (slow and fast gamma range, e.g., Manabe and Mori, 2013). Our main aim was to further elucidate the active properties of GC apical dendrites via direct observation and characterization of AP-evoked Na + signals.

Animals and Slice Preparation
Rats were decapitated under deep anesthesia with isoflurane according to the stipulations of the German law governing animal welfare (Tierschutzgesetz) and according to the EU directive 2010/63/EU, as approved by the Bavarian state government (Regierung von Oberbayern). Brains were removed and horizontal olfactory bulb brain slices (300 µm thick) were prepared of juvenile wild-type (Wistar) of either sex (postnatal days 11 -18). The slices were incubated in artificial cerebrospinal fluid (ACSF, composition: 125 mM NaCl, 26 mM NaHCO 3 , 1.25 mM NaH 2 PO 4 , 20 mM glucose, 2.5 mM KCl, 1 mM MgCl 2 , and 2 mM CaCl 2 ) infused with carbogen gas (95% O 2 , 5% CO 2 ) in a heated water bath at 33 • C for 30 min and then kept at room temperature (22 • C) until experimentation.
All experiments were performed at room temperature (22 • C). The patched MCs and GCs were held in current-clamp mode near their resting potential of −60 and −70 mV, respectively, and the access resistance was monitored. Polarizing step pulses for 500 ms each (first step −90 pA, increased by +50 pA for 10 steps, ending at +180 pA) were applied in order to identify the targeted cell type via its firing pattern. During Na + imaging experiments, cells with a holding current value above −25 pA for GCs and above −50 pA for MCs were rejected.
For Na + imaging, the excitation wavelength of the laser was set to 800 nm and cells were loaded with 1 mM SBFI via the patch pipette. After sufficient filling of the dendritic tree (at least 15 min past establishment of the whole-cell configuration), structures of interest were imaged in free line-scanning mode with a temporal resolution of ∼0.5 ms. At a given dendritic location, several consecutive focal line-scans were recorded during either somatically evoked single action potentials (sAPs) (evoked by an injected current step of 1000 pA) or action potential trains (20 steps at 50 Hz or 80 Hz). To ensure reliable AP induction, current steps had a duration of 5 ms. The duration of the scanning was set to 3.5-5.0 s depending on the decay time course of the Na + transients. Transients recorded sequentially at the same location were averaged and smoothed post hoc.
Post hoc data analysis was performed using custom macros written in IGOR Pro 7 (Wavemetrics, Lake Oswego, OR, USA) and OriginPro (OriginLab Corporation, Northampton, MA, USA). To correct for dye bleaching, Na + imaging was periodically performed without stimulation, during phases of low postsynaptic spontaneous activity. An exponential decay curve was fitted to these data and then subtracted from the traces during which a stimulation was performed.
Dendritic Na + transients, reflected by decreases in fluorescence emission F/F, were analyzed relative to the resting fluorescence F 0 , with their decay measured in terms of half-duration τ 1/2 (Egger et al., 2003). τ 1/2 values were capped at 3 s, because a higher value could not be reliably extrapolated. Rise times of transients were analyzed in terms of the time interval from 20 to 80% of the maximal F/F amplitude.
At the end of individual experiments, the neurites were imaged in a z-stack both in the fluorescence and the infrared channel of the TPLSM. Distance measurements were performed using Fiji's Neurite Simple Tracer plugin (Longair et al., 2011) from z-stack scans of patched cells. The first initial branching of the neurite from the soma was selected as the starting point for the measurement. From that point on, the dendrite was traced in 3D up to the imaged scan line or its first branchpoint and the length of this tracing was determined. The distance of GC somata from the lower border of the MCL was determined in scans imaged in the infrared channel of the TPLSM.

Simulations
The NEURON simulation tool (Release 7.4, Carnevale and Hines, 2006) was used with Python (Hines et al., 2009) to construct a simple compartmental model of the GC apical dendrite for a simulation of Na + entry. Geometric dimensions were based on previous data (mean diameter of distal GC dendrite 1.2 µm, Egger and Stroh, 2009) and on two-photon fluorescence scans of the neurons in this data set (diameter of proximal segment for decreasing cells 2.7 ± 0.2 µm, n = 6 GCs; for increasing cells 2.7 ± 0.4 µm, n = 4). The model consists of 11 connected cylindrical compartments, with the first representing the soma (10 µm diameter and length) and the 10 others a dendrite with linear taper from 2.7 to 1.2 µm or 2.0 to 1.2 µm depending on the cell subtype and 200 µm total dendritic length.
In addition to the passive parameters (R m = 5000 cm 2 , R i = 100 cm, C m = 1 µF/cm 2 ) the model was equipped with two active conductances for Na + and K + channels (g Na and fast g K as in Zylbertal et al., 2015), and radial and longitudinal diffusion for Na + (as in Zylbertal et al., 2015). The stimulation was implemented as the experimental 50 Hz train of 1000 pA current injections. The basal sodium concentration was set to 15 mM.
Core assumptions of the model were (i) faithful propagation of AP trains into proximal and distal dendrite which is in part supported by Ca 2+ imaging data (Egger et al., 2003, but see Discussion). (ii) g Na and g K models as above, which are unlikely to perfectly capture physiological Na + entry itself; both the precise channel kinetics in general and the detailed Na v /K v subtype composition in GC dendrites are not known, therefore we used the model to estimate current densities rather than densities of fictive channels. (iii) a fixed ratio of g Na /g K throughout all compartments (as in Zylbertal et al., 2015). (iv) linear gradients in g Na density along the dendrite.
The source code for this model is available from ModelDB 1 as model entry 225086.

Statistical Analysis
Statistics were performed with the VassarStats online software 2 . The non-parametric Mann-Whitney test was used for all comparisons between two groups except for paired data [e.g., for ( F/F) 50Hz and ( F/F) 80Hz at the same location, or spines and their adjacent dendrite] where the Wilcoxon pair test was used. For correlations, a linear regression analysis was utilized to determine r-values. A cluster analysis of GC properties was performed using the online tool ClustVis (Metsalu and Vilo, 2015), after normalizing the individual parameter values for each property onto the range 0 -1. Data are presented as mean values of parameters ± standard deviation (SD), unless indicated otherwise.

Calibration and General Approach
To allow for the conversion of fluorescence decrements detected by our TPLSM into absolute changes in intracellular [Na + ], we calibrated SBFI fluorescence signals from cells in olfactory bulb slices following a procedure described earlier (Rose et al., 1999;Meier et al., 2006). Such an in situ calibration is required because the properties of SBFI change in intracellular environments, most likely because of increased viscosity. For example, SBFI fluorescence absorption spectra show a significant blue shift when loaded inside cells (e.g., Langer and Rose, 2009). In addition, its FIGURE 1 | Calibration of SBFI fluorescence. (A) Calibration of the Na + sensitivity of SBFI in GC somata in acute tissue slices of the olfactory bulb. Cells were bolus-loaded with the AM-form of SBFI and afterwards subjected to different calibration solutions containing ionophores, the Na + /K + -ATPase inhibitor ouabain and different Na + concentrations as indicated. Stepwise changes in the extracellular [Na + ] from 0 to 50 mM for 15 min each and back caused stepwise changes in the fluorescence of SBFI. The trace represents an average from 12 cells obtained in one experiment. (B) Relationship between changes in the fluorescence of SBFI and intracellular sodium concentration [Na + ] i , normalized to the fluorescence level in Na + -free saline (0 mM). Shown are mean values ± SEM (n = 54 cells, five slices, three animals). The orange line represents a Michaelis-Menten fit of all the data (v max = 73.5%, K D = 20.2 mM). The blue line represents a linear fit of the data obtained between 15 and 50 mM [Na + ] i . Within this range, a 10% change in fluorescence emission of SBFI corresponds to a change of 22.3 mM in [Na + ] i . K + sensitivity is significantly decreased as compared to in vitro calibrations (e.g., Rose and Ransom, 1996;Meier et al., 2006).
For calibration, cells were loaded by bolus injection with the membrane-permeable form of SBFI (SBFI-AM) into the GC layer and MC layer (MCL) as reported earlier for other brain regions (Langer and Rose, 2009;Karus et al., 2015). Next, equilibration of intra-and extracellular Na + concentration was achieved by perfusing SBFI-loaded slices with calibration solutions containing ionophores (3 µM gramicidin D and 10 µM monensin) and the Na + /K + -ATPase blocker ouabain (100 µM).
Slices were then perfused with calibration salines containing different Na + concentrations and resulting changes in SBFI fluorescence were recorded (n = 54 cells; five slices, three animals; cf. Figures 1A,B). As reported earlier (Rose et al., 1999;Meier et al., 2006;Mondragao et al., 2016), changes in relative fluorescence levels were highly correlated to changes in [Na + ] i concentrations, with increasing [Na + ] causing a decrease in fluorescence emission of SBFI.
The changes in relative fluorescence (F/F o ) against [Na + ] i could be fit by a Michaelis-Menten relationship (r = 0.98), revealing an apparent K D of 20.2 mM and a maximal change in fluorescence of v max = 73.5%. This observation indicates that SBFI fluorescence starts to saturate at levels well beyond 50 mM [Na + ] i , and is thus well suited to report [Na + ] i changes within the physiological range.
Pipette [Na + ] was 15 mM in our experiments. Under the assumption that this concentration also represented the baseline [Na + ] i of patch-clamped neurons, a linear correlation was extrapolated for [Na + ] i between 15 and 50 mM (r = 0.99, P < 0.001), with −10% F/F corresponding to an increase [Na + ] i of 22.3 mM. Because our study only reports increases in [Na + ] i rising from the assumed baseline of 15 mM, we chose this simple linear correlation (the r-value of which was even slightly higher than that of the Michaelis-Menten plot), to convert changes in SBFI fluorescence to changes in [Na + ] i . Importantly, and in contrast to the situation if a Michaelis-Menten relationship is applied, this linearization renders data conversion independent from the actual baseline levels, whichin active neurons undergoing axonal and dendritic sodium influx -might be higher that the presumed baseline of 15 mM.
Throughout the remainder of this study, neurons were filled with 1 mM of SBFI in the whole-cell configuration and stimulated with somatic current injections to generate either single APs or trains of 20 APs at 50 Hz or 80 Hz (train duration 400 ms and 250 ms, respectively; see Figure 2A).

Mitral Cell Axons and Dendrites
To establish measurements of Na + transients, we first imaged the initial segment of MC axons (Figures 2A,B). Single APs elicited detectable axonal F/F changes in 6 out of 11 MCs tested (responding axons: mean ( F/F) AP = 4.7 ± 1.2%, corresponding to roughly 10 mM [Na + ] i ). In our set of experiments the detection threshold was ∼2% F/F, corresponding to 4 mM [Na + ] i . All axons showed substantial fluorescence changes in response to trains of 50 Hz (n = 25 cells, mean ( F/F) 50 Hz = 15.3 ± 8.1%, corresponding to 34 mM [Na + ] i ) and 80 Hz (n = 24, mean ( F/F) 80 Hz = 17.8 ± 8.6%, corresponding to 40 mM [Na + ] i ), at a mean distance from the soma of 34 ± 14 µm. We chose to image the axon not too close to the MC soma since we observed that more proximal measurements yielded smaller F/F measurements (n = 3 MCs). While peak amplitudes of Na + transients measured at the same axonal location increased significantly at 80 Hz as compared to the 50 Hz stimulation (n = 22, mean ratio 80/50 = 1.25 ± 0.33, P < 0.005, Wilcoxon test, Figure 2B), the transients decayed with indistinguishable half-durations of τ 1/2_50 Hz = 0.61 ± 0.45 s and τ 1/2_80 Hz = 0.60 ± 0.51 s (not shown). We also investigated Na + transients in MC lateral dendrites because of their involvement in the reciprocal microcircuits with GCs (Figures 2A,C). The measurement locations were on average 66 ± 34 µm distal from the soma (n = 12). Single APs did not evoke detectable changes in SBFI fluorescence. Trains of APs induced small, but reliable fluorescence changes in the majority of cells (50 Hz: n = 10 out of 12 cells, mean ( F/F) 50 Hz = 3.2 ± 1.1% or 7 mM [Na + ] i ; 80 Hz: n = 9 out of 10, mean ( F/F) 80 Hz = 4.4 ± 1.7% or 10 mM [Na + ] i ). Again, peak amplitudes significantly increased between 50 and 80 Hz at the same location (n = 8, mean ratio 80/50 = 1.43 ± 0.35, P < 0.01, Figure 2C). The transients decayed with similar half-durations of τ 1/2_50 Hz = 0.91 ± 0.64 s and τ 1/2_80 Hz = 0.88 ± 0.57 s. Because of the small size of transients we did not explore the propagation of signals along the lateral dendrites.

Granule Cell Apical Dendrites
In GC apical dendrites, responses to single APs could not be resolved (tested in five cells), while fluorescence changes upon 50 Hz trains could be detected in almost all locations tested (in 62 out of 65 locations in 37 GCs; mean ( F/F) 50 Hz = 6.9 ± 3.3%, corresponding to 15 mM [Na + ] i ; Figure 3), at an average distance from the soma of 75 ± 42 µm.
It has been reported before that GC Ca 2+ signals in response to single APs increase with distance from soma and show a plateau within the EPL (Egger et al., 2003). GC Na + signals in response to 50 Hz trains showed similar but less pronounced effects, since ( F/F) 50 Hz amplitudes were weakly positively correlated to the distance of the measurement location from the MCL border ( Figure 3C, r = 0.33, P < 0.005; see also below).
Granule cell dendritic Na + transients ( F/F) 50 Hz decayed with a mean half duration of τ 1/2_50 Hz = 1.02 ± 0.59 s (n = 54 locations). Several distal locations showed half durations of ≥3 s (n = 12), which could not be properly determined due to the limited scan time and/or noise (Methods) and thus were not included in the above mean value. Overall, a correlation with distance was observed, with a clear trend for slow transients to occur in the distal parts of the dendrites (Figure 3C, r = 0.43, P < 0.001; see also GC in Figure 4B). The average rise time of ( F/F) 50 Hz was t rise = 0.30 ± 0.12 s (n = 61 locations) and uncorrelated with distance from the soma (r = 0.24, P = 0.06), indicating that later APs within 50 Hz trains did not fail to invade distal parts of the dendrite (data not shown).
While the pooled data across all GCs and locations ( Figure 3C) suggested a rather weak influence of dendritic location on Na + peak amplitude, analysis of distance versus ( F/F) 50 Hz amplitude within individual GCs revealed a more refined picture (Figure 4). Here, we included all GCs with at least two dendritic measurement locations at a distance of ≥ 50 µm between the first and last location.
We observed two different behaviors, with a subset of GCs showing a consistent decrease (n = 6 of 12 cells) and the remainder (n = 6 cells) a consistent increase of peak amplitudes, as shown in Figure 4C. The average linear slopes of each subtype's pattern were similar except for their sign (−0.07 ± 0.04% F/F per µm vs. 0.06 ± 0.04% F/F per µm). Interestingly, all the GCs with the increasing pattern except for one showed smaller Na + transients within the most proximal apical dendrite measurements (with respect to the MCL; < 5% F/F or 11 mM [Na + ] i ) whereas all GCs with the decreasing pattern showed larger proximal Na + transients (> 5% F/F; P < 0.025; Figure 4C).
Since absolute concentration changes are tightly linked to dendritic SVRs and thus to dendritic morphology, we explored possible correlations between morphology and the observed Na + transient patterns. Because the morphologies of individual GCs could not always be fully recovered at the end of experiments, our analysis was restricted to the parameters that were available for all or most GCs. While we could not detect any significant link to branching itself, i.e., a systematic drop or increase in ( F/F) 50Hz for measurement locations in front of versus beyond individual branch points, there were correlations between the slope of F/F per µm and the distance of the first branchpoint from the MCL and also between the slope and the diameter of the proximal dendrite ( Figure 4D, r = 0.60, P < 0.025 and r = 0.89, P < 0.005). The location of the first branchpoint is indicative of the anatomical subtype of GC: small, compact GCs like the one in Figure 4A with early branching preferentially innervate the deep EPL, whereas GCs with dendrites innervating mostly higher levels of the EPL also usually show more distal first branch points (Mori et al., 1983;Orona et al., 1983; see Discussion). F/F half duration τ 1/2 . Lines: Linear fits of the data. Weak positive correlations were observed for both the amplitude and the half duration. * * * P < 0.005.
GCs with negative slope had their first branchpoint on average 20 ± 22 µm above the MCL, whereas the first branchpoints of GCs with positive slope were located at 70 ± 45 µm (P < 0.025). The open black square denotes the GC that was not clustered with the other cells with positive slope (see Results). (D) Top panel: Plot of the cell-wise slope of ( F/F) 50Hz /distance from MCL versus the distance of the respective GC's first branch point from MCL. Bottom panel: Plot of the cell-wise slope of ( F/F) 50Hz /distance from MCL versus the diameter of the proximal dendritic segment. Again, the open black squares denote the GC that was not clustered with the other cells with positive slope. (E) Simulation of [Na + ] i (top panels) and the associated Na v current density I Na (bottom panels) in a simple compartmental model of the apical dendrite, with an initial diameter of 2.0 µm for the decreasing and of 2.7 µm for the increasing cell type, both tapering linearly to 1.2 µm over 200 µm. I Na was adjusted to match the average experimental [Na + ] i for the respective cell type from panel C above (see Materials and Methods). * P < 0.05; * * * P < 0.005.
A cluster analysis of GCs with respect to the four parameters (slope of transients, amplitude of most proximal F/F measurement, distance of first branch point and diameter of proximal dendrite, see Materials and Methods) indeed resulted in a clustering into two groups according to slope sign, except for the increasing cell with the high initial F/F marked in the left panel of Figure 4C. However, since this data point is already rather remote from the MCL, a more proximal measurement in the same cell might have been of considerably smaller amplitude.
To provide a qualitative explanation of the observed GC subtype behavior, we used a simple compartmental model (see Materials and Methods, Figure 4E) in which we varied the Na + current density along the dendrite and matched its values to the observed [Na + ] i . Based on the key assumption that AP trains faithfully propagate into the dendrite, these simulations yielded a negative gradient for the decreasing pattern that would be about two times stronger than the positive gradient for the increasing pattern. This difference as well as the plateau in the proximal part of the decreasing cell pattern can be explained by the effect of tapering, which by itself creates an increase in [Na + ] i with distance because of the increased SVR. The model was robust with respect to branching and effects of Na + diffusion, both in line with our data.

Granule Cell Large Spines and Adjacent Dendritic Shafts
The large GC spines on the apical dendrite are known to contain reciprocal synapses and to exhibit voluminous spine heads and long spine necks (Woolf et al., 1991). Figure 5A shows representative responses to 50 Hz trains within putative reciprocal GC spines and their adjacent dendritic shafts. The 16 spines in our sample were located on average 98 ± 31 µm away from the soma and responded with a mean ( F/F) 50 Hz = 7.1 ± 2.7% corresponding to 16 mM [Na + ] i . Again, the decay of GC spine Na + transients was mostly on the order of 1 s with a mean half duration τ 1/2_50 Hz = 1.10 ± 0.76 s (n = 12 spines) with a subset of four spines showing half durations of ≥ 3 s that could not be precisely measured as described above for the GC dendrite.
Compared to the adjacent dendrite, there were no consistent differences between the transients of a given spine and its dendrite with respect to amplitude, half duration and rise time, as shown in Figures 5A,B (ratio S/D F/F 50 Hz = 1.02 ± 0.24, n = 16, P = 0.76; ratio S/D τ 1/2_50 Hz = 1.29 ± 1.27, n = 15, P = 0.50; ratio S/D t rise_50 Hz = 1.09 ± 0.44, n = 14, P = 0.91). Spine transients as a whole set also showed no correlation with distance from the MCL, for both amplitudes ( Figure 5C) and half durations (not shown).

DISCUSSION
Technical Aspects: Suitability of Na + Imaging to Explore MC and GC Function In this study, we have established two-photon Na + imaging with the fluorescent indicator dye SBFI to explore propagation of action potentials based on Na + influx through voltagegated channels in mitral and GC neurites. Compared to commonly used chemical calcium indicators such as Fura-2 or Oregon-Green-Bapta-1, SBFI has a low quantum efficiency and exerts a high K D (Schreiner and Rose, 2012). Moreover, both relative changes in [Na + ] i over baseline and relative changes FIGURE 5 | Na + signals in large GC spines and their adjacent dendrites. (A) Two-photon scans of four representative large GC spines with averaged fluorescence transients in response to 50 Hz trains imaged at the spine heads (S, red traces) and the adjacent dendrites (D, black traces). Train durations are indicated by the gray rectangles. (B) Cumulative data for all spine-dendrite pairs for ( F/F) 50Hz amplitudes and half durations. (n = 16 pairs). Open symbols: individual data, solid black symbols: averaged data ± SD. Data from the same spine-dendrite pair are interconnected. Because F/F in more than one location showed half durations ≥3 s, the respective data are accumulated in the 3 s data points and more than one pair can be connected to these points. (C) Cumulative data for all spine-dendrite pairs for ( F/F) 50Hz amplitudes versus their distance from the MCL border. There is no significant trend, and the linear fits for spine and dendrite data each (red and black dotted line) are highly similar.
in dye emission are rather small during physiological activity, necessitating the use of high dye concentrations to obtain an acceptable signal-to-noise-ratio (Rose et al., 1999). Importantly, and again in contrast to high-affinity Ca 2+ indicator dyes (Neher and Augustine, 1992;Maravall et al., 2000), SBFI does not buffer Na + and, therefore, does not distort Na + signals even at a dye concentration of 1 mM as used in the present study (Sabatini et al., 2001;Mondragao et al., 2016).
Similar to earlier reports (e.g., Mondragao et al., 2016), our in situ calibrations revealed an acceptable linear correlation between changes in SBFI emission ( F/F) and changes in [Na + ] i between 15 and 40-50 mM [Na + ] i . This enabled a calculation of absolute changes in [Na + ] i without the necessity of knowing baseline [Na + ] i . Such quantitative measurements with two-photon microscopes are far more tedious for Ca 2+ concentration changes, where basal and maximal fluorescence has to be established in every cell of interest (e.g., Maravall et al., 2000;Egger and Stroh, 2009). Taken together, these considerations illustrate that SBFI is an indicator well-suited for quantitative measurement of activity-induced [Na + ] i transients.
Large [Na + ] i transients were found in MC axons following somatic current injection and induction of APs. Even single APs evoked detectable [Na + ] i in half of the axons, amounting to an average increase by 10 mM, while trains of APs resulted in [Na + ] i transients by more than 30 mM. Activity-related [Na + ] i signals in this range have also been reported from dendrites and spines of CA1 pyramidal neurons in response to suprathreshold synaptic stimulation and opening of ionotropic glutamate receptors (Rose and Konnerth, 2001). Interestingly, peak amplitudes were largest at a distance of ∼35 µm from the soma. While we have not investigated this observation in greater detail, it indicates that the spike initiation zone is not directly located at the cell body but at some distance from it, as demonstrated before based on Na + imaging or patch recordings from axon initial segments, e.g., in cortical pyramidal neurons (Kole et al., 2008;Hu et al., 2009;Fleidervish et al., 2010), cochlear interneurons (Bender and Trussell, 2009) and dentate gyrus GCs (Scott et al., 2014). Since in a few MCs the train-evoked transients reached amplitudes beyond the linear regime of calibration, we would recommend to use shorter AP trains for future studies of MC axonal Na + signaling.
Action potential-induced axonal Na + transients rapidly decayed with half-durations of ∼0.6 s. This is in the range of decay times determined for layer five pyramidal neurons (∼0.5 s; Fleidervish et al., 2010). The latter study also provided evidence that the decay of axonal Na + transients is largely governed by fast lateral diffusion. The same phenomenon was described for local, glutamate-induced Na + transients in dendrites of CA1 pyramidal neurons, where decay to baseline was dominated by diffusion, resulting in a more than 10-fold higher rate of recovery as compared to extrusion through the Na + /K + -ATPase (Mondragao et al., 2016).
Na + imaging also enabled detection of Na + APs backpropagating into MC lateral and apical dendrites. In both dendrite types, Na + influx in response to single APs could not be resolved, while trains of 20 APs resulted in small signals as observed earlier (Rose et al., 1999). However, on our TPLSM signals from both compartments were close to the detection limit, especially in the apical dendrite. Thus we would like to suggest that in MCs SBFI-based Na + imaging will mainly be helpful to explore axonal function or to study activity in the highly excitable apical dendritic tuft (as performed by Zylbertal et al., 2015; not investigated here), whereas it is less suitable for studies of propagation into MC dendrites, which might be better studied via dendritic patching (Chen et al., 1997;Christie and Westbrook, 2003).
In GC apical dendrites we observed substantially larger Na + transients as compared to MC lateral dendrites in response to 50 Hz trains which allowed us to quantitatively compare signals along the dendrite of individual GCs and also across spines and their adjacent dendrites. In general, the observed GC ( F/F) 50Hz signals were within the linear regime of our calibration (range 20-45 mM including baseline [Na] i = 15 mM, Figure 1B).

GC Apical Dendrite Na + Signals vs. Ca 2+ Signals
In an earlier study, we have observed that single sAPs result in Ca 2+ entry throughout GC apical dendrites (Egger et al., 2003). These sAP-mediated Ca 2+ signals did not decrease with distance from the soma but rather showed a plateau within the EPL, which we interpreted as a property tied to dendritic transmitter release onto MCs and TCs within the EPL. Here, we found that sAP-train-mediated Na + signals did either rise or decline with distance from the MCL, in a manner that was correlated with morphological features of the respective cells.
Since such differential patterns might have been obscured in the initial Ca 2+ imaging study due to averaging across cells, we have now performed a cell-wise reanalysis of the Ca 2+ imaging data set with respect to morphology and slope in ( F/F) sAP but could not reveal any significant correlations between Ca 2+ transient amplitude changes along the apical GC dendrite and branchpoint location or proximal diameter (n = 19 cells; only 4 of these showed negative slopes in ( F/F) sAP in the first place). If we assume that there is little failure of invasion by late APs in the train (as indicated by the lack of a negative correlation between the rise time of ( F/F) 50Hz and the distance from soma in the Na + imaging data set, see results), the difference between Na + and Ca 2+ signaling might be explained by the fact that a substantial fraction of Ca 2+ is carried by low-voltage activated Ca 2+ channels (Egger et al., 2003;Pinato and Midtgaard, 2005) and thus a slight attenuation or increase in Na + entry with distance might have only a limited effect on the sAP-associated Ca 2+ signals. Moreover, if the observed Na + signaling patterns can indeed be explained by variations in Na v current density, Ca v current density in turn is unlikely to systematically covary with Na v current density, which also has not been observed in other neuron types (e.g., Magee, 2008). Also, if Na + imaging would simply mirror GC dendritic excitability with respect to Ca 2+ entry, Na + imaging in GCs would be of little extra use compared to Ca 2+ imaging since the Na + -dye properties do not allow to resolve single APs so far.
To shed light on local correlations between Na + and Ca 2+ entry, simultaneous dendritic Na + and Ca 2+ imaging would be required, as performed, e.g., by Jaffe et al. (1992), who observed in CA1 pyramidal cell dendrites sAP-train-mediated Na + and Ca 2+ signals at distances at least up to 200 µm from the soma. Moreover, they found that distal Na + transients were smaller than proximal signals and that their rise was restricted to the initial phase of sAP trains, probably due to an invasion failure of later APs. Ca 2+ transients in the same cells were found to mirror this behavior, leading the authors to the conclusion that AP invasion was required for Ca 2+ influx to occur. In our studies in GCs there was no such indication of propagation failures during trains, probably due to the high excitability of the GC axo-dendrite.

Functional Aspects: GC Morphology and Dendritic Excitability
The distance of the first branchpoint of the apical GC dendrite from the MCL is indicative of the morphological subtype of GC, as can be gathered from the GC anatomies shown in Mori et al. (1983) and Orona et al. (1983): GCs with early branching preferentially innervate the deep EPL or possibly the entire EPL, while GCs with more distal first branch points innervate mostly higher levels of the EPL. Mori et al. differentiated the subtypes I, II, III depending on the innervated sublayer(s) of the EPL. While our sample of GCs is biased with respect to the distance of the soma from the MCL (all closer than 100 µm), all three subtypes might have been picked up, since the somata of subtype I and II can be found throughout the GCL and subtype III occurs preferentially close to or within the MCL (Mori et al., 1983). In addition, we also might have sampled type V shrub GCs (Merkle et al., 2014;Nagayama et al., 2014). The GCs with positive slope in F/F would thus correspond mostly to subtype III, projecting into the superficial EPL and thus interacting preferably with TCs. On the other hand, the GCs with decreasing F/F would belong to subtype II and perhaps also I, interacting rather with MCs. With respect to the GCs in our data set, an unequivocal assignment to these subtypes would require both a complete recovery of GC morphology, including their extent relative to the EPL width, and a more quantitative description of GC morphotypes in the literature.
Nevertheless, our findings for the first time reveal functional differences in active dendritic properties of GCs innervating different sublayers of the EPL. Together with the observation that deep-branching cells start out with larger Na + signals compared to the higher branching cells, these patterns might relate to different Na v channel densities and gradients along the dendrites as implied by our simulations of Na + current densities ( Figure 4E). Our model also predicts that the increasing pattern is not caused by decreasing dendritic diameters alone (due to the concomitant increase in SVR along the dendrite) but could involve a roughly twofold increase in current density along the dendrite. Similarly diverse types of Na v channel or Na v current distributions have been reported from several other cell types, with a negative gradient of Na v 1.6 expression in CA1 pyramidal cells (Lorincz and Nusser, 2010), a constant high density in a subtype of hippocampal interneuron (Martina et al., 2000) and an increasing Na + current density in CA3 pyramidal cells (Kim et al., 2012).
In GCs the observed differential patterns might serve to foster AP propagation for GC subtypes with larger dendritic trees or at least might compensate for an increased likelihood of propagation failures because of the extended morphology, while the higher current density in the less wide proximal segment of the smaller cells might support dendrosomatic AP conduction.
If the different subclasses of GCs indeed dispose of different Na v channel densities and different dendritic gradients, these differences could also impact synaptic integration and thus GC subtype interaction with the respective mitral/TC subnetworks via their dendrodendritic output. Recent studies indicate substantial functional differences between MCs and TCs with respect to temporal processing of olfactory stimuli: TCs are firing earlier in the respiration cycle and are thought to encode the detection of the presence of an odor, whereas MCs are firing later and might contribute to difficult odor discriminations (Fukunaga et al., 2012;Igarashi et al., 2012). Moreover, TCs are involved in early onset fast gamma oscillations while MCs contribute to late onset slow gamma oscillations (Manabe and Mori, 2013). Thus the dendritic excitability of their GC partners in the respective subnetworks could be specifically tuned toward these temporal processing requirements. Since increasing cells are preferentially interacting with TCs and are predicted to dispose of a higher current density in their distal dendrites (than the decreasing cells that interact with MCs), it is tempting to speculate that increasing GCs can generate APs (or Na + spikelets, see Introduction) in the distal dendrite and thus lateral output within the GC-TC network more rapidly than decreasing GCs, both in line with a contribution to fast gamma and to fast TC output. Decreasing GCs on the other hand can afford to provide lateral inhibition on a slower time scale, potentially serving to decorrelate MC activity during difficult odor discriminations (Friedrich and Wiechert, 2014).

GC Spine Na + Signaling
As observed earlier for other cell types such as CA1 pyramidal neurons (Rose et al., 1999), trains of backpropagating APs not only induced Na + signals in GC dendrites, but also in their large spines. GC spine signals did not differ from dendritic Na + signals in respect to their peak amplitude, rise or decay times. These data, therefore, do not allow to draw a clear conclusion about the presence or absence of voltage-gated Na v in spine heads. A recent study in CA1 pyramidal neurons showed that diffusion of Na + ions along dendrites is considerably slower than previously thought (320 µm 2 /s; Mondragao et al., 2016). This value, is, however, still about 20-fold faster than that of buffered diffusion of Ca 2+ ions (∼15 µm 2 /s; Yuste et al., 2000), indicating that any differences in Na + concentrations between spines and dendrites upon differential Na + influx will be equilibrated rapidly.
On the other hand, there is substantial indirect evidence for the presence of Na v in GC spine heads from Ca 2+ imaging (Bywalez et al., 2015) where quasi-synaptic stimulation of GC spines via uncaging of glutamate resulted in a TTX-sensitive Ca 2+ signal that can be explained by the existence of a local postsynaptic Na + spike within the spine head. GC Ca 2+ signals due to propagation of sAPs or due to synaptically evoked global APs were also found to be very similar in spines and their adjacent dendrites (Egger et al., 2003;Egger, 2008). Moreover, both Ca 2+ handling within and the SVRs of these two compartments were found to be highly similar (Egger and Stroh, 2009). Taken together, these findings argue for a similar density of Na v channels and low-voltage activated Ca 2+ channels within GC spines and their adjacent dendrites. To fully validate this conclusion the precise Na + extrusion/diffusion mechanisms in GC dendrites need to be determined.

CONCLUSION
In summary, we have shown that Na + imaging can reveal distinct functional aspects of the active GC dendrites that have not been accessible to Ca 2+ imaging. In the future, we expect that such Na + imaging data will allow us to further constrain simulations of dendritic integration in GCs or of post-synaptic GC spine signals (as in Bywalez et al., 2015), and yield insights into the role of Na + currents in reciprocal processing at the dendrodendritic MC-GC or TC-GC synapse.

AUTHOR CONTRIBUTIONS
All authors take responsibility for the integrity of the data and the accuracy of the data analysis. Study concept and design: CR and VE. Calibration: NG and TO-J. Acquisition of other data: TO-J. Simulations: SS. Analysis and interpretation of data: all.

FUNDING
The instrumentation used in this work was funded by the BMBF (FKZ 01GQ1502; VE) with additional equipment funding by DFG-SFB 870 (EG135/5-1) and GSN-LMU. NG and CR were supported by the DFG (Ro2327/6-1). SS was supported by a fellowship from the "I.R. Iran Cognitive Sciences and Technologies Council."