Odor identity coding by distributed ensembles of neurons in the mouse olfactory cortex

  1. Benjamin Roland
  2. Thomas Deneux
  3. Kevin M Franks
  4. Brice Bathellier  Is a corresponding author
  5. Alexander Fleischmann  Is a corresponding author
  1. Center for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS UMR 7241, INSERM U1050, France
  2. Centre National de la Recherche Scientifique, UPR 3293, France
  3. Duke University, United States

Abstract

Olfactory perception and behaviors critically depend on the ability to identify an odor across a wide range of concentrations. Here, we use calcium imaging to determine how odor identity is encoded in olfactory cortex. We find that, despite considerable trial-to-trial variability, odor identity can accurately be decoded from ensembles of co-active neurons that are distributed across piriform cortex without any apparent spatial organization. However, piriform response patterns change substantially over a 100-fold change in odor concentration, apparently degrading the population representation of odor identity. We show that this problem can be resolved by decoding odor identity from a subpopulation of concentration-invariant piriform neurons. These concentration-invariant neurons are overrepresented in piriform cortex but not in olfactory bulb mitral and tufted cells. We therefore propose that distinct perceptual features of odors are encoded in independent subnetworks of neurons in the olfactory cortex.

https://doi.org/10.7554/eLife.26337.001

Introduction

All sensory systems must be able to unambiguously determine stimulus identity in the face of variable stimulus intensity. In vision, for example, the perception of colors is stable across a wide range of luminance, despite the strong dependency of photoreceptor activation on wavelength and light intensity (Nunn et al., 1984). This problem is of particular importance in olfaction, given the massive and rapid fluctuations in odorant concentration, for example encountered in odor plumes in the environment (Murlis and Jones, 1981).

Molecular features of odorants are detected by odorant receptors. Odorant receptors are expressed on the dendrites of sensory neurons in the olfactory epithelium, and odorant receptors are broadly tuned such that each odorant typically activates multiple receptors (Jiang et al., 2015; Malnic et al., 1999). A given olfactory sensory neuron expresses one of a large repertoire of odorant receptors (Buck and Axel, 1991; Zhang and Firestein, 2002), and neurons expressing a given receptor project to two spatially invariant glomeruli in the olfactory bulb (Mombaerts, 2001). Thus, the molecular features of an odorant are represented as a discrete and stereotyped map of glomerular activity (Bozza et al., 2004; Ma et al., 2012; Soucy et al., 2009; Uchida et al., 2000). Odor information encoded in patterns of glomerular activity must then be integrated at higher olfactory centers in the brain to generate unified odor objects, defined by perceptual features such as odor identity and intensity (Gottfried, 2010; Wilson and Sullivan, 2011; Wojcik and Sirotin, 2014). The piriform cortex has been suggested to serve as such a site of integration.

The piriform cortex is a simple, three-layered paleocortical structure, which receives dense projections from mitral and tufted cells, the main output neurons of the olfactory bulb. Mitral and tufted cells extend an apical dendrite into a single glomerulus, and thus only receive direct excitatory input from sensory neurons expressing a single odorant receptor. Odor information encoded in the spatio-temporal patterns of mitral and tufted cell activity is then transmitted to higher olfactory centers, including the piriform cortex. Mitral and tufted cells project axons to large areas of the piriform cortex, without identifiable topographic organization (Ghosh et al., 2011; Igarashi et al., 2012; Nagayama et al., 2010; Sosulski et al., 2011). Individual piriform neurons receive inputs from multiple and broadly distributed glomeruli, thus providing an opportunity for molecular feature integration (Apicella et al., 2010; Miyamichi et al., 2011). Consistent with this model, optical stimulation of the olfactory bulb suggests that piriform neurons respond to combinations of co-active glomeruli (Davison and Ehlers, 2011; Haddad et al., 2013). Calcium imaging and electrophysiological recordings show that odors activate sparse ensembles of piriform neurons, which are distributed across the piriform cortex without apparent spatial organization (Poo and Isaacson, 2009; Rennaker et al., 2007; Stettler and Axel, 2009; Tantirigama et al., 2017). However, how information about the identity and intensity of an odor is encoded in the response patterns of piriform ensemble activity remains poorly understood.

The ability to unambiguously identify odors across a wide range of concentrations is essential for olfactory perception and behavior (Cleland et al., 2011; Stopfer et al., 2003; Wojcik and Sirotin, 2014). Indeed, rats can identify odorants with consistently high accuracy over a greater than 50,000-fold range in concentration (Homma et al., 2009). However, the specificity of odorant - receptor binding steeply depends on odorant concentration (Jiang et al., 2015; Malnic et al., 1999), and consequently, patterns of odor-evoked glomerular activity change with changing odorant concentration (Bozza et al., 2004; Meister and Bonhoeffer, 2001; Rubin and Katz, 1999; Spors and Grinvald, 2002). Similarly, the spatio-temporal patterns of mitral and tufted cell strongly depend on odorant concentration (Banerjee et al., 2015; Bathellier et al., 2008; Fukunaga et al., 2012; Kato et al., 2013; Miyamichi et al., 2013; Sirotin et al., 2015). Therefore, we sought to determine if the piriform cortex is capable of forming concentration-invariant representations of odor identity from sensory input that is concentration-dependent.

Here, we have used in vivo two-photon calcium imaging in anesthetized mice to record odor responses from large, unbiased ensembles of piriform neurons. We find that odor identity can accurately be decoded from the spatial patterns of local piriform odor responses. However, we also observe that piriform odor representations change substantially across a 100-fold range in odor concentration, degrading information about the identity of the odor. We propose a solution for this potential confound by identifying a subpopulation of concentration-invariant piriform neurons, which accurately encodes odor identity across a broad range of odor concentrations. These concentration-invariant neurons are present at numbers significantly above chance in piriform ensembles but not in olfactory bulb mitral and tufted cells, indicating that the ability to form concentration-independent representations of odor identity in functionally distinct neural subnetworks is an emergent property of piriform cortex.

Results

Representations of odor identity in the piriform cortex

We stereotaxically injected adeno-associated virus (AAV) expressing the calcium indicator GCaMP6s (Chen et al., 2013) into the piriform cortex of adult (6- to 10-week-old) mice. Ten days after infection, we surgically exposed the piriform cortex for two-photon imaging under ketamine/xylazine anesthesia. Viral expression of GCaMP6s resulted in dense labeling of layer II piriform neurons (Figure 1a and b, and Materials and methods). To monitor calcium signals in such large and densely packed neural ensembles, we developed an automated cell segmentation algorithm based on calcium signal similarities (Figure 1—figure supplement 1 and Materials and methods). This algorithm operates on the entire data set obtained from individual imaging sites and efficiently identifies individual neurons by iteratively clustering neighboring pixels with high-signal covariance. Under these conditions, we could simultaneously monitor the activity of 100–400 (mean ± SD: 242 ± 105) neurons per imaging site.

Figure 1 with 1 supplement see all
Calcium imaging of odor-evoked activity in piriform cortex.

(a) (Top) Schematic of the experimental protocol. AAV-GCaMP6s was stereotaxically injected into piriform cortex. After 10 days, piriform cortex was surgically exposed and neural activity in response to odors was recorded with two-photon imaging. (Bottom) Coronal section of a brain used for piriform imaging. Scale bar: 100 μm. (b) (Top) Stack average of an imaging site, and (bottom) masks of the regions of interest (ROIs) identified by our clustering algorithm (see also Figure 1—figure supplement 1). Scale bar: 50 μm. (c) Single-trial example traces of four different cells (rows) in response to four different monomolecular odorants (columns). Cells respond with an increase (red traces) or a decrease (blue traces) in fluorescence (gray traces: non responsive trials). Note the different scale of ∆F/F values (y axis) for each neuron. Red bar: odor presentation. (d) Spatial patterns of piriform activity in response to four trials (columns) of three different monomolecular odorants (rows). ∆F/F values are clipped at 100% here and henceforth in all figures. Odorants activate sparse, distributed, and partially overlapping ensembles of piriform neurons. (e) Reliability of activation (red) and suppression (blue), measured as the number of trials each cell-odor pair responded to a given odorant, averaged across imaging sites (n = 6 sites in three mice). Error bars: 95% CI of the mean. (f) Percent of neurons activated (red) or suppressed (blue), averaged across four trials. Dots: single data points from individual imaging sites. Horizontal bars: mean across imaging sites (n = 6). Error bars: 95% CI of the mean. (g) Tuning curve of activation (red) and suppression (blue), averaged across imaging sites (n = 6). Error bars: 95% CI of the mean.

https://doi.org/10.7554/eLife.26337.002

In initial experiments, to establish population coding properties of large ensembles of piriform neurons, we measured responses to a test panel of 13 different monomolecular odorants. Two-second odor pulses of these stimuli (1:10,000 dilution in mineral oil, 0.01% vol./vol.) elicited relatively sparse but partially overlapping activity of piriform neurons, consistent with one previously published report (Stettler and Axel, 2009) (Figure 1c and d). We observed that individual neurons responded selectively with an increase or decrease in fluorescence. Across six imaging sites in three mice (total number of neurons = 1706), 20% (±8.4%) of the neurons responded with an increase in fluorescence, and 11% (±5.7%) of the neurons responded with a decrease in fluorescence (Figure 1f). Most neurons responding with an increase in fluorescence exhibited narrow stimulus tuning, with the exception of a small subpopulation of neurons (8.5 ± 5.4%) that responded to 12 or all 13 odorants of the test panel. Neurons responding with a decrease in fluorescence exhibited similarly selective odor tuning, but only a minimal number of broadly tuned neurons could be observed (Figure 1g). Strikingly, many neurons displayed high trial-to-trial variability in response to the repeated delivery of the same odorant (Figure 1e).

We next sought to quantify the similarity of population responses evoked by different odorants. For each trial, we constructed population activity vectors, defined as the mean temporally deconvolved change in fluorescence of all simultaneously recorded cells of an imaging site over a 4 s time window after stimulus onset (Figure 2a). We then computed pairwise cross-correlations between all single-trial population activity vectors. The results of this analysis for the imaging site shown in Figure 1b are displayed in Figure 2b as a cross-correlation matrix. As indicated by the square boxes along the diagonal (4 × 4 trials), repeated exposure to the same odorant triggered highly correlated population response patterns (intra-odorant cross-trial correlation coefficients, neurons pooled across six imaging sites: 0.67 ± 0.07 (across odorants)). However, this cross-correlation analysis also revealed that patterns elicited by different odorants were fairly similar (mean inter-odorant correlation coefficient: 0.44 ± 0.08). Such overlap may, at least in part, be a consequence of correlated noise at a given imaging site. We therefore pooled neurons across imaging sites and projected this pseudo-population onto the first three principle components in principal component space (Figure 2c). While response patterns for a few odorants appeared to segregate from each other, substantial overlap remained between neural ensembles encoding different odorants.

Figure 2 with 1 supplement see all
Odor identity is encoded in the spatial patterns of piriform activity.

(a) Population response of the imaging site shown in Figure 1a. The mean ∆F/F value after stimulus onset for each trial (columns) in each cell (rows) is color-coded. (b) Similarity matrix obtained by computing the pairwise correlation coefficients between all population response vectors in (a). Squares along the diagonal (4 × 4 trials) represent the similarity of responses to four exposures of the same odorant (intra-stimulus cross-trial correlations). (c) Patterns of piriform activity in response to single odor presentations (dots represent different odorants, color-coded as in (a)) projected onto space of the first three principal components (accounting for ~20% of the total variance). Neurons were pooled across imaging sites. (d) Confusion matrix summarizing the performance of a linear classifier trained to discriminate the odorants in A, summed over imaging sites (n = 6 sites in three mice). (e) Percent of trials correctly identified by a linear classifier, for neurons pooled across all six imaging sites (blackline). Green, red and blue lines represent classification accuracy obtained with a classifier trained on response patterns at 0.27 (green arrowhead), 1 (red arrowhead) and 3 (blue arrowhead) s after odor onset. Gray dashed line: theoretical chance level. Gray square: odor exposure.

https://doi.org/10.7554/eLife.26337.004

This raises the question as to whether the spatial patterns of odor-evoked piriform activity we can observe with calcium imaging, which lacks precise temporal information at small time scales, contain sufficient information to discriminate between different odorants. To address this question, we tested the accuracy with which a linear classifier correctly identify odorants based on single-trial response patterns (see Materials and methods). Despite both trial-to-trial variability and considerable overlap in response patterns to different odorants, the classifier could correctly identify all 13 stimuli in our odorant test panel (Figure 2d). Classification accuracy reached 94% when pooling all six imaging sites (Figure 2e). For individual imaging sites, the average classification accuracy was 71% (±5.8%) (Figure 2—figure supplement 1). Correct classification slowly rises after odor onset to reach high accuracy within 0.5 s. This relatively slow rise can largely be explained by the slow rise time of the calcium indicator GCAMP6s (Chen et al., 2013). Interestingly, classification accuracy remained significantly above chance level for several seconds after odor delivery, despite a marked decline in odor-evoked fluorescence (Figure 2e). This finding is consistent with the observation that odor-specific network configurations persist for extended periods of time in the olfactory bulb (Bathellier et al., 2008; Friedrich and Laurent, 2001). To determine whether this sustained piriform activity is odor-specific, we trained the classifier at defined time points after odor onset (0.27, 1, and 3 s, Figure 2e, colored arrows) and measured classification accuracy over time. We observed that a classifier trained on response patterns 1 s after odor onset accurately predicted odor identity during odor delivery, but classification accuracy readily declined after odor offset. In contrast, a classifier based on response patterns at 0.27 and 3 s after odor onset yielded a lower classification accuracy; however, classification success was more stable over time. This suggests that early and late (post-odor offset) representations are based on weak activity configurations present during odor exposure but masked by another much stronger component, which quickly vanishes after the odor is withdrawn. Thus, odor representations are dynamically rearranged after odor offset to maintain odor information for several seconds.

Previous anatomical and functional experiments suggested that odor information was encoded in randomly distributed ensembles of piriform neurons (Ghosh et al., 2011; Illig and Haberly, 2003; Sosulski et al., 2011; Stettler and Axel, 2009). Therefore, to test whether odorant-selective piriform neurons were clustered or randomly distributed across individual imaging sites, we selected cells that were significantly modulated by odorants (1-way ANOVA, p<0.05) and mapped their odorant preference (i.e. the odorant that triggered the strongest response, see example map in Figure 3a). As a sensitive measure of spatial clustering, we computed for each imaging site a nearest neighbor index (NNI Theodoridis and Koutroumbas, 2009), as the mean distance of a cell preferring a given odor to the nearest cell preferring the same odor. We then computed the NNI distribution for 1000 spatially shuffled maps to estimate the probability p that the observed value can be explained by random spatial organization (null hypothesis). Using this measure, we found no evidence for spatial clustering in any of the six imaging sites (p values of 0.49, 0.81, 0.09, 0.21, 0.91, 0.79). To validate the sensitivity of this test to detect spatial clusters, we again simulated populations but now with subtle inhomogeneities in the spatial distributions of neurons. We divided the simulated site into 16 equally sized sub-areas and allowed preference for a given odorant to occur in only three or four randomly chosen sub-areas (Figure 3b). In this case, comparing observed and simulated NNI resulted in highly significant values for p.

No evidence for topography in patterns of piriform activity.

(a) Odorant preference map for an imaging site, with the odorant preference of individual piriform neurons color coded as in Figure 2. (b) Example of simulated clustering. The imaging site was divided into 16 equally sized sub-areas and preference for an odorant was allowed to occur in only three or four randomly chosen sub-areas. The p value indicates the probability that the computed nearest neighbor index is different from randomly distributed neurons. (c) (Top) Example of an ensemble of piriform neurons (in red) locally constrained around a ‘starter cell’ (in yellow), and used for classification of the stimulus set of Figure 2. (Bottom) Heatmaps of the classification accuracy of different starter cells for piriform ensembles of increasing size.

https://doi.org/10.7554/eLife.26337.006

These observations suggest that information about odor identity is homogeneously distributed across individual imaging sites. To further illustrate this idea, we sequentially defined each cell within an imaging site as a ‘starter cell’, and iteratively built neural ensembles of increasing size by adding neighboring neurons (Figure 3c). We observed that different clusters of the same size encoded odorant identity with similar accuracy, and that no ‘hotspots’ of classification were observed. Finally, we found that classification accuracy was very similar for the six different imaging sites in three mice, which were up to 1 mm apart along the rostro-caudal axis of the piriform cortex (Figure 2—figure supplement 1, and data not shown).

Taken together, our data show that odorants activate spatially distributed ensembles of piriform neurons with significant overlap and variability, and no apparent spatial organization. Furthermore, the spatial patterns of odor-evoked piriform activity contain sufficient information to robustly decode stimulus identity.

Patterns of piriform activity decorrelate with increasing concentrations

The ability to correctly classify individual trials as corresponding to a given odorant reveals that odor response patterns are different, but does not directly address how odor identity is encoded within piriform ensembles. Importantly, perceived odor identity of monomolecular odorants typically remains stable across a large range of odor intensities (Homma et al., 2009), providing an experimental opportunity to more directly address this question. Therefore, to test if odor identity coding is stable across changing odorant concentrations, that is concentration-invariant, we next analyzed piriform neural activity in response to three different odorants over a 100-fold range in concentrations (acetophenone, ethyl acetate, and hexanone; at 1:10,000, 1:1,000, 1:100 vol./vol. dilution in mineral oil). We confirmed, using photoionization detector (PID) measurements, that odorant presentations were reliable and scaled according to volumetric ratios (Figure 4—figure supplement 1). We analyzed 13 imaging sites in 11 mice (total number of neurons = 2935). We found that the fraction of both activated and suppressed neurons increased moderately but statistically significantly with increasing concentrations (Figure 4a and b, mixed effect ANOVA, significant effect of concentration on the fraction of activated neurons: F(2,96) = 52.49, p<0.01, and on the fraction of suppressed neurons: F(2,96) = 28.79, p<0.01). Thus, population sparseness decreases with increasing concentrations. We then asked how changes in odorant concentration impacted the response properties of individual neurons. To follow the evolution of each cell’s odorant selectivity, we calculated lifetime sparseness across the 100-fold range in concentration (see Materials and methods). We observed diverse concentration-induced changes in cell selectivity across the population, indicating that an increase in odorant concentration did not systematically broaden odor tuning (Figure 4c and d). Indeed, when considered across the entire population, the distribution of lifetime sparseness of individual neurons was maintained across stimulus intensities (mean ± SD: 0.01%: 0.35 ± 0.23, 0.1%: 0.37 ± 0.24, 1%: 0.38 ± 0.24, Figure 4e).

Figure 4 with 1 supplement see all
Odor-evoked activity and sparseness of individual piriform neurons depends on odorant concentration.

(a) Spatial patterns of piriform activity in response to acetophenone, ethyl acetate and hexanone at three different concentrations (1:10,000, 1:1,000, 1:100 vol./vol. dilutions). (b) Percent of neurons activated or suppressed by acetophenone, ethyl acetate, and hexanone at three different concentrations (1:10,000, 1:1,000, 1:100 vol./vol. dilutions). Dashed gray lines represent individual imaging sites, thick red (activation) and blue (suppression) lines represent averages across sites (n = 13 sites in 11 mice). Error bars: 95% CI of the mean. (c) Spatial patterns of piriform activity in response to acetophenone, ethyl acetate and hexanone at three different concentrations (1:10,000, 1:1,000, 1:100 vol./vol. dilutions). Only cells responding at least two out of four trials are depicted. Cells responding to multiple odorants are color-coded in white. (d) Matrix of lifetime sparseness across concentrations for the cells in (c), sorted by hierarchical clustering. (e) Distribution of lifetime sparseness of all neurons pooled across imaging sites (n = 13, total number of neurons = 2935).

https://doi.org/10.7554/eLife.26337.007

To visualize changes in the population response patterns with increasing odorant concentration, we next rearranged the population response matrix of the imaging site shown in Figure 3a, using hierarchical clustering (Figure 5a and b). This analysis further supported the observation that neural responses varied with the odorant stimulus, and that response magnitudes could increase or decrease with increasing odorant concentrations. To quantify the similarities of response patterns elicited by odorants at different concentrations, and to evaluate the contributions of these different response profiles to the encoding of the identity of an odorant, we next calculated the cross-trial correlations between individual response vectors of neurons pooled across 13 imaging sites (Figure 5c). We first noted that pair-wise trial correlations increased with concentration, indicating that trial-to-trial variability across the population decreases at higher concentrations (mean ± SD across odorants: 0.01%, 0.67 ± 0.03; 0.1%, 0.72 ± 0.06; 1%, 0.74 ± 0.05). Note that this decrease in trial-to-trial variability partly accounts for the increased number of cells with responses that are statistically significantly different from mineral oil at higher concentrations (Figure 4b). Furthermore, we found that responses to a given odorant at a given concentration were significantly more correlated than responses to a given odorant at different concentrations (mixed-effect ANOVA, effect of concentration F2,96 = 51.39, p<0.01). This concentration-dependent decorrelation was gradual: the average correlation between response patterns elicited by the three odorants at low (0.01% vol./vol.) and intermediate (0.1% vol./vol.) concentrations was 0.53 (±0.12), while the average correlation between response patterns at low and high (1% vol./vol.) odorant concentrations dropped to 0.36 (±0.12) (Figure 5c and d). Thus, the spatial patterns of piriform odor responses were not concentration-invariant. Instead, response patterns to a given odorant changed substantially as odor concentration increased and became as different over a 100-fold change in concentration as responses to a different odorant.

Patterns of piriform activity decorrelate with increasing odorant concentrations.

(a) Population response of the imaging site in Figure 3a to acetophenone, ethyl acetate and hexanone at three different concentrations (1:10,000, 1:1,000, 1:100 vol./vol. dilutions). Cells are sorted by hierarchical clustering. (b) Example response profiles of cells suppressed at higher concentrations for all three odorants (top panel), cells moderately enhanced by increasing concentrations of ethyl acetate or hexanone (middle panel), and cells strongly enhanced by increasing concentrations of acetophenone (bottom). (c) Similarity matrix obtained by computing the pairwise correlation coefficients between all response vectors pooled across imaging sites. Squares along the diagonal (4 × 4 trials) represent the similarities of responses to a single odorant/concentration pair (intra-stimulus cross-trial correlations). Large squares (12 × 12 trials) represent the similarities of responses to an odorant at varying concentrations (intra-odorant cross-trial correlations). The similarities of responses to increasing concentrations of a given odorant (intra-odorant inter-concentrations similarity) are highlighted by dashed line rectangles. (d) Correlation coefficients of the patterns of piriform activity elicited at increasing concentrations with the patterns elicited at low concentrations. Dashed gray lines represent individual imaging sites, thick black lines the average across sites (n = 13 sites in 11 mice). Error bars: 95% CI of the mean. Patterns of piriform activity along increasing concentrations gradually decorrelate from the patterns elicited at low concentration.

https://doi.org/10.7554/eLife.26337.009

A concentration-invariant subpopulation of piriform neurons

Because overall piriform patterns of piriform activity change across odorant concentrations, we next asked if we could identify subnetworks of concentration-invariant piriform neurons that can represent the identity of an odor independent of its intensity. We used a linear regression approach (Rigotti et al., 2013) to determine if individual neurons were present within piriform ensembles whose response profiles could be accounted for solely by the identity of the odorants, irrespective of their concentrations. We performed an analysis of variance of each cell’s response profile, with odorant concentration and odorant identity as the two explanatory variables (Figure 6a). Consistent with our qualitative inspection of response profiles (Figure 5a and b), we found that 37.2 ± 11.7% of cells were selectively responsive to odorant, while 27.8 ± 10.6% of cells showed responses that were modulated by an interaction between odorant identity and concentration (Figure 6a). Crucially, of the cells that were selectively modulated by odorant, 30.1% ± 11.6%, constituting 10.4 ± 4.5% of the total population of cells, exhibited concentration-invariant responses for all three odors, according to the analysis of variance (Figure 6c). Representative response traces of concentration-invariant neurons are shown in Figure 6b. Thus, although many neurons displayed mixed selectivity to the odorant identities and concentrations tested in our experiment, a substantial fraction of concentration-invariant neurons can be identified within piriform neural ensembles.

Figure 6 with 2 supplements see all
A concentration-invariant subnetwork of piriform neurons.

(a) Summary of the analysis of variance of the response profiles of cells in Figure 4a. Two ANOVAs are performed successively. In the first step, neurons with mixed selectivity are excluded from the analysis. Concentration-invariant cells are then identified as cells significantly modulated by odorant identity but not concentration (red arrowheads: example cells). These cells are used for subsequent analyses against the population of all other ‘generic neurons’. Red dashed line: significance threshold, set at p=0.01. (b) Deconvolved response traces of 4 concentration-invariant cells (rows) to acetophenone, ethyl acetate and hexanone at three different concentrations (light grey: 1:10,000, dark grey: 1:1,000, black 1:100 vol./vol. dilutions). ∆F/F values are normalized to each cell’s maximum ∆F/F. Red bar: odor presentation. Shaded area: time interval used to integrate ∆F/F values for the analysis of variance in (a). (c) Response matrix of the population of cells in (a) exclusively modulated by odorant identity. ∆F/F values do not vary significantly with increasing odorant concentrations. (d) Similarity matrix obtained by computing the correlation coefficients between patterns of concentration-invariant neurons pooled across imaging sites (n = 13). Responses of concentration-invariant neurons to different concentrations of a given odorant (intra-odorant inter-concentrations similarity) are as correlated as responses to the same odorant/concentration pair. (e) Patterns of activity of generic neurons (left) or concentration-invariant neurons (right) in response to single-odor presentations (dots) projected onto space of the first three principal components. Neurons were pooled across imaging sites (n = 13). (f) Confusion matrix summarizing the accuracy of the classification of odorant identity by generic neurons (left) or concentration-invariant neurons (right) in a generalization task (see Materials and methods), summed over imaging sites (n = 13). The classifier assigns each odor trial to one of seven stimulus groups. Two concentrations of the tested odorant (including the tested concentration) are excluded from the training data (white boxes). Difficult generalization tasks across a 100-fold change in odorant concentration are highlighted by the red squares.

https://doi.org/10.7554/eLife.26337.010

To verify our selection procedure for concentration-invariant piriform neurons, we quantified the similarities of their response patterns across the 100-fold odorant concentration range. We found that, for a given odorant, responses remained highly correlated across increasing concentrations (mean inter-concentration similarity ± SD across odorants, 0.01% to 0.1%: 0.60 ± 0.09; 0.01% to 1%: 0.57 ± 0.10, Figure 6d, and Figure 6—figure supplement 1). These across-concentration response correlations for concentration-invariant subpopulations approached those for within-concentration responses observed in the general population (0.65 ± 0.10), although a small but significant difference remained (Wilcoxon rank sum test, p<0.05).

To qualitatively display differences in the odor representation formed by the concentration-invariant subnetwork of neurons compared to the rest of the piriform ensemble, we next projected their population response patterns onto the first three principal components in principal component space. This analysis revealed that odorant representations of concentration-invariant neurons clustered irrespective of concentration, while odorant representations of generic neurons clustered by concentrations of the same odorants, but not systematically by odorant groups (Figure 6e and Figure 6—figure supplement 2). To more quantitatively evaluate the odor-coding properties of concentration-invariant piriform subnetworks and to test if concentration-invariant neurons could generalize odorant identity across changing odorant concentrations, we next trained a linear classifier to predict odorant identity based on a single concentration. We then tested the classifier on all other odorants and concentrations (‘generalization learning’, Figure 6f, see Materials and methods). We found that the classification accuracy of odorant identity across a 10-fold change in odorant concentration was similar between the concentration-invariant subpopulation of neurons and the entire population of ‘generic’ neurons. However, for the more difficult generalization tasks across a 100-fold change in odorant concentration (0.01% to 1%, and 1% to 0.01%, Figure 5f, red squared boxes), subnetworks of concentration-invariant neurons were much more accurate in predicting odorant identity than the entire population of ‘generic’ neurons (1% to 0.01% generalization: concentration-invariant cells mean = 63.5 ± 27.1%, Mann-Whitney test n = 13 U=20, Benjamini and Hochberg’s FDR adjusted p<0.01, ‘generic’ neurons mean = 19.2%±17.8%, Mann-Whitney test n = 13 U=82, FDR adjusted p=0.39). Taken together, these analyses demonstrate that for a given set of odorants and concentrations it is possible to identify a subpopulation of piriform neurons that can encode the identity of an odor largely independent of odor intensity.

How stable are concentration-invariant piriform subnetworks with varying stimulus intensity range and complexity? To address this question, we first identified concentration-invariant neurons across a 10-fold instead of a 100-fold range in odorant concentration. Using the same selection criteria as for the original data set, this analysis yielded 19.2% (±7.4% SD across experiments) concentration-invariant neurons. Second, we identified concentration-invariant neurons for pairs of two odorants instead of the three odorants in our test panel. We found that 22% (±8%) of the cells were concentration-invariant for at least one pair of odorants. Of those cells, 49% were identified as concentration-invariant for all three odorants while the response of the other 51% of cells was modulated by the concentration of one of the three odorants. Taken together, these data suggest that concentration-invariant subnetworks of neurons can be modulated by stimulus complexity and concentration, yet remain relatively stable within the range of our stimulus set.

We then tested if concentration-invariant neurons exhibited response profiles that differentiate them from other neurons. We compared response rise time (from 10% to 90% of peak ΔF/F, concentration-invariant neurons: 0.55 ± 0.05 s; generic neurons: 0.52 ± 0.10 s) and duration (width at 50% peak ΔF/F, concentration-invariant neurons: 2.02 ± 0.07 s; generic neurons: 2.01 ± 0.11 s), odor-evoked peak change in fluorescence (deconvolved ΔF/F, concentration-invariant neurons: 10.0 ± 0.7%; generic neurons: 9.2 ± 1.3%), and trial-to-trial variability, but found no significance between the two populations (p>0.05, Wilcoxon ranked sum test, n = 13) (Figure 6—figure supplement 1). We next examined whether concentration-invariant cells were spatially clustered. Visual inspection of the localization of concentration-invariant cells at individual imaging sites did not reveal obvious clustering (Figure 6—figure supplement 1). Furthermore, performing the statistical analysis based on the nearest neighbor index (see Figure 3), we found that the organization of concentration-invariant neurons at 11 out of 13 imaging sites was undistinguishable from the random distribution obtained for shuffled data. In 2 of the 13 imaging sites, spatial distributions were moderately but significantly different from random.

Taken together, our analysis identifies a subpopulation of piriform neurons, with response profiles and spatial distributions that are similar to other odor-responsive neurons, but which encode odor identity independent of concentration.

Concentration-invariant neurons are overrepresented in piriform cortex but not in the olfactory bulb

A representation of odor identity emerges in piriform cortex from the integration of odor-evoked mitral and tufted cell activity from the olfactory bulb. Considerable normalization of odor-evoked neural activity across a range of odorant concentrations has been observed in the olfactory bulb (Banerjee et al., 2015; Kato et al., 2013; Miyamichi et al., 2013; Roland et al., 2016; Zhu et al., 2013), suggesting that concentration-invariant piriform odor responses could be inherited from the olfactory bulb. Alternatively, the formation of segregated, concentration-invariant odor identity representations in subpopulations of piriform neurons may emerge within cortex itself. Therefore, to distinguish between these models, we next analyzed odor-evoked responses of olfactory bulb mitral and tufted cells. We used previously described mitral/tufted cell calcium imaging data (Roland et al., 2016) (Figure 7a), obtained under equivalent experimental conditions but using GCaMP3 instead of GCaMP6 as the calcium indicator. We then performed an analysis of variance of each cell’s response profile, as described above for the piriform imaging data (19 imaging sites in eight mice, total number of cells = 523, Figure 7b).

Figure 7 with 1 supplement see all
Concentration-invariant representations of odor identity emerge in the piriform cortex.

(a) (left) Schematic of the experimental protocol. Rabies-GCaMP3 was injected underneath the lateral olfactory tract. After 5–7 days, the olfactory bulb was surgically exposed and mitral cell activity in response to odorants was recorded with two-photon imaging. (middle) Stack average of an imaging site. Scale bar: 50 μm. (right) Deconvolved response traces of two neurons to acetophenone, ethyl acetate and hexanone at three different concentrations (light grey: 1:10,000, dark grey: 1:1,000, black 1:100 vol./vol. dilutions). Red bar: odor presentation. Shaded area: time interval used to integrate ∆F/F values for the analysis of variance in (b). (b) Summary of the analysis of variance of the response profiles of a mitral and tufted cell imaging site. See also Figure 5a, and Materials and methods for details. (c) Population response of the imaging site in (a) to acetophenone, ethyl acetate and hexanone at three different concentrations (1:10,000, 1:1,000, 1:100 vol./vol. dilutions). Cells are sorted by the p-value of the effect of odorant identity. (d) Percent of concentration-invariant neurons (red line) identified in the olfactory bulb (left, n = 19 imaging sites in eight mice) and in piriform cortex (right, n = 13 imaging sites in 11 mice), overlaid onto the distribution of the percent of concentration-invariant neurons found in the bootstrap samples. (e) p-Values for the number of concentration-invariant neurons identified at each imaging site in the olfactory bulb (left, n = 19) and the piriform cortex (right, n = 13). The number of concentration-invariant neurons is significantly above chance (red dots, p<0.01) in 8 out of 13 imaging sites in the perform cortex, but only in 2 out of 19 imaging sites in the olfactory bulb.

https://doi.org/10.7554/eLife.26337.013

Visual inspection of the response matrix suggested that the majority of mitral cell responses are concentration-dependent (Figure 7c). Consistent with this impression, our analysis of variance indicated that 25.9% ± 14.9% of neurons exhibited odorant-selective responses. Of these, only 20.5 ± 15.1%, constituting only 5.2 ± 5.9% of the total population were concentration-invariant. Thus, many fewer neurons were exclusively modulated by odorant identity in the olfactory bulb compared to the piriform cortex (Figure 7d) (olfactory bulb, 5.2 ± 5.9%, piriform cortex 10.4 ± 4.5%, Mann-Whitney test, nob = 19, npir = 13, U = 59, p<0.01). However, we imaged more piriform neurons (2935) than olfactory bulb mitral cells (523). Therefore, to ensure that differences in the fractions of concentration-invariant neurons did not result from biased sampling, we subsampled piriform cortex to match the numbers of olfactory bulb cells. Furthermore, we tested whether relaxing the significance criterion from p<0.01 to p<0.05 would change our results. We find that the fraction of concentration-invariant cells in piriform cortex is consistently and significantly higher than that observed in the olfactory bulb, independent of sampling size and the significance criterion used in our model (Figure 7—figure supplement 1, Materials and methods). Note also that, as a consequence of the lower dynamic range of GCaMP3 compared to GCaMP6s (Chen et al., 2013), we are likely to overestimate the concentration-invariance of mitral cell responses, so that the difference in the fraction of concentration-invariant cells in piriform cortex compared to olfactory bulb may be greater than this analysis suggests.

Responses in both piriform and mitral cells are highly heterogeneous. We therefore used a bootstrapping analysis to ensure that the population of concentration-invariant neurons, we observed are indeed over-represented and are not merely a consequence of the inherent response variability across population within each dataset. We repeated the analysis of variance on 10,000 shuffled cell-odor pairs, in which cell identities were scrambled across stimuli but the population statistics for each odorant response were preserved (see Materials and methods). Shuffling mitral cell identities across stimuli indeed resulted in the identification of the same percentage of concentration-invariant neurons observed experimentally (bootstrap mean 3.9 ± 0.8%; observed value 5.5%, p=0.073). By contrast, shuffling cell identities from the piriform dataset yielded significantly fewer concentration-invariant neurons than observed experimentally (bootstrap mean 5.3 ± 0.41%; observed value 11.7%, p<0.001, Figure 7e). Finally, we computed the significance of these findings for each individual imaging site. We considered the observed number of concentration-invariant neurons per imaging site to be significantly above chance if it was higher than 99% of the values calculated from the bootstrap samples. We found that the observed number of concentration-invariant neurons in piriform was above chance in 8 out of 13 imaging sites. In contrast, only 2 out of 19 mitral and tufted cell imaging sites contained concentration-invariant neurons above chance level (Figure 7f). Together, these results indicate that the encoding of odor identity independent of concentration is robust in piriform cortex, but not in olfactory bulb mitral and tufted cells.

Discussion

We have examined how the identity of an odor is represented in neural ensemble activity in the piriform cortex using two-photon calcium imaging in anesthetized mice. We found that despite substantial overlap between response patterns evoked by different odorants, odor identity could correctly be predicted from local ensembles of piriform neurons. Furthermore, we observed that piriform response patterns across the population change substantially with increasing odorant concentration, potentially confounding odor identification. However, a substantial fraction of odor-selective piriform neurons exhibit largely concentration-invariant odor responses, and odor identity - independent of intensity - could accurately be decoded from this subpopulation of piriform neurons. Concentration-invariant neurons are overrepresented in piriform cortex but not in olfactory bulb mitral and tufted cells, suggesting that concentration-invariant subnetworks for odor identity emerge in cortical neural circuits for olfaction.

The structure of odor identity-encoding piriform ensembles

Previous immunohistochemical, electrophysiological and imaging experiments have revealed that piriform odor representations are distributed across a large area of the cortex, and that piriform neurons display discontinuous receptive fields (Poo and Isaacson, 2009; Rennaker et al., 2007; Stettler and Axel, 2009; Yoshida and Mori, 2007). These observations led to the speculation that information about odor is encoded by ensembles of coordinately active neurons distributed across piriform cortex without topographic organization (Stettler and Axel, 2009). Here, we explicitly tested this prediction, and we found that odor identity could indeed be decoded from spatially distributed patterns of odor-evoked piriform activity. Our results are consistent with recent data from extracellular recordings in awake rats (Miura et al., 2012) and mice (see co-submitted manuscript by Bolding and Franks), which show that odor identity can be accurately decoded from the firing rates of piriform neural ensembles. However, extracellular recordings cannot reveal the spatial organization of the neurons that participate in odor coding. Our data show that odorant-selective neurons do not cluster in space. Furthermore, odor identity can be decoded with similar accuracy from multiple different imaging sites, and information about odor identity appears to be homogeneously distributed within an individual imaging site. Thus, our results provide robust evidence that odor information is encoded without topographic organization by ensembles of piriform neurons.

Comparison with extracellular recordings in awake mice

Bolding and Franks (2017) have used extracellular recordings to explore how odor identity and intensity is encoded in piriform cortex. It is interesting to compare their results with our results obtained from optical imaging. The majority of piriform neurons are narrowly tuned to odor, exhibit high trial-to-trial variability to repeated presentations of the same odor, and substantial overlap between response patterns elicited by different odors. Piriform odor representations are less sparse than previously suggested, and odor identity can accurately be decoded from piriform ensembles in the absence of precise temporal information, consistent with an earlier report (Miura et al., 2012). Both studies report that piriform odor representations are not concentration-invariant: odorants presented at a range of different intensities (30-fold concentration range in Bolding and Franks, 100-fold in this study) elicit highly dissimilar response patterns, as dissimilar as the response patterns observed for two different odorants. An important difference between the two studies is that, due to the surgical preparation required to expose piriform cortex for two-photon imaging, optical recordings were performed in ketamine/xylazine-anesthetized mice while electrophysiological recordings were obtained from awake, head-fixed mice. Thus, despite the potentially diverse effects of anesthesia on odor sampling and neural physiology, key features of piriform odor responses identified in the two studies are very similar.

Bolding and Franks utilize the high temporal resolution of electrophysiological recordings to propose that information about odor intensity can be encoded by the synchrony of piriform odor responses. Our study, on the other hand, describes the spatial organization of piriform odor representations and provides evidence for a concentration-invariant subnetwork of piriform neurons that encodes odor identity - independent of odor intensity. Thus, the two studies propose complementary coding schemes for non-interfering representations of odor identity and odor intensity in the mouse olfactory cortex. Interestingly, Bolding and Franks also observe that a subpopulation of neurons exhibits short-latency concentration-invariant odor responses. While differences in the experimental design of the two studies preclude a direct comparison, it is tempting to speculate that these represent the same neural subpopulations.

Anesthesia interferes with active sniffing, and awake mice and rats typically exhibit faster sniff rates that can be dynamically modulated by odor (Blauvelt et al., 2013; Carey and Wachowiak, 2011; Wachowiak et al., 2013). In contrast, we do not observe significant modulation of sniff rate in anesthetized mice. Anesthesia has also been shown to modulate the activity of neurons in the olfactory bulb. In awake mice, periglomerular and granule cell inhibitory neurons exhibit higher levels of activity (Cazakoff et al., 2014; Kato et al., 2012; Wachowiak et al., 2013), while some mitral cells exhibit diminished spontaneous and odor-evoked activity (Kato et al., 2012; Kollo et al., 2014; Shusterman et al., 2011). Anesthesia is likely to differentially affect different types of piriform neurons and has recently been shown to modulate baseline neural activity in piriform cortex (Tantirigama et al., 2017). Of note, at elevated odorant concentrations, recordings in awake mice show that the fraction of activated neurons remains stable, while imaging experiments in anesthetized mice show a moderate increase (this study, and Stettler and Axel, 2009). This observation suggests that normalization of neural activity is incomplete under anesthesia, consistent with recent reports on the activity of cortical feedback projections to the olfactory bulb (Boyd et al., 2015; Otazu et al., 2015; Rothermel and Wachowiak, 2014). The activity of cortical feedback projections has been suggested to contribute to signal normalization and is attenuated under anesthesia.

Odor identity and intensity

The ability to accurately determine stimulus identity independent of intensity is critical for olfactory perception and behavior (Cleland et al., 2011; Sirotin et al., 2015; Stopfer et al., 2003). Behavioral experiments have shown that rats can be trained to identify, with high accuracy, monomolecular odorants across a greater than 50,000-fold range in concentration (Homma et al., 2009). On the other hand, odor concentration-invariance is not absolute. Behavioral experiments in humans and insects have shown that the perceived identity of an odor can change with concentration (Bhagavan and Smith, 1997; Gross-Isseroff and Lancet, 1988; Laing et al., 2003; Pelz et al., 1997). While the perceptual boundaries of odor concentration-invariance remain poorly defined, it is clear that odor identity and intensity must be, at least in part, independently represented in the brain. Earlier experiments in insects have suggested that odor identity and intensity information is intermingled in the antennal lobe. Interestingly, however, multidimensional manifolds representing odor identity emerged after non-linear dimensionality reduction of the data, suggesting that downstream structures could extract concentration-invariant information about odor identity from antennal lobe activity (Stopfer et al., 2003). Moreover, experiments in the fish olfactory bulb suggest that temporal multiplexing could be used to independently transmit odor identity and intensity information to higher olfactory centers in the brain (Friedrich et al., 2004) (see also companion manuscript by Bolding and Franks).

We propose an alternative, simple solution for representing these two distinct features of an odor stimulus in the mammalian olfactory cortex; that information about odor identity and odor intensity can be separately represented in distinct subpopulations of piriform neurons. We observed that piriform odor representations change systematically with increasing odorant concentrations, such that responses evoked by an odorant at different concentrations can become as dissimilar as responses evoked by two different odorants. However, we found that ~30% of odor-selective neurons support concentration-invariant odor representations. Such subnetworks can provide a stable representation of odor identity, while information about odor concentration can simultaneously be encoded in other neural ensembles.

What cellular and circuit mechanisms could underlie the generation of concentration-invariant piriform odor responses? The most parsimonious model for such functionally distinct subpopulations of piriform neurons is that these neural subpopulations represent distinct piriform neural cell types. Recent work in acute slice preparations has indeed highlighted the functional diversity of piriform layer II neurons. Piriform layer II cells can be classified into semilunar and superficial pyramidal cells, and cells of intermediate phenotype (Suzuki and Bekkers, 2011; Wiegand et al., 2011). Superficial semilunar cells have higher input resistance and shorter membrane time constants than pyramidal cells in deep piriform layer II, and semilunar cells receive stronger excitatory input from the olfactory bulb, but weaker associational input than pyramidal cells. Such differences in intrinsic properties and synaptic connectivity could underlie some of the heterogeneity in response types we observe. Interestingly, layer II semilunar and pyramidal cells project to distinct piriform target areas (Chen et al., 2014; Diodato et al., 2016), providing an opportunity to selectively transmit distinct features of the odor stimulus to different targets. For example, layer II pyramidal cells send cortical feedback projections to the olfactory bulb, which have been implicated in signal normalization - a function that primarily relies on information about stimulus intensity but may be largely independent of odor identity (Boyd et al., 2015; Otazu et al., 2015). On the other hand, semilunar cells projecting to the cortical amygdala, which has been implicated in the encoding of odor valence (Root et al., 2014), may transmit information about odor identity, independent of intensity. The identification of feature-selective subnetworks in piriform cortex, and advances in the characterization of piriform neural connectivity will open up new possibilities for examining odor information routing in cortical neural circuits for olfaction.

Materials and methods

Surgical preparation for piriform imaging

Request a detailed protocol

Adult (6- to 10-week-old) male mice on a mixed genetic background (C57BL/6; 129Sv) were used for experiments. All experiments were performed according to European and French institutional animal care guidelines (protocol number B750512/00615.02). A total of 400–700 nl of AAV-GCaMP6s (AV-1-PV2824, Penn Vector, Philadephia, PA) were stereotaxically injected at multiple sites into the piriform cortex at 0.5–1 mm posterior to bregma, using manually controlled pressure injection. 10–13 days later, the piriform cortex was surgically exposed, following experimental procedures described in Stettler and Axel (2009). Briefly, mice were anesthetized with ketamine/xylazine (100 mg/kg/10 mg/kg, Sigma Aldrich) and a head-post was glued to the skull. Skin was retracted to expose the masseter muscle, the superficial temporal vein was cauterized, and the zygomatic bone was removed with fine scissors. The upper portion of the lower jawbone was cleared from tendons and cut out. Minor bleeding was stopped with gelatin sponge (Gelfoam). A well was constructed around the surgical site with silicone sealant (WPI, Sarasola, FL). A small craniotomy, typically 1 × 2 mm in size, was then drilled over the piriform, and the thinned bone was removed with fine forceps. After removal of the dura, the silicone well surrounding the craniotomy was filled with artificial cerebral spinal fluid (ACSF; 125 mM NaCl, 5 mM KCl, 10 mM glucose, 10 mM HEPES, 2 mM CaCl2, 2 mM MgSO4) at all times. A small glass coverslip was placed over the craniotomy and sealed in place using 2% agarose. Six imaging sites in three mice (total number of neurons = 1706) were analyzed for the 13 odorant test panel, and at 13 imaging sites in 11 mice (total number of neurons = 2935) for acetophenone, ethyl acetate and hexanone at three different concentrations.

Surgical preparation for mitral cell imaging

Request a detailed protocol

Methods for mitral and tufted cell imaging are described in Roland et al., 2016. Briefly, 3 to 3.5 nL of rabies-GCaMP3 virus was slowly pressure injected underneath the LOT. 5–7 days later, mice were anesthetized using ketamine/xylazine and the skull overlying the olfactory bulb was thinned using a dental drill and removed with forceps, and the dura was peeled back using fine forceps. A small circular glass coverslip was placed over the exposed bulb and sealed in place using 2% agarose. Activity at 19 imaging sites in eight mice (total number of cells = 523) was analyzed.

Functional imaging

Request a detailed protocol

A typical piriform imaging experiment lasted between 2 to 3 hr, and a maximum of three different fields of view were imaged per mouse, at a position between 0 to 1.5 mm posterior to Bregma. Body temperature was maintained at 37°C using a feedback-controlled heating pad (FST). Piriform imaging data were acquired on two different microscopes: a Leica SP5 with a 25x Olympus objective, 256 × 256 pixels for a 347 × 347 μm field of view, and a Scientifica Multiphoton VivoScope with a 20x Olympus objective, 256 × 160 pixels for a 357 × 220 μm field of view. A Mai Tai DeepSee laser (Spectra-Physics, Santa Clara, CA) was tuned to 910 nm. Densely packed piriform layer II neurons at a depth of ~250 µM below the pial surface were imaged. 30 s movie sequences were acquired at a frame rate of ~15 Hz. Mitral cell imaging data were acquired on two different microscopes: Ultima, Prairie Technologies with a 16x objective at 2x zoom or Leica SP5 with a 25x Olympus objective. 25 s movie sequences at 256 × 256 pixels were acquired at a frame rate of 2.53 Hz (Ultima) or 2.9 Hz (Leica SP5).

An odor trial lasted 30 s (8 s of pre-stimulus baseline, 2 s of stimulation, 20 s of post-stimulus acquisition). Odors were delivered at a flow rate of 1 L/min with inter-trial intervals of ~35 s. Odor stimuli for a given experiment consisted of one of two odor sets, delivered through a 16-channel olfactometer (Automate Scientific, Berkley, CA): 13 monomolecular odorants (purchased from Sigma Aldrich at the highest purity available), diluted at 1:10,000 vol./vol. in mineral oil (Sigma Aldrich), and a ‘concentration series’ consisting of acetophenone, ethyl acetate, and hexanone, at 10-fold increasing concentrations (1:100, 1:1000 and 1:10,000 vol./vol. dilutions). Odorants were presented four times each, in pseudo-randomized order to avoid habituation (average inter-stimulus interval for the same odor stimulus = 7 min). A photoionization detector (miniPID 200B, Aurora Scientific, Canada) was used to confirm reliable odor delivery and to verify that odorant concentration scales according to volumetric ratios. No stimulus was presented twice in a row. For the concentration series experiment, odorant identity was changed at every trial, that is different concentrations of the same odorant were never presented in a row to avoid adaptation.

Data analysis

Signal extraction

Request a detailed protocol

Data analysis was conducted in Matlab. Motion artifacts were first corrected by using a subpixel translational-based discrete Fourier analysis. Regions of interest (ROIs) were then manually drawn for mitral cell data. For piriform data, ROIs were selected using a semi-automated hierarchical clustering algorithm based on pixel covariance over time (see detailed method below), and the weighted pixel gray value average inside each ROI was used to estimate the fluorescence of single cells at each time frame. The raw fluorescence trace was then upsampled to match the highest sampling rate of each set of experiments (i.e. ~15 Hz for piriform datasets, ~3 Hz for mitral cell datasets). When needed, we corrected for piriform neuropil contamination using a published method (Kerlin et al., 2010): the neuropil signal Fneuropil (t) surrounding each cell was estimated by averaging the signal of all pixels within a 20 μm circular region from the cell center (excluding all other ROIs). The true fluorescence signal of a cell body was estimated as follows: Ftrue (t) = Fmeasured (t) – r x Fneuropil (t), with r = 0.5.

For each trial, the change in fluorescence (∆F/F0) was calculated as (F-F0)/F0, where F0 is the median value between seconds 4 and 8 of the pre-stimulus period. We estimated the baseline fluctuation for a given trial as the standard deviation (SD) of ∆F/F0 during the baseline period. Odor responses were assessed over a 4-s period following odor onset. A ROC analysis (including blank trials consisting of 30 s of recording with no stimulation for the evaluation of the false-positive rate) was used to determine a threshold with the best sensitivity/specificity ratio. A cell was deemed responsive if it reached and remained above threshold (two times the standard deviation of the baseline) for 21 (activation) or 19 (suppression) frames during this response window. These criteria yielded a detection accuracy (ACC) of 0.92 and 0.915, and a true positive rate (TPR) of 0.89 and 0.87, respectively.

To estimate the time course of firing rate, the calcium signal was temporally deconvolved using the following formula: r(t) = f’(t) + f(t) / τ in which f’(t) is the first derivative of f(t) and τ the decay constant set to 2 s for GCaMP6s (as estimated from the decay of the GCAMP6s fluorescent transients), and 0.5 s for GCaMP3. This signal was low-pass filtered using a four-pole Butterworth filter with a cutoff frequency of 2.5 Hz and used for all subsequent analysis, except for example traces in Figure 1, for which the raw signal was plotted.

Automated cell segmentation

Request a detailed protocol

We developed an original method for automatically detecting neurons in the recording region, based on activity time courses. As in other existing methods (Mukamel et al., 2009; Pnevmatikakis et al., 2016), our method uses temporal activity patterns to automatically segment neurons, including the cells undetectable with contrast-based methods because of low fluorescent baseline. But while other methods focus only on the activity time-courses, our method also takes connexity into account, i.e. that pixels belonging to the same neuron are neighboring each other. We perform iterative clustering of all pixels in the image by merging together neighboring pixels or regions whose activities are correlated.

The clustering procedure is schematized in Figure 1—figure supplement 1. At each step, correlations are calculated for each pair of connected pixels in the entire image (with connections only allowed in the horizontal and vertical directions), and the two pixels with the highest correlation are merged together (Figure 1—figure supplement 1a). The averaged signal of the merged region is then computed and connections and correlations with neighbors are updated. This process continues until all pixels have been merged. We observed that the first pixels merged by this procedure belong to neurons (or axons and dendrites), while pixels belonging to diffuse neuropil regions aggregate at later stages, without altering the segmentation of the neurons (Figure 1—figure supplement 1c). Hence, there is a large range of iteration counts for which we obtain a stable number of clusters corresponding to well segmented neurons and a variable number of regions corresponding to neuropil regions (depending on the stage of the clustering). To obtain a first neuronal segmentation, we therefore selected the clustering output obtained at a particular iteration step. This parameter is user selectable, but the final result is weakly sensitive to it. A value of I = 25,000 was typically chosen for these datasets.

We then extract neurons from this first segmentation step by identifying merged pixel ensembles that have typical neuronal size and shapes (Figure 1—figure supplement 1d). To do so, we compute three measures from each aggregated region and retain only regions falling within a particular range:

  1. Size of the region defined as regsize= number of pixels must be included between a minsize and maxsize parameters.

  2. Region dispersion defined as <| | xcenter | |2>xregionregsize must be smaller than a threshold maxdisp (typically 0.5): This will select round-shaped regions and discard elongated regions.

  3. The ratio between average pixel weights (see below for weights calculation) of the pixels located at the border of the region and of all the pixels should be smaller than a threshold maxborder (typically 0.9): This will select regions whose pixel contributions decrease near the border of the region, which is typical of neuron regions but not of neuropil fragment regions.

A graphical interface permits to adjust these selection parameters if necessary. The final segmentation is obtained after a visual quality check in which the user has the possibility to add or remove agglomerated regions within a dedicated graphical interface that helps accelerating the procedure (Figure 1—figure supplement 1d). Importantly, the correlations used in the algorithm are calculated on the full duration of the dataset, but to reduce computational costs, the data is temporally binned into bins of 30 s. Also, two additional preprocessing procedures are applied to remove correlations with large spatio-temporal scales:

  1. Global slow drifts are removed by high-pass filtering all time courses with a cutoff period of 100 s.

  2. The average signal over all pixels is calculated, and its contribution in each individual pixel (i.e. the projection of the pixel time courses onto this average signal) is subtracted.

To compute the ‘average’ signals of the regions obtained at the end of the first neuronal segmentation step, we assign weights to each individual pixel in order to obtain the best estimate of the signal common to all these pixels. Note that the regular average corresponds to the case in which all these weights are equal to 1/N (N being the number of pixels inside the region), but does not necessarily constitute the best estimate of the common signal, as pixel with low signal level and comparatively high noise contribute as much as pixels with high signal levels. To find the appropriate weights, we proceed as follows:

  1. All weights wi are initialized to 1/N

  2. The signal of the region is calculated as:

    xreg=iwixi=Xw

    (where xi are column vectors of individual pixel signals)

  3. Weights are updated as:

W=X+xreg

(where X+=(XTX)1XT is the pseudo-inverse of X). Steps 2 and 3 are repeated (typically three times) until convergence.

Note that the weights are then visualized on our graphical interface, highlighting the regions of the neurons that contribute with the strongest signals.

Population vector analysis

Request a detailed protocol

To build response vectors, we averaged the deconvolved ∆F/F0 signal of all cells over the 4 s following odor onset. This provides an estimate of a neuron’s response to each odor trial. We obtained a matrix (piriform neurons x odors trials) representing the population response after odor delivery for every trial. To build the cross-correlation matrix of the patterns of activity, we calculated Pearson’s correlation coefficients between every pair of such odor trials.

Lifetime sparseness

Request a detailed protocol

Lifetime sparseness was calculated as:

Sl= 1 [ j=1NrjN ]2j=1Nrj2N1 1N

where rj are the neurons’ responses to individual odors and N is the total number of odors. Lifetime sparseness quantifies the specificity of the neurons’ odor-evoked responses (0: uniformly distributed across odors, 1: highly selective for one odor).

Linear classifier

Request a detailed protocol

To quantify the information contained within patterns of piriform activity, we used a linear classifier to predict stimulus identity based upon the response patterns to single odor trials. We obtained comparable classification performance using one-vs-all Support Vector Machine (SVM) with a linear kernel on the raw ∆F/F data, or linear discriminant analysis (LDA) on the principal components encompassing 95% of variance in the data. For computational efficiency, LDA was used for the analysis in Figure S2. To build the response vectors, we accumulated ∆F/F0 signal over a five frame (~333 ms) sliding window (Figure 2), or used a mean response by accumulating signal over the 4 s following odor onset (Figure 5). These vectors define a multidimensional (1 neuron = 1 dimension) representation of odors and were used to classify single-trial response patterns. To avoid overfitting, we used a leave-one-out cross-validation strategy, whereby the assessed trial is excluded from the calculation of the centroids.

To assess the decoding of odors based on stimulus identity (Figure 2), trials were classified as belonging to one of the stimulus group (13 odors, chance level: 1/13 = 7%). To assess the decoding of odor identity across concentration (Figure 5), we adapted a previously published protocol (Stopfer et al., 2003). For each odor trial tested, only the most distant concentration was kept to train the classifier (when testing trials at 0.1% concentration, only trials at 0.01% concentration were kept). Thus, outcome of this classification was between seven possible groups of stimuli: one group with the tested odorant identity, and 2 odorants x 3 concentrations groups of the other odors. Classification was deemed correct if the trial was assigned to the correct identity group (1 group out of 7 possible, 1/7 = 14% chance level).

Analysis of variance of single neurons and bootstrap methods

Request a detailed protocol

To investigate the source of modulation of individual neural responses, we used a linear regression approach by fitting an analysis of variance (ANOVA) with concentration, odorant identity, and trial number as fixed effects, using type II sum of square. We defined concentration-invariant cells as cells that were only significantly modulated by odorant identity (Test 1, p<0.01); containing information that enables them to identify at least one of the three odorants, but were not significantly modulated by concentration (Test 2, p>0.01) or by interactions between identity and concentration (Test 3, p>0.01). ‘Generic neurons’ are all other cells. It is important to note that in our dataset, only 39.9% of all recorded piriform neurons were significantly modulated by odorant identity. Thus, the 11.7% of piriform cells that were identified as concentration-invariant represent ~30% of the neurons that contained significant odorant information. To evaluate the accuracy our method, we estimated an upper bound on the false-positive rate (FPR) of the statistical test. The expected FPR for the intersection of the three tests is the product of the FPRs of each test. Modulation by odor identity (Test 1) is assessed an FPR of 1%, while the FPR for the absence of modulation by intensity (Test 2) and by intensity-identity interactions (Test 3) are not precisely known as they correspond to the false negative rates of the associated tests. However, given that the FPRs for tests 2 and 3 are bounded by 1, the FPR for the concentration-invariant cells is less than 1%. This indicates that at least 10.7% are true positive for concentration invariance (subtracting the 1% FPR). Similarly, of the 5.5% of concentration-invariant neurons identified in the olfactory bulb, at least 4.5% are true positives. To test for the stability of the concentration-invariant subpopulation with respect to the statistical threshold, we repeated the analysis with an alpha-value of 0.05 instead of 0.01. We identified 12.8% of concentration-invariant neurons and, importantly, all concentration-invariant cells previously identified with the 0.01 threshold were included in this ensemble. Thus, the detection of concentration invariant cells is only marginally affected by the exact statistical parameters.

In a second type of analysis, we tested the hypothesis (different from the one tested above) that concentration invariance was a property arising from a random coding scheme, in which response magnitude is arbitrarily assigned for each neuron and odor-concentration pair. We generated 1000 bootstrapped datasets in which we randomly shuffled the cell identities for each odor-concentration pair. In other words, if the data set is described as a 3D array A[Cell_list,Odor_Conc_list,Trials_list], one surrogate dataset AS is generated by performing a randomization of ‘Cell_list’ for each of the nine items of the Odor_Conc_list. Note that individual trials of a given cell-odor-concentration triplet are kept together. We then submitted the surrogate datasets to the statistical analysis described above for detecting concentration-invariant cells. This enabled us to compute the expected distribution of the fraction of concentration-invariant cells for the ‘random coding scheme’ hypothesis. We found that the fraction of cells observed in cortex is incompatible with this hypothesis (Figure 6d, p<0.001), in contrast to the olfactory bulb (Figure 6d, p=0.073).

Statistics

All descriptive statistics in text and figure legends are mean ± SD. The percent of responding neurons to each stimulus was calculated as the average number of active neurons across four trials. To construct the spatial odor maps, as well as to calculate any parameter implying a thresholding of cell activity (e.g. cell tuning), only cells that responded at least twice out of four trials were included.

References

    1. Meister M
    2. Bonhoeffer T
    (2001)
    Tuning and topography in an odor map on the rat olfactory bulb
    Journal of Neuroscience 21:1351–1360.
    1. Pelz C
    2. Gerber B
    3. Menzel R
    (1997)
    Odorant intensity as a determinant for olfactory conditioning in honeybees: roles in discrimination, overshadowing and memory consolidation
    The Journal of Experimental Biology 200:837–847.
  1. Book
    1. Theodoridis S
    2. Koutroumbas K
    (2009)
    Pattern Recognition
    London: Academic Press.

Decision letter

  1. Upinder Singh Bhalla
    Reviewing Editor; National Centre for Biological Sciences, Tata Institute of Fundamental Research, India

In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.

[Editors’ note: a previous version of this study was rejected after peer review, but the authors submitted for reconsideration. The first decision letter after peer review is shown below.]

Thank you for submitting your work entitled "Odor concentration-invariant subnetworks in the mouse olfactory cortex" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and a Senior Editor. Our decision has been reached after consultation between the reviewers. Based on these discussions and the individual reviews below, we regret to inform you that your work will not be considered for publication in eLife at this time. The reviewers felt that the concerns with the paper were too substantial for acceptance, and would not be feasible to address within a short period. They were unsure if the conclusions would still stand after a re-analysis, but did feel that with substantial changes the authors might be able to address the major issues in a resubmission.

A summary of the key points made by the reviewers is below:

This study looks at large-scale coding using optical recordings and concludesthat a subset of piriform neurons respond in a concentration-invariant way,and can be used for identity coding. The dataset is a potentially useful, broadpopulation readout of piriform cortex activity.

1) The reviewers were concerned that the identification, analysis, andrelevance of the concentration-invariant neurons was not well supported.

The authors should be able to devise a model that would more completelyutilize the data.

2) The authors should provide a strong analysis to support theirassertion that the subsample of cells reliably defined odors in aconcentration-invariant way.

3) There should be a more complete analysis of topography of odor responses.

4) This paper doesn't really draw upon the data from the companion paper,and indeed in some ways doesn't seem to match.

There should be a much better discussion of the effects ofanesthesia, which could compare the two datasets.

The analysis of concentration coding in the two datasets could bedone in a comparative manner, with equivalent analyses in the two cases.

5) The authors should provide a better analysis of the dynamics of responses.

Reviewer #1:

This study looks at large-scale coding using optical recordings and concludesthat a subset of piriform neurons respond in a concentration-invariant way,and can be used for identity coding. These are anesthetized recordings, andit would be difficult to do these in awake animals given the quite invasivesurgical preparation.

I review this in context of the companion manuscript by Bolding and Franks.

What do we learn from the current paper that the Bolding manuscript does not say?

The main point here is to identify a subset of concentration-invariant neurons. I'm not sure that this isn't there in the Bolding dataset, but the authors don't seem to have looked for it. I'm surprised that the classification accuracy is so low (69%), with a few hundred neurons on each trial. I'm also surprised at how slowly the classification accuracy builds up and declines. Possibly this is an artifact of the anesthesia.

Overall, I get the sense that the power of optical recordings hasn't beenfully utilized in this study, and that the analysis could be extended toextract more information from the already obtained data. On its own, and inthe present form, the main result of a few cells being concentration-invariant seems limited.

1) The key point of the paper is that piriform odor representations change withincreasing concentrations, but a 10% subset of neurons stay consistent andappear to code only for identity. The authors compare this with an analysisof previous data on OB M/T cells, which are mostly concentration-dependent.

While this result is interesting to olfactory physiologists, I was looking fora more detailed analysis, or a mechanistic insight. For example, ifone confines oneself to fewer odors, does the same subset of neurons stayconsistent? Can the authors devise a model that would span the range fromfully concentration-invariant to concentration-dependent for all odors?

2) I don't see any attempt to relate the optical findings to the electricalrecordings in the companion paper. If one were to filter out the different time-resolutions, would the optical findings simply fall out of the electrical ones? Does the fraction of concentration-independent neurons match?

3) The frame rate for these recordings is quite good, 15 FPS. I am surprisedthere is no analysis of the dynamics of responses for individual neurons,or to compare them between the concentration dependent and independent cells.

4) Did the authors record respiration? Given that the animals were anesthetized, it seems likely that at this frame rate it should be possible to examine respiration dependence of the piriform responses.

5) There are a couple of findings that the authors don't really follow up.

For example, the authors briefly mention functional heterogeneity of PFneuron responses, but they don't seem to follow this up. Also, there isthe finding that the number of responsive Piriform neurons is roughlyindependent of odor concentration.

Reviewer #2:

In their manuscript, 'Odor concentration-invariant subnetworks in the mouse olfactory cortex', Roland et al., investigate how ensemble activity patterns in the olfactory cortex decode odor identity vs. intensity. The authors employ multiphoton imaging in anaesthetized mice and probe GCaMP6s responses to different odors and concentrations. They find that odor identity can be accurately decoded from the population activity using a linear classifier. Across different concentrations, the population representations of odorants change, thus degrading information about odor identity. As a solution to this issue, the authors propose that a subset of concentration-invariant cells in the olfactory cortex is well suited to decode odor identity in a concentration invariant manner.

The study is topical and the experiments are carefully performed. However, I have several major concerns with the interpretation of the results presented here, which, in my opinion, preclude the publication of the manuscript in the current form in eLife.

1) The authors report ~10% concentration-invariant neurons in the PC vs. ~5% in the OB. They perform a shuffling control and suggest that, unlike the OB, the emergence of concentration-invariant cells in the cortex cannot be attributed to chance sampling. Therefore, these neurons represent an emerging feature of the cortex and can solve the identity problem irrespective of changes in concentration.

The fraction of invariant cells is rather small in both the case of cortex and the bulb. The large majority of neurons in the cortex vary substantially in their responses with changes in concentration (Figure 3). Therefore, to this reviewer, asserting that concentration invariance is solved in particular by these neurons is somewhat more anecdotal rather than based on testable evidence. The authors sample 6 times more neurons in the olfactory cortex compared to the bulb (~3,000 vs. 500). Therefore, it is important to determine whether sub-sampling the olfactory cortex data to match the OB recorded cells, and relaxing the significance criterion from p<0.01 to p<0.05 (Figure 6E) change the results presented.

The cross-concentration correlations of odor representations given by even the neurons identified as 'concentration-invariant' are relatively low (0.56, 0.60). What is the average correlation between the neuronal representations of the same odor across repeats, for each of three concentrations sampled? For example, in Figure 5D, for ethyl acetate and hexanone, the correlations across repeats for the highest concentration appear by visual inspection higher than 0.8. This value is substantially higher than the correlations reported between the low-medium, medium-high or low-high concentrations in the sampled range. These observations put into question the concentration-invariant status of these neurons.

Can the authors provide a testable model on how the concentration invariant cells are preserved as such across a wider concentration range (for example 50,000 fold vs. 100 fold sampled), and also how this information is readout selectively?

In addition, in the companion paper (Bolding and Franks), in awake animals, a different mechanism (at the whole population level) is proposed for decoding odor identity, while concentration is read by changes in response latency. Can the authors comment how these two perspectives relate to one another?

2) Comparison of olfactory cortex and olfactory bulb activity patterns is hard to interpret in anaesthetized mice since recent reports suggest that cortical feedback is strongly affected by anaesthesia (Rothermel et al., 2014, Boyd et al., 2015, Otazu et al. 2015).

3) Figure 3's title is not supported by the data. Both the percentage of active neurons and their identity change across the sampled concentration range (3B, C). At the level of individual cells, there is substantial variability in lifetime sparseness across concentrations (3D), even though the percentage of neurons that show similar lifetime sparseness values is same across concentrations. In my opinion, such observations invalidate the title's claim that the overall levels of piriform cortex activity remain largely stable across changes in concentration.

Reviewer #3:

Concentration invariance is known to be a key property of odor identity coding, and the work of Roland et al. is the first (together with a co-submitted article by Bolding and Franks) to study how piriform cortex solves concentration invariance. Furthermore, the authors use a preparation which allows for simultaneous monitoring of large ensembles of spatially-defined neuronal populations, allowing them to examine how concentration-invariance is implemented at the population level.

The authors demonstrate the presence of concentration-invariant neurons, which comprise a subset of piriform neurons. Additionally, a separate sub-population of piriform neurons change their activity as a function of odor concentration. This is an important observation and is worth publishing. However I have a few comments:

1) One of the novel observations is that the concentration-invariant subnetwork lacks topographic organization in piriform cortex. The authors should elaborate more on this point. The case against topographic organization is made in Figure 2—figure supplement 2 by showing similar classification using a local versus random clustering rule. However, this may not rule out that topography exists. (1) It is possible that any single specific topographic mapping may be washed out in the averaging procedure. The authors should perform an additional analysis: artificially simulating topography and performing the exact same analysis to check that this simulated topography actually yields a classification plot that is different from random. (2) Another possible suggestion to examine topography in piriform representations is to examine the spatial distribution of weights in the classifier.

2) The authors made a point about the trial-to-trial variability of piriform neurons (e.g. Figure 1E). A pertinent question is, what is the trial-to-trial variability for neurons in the concentration invariant subnetwork?

3) The author should discuss the effect of anesthesia in their preparations. They compared their research with electrophysiological study co-submitted by Bolding and Franks. The current study revealed the spatial distributions of identity encoding neurons. The authors should discuss the limitations of their approach and clarify what other conclusions can be safely extrapolated for awake case?

[Editors’ note: what now follows is the decision letter after the authors submitted for further consideration.]

Thank you for submitting your article "Odor concentration-invariant subnetworks in the mouse olfactory cortex" for consideration by eLife. Your article has been favorably evaluated by a Senior Editor and three reviewers, one of whom is a member of our Board of Reviewing Editors. The reviewers have opted to remain anonymous.

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

Summary:

The reviewers felt that the paper was stronger in this resubmission, and that several of the earlier comments had been addressed. Nevertheless, some key points remain.

Essential revisions:

1) The authors must clarify their methods specifically relating to non-randomness of the concentration-invariant responses.

2) The number of cells classified as concentration-invariant seems small and there was some discomfort on whether this population was statistically robust.

3) There must be clarification on what the 'average correlation' of 0.65 refers to: for each odor, or each concentration?

Reviewer #1:

In this revised paper, Roland et al. carry out optical recordings from anterior piriform cortex. They have performed substantial additional analyses of this rich dataset. In doing this they have mostly addressed my concerns from the earlier version of the paper.

The point about subsets of neurons involved in concentration-independent coding has been taken up with further analysis and discussion. The analysis on this is much more thorough.

The current version more completely discusses how the optical findings relate to the electrical recordings by Bolding et al.

The authors don't examine respiration dependence of responses, but do explain that the slow dynamics of GCaMP6s would have made this less useful.

The authors also carry out other useful analyses, including spatial organization of responses, and dynamics of responses.

Overall the study is a valuable characterization of population-wide odor responses in the piriform cortex and reveals interesting features of coding.

Reviewer #2:

In the second submission of the manuscript "Odor concentration-invariant subnetworks in the mouse olfactory cortex", the authors addressed all reviewers' requests and significantly improved the manuscript. I have only one concern.

In the original review, I requested an independent analysis of different features of the of concertation-invariant and generic cell responses in the piriform cortex. The authors presented results on Figure 7—figure supplement 1 Figure. I was hoping that it will reveal some differences between these two classes of cells. However, variability and temporal profiles of concentration-invariant and generic cell responses were undisguisable. This fact raised a concern about statistical significance of the phenomenon. The authors discussed "non-randomness" of concentration invariant responses. This is a very crucial piece of analysis, and the authors should present a better description for their methods, and some discussion of the results and their limitations. For example, if the observed proportion of concertation invariant responses is 11.7% and proportion of the responses due to random sampling is 5.3%, does it mean that the real concentration-invariant pool of cell is significantly smaller than 11.7%? Or what proportion of concentration invariant responses can be missed by chance, etc.?

Reviewer #3:

The authors did perform additional analysis as suggested, but in my opinion, they did not address the main points of criticism. Therefore, I still have concerns regarding the results presented here.

1) The new analysis unfortunately does not fully clarify the points raised related to the concentration invariance. It does show that this group of cells is less concentration variant compared to the entire ensemble of cells imaged. Yet, In the revised version of the manuscript, visual inspection still strongly suggests (Figure 6D) that for a given odor, the correlations across repeats for the same (highest) concentration are higher than correlations across repeats for lower concentrations.

In addition, the correlations across repeats for the same (highest) concentration for a given odor are higher than for the same odor across concentrations. It is unclear whether the average correlation (0.65 +/- 0.10) value given across repeats refers indeed to all odors and all concentrations. Is it the same (0.65 +/- 0.10) for each odor and for each concentration? Plotting side by side, for each odor the correlations across concentration (lowest-highest, lowest-intermediate, intermediate-highest), and respectively the average correlations across repeats for each odor concentration (at low, intermediate and highest concentration) would allow a direct comparison for each stimulus.

2) The small fraction of cells that are classified as concentration invariant, remains a concern, as well as lack of evidence that these cells in the data set indeed constitute a cell type as defined by layer, specific input/output patterns or any other features, besides the difference in observed responses (though proposed as such in the manuscript). It is also unclear what the decoding scheme is for a wider concentration range, except an actual change in perceived odor identity mentioned in the text.

It is indeed important to document that a higher percentage of neurons are called as concentration invariant in the piriform cortex compared to the bulb. However, I'm not convinced that the strong message and title of the manuscript should be focused on a result that summarizes, in the best case scenario, the behavior of ~10% of the responsive neurons in the absence of any additional evidence that these 10% of responsive neurons are doing the job as proposed.

3) The anesthetized vs. awake explanation is not robust. In my opinion, different cells can be differentially affected by the brain state, depending on the local and long range inputs they receive and the strength of corresponding activity patterns.

https://doi.org/10.7554/eLife.26337.016

Author response

[Editors’ note: the author responses to the first round of peer review follow.]

[…] A summary of the key points made by the reviewers is below:

This study looks at large-scale coding using optical recordings and concludesthat a subset of piriform neurons respond in a concentration-invariant way,and can be used for identity coding. The dataset is a potentially useful, broadpopulation readout of piriform cortex activity.

We thank the reviewers for their thoughtful comments and constructive suggestions. The reviewers raise several critical points, which we address by performing extensive additional analyses. These additional tests strongly support and extend our initial conclusions and now provide a much more comprehensive description of our data set.

1) The reviewers were concerned that the identification, analysis, andrelevance of the concentration-invariant neurons was not well supported.

The authors should be able to devise a model that would more completelyutilize the data.

We have performed comprehensive additional tests for the selection of concentration-invariant piriform odor responses. We have substantially extended the analysis of our model and discuss how it can accommodate a larger stimulus range. This new analysis provides strong additional support for the relevance of concentration-invariant subnetworks in piriform cortex.

2) The authors should provide a strong analysis to support theirassertion that the subsample of cells reliably defined odors in aconcentration-invariant way.

We have addressed this concern directly. We now provide additional analysis of the response properties of concentration-invariant piriform neurons, as well as important additional tests to further demonstrate their ability to encode odor identity – independent of odor intensity.

3) There should be a more complete analysis of topography of odor responses.

We have extended our analysis of the topography of piriform odor responses, as suggested by the reviewers. This additional analysis strengthens our initial conclusions.

4) This paper doesn't really draw upon the data from the companion paper,and indeed in some ways doesn't seem to match.

There should be a much better discussion of the effects ofanesthesia, which could compare the two datasets.

The analysis of concentration coding in the two datasets could bedone in a comparative manner, with equivalent analyses in the two cases.

We have added a detailed discussion of the potential effects of anesthesia on piriform odor representations, and we provide additional analysis of concentration coding. Furthermore, we now include an extensive comparison of the two studies and discuss the complementarity of their findings.

5) The authors should provide a better analysis of the dynamics of responses.

We have added a detailed analysis of the dynamics of piriform odor responses in the revised version of the manuscript.

Reviewer #1:

This study looks at large-scale coding using optical recordings and concludesthat a subset of piriform neurons respond in a concentration-invariant way,and can be used for identity coding. These are anesthetized recordings, andit would be difficult to do these in awake animals given the quite invasivesurgical preparation.

I review this in context of the companion manuscript by Bolding and Franks.

What do we learn from the current paper that the Bolding manuscript does not say?

The main point here is to identify a subset of concentration-invariant neurons. I'm not sure that this isn't there in the Bolding dataset, but the authors don't seem to have looked for it. I'm surprised that the classification accuracy is so low (69%), with a few hundred neurons on each trial. I'm also surprised at how slowly the classification accuracy builds up and declines. Possibly this is an artifact of the anesthesia.

Overall, I get the sense that the power of optical recordings hasn't beenfully utilized in this study, and that the analysis could be extended toextract more information from the already obtained data. On its own, and inthe present form, the main result of a few cells being concentration-invariant seems limited.

According to the reviewer’s suggestion we have substantially extended our analysis and Discussion to better bring into focus the novelty and relevance of our findings, and the complementarity of the optical and electrophysiological approaches.

Our study documents, for the first time, the odor coding properties of piriform neurons based on two- photon calcium imaging data. One driving force behind the development of a comprehensive population coding analysis approach was the observation that piriform odor responses exhibited high levels of heterogeneity and trial-to-trial variability. This key characteristic of piriform odor responses had not been recognized in an earlier piriform imaging paper by Stettler and Axel (Neuron 2009), due to technical limitations and the consequent need to average.

As the reviewer points out, an important finding of our study is the identification of functionally distinct piriform neurons with different odor coding properties. The characterization of distinct coding properties in piriform neurons is a major advance over previous studies, which had treated piriform neurons as homogeneous. Furthermore, our study directly addresses a major unresolved and controversial question in the field. Stettler and Axel provocatively stated that odor information is encoded in randomly distributed ensembles of piriform neurons. This claim has remained untested. We have explicitly tested this prediction, and we now provide new additional analysis to support our finding that odor identity can indeed be decoded from distributed piriform response patterns without spatial clustering. Finally, our imaging data analysis required the development of a novel algorithm for automated cell segmentation of GCaMP6-expressing neurons, which we describe here for the first time. While other methods for cell segmentation have been published recently, the lack of a reliable and efficient means to identify and sort densely packed neurons – in piriform and many other brain areas – has remained a major obstacle in the field. We thus believe that, with the critical new analysis we provide in the revised version of this manuscript, the conceptual and technological advances described in this work and the complementarity of our results in comparison with the results obtained from electrophysiological recordings will of broad interest.

1) The key point of the paper is that piriform odor representations change withincreasing concentrations, but a 10% subset of neurons stay consistent andappear to code only for identity. The authors compare this with an analysisof previous data on OB M/T cells, which are mostly concentration-dependent.

While this result is interesting to olfactory physiologists, I was looking fora more detailed analysis, or a mechanistic insight. For example, ifone confines oneself to fewer odors, does the same subset of neurons stayconsistent? Can the authors devise a model that would span the range fromfully concentration-invariant to concentration-dependent for all odors?

We thank the reviewer for this interesting idea. According to the reviewer’s suggestion, we have extended our analysis of concentration invariance. First, we identified concentration-invariant neurons across a 10-fold instead of a 100-fold range in odorant concentration. Using the same selection criteria as for the original data set, this analysis yielded 19.2% ( ± 7.4% SD across experiments) concentration-invariant neurons. Second, we have identified concentration-invariant neurons for pairs of two odorants instead of the three odorants in our test panel. We found that 22% ( ± 8% SD) of the cells were concentration-invariant for at least one pair of odorants. Of those cells, 49% were identified as concentration-invariant for all three odorants while the response of the other 51% of cells was modulated by the concentration of one of the three odorants. We have included this additional analysis in the revised manuscript (Results, subsection “A concentration-invariant piriform subnetwork”, fourth paragraph).

These data suggest that concentration-invariant subnetworks of neurons remain relatively stable within the range of our stimulus set and can provide concentration-invariant odor identity information across a large range of concentration, but that concentration-invariance is not absolute. Concentration-invariant subnetworks depend, to a certain extent, on stimulus complexity and concentration range. This model is consistent with behavioral experiments in humans and insects, which have shown that the perceived odor identity can vary with concentration and stimulus complexity (Bhagavan and Smith, 1997; Engen, 1964; Gross-Isseroff and Lancet, 1988; Laing et al., 2003). We have included a more extensive discussion of odor concentration invariance in the revised manuscript (Discussion, subsection “Odor identity and intensity”, first paragraph). In future experiments, it will be interesting to establish behavioral tests combined with electrophysiological or optical recordings to quantitatively define the boundaries of concentration-invariance. However, these experiments are beyond the scope of this manuscript.

2) I don't see any attempt to relate the optical findings to the electricalrecordings in the companion paper. If one were to filter out the different time-resolutions, would the optical findings simply fall out of the electrical ones? Does the fraction of concentration-independent neurons match?

We thank the reviewer for his/her suggestion to more directly relate our findings to the results obtained from extracellular recordings by Bolding and Franks. We have added additional analyses, and we carefully discuss the similarities and differences of the two studies and their complementarity (Discussion, subsection “Comparison with extracellular recordings in awake mice”).

Bolding and Franks utilize the high temporal resolution of their electrophysiological recordings to propose that information about odor intensity is encoded in the synchrony of piriform odor responses. Our study, on the other hand, focuses on the spatial organization of piriform odor representations and provides evidence for a concentration-invariant subpopulation of piriform neurons, which encodes odor identity – independent of odor intensity. Thus, the two studies propose complementary coding schemes for non-interfering representations of odor identity and odor intensity in the mouse olfactory cortex.

Interestingly, Bolding and Franks report that a subpopulation of ~20% of neurons exhibits concentration-invariant odor response latencies. Differences in the experimental design of the two studies, including the number and concentration range of odorants and the potential bias of extracellular recordings for active neurons preclude a direct comparison of the fraction of concentration-invariant neurons. However, the results of the two studies are compatible and it is tempting to speculate that they indeed represent the same type of piriform neuron.

A difference between the two studies is anesthesia: due to the surgical preparation required to expose piriform cortex for two-photon imaging, optical recordings were performed in ketamine/xylazine-anesthetized mice while electrophysiological recordings were obtained from awake, head-fixed mice. Despite the diverse effects of anesthesia on odor sampling and olfactory processing, many key parameters of piriform odor responses are very similar. The majority of piriform neurons are narrowly tuned to odor, exhibit high trial-to-trial variability to repeated presentations of the same odor, and substantial overlap between response patterns elicited by different odors. Piriform odor representations are less sparse than previously suggested, and odor identity can accurately be decoded from piriform ensembles in the absence of precise temporal information, consistent with a previous report by the Uchida lab (Miura et al., 2012). Importantly, both studies then report that piriform odor representations are not concentration-invariant: odorants presented at a range of different intensities (30-fold concentration range in Bolding and Franks, 100-fold in this study) elicit highly dissimilar response patterns, as dissimilar as the response patterns observed for two different odorants.

3) The frame rate for these recordings is quite good, 15 FPS. I am surprisedthere is no analysis of the dynamics of responses for individual neurons,or to compare them between the concentration dependent and independent cells.

We thank the reviewer for his/her suggestion to explore the dynamics of piriform odor responses. We now include a new analysis of population dynamics over time (Figure 2E). We trained a classifier at defined time points after odor onset, and we measured classification accuracy over time. This analysis shows that odor representations are dynamically rearranged to maintain odor information for several seconds after odor offset.

According to the reviewer’s suggestion, we also provide additional analyses on the dynamics of piriform odor responses. We have characterized response onset, rise time and duration, for concentration-invariant and concentration-dependent cells (Figure 6—figure supplement 1). Our analysis shows that concentration-invariant neurons exhibit response profiles similar to those of other odor-responsive neurons. However, it is important to note that the response dynamics of GCaMP6s (rise time to peak: ~0.5 s, half decay time: >1 s; Chen at el., Nature 2013) are likely to be the limiting factor in our analysis of response dynamics analysis, rather than the frame rate of our recordings.

4) Did the authors record respiration? Given that the animals were anesthetized, it seems likely that at this frame rate it should be possible to examine respiration dependence of the piriform responses.

We obtained reliable recordings of respiration frequency in the majority of our experiments (11 of 14 mice, 15 of 19 imaging sites). We observe an average breath period of 370 ms (Figure 4—figure supplement 1). This breath frequency is very similar to that observed by Bolding and Franks in awake mice. However, while Bolding and Franks observe a small increase of breath frequency at increasing odor concentrations, breath frequency in anesthetized mice did not change following odor delivery at any concentration (breath period, mean ± SD: 0.01%: 0.36 ± 0.05, 0.1%: 0.37 ± 0.06, 1%:

0.37 ± 0.05). We have included this additional information in Figure 4—figure supplement 1 of the revised manuscript.

We also provide additional analysis of population dynamics over time (Figure 2E). However, due to the slow temporal dynamics of the GCaMP6s calcium indicator, and given that respiration in our experiments is independent from odor exposure, we did not examine respiration dependence of piriform odor representations.

5) There are a couple of findings that the authors don't really follow up.

For example, the authors briefly mention functional heterogeneity of PFneuron responses, but they don't seem to follow this up. Also, there isthe finding that the number of responsive Piriform neurons is roughlyindependent of odor concentration.

According to the reviewer’s suggestion, we have added a detailed characterization of the response properties of piriform neurons (temporal dynamics, change in fluorescence, trial-to-trial variability), which further document their functional heterogeneity. We now include additional analysis to describe the concentration-dependence of piriform odor representations and its implications for odor coding in the revised manuscript. Please note that while the number of responsive piriform neurons only moderately changes with odor concentration, the individual response properties of the majority of neurons are odor concentration-dependent. We discuss this in the manuscript.

Reviewer #2:

In their manuscript, 'Odor concentration-invariant subnetworks in the mouse olfactory cortex', Roland et al., investigate how ensemble activity patterns in the olfactory cortex decode odor identity vs. intensity. The authors employ multiphoton imaging in anaesthetized mice and probe GCaMP6s responses to different odors and concentrations. They find that odor identity can be accurately decoded from the population activity using a linear classifier. Across different concentrations, the population representations of odorants change, thus degrading information about odor identity. As a solution to this issue, the authors propose that a subset of concentration-invariant cells in the olfactory cortex is well suited to decode odor identity in a concentration invariant manner.

The study is topical and the experiments are carefully performed. However, I have several major concerns with the interpretation of the results presented here, which, in my opinion, preclude the publication of the manuscript in the current form in eLife.

1) The authors report ~10% concentration-invariant neurons in the PC vs. ~5% in the OB. They perform a shuffling control and suggest that, unlike the OB, the emergence of concentration-invariant cells in the cortex cannot be attributed to chance sampling. Therefore, these neurons represent an emerging feature of the cortex and can solve the identity problem irrespective of changes in concentration.

The fraction of invariant cells is rather small in both the case of cortex and the bulb. The large majority of neurons in the cortex vary substantially in their responses with changes in concentration (Figure 3). Therefore, to this reviewer, asserting that concentration invariance is solved in particular by these neurons is somewhat more anecdotal rather than based on testable evidence. The authors sample 6 times more neurons in the olfactory cortex compared to the bulb (~3,000 vs. 500). Therefore, it is important to determine whether sub-sampling the olfactory cortex data to match the OB recorded cells, and relaxing the significance criterion from p<0.01 to p<0.05 (Figure 6E) change the results presented.

The cross-concentration correlations of odor representations given by even the neurons identified as 'concentration-invariant' are relatively low (0.56, 0.60). What is the average correlation between the neuronal representations of the same odor across repeats, for each of three concentrations sampled? For example, in Figure 5D, for ethyl acetate and hexanone, the correlations across repeats for the highest concentration appear by visual inspection higher than 0.8. This value is substantially higher than the correlations reported between the low-medium, medium-high or low-high concentrations in the sampled range. These observations put into question the concentration-invariant status of these neurons.

Can the authors provide a testable model on how the concentration invariant cells are preserved as such across a wider concentration range (for example 50,000 fold vs. 100 fold sampled), and also how this information is readout selectively?

In addition, in the companion paper (Bolding and Franks), in awake animals, a different mechanism (at the whole population level) is proposed for decoding odor identity, while concentration is read by changes in response latency. Can the authors comment how these two perspectives relate to one another?

We thank the reviewer for bringing up several important points.

First, the reviewer suggests crucial additional tests to strengthen our assertion that concentration- invariant neurons are overrepresented in piriform cortex compared to olfactory bulb. According to the reviewer’s suggestion, we subsampled piriform cortex to balance the numbers of neurons. Furthermore, we tested our model with a relaxed significance criterion of p < 0.05 (instead of p < 0.01). The results strongly support our initial conclusions. We performed 1000 random iterations of subsampling of piriform neurons and obtained in all cases a fraction of concentration-invariant neurons larger than the one observed for the olfactory bulb. We also performed a bootstrap analysis as in Figure. 7D and found in all subsamplings that the observed fraction of concentration invariant neurons was significantly larger than chance (p<0.01). Furthermore, relaxing the selection criterion to p < 0.05 yields 12.7% ± 5.6% of concentration-invariant piriform neurons, significantly more than in the olfactory bulb (6.7 ± 6.1%, Wilcoxon test, p=0.0076). This additional analysis is now presented in Figure 7 – —figure supplement 1 of the revised manuscript.

The reviewer points out that correlations of odor representations by concentration-invariant neurons appear lower across concentrations than across individual trials. The average correlations for the three odorants across individual trials for each of the three concentrations are 0.65 ( ± 0.10). The average correlations for the three odorants across the three concentrations are 0.6 ( ± 0.09) (low to intermediate) and 0.56 ( ± 0.1) (intermediate to high) for the concentration-invariant neurons, compared to 0.53 ( ± 0.12) and 0.36 ( ± 0.12) for the entire piriform ensemble. We now include this additional information and a more detailed discussion of concentration-invariance in the revised version of the manuscript (Results: subsection “A concentration-invariant piriform subnetwork”, Discussion: subsection “Odor identity and intensity”).

Our analysis suggests that odor representations comprised of concentration-invariant neurons are significantly more similar across the 100-fold range of odor concentration than the entire piriform ensemble, and that they provide a computational advantage to identify odors independent of concentration (Figure 5E and F). However, the concentration-invariant subnetworks we identify are not fully concentration-invariant across the tested range in odorant concentration (see also reviewer 1, point 1). The ability to identify an odor across a range of concentrations is essential for olfactory- driven behaviors. However, behavioral studies in humans and insects have also shown that some odors can change perceived odor identity with changing odor concentrations (Bhagavan and Smith, Physiol Behav. 1997; Engen, Ann N Y Acad Sci. 1964; Gross-Isseroff and Lancet, Chem. Senses 1988; Laing et al., 2003; Peltz et al., 1997). Thus, odor representations must provide representations of odor identity that are typically robust to changes in odorant concentration, but concentration- invariance is not absolute. It will be interesting to establish behavioral tests combined with electrophysiological or optical recordings to further examine the boundaries of concentration- invariance. However, these experiments are beyond the scope of this manuscript.

Finally, according to the reviewer’s suggestion, we now discuss how concentration-invariant information encoded in piriform neural subnetworks may be read out by downstream brain areas, and how our model complements the findings by Bolding and Franks, who propose a multiplexed spatial and temporal coding strategy for representing distinct odor features. In brief, we propose that functionally distinct subpopulations of piriform neurons represent distinct piriform neural cell types. Recent work in acute slice preparations has indeed highlighted the diversity of piriform layer II cells. Piriform layer II is composed of semilunar and superficial pyramidal cells, and cells of intermediate phenotype (Suzuki and Bekkers, 2011; Wiegand et al., 2011). Differences in intrinsic properties and synaptic connectivity could explain how different piriform neural cell types can represent different stimulus features. Moreover, layer II semilunar and pyramidal cells project to distinct piriform target areas (Chen et al., 2014; Diodato et al., 2016), providing an opportunity to selectively read out stimulus information. For example, layer II pyramidal cells send cortical feedback projections to the olfactory bulb, which have been implicated in signal normalization – a function that primarily relies on information about stimulus intensity but may be largely independent of odor identity (Boyd et al., 2015; Otazu et al., 2015). On the other hand, semilunar cells projecting to the cortical amygdala, which has been implicated in the encoding of odor valence (Root et al., 2014), may transmit information about odor identity, independent of intensity.

2) Comparison of olfactory cortex and olfactory bulb activity patterns is hard to interpret in anaesthetized mice since recent reports suggest that cortical feedback is strongly affected by anaesthesia (Rothermel et al., 2014, Boyd et al., 2015, Otazu et al. 2015).

We agree with the reviewer’s point that anesthesia affects the activity cortical feedback projections to the olfactory bulb. An important finding of our study is that odor identity independent of odor intensity can be decoded from a functionally distinct subnetwork of piriform neurons, and that concentration- invariant cells are overrepresented in piriform cortex compared to the olfactory bulb. The most parsimonious model for concentration-invariant subnetworks in piriform cortex – and for functionally distinct subpopulations of piriform neurons more generally – is that such neural subpopulations represent distinct piriform neural cell types. Work in acute slice preparations has indeed highlighted the diversity of intrinsic properties and synaptic connectivity of piriform layer II cells (Choy et al., 2015; Johenning et al., 2009; Large et al., 2016; Suzuki and Bekkers, 2011).

We propose that the functional differences of odor responses we observe reflect the response properties of distinct piriform neural cell types. While anesthesia can modulate neural odor responses, it is unlikely that anesthesia imposes functional differences between different neural cell types. Future experiments will be required to directly test the extent to which our model is valid across brain states.

3) Figure 3's title is not supported by the data. Both the percentage of active neurons and their identity change across the sampled concentration range (3B, C). At the level of individual cells, there is substantial variability in lifetime sparseness across concentrations (3D), even though the percentage of neurons that show similar lifetime sparseness values is same across concentrations. In my opinion, such observations invalidate the title's claim that the overall levels of piriform cortex activity remain largely stable across changes in concentration.

We thank the reviewer for pointing out this potential confound. We have changed the figure caption to: Figure 4. Odor-evoked activity and sparseness of individual piriform neurons depends on odorant concentration.

Reviewer #3:

Concentration invariance is known to be a key property of odor identity coding, and the work of Roland et al. is the first (together with a co-submitted article by Bolding and Franks) to study how piriform cortex solves concentration invariance. Furthermore, the authors use a preparation which allows for simultaneous monitoring of large ensembles of spatially-defined neuronal populations, allowing them to examine how concentration-invariance is implemented at the population level.

The authors demonstrate the presence of concentration-invariant neurons, which comprise a subset of piriform neurons. Additionally, a separate sub-population of piriform neurons change their activity as a function of odor concentration. This is an important observation and is worth publishing. However I have a few comments:

1) One of the novel observations is that the concentration-invariant subnetwork lacks topographic organization in piriform cortex. The authors should elaborate more on this point. The case against topographic organization is made in Figure 2—figure supplement 2 by showing similar classification using a local versus random clustering rule. However, this may not rule out that topography exists. (1) It is possible that any single specific topographic mapping may be washed out in the averaging procedure. The authors should perform an additional analysis: artificially simulating topography and performing the exact same analysis to check that this simulated topography actually yields a classification plot that is different from random. (2) Another possible suggestion to examine topography in piriform representations is to examine the spatial distribution of weights in the classifier.

We thank the referee for bringing up this important point. To more quantitatively test the spatial organization of odor-responsive piriform neurons we have developed a statistical test to detect spatial clustering of functionally distinct groups of cells. We first define functionally distinct neurons by their odor tuning, and we compute the nearest neighbor index (NNI) as the mean distance of each cell to its nearest neighbor within the same group. We then test whether the observed NNI is significantly smaller than for neurons whose positions are shuffled to generate the expected NNI distribution for randomly distributed ensembles. Furthermore, we simulate subtle inhomogeneities in the spatial distribution of neurons to validate the sensitivity of this test to detect spatial clustering. This additional analysis is now presented in Figure 3 and fails to detect spatial organization (p > 0.05) in piriform response patterns to 13 different odorants for all 6 imaging sites.

We then use the same model to test the spatial organization of concentration-invariant piriform subnetworks. A random distribution of neurons (p > 0.05) was again observed for 11 out of 13 imaging sites, and only weak clustering was observed for 2 imaging sites. To illustrate these findings we now show results of the statistical tests and the distributions of concentration-invariant neurons for each imaging sites in Figure 6 —figure supplement 1.

We favor to keep our analysis of odor classification performance across increasingly large patches of cells as a complementary, more qualitative illustration of the fact that odor identity information is broadly distributed across piriform ensembles. Please note that homogenous classification performance can be obtained even with spatial clustering of odor identity preference, thus we do not suggest this analysis to be a proof of the lack of spatial organization. Finally, we considered the idea of using classifier weights as a test for statistical significance. However, with 13 odorants and one classifier for each odorant pair it is difficult to meaningfully interpret the distribution of the 13 corresponding sets of weights.

2) The authors made a point about the trial-to-trial variability of piriform neurons (e.g. Figure 1E). A pertinent question is, what is the trial-to-trial variability for neurons in the concentration invariant subnetwork?

According to the reviewer’s suggestion we have determined trial-to-trial variability of concentration- invariant neurons (Figure 6—figure supplement 1). We do not observe a significant difference in trial- to-trial variability between concentration-invariant and ‘generic’ neurons.

3) The author should discuss the effect of anesthesia in their preparations. They compared their research with electrophysiological study co-submitted by Bolding and Franks. The current study revealed the spatial distributions of identity encoding neurons. The authors should discuss the limitations of their approach and clarify what other conclusions can be safely extrapolated for awake case?

We thank the reviewer for this suggestion. We have included a detailed discussion of potential anesthesia effects, and a direct comparison with the results obtained from electrophysiological recordings in awake mice by Bolding and Franks in the revised manuscript.

As the reviewer points out, an important difference between the two studies is anesthesia: due to the surgical preparation required to expose piriform cortex for two-photon imaging, optical recordings were performed in ketamine/xylazine-anesthetized mice while electrophysiological recordings were obtained from awake, head-fixed mice. Despite the diverse effects of anesthesia on odor sampling and olfactory processing (see below), many key parameters characteristic of piriform odor responses are very similar. Piriform odor representations are relatively sparse, and the majority of piriform neurons are narrowly tuned to odor. Importantly, piriform odor representations exhibit high trial-to-trial variability to repeated presentations of the same odor and substantial overlap between response patterns elicited by different odors. Finally, and consistent with an earlier report (Miura et al., 2012), odor identity can accurately be decoded from piriform ensembles in the absence of precise temporal information. Both studies then report that piriform odor representations are not concentration-invariant: odorants presented at a range of different intensities (30-fold concentration range in Bolding and Franks, 100-fold in this study) elicit highly dissimilar response patterns, as dissimilar as the response patterns observed for two different odorants.

Bolding and Franks utilize the high temporal resolution of their electrophysiological recordings to propose that information about odor intensity is encoded in the synchrony of piriform odor responses. Our study, on the other hand, focuses on the spatial organization of piriform odor representations and provides evidence for a concentration-invariant subnetwork of piriform neurons that encodes odor identity – independent of odor intensity. Thus, the two studies propose complementary coding schemes for non-interfering representations of odor identity and odor intensity in the mouse olfactory cortex. Interestingly, Bolding and Franks also observe a subpopulation of neurons that exhibit short- latency concentration-invariant odor responses. While differences in the experimental design of the two studies preclude a direct comparison of these neural subpopulations, it is tempting to speculate that they represent the same type of piriform neuron.

Anesthesia interferes with active sniffing, and awake mice and rats generally exhibit faster sniff rates that can be dynamically modulated by odor (Blauvelt et al., 2013; Carey and Wachowiak, 2011; Wachowiak et al., 2013). In contrast, we do not observe significant modulation of sniff rate in ketamine/xylazine-anesthetized mice. Anesthesia has also been shown to modulate the activity of neurons in the olfactory bulb. In awake mice, periglomerular and granule cell inhibitory neurons exhibit higher levels of activity (Cazakoff et al.et al., 2014; Kato et al., 2012; Wachowiak et al., 2013), while some mitral cells exhibit diminished spontaneous and odor-evoked activity (Kato et al., 2012; Kollo et al., 2014; Shusterman et al., 2011). Tantirigama et al. (2017) report that ketamine/xylazine anesthesia affects baseline activity of piriform neurons. In contrast, unpublished data by Bolding and Franks (manuscript in preparation) suggests that ketamine/xylazine anesthesia has no significant effect on the fraction of odor-responsive piriform neurons or their odor tuning (see Author response image 1). However, anesthesia does alter sniff coupling and the temporal response properties of piriform neurons (Bolding and Franks, companion paper and Author response image 1). Of note, at elevated odorant concentrations, extracellular recordings in awake mice show that the fraction of activated neurons remains stable, while optical recordings in anesthetized mice show a moderate increase. This observation suggests that normalization of neural activity is incomplete under anesthesia, consistent with recent reports on the activity of cortical feedback projections to the olfactory bulb (Boyd et al., 2015; Otazu et al., 2015; Rothermel and Wachowiak, 2014). The activity of cortical feedback projections has been suggested to contribute to signal normalization and is attenuated under anesthesia. We now discuss this distinction in our paper.

Author response image 1
Odor responses are similar in awake and anesthetized mice.

(A) Example showing respiration and spontaneous activty in a mouse before (black) and during ketamine/xylazine (k/x, cyan) anesthesia. (B) Example PSTHs in response to 3 odors before and during anesthesia. Each row is a different cell. Responses are aligned to inhalation. (C) Population and lifetime sparseness, averaged across experiments, are similar in awake and during k/x anesthesia. (D) Distribution of peak firing rates are similar in awake and anesthetized mice. (E) Responses awake mice occur earlier than under anesthesia. (F) Response durations are similar in awake and anesthetized mice.

https://doi.org/10.7554/eLife.26337.015

The most parsimonious model for concentration-invariant subnetworks in piriform cortex – and for functionally distinct subpopulations of piriform neurons more generally – is that such neural subpopulations represent distinct piriform neural cell types. Work in acute slice preparations has indeed highlighted the diversity of intrinsic properties and synaptic connectivity of piriform layer II cells (Choy et al., 2015; Johenning et al., 2009; Large et al., 2016; Suzuki and Bekkers, 2011). However, experiments investigating the in vivo odor response properties of piriform neurons have treated these different cell types as homogeneous (Miura et al., 2012; Poo and Isaacson, 2009; Rennaker et al., 2007; Stettler and Axel, 2009), and the response properties of distinct piriform principle cells to odor remain unknown We favor the model that the functional differences we observe reflect the response properties of distinct piriform neural cell types and are thus unlikely to be imposed by anesthesia. Future experiments will be required to directly test the extent to which this model is valid across brain states.

[Editors' note: the author responses to the re-review follow.]

Essential revisions:

1) The authors must clarify their methods specifically relating to non-randomness of the concentration-invariant responses.

We have added a comprehensive description of the statistical analysis of concentration-invariant neurons to the Methods section of the revised manuscript. We have also included a detailed discussion of the results and limitations of this analysis.

2) The number of cells classified as concentration-invariant seems small and there was some discomfort on whether this population was statistically robust.

There are essentially two concerns here: the seemingly small size of the concentration-invariant cells and the robustness of the statistics used to identify them.

Regarding the size of the concentration-invariant subpopulation, it is important to bear in mind that in our data set, ~40% of the cells responded selectively to any of the three odorants. Thus, the ~10% of cells that have concentration-invariant responses represent almost one third of all odor-selective cells. While this was described in the Methods section and further explored in the additional analysis of concentration-invariant subpopulations requested by reviewer 1, it was not explicitly stated in the text. We have now clarified this point in the Results section and Discussion of the revised manuscript.

We have also extended the analysis and discussion of the statistics used for the identification of concentration-invariant piriform neurons. Importantly, cross-correlations, PCA and linear classification analysis consistently support the robustness of our conclusions.

3) There must be clarification on what the 'average correlation' of 0.65 refers to: for each odor, or each concentration?

We apologize for confusion here. We now provide a detailed description of the cross-correlation analysis in Figure 6—figure supplement 1. Following reviewer 3’s recommendation, we have plotted, side by side, the cross-correlation values within and across odorant concentrations for each odorant. We see how this helps clarify the issue and thank the reviewer for this suggestion.

Reviewer #2:

In the second submission of the manuscript "Odor concentration-invariant subnetworks in the mouse olfactory cortex", the authors addressed all reviewers' requests and significantly improved the manuscript. I have only one concern.

In the original review, I requested an independent analysis of different features of the of concertation-invariant and generic cell responses in the piriform cortex. The authors presented results on Figure 7—figure supplement 1. I was hoping that it will reveal some differences between these two classes of cells. However, variability and temporal profiles of concentration-invariant and generic cell responses were undisguisable. This fact raised a concern about statistical significance of the phenomenon. The authors discussed "non-randomness" of concentration invariant responses. This is a very crucial piece of analysis, and the authors should present a better description for their methods, and some discussion of the results and their limitations. For example, if the observed proportion of concertation invariant responses is 11.7% and proportion of the responses due to random sampling is 5.3%, does it mean that the real concentration-invariant pool of cell is significantly smaller than 11.7%? Or what proportion of concentration invariant responses can be missed by chance, etc.?

We appreciate the reviewer’s suggestion to provide a comprehensive description of the statistical methods used to characterize concentration-invariant piriform neurons. We have added a detailed explanation of the bootstrap analysis used for the characterization of concentration-invariant subpopulations of piriform neurons in the Methods section of the revised manuscript. We also discuss the results obtained from this analysis and its limitation (subsection “Analysis of variance of single neurons and bootstrap methods”).

Reviewer #3:

The authors did perform additional analysis as suggested, but in my opinion, they did not address the main points of criticism. Therefore, I still have concerns regarding the results presented here.

1) The new analysis unfortunately does not fully clarify the points raised related to the concentration invariance. It does show that this group of cells is less concentration variant compared to the entire ensemble of cells imaged. Yet, In the revised version of the manuscript, visual inspection still strongly suggests (Figure 6D) that for a given odor, the correlations across repeats for the same (highest) concentration are higher than correlations across repeats for lower concentrations.

In addition, the correlations across repeats for the same (highest) concentration for a given odor are higher than for the same odor across concentrations. It is unclear whether the average correlation (0.65 +/- 0.10) value given across repeats refers indeed to all odors and all concentrations. Is it the same (0.65 +/- 0.10) for each odor and for each concentration? Plotting side by side, for each odor the correlations across concentration (lowest-highest, lowest-intermediate, intermediate-highest), and respectively the average correlations across repeats for each odor concentration (at low, intermediate and highest concentration) would allow a direct comparison for each stimulus.

We thank reviewer 3 for pointing out additional points for clarification.

Regarding the correlation matrices presented in Figures 5 and 6, it is important to keep in mind two principles. First, this presentation of the data is not aimed at assessing the discriminability of representations (this must be done with cross-validated classifiers, as we did in Figure 6F) but rather is presented as a way to qualitatively visualize the variability of population responses. Second, as a rule of thumb (for details, see Bathellier et al., Neuron 2012), one can consider the response patterns for two stimuli significantly distinct from each other if their correlations are lower than the trial-to-trial correlations of repeated presentations of either stimulus alone.

As suggested by the reviewer, we have plotted, side by side, the cross-correlations across trials and concentrations for each odorant (Figure 6—figure supplement 1). Indeed, this presentation helps clarify our conclusion that, for all three odorants, responses of concentration-invariant piriform neurons are highly correlated across trials and concentrations.

Despite an apparent trend (which disappears when using concentration-invariant subpopulations identified with p < 0.05 as the significance criterion), there is no statistically significant difference (Wilcoxon ranked sum test p > 0.05, see new Figure 6—figure supplement 1) between trial-to-trial correlations for responses to odorants at low (0.01%) versus high (1%) concentrations. Furthermore, we tested if trial-to-trial correlations at low odorant concentrations were different from the low versus high odorant concentration correlations. We found no significant differences for acetophenone and hexanone, and only a small but significant difference for ethyl acetate (0.57 ± 0.04% SD versus 0.65 ± 0.05% SD, p = 0.007).

This does not alter the main message of Figure 6D, which is that the response correlations across concentrations in concentration-invariant subpopulations are strikingly different from those of the entire ensemble, which strongly decorrelate with increasing concentrations (Figure 5C and D). We agree that there might be small differences with some odors between the representations of low and high concentrations. Crucially, these differences are too small to make the responses for low and high concentration discriminable on a single trial basis, as we show in Figure 6F. We have included a more detailed description of these data in Results section of the revised manuscript (subsection “A concentration-invariant subpopulation of piriform neurons”, second paragraph, and Figure 6—figure supplement 1).

2) The small fraction of cells that are classified as concentration invariant, remains a concern, as well as lack of evidence that these cells in the data set indeed constitute a cell type as defined by layer, specific input/output patterns or any other features, besides the difference in observed responses (though proposed as such in the manuscript). It is also unclear what the decoding scheme is for a wider concentration range, except an actual change in perceived odor identity mentioned in the text.

It is indeed important to document that a higher percentage of neurons are called as concentration invariant in the piriform cortex compared to the bulb. However, I'm not convinced that the strong message and title of the manuscript should be focused on a result that summarizes, in the best case scenario, the behavior of ~10% of the responsive neurons in the absence of any additional evidence that these 10% of responsive neurons are doing the job as proposed.

Given that only 10% of cells (out of ~1 million cells in piriform cortex) are identified as concentration-invariant and we have not yet identified additional identifying features of these cells, the reviewer is uncomfortable with us using a title that focuses exclusively on the concentration-invariant subpopulation. In fact, we agree that this is just one of multiple novel results that we present in this manuscript. Please note we that here we also provide for the first time, a detailed analysis of trial-to-trial variability within and across concentration for the general population, quantify normalization of piriform responses across concentrations, and demonstrate that there is no discernible topography in piriform representations. We have therefore made the title of this manuscript more generic (“Odor identity coding by distributed ensembles of neurons in the mouse olfactory cortex”), per reviewer 3’s suggestion.

3) The anesthetized vs. awake explanation is not robust. In my opinion, different cells can be differentially affected by the brain state, depending on the local and long range inputs they receive and the strength of corresponding activity patterns.

We agree with the reviewer that different piriform cell types could differentially be affected by anesthesia. We have extended our discussion of potential anesthesia effects to include this point (subsection “Comparison with extracellular recordings in awake mice”, first and last paragraphs).

https://doi.org/10.7554/eLife.26337.017

Article and author information

Author details

  1. Benjamin Roland

    Center for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS UMR 7241, INSERM U1050, Paris, France
    Contribution
    BR, Data curation, Formal analysis, Investigation, Methodology, Writing—review and editing
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3413-1044
  2. Thomas Deneux

    Unité de Neuroscience, Information et Complexité, Centre National de la Recherche Scientifique, UPR 3293, Gif-sur-Yvette, France
    Contribution
    TD, Data curation, Formal analysis, Methodology, Writing—review and editing
    Competing interests
    The authors declare that no competing interests exist.
  3. Kevin M Franks

    Department of Neurobiology, Duke University, Durham, United States
    Contribution
    KMF, Conceptualization, Supervision, Writing—review and editing
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6386-9518
  4. Brice Bathellier

    Unité de Neuroscience, Information et Complexité, Centre National de la Recherche Scientifique, UPR 3293, Gif-sur-Yvette, France
    Contribution
    BB, Conceptualization, Data curation, Formal analysis, Methodology, Writing—review and editing
    For correspondence
    bathellier@unic.cnrs-gif.fr
    Competing interests
    The authors declare that no competing interests exist.
  5. Alexander Fleischmann

    Center for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS UMR 7241, INSERM U1050, Paris, France
    Contribution
    AF, Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing—original draft, Writing—review and editing
    For correspondence
    alexander.fleischmann@college-de-france.fr
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7956-9096

Funding

European Molecular Biology Organization (ASTF 395 - 2014)

  • Benjamin Roland

LabEx Memolife

  • Benjamin Roland

National Institute on Deafness and Other Communication Disorders (DC009839)

  • Kevin M Franks

National Institute on Deafness and Other Communication Disorders (DC015525)

  • Kevin M Franks

Agence Nationale de la Recherche (SENSEMAKER)

  • Brice Bathellier

Human Frontier Science Program (CDA-0064-2015)

  • Brice Bathellier

Marie Curie Program (CIG 334581)

  • Brice Bathellier

Marie Curie International Reintegration Grant (IRG 276869)

  • Alexander Fleischmann

Fondation pour la Recherche Médicale (AJE201106)

  • Alexander Fleischmann

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Ludovic Cacheux for help with imaging experiments and data analysis, Dmitri Bryzgalov and Simon Daste for help with data analysis, Yves Dupraz for his work on the in vivo imaging set-up, and Jérémie Teillon and Philippe Mailly for help with imaging and imaging data pre-processing. We thank Kevin Bolding, Thomas Preat, Andreas Schaefer, German Sumbre and Jonathan Touboul for critical comments of the manuscript. This work was supported by a Marie Curie International Reintegration grant (IRG 276869), and the ‘Amorçage de jeunes équipes’ program (AJE201106) of the Fondation pour la Recherche Médicale (to AF), an EMBO short term fellowship (ASTF 395–2014) and a postdoctoral fellowship by the LabEx ‘MemoLife’ (to BR), by grants from the NIDCD (DC009839 and DC015525) to KMF), grants from the Agence Nationale pour la Recherche (ANR ‘SENSEMAKER’), the Marie Curie Program (CIG 334581), and The International Human Frontier Science Program Organization (CDA-0064–2015) (to BB).

Ethics

Animal experimentation: This study was performed in strict accordance with French National and INSERM animal care and use committee guidelines (#B750512/00615.02). All surgery was performed under ketamine/xylazine anesthesia.

Reviewing Editor

  1. Upinder Singh Bhalla, National Centre for Biological Sciences, Tata Institute of Fundamental Research, India

Publication history

  1. Received: February 24, 2017
  2. Accepted: April 29, 2017
  3. Accepted Manuscript published: May 10, 2017 (version 1)
  4. Version of Record published: May 19, 2017 (version 2)

Copyright

© 2017, Roland 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

  • 4,874
    Page views
  • 1,064
    Downloads
  • 72
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Benjamin Roland
  2. Thomas Deneux
  3. Kevin M Franks
  4. Brice Bathellier
  5. Alexander Fleischmann
(2017)
Odor identity coding by distributed ensembles of neurons in the mouse olfactory cortex
eLife 6:e26337.
https://doi.org/10.7554/eLife.26337

Further reading

    1. Neuroscience
    Miriam Menzel, David Gräßel ... Marios Georgiadis
    Research Article Updated

    Disentangling human brain connectivity requires an accurate description of nerve fiber trajectories, unveiled via detailed mapping of axonal orientations. However, this is challenging because axons can cross one another on a micrometer scale. Diffusion magnetic resonance imaging (dMRI) can be used to infer axonal connectivity because it is sensitive to axonal alignment, but it has limited spatial resolution and specificity. Scattered light imaging (SLI) and small-angle X-ray scattering (SAXS) reveal axonal orientations with microscopic resolution and high specificity, respectively. Here, we apply both scattering techniques on the same samples and cross-validate them, laying the groundwork for ground-truth axonal orientation imaging and validating dMRI. We evaluate brain regions that include unidirectional and crossing fibers in human and vervet monkey brain sections. SLI and SAXS quantitatively agree regarding in-plane fiber orientations including crossings, while dMRI agrees in the majority of voxels with small discrepancies. We further use SAXS and dMRI to confirm theoretical predictions regarding SLI determination of through-plane fiber orientations. Scattered light and X-ray imaging can provide quantitative micrometer 3D fiber orientations with high resolution and specificity, facilitating detailed investigations of complex fiber architecture in the animal and human brain.

    1. Computational and Systems Biology
    2. Neuroscience
    Jamie D Costabile, Kaarthik A Balakrishnan ... Martin Haesemeyer
    Research Article

    Brains are not engineered solutions to a well-defined problem but arose through selective pressure acting on random variation. It is therefore unclear how well a model chosen by an experimenter can relate neural activity to experimental conditions. Here we developed 'Model identification of neural encoding (MINE)'. MINE is an accessible framework using convolutional neural networks (CNN) to discover and characterize a model that relates aspects of tasks to neural activity. Although flexible, CNNs are difficult to interpret. We use Taylor decomposition approaches to understand the discovered model and how it maps task features to activity. We apply MINE to a published cortical dataset as well as experiments designed to probe thermoregulatory circuits in zebrafish. MINE allowed us to characterize neurons according to their receptive field and computational complexity, features which anatomically segregate in the brain. We also identified a new class of neurons that integrate thermosensory and behavioral information which eluded us previously when using traditional clustering and regression-based approaches.