$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
In this study, we conducted a detailed analysis of the respiratory patterns of mouse coughing and sneezing behaviors using whole-body plethysmography (WBP) and ultrasound recording techniques. The results showed that mice subjected to tracheal or nasal surgery demonstrate comparable cough and sneeze responsiveness to naive controls. However, capsaicin (chemical) and mechanical stimulation effectively evoke cough and sneeze reflexes in these surgically prepared animals (Figure 2A,B). WBP results revealed that both behaviors exhibited single- and double-peak patterns, with their sound signals appearing during the expulsive stage of the respiratory pattern (Figure 3A,C,E,G). Further ultrasound analysis revealed distinct acoustic profiles: Cough sounds exhibited abrupt onset with peak intensity occurring instantaneously, predominantly within 0-30 kHz frequency range, and shorter duration (Figure 3B,D; Supporting Video S1). Sneeze sounds demonstrated progressive intensity buildup to maximum amplitude, broadband frequency distribution (0-80 kHz), and prolonged duration (Figure 3F,H; Supporting Video S2). Statistical analysis revealed that the average compression phase of single-peak coughing behavior was 32.39 ± 0.73 ms, significantly shorter than that of double-peak coughing behavior at 53.17 ± 2.29 ms (Figure 4A,B,C). The average compression phase of single-peak sneezing behavior was 44.03 ± 1.84 ms, lower than the average value of double-peak coughing behavior, while the average compression phase of double-peak sneezing behavior was 74.24 ± 1.29 ms, significantly longer than that of single-peak sneezing behavior (Figure 4E,F,G). Additionally, the occurrence ratio of single-peak pattern in coughing behavior was 77.55 ± 2.10% (Figure 4D), while in sneezing behavior, the single-peak pattern accounted for only 32.42 ± 5.88%, with the double-peak pattern making up 67.58 ± 5.88% (Figure 4H). The statistical analysis of audio duration indicated significant differences between coughing and sneezing audio times, but there was considerable overlap between the two (Figure 5A). Through spectral analysis of the ultrasound audio of coughing and sneezing, we found that the frequency spectrum of coughing audio was mainly concentrated in the 0-30 kHz range, while sneezing audio had a higher proportion in the 0-60 kHz range (Figure 5B). Although coughing audio also had frequencies in the 30-120 kHz range, the proportion of low frequencies was always higher than that of high frequencies. Therefore, we used the slope of sound intensity from the starting point to the maximum value (Figure 5C) as the vertical coordinate and the energy ratio of 30-60 kHz to 0-30 kHz (Figure 5D) as the horizontal coordinate to establish a coordinate system, which could effectively distinguish coughing and sneezing audio. In addition, there were no significant differences between single- and double-peak coughing as well as single- and double-peak sneezing in this coordinate system (Figure 5E).
Next, we employed the aforementioned classification methods to categorize the coughing and sneezing behaviors induced by capsaicin nebulization in mice. We found that the coughing and sneezing behaviors induced by capsaicin nebulization overlap with those evoked by specific stimulation on the same coordinate positions (Figure 5F,G). Following transection of nasal sensory afferent nerves (anterior ethmoidal and infraorbital nerves), capsaicin nebulization significantly reduced sneeze-associated acoustic signals while cough-related audio profiles remained unaffected (Figure 5H,I).

Figure 1: Schematic diagram of the operation in the stimulated trachea and nasal cavity model. (A) Schematic representation of the catheter implanted into the trachea.(B) Frontal (left) and coronal (right) planes show the position of the catheter implanted into the trachea of the mouse. (C) The sagittal plane shows the position of catheter implantation in the nasal cavity and the position of filament fiber stimulation. Please click here to view a larger version of this figure.

Figure 2: Quantification of cough- and sneeze-like behaviors elicited by site-specific stimulation. (A) Sneeze quantification of the normal control group, sham group, and mechanical nasal stimulation group. The sham group exhibited comparable sneeze frequency to controls. Mechanical nasal stimulation evoked a significant increase in sneeze events. (B) Cough quantification of the normal control group, the sham group, and the capsaicin tracheal challenge group. The sham group showed no significant difference in cough frequency compared to controls. Capsaicin tracheal stimulation significantly increased cough counts. Data are represented as the mean ± SEM, ***P < 0.001, ****P < 0.0001. Please click here to view a larger version of this figure.

Figure 3: The respiratory traces and sound oscillogram of coughing and sneezing. (A) Representative images showing the single peak respiratory trace and audio recording during a cough induced by capsaicin. (B) The sound oscillogram (top) and spectrogram (bottom) of the cough audio recording in A. (C) Representative images showing the double peaks respiratory trace and audio recording during a cough induced by capsaicin. (D) The sound oscillogram (top) and spectrogram (bottom) of cough audio recording in C. (E) Representative images showing the single peak respiratory trace and audio recording during a sneeze induced by mechanical stimulation. (F) The sound oscillogram (top) and spectrogram (bottom) of the sneeze audio recording in E. (G) Representative images showing the double peaks respiratory trace and audio recording during a sneeze induced by mechanical stimulation. (H) The sound oscillogram (top) and spectrogram (bottom) of the sneeze audio recording in G. Please click here to view a larger version of this figure.

Figure 4: The single and double peaks of the respiratory traces of cough-like and sneeze-like behaviors. (A) Average (red) and individual (gray) and average (red) respiratory traces of cough-like behaviors with a single respiratory peak (n = 52).(B) Average (red) and individual (gray) and average (red) respiratory traces of cough-like behaviors with double respiratory peaks (n = 29).(C) Duration of single and double respiratory peaks compressive phases for coughs induced by capsaicin (single peak n = 52, double peaks n = 29). (D) The average ratio of single and double respiratory peak traces in cough (n = 6 mice). (E) Average (red) and individual (gray) and average (red) respiratory traces of sneeze-like behaviors with single respiratory peaks (n = 40). (F) Average (red) and individual (gray) and average (red) respiratory traces of sneeze-like behaviors with double respiratory peaks (n = 74). (G) Duration of single and double respiratory peaks compressive phases for sneezes induced by mechanical stimulation (single peak n = 40, double peaks n = 74). (H) The average ratio of single and double respiratory peak traces in sneeze (n = 6 mice). Data are represented as the mean ± SEM, ****P < 0.0001. Please click here to view a larger version of this figure.

Figure 5: The frequency and waveform classification of coughing and sneezing sound oscillogram. (A) Quantification of sonogram durations of coughs and sneezes.(B) The average spectrogram analysis of mechanical stimulation-induced sneezes (n = 69, black) and capsaicin-induced coughs (n = 140, red).(C) A representative schematic diagram for calculating the slope of a sound oscillogram. Comparative analysis of the smoothing spline reveals that at 25 points, excessive smoothing obscures original data features, whereas values exceeding 50 points (e.g., 75 points) show substantial overlap with the 50-spline curve but introduce additional inflection points.(D) A representative schematic diagram of area ratio calculation, using spectral analysis results to calculate the area under different frequency curves in GraphPad Prism. The area ratio is equal to the area under the 30-60 kHz curve and the area under the 0-30 kHz curve.(E) The frequency and waveform classification of mechanical stimulation-induced sneezes and capsaicin-induced coughs sound oscillogram, with the y-axis representing the exponent of the slope of the sound waveform from baseline to peak, and the x-axis representing the exponent of the ratio of sound waves between 30-60 kHz and 0-30 kHz.(F) The average spectrogram analysis of capsaicin nebulization-induced sneezes (n = 36, black) and capsaicin nebulization-induced coughs (n = 42, red). (G) The frequency and waveform classification of capsaicin nebulization-induced sneezes (n = 36, black) and capsaicin nebulization-induced coughs (n = 42, red) sound oscillogram, with the y-axis representing the exponent of the slope of the sound waveform from baseline to peak, and the x-axis representing the exponent of the ratio of sound waves between 30-60 kHz and 0-30 kHz. (H) Statistical data on capsaicin-induced sneeze reflex in groups with transected anterior ethmoidal and infraorbital nerves versus the sham surgery group. (I) Statistical data on capsaicin-induced cough reflex in groups with transected anterior ethmoidal and infraorbital nerves versus the sham surgery group. Data are represented as the mean ± SEM, **P < 0.01, ****P < 0.0001. Please click here to view a larger version of this figure.
Supplementary Video 1: Video of cough-like behavior synchronized with its audio spectrogram, displayed at normal and 0.1x speed. Please click here to download this File.
Supplementary Video 2: Video of sneeze-like behavior synchronized with its audio spectrogram, displayed at normal and 0.1x speed. Please click here to download this File.