Fast and slow feedforward inhibitory circuits for cortical odor processing
Abstract
Feedforward inhibitory circuits are key contributors to the complex interplay between excitation and inhibition in the brain. Little is known about the function of feedforward inhibition in the primary olfactory (piriform) cortex. Using in vivo two-photon-targeted patch clamping and calcium imaging in mice, we find that odors evoke strong excitation in two classes of interneurons – neurogliaform (NG) cells and horizontal (HZ) cells – that provide feedforward inhibition in layer 1 of the piriform cortex. NG cells fire much earlier than HZ cells following odor onset, a difference that can be attributed to the faster odor-driven excitatory synaptic drive that NG cells receive from the olfactory bulb. As a result, NG cells strongly but transiently inhibit odor-evoked excitation in layer 2 principal cells, whereas HZ cells provide more diffuse and prolonged feedforward inhibition. Our findings reveal unexpected complexity in the operation of inhibition in the piriform cortex.
Editor's evaluation
Feedforward inhibition (FFI) typically exerts a powerful effect shaping neural activity. In this paper, Suzuki et al., use a combination of in vivo and in vitro experiments to characterize, for the first time, responses in the two main classes of FFIs in the mouse olfactory cortex, neurogliaform cells (NG) and horizontal cells (HZ). They find that these two cell types have different responses and different connectivity, which partially explains their different responses. This paper also helps resolve a previously perplexing result from a recent publication proposing that FFI in the mouse olfactory cortex plays a negligible role in shaping cortical odor responses, presumably because those authors were only recording from HZ, but not NG, cells.
https://doi.org/10.7554/eLife.73406.sa0Introduction
Synaptic inhibition in the cortex is an intricate process with many interacting parts. Cortical interneurons are classified into dozens of subtypes with distinctive morphological, transcriptomic, and electrophysiological properties (Kawaguchi and Kubota, 1997; Tremblay et al., 2016; Zeng, 2017; Gouwens et al., 2020). These diverse interneurons are embedded in at least three types of cortical inhibitory circuits – feedforward inhibition, feedback inhibition, and disinhibition (Kullmann and Lamsa, 2007; Isaacson and Scanziani, 2011; Denève and Machens, 2016; Tremblay et al., 2016; Lourenço et al., 2020). Within these three categories there exist numerous variants that implement a variety of neural computations. For example, feedforward and feedback inhibition exert different effects depending on whether the effector interneurons innervate the perisomatic or distal dendritic regions of their targets (Bruno and Sakmann, 2006; Isaacson and Scanziani, 2011; Pardi et al., 2020), and disinhibition can depend on the subtypes of interneurons involved and the origin of their inputs (Pi et al., 2013; Fu et al., 2014; Tremblay et al., 2016). In this paper we focus on feedforward inhibition in the primary olfactory (piriform) cortex, and report that a surprising complexity of inhibitory processing is also a feature of this type of inhibition in the paleocortex.
Feedforward inhibition is conventionally understood as an input-tracking mechanism that does not depend on local circuit activity (Tremblay et al., 2016). In the primary neocortex, feedforward inhibition is often mediated by parvalbumin-expressing fast-spiking basket cells localized in the deeper layers that receive direct afferent input from first-order thalamic nuclei (Beierlein et al., 2003; Gabernet et al., 2005; Cruikshank et al., 2007). In this role feedforward inhibition improves sensory discrimination by favoring coincidence detection, but in addition it can provide gain modulation through input normalization (Bruno and Sakmann, 2006; Carandini and Heeger, 2011). More recently it has been reported that feedforward inhibition is also prominent in layer 1 of the neocortex, where it can engage afferent input from other cortical regions and from higher-order thalamic nuclei (Pardi et al., 2020; Anastasiades et al., 2021). However, the properties and functions of layer 1 neocortical interneurons have received relatively little attention and remain enigmatic (Schuman et al., 2019; Fan et al., 2020; Lourenço et al., 2020).
The piriform cortex (PCx) is unusual in that it receives all of its primary afferent input from the olfactory bulb into the superficial part of layer 1 (layer 1a; Neville and Haberly, 2004; Bekkers and Suzuki, 2013). Layer 1a mostly contains the apical dendrites of glutamatergic principal cells located in layers 2 and 3, as well as the somas, axons, and dendrites of two classes of GABAergic interneurons, namely, neurogliaform (NG) cells and horizontal (HZ) cells (Suzuki and Bekkers, 2010a). NG and HZ cells in layer 1a are thus ideally situated to provide feedforward inhibition onto the distal dendrites of layer 2/3 principal cells. Classic work (Biedenbach and Stevens, 1969; Haberly, 1973) as well as more recent studies (e.g. Luna and Schoppa, 2008; Stokes and Isaacson, 2010; Suzuki and Bekkers, 2012; Sheridan et al., 2014; Stokes et al., 2014; Large et al., 2016a; Large et al., 2016b) have used electrical stimulation to explore some of the basic properties of feedforward inhibition in the PCx. However, less is known about how particular types of interneurons respond to odor stimulation. The existence and importance of odor-evoked inhibition in the PCx has been established (Poo and Isaacson, 2009; Zhan and Luo, 2010; Franks et al., 2011; Poo and Isaacson, 2011; Sturgill and Isaacson, 2015; Bolding and Franks, 2017; Tantirigama et al., 2017; Bolding and Franks, 2018) but the identities and properties of neurons responsible for the different kinds of inhibition are, for the most part, uncertain.
Here, we approach this question by using two-photon-targeted whole-cell patch clamping and functional Ca2+ imaging to record from visually identified NG and HZ cells in layer 1a of the PCx in vivo. We find that the odor-evoked feedforward inhibition provided by these two types of interneurons is strikingly different: NG cells express a powerful and transient inhibition that begins rapidly after odor onset, whereas HZ cells provide a more diffuse and delayed form of feedforward inhibition. Thus, two distinctive feedforward inhibitory circuits exist in the superficial layer of the PCx, where they are well placed to generate dynamically complex patterns of inhibition in the distal dendrites of principal cells.
Results
Interneurons that provide feedforward inhibition in the PCx can be targeted ‘in vivo’
In this study we took advantage of the simple architecture of the PCx to examine in isolation only those neural circuits that generate feedforward synaptic inhibition. Afferent input from the olfactory bulb via the lateral olfactory tract (LOT) to the PCx is strictly confined to layer 1a (L1a) of the PCx (Figure 1A; Neville and Haberly, 2004). Hence, only those interneurons with dendrites concentrated in L1a, where they are able to intercept axons from the LOT, are able to generate feedforward inhibition (red cells, Figure 1A). Two such interneuron types have been identified in the PCx, HZ cells and L1a NG cells (Suzuki and Bekkers, 2010a). Conversely, feedback inhibition is largely mediated by interneurons with dendrites concentrated in deeper associational layers, two prominent types being fast-spiking (FS) and regular-spiking (RS) interneurons (blue cells, Figure 1A). In this paper we focus on the feedforward inhibitory neurons, HZ cells and L1a NG cells (here called ‘NG cells’ for short), located in the anterior PCx.
Pair recordings in slices show that all types of interneurons in the anterior PCx synaptically inhibit principal cells (Figure 1A; Suzuki and Bekkers, 2010a; Suzuki and Bekkers, 2012). Inhibition is commonly observed at the soma of principal cells in response to odors in vivo (e.g. Poo and Isaacson, 2009; Zhan and Luo, 2010; Poo and Isaacson, 2011; Sturgill and Isaacson, 2015; Bolding and Franks, 2017; Tantirigama et al., 2017; Bolding and Franks, 2018), although it is not immediately clear whether this inhibition is dominated by feedforward or feedback inhibitory circuits. We decided to explore feedforward inhibition in this study because the relevant interneurons (NG and HZ cells) are close to the cortical surface and more readily accessible to two-photon-targeted patch clamping in vivo.
In GAD67-GFP (Δneo) mice, NG and HZ cells can readily be distinguished under the two-photon microscope: HZ cells are only found close to the LOT whereas NG cells are distributed throughout L1a; HZ cell somata have an elongated shape whereas NG somata are spherical; and HZ cells express much lower levels of GFP than NG cells (Figure 1B, top; Suzuki and Bekkers, 2010a; Suzuki and Bekkers, 2010b). After whole-cell recording, their identities could be confirmed by their dendritic morphology (NG: short, thin, highly branched dendrites; HZ: longer, less branched, often spiny; Figure 1B, bottom) and distinctive intrinsic electrical properties (Figure 1C and D), similar to in vitro (Suzuki and Bekkers, 2010a). Thus, we were confident that we could record from identified NG and HZ cells in vivo.
NG and HZ cells respond strongly to odors
We applied up to 15 odorants from a variety of chemical functional groups to urethane-anesthetized mice and measured the voltage responses of identified NG and HZ cells in the anterior PCx (Figure 1E and F). Three features were apparent in these responses: first, both cell types responded to many different odors (i.e. they were broadly tuned); second, regular oscillations in the membrane potential (Vm) often became larger in the presence of odor; and third, HZ cells appeared to respond more slowly to odors than did NG cells. Each of these features is examined in the following sections.
NG and HZ cells are broadly excited by odors
Our palette of odorants was drawn from a variety of chemical functional groups intended to span a large part of ‘odor space’ (although this term is difficult to define; Pashkovski et al., 2020; Ravia et al., 2020). We measured odor responsiveness by median-filtering each trace to remove action potentials (APs) then testing whether the z-scored Vm amplitude exceeded a positive threshold during the period of odor application (Experimental procedures; Figure 1—figure supplement 1). These data were analyzed in two different ways: (i) by calculating the fraction of odors that each cell responded to, then averaging across cells (‘cell-averaged index’); and (ii) by calculating the fraction of cells each odor activated, then averaging across odors (‘odor-averaged index’). NG and HZ cells both had large cell-averaged indices that were not significantly different from each other (NG: 0.78 ± 0.04, n = 28 cells from 19 mice; HZ: 0.66 ± 0.12, n = 9 cells from 9 mice; p = 0.34, Welch’s two-sample unpaired t-test). The odor-averaged index was also large for both cell types but was significantly smaller for HZ cells (NG: 0.75 ± 0.03, n = 15 odors; HZ: 0.60 ± 0.06, n = 15; p = 0.025, Welch’s two-sample unpaired t-test). Similar results were found by using an alternative pair of measures, the lifetime and population sparseness, which are analogous to the cell- and odor-averaged indices, respectively (Willmore and Tolhurst, 2001; Poo and Isaacson, 2009).
Respiration-locked oscillations in Vm are prominent in NG and HZ cells
Subthreshold oscillations in Vm were frequently seen in both NG and HZ cells (e.g. Figure 2A and B, top, show the same cells as in Figure 1E and F responding to ethyl-n-butyrate; gray traces are respiration, with upward transients indicating onset of exhalation). Expanding the traces in windows before (b), during (d), and after (a) odor application revealed that oscillations in Vm were synchronized to respiration and, at least in the NG cell, the amplitude of Vm oscillations appeared to increase during the odor (Figure 2A and B, bottom). These observations were quantified by excising the segments of Vm that lay between successive positive peaks of the respiration trace, linearly warping them to have the same time axis, then averaging together all such segments within each of the three windows (b, d, a) for each odor. The peak-to-peak amplitude of the average Vm and the time of the positive peak of average Vm, expressed as a phase of the respiration cycle, were plotted for each odor application (Figure 2C; each triplet of connected points is from each odor; only the data from the NG cell in Figure 2A are shown). For this particular NG cell, the peak-to-peak Vm amplitude increased significantly during the odor (b: 7.7 ± 0.8 mV; d: 11.1 ± 1.2 mV; n = 12, p = 0.002, one-way ANOVA with Tukey’s contrasts), whereas the phase of the peak was unchanged (b: 0.49 ± 0.02; d: 0.45 ± 0.03, both expressed as a fraction of the respiration cycle; p = 0.47).
A similar analysis was done for all NG and HZ cells in our dataset and the summary is shown in Figure 2D. The mean peak-to-peak Vm amplitude increased significantly during odors in NG cells but not in HZ cells (Figure 2D, left; n = 232 trials in 22 NG cells, p < 0.0002; n = 98 trials in 9 HZ cells, p = 0.78; one-way ANOVA with Tukey). In contrast, the phase of the peak Vm was unchanged by odors in both NG and HZ cells (Figure 2D, right). The same result was obtained if we used a different measure of phase-locking between Vm and respiration, the cross-covariance (Figure 2—figure supplement 1). Thus, while NG and HZ cells both show strong respiration-coupled oscillations in their subthreshold Vm, odors affected these cells differently: in NG cells the oscillation amplitude increased, whereas in HZ cells Vm tended to depolarize without a change in the oscillation.
We wondered if the oscillations in Vm we found in NG (and HZ) cells in the PCx were also seen in NG cells in other cortical areas and, if so, whether they were synchronized to respiration. We made targeted whole-cell recordings from NG cells in the upper layers of primary somatosensory cortex in vivo and found that, although Vm oscillations were strongly present, their temporal structure was different from those in PCx and they were not synchronized to respiration (Figure 2—figure supplement 2; see also Suzuki et al., 2014).
NG and HZ cells tend to fire early and late, respectively, following odor onset
Next, we turned to the apparent difference in the kinetics of the odor response, with NG cells appearing to be excited more quickly following odor application (Figure 1E and F). We confirmed this impression by constructing AP raster plots and peristimulus time histograms (Figure 3A): spiking in NG cells reached a peak at 0.43 ± 0.14 s after odor onset (n = 7 cells), whereas in HZ cells the peak occurred at 1.92 ± 0.24 s (n = 5 cells; significantly different, p = 0.0013, Welch’s two-sample unpaired t-test). Control experiments confirmed that none of these responses were limited by the rate of delivery of odorants by our olfactometer (Figure 3—figure supplement 1).
To test whether this difference in dynamics could be observed using a less invasive approach, we turned to two-photon Ca2+ imaging with the red-shifted indicator Cal-590. NG and HZ cells could be distinguished as before from their soma location, soma shape and GFP fluorescence (Figure 3B, top), and odor-evoked spiking could be resolved from changes in ΔF/F0 (Figure 3B, bottom). Consistent with whole-cell patch clamping (Figure 3A), Ca2+ imaging showed that NG cells fired quickly after odor onset and HZ cells fired with a delay (Figure 3B; traces are averages of n = 6 NG cells and n = 5 HZ cells; individual traces were too noisy to reliably determine the average time to peak).
We also took advantage of the Ca2+ imaging approach to examine the effect of anesthetics on this neural circuit. (Because the surgery to expose the PCx is so invasive, our experiments could not use awake animals.) All of the above experiments used urethane (0.7 g/kg). We repeated the imaging experiment in Figure 3B using fentanyl plus medetomidine, which has been used to induce a more awake-like state of anesthesia (Constantinople and Bruno, 2011; Altwegg-Boussac et al., 2014). The result was the same (Figure 3—figure supplement 2): in response to odors, NG cells fired early, HZ cells late. These findings, together with a report that urethane at the relatively low concentration we used here has little effect on ligand-gated synaptic receptors (Hara and Harris, 2002), suggest that the different odor-response dynamics of NG and HZ cells are not related to anesthesia.
Different kinetics of odor-evoked EPSPs in NG and HZ cells
Having found that NG and HZ cells differ in their AP responses to odors, we next asked whether similar differences could be observed in the underlying EPSPs. Odor-evoked EPSPs were median-filtered to remove APs and notch-filtered at 2–4 Hz to remove respiration-coupled oscillations, then all EPSPs measured in the same cell for different odors were averaged together (responses to 4–15 odors per cell; Figure 4—figure supplement 1). Because these averages were still relatively noisy, a smooth curve (Materials and methods) was fitted to each averaged EPSP and the key parameters were measured from this curve (Figure 4—figure supplement 1). The analysis showed that odor-evoked EPSPs reached their peak about twofold faster in NG cells than in HZ cells (Figure 4A; NG: 0.53 ± 0.05 s, n = 16 cells; HZ: 1.26 ± 0.26 s, n = 6 cells; p = 0.038, Welch’s two-sample unpaired t-test). Thus, the slower-rising EPSPs in HZ cells (Figure 4A) can explain the longer delay to firing of these cells (Figure 3).
To confirm this result, we used a different method that focused on the rising phase of individual, unaveraged EPSPs. Traces were again median-filtered to remove APs, then the z-scored mean Vm for each respiration cycle was plotted versus time and a straight line was fitted over the period of odor application (Figure 4B; typical NG cell on left, HZ cell on right). The slope of this line tended to be negative for NG cells and positive for HZ cells (red and blue line, respectively; Figure 4B, bottom). These data were quantified for all odors and cells by plotting the peak z-scored Vm amplitude during the odor period versus the fitted slope (Figure 4C, top; NG cells in red, n = 310 trials; HZ cells in blue, n = 188 trials). Trials that did not give an odor response (i.e. with peak z-scored Vm values ≤ 2.5) are grayed out (Figure 4C, top). Histograms of all the remaining trials confirm that the EPSP slopes for NG and HZ cells were significantly different (Figure 4C, bottom; NG slope, –0.27 ± 0.02, n = 232 trials above the z-score threshold; HZ slope, 0.21 ± 0.04, n = 84 trials above threshold; p < 0.001, Kolmogorov-Smirnov [KS] test). Thus, during odors, EPSPs in NG cells tend to decline from an early peak while EPSPs in HZ cells tend to rise to a later peak. These behaviors are consistent with the odor-evoked AP firing observed in NG and HZ cells (Figure 3).
Differences in excitatory synaptic input can explain the odor response differences between NG and HZ cells
What cellular mechanisms might explain the different odor response dynamics of NG and HZ cells? Given the slow rise times of odor-evoked EPSPs in these cells (hundreds of milliseconds, much slower than electrically evoked EPSPs measured in slices; Suzuki and Bekkers, 2010a; Suzuki and Bekkers, 2012), we looked for ways in which slowly rising compound EPSPs in vivo could be constructed from trains of unitary EPSPs. Perhaps the simplest explanation is that NG cells receive depressing EPSPs from the olfactory bulb via the LOT, whereas HZ cells receive facilitating EPSPs (Figure 5A). To test this idea, we first recorded from NG and HZ cells in slices and applied patterned electrical stimuli to the LOT while pharmacologically blocking GABAA receptors. For brief trains of stimuli (6 pulses at 40 Hz, modeled on in vivo patterns; Suzuki and Bekkers, 2006), the response was the exact opposite of what was required: NG cells facilitated and HZ cells depressed (Figure 5B, top). We next tried a more realistic model of odor-evoked afferent excitation, that is, a train of 40 Hz stimuli repeating at 3 Hz (the approximate respiration frequency in mice; Poo and Isaacson, 2009). This stimulus produced a modestly depressing envelope of EPSPs in both NG and HZ cells (Figure 5B, bottom). Thus, the short-term dynamics of the LOT afferents could explain the declining response during odors in NG cells but not the facilitating response in HZ cells.
To further explore this question, we returned to in vivo recordings but now used whole-cell voltage clamp to isolate odor-evoked EPSCs at a holding potential of –70 mV (Figure 5C). Similar to the result for EPSPs (previous section), the time between odor onset and EPSC peak was about twofold faster in NG cells than in HZ cells (NG: 0.72 ± 0.09 s, n = 30 cell-odor pairs; HZ: 1.24 ± 0.12 s, n = 21 cell-odor pairs; p = 0.0014, Welch’s two-sample unpaired t-test; Figure 5C inset, bottom). Thus, given that the membrane time constant is fast (~4 ms) and not significantly different between NG and HZ cells (Suzuki and Bekkers, 2010a), the slower EPSC rise time in HZ cells is consistent with the slower-rising EPSP (Figure 4A) and the delayed onset of odor-evoked spiking (Figure 3) in HZ cells.
In summary, despite the in vitro finding that HZ cells do not receive facilitating EPSPs from the LOT (Figure 5B), both current clamp (Figure 4A) and voltage clamp (Figure 5C) experiments in vivo show that HZ cells receive slower-rising synaptic inputs than NG cells. Thus, the dynamics of excitatory synaptic inputs at least partially determine the delayed odor responses of HZ cells, although the origin of these slower dynamics remains unclear (see Discussion).
Synaptic inhibition also contributes to the odor responses of NG and HZ cells
Although we have so far focused on excitatory synaptic inputs, a likely contribution of inhibitory inputs cannot be excluded. NG and HZ cells may engage in mutual inhibition, thereby modifying the odor-evoked excitation both receive from the LOT. We addressed this possibility by making in vivo whole-cell voltage clamp recordings to look for odor-evoked IPSCs in isolation at a holding potential of +10 mV (Poo and Isaacson, 2009). Such IPSCs were, indeed, present in both HZ and NG cells (Figure 6A). Their properties were not significantly different between the two cell types (mean amplitude: 50.3 ± 4.3 pA versus 49.9 ± 11.0 pA; time to peak: 0.58 ± 0.12 s versus 0.35 ± 0.05 s; data for NG and HZ cells, respectively; all p > 0.1, n = 14 and 5 cell-odor pairs, respectively). Current clamp experiments confirmed that odor-evoked IPSPs could also be observed in HZ cells when they were depolarized by somatic current injection (Figure 6—figure supplement 1).
What might be the origin of these inhibitory inputs? We addressed this question by making pair recordings in slices. These experiments revealed that NG→HZ cell connections were frequent and strong (Figure 6B; 14/43 = 32.6% connectivity; mean connection conductance 0.87 ± 0.41 nS, range 0.10–5.86 nS, n = 14 pairs) whereas the reverse connection was never seen (0/13 = 0%; significantly different, p = 0.02, chi-square 2 × 2 contingency test). NG→NG cell connectivity (2/17 = 11.8%) was not significantly different from NG→HZ connectivity (p = 0.1, chi-square 2 × 2 contingency test; Figure 6B, right) although the n-values are modest. In contrast, connectivity between HZ cells was weaker (1/19 = 5.3%; Figure 6B, right). To sum up, these connectivity experiments suggest that odor-evoked synaptic inhibition of both NG and HZ cells most likely comes from NG cells after they are excited by the LOT.
If NG cells are responsible for much of the inhibition of layer 1a interneurons, do NG cells also receive privileged input from the LOT? We addressed this question by making dual whole-cell recordings from an NG and HZ cell in a slice while applying a train of minimal electrical stimuli (five at 20 Hz) to the LOT while pharmacologically blocking GABAA receptors (Figure 6C). In some of these experiments (4 out of 7 pairs), a clear plateau region was seen in a plot of EPSC amplitude versus stimulus number as the stimulus strength was progressively increased in small steps (Figure 6C, bottom left, stimulus numbers 20–45); this is suggestive of a unitary input (Stokes and Isaacson, 2010). Notably, the EPSC in the NG cell was always much smaller than that in the HZ cell for the first EPSC in the train (e.g. Figure 6C, bottom left and top right). This result did not depend on which cell was closer to the LOT or stimulating electrode (data not shown). A similar effect was seen in the three experiments in which a plateau was not observed (e.g. Figure 6—figure supplement 2).
Interestingly, because of the characteristic short-term facilitation and depression of LOT inputs to NG and HZ cells, respectively (Suzuki and Bekkers, 2010a), the imbalance in the strength of LOT inputs to these two cells changed during the train. This effect was quantified by calculating the ratio of the EPSC amplitude in the NG cell to that in the HZ cell for each of the five EPSCs in the train (Figure 6C, bottom right). For the first EPSC, the mean NG:HZ amplitude ratio was 0.27 ± 0.03 (mean ± standard error of the mean [SEM], n = 7 cell pairs), but by the third EPSC this ratio had increased to 1.38 ± 0.32 (significantly different, p = 0.0002, n = 7, one-way ANOVA with Tukey; Figure 6C, bottom right). These results show that NG cells do not receive privileged input from the LOT compared to HZ cells; rather, the reverse is the case for the first EPSC in a train. Later in the train, however, the inputs become more equal. Given the known dynamics of afferent input from the olfactory bulb (i.e. modulated bursts of APs; Cang and Isaacson, 2003; Margrie and Schaefer, 2003), these results suggest that NG and HZ cells mostly operate in a regime where they receive similar afferent excitation.
Feedforward inhibition alters the synaptic responsiveness of SP cells
Lastly, we explored the functional consequences of HZ and NG cell inhibition for one of their major targets, layer 2 superficial pyramidal (SP) cells. Because it is difficult to disambiguate these two types of feedforward inhibition in vivo, we conducted the experiments in slices.
We began by eliciting in vivo-like IPSPs in SP cells. Extracellular stimuli were applied to layer 1a in patterns obtained from the in vivo odor-evoked firing patterns of NG and HZ cells (Figure 7A, red traces; gray bars labelled ‘odor period’ represent the period during which the odor was applied in the in vivo experiments). Averaged postsynaptic IPSPs were recorded in SP cells while pharmacologically blocking ionotropic glutamate receptors (Figure 7A, black traces; average of n = 44 or 51 single traces for 13 or 3 different NG or HZ cell stimulus patterns, respectively, while recording from 6 different SP cells).
Two further manipulations were done to make the recordings more in vivo-like. First, NG stimulus recordings were made distant (>400 µm) from the LOT to avoid stimulating HZ cell axons, which are clustered around the LOT (Suzuki and Bekkers, 2010a). On the other hand, HZ stimulus recordings were made near the LOT, where a mixture of HZ and NG cell axons were likely excited. Second, we warped the stimulus patterns so the respiration trace recorded for each stimulus pattern matched a reference respiration trace (shown in gray, Figure 7A, bottom; see Materials and methods; a similar approach was used in Figure 2). This warping of the time base was done to preserve any respiration-synchronized structure in the stimulus patterns (as in Figure 2) when averaging across different patterns. Such synchronization is apparent in the respiration-locked oscillations in the averaged IPSPs (Figure 7A, black traces). These results show that NG cells generate a large, rapid IPSP in SP cells, whereas HZ cells generate a smaller and more diffuse IPSP that persists beyond the end of odor application.
In a final series of experiments we examined the effect of these two types of synaptic inhibition on spiking patterns in postsynaptic SP cells. The method described in the previous paragraph was employed, except that a dynamic clamp was used to replay into the SP cell an odor-evoked excitatory postsynaptic conductance (EPSG) that had previously been recorded from an SP cell in vivo (Figure 7—source code 1). Again, the EPSG and all stimulus patterns were warped to match their respiration traces. When injecting the EPSG alone, an in vivo-like train of APs was evoked in the SP cell (Figure 7B, black traces). When patterned extracellular stimuli were applied at the same time, eliciting IPSPs, fewer APs were evoked in the SP cell by the EPSG (Figure 7B, blue traces; spikes that dropped out are indicated by red asterisks above the black traces). These effects were quantified by calculating normalized PSTHs for APs in the SP cell (Figure 7B, bottom; black, EPSG alone; blue, EPSG plus IPSPs; averages from n = 7 SP cells with NG stimulus patterns, or n = 4 SP cells with HZ stimulus patterns). The results confirm that NG cells strongly suppress early odor-evoked firing in SP cells (p < 0.001, KS test), whereas HZ cells tend to have a weaker, delayed effect (here, not significant; p = 0.81, KS test).
Discussion
In this paper we used whole-cell patch clamping and two-photon Ca2+ imaging in vivo to characterize the odor responses of two types of GABAergic interneurons that provide feedforward inhibition in the input layer (layer 1a) of the anterior PCx. We find that both NG cells and HZ cells are broadly excited by different odors, but the time to reach peak excitation is much slower in HZ cells than in NG cells. This difference can be explained by a twofold slower time to peak of the odor-evoked compound EPSC in HZ cells compared to NG cells. Synaptic inhibition of HZ cells by NG cells may also help to suppress early firing of HZ cells, and suggests how lateral inhibition and feedforward inhibition may interact in the same circuit. In addition, these two cell types differ in their oscillatory response to odors, with NG cells showing larger-amplitude respiration-coupled EPSPs during odor sampling. Lastly, we show that NG and HZ cells have distinctive effects on the excitability of downstream pyramidal neurons: NG cells generate powerful inhibition immediately after odor onset, whereas HZ cells exert more diffuse and prolonged inhibition. Our findings reveal different types of inhibitory responses at the first stage of cortical odor processing, and add to a growing understanding of the role of afferent-driven feedforward inhibition in cortical processing more broadly (Khubieh et al., 2016; Fan et al., 2020; Anastasiades et al., 2021).
Use of anesthesia
Two-photon-targeted patch clamping and Ca2+ imaging as implemented in this study required direct access to the surface of the PCx (Margrie et al., 2003; Tantirigama et al., 2017). Although similar experiments can be done during wakefulness and semi-paralysis of mice (Pashkovski et al., 2020), our animal ethics protocols required us to conduct all surgery and experiments under general anesthesia. We used urethane at the minimum concentration empirically determined to be effective in providing stable anesthesia with complete abolition of reflexes (0.7 g/kg s.c.). Urethane has been widely used in previous in vivo studies of the PCx (e.g. Barnes et al., 2008; Poo and Isaacson, 2009; Poo and Isaacson, 2011; Wesson et al., 2011; Chapuis and Wilson, 2011; Sturgill and Isaacson, 2015), and we have confirmed that urethane at this concentration has no effect on electrical activity in the PCx compared with mice anesthetized with fentanyl/medetomidine (which induces a more awake-like brain state; Figure 3—figure supplement 2; Constantinople and Bruno, 2011; Tantirigama et al., 2017). We avoided using ketamine/xylazine anesthesia which has been reported to alter PCx activity through its blockade of NMDA receptors (Fontanini and Bower, 2005; Tantirigama et al., 2017). Despite these considerations, however, it remains possible that our findings are affected by urethane and further study is warranted.
Comparison with previous findings
Most previous work has used electrical stimulation in slices to study feedforward inhibition in the PCx (Luna and Schoppa, 2008; Stokes and Isaacson, 2010; Suzuki and Bekkers, 2010a; Suzuki and Bekkers, 2012; Sheridan et al., 2014; Stokes et al., 2014; Large et al., 2016a; Large et al., 2016b). Only two reports, to our knowledge, have explicitly examined this feedforward circuit using odor stimulation in vivo. Poo and Isaacson, 2009, made blind whole-cell recordings from neurons in layer 1, while Bolding and Franks, 2018, made unit recordings from optogenetically identified GABAergic neurons in layer 1. However, neither study distinguished interneuron subtypes.
We found broad tuning of odor-evoked EPSPs in NG and HZ cells, consistent with previous findings (Poo and Isaacson, 2009; odor-averaged index, 0.50 ± 0.04 cf 0.75 and 0.60 for NG and HZ cells, respectively, found here). Interestingly, broad odor tuning has also been reported for some (Zhan and Luo, 2010; Poo and Isaacson, 2011; Bolding and Franks, 2017) but not all (Sturgill and Isaacson, 2015) classes of feedback inhibitory neurons located in deeper layers of the PCx. It is likely that this diversity in stimulus tuning reflects the number and variety of functional inputs from upstream neurons. Indeed, it has been shown that layer 1 interneurons receive a higher convergence of afferent input from the olfactory bulb than do principal neurons (Poo and Isaacson, 2009; Miyamichi et al., 2011). Future work could repeat these experiments while distinguishing NG and HZ cells. Differences might arise because HZ cells, unlike NG cells, are clustered around the LOT where they may encounter a higher density of afferents (Suzuki and Bekkers, 2010a).
We often observed oscillations in subthreshold Vm that were phase-locked to respiration, consistent with previous reports that such oscillations are ubiquitous in the olfactory system (Fontanini and Bower, 2006; Kay et al., 2009; Poo and Isaacson, 2009; Wilson, 2010; Oswald and Urban, 2012; Kay, 2014; Jiang et al., 2017). Our novel finding was that respiratory oscillations in Vm often increased in amplitude during odor sampling, but only in NG cells. These odor-dependent oscillations were functionally relevant because they generated respiratory-patterned IPSPs in target neurons (Figure 7A) and provided rhythmic inhibition of similarly patterned EPSPs (Figure 7B; see also Poo and Isaacson, 2009). This difference between NG and HZ cells suggests that NG cells provide fast, phase-critical feedforward inhibition whereas HZ cells provide a slower, tonic form of inhibition. Interestingly, Bolding and Franks, 2018, reported only the latter kind of feedforward inhibition in multi-unit recordings from awake mice.
Difference in odor-evoked EPSP kinetics
Our main finding is that odor-evoked firing in HZ cells has a delayed onset, contrasting with the rapid onset in NG cells. We showed that this effect can be at least partly explained by the slower time-to-peak of the odor-evoked compound EPSC in HZ cells (Figure 5C). Surprisingly, however, these kinetic differences cannot be explained by the properties of short-term synaptic plasticity at LOT synapses, assayed in slices (Figure 5B). What might be alternative explanations?
The LOT is a heterogeneous fiber tract that contains the axons of two distinct types of projection neurons in the olfactory bulb, mitral and tufted cells. Mitral cells have been reported to respond about twofold more slowly to odor stimulation than tufted cells (Fukunaga et al., 2012; Igarashi et al., 2012) because of differences in synaptic inhibition in the bulb (Fukunaga et al., 2012; Geramita and Urban, 2017). An appealing possibility is that mitral and tufted cell axons preferentially target HZ and NG cells, respectively; that is, the delayed odor response of HZ cells may originate in the olfactory bulb rather than at synaptic terminals in the PCx. A difficulty with this idea is that the axons of tufted cells in the dorsal olfactory bulb do not extend far from the LOT and only reach about halfway along the length of the anterior PCx, that is, to ~1 mm anterior to Bregma (Paxinos and Franklin, 2001; Igarashi et al., 2012). Our in vivo recordings were made close to the LOT but more caudally, around 0.6 mm anterior to Bregma. However, it remains possible that tufted cells in other parts of the olfactory bulb (not just in the dorsal bulb) have a more caudal projection. This hypothesis could be tested in future work.
A second possible explanation is that HZ cells may receive delayed bulbar input via an intermediary, such as the anterior olfactory nucleus (AON; Haberly and Price, 1978; McGinley and Westbrook, 2011; Kay and Brunjes, 2014). Although AON axons are sparser in layer 1a than in deeper layers, they are consistently seen there (Haberly and Price, 1978; Hagiwara et al., 2012; Russo et al., 2020). This hypothesis could be tested anatomically or by silencing the AON during odor application.
A third possibility is that lateral inhibition from NG cells counteracts early excitation in HZ cells, delaying depolarization. In support of this idea, NG cells profusely innervate HZ cells (Figure 6B) and an odor-evoked hyperpolarization with the right time course can be observed in HZ cells (Figure 6—figure supplement 1). In addition, NG cells receive facilitating excitatory input from the LOT, driving them at least as strongly as HZ cells during the bursts of spikes that typify odor-evoked output from the olfactory bulb (Figure 6C; Cang and Isaacson, 2003; Margrie and Schaefer, 2003). On the other hand, NG cells also receive odor-evoked synaptic inhibition, presumably from other NG cells (Figure 6A), yet do not exhibit delayed odor responses like HZ cells. It is plausible that inhibition from NG cells is larger in HZ cells but it was not apparent here because the number of experiments was relatively small.
It should be kept in mind that both NG and HZ cells may also receive feedback inhibition from interneurons located in deeper layers, particularly from layer 3 somatostatin-positive RS interneurons with axons that ramify in L1 (Suzuki and Bekkers, 2010a). HZ cells can extend their dendrites into L1b (Suzuki and Bekkers, 2010a) so might be especially susceptible to feedback inhibition of this kind. Thus, depending on patterns of activity in afferent and associational circuits, feedforward and feedback inhibition may interact.
Functional significance of feedforward inhibition in the PCx
Feedforward inhibition is generally recognized as an input-tracking mechanism which, in the hippocampus and neocortex, can synchronize spike timing (Pouille and Scanziani, 2001; Gabernet et al., 2005) and modulate gain through input normalization (Pouille et al., 2009). The PCx is unusual in that it receives its feedforward inhibitory input onto the distal apical dendrites of principal cells, rather than onto perisomatic regions as in other cortices (Lourenço et al., 2020). What are the consequences of this atypical architecture for the role of feedforward inhibition in the PCx?
One slice study has shown that feedforward inhibition works together with feedback inhibition to provide spatially and temporally modulated inhibition during bursts of olfactory input (Stokes and Isaacson, 2010). Intriguingly, the shift in inhibition reported by Stokes and Isaacson, 2010, is in the opposite direction to that found in other brain regions (Pouille and Scanziani, 2004; Silberberg and Markram, 2007). These results suggest that feedforward inhibition in the PCx is functionally important but has unusual features.
More recently, an in vivo unit recording study in awake mice showed that odor-evoked spiking in feedforward inhibitory neurons in the PCx increases only slowly and weakly after inhalation, suggesting that these neurons provide tonic inhibition and do not play a major role in sculpting odor-evoked responses (Bolding and Franks, 2018). Modeling confirmed that feedforward inhibition provides modest subtractive normalization in the PCx, in contrast to the divisive normalization it provides in other cortical regions (Stern et al., 2018). An intuitive explanation is that PCx uses a temporal or rate-based code that is more susceptible to ‘recurrent normalization’ provided by feedback inhibition (Sturgill and Isaacson, 2015; Bolding and Franks, 2018; Stern et al., 2018; Bolding et al., 2020; Pashkovski et al., 2020). Hence, the evidence of Bolding and colleagues suggests that feedforward inhibition plays a minor role in odor processing under the conditions of their experiments.
The results of Bolding and colleagues seem at variance with our finding that many layer 1a interneurons, especially NG cells, respond vigorously to odors. The discrepancy may lie in the different recording conditions, for example, whole-cell versus unit recordings, placement of the electrodes (close to or far from the LOT), or our use of anesthesia. Further work is required. However, if feedforward inhibition in the PCx is, indeed, weakly engaged in normal olfactory processing, what other functions might it serve?
The key property of feedforward inhibition in the PCx may be its dendritic localization. By inhibiting the distal apical dendrites of layer 2/3 principal cells, NG and HZ cells have the potential to dampen passive propagation of afferent EPSPs to the soma and to veto dendritic electrogenesis (Larkum et al., 1999; Palmer et al., 2014; Pardi et al., 2020; Anastasiades et al., 2021). Indeed, it has been directly shown in vitro that a single layer 1 interneuron can inhibit Ca2+ signaling in the distal dendrites of PCx principal cells in a branch-specific fashion (Stokes et al., 2014). Feedforward inhibition might then provide a mechanism for regulating processes that involve dendritic electrogenesis and plasticity – including burst-firing (Tseng and Haberly, 1989; Protopapas and Bower, 2001), spike timing-dependent plasticity (Kanter et al., 1996; Johenning et al., 2009; Cassenaer and Laurent, 2012), and NMDA spikes (Kumar et al., 2018; Kumar et al., 2021) – and which may only become apparent during experimental paradigms that engage olfactory learning (Wilson and Stevenson, 2006; Ghosh et al., 2015; Shakhawat et al., 2015; Meissner-Bernard et al., 2019). Future work would need to explore this possibility.
Two types of feedforward inhibition
We have previously shown in slices that layer 1a NG cells and HZ cells generate slow-rising (2–3 ms) and fast-rising (1 ms) feedforward unitary IPSPs, respectively, in layer 2 principal neurons, leading us to suggest that NG cell-mediated feedforward inhibition is slower and more diffuse than that provided by HZ cells (Suzuki and Bekkers, 2012). In the present report, after taking into account the much slower odor-evoked excitatory drive received by HZ cells, we reach the opposite conclusion. Our results are a reminder that the unitary properties of synapses, measured in vitro, can be less important than the patterns of concerted synaptic activity occurring in vivo.
What might be the value of delivering fast (NG) and slow (HZ) feedforward inhibition? We have shown in a slice simulation (Figure 7) that NG cells can strongly inhibit firing in principal cells immediately after a step application of odor. Although this effect was not observed in the study of Bolding and Franks, 2018, there nevertheless appears to be scope for a physiological role for this kind of fast inhibition. For instance, NG cells generate GABA transients that spill out of the synaptic cleft and can activate heterosynaptic GABAA and GABAB receptors (Karayannis et al., 2010; Overstreet-Wadiche and McBain, 2015). Future work could explore ways in which this ‘volume transmission’ could enhance computational complexity in the PCx.
The role of HZ cell-mediated feedforward inhibition is more puzzling. Our slice simulation showed that HZ cells generate a delayed, diffuse response with little effect on principal cell excitability (Figure 7B). Trained rodents can discriminate odors within a single sniff (<100 ms; Uchida and Mainen, 2003; Abraham et al., 2004; Rinberg et al., 2006; Resulaj and Rinberg, 2015). What could be the relevance of late-firing HZ cells in this scenario? Two unique features of HZ cells suggest that they are potentially important. First, HZ cells are only found close to the LOT (<~200 µm), giving them privileged inhibitory influence over a spatially restricted subset of principal cells (Suzuki and Bekkers, 2010a; Suzuki and Bekkers, 2012). Second, they are the only class of interneuron in the PCx, and one of few in the cerebral cortex, that are profusely spiny (Haberly, 1983). Given the importance of spines for synaptic plasticity (Engert and Bonhoeffer, 1999; Maletic-Savatic et al., 1999; Toni et al., 1999; Matsuzaki et al., 2004; Tonnesen et al., 2014), HZ cells are poised to receive learning-dependent synaptic excitation, which in turn may allow them to exert adjustable inhibitory control over neurons in their neighborhood.
In summary, we have shown that two forms of odor-evoked feedforward inhibition, fast and slow, are expressed in layer 1a of the PCx. Our work reveals previously unknown circuit elements in the PCx, and adds to a growing understanding of the role of neurons in layer 1 of the cerebral cortex.
Materials and methods
Animals and surgery for ‘in vivo’ experiments
Request a detailed protocolAll experimental procedures were approved by the Animal Experimentation Ethics Committee of the Australian National University and conform to the Australian Code for the Care and Use of Animals for Scientific Purposes, published by the National Health and Medical Research Council of Australia.
Experiments used heterozygous GAD67-GFP (Δneo) mice, which have the EGFP gene targeted to the Gad1 locus (Tamamaki et al., 2003). The colony was maintained on a C57BL6/J background. Animals were aged 35–75 days and weighed 17–24 g. For surgery, an animal was sedated with chlorprothixene (5 mg/kg i.p.) then injected with urethane (0.7 g/kg s.c.) for general anesthesia plus atropine (0.2–0.3 mg/kg s.c.) to reduce secretions. The correct level of anesthesia was confirmed by observing regular respiration and the complete absence of a pinch reflex; a topup dose of urethane was sometimes required after 5–6 hr. A local anesthetic (prilocaine, 0.2 mg/kg) was applied topically to incision sites during surgery. Access to the PCx was via the cheek and upper mandible, as previously described (Stettler and Axel, 2009; Tantirigama et al., 2017). Briefly, the following surgical procedure was used (total duration ~2 hr). After retracting the skin, superficial blood vessels of the cheek were cauterized, then the temporalis muscle was carefully detached and retracted toward its base near the eye, revealing the temporal aspect of the skull. The zygomatic bone and the upper sections of the mandible, including the coronoid and the condyloid processes, were removed, followed by the upper section of the masseter muscle, exposing the basolateral surface of the skull. At this point, the PCx was visible under the translucent skull, identifiable using anatomical landmarks such as the LOT and the middle cerebral artery (MCA), which are roughly perpendicular to each other. A metal head post was then glued on top of the skull and the head was stabilized. A craniotomy (~2 mm2) was made with a fine drill (Osada Electric, Nagoya, Japan, or Ram Products Inc, Dayton, NJ) just caudal to the MCA and close to the LOT. This placement of the craniotomy ensured that all recordings were made in the anterior PCx at approximately +0.6 mm from Bregma. In many experiments the dura was left intact, but in some the dura was carefully removed using a needle and fine forceps. After completion of the surgery, a small chamber made from a plastic weighboat and dental cement was constructed around the site. To keep the area hydrated and allow immersion of the microscope objective, the chamber was filled with a Ringer’s solution containing (mM) 135 NaCl, 5.4 KCl, 1.8 CaCl2, 1 MgCl2, 5 HEPES at pH 7.4. For all the above procedures, as well as during recordings, the animal was placed on an electrically heated surface at ~37°C and was kept hydrated by periodic s.c. injections of normal saline with 2% dextrose.
Two-photon targeted patch clamping
Request a detailed protocolGFP-positive interneurons in GAD67-GFP (Δneo) mice were visualized using a two-photon MOM microscope (Sutter Instrument Company, Novato, CA) with a 40×/0.8 NA water immersion objective (Olympus, Tokyo, Japan) and a Chameleon Ultra Ti:Sapphire laser (Coherent, Santa Clara, CA) tuned to 800–820 nm (Suzuki et al., 2014). Frames (512 × 512 pixels) were acquired simultaneously through a red and green channel at ~5 Hz under the control of ScanImage (Vidrio Technologies, Ashburn, VA). Patch pipettes were pulled with a longer taper than usual and had resistances of 5–8 MΩ when filled with internal solution comprising (in mM) 135 KMeSO4, 7 NaCl, 0.1 EGTA, 2 Na2ATP, 2 MgCl2, 0.3 GTP, 10 HEPES at pH 7.2, supplemented with 0.2–0.4% biocytin (295–300 mOs/kg). This solution had a Cl- concentration of 11 mM and a measured junction potential of –7 mV. For voltage clamp experiments CsMeSO3 replaced KMeSO4. These solutions also contained Alexa Fluor 594 (20 µM) for visualizing the electrode in the red channel. Patch electrodes were positioned using a micromanipulator (MP-285, Sutter Instrument, Novato, CA) and electrical recordings were obtained with a MultiClamp 700B amplifier (Molecular Devices, San Jose, CA). Data were filtered at 10 kHz and sampled at 20–50 kHz using an Instrutech ITC-18 digitizing interface (HEKA, Ludwigshafen, Germany) under the control of Axograph X (Axograph Scientific, Sydney, Australia). The reference electrode was a Ag/AgCl wire inserted under the skin. The patch electrode was advanced rapidly to penetrate the dura, then more slowly to approach the selected cell and obtain a gigaseal whole-cell recording in the usual way (Margrie et al., 2002). For current clamp recordings, bridge balance and capacitance neutralization were adjusted and the cell was allowed to remain at its resting potential. For voltage clamp recordings, series resistance compensation was not used. Cells were included in the dataset if they had a mean resting potential more hyperpolarized than –50 mV and were stable enough to allow the recording of responses to at least five odors. In addition, cells had to be unambiguously identified as either NG or HZ cells according to the criteria given in the Results. At the end of the recording an image stack of the cell was acquired in both the red and green channels.
Blind ‘in vivo’ patch clamping
Request a detailed protocolThis method was used to measure the EPSGs used in the in vitro simulation in Figure 7B. The dura was removed but a coverslip was not used. Patch electrodes were prepared and filled as for targeted patch clamping, then an electrode was positioned at the surface of the anterior PCx using a dissection microscope. The micromanipulator (MP-285, Sutter Instrument) was zeroed at the surface, then the electrode was advanced rapidly to the search depth (150–280 µm, corresponding to layer 2) while applying high pressure (25 kPa), after which the pressure was reduced (4–8 kPa) and a whole-cell recording obtained using standard techniques (Margrie et al., 2002; Poo and Isaacson, 2009). Data were acquired as described above. EPSGs were calculated from EPSCs recorded under voltage clamp at a holding potential of –70 mV, close to the chloride reversal potential for these solutions. Layer 2 principal neurons (SL and SP cells) were identified by their intrinsic electrical properties (Suzuki and Bekkers, 2006) and by the recording depth (SL: 150–200 µm; SP: 200–280 µm). Cell identity was also confirmed by fixing the brain at the end of the experiment and recovering the morphology of the recorded neuron as previously described (Suzuki et al., 2014).
‘In vivo’ functional Ca2+ imaging
Request a detailed protocolImaging used the red-shifted Ca2+ indicator Cal-590 AM (AAT Bioquest, Sunnyvale, CA), which was prepared and injected as previously described (Tischbirek et al., 2015). Briefly, the dura was removed and dye (1 mM) was pressure-injected into the PCx at a depth of ~200 µm using a glass pipette (tip diameter ~10 µm). A coverslip was glued over the PCx and imaging commenced >1 hr after injection. Imaging frames were acquired at 30 Hz using a custom-modified B-scope two-photon microscope (Thorlabs, Newton, NJ) with a 16×/0.8 NA water immersion objective (Nikon, Tokyo, Japan), resonance-galvanometer scanners and a Chameleon Ultra Ti:Sapphire laser (Coherent, Santa Clara, CA) tuned to 800–820 nm. Cells were included in the dataset if they unambiguously satisfied the fluorescence and soma morphology criteria given in Results, and if they exhibited clear odor responses according to the criteria under Data analysis, below.
Odor presentation
Request a detailed protocolA custom-built flow-dilution olfactometer was used to deliver up to 15 odors which were diluted to 10% of their saturated vapor pressure in charcoal-filtered medical air (flow rate 1 L/min; Bozza et al., 2004). The odors used in this study were: acetophenone, anisole, benzaldehyde, butyric acid, ethanol, ethyl-n-butyrate, eugenol, 1-heptanal, 2-heptanone, lavender oil, limonene, 1-pentanol, propionic acid (Merck/Sigma-Aldrich, St Louis, MO), amyl acetate and isoamyl acetate (MP Biomedicals, Sydney, Australia). Odors were presented for 3 s at 60 s intervals. Control experiments used a miniPID photo-ionization device (Aurora Scientific, Aurora, Canada) to confirm that odors were presented in a step-like manner (20–80% rise time in 46 ± 2 ms, n = 9 odors; Figure 3—figure supplement 1). However, there was a consistent delay of 242 ± 6 ms (n = 9) between the switching time of the final valve and the arrival of the odorant at the detector inlet (Figure 3—figure supplement 1). All odor arrival times in the paper have been corrected for this delay. The mouse was freely breathing and its respiration was recorded using a piezoelectric strap (Kent Scientific, Torrington, CT) around the abdomen. Control experiments confirmed that the onset of inhalation coincided with the start of the downward step visible in some respiration traces (e.g. Figure 1—figure supplement 1) and the beginning of exhalation corresponded to the peak of the upward spike. Because the upward spike was a more reliable feature, we estimated inhalation onset by reference to the exhalation. In a subset of experiments in which both features were clear, we measured the mean latency from the peak of the exhalation spike to the start of the following inhalation (239 ± 8 ms; mean time between exhalation spikes, 365 ± 12 ms; n = 10 mice). Hence, expressed as a fraction of a respiration cycle, inhalation onset occurred at 0.654 ± 0.007 after the exhalation spike; this was used to estimate the time of onset of inhalation for all experiments. Odor stimulus onset was defined as the time of the first estimated onset of inhalation that occurred following the corrected odor arrival time. To avoid habituation, each odor was presented only once. Cross-habituation between different odors was not observed.
Slice experiments
Request a detailed protocolCoronal slices (300 µm thick) were obtained from the anterior PCx of GAD67-GFP (Δneo) mice aged 20–30 days, as previously described (Suzuki and Bekkers, 2006; Suzuki and Bekkers, 2012). Briefly, slices were prepared on a tissue slicer (Campden Instruments, Loughborough, UK) in ice-cold high-Mg2+ cutting solution comprising (in mM) 125 NaCl, 3 KCl, 0.5 CaCl2, 6 MgCl2, 25 NaHCO3, 1.25 NaH2PO4, 2 ascorbate, 3 pyruvate, and 10 glucose (osmolarity 305 mOs/kg), bubbled with 5% CO2/95% O2 (carbogen). The slices were incubated for 40 min at 34°C in carbogen-bubbled artificial cerebrospinal fluid (ACSF; composition below) then were held at room temperature until required.
Whole-cell patch-clamp recordings were made from visually identified GFP-positive interneurons using an Olympus BX51WI microscope equipped with infrared differential interference contrast and wide-field fluorescence, as described previously (Suzuki and Bekkers, 2012). Slices were superfused with warmed ACSF containing (in mM) 125 NaCl, 3 KCl, 2 CaCl2, 1 MgCl2, 25 NaHCO3, 1.25 NaH2PO4, and 25 glucose (310 mOs/kg), bubbled with carbogen and maintained at 33°C ± 1°C. For the experiments in Figures 6B and 7 the bath solution contained 6,7-dinitroquinoxaline-2,3-dione (DNQX, 10 µM) and D-aminophosphonovaleric acid (D-AP5, 50 µM) to block ionotropic glutamate receptors (Tocris, Abingdon, UK). For the experiments in Figure 5B the bath solution instead contained 100 µM picrotoxin (Merck/Sigma-Aldrich), while in Figure 6C and Figure 6—figure supplement 2 it contained 20 µM SR 95531 (gabazine; HelloBio, Bristol, UK), in both cases to block GABAA receptors. Patch electrodes had resistances of 4–6 MΩ when filled with KMeSO4- or CsMeSO3-based internal solution (same as used in vivo). Unless stated otherwise, compounds were obtained from Merck/Sigma-Aldrich.
Recordings were obtained using the same instrumentation and software as for the in vivo experiments. The stimulating electrode was made from a low-resistance patch electrode (~1 MΩ) filled with ACSF and coated with conductive paint. For connected pair recordings (Figure 6B), the pre- and postsynaptic electrodes were filled with KMeSO4- and CsMeSO3-based internal solutions, respectively, allowing the postsynaptic neuron to be voltage clamped at a holding potential of +50 mV while the presynaptic neuron was allowed to remain at its resting potential in current clamp mode. For the dual recordings (Figure 6C, Figure 6—figure supplement 2), both recording electrodes were filled with the KMeSO4-based internal solution, and trains of extracellular stimuli (5 × 200 µs at 20 Hz) were delivered by a constant-voltage stimulator (Digitimer, Welwyn Garden City, UK) over the range 0–85 V in 0.1–5 V increments. Cells were included in the dataset if they had a stable resting potential more hyperpolarized than –60 mV (current clamp) or a stable holding current <100 pA when held at –70 mV (voltage clamp). At the end of the experiment the slice was fixed and processed with an ABC kit (Vector Laboratories, Burlingame, CA), allowing recovery of the morphologies of the recorded neurons.
Patterned stimulation and dynamic clamp experiments in slices
Request a detailed protocolIn vivo-like patterns of inhibitory synaptic stimulation in slice experiments (Figure 7) were adjusted to a common respiration time base as follows. First, a subset of odor-evoked firing patterns from NG and HZ cells, together with their associated respiration traces, was randomly selected from the full dataset for each cell type (sample data in Figure 7—source code 1). For each subset, one firing/respiration combination was chosen as a reference, and the first upward peak in the reference respiration trace after odor onset was defined as t0. Every other odor-evoked firing pattern in that subset was translated in time to align its first respiration peak after odor onset to t0. Working forward and backward from t0, for each respiration interval the respiration trace (and associated firing pattern) of each other firing pattern was excised from the original data recording and linearly stretched (‘warped’) so the duration of that respiration interval matched the corresponding interval in the reference respiration trace. Finally, all of these warped segments were concatenated in their original order, and the resultant stimulus patterns (examples in Figure 7A, red traces) were used as the trigger to the extracellular stimulator in slice experiments. The same method was used to align the EPSG to the inhibitory stimulation patterns (Figure 7B, Figure 7—source code 1). In this case the respiration trace for the EPSG was used as the reference time base.
The EPSG used in Figure 7B was obtained from an EPSC recorded blind in vivo from an SP cell in response to a 3-s-long application of ethyl-n-butyrate (data in Figure 7—source code 1). Although the in vivo EPSC was measured near the chloride reversal potential to minimize contamination due to direct inhibition onto that cell, the EPSC will still be the result of excitatory input from all other neurons in the circuit, including other SP cells that have received their own feedforward inhibition. Furthermore, injecting an excitatory conductance at the soma does not accurately replicate excitatory synaptic inputs distributed across the dendrites. These details were disregarded for the purpose of this simple in vitro simulation. For each neuron the conductance magnitude was adjusted to produce a similar firing rate with the extracellular stimulus switched off, then the stimulator was switched on to record the effect of synaptic inhibition. To provide a reference firing rate for normalizing the PSTH (Figure 7B), a fixed conductance stimulus was inserted at the end of the EPSG, well past the odor period (not visible in Figure 7B). The dynamic clamp was implemented using Igor Pro (Wavemetrics, Lake Oswego, OR; Figure 7—source code 1).
Data analysis
Request a detailed protocolAll analysis was done using Igor Pro, Axograph X, Matlab (MathWorks, Natick, MA), or R (running under RStudio, Boston, MA).
AP properties (Figure 1D) were measured as previously described for slice experiments (Suzuki and Bekkers, 2010a). Latency to the first AP was the delay from the beginning of the current step to the first AP at rheobase. AP height was the difference between the peak of the first AP at rheobase and its threshold voltage (defined as the membrane potential [Vm] at which dVm/dt first exceeded 15 V/s). AP halfwidth was the width of the first AP at rheobase measured at half its height. Respiration-correlated oscillations in Vm (Figure 2) were characterized as follows. Respiration traces were band-pass filtered at 1–50 Hz and normalized to oscillate between 0 and 1. The time between successive exhalation spikes was used to calculate the mean respiration frequency in three 3-s-long windows (−6 to –3 s, 0–3 s, and 16–19 s, all with respect to odor onset). The same three windows were used for the analysis of oscillations in Vm. First, Vm was median-filtered to remove APs. Mean peak-to-peak amplitude of the oscillations in Vm was found by excising segments of Vm between successive exhalation spikes, linearly stretching (‘warping’) them so they had the same time axis, then averaging together all segments within each of the three analysis windows. Mean peak-to-peak amplitude was taken as the difference between the maximum and minimum amplitude of the averaged segment of Vm. For cross-covariance analysis (Figure 2—figure supplement 1, Figure 2—figure supplement 2), the respiration trace was replaced by a normalized series of Gaussians representing the upward spikes in the respiration cycle; this was done to eliminate the effect of differences in the amplitude of the recorded respiration trace. The amplitude of the largest positive peak in the covariance was found, as well as the location of this peak expressed as a fraction of the mean time between exhalation spikes for that window.
The odor response indices were found by: (i) calculating the fraction of odors that each cell responded to, then averaging across cells (‘cell-averaged index’) and (ii) calculating the fraction of cells each odor activated, then averaging across odors (‘odor-averaged index’). Odor responsiveness was determined using a z-score criterion: median-filtered Vm segments demarcated by exhalation spikes were excised as above, and the mean Vm in each of these segments was calculated for the entire odor trial. These values were converted to a z-score by subtracting the mean and dividing by the standard deviation of the list of mean Vm values during the 8-s-long pre-odor baseline period. An odor response was said to occur if this z-score value exceeded 2.5 for any respiration segment during the 3-s-long odor application period (Figure 1—figure supplement 1). Lifetime and population sparseness were calculated with modified expressions (Pashkovski et al., 2020) that used the peak z-scored response (without thresholding) as the variable.
Individual odor-evoked EPSPs (Figure 4A), EPSCs (Figure 5C), and IPSCs (Figure 6A) were notch-filtered at 2–4 Hz, if required, to remove respiration-associated oscillations prior to averaging (see also Figure 4—figure supplement 1). To estimate the latency to peak and halfwidth, each individual EPSC was fitted to the equation I(t) = a(1 − e-t/τ1)2e-t/τ2, where a is amplitude and τ1 and τ2 are the rising and falling time constants, respectively. The latency and halfwidth were then measured from this fitted curve.
In vivo Ca2+ imaging experiments (Figure 3, Figure 3—figure supplement 2) were analyzed by manually drawing regions of interest around the somas of identified neurons in the green channel then measuring mean somatic fluorescence, F, in the red channel for each frame. Baseline fluorescence (F0) for each neuron was defined as the median of the lower 80% of all values of F measured in all frames acquired prior to odor application. Plots show ΔF/F0 = (F − F0)/F0.
Connected-pair recordings (Figure 6B) and dual recordings (Figure 6C) in slices were analyzed as follows. The peak postsynaptic current was found within a 5- or 10-ms-long window starting 2 ms after the time of the presynaptic AP or extracellular stimulus, then the mean peak amplitude (apk) was found by averaging over a 1 ms long window around the peak and subtracting the averaged baseline over a 5- or 10-ms-long window ending 2 ms before the time of the presynaptic AP or extracellular stimulus. For the connected-pair recordings (Figure 6B) the same measurement procedure was repeated for 200 randomly-chosen times during the 500-ms-long baseline period preceding the time of the presynaptic AP, and the standard deviation (sb) of these baseline mean amplitudes was calculated. A synaptic connection was identified if |apk| > |3 sb|. For the dual recordings (Figure 6C) the mean NG:HZ EPSC amplitude ratio was calculated as follows. For each stimulator setting, all evoked EPSC amplitudes at that setting (typically n = 5) were averaged together and the ratio was calculated. In the case of experiments with a clear plateau (e.g. Figure 6C, bottom left), the individual ratios for EPSCs on the plateau were averaged together to yield the overall mean NG:HZ ratio for that experiment. In the case of experiments without a clear plateau (e.g. Figure 6—figure supplement 2), the overall mean only included individual ratios for which both EPSCs were >10 pA.
Statistical analysis
Request a detailed protocolAll statistical analysis was done using R (version 3.6.0) running under RStudio. Sample sizes were not predetermined using a statistical test; we established that our sample sizes were sufficient from the size and statistical significance of the results, and our sizes are similar to those commonly used in the field. Data collection and analysis were not blinded or randomized, but analysis was automated whenever possible. Results are presented as mean ± standard error of the mean (SEM) with associated exact p value (n = number of cells, cell pairs, or cell-odor pairs, as indicated). Pairwise comparisons were done using Welch’s unpaired two-tailed t-test (t.test function in R). Multiple paired comparisons used one-way ANOVA with Tukey’s contrasts (lmer and glht functions in R) with blocking by cell and odor number. Visual inspection of residual plots did not reveal any large deviations from homoscedasticity. Distributions were compared using the KS test (ks.test). Significance is indicated on the figures by ns (not significant) or by one or more asterisks, following the usual convention in R.
Data availability
Source data files have been provided for Figures 1, 2 and 4. Source code and sample data have been provided for Figure 7.
References
-
Olfactory perceptual stability and discriminationNature Neuroscience 11:1378–1380.https://doi.org/10.1038/nn.2217
-
Two dynamically distinct inhibitory networks in layer 4 of the neocortexJournal of Neurophysiology 90:2987–3000.https://doi.org/10.1152/jn.00283.2003
-
Neurons and circuits for odor processing in the piriform cortexTrends in Neurosciences 36:429–438.https://doi.org/10.1016/j.tins.2013.04.005
-
Synaptic organization of cat olfactory cortex as revealed by intracellular recordingJournal of Neurophysiology 32:204–214.https://doi.org/10.1152/jn.1969.32.2.204
-
Recurrent cortical circuits implement concentration-invariant odor codingScience (New York, N.Y.) 361:1–12.https://doi.org/10.1126/science.aat6904
-
Cortex is driven by weak but synchronously active thalamocortical synapsesScience (New York, N.Y.) 312:1622–1627.https://doi.org/10.1126/science.1124593
-
In vivo whole-cell recording of odor-evoked synaptic transmission in the rat olfactory bulbThe Journal of Neuroscience 23:4108–4116.https://doi.org/10.1523/JNEUROSCI.23-10-04108.2003
-
Normalization as a canonical neural computationNature Reviews. Neuroscience 13:51–62.https://doi.org/10.1038/nrn3136
-
Bidirectional plasticity of cortical pattern recognition and behavioral sensory acuityNature Neuroscience 15:155–161.https://doi.org/10.1038/nn.2966
-
Efficient codes and balanced networksNature Neuroscience 19:375–382.https://doi.org/10.1038/nn.4243
-
Variable coupling between olfactory system activity and respiration in ketamine/xylazine anesthetized ratsJournal of Neurophysiology 93:3573–3581.https://doi.org/10.1152/jn.01320.2004
-
Slow-waves in the olfactory system: an olfactory perspective on cortical rhythmsTrends in Neurosciences 29:429–437.https://doi.org/10.1016/j.tins.2006.06.013
-
Persistent CaMKII activation mediates learning-induced long-lasting enhancement of synaptic inhibitionThe Journal of Neuroscience 35:128–139.https://doi.org/10.1523/JNEUROSCI.2123-14.2015
-
Unitary analysis of opossum prepyriform cortexJournal of Neurophysiology 36:762–774.https://doi.org/10.1152/jn.1973.36.4.762
-
Association and commissural fiber systems of the olfactory cortex of the rat. II. Systems originating in the olfactory peduncleThe Journal of Comparative Neurology 181:781–807.https://doi.org/10.1002/cne.901810407
-
Structure of the piriform cortex of the opossum. I. Description of neuron types with Golgi methodsThe Journal of Comparative Neurology 213:163–187.https://doi.org/10.1002/cne.902130205
-
Optophysiological analysis of associational circuits in the olfactory cortexFrontiers in Neural Circuits 6:1–17.https://doi.org/10.3389/fncir.2012.00018
-
The anesthetic mechanism of urethane: the effects on neurotransmitter-gated ion channelsAnesthesia and Analgesia 94:313–318.https://doi.org/10.1097/00000539-200202000-00015
-
Dendritic compartment and neuronal output mode determine pathway-specific long-term potentiation in the piriform cortexThe Journal of Neuroscience 29:13649–13661.https://doi.org/10.1523/JNEUROSCI.2672-09.2009
-
Olfactory oscillations: the what, how and what forTrends in Neurosciences 32:207–214.https://doi.org/10.1016/j.tins.2008.11.008
-
Circuit oscillations in odor perception and memoryProgress in Brain Research 208:223–251.https://doi.org/10.1016/B978-0-444-63350-7.00009-7
-
Diversity among principal and GABAergic neurons of the anterior olfactory nucleusFrontiers in Cellular Neuroscience 8:111.https://doi.org/10.3389/fncel.2014.00111
-
Long-term synaptic plasticity in hippocampal interneuronsNature Reviews. Neuroscience 8:687–699.https://doi.org/10.1038/nrn2207
-
Inhibition by somatostatin interneurons in olfactory cortexFrontiers in Neural Circuits 10:1–17.https://doi.org/10.3389/fncir.2016.00062
-
GABAergic circuits control input-spike coupling in the piriform cortexThe Journal of Neuroscience 28:8851–8859.https://doi.org/10.1523/JNEUROSCI.2385-08.2008
-
Theta oscillation coupled spike latencies yield computational vigour in a mammalian sensory systemThe Journal of Physiology 546:363–374.https://doi.org/10.1113/jphysiol.2002.031245
-
Membrane and synaptic properties of pyramidal neurons in the anterior olfactory nucleusJournal of Neurophysiology 105:1444–1453.https://doi.org/10.1152/jn.00715.2010
-
Encoding of odor fear memories in the mouse olfactory cortexCurrent Biology 29:367–380.https://doi.org/10.1016/j.cub.2018.12.003
-
BookOlfactory cortexIn: Shepherd GM, editors. The Synaptic Organization of the Brain. New York: Oxford University Press. pp. 100–736.
-
Interactions between behaviorally relevant rhythms and synaptic plasticity alter coding in the piriform cortexThe Journal of Neuroscience 32:6092–6104.https://doi.org/10.1523/JNEUROSCI.6285-11.2012
-
Neurogliaform cells in cortical circuitsNature Reviews. Neuroscience 16:458–468.https://doi.org/10.1038/nrn3969
-
NMDA spikes enhance action potential generation during sensory inputNature Neuroscience 17:383–390.https://doi.org/10.1038/nn.3646
-
Input normalization by global feedforward inhibition expands cortical dynamic rangeNature Neuroscience 12:1577–1585.https://doi.org/10.1038/nn.2441
-
Spike coding in pyramidal cells of the piriform cortex of ratJournal of Neurophysiology 86:1504–1510.https://doi.org/10.1152/jn.2001.86.3.1504
-
Novel behavioral paradigm reveals lower temporal limits on mouse olfactory decisionsThe Journal of Neuroscience 35:11667–11673.https://doi.org/10.1523/JNEUROSCI.4693-14.2015
-
Sparse odor coding in awake behaving miceThe Journal of Neuroscience 26:8857–8865.https://doi.org/10.1523/JNEUROSCI.0884-06.2006
-
Synaptic organization of anterior olfactory nucleus inputs to piriform cortexThe Journal of Neuroscience 40:9414–9425.https://doi.org/10.1523/JNEUROSCI.0965-20.2020
-
Four unique interneuron populations reside in neocortical layer 1The Journal of Neuroscience 39:125–139.https://doi.org/10.1523/JNEUROSCI.1613-18.2018
-
Single dendrite-targeting interneurons generate branch-specific inhibitionFrontiers in Neural Circuits 8:139.https://doi.org/10.3389/fncir.2014.00139
-
Neural coding by two classes of principal cells in the mouse piriform cortexThe Journal of Neuroscience 26:11938–11947.https://doi.org/10.1523/JNEUROSCI.3473-06.2006
-
Inhibitory neurons in the anterior piriform cortex of the mouse: classification using molecular markersThe Journal of Comparative Neurology 518:1670–1687.https://doi.org/10.1002/cne.22295
-
Microcircuits mediating feedforward and feedback synaptic inhibition in the piriform cortexThe Journal of Neuroscience 32:919–931.https://doi.org/10.1523/JNEUROSCI.4112-11.2012
-
Green fluorescent protein expression and colocalization with calretinin, parvalbumin, and somatostatin in the GAD67-GFP knock-in mouseThe Journal of Comparative Neurology 467:60–79.https://doi.org/10.1002/cne.10905
-
Spine neck plasticity regulates compartmentalization of synapsesNature Neuroscience 17:678–685.https://doi.org/10.1038/nn.3682
-
Speed and accuracy of olfactory discrimination in the ratNature Neuroscience 6:1224–1229.https://doi.org/10.1038/nn1142
-
Sensory network dysfunction, behavioral impairments, and their reversibility in an Alzheimer’s beta-amyloidosis mouse modelThe Journal of Neuroscience 31:15962–15971.https://doi.org/10.1523/JNEUROSCI.2085-11.2011
-
BookLearning to Smell: Olfactory Perception from Neurobiology to BehaviorBaltimore, MD: The Johns Hopkins University Press.
-
Single-unit activity in piriform cortex during slow-wave state is shaped by recent odor experienceThe Journal of Neuroscience 30:1760–1765.https://doi.org/10.1523/JNEUROSCI.5636-09.2010
-
Neuronal cell-type classification: challenges, opportunities and the path forwardNature Reviews. Neuroscience 18:530–546.https://doi.org/10.1038/nrn.2017.85
-
Diverse patterns of odor representation by neurons in the anterior piriform cortex of awake miceThe Journal of Neuroscience 30:16662–16672.https://doi.org/10.1523/JNEUROSCI.4400-10.2010
Article and author information
Author details
Funding
National Health and Medical Research Council (1009382)
- John M Bekkers
National Health and Medical Research Council (1050832)
- John M Bekkers
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
This work was performed in strict accordance with the guidelines contained in the Australian Code for the Care and Use of Animals for Scientific Purposes, published by the National Health and Medical Research Council of Australia. All experimental procedures were approved by the Animal Experimentation Ethics Committee of the Australian National University (Protocols A2014/11, A2017/27, A2018/52). All surgeries and experiments were performed under anesthesia with either urethane or fentanyl/medetomidine, with supplementary isoflurane as required, and every effort was made to minimize suffering.
Copyright
© 2022, Suzuki et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,874
- views
-
- 212
- downloads
-
- 4
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Neuroscience
Otolith organs in the inner ear and neuromasts in the fish lateral-line harbor two populations of hair cells oriented to detect stimuli in opposing directions. The underlying mechanism is highly conserved: the transcription factor EMX2 is regionally expressed in just one hair cell population and acts through the receptor GPR156 to reverse cell orientation relative to the other population. In mouse and zebrafish, loss of Emx2 results in sensory organs that harbor only one hair cell orientation and are not innervated properly. In zebrafish, Emx2 also confers hair cells with reduced mechanosensory properties. Here, we leverage mouse and zebrafish models lacking GPR156 to determine how detecting stimuli of opposing directions serves vestibular function, and whether GPR156 has other roles besides orienting hair cells. We find that otolith organs in Gpr156 mouse mutants have normal zonal organization and normal type I-II hair cell distribution and mechano-electrical transduction properties. In contrast, gpr156 zebrafish mutants lack the smaller mechanically evoked signals that characterize Emx2-positive hair cells. Loss of GPR156 does not affect orientation-selectivity of afferents in mouse utricle or zebrafish neuromasts. Consistent with normal otolith organ anatomy and afferent selectivity, Gpr156 mutant mice do not show overt vestibular dysfunction. Instead, performance on two tests that engage otolith organs is significantly altered – swimming and off-vertical-axis rotation. We conclude that GPR156 relays hair cell orientation and transduction information downstream of EMX2, but not selectivity for direction-specific afferents. These results clarify how molecular mechanisms that confer bi-directionality to sensory organs contribute to function, from single hair cell physiology to animal behavior.
-
- Neuroscience
Animals navigate by learning the spatial layout of their environment. We investigated spatial learning of mice in an open maze where food was hidden in one of a hundred holes. Mice leaving from a stable entrance learned to efficiently navigate to the food without the need for landmarks. We developed a quantitative framework to reveal how the mice estimate the food location based on analyses of trajectories and active hole checks. After learning, the computed ‘target estimation vector’ (TEV) closely approximated the mice’s route and its hole check distribution. The TEV required learning both the direction and distance of the start to food vector, and our data suggests that different learning dynamics underlie these estimates. We propose that the TEV can be precisely connected to the properties of hippocampal place cells. Finally, we provide the first demonstration that, after learning the location of two food sites, the mice took a shortcut between the sites, demonstrating that they had generated a cognitive map.