Skip Navigation


Chemical Senses Advance Access originally published online on August 16, 2006
Chemical Senses 2006 31(9):795-806; doi:10.1093/chemse/bjl021
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
31/9/795    most recent
bjl021v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (3)
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Johnson, B. N.
Right arrow Articles by Sobel, N.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Johnson, B. N.
Right arrow Articles by Sobel, N.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2006. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oxfordjournals.org

A Comparison of Methods for Sniff Measurement Concurrent with Olfactory Tasks in Humans

Bradley N. Johnson1,2, Cheryl Russell3, Rehan M. Khan2 and Noam Sobel1,2,4

1 Joint Graduate Group in Bioengineering, UC Berkeley & UC San Francisco 2 Helen Wills Neuroscience Institute, 3 Department of Bioengineering, UC Berkeley 4 Department of Psychology, UC Berkeley

Correspondence to be sent to: Bradley N. Johnson, Helen Wills Neuroscience Institute, UC Berkeley, Berkeley, CA 94720, USA e-mail: bnjohnso{at}socrates.berkeley.edu


    Abstract
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Appendix
 Acknowledgements
 References
 
There is a growing appreciation for the role of sniffing in the formation of the olfactory percept. With this in mind, monitoring and measurement of sniffing is an important aspect of olfactory experiments. There are several methods for measuring human sniffs concurrent with odor delivery in olfactory experiments. Here, we set out to compare the temporal sensitivity and power of these different methods by applying them all simultaneously with an olfactory task. We discuss the advantages and disadvantages of each method and conclude in recommending the use of a nasal cannula linked to a pressure sensor whenever possible.

Key words: psychophysics, spirometer, respitrace, temperature, plethysmograph


    Introduction
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Appendix
 Acknowledgements
 References
 
Researchers have long debated the proper management of sniffing during olfactory experiments. Whereas several early studies argued for the elimination of sniffing within olfactory experiments (Gamble 1899Go; Elsberg and Levy 1935Go), others have highlighted the importance of allowing near-natural olfactory behavior during olfactory tasks (Mainland and Sobel 2006Go). Indeed, the important role of sniffing in the formation of the olfactory percept has been highlighted in a recent series of reviews in this journal (Buonviso et al. 2006Go; Kepecs et al. 2006Go; Mainland and Sobel 2006Go; Schoenfeld and Cleland 2006Go). Considering the influence of sniffing on the resultant olfactory percept (Le Magnen 1945Go; Laing 1983Go) and on patterns of neural activity throughout the olfactory system (Tucker 1963Go; Sobel et al. 1998Go; Laurent 2002Go; Spors and Grinvald 2002Go; Kepecs et al. 2006Go), accurate measurement of sniff parameters is paramount. Although one can point to multiple sniff parameters of interest (Youngentob et al. 1987Go), the most important sniff measures are related to the temporal dynamics of the sniff and the resultant airflow volumetrics that deliver the odorant to the olfactory epithelium.

One can posit 2 reasons for measuring airflow in an olfaction study. The first reason is to know exactly how much odorant-laden air arrived at the sensory surface, namely, the olfactory epithelium. Experiments in cadavers (Paulsen 1882Go) and computational fluid dynamic models using computed tomography or magnetic resonance imaging generated from human nasal structures (Hahn et al. 1994Go; Keyhani et al. 1995Go; Zhao et al. 2004Go, 2006Go) have shown that the amount of air reaching the sensory epithelium is approximately 5–10% of the air entering the nostrils. However, considering individual variation in nasal structure, possibly enhanced by pathology, one cannot simply safely assume that 5–10% of what enters a nostril always reaches the epithelium. Due to this disparity, one would ideally measure the flow adjacent to the sensory epithelium. However, no such device for measuring airflow in the upper airway concurrent with an olfaction task is currently available for human studies.

A second reason for measuring sniffing during an olfaction task is to assure that sniffing behavior was held constant across different olfactory conditions or across repeated trials of the same condition. For example, one may be measuring the neural response to 2 different odors. In order to assure that variance in neural activity is reflecting the variance in odorants and not variance in sniffing behavior, one must assure that sniffing was constant. In this case, individual variability leading to differing epithelial exposure (e.g., polyps, septal deviation) across subjects is irrelevant because it remains constant within a subject across trials and conditions. It is this second reason for measuring sniffing that is addressed in this article.

Various methods have been employed in order to measure the external airflow sniff dynamics in humans. These include 1) respiratory belts placed around the abdomen, chest, or both that measure respiratory sniff airflow–related expansion and contraction, 2) thermistors placed at the nares that measure the sniff airflow–related changes in temperature, 3) cannulas places at the nares and linked to pressure sensors that measure sniff airflow–related changes in relative pressure, and 4) pneumotachometers placed at the inlet of a nasal mask where they similarly measure sniff airflow–related changes in relative pressure. These methods are not all equally effective in measuring sniffing within an olfactory task. At times, the only information necessary is whether the subjects sniffed and when. For such gross information, it is likely that any of these measures is sufficient. However, some experiments call for more accurate measurements, and it is unknown which of these techniques effectively captures sniff dynamics.

For example, sniffs are odorant specific. Whereas intense and unpleasant odorants are sampled with moderate sniffs, mild and pleasant odorants are sampled with vigorous sniffs (Laing 1983Go; Warren et al. 1994Go; Sobel et al. 2001Go; Walker et al. 2001Go). This inverse relationship between sniff vigor and odorant magnitude is evident within as little as ~200 ms following sniff onset (Johnson et al. 2003Go) and is sufficiently stereotyped so as to serve as an olfactory clinical diagnostic tool (Frank et al. 2003Go). When studying such odorant-dependent sniffing, one would want increased temporal sensitivity in sniff measurement combined with minimal interference of the olfactory task. Here we set out to compare the sensitivity of different methods for sniff measurement in the context of odorant-specific sniffing. Subjects sniffed either an intense unpleasant odorant or a mild pleasant odorant concurrent with sniff measurement using all the 4 above-mentioned methods simultaneously.


    Materials and methods
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Appendix
 Acknowledgements
 References
 
Subjects

Twenty-one college students (8 females), aged 18–23, participated in the experiment after giving informed consent to procedures that conformed to National Institutes of Health guidelines and were approved by the University of California, Berkeley, Committee on Protection of Human Subjects. Subjects were screened for no history of nasal trauma or neurological disease, no regular use of medication, and self-reported intact olfactory acuity.

Testing room

Testing was performed in a designated 9 x 11–ft room that was completely coated with stainless steel to prevent odor contamination. The room was equipped with high efficiency particulate air and carbon high-rate filtration to further minimize odorant contamination. Subjects were comfortably seated in a dental chair centered in the room. During experimentation, subjects were alone in the room, monitored via a one-way mirror window, and video monitor in the adjacent control room.

Sniff measurement

Subjects were fitted with plethysmograph belts (respitrace), nasal temperature sensors, nasal pressure sensors, and nasal mask pneumotachometer, allowing simultaneous recording of all measures (Figure 1).


Figure 1
View larger version (51K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 1 Subject setup. Panel A shows the cannula (CAN) that is attached to the piezoelectric pressure sensor (not shown) and the 2 thermistors (TMP). Small zip-ties are used to fasten the thermistors to each orifice of the cannula. Panel B shows the CAN and TMP apparatus attached to a subject. The remaining measurement devices are shown in panel C. The thoracic belt (THX) is secured below the arms, whereas the abdomen belt (ABD) is secure around the midsection at the level of the naval. The pneumotachometer flow head (TAC) is attached to the mask via a short flexible tube.

 
Respitrace
We used 2 piezoelectric respiratory belt transducers (MLT1132, ADInstruments, Grand Junction, CO) for the respiratory plethysmograph recording. One belt was placed at the level of the abdomen, and the other was placed at the level of the chest. Each belt contains a piezoelectric element that responds to length changes, requires no driving circuitry, and connects directly to an instrumentation amplifier.

Temperature
We used 2 thermistors that were designed specifically for the measurement of nasal airflow–related changes in temperature (MLT415/A, ADInstruments). One thermistor was placed at the opening of each nostril. Each thermistor was attached to a signal conditioner allowing resolution of 0.01 °C (Thermistor POD MC309, ADInstruments) and then to an instrumentation amplifier. Critically, when placing these sensors, we made sure that the sensing element was not in direct contact with the skin but rather suspended in the nostril orifice. Direct contact with the skin acts as a heat sink that significantly reduces sensitivity.

Pneumotachometer
We used a pneumotachometer (high-sensitivity flowmeter model #4719, Hans Rudolph, Inc., Kansas City, MO) that was attached in-line with the vent port of the nasal mask. The pneumotachometer differential pressure was measured and converted to a voltage signal using a spirometer (ML141, ADInstruments) that delivered the voltage to an instrumentation amplifier.

Pressure
We used a respiratory nasal cannula (#1600, Salter Labs, Avin, CA, http://www.salterlabs.com) attached to a small piezoelectric silicon device (30INCH-D-MV, AllSensors, http://www.allsensors.com) with adjacent signal conditioning circuitry (described in detail in Appendix) for the millivolt output from the transducer to an instrumentation amplifier. The pressure transducer was calibrated with a manometer to allow conversion from volts to millimeters of H2O.

All measurements were simultaneously acquired and sampled at 1 kHz using an ADInstruments PowerLab 16SP instrumentation amplifier and Chart v5.0 software for digital recording. All channels were low-pass filtered (100 Hz) by the instrumentation amplifier.

Olfactometer

The olfactometer was based on previously described theoretical designs (Benignus and Prah 1980Go; Kobal 1985Go; Prah et al. 1995Go). All wetted parts of the olfactometer were made of stainless steel or Teflon. The olfactometer was driven by custom software written in LabVIEW (National Instruments, Austin, TX) that controlled 1) odorant identity, 2) odorant concentration, 3) gas flow rate, 4) flow temperature, 5) flow humidity, and 6) odorant delivery timing. Odorant identity was achieved by flowing a carrier gas, medical-grade breathing air, over undiluted liquid of the odorant placed in a specially designed stainless steel canister with 510-mm2 surface area. The carrier gas pulled away the vapor of the evaporating liquid. Concentration range was accomplished by mixing clean and odorized air using electronic mass flow controllers (MFCs, M100B, MKS Instruments, Methuen, MA). At the olfactometer output, MFCs controlled a clean air line, an odorized air line, and 2 vacuum lines, all joining at a "railroad switch" under the subject's nose. Switching between odorant and no-odorant conditions was accomplished by activating a 3-way solenoid that pulled vacuum from only one of the 2 vacuum lines. This resulted in removal of either the odorized or the clean air from the railroad switch, allowing only the other air source to reach the nose. The outcome was a seamless shift between odorant and no-odorant conditions with no visual, auditory, tactile, humidity, or thermal cues as to the alteration. For more details on the construction and performance of the device used, see Johnson et al. (2003)Go.

Odorants

Two odorants were used in the study, phenethyl alcohol (99%+, P1 360-6, Sigma-Aldrich, http://sigmaaldrich.com) and valeric acid (99%+, 129542, Sigma-Aldrich, http://sigmaaldrich.com). These odorants were chosen because they are commonly reported as pleasant (phenethyl alcohol) and unpleasant (valeric acid) (Doty 1975Go). We used 1 concentration of each odorant that corresponded to 100% of full maximum at 6 LPM over 8 ml of liquid odorant in the stainless steel canister. In this article, the phenethyl alcohol trials are identified PEA and the valeric acid trials as VAL. In addition to the PEA and VAL conditions, we used a third condition of 6 LPM clean air and identified it as AIR.

In order to estimate the absolute concentration of the stimuli used, we measured the depletion rate of liquid odorant as a function of carrier gas flow rate (Allison and Katz 1919Go; Johnston 1967Go). Odorant was loaded into the olfactometer after measuring the initial mass with an analytical balance (L-Series, 0.1 mg precision, Acculab, Bradford, MA). Odorized air at 6 LPM carrier gas flow rate was produced by the olfactometer for a 20- to 30-min period, after which the remaining mass was measured. The total volume of air was measured by integrating the flow measurement by the olfactometer. The mass loss was converted to moles (mass/molecular weight) and then to molar concentration (moles of odor vapor per million moles of carrier gas at 1 atm and 37 °C). The obtained values were then used to assess the diluted concentration for the flow rate used in this experiment. Based on these measurements, we estimated the final sniffed concentrations to be about 0.15 ppm for PEA and 0.23 ppm for VAL.

Experimental design

Each experiment contained 30 trials with a 45-s intertrial interval, 10 trials for each condition (VAL, PEA, and AIR). For each trial, the following digitized audio and visual command queue was given: "please prepare to sniff at the tone," "3," "2," "1," tone. At the tone, the subject took one sniff that contained PEA, VAL, or AIR delivered in a counter-balanced order. Following the sniff, the subject used a trackball to rate the odorant intensity and pleasantness on an on-screen visual analog scale. All experimental components, including task instructions, odorant delivery and triggering, and acquisition of both respiratory and perceptual data, were computer controlled and administered via digitized voice, thus preventing any unintentional experimenter-related cues or bias.

Data reduction

Odorant perceptual data
Perceptual responses were averaged across repeats of the same odorants within each subject, graphed, and tested for significant differences.

Sniff data
We first eliminated one thermistor trace and added one computed respitrace measure as follows: There is often a difference across nostrils in the rate of airflow (Principato and Ozenberger 1970Go; Hasegawa and Kern 1977Go). Because the thermistor reflects airflow-dependent changes in temperature, it provides increased signal from the nostril with greater airflow (indeed, sometimes one nostril had almost no sniff-related changes in temperature). Thus, we first selected for each subject the thermistor trace with greater signal and used only that for further analysis. In turn, we added the commonly used computed respitrace measure that is the weighted sum of 2 times thoracic and abdominal measures (Banzett et al. 1995Go).

Sniffing data was first converted from Chart binary files into ASCII text files that were then loaded into MATLAB 6.5 (Mathworks, Inc., Natick, MA). MATLAB was used for all subsequent analysis. Each data trace was then preprocessed in an identical way—data was low-pass filtered, third-order Butterworth Fc = 10 Hz, and then phase corrected to compensate for the temporal shift produced by the filter. All trials were then aligned in time, DC offset was removed, and then data was scaled to the maximum (minimum for temperature) value collected for that measure in that subject. Finally, for each trial, the data trace was truncated at 4000 ms.

To parameterize the time series data for statistical analysis, we calculated 3 summary variables: mean value, max value, and the integral (e.g., area under the curve). For each of these summary variables of each measurement technique, we calculated an F-statistic using a 1-way repeated measures analysis of variance (ANOVA). Additionally, we calculated the F-statistic versus time. We did this in 2 slightly different ways. In the first method, the input parameter to the ANOVA was the value from the time series at that particular time point. In essence, we created an "instantaneous" estimate of the F-statistic that was independent from prior ANOVAs in as much as the time series data are independent. In contrast, the second method calculated the F-statistic for successively longer time periods, where we used the integral within this growing time period as the measured parameter for the ANOVA. That is, each successive ANOVA computation was the same as the previous one, except that we added another 10 ms of data. This procedure created a "running" estimate of the F-statistic versus time that was dependent on all data back to time zero.

To compare the overall statistical sensitivity of each measure, we used bootstrapping (Efron and Tibshirani 1993Go) to create distributions of the F-statistic for the 3 summary variables (mean, max, and integral) for each measurement technique and then used these distributions to compare the various measurement techniques. The bootstrapping was done as follows. Within each subject and condition, we randomly selected 10 repeats (e.g., trials), allowing individual repeats to be selected more than once (e.g., random with replace). These 10 repeats were averaged and thus formed a "resampled" version of a particular subject's response to condition X. After resampling all subjects and the 3 conditions, we computed the 3 summary variables (mean, max, and integral) and then computed the F-statistic from these values. This technique allowed us to estimate the variance of the F-statistics determined from the original data.


    Results
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Appendix
 Acknowledgements
 References
 
Data from 7 subjects was eliminated from the subsequent data analysis. One subject was excused because of technical problems (communication between the olfactometer and acquisition software) that could not be resolved during the allotted experiment time period, one subject was excused during the training period because they were uncooperative, three subjects were eliminated for missing more than 3 trials during the experiment, and two subjects were eliminated for failure of the nasal cannula measurement channel. The 2 nasal cannula failures reflected the various hardware (wire, cabling, tubing, etc) of one measurement interfering with the hardware of another measure. Typically, movement of the nasal mask or movement of the thermistor cable pulled out the nasal cannula. Because the simultaneous recording of all measures was critical for this study, subjects where any measure was eliminated were then completely eliminated from analysis. This left 16 subjects for further analysis. It is noteworthy that this attrition rate does not reflect a realistic attrition rate of the nasal cannula but rather is related to the specific design of this study where placement of one measure complicated the placement of another. In a typical olfaction study, one would use only one of these devices, and from our experience, when used alone, the attrition rate of these devices is negligible.

We conducted 4 separate analyses of our data. The first analysis was to determine the perceptual differences between the 3 stimulus conditions. In the second analysis, we assessed each techniques' ability to reveal odorant-specific sniffing. In the third analysis, we directly compared the signal-to-noise ratios for each measurement technique. Finally, in the fourth analysis, we examined the temporal response of each measurement technique.

Odorant perceptual data

We first examined the intensity and pleasantness responses to ensure that the chosen odorants were perceived differently. The group mean perceptual results shown in Figure 2 demonstrate that we successfully generated 2 very different odorants. There was a significant effect of condition on perceived intensity (F(2, 30) = 38.9, P < 4.6 x10–9) and perceived pleasantness (F(2, 30) = 22.27, P < 1.2 x10–6). Post hoc analysis using multiple comparisons–corrected t-tests showed a significant intensity difference between AIR and VAL (t(15) = –7.9371, P < 1 x10–6), PEA and AIR (t(15) = 2.3245, P < 0.035), and PEA and VAL (t(15) = –5.6064, P < 1 x10–4) and a significant pleasantness difference for AIR and VAL (t(15) = 5.2021, P < 1 x10–4) and PEA and VAL (t(15) = 4.6485, P < 3 x10–4) but not for PEA and AIR (t(15) = –1.6439, P = 0.121). In summary, VAL differed from PEA in both intensity and pleasantness. See Table 1 for mean and variance estimates for the ratings of each condition.


Figure 2
View larger version (14K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 2 Perceptual results. Group mean estimates for intensity (left) and pleasantness (right) for the 3 conditions: valeric acid (VAL), clean air (AIR), and phenethyl alcohol (PEA). Errors bars are ±1 SE but are not representative of actual within-subject differences. There was a significant effect for condition on perceived intensity (F(2, 30) = 38.9, P < 4.6 x 10–9) and perceived pleasantness (F(2, 30) = 22.27, P < 1.2 x 10–6).

 

View this table:
[in this window]
[in a new window]

 
Table 1 Consolidated data from ANOVA tests

 
It is important here to remind that if one were setting out to test which particular olfactory aspects most influence sniffing, the above olfactory perceptual space would be suboptimal in that it would be hard to dissociate whether odorant intensity or pleasantness influenced sniffing patterns. In this study, however, we set out to create 2 stimuli that differed in a manner that was sure to be accompanied by different sniffs, so that we could address our question of interest, namely, the sensitivity of different sniff measurement techniques.

Summary variables revealing odorant-induced changes in sniffing

To begin analyzing the sniffing data, we determined which technique showed differences between our 3 odorant conditions. We created the various summary variables (mean, max, integral) as detailed in Data Reduction and then calculated a 1-way repeated measures ANOVA to determine which summary variables were different between the 3 odorant conditions. The group mean time series and group mean summary variables for each measurement technique are shown in Figure 3. Overall, we found that every measure but one was able to discriminate odorant-induced sniffing changes to some degree. Only the thorax belt (and combo thorax + abdominal) showed no difference for any of the 3 summary variables (volume, mean, or max). However, the abdomen belt showed significant differences for mean (F(2, 30) = 3.534, P < 0.042). The pneumotachograph showed significant differences for volume (F(2, 30) = 4.93, P < 0.014) and mean (F(2, 30) = 4.06, P < 0.023). The nasal cannula showed significant differences for volume (F(2, 30) = 8.18, P < 0.002), max (F(2, 30) = 4.38, P< 0.02), and mean (F(2, 30) = 7.36, P < 0.003). The temperature sensor showed significant differences for volume (F(2, 30) = 7.31, P < 0.003), max (F(2, 30) = 12.11, P < 0.0001), and mean (F(2, 30) = 8.88, P < 0.0009). Post hoc analysis to determine individual condition differences consistently showed separation between VAL and AIR and less frequently between AIR and PEA or VAL and PEA. See Table 1 for mean and variance estimates for each summary variable.


Figure 3
View larger version (42K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 3 Group mean time series (left most panel) and group mean summary variables (right 3 panels) for each measurement technique. Solid lines and bars are for valeric acid (VAL), dotted lines and bars for clean air (AIR), and dashed lines and bars for phenethyl alcohol (PEA). The measurements are abbreviated: thorax belt (THX), abdomen belt (ABD), respitrace combination (CMB), pneumotachograph (TAC), nasal cannula (CAN), and temperature (TMP). The errors bars in the summary variable panels represent ±1 SE and represent the between-subject error and not the within-subject as tested by the ANOVA. Panels with a significant ANOVA (P < 0.05) are highlighted by the gray background shading.

 
Rank ordering of measurement signal-to-noise

We next set out to determine which measurement technique showed the greatest differences between the 3 odorant conditions. We used the F-statistic from the ANOVA as an estimate of each techniques' signal-to-noise ratio, that is, how sensitive it was to differences among the different sniff conditions.

The F-statistic for each summary variable is shown in Figure 4 and are the same F-statistics reported above (summary variables revealing odorant-induced changes in sniffing). The dotted line in Figure 4, at F = 3.3158, represents the critical value for an alpha = 0.05 level of significance. The bar graphs convey the picture that the spirometer, pressure, and temperature have a generally better signal-to-noise ratio than the respiratory band methods. In order to compare the measures statistically, we need to know the error distribution of the F-statistic in each condition. We used bootstrapping, a resampling procedure (Efron and Tibshirani 1993Go), to construct error distributions for the F-statistic in each condition.


Figure 4
View larger version (16K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 4 F values of the 1-way repeated measures ANOVAs on the actual data, for each summary parameter: integral (left), max value (middle), and mean value (right). The 6 measurement techniques are represented within each bar graph as follows: thoracic belt (THX), abdominal belt (ABD), combined belts (CMB), pneumotachograph (TAC), cannula (CAN), and temperature (TMP). The dotted line (F = 3.3158) is alpha = 0.05 (degrees of freedom, df = 2,30). The error bars represent the standard deviation of the bootstrapped F-statistic distribution and allow comparison between individual bars to assess if one measure is significantly better than the second (e.g., temperature max vs. cannula max) and not if the measures themselves are significant (as both temperature max and cannula max are significant).

 
For every bootstrap iteration, we obtained an estimated F-statistic for each of the 18 variables shown in Figure 4. These estimates allowed us to generate distributions for each of these 18 variables. We then calculated the standard deviation for each distribution and displayed them as the error bars in Figure 4. The means of these distributions generally agreed with, although were slightly less than, the F-statistics calculated from the actual data (e.g., that shown in Figure 4). These variance estimates allow for visual comparison of the 18 variables.

To statistically compare the 18 variables, we again used the bootstrap analysis. To determine if these differences in the F-statistic are significant, we pairwise compared the measurement techniques that were significant, abdomen belt (ABD), pneumotachograph (TAC), nasal cannula (CAN), and temperature (TMP), for each summary variable (volume, max, and mean) for each bootstrap iteration. That is, for each of the 1,000 iterations of the bootstrapping, we obtained a number that was the difference in F-statistic value for a particular summary variable between 2 measurements, for example, pneumotachograph and nasal cannula (TAC vs. CAN). If there were no difference between the F values produced by the measurement techniques, then the pairwise subtraction of the F-statistic values would yield a distribution centered around zero. That is, for every bootstrap iteration, there is an equal probability that the pairwise subtraction will fall above or below zero. Thus, we estimated the probability to observe only N samples above (or below) zero as N/500 or alpha = N/500. For the volume summary variable, we found significantly higher signal-to-noise (F value) for the TAC versus ABD (P < 0.001), for the CAN versus ABD (P < 0.001), and for the TMP versus ABD (P < 0.004) but not for CAN versus TAC (P ~ 0.09), TMP versus TAC (P ~ 0.61), or TMP versus CAN (P ~ 0.58). For the max value summary variable, we found significant differences for TMP versus ABD (P < 0.02) and TMP versus TAC (P < 0.04) but not for TAC versus ABD (P ~ 0.77), TAC versus ABD (P ~ 0.77), CAN versus ABD (P ~ 0.27), CAN versus TAC (P ~ 0.20), or TMP versus CAN (P ~ 0.22). Using the mean parameter, there was no significant difference between the measurement parameters.

Comparing temporal response of each measurement

In our last analysis, we examined the behavior of the F-statistic as a function of time (e.g., from sniff onset to the end of the sniff) to compare the temporal responses of the measurement techniques. To examine the behavior of the F-statistic over time, we plotted the instantaneous and running F value for our actual data (not the resampled bootstrap iterations) in Figure 5. There is a very small time window where the abdomen belt is significant (~1000 ms, P < 0.05); however, this separation is not sustained as in the other measurements. The spirometer is significant past ~500 ms (instantaneous) or ~1000 ms (running). The pressure is significant past ~400 ms (instantaneous) or ~600 ms (running). Finally, the temperature sensor did not reach significance until ~900 ms (instantaneous) or ~1100 ms (running).


Figure 5
View larger version (14K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 5 P values, as determined by the F-statistic, as a function of time for (top to bottom) the abdomen belt, spirometer, pressure, and temperature sensors. Solid lines are for the 10 ms "instantaneous" estimate and dashed lines for the "running" estimate of the F-statistic. The abdomen belt showed only a very small window where the instantaneous estimate was significant. Conversely, the spirometer, pressure, and temperature sensors all showed sustained discrimination of the odorant conditions.

 
These findings are consistent with the previous analysis using volume, max, and mean, but the added temporal sensitivity has allowed us to locate the obvious temporal performance differences between the spirometer, pressure, and temperature sensors. Overall, the temperature sensor had the highest F-statistic, provided it had access to the entire sniff. However, when temperature was examined for the first second, not only did it perform worse than the pressure sensor but also the F-statistic for the actual data did not produce an alpha < 0.05 until after ~1 s, whereas both the pressure and spirometer did so at around 500–600 ms (Figure 5).


    Discussion
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Appendix
 Acknowledgements
 References
 
In this study, we sought to determine which of the most commonly used sniff-monitoring techniques are suitable for use in olfactory experiments. We were particularly interested in techniques that were sensitive enough to show odorant-related sniff modulation. We found that a thoracic respiratory belt did not provide information sufficient for dissociating odorant-specific sniffs. Furthermore, a commonly used computed combination of thoracic and abdominal respiratory belts was also not sufficiently sensitive. This information is important because often these very measures have been used to assure that sniffing was constant across different olfactory conditions within a study. Our results suggest that in the future other methods of measuring sniffing should be applied in order to substantiate such statements. In contrast, an abdominal respiratory belt, temperature sensors in the nostrils, and measures of sniff-induced pressure differences measured either in the nose by cannula or around the nose through nasal mask were all sufficiently sensitive. If one is interested in an overall odorant-dependent sniff difference, the rank ordering of these measurements from most to least powerful is as follows: temperature, nasal cannula (attached to pressure sensor), pneumotachometer flow head (attached to spirometer pressure sensor), and lastly abdominal belt.

In contrast, if one is interested in the temporal dynamics of the sniff, in other words, not only whether odorant-related sniff modulation occurred but also when it first emerged, then the rank ordering of these measurements from most to least time sensitive was as follows: nasal cannula (attached to pressure senor), pneumotachometer flow head (attached to spirometer pressure sensor), and lastly abdominal belt.

It addition to their rank ordering in power and temporal sensitivity, each of these measures enjoys particular advantages and disadvantages that may render them particularly appropriate or inappropriate for a given experiment. Overall, each of the techniques disappoint because they do not directly measure the flow of odorant-laden air that reaches the olfactory sensory surface. Ideally, the "gold-standard" for which we had based our judgments of these 5 external measurement techniques would be a simultaneous recording of the airflow in the upper airway. Because no such device is currently available, we instead judged the techniques on their ability to show odorant-related sniff modulation.

What we measured with the external devices was the rapid change in sniff vigor that was related to the content of the odorant-laden air that reached the olfactory epithelium. Hence, because we did not directly measure the upper airway flow, we did not directly measure the flow that induces the sniff modulation. This sniff modulation manifests as a change in inspiratory drive that is intended to change the amount of odorant-laden air that reaches the sensory surface. The modulation of inspiratory drive changes not just the airflow in the upper airway but also the total air entering the nostril. It is this change in the total air entering the nostril that we measured here. However, in some experiments, particularly those designed to directly study the dynamics of the modulation-inducing airflow, all the techniques used here may be inadequate.

Abdominal respiratory belt

The abdominal respiratory belt was the least powerful of the methods still capable of dissociating odorant-specific sniffing. Respiratory plethysmograph was originally designed as a nonintrusive measure of tidal volume. There are many reports dedicated to assessing this technique's accuracy (Konno and Mead 1967Go; Milledge and Stott 1977Go; Chadha et al. 1982Go; Gonzalez et al. 1984Go; Watson et al. 1988Go) and its validity toward diagnosis of respiratory function, especially when using the correlation and phase relationship (the so called Konno–Mead plot) between the abdominal and thoracic expansion (Tripp and Bolton 1998Go; Willis et al. 2004Go). That said, all these studies apply to data obtained during "quiet breathing." In contrast, olfactory sniffing constitutes "active breathing" where the normal phase relationship between thoracic and abdominal expansion breaks down. Here we observed that nearly every subject had a different phase relationship between the 2 belts, and even within a single subject, this relationship changed through the course of the experiment. These disadvantages of using respiratory belts to measure sniffing stem from what is also their major redeeming feature—they do not record from the nose. Thus, one advantage of this method is its applicability to experiments where the nose is engaged within some sort of apparatus that does not permit application of direct measurement devices. A final advantage of this method is its availability. For example, respiratory belts often come "built-in" to MRI scanners and are thus easily applied in fMRI of olfaction experiments. However, the lesson learned from our results is that the belt should be placed at the abdomen and not at the thorax during such experiments.

Temperature sensors

It was surprising to us that the temperature sensor performed so well. Thermistors are a popular technique to measure sniffing (Utenick 1971Go; Laing 1982Go, 1983Go, 1985Go, 1986Go; Uchida and Mainen 2003Go), but they are difficult to calibrate and have poor inherent temporal resolution (temporal resolution can be improved if they are heated above atmospheric temperature, but this can be inconvenient). Indeed, one possible explanation for the marginally higher statistical power of the temperature sensor is the inherent temporal filtering of the device. It is very likely that we could increase the statistical power of the cannula and pneumotachometer by applying a more aggressive low-pass filter. In our analysis, we used a third-order Butterworth low pass with Fc at 10 Hz to prevent altering the temporal dynamics of the measurements. Had we used Fc = 1 Hz or even Fc = 1/3 Hz, it is likely we would have much higher overall discrimination power at the cost of earlier temporal discrimination. On this basis, we repeated our analysis after applying a third-order Butterworth low pass with Fc = 1 Hz. We observed an improvement in the cannula and pneumotachometer F-statistic for only the max value summary parameter. We found that the more aggressive low-pass filter "rounded off" the top of each sniff peak, thus making the max value more consistent, yet maintained the overall shape and duration of the sniff, thus not affecting the integral or mean value summary variables.

Finally, despite the statistical power of the temperature sensors, 2 aspects reduce its convenience. First, application of thermistors in demanding environments such as the MRI is restricted. Second, because one nostril can have dramatically lower airflow rates than the other, it is critical to measure using 2 thermistors, one in each nostril, in order to assure sufficient signal.

Pneumotachometer

The pneumotachometer linked to spirometer provided both good statistical power and good temporal resolution. Furthermore, although we used a high-end device, precalibrated USB-based pneumotachometer-spirometer–based instruments are readily available, accurate, and moderately priced. An additional advantage of this method is that the sensor, namely the pneumotachometer, can be linked to the spirometer by pressure tubing alone, with no wires. Considering the availability of pneumotachometers with nylon rather than steel meshes, this allows easy application of the pneumotachometer in demanding environments such as an MRI. In turn, a drawback to pneumotachometers coupled to spirometers is that they require a nasal mask to capture and direct airflow through the device orifice, and breathing through a mask can itself alter breathing patterns (Gilbert et al. 1972Go; Perez and Tobin 1985Go; Rodenstein et al. 1985Go; Paek and McCool 1992Go). In olfaction experiments, this concern may be unavoidable when a nasal mask is a necessary part of the odorant delivery apparatus. However, if a nasal mask is not necessary for odorant delivery, adding such a mask just for the sake of spirometry is a disadvantage. Finally, an additional concern one may raise regarding the use of a pneumotachometer is that a subject's breathing pattern may be influenced by their awareness of monitoring (Han et al. 1997Go). However, in our experience, subjects always assume that the nasal mask is only there to deliver odors, and they practically never think that it is also measuring their respiration (Bensafi et al. 2003Go).

Nasal cannula

The pressure measurement via nasal cannula provided good statistical power and high temporal resolution. Like the pneumotachometer linked to spirometer, the cannula linked to pressure sensor measures differences in pressure resulting from the sniff. However, this device has several advantages over the pneumotachometer. Because it is placed directly in the nares, its temporal resolution is not dampened by the air volume of a nasal mask, as is the case for the pneumotachometer. This is what probably accounts for the slight advantage this method offers in temporal resolution. Furthermore, the geometry of the front end, namely the cannula, enables measurement with or without a nasal mask. Furthermore, like the pneumotachometer and in contrast to thermistors, this method does not entail leading wires into the nose. Use of pressure tubing alone provides maximal flexibility in environments such as the MRI. Last but not least, pressure measurement via nasal cannula is cheap. The cost of the current device was approximately US$120. Indeed, considering these advantages, it is not surprising that this type of method has been applied to devices aimed at measuring sniffing for diagnostic purposes (Frank et al. 2003Go). The only drawbacks we see to this method are that it is difficult to calibrate and it calls for some care when securing the cannula in place within the nares. Data from 2 subjects were lost in this study due to the cannula slipping out of the nares mid study, albeit these instances were instigated by the numerous other hardware weighing on the cannula tubing. These drawbacks not withstanding, we conclude in recommending pressure measurement via nasal cannula as the best method for measuring sniffing concurrent with olfactory tasks in humans.


    Appendix
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Appendix
 Acknowledgements
 References
 
Formula


    Acknowledgements
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Appendix
 Acknowledgements
 References
 
This work was supported by a grant from the Army Office of Research 46666-LS. We thank Arak Elite.


    References
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Appendix
 Acknowledgements
 References
 
Allison VC and Katz SH. (1919) An investigation of stenches and odors for industrial purposes. J Ind Eng Chem-Us 11:336–8.[CrossRef]

Banzett RB, Mahan ST, Garner DM, Brughera A, Loring SH. (1995) A simple and reliable method to calibrate respiratory magnetometers and respitrace. J Appl Physiol 79:2169–76.[Abstract/Free Full Text]

Benignus VA and Prah JD. (1980) A computer-controlled vapor-dilution olfactometer. Behav Res Methods 12:535–40.[Web of Science]

Bensafi M, Porter J, Pouliot S, Mainland J, Johnson B, Zelano C, Young N, Bremner E, Aframian D, Khan R, et al. (2003) Olfactomotor activity during imagery mimics that during perception. Nat Neurosci 6:1142–4.[CrossRef][Web of Science][Medline]

Buonviso N, Amat C, Litaudon P. (2006) Respiratory modulation of olfactory neurons in the rodent brain. Chem Senses 31:145–54.[Abstract/Free Full Text]

Chadha TS, Watson H, Birch S, Jenouri GA, Schneider AW, Cohn MA, Sackner MA. (1982) Validation of respiratory inductive plethysmography using different calibration procedures. Am Rev Respir Dis 125:644–9.[Web of Science][Medline]

Doty RL. (1975) An examination of relationships between the pleasantness, intensity and concentration of 10 odorous stimuli. Percept Psychophysiol 17:492–6.[Web of Science]

Efron B and Tibshirani R. (1993) An introduction to the bootstrap. (Chapman and Hall, London).

Elsberg CA and Levy I. (1935) The sense of smell. I. A new and simple method of quantative olfactometry. Bull Neurol Inst NY 4:5–19.

Frank RA, Dulay MF, Gesteland RC. (2003) Assessment of the sniff magnitude test as a clinical test of olfactory function. Physiol Behav 78:195–204.[CrossRef][Medline]

Gamble RC. (1898) The applicability of Weber's Lab to smell. Am J Psychol 10:82–142.[CrossRef]

Gilbert R, Auchincloss JH, Brodsky J, Boden W. (1972) Changes in tidal volume, frequency and ventilation induced by their measurement. J Appl Physiol 33:252–4.[Free Full Text]

Gonzalez H, Haller B, Watson HL, Sackner MA. (1984) Accuracy of respiratory inductive plethysmograph over wide range of rib cage abdominal compartmental contributions to tidal volume in normal subjects and in patients with chronic obstructive pulmonary disease. Am Rev Respir Dis 130:171–4.[Web of Science][Medline]

Hahn I, Scherer PW, Mozell MM. (1994) A mass transport model of olfaction. J Theor Biol 167:115–28.[CrossRef][Web of Science][Medline]

Han JN, Stegen K, Cauberghs M, Van de Woestijne KP. (1997) Influence of awareness of the recording of breathing on respiratory pattern in healthy humans. Eur Respir J 10:161–6.[Abstract]

Hasegawa M and Kern EB. (1977) The human nasal cycle. Mayo Clin Proc 52:28–34.[Web of Science][Medline]

Johnson BN, Mainland JD, Sobel N. (2003) Rapid olfactory processing implicates subcortical control of an olfactomotor system. J Neurophysiol 90:1084–94.[Abstract/Free Full Text]

Johnston JW. (1967) Quantification of olfactory stimuli. In Brown CC (Ed.). Methods in psychophysiology(The Williams & Wilkins Company, Baltimore, MD).

Kepecs A, Uchida N, Mainen ZF. (2006) The sniff as a unit of olfactory processing. Chem Senses 31:167–79.[Abstract/Free Full Text]

Keyhani K, Scherer PW, Mozell MM. (1995) Numerical simulation of airflow in the human nasal cavity. J Biomech Eng 117:429–41.[Web of Science][Medline]

Kobal G. (1985) Pain-related electrical potentials of the human nasal mucosa elicited by chemical stimulation. Pain 22:63.

Konno K and Mead J. (1967) Measurement of separate volume changes of rib cage and abdomen during breathing. J Appl Physiol 22:407–22.[Free Full Text]

Laing DG. (1982) Characterisation of human behaviour during odour perception. Perception 11:221–3.[Web of Science][Medline]

Laing DG. (1983) Natural sniffing gives optimum odour perception for humans. Perception 12:99–117.[Web of Science][Medline]

Laing DG. (1985) Optimum perception of odor intensity by humans. Physiol Behav 34:569–74.[CrossRef][Medline]

Laing DG. (1986) Identification of single dissimilar odors is achieved by humans with a single sniff. Physiol Behav 37:163–70.[CrossRef][Medline]

Laurent G. (2002) Olfactory network dynamics and the coding of multidimensional signals. Nat Rev Neurosci 11:884–95.

Le Magnen J. (1945) Etude des facteurs dynamiques de l'excitation olfactive. L'Année Psychologique 6:77–89.

Mainland J and Sobel N. (2006) The sniff is part of the olfactory percept. Chem Senses 31:181–96.[Abstract/Free Full Text]

Milledge J and Stott F. (1977) Inductive plethysmograph—a new respiratory transducer. J Physiol 267:4–5.

Paek D and McCool FD. (1992) Breathing patterns during varied activities. J Appl Physiol 73:887–93.[Abstract/Free Full Text]

Paulsen E. (1882) Experimentelle untersuchungen uber Stromungen der Luft in der Nasenhohle. Sitzungsber Akad Wiss.

Perez W and Tobin MJ. (1985) Separation of factors responsible for change in breathing pattern induced by instrumentation. J Appl Physiol 59:1515–20.[Abstract/Free Full Text]

Prah D, Sears S, Walker JC. (1995) Modern approaches to air dilution olfactometry. In Doty RL (Ed.). Handbook of olfaction and gustation(Marcel Dekker, Inc, New York) pp. 227–57.

Principato JJ and Ozenberger JM. (1970) Cyclical changes in nasal resistance. Arch Otolaryngol 91:71–7.[Abstract/Free Full Text]

Rodenstein DO, Mercenier C, Stanescu DC. (1985) Influence of the respiratory route on the resting breathing pattern in humans. Am Rev Respir Dis 131:163–66.[Web of Science][Medline]

Schoenfeld TA and Cleland TA. (2006) Anatomical contributions to odorant sampling and representation in rodents: zoning in on sniffing behavior. Chem Senses 31:131–44.[Abstract/Free Full Text]

Sobel N, Prabhakaran V, Desmond JE, Glover GH, Goode RL, Sullivan EV, Gabrieli JD. (1998) Sniffing and smelling: separate subsystems in the human olfactory cortex. Nature 392:282–6.[CrossRef][Medline]

Sobel N, Thomason ME, Stappen I, Tanner CM, Tetrud JW, Bower JM, Sullivan EV, Gabrieli JD. (2001) An impairment in sniffing contributes to the olfactory impairment in Parkinson's disease. Proc Natl Acad Sci USA 98:4154–9.[Abstract/Free Full Text]

Spors H and Grinvald A. (2002) Spatio-temporal dynamics of odor representations in the mammalian olfactory bulb. Neuron 34:301–15.[CrossRef][Web of Science][Medline]

Tripp H and Bolton J. (1998) Phrenic nerve injury following cardiac surgery: a review. J Cardiothorac Surg 31:1265–7.

Tucker D. (1963) Physical variables in the olfactory stimulation process. Gen Physiol 46:453–89.[Abstract/Free Full Text]

Uchida N and Mainen ZF. (2003) Speed and accuracy of olfactory discrimination in the rat. Nat Neurosci 6:1224–9.[CrossRef][Web of Science][Medline]

Utenick MR. (1971) Design of a hot-wire anemometer. ISA Trans 10:21–8.[Medline]

Walker JC, Kendal-Reed M, Hall SB, Morgan WT, Polyakov VV, Lutz RW. (2001) Human responses to propionic acid. II. Quantification of breathing responses and their relationship to perception. Chem Senses 26:351–8.[Abstract/Free Full Text]

Warren DW, Walker JC, Drake AF, Lutz RW. (1994) Effects of odorants and irritants on respiratory behavior. Laryngoscope 104:623–6.[Web of Science][Medline]

Watson HL, Poole DA, Sackner MA. (1988) Accuracy of respiratory inductive plethysmographic cross-sectional areas. J Appl Physiol 65:306–8.[Abstract/Free Full Text]

Willis BC, Graham AS, Wetzel R, Newth CJL. (2004) Respiratory inductance plethysmography used to diagnose bilateral diaphragmatic paralysis: a case report. Pediatr Crit Care Med 5:399–402.[CrossRef][Medline]

Youngentob SL, Mozell MM, Sheehe PR, Hornung DE. (1987) A quantitative analysis of sniffing strategies in rats performing odor detection tasks. Physiol Behav 41:59–69.[CrossRef][Medline]

Zhao K, Dalton P, Yang GC, Scherer PW. (2006) Numerical modeling of turbulent and laminar airflow and odorant transport during sniffing in the human and rat nose. Chem Senses 31:107–18.[Abstract/Free Full Text]

Zhao K, Scherer PW, Hajiloo SA, Dalton P. (2004) Effect of anatomy on human nasal air flow and odorant transport patterns: implications for olfaction. Chem Senses 29:365–79.[Abstract/Free Full Text]

Accepted 24 July 2006


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
Chem SensesHome page
A. Arzi, L. Sela, A. Green, G. Givaty, Y. Dagan, and N. Sobel
The Influence of Odorants on Respiratory Patterns in Sleep
Chem Senses, January 1, 2010; 35(1): 31 - 40.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
A. Koulakov, A. Gelperin, and D. Rinberg
Olfactory Coding With All-or-Nothing Glomeruli
J Neurophysiol, December 1, 2007; 98(6): 3134 - 3142.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
C. Zelano, J. Montag, B. Johnson, R. Khan, and N. Sobel
Dissociated Representations of Irritation and Valence in Human Primary Olfactory Cortex
J Neurophysiol, March 1, 2007; 97(3): 1969 - 1976.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
31/9/795    most recent
bjl021v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (3)
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Johnson, B. N.
Right arrow Articles by Sobel, N.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Johnson, B. N.
Right arrow Articles by Sobel, N.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?