- 1Interdepartmental Program in Neuroscience, University of Utah, Salt Lake City, UT, United States
- 2 Department of Neurobiology, University of Utah School of Medicine, Salt Lake City, UT, United States
In the mammalian olfactory bulb (OB), mitral/tufted (MT) cells respond to odorant inhalation with diverse temporal patterns that are thought to encode odor information. Much of this diversity is already apparent at the level of glutamatergic input to MT cells, which receive direct, monosynaptic excitatory input from olfactory sensory neurons (OSNs) as well as a multisynaptic excitatory drive via glutamatergic interneurons. Both pathways are also subject to modulation by inhibitory circuits in the glomerular layer of the OB. To understand the role of direct OSN input vs. postsynaptic OB circuit mechanisms in shaping diverse dynamics of glutamatergic drive to MT cells, we imaged glutamate signaling onto MT cell dendrites in anesthetized mice while blocking multisynaptic excitatory drive with ionotropic glutamate receptor antagonists and blocking presynaptic modulation of glutamate release from OSNs with GABAB receptor antagonists. GABAB receptor blockade increased the magnitude of inhalation-linked glutamate transients onto MT cell apical dendrites without altering their inhalation-linked dynamics, confirming that presynaptic inhibition impacts the gain of OSN inputs to the OB. Surprisingly, blockade of multisynaptic excitation only modestly impacted glutamatergic input to MT cells, causing a slight reduction in the amplitude of inhalation-linked glutamate transients in response to low odorant concentrations and no change in the dynamics of each transient. The postsynaptic blockade also modestly impacted glutamate dynamics over a slower timescale, mainly by reducing adaptation of the glutamate response across multiple inhalations of odorant. These results suggest that direct glutamatergic input from OSNs provides the bulk of excitatory drive to MT cells, and that diversity in the dynamics of this input may be a primary determinant of the temporal diversity in MT cell responses that underlies odor representations at this stage.
Introduction
In the mammalian olfactory system, the neural representation of olfactory information is inherently dynamic, with respiration and active odor sampling (i.e., sniffing) driving inhalation-linked bursts of activity in olfactory sensory neurons (OSNs) that are passed on to higher-order neurons, including mitral and tufted (MT) cells, the main output neurons of the mammalian olfactory bulb (OB; Schaefer and Margrie, 2007; Wachowiak, 2011). At the same time, OSNs and MT cells show slower temporal patterning across respiration cycles (Patterson et al., 2013; Eiting and Wachowiak, 2020). Temporal patterning at each of these timescales is hypothesized to play important roles in encoding odor information at the level of the OB and in the piriform cortex, the main target of MT cell projections (Schaefer and Margrie, 2007; Stern et al., 2018; Chong et al., 2020). The cellular and circuit mechanisms underlying temporal patterning at the level of MT cell output from the OB remain unclear.
Potential contributions towards the temporal patterning of MT cell activity include diverse patterns of OSN activation (Spors et al., 2006; Short and Wachowiak, 2019), presynaptic inhibition of glutamate release from OSN terminals (McGann et al., 2005; Wachowiak et al., 2005), or neural circuits of the OB. Potential OB circuit mechanisms for shaping MT cell output dynamics have been well-characterized using OB slice preparations, and include modulation of glutamate release from ET cells onto MT cells (Hayar et al., 2004a; De Saint Jan et al., 2009; Gire et al., 2012) as well as feedforward inhibition via several distinct inhibitory interneuron pathways (Murphy et al., 2005; Gire and Schoppa, 2009; Shao et al., 2012). In a recent study using glutamate and Ca2+ imaging from the mouse OB in vivo (Moran et al., 2021), we found that glutamatergic signaling onto MT cells showed complex dynamics both within and across inhalations and that the slower temporal patterning was well-correlated with that of MT cell Ca2+ signals, highlighting the importance of excitatory pathways in generating MT cell patterning.
Here, we sought to dissect the contribution of OSNs vs. second-order OB circuits to generating dynamic glutamatergic signaling onto MT cells. We used in vivo pharmacology to examine the contributions of these different sources to the dynamics of MT cell glutamatergic input in anesthetized mice. We directly imaged odorant-evoked glutamate signaling onto MT cells while pharmacologically blocking multisynaptic excitatory drive with ionotropic glutamate receptor antagonists and blocking presynaptic modulation of glutamate release from OSNs with GABAB receptor antagonists. Neither manipulation substantially impacted the inhalation-linked temporal dynamics of glutamate signaling onto MT cells. Furthermore, blocking multisynaptic excitation only weakly reduced the magnitude of glutamatergic excitation, and modestly impacted glutamate dynamics over a slower timescale spanning multiple inhalations of odorant. Overall, these results suggest that direct glutamatergic input from OSNs provides the bulk of excitatory drive to MT cells and that this direct OSN–MT cell pathway may be the primary determinant of inhalation-linked temporal patterning of MT cell activity.
Materials and Methods
Animals
Experiments were performed on male and female mice expressing Cre recombinase (Cre) in defined neural populations. Mouse strains used were: Pcdh21-Cre [Tg(Pcdh21-cre)BYoko], Gensat Stock #030952-UCD; OMP-Cre [Tg(Omp-tm4-Cre)Mom], JAX Stock #006668, Tbet-Cre [Tg(Tbx21-cre)1Dlc], JAX Stock #024507, CCK-IRES-Cre [Tg(CCK-IRES-Cre)Zjh], JAX Stock #012706 (Haddad et al., 2013), and Thy1-jRGeCO1a Tg(Thy1-jRGECO1a)GP8.31Dkim/J, JAX Stock #030526 (Dana et al., 2018). Mice ranged from 3 to 8 months in age. Mice were housed up to four/cage and kept on a 12/12 h light/dark cycle with food and water available ad libitum. All procedures were carried out following the National Institutes of Health Guide for the Care and Use of Laboratory Animals and were approved by the University of Utah Institutional Animal Care and Use Committee.
Viral Vector Expression
Viral vectors were obtained from the University of Pennsylvania Vector Core (AAV1 or 5 serotype, AAV.hSynap-FLEX.iGluSnFR, AAV1 serotype AAV.hSynap-FLEX.GCaMP6f), HHMI Janelia Campus or Vigene (AAV1 or 5 serotype, pAAV.hSynap-FLEX.SF-iGluSnFR.A184V, pAAV.hSynap-FLEX.SF-iGluSnFR.A184S). Virus injection was done using pressure injections and beveled glass pipettes, as described previously (Rothermel et al., 2013; Wachowiak et al., 2013; Short and Wachowiak, 2019). After injection, mice were given carprofen (Rimadyl, S.C., 5 mg/kg; Pfizer) as an analgesic and enrofloxacin (Baytril, I.M., 3 mg/kg; Bayer) as an antibiotic immediately before and 24 h after surgery. Mice were singly housed after surgery on ventilated racks and used 21–35 days after virus injection. In some mice, viral expression was characterized with post hoc histology using native fluorescence.
In vivo Two Photon Imaging
Two-photon imaging in anesthetized mice was performed as described previously (Wachowiak et al., 2013; Economo et al., 2016). Mice were initially anesthetized with pentobarbital (50–90 mg/kg) then maintained under isoflurane (0.5–1% in O2) for data collection. Body temperature and heart rate were maintained at 37°C and ~400 beats per minute. Mice were double tracheotomized and isoflurane was delivered passively via the tracheotomy tube without contaminating the nasal cavity (Eiting and Wachowiak, 2018). Two-photon imaging occurred after removal of the bone overlying the dorsal olfactory bulb.
Imaging was carried out with a two-photon microscope (Sutter Instruments or Neurolabware) coupled to a pulsed Ti:Sapphire laser (Mai Tai HP, Spectra-Physics; or Chameleon Ultra, Coherent) at 920–940 nm and controlled by either Scanimage (Vidrio) or Scanbox (Neurolabware) software. Imaging was performed through a 16×, 0.8 N.A. objective (Nikon) and emitted light detected with GaAsP photomultiplier tubes (Hamamatsu). Fluorescence images were acquired using unidirectional resonance scanning at 15.2 or 15.5 Hz. For dual-color imaging, a second laser (Fidelity-2; Coherent) was utilized to optimally excite jRGECO1a (at 1,070 nm) and emitted red fluorescence collected with a second PMT, as described previously (Short and Wachowiak, 2019; Moran et al., 2021).
In vivo Pharmacology
in vivo pharmacology was carried out after removing the bone and dura overlying the dorsal olfactory bulb, using protocols described previously (Pírez and Wachowiak, 2008; Brunert et al., 2016). For experiments where both OBs were imaged sequentially in a single mouse, we performed pharmacology on one bulb while the other bulb remained encased in bone. After the experiment was complete on the initial OB, the bone overlying the contralateral OB was repeatedly rinsed with Ringers solution, thinned then craniotomized, the underlying dura removed, and the experiment repeated. Drug solutions (1 mM CGP35348 and 0.5 mM NBQX/1 mM APV) were dissolved in Ringers solution, pre-warmed on a heating block, and applied in bulk to the dorsal OB without the use of agarose or coverslip. For higher concentration experiments, we utilized 2.5 mM NBQX/5 mM APV dissolved in Ringers immediately prior to use. We waited at least 10 min after drug application to allow for absorption and temperature equilibration across the tissue. For the experiments where we applied ionotropic glutamate blockers following GABAB blockade, NBQX + APV was applied immediately after CGP35348 without an intervening control wash. Due to the extremely slow washout of CGP35348 in vivo, we considered CGP35348 to still be present during the subsequent NBQX + APV application.
Odorant Stimulation
In most experiments, odorants were presented as precise dilutions from saturated vapor (s.v.) in clean, desiccated air using a custom olfactometer under computer control, as described previously (Bozza et al., 2004; Economo et al., 2016). Odorants were presented for durations ranging from 2 to 8 s for most experiments, with inhalation rates ranging from 1 to 2 Hz, as specified; for experiments used to obtain inhalation-triggered averages (ITAs), odorant duration was 70 s with 0.25 Hz inhalations. Inter-trial intervals were 20 s for 2-s odorant presentations and 36 s for 4- or 8-s presentations. Clean air was passed across the nostrils in between trials to avoid contribution from extraneous odorants in the environment. Odorants were pre-diluted in solvent (1:10 or 1:25 in mineral oil or medium chain triglyceride oil) to allow for lower final concentrations and then diluted to concentrations ranging from 1% to 3% s.v. Relative increases in concentration were confirmed with a photoionization detector (miniRAE Lite, PGM-7300, RAE Systems) 3 cm away from the flow dilution port. Estimated final concentrations of odorants used ranged from 0.03 to 48 ppm, depending on vapor pressure and s.v. dilution (Table 1). For experiments testing a larger panel of 12 odorants (e.g., Figure 6), we used a novel olfactometer design that allowed for rapid switching between odorants with minimal cross-contamination (Burton et al., 2019). Here, odorants were presented for 3 s, in random order, using 10 s interstimulus intervals (Table 1) and repeated three times each. The odorant presentation was as described in a previous publication (Burton et al., 2019), using an eductor nozzle for additional mixing in a carrier stream of filtered air. The end of the eductor was placed 5–7 cm from the nose. With the configuration used, estimated dilutions of odorant were approximately 1.5% s.v.; odorants were prediluted to achieve relatively sparse activation of dorsal glomeruli (Burton et al., 2019). Estimated final concentrations ranged from 0.03 to 15 ppm.
Figure 1. GABAB-mediated presynaptic inhibition regulates the strength, but not the timing, of inhalation-driven glutamatergic signaling onto mitral/tufted cells. (A) Baseline fluorescence (mean F) and ΔF maps (OMP-Cre × Rosa26-GCaMP6f) of calcium responses to methyl valerate before (Predrug) and after application of CGP35348. (B) Traces showing mean fluorescence of all glomerular regions of interest (ROIs) in one olfactory bulb (OB) to the odorant presentation (ethyl butyrate) before and after application of CGP35348. Trace showing artificial inhalation pressure is below the response traces (2 Hz). Lower traces show expansion of responses to the first and second inhalation of odorant scaled to the same maximum. (C) Peak OSN—GCaMP6f responses (ΔF/F) in each measured ROI, normalized to the maximal predrug response in each experiment. Each color is an individual experiment. Summary statistics for each experiment (median CGP35348/Predrug ratio, Wilcoxon Signed Ranks Test, p-value corrected for multiple comparisons): Exp. 1: 1.60, Z = −2.89, p = 0.004, Exp. 2: 1.91, Z = −3.60, p = 3.21 × 10−4, Exp 3: 1.76, Z = −3.70, p = 2.14 × 10−4, Exp. 4: 2.13, Z = −3.80, p = 1.43 × 10−4, Exp. 5: 2.46, Z = −3.02, p = 0.003. (D) Odorant-evoked ΔF maps showing glutamate signals (SF-iGluSnFR.A184V) on MT cells before and after CGP35348 application. (E) Mean iGluSnFR fluorescence across all glomerular regions of interest (ROIs) in one olfactory bulb (OB) to the odorant presentation (ethyl butyrate) before and after CGP35348 application. Odorant was delivered using 1 Hz artificial inhalation. Bottom traces show expanded responses to the first two inhalations of odorant, scaled to the peak of the first inhalation response. (F) Peak MT Cell—iGluSnFR ΔF/F responses (ΔF/F) in each measured ROI before and after CGP35348 application, normalized to the maximal Predrug response in each experiment. Summary statistics, reported as in (C): CGP35348/Predrug (first inhal.; mean ± SD): Exp. 1: 1.2 ± 0.06, n = 19, p = 2.9 × 10−10; Exp. 2: 0.97 ± 0.17, p = 1.0, n = 20; Exp. 3: 1.24 ± 0.15, n = 16, p = 2.6 × 10−5; Exp 4: 1.4 ± 0.2, n = 14, p = 2.5 × 10−5. CGP35348/Predrug (second inhal): Exp. 1: 1.34 ± 0.11, p = 3.0 × 10−10, Exp 2: 1.15 ± 0.14, p = 3.9 × 10−4, Exp 3: 1.34 ± 0.20, p = 1.7 × 10−5, Exp 4: 1.48 ± 0.13, p = 2.1 × 10−8. Asterisks indicate Experiments with p <0.05 after correction for multiple comparisons). (G) Change in onset and peak latencies of iGluSnFR response to the first inhalation of odorant after application of CGP35348. Same mice and glomeruli as in (E,F). Each point is an ROI. Top: Change in latencies relative to predrug values. Boxes indicate 25th and 75th percentiles, line indicates median, whiskers denote outliers with a coefficient of 1.5. Bottom: Absolute latencies relative to inhalation start, before and after CGP35348. CGP35348 had little to no effect on the onset and peak latencies of response to first inhalation (Wilcoxon signed ranks test, medians and p-values per experiment): Δ onset latencies, Exp. 1: −6.3 ms, p = 0.04; Exp. 2: 3.5 ms, p = 0.2; Exp. 3: −3.1 ms, p = 0.016; Exp. 4: 4.9 ms, p = 0.16. Δ peak latencies, Exp 1: 3.3 ms, p = 0.16; Exp. 2: 8.4 ms, p = 8 × 10−4; Exp. 3: −2 ms, p = 0.32 ; Exp. 4: 2.8 ms, p = 0.52.
Figure 2. GABAB-mediated presynaptic inhibition has a heterogeneous impact on mitral/tufted cell calcium responses. (A) Baseline fluorescence (mean F) and ΔF maps of odorant-evoked GCaMP6f signals in mitral/tufted cells, imaged with two-photon microscopy; z-plane is in the deep glomerular layer/outer external plexiform layer of the dorsal OB. Maps show responses to the same odorant (ethyl butyrate) before and after application of CGP35348, and after subsequent application of APV+NBQX. Pseudocolor scale normalized to Predrug response levels. APV+NBQX strongly reduced excitatory responses and eliminated suppressive responses. (B) Response traces for the three glomeruli indicated in (A), (2 are excited, 1 is suppressed) across each condition. (C) Comparison of mitral/tufted (MT) cell GCaMP6f responses for Predrug vs. CGP35348 conditions (12 mice; left plots) and CGP35348 and subsequent APV+NBQX (four mice; right plots) conditions. Colors indicate individual mice. Top row: peak response amplitudes to the first inhalation of odorant, normalized to the maximum Predrug response in each mouse. Bottom row: mean response across odorant presentation. See Text for CGP35348 summary statistics across mice. For APV+NBQX/CGP35348 statistics, ratios per experiment, reported as medians [quartile range] were: Exp. 1: 0.47 [0.37–0.9], n = 6 ROIs; Exp 2: 0.39 [0.25–0.52], n = 11; Exp 3: 0.38 [0.19–0.57], n = 5; Exp 4: 0.19 [0.08–0.27], n = 6; Full duration: Exp 1: 0.66 [0.46–1.12], n = 6; Exp. 2: 0.4 [0.31–0.51]; Exp 3: 0.37 [0.27–0.50], n = 4; Exp 4: 0.12 [0.05–0.18], n = 7. ROIs were pooled across the four mice for statistical comparison due to the small number of ROIs in each experiment.
Figure 3. Near-completeblockade of odorant-evoked excitation of CCK+ superficial tuftedcells by APV+NBQX. (A) Odorant response maps of GCaMP6fsignals imaged from CCK+ superficial tufted cells, taken acrossPredrug, CGP35348, and APV+NBQX conditions, scaled to Predrugresponse levels. (B) Response traces for the two glomeruliindicated in (A), showing enhanced response amplitudes afterapplication of CGP35348 and either elimination (top) or sharpreduction (bottom) in responses after subsequent application ofAPV+NBQX. Note that in ROI 1, the response after APV+NBQX iseliminated whereas, in ROI 2, the remaining response is present butdelayed. Note also the two-fold difference in scale for ROI2 due to differences in overall response magnitude. (C)Example of high-zoom imaging of GCaMP6f signals imaged fromsuperficial CCK+ neurons innervating a glomerulus, with response mapsin Predrug, CGP35348, and APV+NBQX conditions. (D) Tracesshowing odorant-evoked response from glomerular neuropil in(C). CGP35348 increases response amplitudes. APV+NBQXdecreases response amplitude, with delay in remaining signal relativeto inhalation. Lower traces show the expansion of response to thefirst two inhalations of odorant, highlighting delayed GCaMP6fresponse after APV+NBQX. (E) Comparison of CCK+ GCaMP6fresponses across CGP35348 and subsequent APV+NBQX conditions, imagedfrom 14 glomeruli in three OBs (two mice). Scatter plots show peakresponse to first inhalation of odorant in Predrug compared toCGP35348 condition (left) or CGP35348 and APV+NBQX conditions(right). Shaded regions in APV+NBQX vs. CGP35348 plot (upper right)indicate two groups of glomeruli defined by the magnitude of responsereduction by APV+NBQX (see Text). Group 1 defines glomeruli whoseresponses persisted following APV+NBQX; group 2 defines glomeruli whose responses were eliminated. (F) Comparison of CCK+ GCaMP6f responses before (Predrug) and after higher concentrations of APV+NBQX (5 mM/2.5 mM), and without prior application of CGP35348. Separate experiments from (E). Responses are fully blocked in 11 of 13 glomeruli, (three experiments). (G) Change in onset and peak latencies of CCK+ GCaMP6f response to the first inhalation of odorant after application of CGP35348 and APV+NBQX, for the six group one glomeruli shown in (E). Response latencies for group 1 glomeruli are unchanged by CGP35348 application (left), but increase significantly after APV+NBQX application (right). Summary statistics, mean ± SD, paired t-test: CGP35348 vs. Predrug Δ onset latency: −0.5 ± 12 ms, p = 0.91; Δ peak latency: 19 ± 45 ms, p = 0.34; APV+NBQX vs. CGP35348 Δ onset latency: 255 ± 79 ms, p = 5 × 10−4; Δ peak latency: 257 ± 51 ms, p = 6 × 10−5.
Figure 4. Blockade of postsynaptic excitation does not impact the magnitude or latency of inhalation-linked glutamate signaling onto piriform-projecting mitral/tufted (pcMT) cells. (A) Odorant-evoked SF-iGluSnFR.A184V response maps across Predrug, CGP35348, and APV+NBQX conditions in an aPC-injected Tbet-Cre mouse. Maps are scaled to max of Predrug condition. (B) Top: Traces showing responses of ROI 1 (from A) across conditions. Inhalation frequency, 1 Hz. inset shows the expansion of response to the first inhalation of odorant in each condition, showing a slight increase in latency to peak after APV+NBQX application. Bottom: Example traces from ROI in second preparation (Exp. 2, imaged with epifluorescence), showing similar drug effects, although with some change in tonic (non-inhalation-linked) signal. Inset shows expansion of response to the fourth inhalation, with traces shifted to start at same pre-inhalation baseline. Traces low-pass filtered at 5 Hz for display. (C) Comparison of peak SF-iGluSnFR.A184V responses (1st inhalation) across CGP35348 (left) and APV+NBQX (right) conditions. Points indicate glomeruli from two mice (red and black). for the two mice. Statistics for each of two experiments: CGP35348/Predrug ratios (mean ± SD, t-test compared to ratio of 1): Exp. 1, 1.35 ± 0.097, p = 5 × 10−4, n = 6 glomeruli; Exp. 2, 1.90 ± 0.3, p = 1.6 × 10−3, n = 6 glomeruli; APV+NBQX/CGP35348 ratios: Exp. 1, 1.02 ± 0.07, p = 0.94; Exp. 2, 0.91 ± 0.1, p = 0.2. (D) Little change in onset and peak latencies of SF-iGluSnFR.A184V response to the first inhalation of odorant after application of CGP35348 and APV+NBQX, for same glomeruli and animals as in (C). Top: Change in onset or peak latency relative to pre-drug condition. Asterisks indicate a significant change in latency after drug treatment (see Text for statistics). Bottom: Absolute onset and peak latencies before (CGP) and after APV+NBQX application for each ROI. Latency measurements were pooled across the two experiments due to the lower number of ROIs supporting reliable measurements.
Figure 5. Dual-colorimaging confirms small impact of multisynaptic excitation onglutamatergic input to MT cells. (A) Top: Mean fluorescenceand ΔF response maps showing odorant-evoked (ethyl tiglate,0.7 ppm) glutamate signals in several glomeruli afterSF-iGluSnFR.A184S (A148S) expression in pcMT cells. Bottom:Inhalation-triggered average (ITA) responses from four of theglomeruli, imaged before and after application of APV+NBQX(5 mM/2.5 mM). Response amplitudes in each glomerulus are modestly reduced, but relative latencies and durations are unchanged. Lower traces show ITAs for each of four glomeruli, before and after APV+NBQX application, scaled to the same maximum (gray trace is post-drug). (B) Dual imaging of glutamate and Ca2+ from overlapping subsets of MT cells. Left: post hoc histology showing expression of SF-iGluSnFR.A184S in pcMTs (green) and Thy1-driven jRGeCO1a expression in MT cells (magenta), in a Tbet-Cre × Thy1-jRGECO1a mouse cross. Insets highlight expression in the external plexiform layer (EPL) and mitral cell layer (MCL). Note that SF-iGluSnFR is expressed preferentially in mitral cells, while jRGECO1a is expressed in both mitral cells and tufted cells in the EPL. (C) Left: Mean SF-iGluSnFR.A184S fluorescence and ΔF response maps for SF-iGluSnFR.A184S and jRGECO1a, imaged simultaneously in response to odorant stimulation (vinyl butyrate, 3 ppm). Right: ITA traces for the SF-iGluSnFR.A184S (green) and jRGECO1a (magenta) signals imaged from the same glomerulus, comparing Predrug (bold line) and APV+NBQX conditions. (D) Top: ITA response amplitudes before and after APV+NBQX (5 mM/2.5 mM) application in pcMT cells from six experiments. In each experiment, odorant was presented at two concentrations varying by a factor of 2.5–10; each point indicates a particular glomerulus-concentration response, before and after drug application. Responses are normalized to the maximal pre-drug response in each experiment. Asterisks indicate experiments with significant change in ITA response amplitude. Summary statistics per experiment (mean ± SD APV+NBQX/Predrug ratios, paired t-tests with Bonferroni correction, n = glom − odor concentration pairs): Exp. 1, 1, 07 ± 0.35, p = 1.0, n = 13; Exp. 2, 0.83 ± 0.41, p = 1.0, n = 16; Exp. 3, 0.45 ± 0.15, p = 0.012, n = 8; Exp. 4, 0.58 ± 0.14, p = 3.7 × 10−4, n = 14; Exp 5, 0.97 ± 0.29, p = 1.0, n = 4; Exp 6, 0.53 ± 0.13, p = 8.2 × 10−3, n = 7. (E) Little change in onset latencies and peak times of inhalation-linked responses after APV+NBQX application. Left plot: change in latency measurements from CCK+ GCaMP6f datasets (persistent glomeruli, same data as Figure 3G), pcMT SF-iGluSnFR.A184V with 1 mM/0.5 mM APV+NBQX (data in Figure 4D), and pcMT SF-iGluSnFR.A148S with 5 mM/2.5 mM APV+NBQX. Onset latencies were not measured for the final dataset due to the relatively low signal-to-noise ratio of these responses. Right plot: Peak latencies relative to inhalation for each ROI in 5 mM/2.5 mM APV+NBQX dataset. Note no change in latencies for iGluSnFR signals. Boxes indicate 25th and 75th percentiles, line indicates median, and whiskers denote outliers with a coefficient of 1.5.
Figure 6. Impact ofmultisynaptic excitation on slow temporal dynamics of glutamatesignaling across inhalations.(A) Left: ΔF responsemaps showing MT cell glutamate responses to two odorants beforeapplication of APV+NBQX (1 mM/0.5 mM). Signals are fromSF-iGluSnFR.A184S injected into the OB of a Tbet-Cre mouse.Right: Traces (mean of four presentations) from two glomeruli (top,bottom) showing diverse temporal responses to different odorants, before and after APV+NBQX application. Response to cyclohexylamine (bottom) shows a loss of the initial transient and replacement by tonic signal. Response to vinyl butyrate (top) shows an increase in initial response magnitude and apparent loss of adaptation. Temporally distinct response patterns to other odorants persist after drug application. T1, T2, Tmax indicate time-points used to quantify slow temporal dynamics, as in (C). (B) Box plots (as in Figure 5) showing no overall impact of APV+NBQX on mean and peak response amplitudes (each measured across the duration of odorant presentation) across all glomerulus-odorant pairs (n = 83, three mice). (C) Scatter plot of T2 − T1/Tmax values before and after APV+NBQX for all responsive glomerulus-odor pairs across three experiments (each experiment indicated by color), showing decrease in adaptation for many pairs after application of APV+NBQX. Summary statistics: Paired t-tests, Exp. 1, n = 7 glomerulus-odor pairs, p = 0.004; Exp. 2: n = 34, p = 0.004; Exp. 3: n = 42, p = 1.63 × 10−5.
Data Analysis
Image analysis was performed using custom software written in MATLAB (Mathworks). For display, odorant response maps were displayed using ΔF values rather than ΔF/F to minimize noise from nonfluorescent regions. Activity maps were scaled as indicated in the figure and were kept to their original resolution (512 × 512 pixels) and smoothed using a Gaussian kernel with σ of 1 pixel. For time-series analysis, regions of interest (ROIs) were chosen manually based on the mean fluorescence image and were further refined based on odorant response ΔF maps, then all pixels averaged within an ROI to generate time-series for analysis. Time-series were computed and displayed as ΔF/F, with F defined as the mean fluorescence in the 1–2 s prior to odorant onset, and upsampled to 150 Hz for analysis using the MATLAB pchip function.
To analyze response patterns, traces were averaged across four-eight presentations of each odorant. Responses were classified as having significant excitatory and/or suppressive components (excitatory for Figures 1–6, suppressive only for Figure 2) as follows. First, each averaged response in the Predrug condition was z-scored using a baseline from 0.5 to 2 s before odorant onset, with z defined as the SD of the baseline period concatenated for each odorant response spanning each ROI. Peak excitation was measured as 95% of the maximum signal during the duration of the odorant presentation (4–8 s; when applicable, suppressive responses (i.e., Figure 2) were measured as the 15th percentile of all values in a time window from odorant onset to 500 ms after odorant offset. We then used a conservative criterion for the significance of z = ± 4 SD for identifying significant excitatory or suppressive responses for further analysis.
For analysis of inhalation-linked dynamics, in most experiments, responses to the first inhalation after odorant onset were analyzed, after averaging across 4–8 presentations. In other experiments (e.g., Figure 5) inhalation-triggered average (ITA) responses were generated by averaging each inhalation, repeated at 0.25 Hz, over a 70 s odorant presentation (17 inhalations averaged in total), as described previously (Short and Wachowiak, 2019). Onset and peak latencies, relative to the start of inhalation, were calculated from the ΔF/F time-series using the “risetime” function in MATLAB to fit a rising sigmoid to the baseline and peak ΔF/F levels and return the time to 10% of the rise from baseline to peak (onset latency) and the time to 90% of the rise from baseline to peak (peak latency). Latency values were calculated after upsampling but with no additional temporal filtering.
For analysis of odorant-evoked dynamics across multiple sniffs (i.e., Figure 6), responses to 4–8 presentations of odorant were averaged before analysis. Changes in response amplitude over time, or T2 − T1/Tmax, were calculated as the difference in amplitude between the peak ΔF/F following the first (T1) and the second-to-last (T2) inhalation during the odorant presentation, divided by the maximum ΔF/F during the 4-s odor presentation.
For statistical analysis of drug effects, each experiment was treated as an independent observation, since drug treatment was applied across all ROIs in an experiment. For datasests with five or more experiments, summary statistics are reported as the mean (or median, for non-normally-distributed data) ± standard deviation (SD) effect across all imaged ROIs, and comparisons made on these mean/median values. For datasets with fewer than five experiments, mean (or median) effects across ROIs are reported for each experiment and, when possible given the number of ROIs, tests for significance performed on each experiment, using ROIs as the independent measure, and p-values adjusted for multiple comparisons based on the number of experiments per dataset. Statistical tests were performed in Origin (OriginLab Corp.). All measurement of response parameters was done using analysis code that was independent of treatment or comparison condition.
Data and analysis code underlying the results of this study are available from the corresponding author upon request.
Results
Presynaptic Inhibition Does Not Impact Inhalation-Linked Dynamics of Glutamate Signaling
Neurotransmitter release from OSNs can be modulated by presynaptic dopamine (D2) and GABAB receptors, which modulate Ca2+ influx into the presynaptic terminal and subsequent release of the neurotransmitter (Ennis et al., 2001; Murphy et al., 2005; Wachowiak et al., 2005; Vaaga et al., 2017). GABAB-mediated presynaptic inhibition is activated both tonically and via sensory-evoked, intraglomerular circuits that have been proposed to mediate gain control in vivo in a manner independent of inhalation frequency (McGann et al., 2005; Pírez and Wachowiak, 2008; Shao et al., 2009; Vaaga and Westbrook, 2017). Thus we first sought to characterize how GABAB-mediated presynaptic inhibition modulates the gain and dynamics of the glutamatergic drive onto MT cells in vivo.
First, we confirmed effective blockade of GABAB-mediated presynaptic inhibition by imaging Ca2+ influx into OSN terminals in dorsal OB glomeruli using the genetically-encoded calcium reporter GCaMP6f expressed exclusively in OSNs (Wachowiak et al., 2013; Rothermel and Wachowiak, 2014; Short and Wachowiak, 2019) and comparing responses before and after topical application of the GABAB receptor antagonist CGP35348 (1 mM) to the dorsal OB, as done previously (Wachowiak et al., 2005; Brunert et al., 2016). Odorant was presented at suprathreshold concentrations (17–48 ppm; Table 1) to anesthetized mice using an artificial inhalation paradigm and evoked signals imaged from the dorsal OB under epifluorescence. Robust odorant-evoked presynaptic Ca2+ signals were apparent in multiple foci representing glomeruli across the dorsal OB (Figure 1A). Temporal response patterns during 2 Hz inhalation consisted of inhalation-linked transients apparent atop a sustained signal (Figure 1B). CGP35348 application caused an increase in the odorant-evoked Ca2+ signal that was apparent as early as the first inhalation of odorant and which persisted throughout the odorant presentation (Figure 1B). The magnitude of the increase varied across glomeruli and in different mice (Figure 1C), with a mean increase of 97% across the five mice tested (mean ± SD of median CGP35348/Predrug ratio per experiment: 1.97 ± 0.34; p = 0.003, t-test on median ratios with ratio = 1 as null hypothesis; df = 4). CGP35348 did not appear to change the dynamics of inhalation-linked GCaMP6f transients (Figure 1B, bottom). These results confirm previous findings using calcium-sensitive dyes (Wachowiak et al., 2005; Pírez and Wachowiak, 2008).
Next, we used the second-generation iGluSnFR, SF-iGluSnFR.A184V, to directly monitor the impact of GABAB receptor blockade on odorant-evoked glutamate signaling onto MT cells. As expected from an earlier characterization using first-generation iGluSnFR (Brunert et al., 2016), CGP35348 led to an increase in the amplitude of inhalation-driven glutamate transients (Figures 1D,E). Here, in contrast to the GCaMP6f imaging results, we observed a difference in the magnitude of the effect of GABAB receptor blockade on responses to the first and second inhalations of odorant. Overall, CGP35348 caused a significant increase in response to the second inhalation in all four mice but increased responses to the first inhalation of odorant in only three of the four mice (Figure 1F; see Figure 1 legend for statistical tests). The impact of CGP35348 on responses to the third and fourth inhalations appeared similar to that on the second inhalation (not shown).
However, GABAB receptor blockade had little to no impact on the inhalation-linked dynamics of glutamatergic signaling onto MT cells. Neither onset latencies nor peak latencies (defined as time to 10% of peak and time to 90% of peak, respectively) of the iGluSnFR transient evoked by the first inhalation of odorant changed substantially after CGP35348 application (Δ onset latency, −0.3 ± 2.6 ms; Δ peak latency, 3.1 ± 4.2 ms; mean ± SD of median values from four mice; Figure 1G). While two of the four mice showed a significant change in onset or peak latency (see Figure 1G legend), the magnitude of the change in these experiments (<10 ms) was negligible compared to the range of onset and peak latencies seen across different glomeruli in the same preparation, which ranges over a span of 200–250 ms (Moran et al., 2021). These results confirm that GABAB-mediated inhibition regulates the strength of glutamatergic signaling from presynaptic sources onto MT cells in vivo (McGann et al., 2005; Wachowiak et al., 2005; Brunert et al., 2016), and suggest that presynaptic inhibition contributes little to the inhalation-linked timing of these signals.
To assess the impact of this GABAB-mediated modulation of glutamatergic drive on postsynaptic responses of MT cells in vivo, we returned to Ca2+ imaging, expressing GCaMP6f in MT cells and comparing odorant-evoked responses before and after CGP35348 application. GCaMP6f signals were imaged from the apical tufts of MT cells across the glomerular layer using two-photon imaging, as done previously (Economo et al., 2016; Short and Wachowiak, 2019; Figures 2A,B). CGP35348 had variable effects on MT cell Ca2+ responses in different glomeruli, with increased responses in some glomeruli and decreased responses in others (Figures 2B,C). Suppressive responses reflected as a decrease in GCaMP6f fluorescence and presumably reflecting inhibition of ongoing MT cell spiking (Economo et al., 2016), persisted after application of CGP35348 (Figures 2B,C) indicating that MT cell suppression is largely not mediated by presynaptic inhibition of glutamatergic drive from OSNs. Overall, there was no significant net change in response magnitude, either for peak responses to the first inhalation of odorant, or averaged across the duration of odorant presentation (CGP35348/Predrug ratios, first inhalation of odorant (median [quartile range] of medians: 0.75 [0.64 ± 1.05], p = 0.11, one-sample Wilcoxon Signed Rank test, n = 12; odor duration: 0.91 [0.77 ± 1.33], p = 0.85). Thus, increasing the magnitude of excitatory input to the glomerulus by the enhanced release of glutamate from OSN terminals does not uniformly translate into increased excitation of MT cells.
Limited Contribution of Multisynaptic Pathways to Odorant-Evoked Glutamate Signaling Onto MT Cells
Sensory-evoked glutamatergic input onto MT cells can arise from OSNs or from multisynaptic excitation involving dendritic glutamate release from ET cells and, possibly, superficial tufted cells (sTCs; Hayar et al., 2004a; De Saint Jan et al., 2009; Gire et al., 2012; Vaaga and Westbrook, 2016; Sun et al., 2020). In addition, inhibitory circuits can modulate multisynaptic glutamate signaling (Gire and Schoppa, 2009; Shao et al., 2012; Gire et al., 2019). Glomerular glutamate signals may also arise from glutamate release by MT cell dendrites (Isaacson, 1999; Najac et al., 2015), and thus at least partially reflect MT cell excitation itself. To isolate the direct contribution of OSN inputs to MT cell excitation dynamics, we compared evoked glutamate dynamics before and after pharmacological blockade of postsynaptic activity with ionotropic glutamate receptor antagonists (Gurden et al., 2006; Pírez and Wachowiak, 2008; Lecoq et al., 2009). Because iGluSnFRs are insensitive to these antagonists (Marvin et al., 2013, 2018), we reasoned that this approach would largely prevent OSN-driven excitation of postsynaptic neurons and allow us to image evoked MT cell glutamate signals arising solely from OSN inputs.
To confirm our ability to block ionotropic glutamate receptors with the in vivo drug application protocol, we first continued the GCaMP6f imaging from MT cells in the same mice used for the GABAB blockade experiments, applying a cocktail of APV+NBQX (1 mM/0.5 mM) following the initial CGP35348 treatment (n = 4 mice). APV+NBQX sharply reduced odorant-evoked GCaMP6f signals, such that distinct inhalation-linked transients were eliminated in most glomeruli and replaced with a weak tonic glutamate signal (Figure 2B); this signal might reflect modulation by metabotropic glutamate receptors (Matsumoto et al., 2009; Dong and Ennis, 2014), which we did not attempt to block in these experiments. Excitatory responses to both the first inhalation and over the duration of odorant presentation were reduced to 36 ± 20% (p = 3.5 × 10−10, n = 29 ROIs from four mice, paired t-test); and 36 ± 15% (p = 1.2 × 10−6, n = 23 ROIs, paired Wilcoxon Signed Ranks) respectively, of their baseline values (Figures 2B,C; see Figure legend for summary statistics). Suppressive responses were eliminated after APV+NBQX application and were replaced with a slow, tonic increase (Figure 2B).
We next tested the ability of APV+NBQX to block activation of cholecystokinin (CCK)-expressing sTCs, as these, like ET cells, are strongly driven by monosynaptic OSN input (Isaacson, 1999; Najac et al., 2015). We imaged activation of CCK+ sTCs at their dendrites in CCK-IRES-Cre mice crossed to a Cre-dependent GCaMP6f reporter line (n = 3 experiments from two mice, see “Materials and Methods” section). We first blocked GABAB-mediated presynaptic inhibition with CGP35348 to further enhance transmitter release from OSNs. As expected, CGP35348 led to a modest increase in the peak CCK+ GCaMP6f response, with median increases of 33, 9, and 24%, respectively, in the three experiments (Figures 3A–E). Subsequent application of APV+NBQX (1 mM/0.5 mM) had mixed effects in different glomeruli, with responses completely or nearly eliminated in 8 of 14 glomeruli (97 ± 3% median reduction in peak response to 1st inhalation) but only partially reduced in the remaining six (38 ± 16% median reduction; Figures 3A–E). In these six glomeruli, responses to successive inhalations after odorant onset were reduced even further or eliminated (e.g., Figure 3E). We performed additional experiments using a 5× higher concentration of antagonist (5 mM APV/2.5 mM NBQX) and without preapplication of CGP35348. In these, APV+NBQX completely blocked postsynaptic activation in 11 of 13 glomeruli (from three experiments), with strongly reduced responses (66 and 75% reduction) in the remaining two (Figure 3F).
Importantly, in glomeruli with persistent (but weakened) responses after APV+NBQX application, inhalation-linked Ca2+ transients were significantly delayed (Figures 3B,D), with a mean increase in onset and peak latency of 255 ± 79 ms and 257 ± 51 ms, respectively (Figure 3G). There was no change in latencies after CGP35348 application alone (Figure 3G; see Figure legend for summary statistics). An explanation for this substantial delay is that, in glomeruli receiving the strongest sensory input, glutamate concentrations are eventually able to overcome APV+NBQX blockade sufficiently to trigger spike bursts in sTCs. Overall, these results demonstrate that APV+NBQX blocks postsynaptic activation in glomeruli receiving all but the strongest OSN inputs, and even in those glomeruli, substantially weakens and delays activation of monosynaptically-driven sTCs.
We next used in vivo pharmacology to test the contribution of mono- vs. di- or polysynaptic glutamate signaling to inhalation-linked glutamate transients onto MT cells. We focused on MT cells projecting to the piriform cortex (pcMTs) as the multisynaptic excitatory drive is thought to predominately impact this subpopulation (Fukunaga et al., 2012; Gire et al., 2012; Igarashi et al., 2012). We selectively expressed SF-iGluSnFR.A184V in pcMTs via retrograde viral infection by targeting virus injections to the anterior piriform cortex of Tbet-Cre mice (n = 2 mice), as described previously (Rothermel et al., 2013). We have established in a recent study that odorant- and inhalation-evoked glutamate transients in pcMTs are indistinguishable from those measured from sTCs or from the general MT cell population (Moran et al., 2021). As expected, inhalation-linked glutamate transients on pcMTs were increased by CGP35348 (Figures 4A,B), with mean increases of 35% and 90% in the two mice, and no impact on inhalation-linked glutamate dynamics (Δ onset and peak latencies, mean ± SD, paired t-tests, n = 10 ROIs pooled across two mice: Δ onset, 8 ± 11 ms, p = 0.07; Δ peak, 8 ± 15 ms, p = 0.2). Surprisingly, subsequent application of APV+NBQX (1 mM/0.5 mM) had no significant impact on response amplitudes (Figures 4A–C; see Figure 4C legend for summary statistics). iGluR blockade did cause a very small, but statistically significant, increase in the latency of the inhalation-linked glutamate transient (Δ onset latency, 13 ± 12 ms, p = 0.01; Δ peak, 34 ± 16 ms, p = 1.2 × 10−4; Figure 4D); however, the magnitude of this change was very small relative to the overall dynamics of the inhalation-linked transient (Figure 4B).
Multisynaptic glutamatergic excitation may preferentially drive MT cell responses to weak inputs (Najac et al., 2011; Vaaga and Westbrook, 2016; Gire et al., 2019), and the odorant concentrations used in the preceding experiments were relatively high (~20–50 ppm; Table 1). Removal of presynaptic inhibition with CGP35348 may also enhance OSN input strength to physiologically excessive levels. Thus, we next tested the impact of glutamate receptor blockade on inhalation-linked transients in pcMT cell responses to odorants presented at 10–100× lower concentrations (0.1–5 ppm), using a higher concentration of APV+NBQX (5 mM/2.5 mM) and omitting CGP35348. In addition, to allow for greater sensitivity of glutamate detection and to focus on inhalation-linked dynamics of the glutamate signal, we used the higher-affinity SF-iGluSnFR.A184S and modified the artificial inhalation protocol to generate averaged responses to a single inhalation repeated at 0.25 Hz (Figure 5A). This protocol yielded inhalation-triggered average waveforms that showed substantial variation in onset and peak latency as well as duration, as we have reported previously (Short and Wachowiak, 2019; Moran et al., 2021). Finally, in a subset of mice (n = 2), we used dual-color imaging to simultaneously image presynaptic glutamate transients and postsynaptic Ca2+ signals in the apical glomerular tufts of overlapping MT cell populations (Figure 5B), using Tbet-Cre mice crossed with the Thy1-jRGeCO1a reporter line (Dana et al., 2018; see “Materials and Methods” section).
In the two dual-color imaging preparations, APV+NBQX completely blocked postsynaptic MT cell activation in 7 of 11 glomeruli analyzed, as reflected in the peak jRGECO1a Ca2+ signal imaged from the glomerular layer, with a reduction to 9, 12, 12 and 34% of pre-drug levels in the remaining four glomeruli. APV+NBQX had minimal impact on the glutamate transients measured in these same glomeruli (Figure 5C). Overall, across all aPC-injected SF-iGluSnFR.A184V mice in this dataset (n = 6), the impact of APV+NBQX on the magnitude of inhalation-linked glutamate transients was modest and variable, causing a statistically significant reduction in three of six experiments (see Figure 4 legend) and a 27 ± 26% reduction in peak inhalation-triggered average (ITA) amplitude across all experiments (mean ± SD of mean ratios, p = 0.06, t-test, n = 6; Figure 5D). However, APV+NBQX did not impact the dynamics of inhalation-linked glutamate responses, with no significant change in peak times across the dataset (mean ± SD of Δ peak latencies, 9 ± 21 ms, p = 0.34, paired t-test, n = 6; Figure 5E). Overall, these results suggest that multisynaptic pathways modestly contribute to the strength of glutamatergic excitation of MT cells in response to low-intensity stimulation, but do not uniquely shape the inhalation-linked dynamics of this excitation.
Finally, we sought to further examine the role of direct vs. multisynaptic excitation in shaping the dynamics of glutamatergic input to MT cells during repeated odorant sampling. We have shown previously that glutamate signals on MT cells show diverse temporal patterns across repeated inhalations of the odorant that include adaptation, facilitation, and even suppression depending on the odorant identity and concentration (Moran et al., 2021). To capture this diversity, we used a 12-odorant panel (delivered concentrations ranging from 0.03 to 15 ppm) and 3 Hz artificial inhalation, comparing response patterns before and after application of APV+NBQX (1 mM/0.5 mM). This approach yielded significant odorant responses in 83 glomerulus-odorant pairs across three mice and included diverse excitatory response patterns (Figure 6A). Consistent with the earlier experiments, APV+NBQX application alone had little overall impact on the magnitude of evoked glutamate signals: while there was variability in the effect for certain glomerulus-odor pairs, neither mean nor peak excitatory response showed a significant overall change after APV+NBQX application (paired t-tests, p = 0.147 for mean amplitude, p = 0.605 for peak amplitude, n = 83; Figure 6B).
APV+NBQX did, however, impact the temporal patterns of the glutamate signal for a subset of responses. In particular, application of APV+NBQX alone tended to reduce adaptation that occurred following the initial inhalation of odorant (e.g., vinyl butyrate response; Figure 6A), or it enhanced long-lasting, tonic-type responses (e.g., ethyl butyrate response; Figure 6A). We used a simple metric, T2 − T1/Tmax, (Moran et al., 2021) to quantify and compare these adapting or facilitating dynamics before and after drug application (Figure 6C, see also “Materials and Methods” section). APV+NBQX caused an overall increase in T2 − T1/Tmax values (Figure 6C), with significant increases seen in each of the three preparations. These results suggest that multisynaptic pathways contribute to the slower temporal patterning of MT cell activity by modestly suppressing the excitatory drive of MT cells across repeated samples of odorant.
Discussion
Neural circuits in the glomerular layer of the OB mediate the initial synaptic processing of olfactory inputs, yet their role in shaping the dynamic response patterns of MT cells that are thought to encode odor information during natural odor sampling remains unclear. Because of the higher complexity of MT cell response dynamics as compared to those of OSNs—including the emergence of suppressive components—much attention has focused on the role of inhibitory circuits in this process (Shao et al., 2013; Fukunaga et al., 2014; Liu et al., 2016). However, in a preceding study (Moran et al., 2021) we found that diverse temporal patterning is prominent already at the level of excitatory, glutamatergic input to MT cells in vivo, suggesting that feedforward excitatory circuits may play an underappreciated role in generating diverse MT cell responses. Here, we used in vivo pharmacology to more directly probe the contribution of different sources of glutamatergic input onto MT cells. Overall, our results suggest a predominant role for direct synaptic input from OSNs to MT cells in determining the dynamics of sensory-driven excitatory drive in response to odorant inhalation.
OSNs themselves show remarkable temporal diversity in their odorant responses, with activation latencies relative to inhalation spanning a range of ~250 ms, and inhalation-driven response durations varying by a similar amount in an odorant- and glomerulus-specific fashion as measured by presynaptic Ca2+ imaging from OSN axon terminals (Spors et al., 2006; Carey and Wachowiak, 2011; Wachowiak et al., 2013; Short and Wachowiak, 2019). While it is not surprising that this diversity persists at the level of glutamatergic signaling, our results indicate a relatively small contribution of OB circuits to further shaping inhalation-linked MT cell excitatory drive.
In addition to direct OSN inputs, feedforward, disynaptic excitation of MT cells by external tufted (ET) cells has been proposed as a primary driver of MT cell excitation (De Saint Jan et al., 2009; Gire et al., 2012). Studies from OB slices suggest that ET cells provide excitatory drive to MT cells via glutamate spillover from the ET- to periglomerular cell synapse (Najac et al., 2011; Gire et al., 2019), allowing inhibitory circuits to regulate MT cell excitability and temporal dynamics through their inhibitory action on ET cells (Murphy et al., 2005; Gire and Schoppa, 2009; Whitesell et al., 2013; Banerjee et al., 2015; Liu et al., 2016). Here, we found that a near-complete blockade of postsynaptic activation, as confirmed with GCaMP reporters expressed in monosynaptically-driven sTCs (Sun et al., 2020), had no impact on the dynamics of inhalation-linked glutamate transients and only modestly impacted the amplitude of transients in response to low odorant concentrations.
These results do not rule out some models of ET cell-mediated excitation of MT cells, which propose that the ET cell pathway is most important in regimes of weak OSN input (i.e., low odorant concentrations; Najac et al., 2011; Vaaga and Westbrook, 2016; Gire et al., 2019). While the overall conclusion of strong OSN-mediated glutamatergic drive is supported by multiple datasets in the present study, the relatively small number of mice in each cell-type/drug treatment dataset limits quantitative assessments of the concentration- or intensity-dependence of this pathway. Direct OSN-driven currents may also be shunted by gap junctions between MT cells (Gire et al., 2012), which would not be reflected in our iGluSnFR recordings. Disynaptic excitation might also provide tonic glutamatergic drive to MT cells via ET cell bursting in vivo, allowing for modulation of MT cell excitability by inhibitory circuits (Hayar et al., 2004b; Hayar and Ennis, 2007). Overall, our findings illustrate the difficulty in extrapolating from OB slice experiments to OB circuit function in vivo; for example, glutamate spillover may contribute less to MT cell drive in vivo than in OB slices due to differences in glutamate transporter efficacy or in the dynamics of odorant-evoked glutamate release from OSN inputs.
One potential concern in interpreting these results is that the in vivo pharmacological approach did not completely block odorant-evoked glutamatergic transmission. While one previous study used higher concentrations of APV (50 mM) and NBQX (5 mM) to block glutamatergic transmission in vivo (Gurden et al., 2006), most prior studies have used concentrations comparable to or lower than the 5 mM/2.5 mM concentrations used here to block postsynaptic activation in the OB in vivo (Pírez and Wachowiak, 2008; Fletcher et al., 2009; Lecoq et al., 2009; Shang and Xing, 2018). We also found that the 2.5 mM NBQX concentration was near the limit of its solubility in normal Ringers solution and that higher NBQX concentrations introduced severe optical interference with the imaging protocol. Nonetheless, it is conceivable that iGluR blockade failed to impact ET cell activation, and that, in the Ca2+ imaging experiments used to verify effective blockade we failed to observe persistent ET cell activation due to a lack of GCaMP expression in these cells.
Several lines of evidence suggest this is unlikely. First, GCaMP signals imaged from CCK-ergic sTCs, which, like ET cells, receive strong monosynaptic input from OSNs (Sun et al., 2020), were completely blocked in most glomeruli and strongly reduced and delayed in the remainder. Second, despite this delay in the persistent sTC response, we observed no corresponding delay in the iGluSnFR response after APV+NBQX treatment. Third, in simultaneous dual-color SF-iGluSnFR and jRGECO1a imaging experiments, we saw only a modest reduction in SF-iGluSnFR response amplitude and no change in SF-iGluSnFR response latency, despite a complete or near-complete blockade of postsynaptic responses imaged in the same glomeruli. Given these results, it seems unlikely that ET cell activation and subsequent glutamate release could have been unperturbed (and undetected) by the high iGluR blocker concentrations used here.
Blocking presynaptic modulation of transmitter release from OSNs using GABAB receptor antagonists impacted the magnitude, but not the dynamics, of inhalation-linked glutamate signaling onto MT cells. This result is predicted from earlier studies based on Ca2+ imaging from OSN terminals or synaptopHluorin-based measurements of transmitter release, which found that presynaptic inhibition mediates gain control of sensory input to the OB (McGann et al., 2005; Wachowiak et al., 2005). Blocking presynaptic inhibition had no impact on the presence of suppressive responses in MT cells, as measured with Ca2+ imaging, indicating that this suppression is mediated by feedforward inhibition rather than inter- or intraglomerular suppression of ongoing drive from OSNs. In some cases, GABAB receptor blockade had less impact on the first inhalation of odorant than on subsequent inhalations (i.e., Figure 1E), an effect that we have not observed with presynaptic Ca2+ imaging (compare Figures 1B,E; see also Pírez and Wachowiak, 2008). This result is expected if the initial inhalation-driven burst of OSN activity is sufficiently fast and synchronous to drive glutamate release before GABAB-mediated inhibition can be activated, which takes approximately 50 ms to reach maximal strength due to G-protein coupled signaling in the presynaptic terminal (Wachowiak et al., 2005).
Glutamate signaling also shows substantial temporal diversity across repeated inhalations of odorant in both anesthetized and awake mice (Moran et al., 2021), and blocking multisynaptic activation with APV+NBQX altered these temporal patterns, generally by reducing adaptation of the glutamate signal or enhancing its facilitation across successive inhalations. This effect is consistent with disinhibition of glutamatergic drive by APV+NBQX, and could reflect the removal of feedback presynaptic inhibition of OSN terminals (Wachowiak et al., 2005; Shao et al., 2009) or blockade of feedforward inhibition of ET cells; GABAB receptor blockade using CGP35348 prior to APV+NBQX application, as we did with low-frequency inhalations, should distinguish these possibilities. Lastly, we did not attempt to block metabotropic glutamate receptors in these experiments, although they have been implicated in diverse contributions to MT cell and ET cell excitability, including mediating prolonged odorant-evoked MT cell responses in vivo (Matsumoto et al., 2009; Dong and Ennis, 2014, 2017). Despite the need for further dissection of the transmitter pathways involved in glomerular processing, our results suggest that multisynaptic circuits may contribute more to slow temporal patterning of glutamatergic drive to MT cells than to shaping inhalation-linked excitatory dynamics.
MT cells show a striking correspondence between temporal patterns of glutamatergic input to their dendrites and their postsynaptic activity patterns, as measured with simultaneous Ca2+ imaging (Moran et al., 2021). A caveat to this similarity is the possibility that the SF-iGluSnFR signal reflects the dendritic release of glutamate from MT cells themselves (Isaacson, 1999; Najac et al., 2015). However, we found that near-complete pharmacological blockade of MT cell activation enhanced, rather than suppressed, responses across repeated inhalation, suggesting that the SF-iGluSnFR signal largely reflects glutamatergic input to, rather than glutamate release from, MT cells. Taken together, these results suggest that direct glutamatergic input from OSNs provides the bulk of excitatory drive to MT cells, and that diversity in the dynamics of this input may be a primary determinant of the temporal diversity in MT cell responses that underlies odor representations at this stage. Experiments using glutamate or other transmitter reporters, multiplexed with measures of postsynaptic activation, will be important in further unraveling the contributions of OB circuits to shaping OB outputs in vivo, and as a function of odorant sampling in the behaving animal.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, upon reasonable request.
Ethics Statement
The animal study was reviewed and approved by University of Utah IACUC.
Author Contributions
AKM and MW designed the study. AKM, TPE, and MW performed research, analyzed data, and wrote the article. All authors contributed to the manuscript editing and approved the submitted version.
Funding
This work was supported by National Institute of Neurological Disorders and Stroke (NINDS) T32NS076067 (to AKM), and the National Institute of Deafness and Other Communication Disorders (NIDCD) R01DC006441 (to MW), R01DC006441-13S1 (to AKM), and F32DC015389 (to TPE).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank L. Looger and J. Marvin and the GENIE Project from HHMI Janelia Research Campus for sharing the SF-iGluSnFR virus. We also thank G. Vasquez-Opazo, R. Kummer, and J. Ball for technical assistance, T. Rust for data analysis software, as well as S. Burton, S. Short, I. Youngstrom, K. Podgorski, and D. Restrepo for providing helpful feedback and discussion.
References
Banerjee, A., Marbach, F., Anselmi, F., Koh, M. S., Davis, M. B., Garcia Da Silva, P., et al. (2015). An interglomerular circuit gates glomerular output andimplements gain control in the mouse olfactory bulb. Neuron 87, 193–207. doi: 10.1016/j.neuron.2015.06.019
Bozza, T., Mcgann, J. P., Mombaerts, P., and Wachowiak, M. (2004). in vivo imaging of neuronal activity by targeted expression of a genetically encoded probe in the mouse. Neuron 42, 9–21. doi: 10.1016/s0896-6273(04)00144-8
Brunert, D., Tsuno, Y., Rothermel, M., Shipley, M. T., and Wachowiak, M. (2016). Cell-type-specific modulation of sensory responses in olfactory bulb circuits by serotonergic projections from the raphe nuclei. J. Neurosci. 36, 6820–6835. doi: 10.1523/JNEUROSCI.3667-15.2016
Burton, S. D., Wipfel, M., Guo, M., Eiting, T. P., and Wachowiak, M. (2019). A novel olfactometer for efficient andflexible odorant delivery. Chem. Senses 44, 173–188. doi: 10.1093/chemse/bjz005
Carey, R. M., and Wachowiak, M. (2011). Effect of sniffing on the temporal structure of mitral/tufted cell output from the olfactory bulb. J. Neurosci. 31, 10615–10626. doi: 10.1523/JNEUROSCI.1805-11.2011
Chong, E., Moroni, M., Wilson, C., Shoham, S., Panzeri, S., and Rinberg, D. (2020). Manipulating synthetic optogenetic odors reveals the coding logic of olfactory perception. Science 368:eaba2357. doi: 10.1126/science.aba2357
Dana, H., Novak, O., Guardado-Montesino, M., Fransen, J. W., Hu, A., Borghuis, B. G., et al. (2018). Thy1 transgenic mice expressing the red fluorescent calcium indicator jRGECO1a for neuronal population imaging in vivo. PLoS One 13:e0205444. doi: 10.1371/journal.pone.0205444
De Saint Jan, D., Hirnet, D., Westbrook, G. L., and Charpak, S. (2009). External tufted cells drive the output of olfactory bulb glomeruli. J. Neurosci. 29, 2043–2052. doi: 10.1523/JNEUROSCI.5317-08.2009
Dong, H.-W., and Ennis, M. (2014). Activation of group I metabotropic glutamate receptors enhances persistent sodium current andrhythmic bursting in main olfactory bulb external tufted cells. J. Neurophysiol. 111, 641–647. doi: 10.1152/jn.00696.2013
Dong, H. W., and Ennis, M. (2017). Activation of group II metabotropic glutamate receptors suppresses excitability of mouse main olfactory bulb external tufted andmitral cells. Front. Cell. Neurosci. 11:436. doi: 10.3389/fncel.2017.00436
Economo, M. N., Hansen, K. R., and Wachowiak, M. (2016). Control of mitral/tufted cell output by selective inhibition among olfactory bulb glomeruli. Neuron 91, 397–411. doi: 10.1016/j.neuron.2016.06.001
Eiting, T. P., and Wachowiak, M. (2018). Artificial inhalation protocol in adult mice. Bio Protoc. 8:e3024. doi: 10.21769/BioProtoc.3024
Eiting, T. P., and Wachowiak, M. (2020). Differential impacts of repeated sampling on odor representations by genetically-defined mitral andtufted cell subpopulations in the mouse olfactory bulb. J. Neurosci. 40, 6177–6188. doi: 10.1523/JNEUROSCI.0258-20.2020
Ennis, M., Zhou, F.-M., Ciombor, K. J., Aroniadou-Anderjaska, V., Hayar, A., Borrelli, E., et al. (2001). Dopamine D2 receptor-mediated presynaptic inhibition of olfactory nerve terminals. J. Neurophysiol. 86, 2986–2997. doi: 10.1152/jn.2001.86.6.2986
Fletcher, M. L., Masurkar, A. V., Xing, J., Imamura, F., Xiong, W., Nagayama, S., et al. (2009). Optical imaging of postsynaptic odor representation in the glomerular layer of the mouse olfactory bulb. J. Neurophysiol. 102, 817–830. doi: 10.1152/jn.00020.2009
Fukunaga, I., Berning, M., Kollo, M., Schmaltz, A., and Schaefer, A. T. (2012). Two distinct channels of olfactory bulb output. Neuron 75, 320–329. doi: 10.1016/j.neuron.2012.05.017
Fukunaga, I., Herb, J. T., Kollo, M., Boyden, E. S., and Schaefer, A. T. (2014). Independent control of gamma andtheta activity by distinct interneuron networks in the olfactory bulb. Nat. Neurosci. 17, 1208–1216. doi: 10.1038/nn.3760
Gire, D. H., Franks, K. M., Zak, J. D., Tanaka, K. F., Whitesell, J. D., Mulligan, A. A., et al. (2012). Mitral cells in the olfactory bulb are mainly excited through a multistep signaling path. J. Neurosci. 32, 2964–2975. doi: 10.1523/JNEUROSCI.5580-11.2012
Gire, D. H., and Schoppa, N. E. (2009). Control of on/off glomerular signaling by a local GABAergic microcircuit in the olfactory bulb. J. Neurosci. 29, 13454–13464. doi: 10.1523/JNEUROSCI.2368-09.2009
Gire, D. H., Zak, J. D., Bourne, J. N., Goodson, N. B., and Schoppa, N. E. (2019). Balancing extrasynaptic excitation andsynaptic inhibition within olfactory bulb glomeruli. eNeuro 6:ENEURO.0247-19.2019. doi: 10.1523/ENEURO.0247-19.2019
Gurden, H., Uchida, N., and Mainen, Z. F. (2006). Sensory-evoked intrinsic optical signals in the olfactory bulb are coupled to glutamate release anduptake. Neuron 52, 335–345. doi: 10.1016/j.neuron.2006.07.022
Haddad, R., Lanjuin, A., Madisen, L., Zeng, H., Murthy, V. N., and Uchida, N. (2013). Olfactory cortical neurons read out a relative time code in the olfactory bulb. Nat. Neurosci. 16, 949–957. doi: 10.1038/nn.3407
Hayar, A., and Ennis, M. (2007). Endogenous GABA andglutamate finely tune the bursting of olfactory bulb external tufted cells. J. Neurophysiol. 98, 1052–1056. doi: 10.1152/jn.01214.2006
Hayar, A., Karnup, S., Ennis, M., and Shipley, M. T. (2004a). External tufted cells: a major excitatory element that coordinates glomerular activity. J. Neurosci. 24, 6676–6685. doi: 10.1523/JNEUROSCI.1367-04.2004
Hayar, A., Karnup, S., Shipley, M. T., and Ennis, M. (2004b). Olfactory bulb glomeruli: external tufted cells intrinsically burst at theta frequency andare entrained by patterned olfactory input. J. Neurosci. 24, 1190–1199. doi: 10.1523/JNEUROSCI.4714-03.2004
Igarashi, K. M., Ieki, N., An, M., Yamaguchi, Y., Nagayama, S., Kobayakawa, K., et al. (2012). Parallel mitral andtufted cell pathways route distinct odor information to different targets in the olfactory cortex. J. Neurosci. 32, 7970–7985. doi: 10.1523/JNEUROSCI.0154-12.2012
Isaacson, J. S. (1999). Glutamate spillover mediates excitatory transmission in the rat olfactory bulb. Neuron 23, 377–384. doi: 10.1016/s0896-6273(00)80787-4
Lecoq, J., Tiret, P., and Charpak, S. (2009). Peripheral adaptation codes for high odor concentration in glomeruli. J. Neurosci. 29, 3067–3072. doi: 10.1523/JNEUROSCI.6187-08.2009
Liu, S., Puche, A. C., and Shipley, M. T. (2016). The interglomerular circuit potently inhibits olfactory bulb output neurons by both direct andindirect pathways. J. Neurosci. 36, 9604–9617. doi: 10.1523/JNEUROSCI.1763-16.2016
Marvin, J. S., Borghuis, B. G., Tian, L., Cichon, J., Harnett, M. T., Akerboom, J., et al. (2013). An optimized fluorescent probe for visualizing glutamate neurotransmission. Nat. Methods 10, 162–170. doi: 10.1038/nmeth.2333
Marvin, J. S., Scholl, B., Wilson, D. E., Podgorski, K., Kazemipour, A., Muller, J. A., et al. (2018). Stability, affinity andchromatic variants of the glutamate sensor iGluSnFR. Nat. Methods 15, 936–939. doi: 10.1038/s41592-018-0171-3
Matsumoto, H., Kashiwadani, H., Nagao, H., Aiba, A., and Mori, K. (2009). Odor-induced persistent discharge of mitral cells in the mouse olfactory bulb. J. Neurophysiol. 101, 1890–1900. doi: 10.1152/jn.91019.2008
McGann, J. P., Pírez, N., Gainey, M. A., Muratore, C., Elias, A. S., and Wachowiak, M. (2005). Odorant representations are modulated by intra- but not interglomerular presynaptic inhibition of olfactory sensory neurons. Neuron 48, 1039–1053. doi: 10.1016/j.neuron.2005.10.031
Moran, A. K., Eiting, T. P., and Wachowiak, M. (2021). Dynamics of glutamatergic drive underlie diverse responses of olfactory bulb outputs in vivo. eNeuro 8:ENEURO.0110-21.2021. doi: 10.1523/ENEURO.0110-21.2021
Murphy, G. J., Darcy, D. P., and Isaacson, J. S. (2005). Intraglomerular inhibition: signaling mechanisms of an olfactory microcircuit. Nat. Neurosci. 8, 354–364. doi: 10.1038/nn1403
Najac, M., De Saint Jan, D., Reguero, L., Grandes, P., and Charpak, S. (2011). Monosynaptic andpolysynaptic feed-forward inputs to mitral cells from olfactory sensory neurons. J. Neurosci. 31, 8722–8729. doi: 10.1523/JNEUROSCI.0527-11.2011
Najac, M., Sanz Diez, A., Kumar, A., Benito, N., Charpak, S., and De Saint Jan, D. (2015). Intraglomerular lateral inhibition promotes spike timing variability in principal neurons of the olfactory bulb. J. Neurosci. 35, 4319–4331. doi: 10.1523/JNEUROSCI.2181-14.2015
Patterson, M. A., Lagier, S., and Carleton, A. (2013). Odor representations in the olfactory bulb evolve after the first breath andpersist as an odor afterimage. Proc. Natl. Acad. Sci. U S A 110, E3340–E3349. doi: 10.1073/pnas.1303873110
Pírez, N., and Wachowiak, M. (2008). in vivo modulation of sensory input to the olfactory bulb by tonic andactivity-dependent presynaptic inhibition of receptor neurons. J. Neurosci. 28, 6360–6371. doi: 10.1523/JNEUROSCI.0793-08.2008
Rothermel, M., Brunert, D., Zabawa, C., Díaz-Quesada, M., and Wachowiak, M. (2013). Transgene expression in target-defined neuron populations mediated by retrograde infection with adeno-associated viral vectors. J. Neurosci. 33, 15195–15206. doi: 10.1523/JNEUROSCI.1618-13.2013
Rothermel, M., and Wachowiak, M. (2014). Functional imaging of cortical feedback projections to the olfactory bulb. Front. Neural Circuits 8:73. doi: 10.3389/fncir.2014.00073
Schaefer, A. T., and Margrie, T. W. (2007). Spatiotemporal representations in the olfactory system. Trends Neurosci. 30, 92–100. doi: 10.1016/j.tins.2007.01.001
Shang, M., and Xing, J. (2018). Blocking of dendrodendritic inhibition unleashes widely spread lateral propagation of odor-evoked activity in the mouse olfactory bulb. Neuroscience 391, 50–59. doi: 10.1016/j.neuroscience.2018.09.003
Shao, Z., Puche, A. C., Kiyokage, E., Szabo, G., and Shipley, M. T. (2009). Two GABAergic intraglomerular circuits differentially regulate tonic andphasic presynaptic inhibition of olfactory nerve terminals. J. Neurophysiol. 101, 1988–2001. doi: 10.1152/jn.91116.2008
Shao, Z., Puche, A. C., Liu, S., and Shipley, M. T. (2012). Intraglomerular inhibition shapes the strength andtemporal structure of glomerular output. J. Neurophysiol. 108, 782–793. doi: 10.1152/jn.00119.2012
Shao, Z., Puche, A. C., and Shipley, M. T. (2013). Intraglomerular inhibition maintains mitral cell response contrast across input frequencies. J. Neurophysiol. 110, 2185–2191. doi: 10.1152/jn.00023.2013
Short, S. M., and Wachowiak, M. (2019). Temporal dynamics of inhalation-linked activity across defined subpopulations of mouse olfactory bulb neurons imaged in vivo. eNeuro 6:ENEURO.0189-0119.2019. doi: 10.1523/ENEURO.0189-19.2019
Spors, H., Wachowiak, M., Cohen, L. B., and Friedrich, R. W. (2006). Temporal dynamics andlatency patterns of receptor neuron input to the olfactory bulb. J. Neurosci. 26, 1247–1259. doi: 10.1523/JNEUROSCI.3100-05.2006
Stern, M., Bolding, K. A., Abbott, L. F., and Franks, K. M. (2018). A transformation from temporal to ensemble coding in a model of piriform cortex. eLife 7:e34831. doi: 10.7554/eLife.34831
Sun, X., Liu, X., Starr, E. R., and Liu, S. (2020). CCKergic tufted cells differentially drive two anatomically segregated inhibitory circuits in the mouse olfactory bulb. J. Neurosci. 40, 6189–6206. doi: 10.1523/JNEUROSCI.0769-20.2020
Vaaga, C. E., and Westbrook, G. L. (2016). Parallel processing of afferent olfactory sensory information. J. Physiol. 594, 6715–6732. doi: 10.1113/JP272755
Vaaga, C. E., and Westbrook, G. L. (2017). Distinct temporal filters in mitral cells andexternal tufted cells of the olfactory bulb. J. Physiol. 595, 6349–6362. doi: 10.1113/JP274608
Vaaga, C. E., Yorgason, J. T., Williams, J. T., and Westbrook, G. L. (2017). Presynaptic gain control by endogenous cotransmission of dopamine andGABA in the olfactory bulb. J. Neurophysiol. 117, 1163–1170. doi: 10.1152/jn.00694.2016
Wachowiak, M. (2011). All in a sniff: olfaction as a model for active sensing. Neuron 71, 962–973. doi: 10.1016/j.neuron.2011.08.030
Wachowiak, M., Economo, M. N., Díaz-Quesada, M., Brunert, D., Wesson, D. W., White, J. A., et al. (2013). Optical dissection of odor information processing in vivo using gcamps expressed in specified cell types of the olfactory bulb. J. Neurosci. 33, 5285–5300. doi: 10.1523/JNEUROSCI.4824-12.2013
Wachowiak, M., Mcgann, J. P., Heyward, P. M., Shao, Z., Puche, A. C., and Shipley, M. T. (2005). Inhibition of olfactory receptor neuron input to olfactory bulb glomeruli mediated by suppression of presynaptic calcium influx. J. Neurophysiol. 94, 2700–2712. doi: 10.1152/jn.00286.2005
Keywords: imaging, pharmacology, neural dynamics, iGluSnFR, presynaptic inhibition, sniffing
Citation: Moran AK, Eiting TP and Wachowiak M (2021) Circuit Contributions to Sensory-Driven Glutamatergic Drive of Olfactory Bulb Mitral and Tufted Cells During Odorant Inhalation. Front. Neural Circuits 15:779056. doi: 10.3389/fncir.2021.779056
Received: 17 September 2021; Accepted: 06 October 2021;
Published: 27 October 2021.
Edited by:
Mariano Soiza-Reilly, CONICET Institute of Physiology, Molecular Biology and Neurosciences (IFIBYNE), ArgentinaReviewed by:
Claire Cheetham, University of Pittsburgh, United StatesRicardo C. Araneda, University of Maryland, College Park, United States
Copyright © 2021 Moran, Eiting and Wachowiak. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Matt Wachowiak, bWF0dC53YWNob3dpYWtAdXRhaC5lZHU=