article: The olfactory system faces the difficult task of identifying an enormous variety of odors independent of their intensity. Primacy coding, where the odor identity is encoded by the receptor types that respond earliest, might provide a compact and informative representation that can be interpreted efficiently by the brain. In this paper, we analyze the information transmitted by a simple model of primacy coding using numerical simulations and statistical descriptions. We show that the encoded information depends strongly on the number of receptor types included in the primacy representation, but only weakly on the size of the receptor repertoire. The representation is independent of the odor intensity and the transmitted information is useful to perform typical olfactory tasks with close to experimentally measured performance. Interestingly, we find situations in which a smaller receptor repertoire is advantageous for discriminating odors. The model also suggests that overly sensitive receptor types could dominate the entire response and make the whole array useless, which allows us to predict how receptor arrays need to adapt to stay useful during environmental changes. Taken together, we show that primacy coding is more useful than simple binary and normalized coding, essentially because the sparsity of the odor representations is independent of the odor statistics, in contrast to the alternatives. Primacy coding thus provides an efficient odor representation that is independent of the odor intensity and might thus help to identify odors in the olfactory cortex.
  The olfactory system identifies and discriminates odors for solving vital tasks like navigating the environment, identifying food, and engaging in social interactions. These tasks are complicated by the enormous variety of odors, which vary in composition and in the concentrations of their individual molecules. In particular, the olfactory system needs to separately recognize the odor identity (what is there?) and the odor intensity (how much is there?). For instance, the identity is required to decide whether to approach or avoid an odor source, whereas the intensity information is important for localizing it. It is not understood how these two odor properties are separated and how odors are discriminated reliably. Odors are comprised of chemicals that bind to and excite olfactory receptors in the nose in mammals and on antenna in insects. Each receptor responds to a wide range of odors and each odor activates many receptor types. The resulting combinatorial code allows to distinguish odor identities [1–3], but also depends on the odor intensity, since receptors respond stronger to more concentrated molecules [4]. To separate these two properties, the neural signals are processed in the olfactory bulb (antennal lobe in insects) and then forwarded to the olfactory cortex, where odors are identified by comparing to memorized patterns. Indeed, experiments indicate that the olfactory cortex receives a concentration-invariant code [5,6], which allows to identify odors irrespective of their intensity. Consequently, the olfactory bulb can be thought of as a signal processor that removes statistical redundancies in the input to provide a more useful representation to the olfactory cortex. However, so far it is not clear what processing the olfactory bulb performs and how this affects odor representations. The olfactory bulb contains neural clusters called glomeruli, which each receive input from a specific receptor type [7–9]. Each glomerulus excites associated projection neurons, which project into the olfactory cortex. Additionally, the glomeruli are connected by local neurons [10,11], which inhibit the projection neurons [12–18]. These local neurons could mediate a global normalization resulting in an intensity-invariant representation of the odor identity [19,20]. However, we showed that simple normalized representations still depend strongly on the number of ligands in a mixture and might thus not be optimal for solving olfactory tasks [21]. An alternative to these normalized representations is rank coding, where the order in which the receptors are excited is used to encode the odor identity robustly and independently of the odor intensity [22]. Indeed, experiments suggests that odors are encoded robustly by the receptor types that respond within a given time window after sniff onset [23–25]. In particular, the odor identity could be robustly encoded by a fixed number of the receptors that respond first, which is known as primacy coding [24,26]. So far, it is unclear how efficient and useful primacy coding is and how it compares to alternative schemes. In this paper, we consider a simple model of primacy coding and investigate how well it represents complex odors. In particular, we identify how much information is transmitted to the cortex and how well this information can be used to perform typical olfactory tasks, like identifying the addition of a target odor to a background or discriminating odor mixtures. Our statistical approach allows linking parameters of the primacy code to results from typical psychophysical experiments. We show that primacy coding provides a robust and compact representation of the odor identity over a wide range of odors, independent of the odor intensity, and that it outperforms other simple coding schemes. However, this good performance of the olfactory system hinges on tuned receptor sensitivities, which suggests that there is a strong selective pressure to adjust the sensitivities on evolutionary and shorter timescales.
  Odors are detected by an array of receptors in the nasal cavity in mammals and on the antenna in insects. The receptor array consists of NR different receptor types, which each are expressed many times. Typical numbers are NR ≈ 50 in flies [7], NR ≈ 300 in humans [31], and NR ≈ 1000 in mice [32]. The excitations of all receptors of the same type are accumulated in an associated glomerulus in the olfactory bulb in mammals and the antennal lobe in insects [33]. Since this convergence of the neural information mainly improves the signal-to-noise ratio, we here capture the excitation of the receptors at the level of glomeruli; see Fig 1A. The excitation en of glomerulus n can be approximated by a linear function of the odor c [4,34], e n = ∑ i = 1 N L S n i c i, (1) where Sni denotes the effective sensitivity of glomerulus n to ligand i. Note that Sni is proportional to the copy number of receptor type n if the response from all individual receptors is summed [30]. The sensitivity matrix Sni could in principle be determined by measuring the response of each glomerulus to each possible ligand. However, because the numbers of receptor and ligand types are large, this is challenging and only parts of the sensitivity matrix have been measured, e. g. , in humans [35] and flies [36]. Using these data, we showed that the measured matrix elements are well described by a log-normal distribution with a standard deviation λ ≈ 1 of the underlying normal distribution [30]. Motivated by these observations, we here consider random sensitivity matrices, where each element Sni is chosen independently from the same log-normal distribution, which is parameterized by its mean 〈 S n i 〉 = S ¯ and variance var (S n i) = S ¯ 2 (e λ 2 - 1) with λ = 1 [30]. Since these receptor sensitivities are broadly distributed, they might not include specific receptors related to innate behavior [37], but they can collectively discriminate concentration differences of several orders of magnitude [30]. The odor representation on the level of glomeruli excitations en depends strongly on the odor intensity, which is quantified by the total concentration ctot = ∑i ci. This dependency complicates the extraction of the odor identity, which is determined by the ligands which are present and their relative concentrations. A concentration-invariant representation could be achieved by normalizing the excitations by their mean [16], which leads to an efficient neural representation on the level of projection neurons [21]. However, recent experimental data suggest an alternative encoding based on the timing of the glomeruli excitation [24]. The key idea of this primacy coding is that the set of receptor types that are excited first is independent of the total concentration ctot and thus provides a concentration-invariant representation. In our simple model of primacy coding, odors are encoded by the identity of the NC glomeruli that respond first. For simplicity, we here neglect the order in which they respond, in contrast to rank coding [22], and we also consider the simple situation where ligands binding to receptors only affect the magnitude of the receptor output, but not the signaling dynamics. In this case, receptors that respond first are also the ones with the largest excitation, so that the primacy code is given by the identity of the NC glomeruli with the largest excitation, which is known as the primacy set [38]. The primacy set can be represented by a binary activity vector a = (a 1, a 2, …, a N R), where an = 1 implies that glomerulus n belongs to the primacy set and is active, while an = 0 denotes an inactive glomerulus not belonging to the primacy set. Since the active glomeruli have the highest excitation, they can be identified using an excitation threshold γ; see Fig 1B. Consequently, the activities are given by a n = { 1 e n > γ (e) 0 e n ≤ γ (e). (2) Physiologically, the activities an could be encoded by projection neurons in insects and mitral and tufted cells in mammals. These neurons receive excitatory input from one glomerulus [39] and are inhibited by a local network of granule cells [20,33]. These granule cells basically integrate the activity of all glomeruli [40] and could inhibit the glomeruli once a threshold is reached. Taken together, this would implement primacy coding since only the glomeruli that respond earliest would be activated. For simplicity, we consider the case where the number NC of active glomeruli is fixed and does not depend on the odor c. The associated constraint N C = ∑ n = 1 N R a n (3) determines the threshold γ. Note that the activity pattern a is sparse since only a fraction NC/NR of all glomeruli is activated. Moreover, a is independent of S ¯ and ctot, implying concentration-invariance. This is because multiplying the concentration vector c by a factor changes the excitations en and the threshold γ by the same factor, so that a given by Eq (2) is unaffected. In essence, only relative excitations are relevant for our model of primacy coding. In the binary representation given by Eq (2), each receptor type can at most contribute 1 bit of information to the odor representation. This worst-case scenario corresponds to large processing noise, such that intermediate excitations cannot be identified in the downstream processing. In fact, the concentration range over which receptors are sensitive is typically small compared to the expected range of odorant concentrations [41]. Consequently, receptors will be activated either very little or very strongly for natural odors, suggesting a binary picture. Moreover, there is evidence that the identity of active neurons can robustly encode odor identity and is actually used in the olfactory system [5]. To see whether primacy coding encodes odor information efficiently [42], we quantify the amount of information I that can be learned about the odor c by observing the binary activity pattern a with given sparsity NC/NR. Since our model is deterministic, I is given by the entropy I = - ∑ a P (a) log 2 P (a), (4) where the probability P (a) of observing an output a depends on the odor environment Penv (c) as well as the properties of the olfactory system, which in our model are quantified by NC, NR, and λ. Since further processing in the downstream regions of the brain can only reduce the amount of information, Eq (4) provides an upper bound for the information that the brain receives about odors when primacy coding as described by Eqs (1) – (3) is used. In an optimal receptor array, each output a occurs with equal probability when encountering odors distributed according to Penv (c) [30]. In our model, only outputs with exactly NC active receptor types are permissible. The resulting representations would be optimal if each receptor type was activated with a probability 〈an 〉 = NC/NR and all types were uncorrelated, cov (an, am) = 0 for n ≠ m. The associated information I max (N C, N R) = log 2 (N R N C) ≈ N C - 1 ln 2 + N C log 2 N R N C (5) provides an upper bound for I given by Eq (4). Here, the approximation on the right hand side is obtained using Stirling’s formula for large receptor repertoires (NR ≫ NC). We start by analyzing the information I transmitted by the primacy code using numerical ensemble averages of Eqs (1) – (4); see Methods and Models. Fig 2A shows that I is very close to the maximal information Imax given by Eq (5), which is obtained when all receptor types have equal activity and are uncorrelated [30]. This indicates that the primacy code uses the different receptor types with similar frequency and that correlations between them are negligible. The expression for Imax implies that the information grows linearly with the primacy dimension NC, but only logarithmically with the number NR of receptor types. Consequently, the number of distinguishable signals, NS = 2I, grows strongly with NC, but the dependence on the repertoire size is weaker; see Fig 2B. Given equal NC, our model thus predicts that the transmitted information in mice is only twice that of flies, although mice possess about 20 times as many receptor types. However, the number of discriminable signals changes by many orders of magnitudes, since it scales exponentially with I. The logarithmic scaling of the transmitted information I with the receptor repertoire size NR could explain why the ability of rats to discriminate odors is not significantly affected when half the olfactory bulb is removed in lesion experiments [43,44]. If this operation removes half the receptor types, our model implies that the transmitted information I is lowered by NC bits; see Eq (5). This corresponds to a reduction of I by about 10% in rats where NR ≈ 1000; see Fig 2C. Conversely, the transmitted information decreases by almost 50% in flies, which have a much smaller receptor repertoire of NR ≈ 50. Our model thus predicts that lesion experiments have a much more severe effect on the performance of animals with smaller receptor repertoires. Taken together, this first analysis already suggests that the primacy code provides a robust odor representation, which is sparse, concentration-invariant, and depends only weakly on the details of the receptor array. However, for this representation to be useful to the animal, it needs to allow solving typical olfactory tasks. Typical olfactory tasks include detecting a ligand in a distracting background, detecting the addition of a ligand to a mixture, as well as discriminating similar mixtures. All these tasks involve discriminating odors with common ligands, implying that the associated primacy sets are correlated. This correlation can be quantified by the expected Hamming distance d between the primacy sets, which counts the number of glomeruli with different activities. The probability η that this distance is larger than 0, so that the two odor representations can be discriminated in principle, is given by η (d) ≈ 1 - (1 - d 2 N C) N C; (6) see Methods and Models. In particular, discriminating similar odors will be impossible (η = 0) if their primacy sets are identical (d = 0). We showed that primacy coding contains sufficient information to perform typical olfactory tasks with experimentally measured accuracy. Although this provides some support for primacy coding, alternative encoding schemes might also be consistent with experimental data. To elucidate this, we next compare primacy coding to two alternatives, which are also based on the simple model described by Eqs (1) and (2). The first alternative is binary coding, where glomeruli become active when their excitation exceeds a constant threshold γ [30,51]. The second alternative is normalized coding, where the threshold is proportional to the mean excitation, γ = α N R - 1 ∑ n e n, and the inhibition strength α determines how many glomeruli are active on average [21]. To see how binary and normalized coding compare to primacy coding, we calculate the probability η that adding a ligand to a mixture of s ligands can be detected; see Fig 6A. The binary code strongly depends on the overall concentration of the presented odor (or, equivalently, the imposed threshold γ). This implies that there is only a narrow region of mixture sizes s where the binary code allows detecting the addition of a ligand. Conversely, the normalized code is concentration-invariant and could thus in principle discriminate odors at all intensities. However, we showed in Ref. [21] that the encoded information and the discriminability still depend strongly on the mixtures size s in this model. Consequently, normalized codes can only discriminate mixtures of realistic sizes when the inhibition strength α is very low and thus many glomeruli get activated on average. The example of the normalized code shows that it is not sufficient to study how well different coding schemes can solve olfactory tasks, but one also needs to consider how useful this code is to the downstream decoder. Without modeling the decoder in detail, we here just propose that sparser codes are preferable since they imply fewer firing neurons, which saves energy and simplifies the downstream processing. In fact, sparse coding is typical for sensory information [19,52]. In our model, the sparsity is given by the fraction 〈an〉 of activated glomeruli and Fig 6B shows this quantity as a function of the mixture size s. Since larger mixtures imply a higher odor intensity, this number increases quickly in binary coding and makes this code inefficient. Conversely, the number of active glomeruli decreases strongly in normalized coding [21], which explains the poor discriminatory performance for large mixtures. In contrast, primacy coding has a constant sparsity, because it is directly controlled by the primacy count NC. Taken together, primacy coding outperforms both binary coding and normalized coding essentially because the sparsity of the representation is independent of the presented odors and can thus be adjusted to be useful and efficient over the whole range of possible odors. The three models discussed here differ in how the statistics of the odor c affect the statistics of the output a. In the binary model, the odor intensity given by the mean concentration ctot affects the mean excitation 〈en〉 and therefore the sparsity 〈an〉. This clearly prevents the response from being useful over a wide concentration range. This dependence on the odor intensity is removed in normalized coding, but the variance of the excitations en still depends on the odor statistics, e. g. larger mixtures imply smaller variations in en. This is problematic since it implies the excitations of fewer glomeruli exceed the fixed threshold in normalized coding, so the sparsity 〈an〉 and the usefulness decrease [21]. In primacy coding, however, the mean activity 〈 a n 〉 = N C N R is independent of the odor statistics, so the system is useful in all situations. In fact, primacy coding can be interpreted as normalized coding with an inhibition strength α that depends on the non-dimensional width of the concentration distribution; see Methods and Models. Primacy coding is thus an example for global inhibition with instantaneous adaptation, which displays better performance than a simple fixed threshold γ. Taken together, this simple model comparison indicates that the mean response of the olfactory systems needs to be controlled and that simple normalization is not sufficient for this. So far, we calculated the transmitted information and tested the discrimination performance of primacy coding under the assumption that all receptor types behave similarly. In fact, we established that the maximal information is achieved when all receptor types are activated with equal probability NC/NR. However, neither the receptor sensitivities nor the odors themselves are distributed equally in realistic situations. Variations in these quantities affect the transmitted information and thus the usefulness of the primacy code. For instance, the transmitted information decreases if a single receptor is activated less often than all the others; see Fig 7A. This effect is small, since in the worst case the receptor is never active and the transmitted information thus corresponds to an array with this receptor removed. Conversely, having a receptor that is active more often than all others can have a much more severe effect; see Fig 7A. In fact, if the receptor type is more than three times as active, the transmitted information I is lower than if the receptor type was remove completely; see Methods and Models. This indicates that receptors can shadow the response of other receptors and thus be detrimental to the overall array when they are overly sensitive. The effect of varying receptor sensitivities can be studied in our model of primacy coding by discussing more general sensitivities matrices. We consider S n i = ξ n S n i iid, where each receptor type can have a different sensitivity factor ξn, which modulates the uniform sensitivity matrix S n i iid where each entry is independently chosen from the same log-normal distribution. The case of homogeneous sensitivities that we discusses so far thus corresponds to ξn = 1 for n = 1,2, …, NR. To investigate the effect of heterogeneous sensitivities, we start by varying the sensitivity factor of one receptor type while keeping all others untouched, i. e. , we change ξ1 while keeping ξn = 1 for n ≥ 2. There are three simple limits that we can discuss immediately. For ξ1 = 0, the first receptor type will never become active, the array behaves as if this type was not present, and the transmitted information is approximately Imax (NC, NR − 1). This value is lower than the maximally transmitted information Imax (NC, NR) reached for the symmetric case ξ1 = 1. However, the associated information loss ΔI = Imax (NC, NR) − Imax (NC, NR − 1) ≈ NC/ (NR ln 2) is relatively small in large receptor arrays (NR ≫ NC); see Fig 7B. Conversely, the transmitted information can be affected much more severely if the sensitivity of the first receptor type is increased beyond ξ1 = 1 and the receptors will thus be active more often than the others. In the extreme case of ξ1 → ∞, the first receptor type will always be active and thus not contribute any information. Since this receptor type would always be part of the primacy set, the information transmitted by the remaining receptor types is approximately Imax (NC − 1, NR − 1), which is smaller than Imax (NC, NR − 1) in the typical case NR ≫ NC. Consequently, an overly active receptor type can be worse than not having this type at all under primacy coding. The fact that overly sensitive receptors are detrimental to the transmitted information is also visible in numerical simulations. Fig 7B shows ensemble averages of the information I transmitted by receptor arrays as a function of the sensitivity factor ξ1. As qualitatively argued above, I is maximal for ξ1 = 1 and it is slightly lower for smaller ξ1 since the receptor type is active less often. In contrast, for ξ1 > 1, I decreases dramatically and falls below the value of ξ1 = 0 for ξ1 ≳ 1. 5. These data suggest that it would be better to remove receptor types that exhibit a 50% higher sensitivity than the other types. To see whether overly sensitive receptor types are also detrimental when all types have varying sensitivities, we next considering sensitivity factors ξn distributed according to a lognormal distribution. Numerical results shown in Fig 7C indicate that the transmitted information indeed decreases with increasing variance var (ξn) of the sensitivity factors. In fact, a variation of var (ξn) /〈ξn〉2 = 0. 5 already implies a reduction of the transmitted information by almost 50% for small concentration variations σ/μ = 1. If the odor concentrations vary more, the information degradation is less severe, but the same trend is visible. Interestingly, rescaling the information by the maximal information Imax given in Eq (5) collapses the curves for all dimensions NC and NR, suggesting that this analysis also holds for realistic receptor repertoire sizes. Note that the reduced transmitted information also implies poorer odor discrimination performance; see Fig 7D. Taken together, this provides a strong selective pressure to limit the variability of the receptor sensitivities so overly sensitive receptors do not dominate the whole array.
  We analyzed a simple model of neural representations of olfactory stimuli, where odors are identified by the NC strongest responding receptor types. This version of primacy coding provides a sparse representation of the odor identity, which is independent of the odor intensity. We showed using numerical simulations and a statistical model that the primacy dimension NC strongly affects the transmitted information and the discriminability of odors. Interestingly, already for small values of NC ≲ 10, the typical olfactory discrimination tasks can be carried out with performances close to experimentally measured ones. Conversely, the number NR of receptor types does not strongly affect the coding capacity and the discriminability of similar odors, in accordance with lesion experiments. Our model even indicates that lowering NR can improve the identification of a target ligand in a background. The advantage of our simple model is that we can analyze its behavior in depth and explicitly link the statistical properties of the olfactory system to data from psycho-physical experiments. In particular, we predict how likely two different odors drawn from a particular statistics can be distinguished. For instance, our model implies that target odors are easier to detect if disturbing backgrounds consist of many ligands. We generally find that representations are sensitive to the relative concentration of ligands in mixtures and that dilute components are basically completely shadowed. Conversely, for fixed ligand concentrations mixtures can typically be discriminated very well. However, identifying the individual ligands in mixture is only possible for mixtures with few components. In any case, our results suggest that the primacy code formed in the olfactory bulb is more useful to identify odors in the subsequent olfactory cortex than simple alternatives, essentially because the statistics of the representations are independent of the odor statistics. Our model predicts that receptors are only useful if their likelihood to respond to incoming odors is similar. This is because receptor types that are overly sensitive and respond strongly to many odors could dominate other types and thus degrade the total information. In fact, having a receptor type that is 50% more sensitive than others, and thus responds about three times as often, can lead to less transmitted information than when this type is absent. This observation is related to the primacy hull discussed in [38], which also predicts strong restrictions on the receptor sensitivities stemming from primacy coding. Various strategies could play a role in keeping the activity of the receptor types similar [53]: On timescales as short as a single sniff, the inhibition strength could be adjusted to regulate the relative importance of receptor excitations [54]. On longer timescales of several weeks, there are changes of the receptor copy number that directly affect the sensitivity of the glomeruli [55–57] and the processing neurons in the olfactory bulb [58,59]. Receptor copy number adaptations influence the signal-to-noise ratio at the receptor level, so the copy number could be increased to improve the detection of frequently appearing odors [60]. In contrast, we predict a decrease of the copy number of overly sensitive receptor types that respond often. Combining the two alternatives, receptor copy numbers could be controlled such that noise is suppressed sufficiently while ensuring that single receptor types do not dominate the array. Finally, receptor sensitivities can also be adjusted by genetic modifications on evolutionary timescales [61,62]. Moreover, direct feedback from higher regions of the brain could modify the processing of olfactory signals, e. g. , in response to the behavioral state [7]. Although our work shows that the activities of the receptors need to be balanced, the actual distribution of the sensitivities matters much less. For instance, log-uniform distributions, which have been suggested to describe realistic receptor arrays [51,63], lead to similar odor discriminability as log-normally distributed sensitivities; see S2 Fig. Our results raise the question why mice have 20 times as many receptor types than flies although the transmitted information under primacy coding is only increased by a factor of 2 (see Eq (5) ) and the odor discriminability is hardly affected by the receptor repertoire size (see Fig 4). The apparent usefulness of large receptor repertoires hints at roles of the olfactory system beyond transmitting the maximal information and discriminating average odors. For instance, having many receptor types might help to hardwire innate olfactory behavior when receptors are narrowly tuned to odors. In this case, our model would only apply to the fraction of the receptor types that are broadly tuned and are not connected to innate behavior. Alternatively, having many receptor types might be advantageous to discriminate very similar odor mixtures, to cover a larger dynamic range in concentrations of individual ligands, or to allow for a larger variation in average sensitivities, enabling quick adaptation to new environments. Finally, biophysical constraints of the receptor structure might imply that many receptors are required to cover a large part of chemical space. We discussed the simplest version of primacy coding with a minimal receptor model and a constant primacy dimension NC implemented by a hard threshold. This model neglects the complex interactions of ligands at the olfactory receptors, which can affect perception [64]. In particular, antagonistic effects can already provide some normalization at the level of receptors [65]. Generally, it is likely that many mechanisms contribute to the overall normalization of the receptor response [66]. A more realistic model of primacy coding might also consider a softer threshold, where receptor types with larger excitation are given higher weight in the downstream interpretation, which is related to rank coding [22]. In this case, information from fewer glomeruli might be sufficient to identify odors, since the rank carries additional information. Realistic olfactory systems could also use a timing code, taking into account more and more receptor types (with decreasing excitation) until an odor is identified confidently. Such a system could explain that the response dynamics in experiments depend on the task [67,68]. Generally, a better understanding of the temporal structure of the olfactory code [8,69–73] might allow to derive more detailed models. These could rely on attractor dynamics that are guided by the excitations and thus respond stronger to the early and large excitations [74,75].
  All numerical simulations are based on ensemble averages over sensitivity matrices Sni. The elements of Sni are drawn independently from a log-normal distribution with var (S n i) / S ¯ 2 = 1. 72 corresponding to λ = 1. In Figs 2A and 7B–7D, an additional ensemble average over odors c is performed using the distribution Penv (c). Here, odors c are chosen by first determining which of the NL ligands are present using a Bernoulli distribution with probability p = s/NL and then independently drawing their concentration from a log-normal distribution with mean μ and standard deviation σ. In all simulations the primacy set a corresponding to c is given by the NC receptors with the highest excitation calculated from Eq (1). Statistics of a and the transmitted information I given by Eq (4) are determined by repeating this procedure 105 and 107 times, respectively. In order to obtain deeper insights into the numerical results, we also develop analytical approximations using a statistical description of all involved quantities, which is based on accounting for the means and variances of the respective distributions. For instance, the statistics of the output a given by Eqs (1) – (3) can be estimated using ensemble averages of sensitivity matrices for different odors c, similar to our treatment presented in [21] and [30]. In particular, Eq (1) implies that the effects of different ligands are additive. Since the log-normal distribution describing the sensitivities is narrow (λ = 1), the excitations en are also well approximated by a log-normal distribution with mean 〈 e n 〉 S = S ¯ ∑ i c i and variance and var S (e n) = var (S n i) ∑ i c i 2 [76], whereas correlations are negligible [21]. The probability that the excitation en exceeds the threshold γ, and the associated receptor type is thus part of the primacy set, reads ⟨ a n ⟩ S = 1 - G (γ (c) ⟨ e n ⟩ S; ζ (c) ) (7) with G (x; ζ) = 1 2 + 1 2 erf (ζ + log (x) 2 ζ 1 2) (8) being the cumulative density function of a log-normal distribution with 〈x〉 = 1 and var (x) = exp (2ζ) − 1. The width of the distribution is determined by the positive parameter ζ = 1 2 ln (1 + var (e n) / 〈 e n 〉 2), which reads ζ (c) = 1 2 ln [ 1 + (e λ 2 - 1) ∑ i c i 2 (∑ i c i) 2 ] (9) for an ensemble average over sensitivities. Note that ζ is concentration-invariant, since it does not change when the concentration vector c is multiplied by a constant factor. In the simple case of ligands that are distributed according to Penv (c), we find 〈 (∑ i c i 2) (∑ i c i) - 2 〉 c = s - 1 (1 + σ 2 / μ 2). Consequently, the distribution width ζ is large for broadly distributed sensitivities (large λ), few ligands in an odor (small s), and wide concentration distributions (large σ/μ). The constraint Eq (3) implies 〈an〉 = NC/NR, so that the mean threshold reads ⟨ γ ⟩ = ⟨ e n ⟩ S · G - 1 (1 - N C N R; ζ), (10) where G−1 is the inverse function of G defined in Eq (8). Using this expression as an estimate for γ in Eq (7) results in concentration-invariant activities an, since 〈γ〉 is proportional to the excitation 〈en〉. This situation is comparable to simple normalized representations resulting from the threshold γ = α〈en〉, where α is a constant inhibition strength [21]. In fact, primacy coding can be interpreted as global inhibition with an inhibition threshold depending on the width of the excitation distribution, α = G - 1 (1 - N C N R - 1; ζ).
summary: Humans can identify odors independent of their intensity. Experimental data suggest that this is accomplished by representing the odor identity by the earliest responding receptor types. Using theoretical modeling, we here show that such a primacy code outperforms alternative encodings and allows discriminating odors with close to experimentally measured performance. This performance depends strongly on the number of receptors considered in the primacy code, but the receptor repertoire size is less important. The model also suggests a strong evolutionary pressure on the receptor sensitivities, which could explain observed receptor copy number adaptations. By predicting psycho-physical experiments, the model will thus contribute to our understanding of the olfactory system.
keywords: medicine and health sciences
  brain
  social sciences
  neuroscience
  probability distribution
  mathematics
  computational neuroscience
  coding mechanisms
  sensory physiology
  animal cells
  sensory receptors
  olfactory receptor neurons
  probability theory
  signal transduction
  piriform cortex
  cellular neuroscience
  psychology
  olfactory bulb
  cell biology
  anatomy
  physiology
  neurons
  olfactory system
  biology and life sciences
  cellular types
  afferent neurons
  sensory perception
  sensory systems
  physical sciences
  computational biology
title: Primacy coding facilitates effective odor discrimination when receptor sensitivities are tuned