Detection threshold for inspiratory resistive loads and respiratory-related evoked potentials
Abstract
The relationship between detection threshold of inspiratory resistive loads and the peaks of the respiratory-related evoked potential (RREP) is unknown. It was hypothesized that the short-latency and long-latency peaks of the RREP would only be elicited by inspiratory loads that exceeded the detection threshold. The detection threshold for inspiratory resistive loads was measured in healthy subjects with inspiratory-interruption or onset load presentations. In a separate protocol, the RREPs were recorded with resistive loads that spanned the detection threshold. The loads were presented in stimulus attend and ignore sessions. Onset and interruption load presentations had the same resistive load detection threshold. The P1, Nf, and N1 peaks of the RREP were observed with loads that exceeded the detection threshold in both attend and ignore conditions. The P300 was present with loads that exceeded the detection threshold only in the attend condition. No RREP components were elicited with subthreshold loads. The P1, Nf, and P300 amplitudes varied with resistive load magnitude. The results support the hypothesis that there is a resistive load threshold for eliciting the RREPs. The amplitude of the RREP peaks vary as a function of load magnitude. The cognitive P300 RREP peak is present only for detectable loads and when the subject attends to the stimulus. The absence of the RREP with loads below the detection threshold and the presence of the RREP elicited by suprathreshold loads are consistent with the gating of these neural measures of respiratory mechanosensory information processing.
the respiratory-related evoked potential (RREP) has been elicited by inspiratory occlusion and graded, easily detected inspiratory resistive (R) loads (10, 14, 15, 17, 28, 35, 44, 51). A relationship between R load magnitude, magnitude estimation, and the short-latency P1 peak amplitude has been demonstrated (23, 24). There is also a report of a minimum threshold for eliciting the long-latency peaks [N2, P2, and P300 peaks (3)], but the authors were unable to record the short-latency P1 and N1 peaks because they had a limited unilateral scalp electrode montage, used only loads presented at the onset of the inspiration, and were unaware of the existence of the frontal Nf peak (14, 28). Thus there is no report of an inspiratory load threshold for eliciting the primary initial sensory peaks, P1 and Nf, of the RREP. However, it is known that presentation of no load does not elicit any peaks of the RREP (10, 14, 35, 51). This suggests a threshold for eliciting these early peaks of the RREP. Although it is clear that there is a relationship between the magnitude estimation of inspiratory resistive loads and the early peaks (Nf, P1, and N1) of the RREP (24), i.e., load perception and neural measures of cortical neural processing of load information, it is unknown whether the load must exceed the cognitive detection threshold for elicitation of these RREP measures of neural sensory processing. The central hypothesis for this study is, therefore, that the short-latency P1 and Nf sensory peaks of the RREP can be elicited only by R loads that exceed the R load detection threshold [the added external R load that has a 50% detection probability (ΔR50)].
The RREP cognitive, attention-related peak most commonly studied was the positive peak that occurs ∼300 ms after the application of the load, P300 (2, 3, 41, 43–45), and it was present when the subject attended (attempted to detect) to the load. It has been reported that the latency and amplitude of the P300 RREP component were correlated with cognitive attention and the stimulus magnitude (2, 3, 42). It has been further reported that the cognitive-related P300 peak was not elicited by loads below the detection threshold and no load (3). No-load control RREP recordings have been reported to not elicit the P1, Nf, or N1 peaks of the RREP (2, 3, 10, 20, 23, 24, 42). It is likely that there is a R load threshold for eliciting both the short-latency and long-latency peaks of the RREP. If the P1 and Nf peaks reflect the arrival of the respiratory load-related mechanosensory information and the P300 is related to the cognitive cortical processing of afferent information, then it was reasoned that the P300 peak would be present only if the P1 peak was also present. This means that the P300 peak should be present only for R loads above the ΔR50 threshold and if the subject attended to the load. We hypothesized that the P300 peak would be present when the P1 peak was also present with loads that exceeded the ΔR50 detection threshold during attend conditions.
Humans can easily detect mechanical loads added to inspiration. The detection threshold for R loads is defined as ΔR50 (45). It has been reported that the probability for R load detection is a function of the background resistance (Ro), i.e., the resistance of the airways plus the resistance of the breathing apparatus (46). The Weber-Fechner law predicts that the detection of a stimulus is a constant fraction of the baseline stimulus intensity (46). For R loads, the Weber fraction (ΔR50/R0) is ∼0.3 (46). Previous R load detection studies have delivered the loads at the onset the inspiration (8, 26, 46, 48). The standard protocol requires the subject to inspire against the load at the beginning of the inspiratory effort and continue to breathe against the load throughout the inspiration (7, 8, 16, 26, 48, 49). The subject signals detection of the load if the load is of sufficient magnitude for the subject to detect an added load. Zechman and Davenport (48) reported that the detection of the load is related to the breathing pattern with R load detection correlated with the inspiratory airflow (V̇i). Although it is known that R load detection has a Weber fraction of ∼0.3 and that near-threshold load detection is related to the airflow pattern, it is unknown whether R load detection thresholds are affected by applying the load as an interruption of the inspiration. The application of inspiratory loads by interruption is essential for reliably eliciting the RREP (35). This method of load presentation has been used to demonstrate a relationship between magnitude estimation of the load and the amplitude of the primary sensory peak (P1) amplitude (23, 24). In a previous study of R load magnitude estimation (25), the slopes of log-log relationship between R load magnitude and magnitude estimation were the same when the loads were presented either at the onset of an inspiration or as an interruption of an inspiration. Although the slopes were unaffected by the presentation method, the loads applied as an interruption of inspiration were estimated as smaller, shifting the relationship to the right and increasing the R load intercept (25). This suggests that the detection threshold increased with the presentation of inspiratory loads as an interruption of inspiration. It was essential to demonstrate the relationship between the load presentation methods and the load detection threshold to use the inspiratory-interruption method for determining the relationship between the detection threshold and the RREP. It was hypothesized ΔR50 would increase with a detection paradigm using R loads presented as an interruption of inspiration.
The present study was designed to initially determine the ΔR50 with loads presented at the onset of the inspiration and as an interruption of inspiration. Loads presented as an inspiratory interruption were then used to elicit the RREP. The RREP was recorded using R loads that were above and below the ΔR50 detection threshold determined by interruption of inspiration for each subject. The RREP was recorded under both ignore and attend conditions. The short-latency (P1, Nf, N1) and cognitive P300 peaks of the RREP were identified, and their amplitudes and latencies were determined for loads above and below the inspiratory interruption ΔR50.
METHODS
The protocol for this experiment was reviewed and approved by the Institutional Review Board of the University of Florida.
Subjects
Twelve healthy nonsmoking subjects were tested (5 men and 7 women). The mean age of subjects was 24.6 ± 4.5 yr (19–33 yr). Mean height and weight were 171.3 ± 12.8 cm and 68.7 ± 13.2 kg, respectively. The subjects were in good general health based on self-reported medical history of no chronic or acute neurological or respiratory disease. The subjects refrained from consuming alcohol and caffeinated beverages for at least 24 h before the study. Each subject was brought to the laboratory and informed of the nature of the study. The general nature of the experiment was described, and each subject's consent was obtained. The subjects participated in two experimental protocols on a single day. Protocol 1 consisted of the detection of R loads using two different load presentation methods: 1) R loads presented at the onset of inspiration and 2) R loads presented as an interruption of inspiration. Protocol 2 consisted of the recording of the RREP using R load magnitudes that were above and below the ΔR50 threshold. The loads used in Protocol 2 were determined from the R detection trial with loads presented as an interruption of inspiration. Protocol 1 preceded Protocol 2 for all subjects.
Protocol 1: R Load Detection Threshold
Subject preparation.
Pulmonary function was measured. The total respiratory resistance was measured with the impulse oscillometry method (Jaeger MS-IOS, Sensormedics). Forced vital capacity (FVC) and forced expiratory volume in 1 s (FEV1) were recorded for all subjects with same apparatus. Three FVC measurements were performed on each subject, and the highest value was used for analysis. A FEV1 above 70% of predicted value was required to admit into the study. All subjects met this inclusion criterion. Pulmonary function was within the normal limits for all subjects. The results of pulmonary function tests are shown in Table 1.
| Measured Value | %Predicted | |
|---|---|---|
| FEV1.0, liters | 3.62±0.25 | 96.45±4.00 |
| FVC, liters | 4.28±0.37 | 96.93±4.41 |
| FEV1.0/FVC | 84.42±1.16% | |
| Raw, cmH2O·l−1·s | 4.06±0.86 |
The subject was then seated in a chair with the head, neck, and back supported. The subject respired through a mouthpiece and nonrebreathing valve (2600 series, Hans Rudolph). Care was taken to suspend the valve to eliminate the need for the subject to contract mouth muscles yet maintain an airtight seal. The subject wore a nose clip. The R loads were sintered bronze disks placed in series in a Plexiglas tube (loading manifold), with stoppered ports between disks. The loading manifold was connected to a pneumotachograph, which was connected by reinforced tubing to the inspiratory port of the nonrebreathing valve. Resistances were selected by removing the stopper and allowing the subject to inspire through the selected port. The background resistance of the system (Rapp) was 1.85 cmH2O·l−1·s. The loading manifold was hidden from the subject's view. Mouth pressure (Pm) was recorded from a port in the center of the non-rebreathing valve using a differential pressure transducer and a signal conditioner. V̇i was recorded using a differential pressure transducer connected to the pneumotachograph. V̇i was integrated to obtain the inspired volume (Vi). Pm, V̇i, and Vi were digitalized at 200 Hz (PowerLab/8s, ADInstruments, Mountain View, CA) and stored on the computer for later analysis. The subjects were monitored with a video camera through the study.
Protocol.
The subject was seated in a sound isolated chamber, separated from the investigator and experimental apparatus. Similar to our laboratory's previous R load detection studies (16, 48, 49), a series of suprathreshold and near-threshold R loads were presented in practice session to familiarize the subject with the load sensation and the range of the loads. The subjects were instructed to press the signal button if they detected a load. The signal button was held in their dominant hand and pressed with their thumb.
During the experimental trials, the subjects listened to the music of their choice, which masked experimental sounds. Two methods of presentation were used: 1) R loads presented at the onset of inspiration and 2) R loads presented as an interruption of inspiration. With inspiratory interruption, the R loads were applied after the onset of inspiratory phase, at the beginning of rising phase of V̇i as indicated by the analog display of the airflow signal on the oscilloscope. Onset R load presentations were applied by selecting the load during expiration resulting in the following inspiration against the load beginning at the onset of inspiration. The sequence of load presentation methods was randomized. Eight R load magnitudes (0.20, 0.80, 1.24, 1.64, 2.48, 3.26, 6.95, and 11.46 cmH2O·l−1·s) were used. Each R load was presented 10 times in a randomized block order, with each individual load presentation separated by 3–6 control breaths. The detection signal was recorded on the computer polygraph system indicating detection of the load. Thus there were 2 experimental trials with 80 loads presented in each trial. The trials differed in the method of load presentation. The subject was allowed a minimum of 5-min rest period off the breathing apparatus between the two trials. The total duration of the Protocol 1 experiment, including preparation time, was ∼1 h.
Data analysis.
All experimental data are presented as means ± SD. The detection threshold was determined separately for each experimental trial. The number of detected loads was summed for each load magnitude. The number of detections was divided by the number of presentations, and this percent detection was plotted against R load magnitude (ΔR). The detection threshold was defined as the R load magnitude corresponding to 50% detection (ΔR50). The ΔR50 was then divided by the total background resistance (Ro), which was the sum of the resistance of the breathing apparatus (Rapp) plus the subject's total airway resistance. The resultant ΔR50/Ro was determined for each subject. The ΔR50 and the ΔR50/Ro were also measured for each load presentation method. Comparison of ΔR50 and the ΔR50/Ro between the two trials was performed by using a two-tailed paired t-test. Pearson correlation was performed to test the correlation between two threshold measurements. The significance criterion was set at P < 0.05.
The peak Pm for each R load magnitude was measured from the computer polygraph recordings. The peak Pm for the two load presentation methods was analyzed with Tukey's honestly significant difference test. The significance criterion was set at P < 0.05.
Protocol 2: RREP
Subject preparation.
Protocol 2 was initiated after the completion of the detection study. The methods used in this part of the study were similar to previous RREP recordings with graded R loads (23). The subject was allowed at least a 15-min rest after finishing Protocol 1. The subject was tested in a separate room configured for RREP recordings. Tin electrodes were placed on both ears. The circumference of the head was measured, and an appropriately sized electrode cap with integral tin electrodes was placed on the head. EEG activity was recorded from 12 scalp sites: F3, FZ, F4, C3, CZ, C4, C3′, CZ′, C4′, P3, PZ, and P4 based on the International 10–20 System. EEG activity was referenced to the joined ear lobes. Vertical electrooculogram (EOG) activity was monitored with bipolar electrodes placed lateral to the canthus of the right eye. The impedance level for each electrode is checked and maintained below 5 kΩ. The EEG activity was monitored with an oscilloscope monitor. The electrode cap was connected to an electroencephalograph system (Neurodata 12, Teledyne). The EEG activity was band-pass filtered at 0.3 Hz to 1 kHz, amplified 50 k, led into an online signal averaging computer system, and digitized at 2.5 kHz (model 1401, Cambridge Electronics Design). The Pm signal was recorded by a differential-pressure transducer (model MP-45, Validyne Engineering), digitized, and led into the online signal-averaging computer system. The Pm was also led to an oscilloscope and monitored by the investigator. Pm was monitored continuously and used to time the presentation of the R loads.
A mouthpiece and non-rebreathing valve were connected to a breathing apparatus similar to that used for the detection trials. The inspiratory port of the non-rebreathing valve was connected to the loading manifold by reinforced tubing. A custom-designed gas-activated occlusion valve was connected to the first port of the loading manifold. Activation of the occlusion valve shunted V̇i through the selected R load. Three R loads were used for testing that spanned the detection threshold. The R1 load (0.80 cmH2O·l−1·s) was below the ΔR50 for all but one subject. This subject had a ΔR50 of 0.60 cmH2O·l−1·s, and the R1 load for this subject was 0.20 cmH2O·l−1·s. The R1 loads for all subjects were ∼30% of the ΔR50. The R2 load (6.95 cmH2O·l−1·s) was above the ΔR50, and >50% of the R2 load presentations were detected by all subjects in the inspiratory-interruption detection trials. The R3 load (11.46 cmH2O·l−1·s) was easily detected with 100% detection for inspiratory-interruption presentations in all subjects. In addition, Ro was presented. The subject wore a nose clip when respiring on the breathing apparatus. The loading manifold was hidden from the subject's view. Pm was recorded from a port in the center of the non-rebreathing valve. During the experiment, the inspiratory load was presented by silently inflating the occlusion valve after the onset of the inspiration, after the onset of the rising phase of inspiratory airflow as indicated by the analog display of the airflow signal on the oscilloscope. Each loaded breath was separated by three to six unloaded breaths. A transistor-transistor logic pulse generated by the inspiratory occlusion valve controller triggered the collection of 50 ms of pretrigger and 950 ms of posttrigger EEG and Pm data. The duration of the load presentation was ∼500 ms.
Protocol.
Subjects were seated comfortably in a reclining lounge chair with their back, neck, and head supported. There were two conditions: ignore and attend presented with a protocol similar to our laboratory's previous studies (51). For each condition, the four inspiratory load levels (Ro, R1, R2, and R3) were presented in an order to minimize temporal, order, and sequence effects as described previously (23, 24). Subjects received 1 trial of 200 presentations for each condition, with a rest period between each trial. Within each trial, load presentations were tested using 40 blocks of 5 presentations each, with the same load presented within a block to minimize manifold manipulations. Ordering of the different blocks was randomly determined with the restriction that each successive set of four blocks contained only one block of each load level. Thus, over the trial of 200 presentations, each block of 5 presentations of a given load was presented 10 times for 50 presentations of each load magnitude. There was 1 trial for each condition, making 50 presentations of each load magnitude (and no load) available for signal averaging.
The order of the two conditions, ignore and attend, was randomized. During the ignore condition, the subject watched a videotape movie. During the attend condition, the subject listened to music and signaled detection of a load similar to Protocol 1. With the attend condition, the subject was asked to signal the detection of the R load to ensure that attention to sensing the load was maintained. However, because this task was designed to maintain the attention rather than assess the detection performance, these detection data were not analyzed and were discarded. Including subject preparation time, the entire study session for Protocol 2 lasted ∼3 h.
Data analysis.
For each load presentation, a 1,000-ms epoch of EEG activity and Pm was stored on disk for subsequent computer signal averaging (Signal 2, Cambridge Electronic Design). Each load magnitude was averaged separately. The computer stored the individual load trials in the order of presentation. The averaged response for an individual load magnitude was obtained as described previously (23). An individual presentation was included in the average if there was 1) a stable prestimulus EEG activity baseline, 2) no EEG transient activity exceeding ±50 μV, 3) no EOG eye blink, and 4) the onset of the Pm change related to the load aligned with previous presentations. This analysis was repeated for each load magnitude and the no-load control. There were no peaks in the no-load control average that corresponded to the peaks observed with the load-elicited RREP. The EEG activity was averaged independently for each load magnitude for each condition. The presence, latency, and amplitude of RREP components (P1, Nf, N1, and P300) were determined from the averaged EEG traces. Based on previous peak localization studies (14, 28, 41, 42), the Nf peak was analyzed in the F3 and F4 electrodes, the P1 peak was analyzed from the C3′ and C4′ electrodes, the N1 peak was analyzed from the CZ electrode, and the P300 peak was analyzed from the CZ and PZ electrodes. Peak latencies were measured as the time from the onset of the occlusion (indicated by the change in Pm) to the EEG peak. The zero-to-peak amplitude was recorded at the peak of each component. The nomenclature for the peaks is based on previous reports (10, 14, 23, 24, 28, 41–43). The RREP components were identified as follows: Nf was the negative peak in the frontal region occurring within a latency range of 25–45 ms; P1 was the positive peak in the central region occurring within a latency range of 45–65 ms; N1 was the negative peak in the vertex (CZ) occurring within a latency range of 85–125 ms; and P300 was the positive peak in the CZ and PZ scalp locations occurring within a latency range of 250–350 ms. Averaged RREPs were obtained for each load magnitude, each condition, and each subject.
The descriptive statistics of all the variables were calculated and are expressed as means ± SD. The differences in latencies and amplitude of short-latency RREP components were compared by using repeated-measures ANOVA and Tukey's honestly significant difference test to determine the effects of load magnitude and condition. Because no P300 peak was present under ignore condition, the repeated-measures ANOVA was performed to determine the difference in amplitude and latency between different load magnitudes. Significance level was set at P < 0.05.
RESULTS
V̇i decreased, whereas Pm became more negative, compared with control breathing when the R loads were delivered before the onset of inspiration. The pressure became more negative as the load magnitude increased (Fig. 1). There were no significant differences in the peak Pm between onset and interruption load presentations. The mean detection threshold (Fig. 2) for the onset-inspiration method was 2.64 ± 0.89 cmH2O·l−1·s (range = 1.24–4.35 cmH2O·l−1·s). The mean ΔR50/Ro for loads presented at the onset of the breath was 0.45 ± 0.13. For inspiratory-interruption delivered loads, V̇i decreased sharply after application of the load, and the Pm had a corresponding rapid decrease. The magnitude of the pressure change increased with increasing load magnitude (Fig. 1). The mean detection threshold for inspiratory interruption presentations (Fig. 2) was 2.66 ± 1.22 cmH2O·l−1·s (range = 0.60–5.10 cmH2O·l−1·s). The mean ΔR50/Ro for loads presented interruption of the breath was 0.44 ± 0.18. There was no significant difference in detection threshold between two load presentation methods (P > 0.05). A significant correlation was found in R load detection threshold between the two load presentation methods (r = 0.67, P < 0.05). There was a significant difference for percent detection among the different load magnitudes for both presentation methods (P < 0.001); as the load magnitude increased, the percent detection increased (Fig. 2). However, no significant interaction was found between presentation method and load magnitude percent detection.

Fig. 1.Relationship between resistive load magnitude, presentation method, and peak mouth pressure (Pm). ○, Group means and SDs for peak Pm for resistive loads presented as an interruption of inspiration; •, group means and SDs for peak Pm for resistive loads presented at the onset of inspiration. There were no significant differences between the peak Pm and presentation method.

Fig. 2.Relationship between % detection, presentation method, and resistive load magnitude. Solid line and data points, group means and SDs for the % detection of each inspiratory resistive (R) load magnitude presented at the onset of the inspiration; dashed line and data points, group means and SDs for the % detection of each R load magnitude presented by interruption of inspiration. Detection threshold [the added external R load that has a 50% detection probability (ΔR50)] is indicated by the arrow for both conditions, and there was no significant difference.
R Load RREP Threshold
The RREP was elicited with the R2 and R3 loads, which were above the ΔR50 of all subjects, during both ignore and attend trials. The R1 and no-load Ro did not elicit RREP peaks in any of the subjects (Fig. 3).

Fig. 3.Example of respiratory-related evoked potential (RREP) Nf peak response to R loads and no-load. The F3 traces are identical to the other frontal electrode recordings. The 4 RREP traces are 50 presentations for each load magnitude averaged for 1 subject; negative voltage is plotted up the y-axis. Loads were presented under ignore conditions. Nf peak is indicated in the R3 load trace. Pm traces for each load magnitude are presented with the first Pm trace Ro (no load), the second Pm trace the R1 load, the third Pm trace the R2 load, and the fourth Pm trace the R3 load.
RREP ignore condition.
The R2 and R3 loads elicited the Nf peak recorded in the F3 and F4 electrodes (Fig. 3). The mean Nf amplitudes for R2 and R3 are summarized in Table 2. The Nf amplitudes with the R3 load for F3 and F4 were significantly greater (P < 0.01) than for R2 (Fig. 4). There was no significant difference in the Nf latencies (Table 3) between R2 and R3 (Fig. 5). The P1 peak was elicited by the R2 and R3 loads and recorded in the C3′ and C4′ electrodes (Fig. 6). The mean P1 amplitudes are summarized in Table 2. The P1 amplitudes (Fig. 4) with the R3 load for C3′ and C4′ were significantly greater (P < 0.05) than for R2. There was no significant difference in the P1 latencies (Table 3) between R2 and R3 (Fig. 5). The R2 and R3 loads elicited the N1 peak recorded in the CZ electrode (Fig. 6). The mean N1 amplitudes (Table 2) were not significantly different between R2 and R3 (Fig. 4). There was also no significant difference in the N1 latencies (Table 3) between R2 and R3 (Fig. 5). The P300 peak was not observed with any of the R loads in the ignore condition.

Fig. 4.Relationship between R load magnitudes and RREP peak amplitudes in the ignore condition. Group means and SDs for peak amplitudes are presented for the electrode locations for each peak. ★ Significant difference between R2 and R3 peak amplitude, P < 0.05.

Fig. 5.Relationship between R load magnitudes and RREP peak latencies in the ignore condition. Group means and SDs for peak latencies are presented for each peak. Nf peak latency is from recording sites F3 and F4, P1 latencies are from recording sites C3′ and C4′ and N1 latency is from the CZ recording site. There were no significant differences between the R2 and R3 peak latencies.

Fig. 6.Example of RREP P1 and N1 peak response to R loads and no load. C3′ traces are identical to the other central electrode recordings. The 4 RREP traces are 50 presentations for each load magnitude averaged for 1 subject; negative voltage is plotted up the y-axis. Loads were presented under ignore conditions. P1 and N1 peaks are indicated in the R3 load trace. Pm traces for each load magnitude are presented with the first Pm trace Ro (no load), the second Pm trace the R1 load, the third Pm trace the R2 load, and the fourth Pm trace the R3 load.
| Nf | P1 | N1 | P300 | |||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| F3 | F4 | C3′ | C4′ | CZ | PZ | |||||||||||||||
| R2 | R3 | R2 | R3 | R2 | R3 | R2 | R3 | R2 | R3 | R2 | R3 | |||||||||
| Ignore, μV | −1.34±1.53 | −3.12±1.72 | −1.86±1.11 | −3.34±1.67 | 1.38±0.88 | 1.91±0.87 | 1.46±0.83 | 3.51±1.84 | −2.08±1.38 | −2.51±1.47 | ||||||||||
| Attend, μV | −0.94±0.65 | −2.57±1.05 | −1.30±0.94 | −2.12±0.87 | 1.17±0.88 | 1.96±1.21 | 1.51±0.78 | 3.28±0.96 | −2.60±1.82 | −4.42±2.38 | 7.42±4.78 | 11.12±5.31 | ||||||||
| Nf | P1 | N1 | P300 | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| R2 | R3 | R2 | R3 | R2 | R3 | R2 | R3 | |||||
| Ignore, ms | 40.6±7.7 | 42.7±10.2 | 68.9±12.8 | 69.4±17.2 | 119.5±24.2 | 114.5±26.4 | ||||||
| CV | 0.19 | 0.24 | 0.19 | 0.25 | 0.20 | 0.23 | ||||||
| Attend, ms | 43.2±9.8 | 42.3±5.8 | 73.6±25.4 | 73.1±24.1 | 120.5±32.2 | 131.0±29.4 | 312.1±21.5 | 306.1±22.7 | ||||
| CV | 0.23 | 0.14 | 0.34 | 0.33 | 0.27 | 0.22 | 0.07 | 0.07 | ||||
RREP attend condition.
Attention to the load resulted in a change in the RREP waveform with the presence of a large amplitude P300 peak in the attend condition (Fig. 7). The R2 and R3 loads elicited the Nf peak recorded in the F3 and F4 electrodes. The mean Nf amplitudes for R2 and R3 are summarized in Table 2. The Nf amplitudes with the R3 load for F3 and F4 were significantly greater (P < 0.05) than for R2 (Fig. 8) There was no, however, significant difference for the Nf amplitudes between the attend and ignore conditions (Table 2). There was also no significant difference in the Nf latencies (Table 3) between R2 and R3 (Fig. 9) and no significant difference for the Nf latencies between the attend and ignore conditions (Table 3). The R2 and R3 loads elicited the P1 peak recorded in the C3′ and C4′ electrodes. The P1 amplitudes (Table 2) with the R3 load for C3′ and C4′ were significantly greater (P < 0.05) than for R2 (Fig. 8). There was also no significant difference for the P1 amplitudes between attend and ignore conditions (Table 2). There was no significant difference in the P1 latencies (Table 3) between R2 and R3 (Fig. 9) and no significant difference for the P1 latencies between attend and ignore conditions (Table 3). The N1 peak was elicited by the R2 and R3 loads, recorded in the CZ electrode. The N1 amplitude (Table 2) with the R3 load for CZ was significantly greater (P < 0.01) than for R2 (Fig. 8). The N1 amplitude for the R2 load in the attend condition was not significantly different from the N1 amplitude for the R2 load in the ignore condition; however, the N1 amplitude (Table 2) was significantly (P < 0.01) greater for R3 in the attend condition than for the same load in the ignore condition (Fig. 8). There was no significant difference in the N1 latencies (Table 3) between R2 and R3 (Fig. 9) and no significant difference for the N1 latencies between attend and ignore conditions.

Fig. 7.Example of RREP P300 peak response elicited by an R3 load between ignore and attend conditions. PZ traces are from the same subject; negative voltage is plotted up the y-axis. The first trace is the ignore condition, and the second trace is the attend condition with the P300 peak indicated. RREP traces are 50 presentations for each load magnitude averaged for 1 subject for each condition. Pm was identical for both conditions and is presented in the third trace.

Fig. 8.Relationship between R load magnitude and RREP peak amplitudes in the attend condition. Group means and SDs for peak amplitudes are presented for the electrode locations for each peak. ★ Significant difference between R2 and R3 peak amplitude, P < 0.05.

Fig. 9.Relationship between R load magnitudes and RREP peak latencies in the attend condition. Group means and SDs for peak latencies are presented for each peak. Nf peak latency is from recording sites F3 and F4, P1 latencies are from recording sites C3′ and C4′, N1 latency is from the CZ recording site, and P300 latency is from the PZ recording site. There were no significant differences between the R2 and R3 peak latencies.
The P300 was observed with R2 and R3 loads only in the attend condition (Fig. 10). The R2 load elicited the P300 peak recorded in the PZ electrode. The P300 amplitudes (Table 2) with the R3 load for PZ was significantly greater (P < 0.01) than for R2 (Fig. 8). There was no significant difference in the P300 latencies (Table 3) between R2 and R3 (Fig. 9).

Fig. 10.Example of RREP P300 peak responses to R loads and no load. The 4 PZ RREP traces are 50 presentations for each load magnitude averaged for 1 subject; negative voltage is plotted up the y-axis. Loads were presented under attend conditions. P300 peak is indicated in the R3 load trace. Pm traces for each load magnitude are presented with the first Pm trace Ro (no load), the second Pm trace the R1 load, the third Pm trace the R2 load, and the fourth Pm trace the R3 load.
DISCUSSION
The results of this study demonstrate an R load threshold for eliciting the short- and long-latency peaks of the RREP. The RREP is only elicited with R loads that exceed the detection threshold. This study further demonstrated that the detection threshold and Weber ratio for R loads presented at the onset of the breath or as an interruption of inspiration are the same and rejects the hypothesis that the increase in the magnitude estimation intercept with inspiratory interruption loads reflects an increase in the detection threshold. This study supports the hypothesis that there is mechanical threshold for respiratory afferents to elicit cognitive awareness of an R load and that this cognitive threshold must be exceeded for the RREP neural measure of cerebral cortical activity to be present. The results further suggest that cognitive awareness of inspiratory loads is mediated by a central neural threshold gating process.
R Load Detection
R load detection is airflow dependent with the detection of near-threshold loads occurring at the point of peak airflow (26, 48). The generation of the inspiratory driving pressure creates the pressure gradient for air to flow into the lung. R loads oppose the movement of air and result in a reduction in airflow (30, 47). Pleural pressure becomes more negative, and the transdiaphragmatic pressure (Pdi) is augmented (49). The pattern of these mechanical changes is directly related to the airflow and the load magnitude. Resistive load magnitudes near the threshold of detection, presented at the onset of the breath (48), produce very small changes in airflow (from control), but they are associated with a more negative Pm and augmented Pdi (49). The Pdi becomes augmented above control progressively earlier in the inspiration as the load magnitude increases. This is also associated with a progressive increase in the detectability of the loads and decrease in the detection latency. This study (49) demonstrated that the detection of an R load was primarily dependent on the change in Pdi. There was a Pdi threshold which, when exceeded, resulted in signaling detection of the load at a constant latency. Thus the change in latency with the change in load magnitude was due to differences in the time required for Pdi to exceed the mechanical threshold in addition to a relatively constant central neural processing time for signaling detection. If it is the change in Pdi that is a major determinant of the detection of the load, then this may explain why the detection threshold for onset and interruption load presentations were the same. A specific load magnitude at a constant respiratory drive will generate the same change in Pdi. The load interrupting the inspiratory effort decreases the airflow and augments inspiratory Pm as a function of load magnitude not timing of the load presentation. The inspiratory interruption loads thus produce a similar magnitude of change in mechanical impediment to inspiration, resulting in the same change in Pm and Pdi, resulting in a similar detection threshold.
The detection of inspiratory loads is dependent on the background mechanical status of the subject (46). The 50% detection threshold for R loads is a constant ratio of the background load: ΔR50/R0 = k (where k is a constant) (46). The results of this study are consistent with previous reports of this relationship, and the ΔR50/R0 is not altered by the load presented at the onset of the breath or as an interruption of inspiration. In a previous study, however, R loads presented by inspiratory interruption reduced magnitude estimation of R load compared with the same loads presented at the onset of inspiration (25). This difference was a change in the magnitude estimation of all loads, yet the slope of log-log relationship between magnitude estimation and load magnitude was the same (25). This suggested that, whereas the absolute magnitude of the load was decreased with inspiratory interruption, the sensitivity of the subject to the load did not depend on load presentation method (25). This also means that the effect of inspiratory load pattern is different for detection and magnitude estimation of the load. Of importance to the present study, the R threshold for eliciting the RREP is not a function of the presentation method but due to the load detection threshold.
RREP Detection Threshold
The RREP has been reported to be elicited by inspiratory occlusion and R loads (2, 3, 10, 14, 17, 22–24, 28, 35, 37, 41–43, 51). The short-latency components P1, Nf and N1 reflect the activation of cortical neurons by afferent signals ascending somatosensory pathways to the sensory-motor cortex (11, 13, 17, 28, 50) and the supplementary motor area (28). These peaks of the RREP are determined mainly by the physical characteristics of the stimuli. In the present study, the RREP was only present with the presentation of suprathreshold R loads similar to the effect of mechanical tactile stimulation somatosensory-evoked potentials (36). The correlation between the load magnitude necessary to elicit detection and the load magnitude necessary to elicit the peaks of the RREP suggests that these peaks are neural indicators of cortical sensory processing of inspiratory mechanical information.
The first peak of the RREP, P1, was a positive voltage that was due to the dipole that occurs when a cerebral cortical column was depolarized by the arrival of activity from a population of afferents that were activated by the occlusion stimulus (10, 14, 17, 28). This first peak signals the arrival of the afferent information at the somatosensory cortex. Previous studies have confirmed that the amplitude of P1 peak is significantly correlated with R load magnitude and magnitude estimation (24, 42). This was confirmed in the present study with the amplitude of the P1 peak greater for R3 than the R2 load. The detection threshold for the P1 peak has not been previously reported. Bloch-Salisbury et al. (3) were unable to identify the P1 peak in their study of the R load threshold for the RREP. They used unilateral recording sites referenced to the mastoid processes. The N1 and P300 peaks can be reliably recorded with their electrode montage and reference, but it is very difficult to observe the P1 at their F3′, CZ, and P3′. In addition, they used only onset load presentations, which elicit smaller amplitude RREP peaks for P1 and N1. Onset load presentations have variable peak latencies that depend on the subject's respiratory drive (10). The peak latency should increase with decreasing load magnitude (48), making is very difficult for identifying and analyzing the onset P1 peak (3). Thus the present study extends their study (3) by using a recording electrode montage that will reliably express the P1 peak and a load presentation method (inspiratory interruption) that controls for the time required for the load to exceed the mechanoreceptor threshold. The results from this study also demonstrate that once the threshold for eliciting the P1 peak is exceeded, similar to our laboratory's previous results (23, 24), there is a direct relationship between the magnitude of the load and P1 peak amplitude, R3 was greater that R2. Thus, the first stage of neural processing is eliciting mechanosensory information of sufficient magnitude to activate somatosensory cortical neurons. When the loads exceeded the detection threshold, the P1 peak amplitude then correlates with the load magnitude and magnitude estimation of the load (23, 24).
The Nf peak was found in the frontal region (14, 28). The Nf peak was not reported with cephalic referenced RREPs (10, 35) nor with electrodes referenced to the mastoids (3). Logie et al. (28) found the source localization for the frontal Nf appears to be elicited by a sensory projection to the lateral premotor cortices. The Nf peak has received little specific analysis. It has been reported that the Nf peak is the result of a neural pathway that is a parallel activation to the somatosensory P1 activation (14). The Nf peak was elicited only with loads that exceeded the detection threshold. When the Nf peak was found with suprathreshold loads, the amplitude of the Nf peak increased with increasing load magnitude. This suggests that the Nf peak amplitude, similar to the P1 peak, is a function of the stimulus intensity. The role of this frontal cortical activation in R load detection and information processing, however, remains unknown.
The relationship between R load magnitude, N1 peak threshold, and N1 peak amplitude relationship to load magnitude has not been reported. A previous RREP threshold study (3) was unable to analyze the N1 peak because they could not reliably record the N1 peak using onset load presentations. In the present study, we demonstrate that the N1 peak was present with loads that exceed the detection threshold. In the absence of attention to the load (ignore condition), however, the amplitude of the N1 peak did not vary with suprathreshold load magnitudes. Localization studies have shown the greatest N1 amplitude at the vertex, CZ (40, 42). Harver et al. (22) reported that N1 was unaffected by attention in young subjects. However, Webster and Colrain (42) reported that the N1 peak amplitude increased with attention to inspiratory occlusions. These latter investigators suggested that the N1 peak is a compound peak consisting of two negative components, one unassociated with attention and a second attention related (41, 42). This is supported by reports that the auditory N100 peak is a composite of two negative components, one related to attention and the other insensitive to attention (21, 31–33, 39). The present study is consistent with a dual source for the N1 peak. There was a significant increase in N1 amplitude for the R3 load with subjects attending to the load. The R2 attend N1 amplitude was nonsignificantly greater than for R2 with the ignore condition. Thus the relationship between load magnitude and stimulus magnitude for the ignore condition did not reach significance primarily because of a smaller R3 elicited N1 peak in the ignore condition. When the subject attended to the loads, the N1 peak amplitude increased, and a significant load magnitude relationship was evident. These results suggest that the RREP N1 peak is a function of the stimulus parameters and is modulated by subject attention to the stimulus. These results are also consistent with the suggestion (32, 33) that the N100 peak is indicative of a triggering process related to the perception of a nonspecific stimulus.
The P300 peak is dependent on the subject's attention to the stimulus (41). The present study demonstrated that the P300 peak is absent with undetectable loads and if the P1 peak is not present. When the P300 peak was present with suprathreshold loads, there was a direct relationship between R load magnitude and P300 amplitude; R3 P300 amplitude was greater than for R2. These results are consistent with previous reports on the relationship between RREP P300 amplitude and R load magnitude (2, 3, 42). The P300 displayed a topographic distribution (9, 44) showing centroparietal dominance identical (29) to that seen in the same subjects in response to target stimuli in an auditory oddball task suggesting that the RREP P300 is similar in origin to the auditory P300. The present results support the load magnitude effect on the amplitude of RREP P300 peak. The R3 load elicited significantly greater P300 peak amplitude compared with the R2 load. The P300 was not present with loads below the detection threshold despite the fact that the subjects were attending to the experimental trial. Thus the RREP P300 peak is observed only when the R load is of sufficient magnitude to elicit somatosensory cortical activity and the subject attends to detecting the load. This result suggested that the presence of P300 peak is correlated with the presence of short-latency RREP components.
This relationship between short-latency RREP components and long-latency RREP components would be expected because respiratory mechanoreceptor information arrives at cerebral cortex, eliciting short-latency stimulus dependent neural activity. This initial cortical activation is then followed by attention related neural activity. The N1 peak may represent the triggering or gating process related to subjects attending to the stimulus, which is then followed by cortical and subcortical neural activity related to cognitive processing of the load information (P300). The pattern of RREP peaks, therefore, represents neural indicators of the temporal sequence of respiratory sensation.
Eupneic breathing normally is not consciously perceived. This suggests that during normal ventilation, respiratory sensory information is gated-out of cognitive centers. When ventilation is obstructed, stimulated, challenged, or attended to, cognitive awareness of breathing occurs (7, 12, 16, 34, 38, 46). The perception of an intrinsic change in respiratory mechanics is an essential component of the cognitive awareness of the onset of an obstructive event. The relationship between the stimulus threshold and central neural activation elicited by this respiratory mechanical stimulation suggests a gated process. Sensory gating has been demonstrated in auditory, visual and somatosensory modalities (1, 5, 6, 19, 27). Modality-specific activation of cortical neural processing centers depends on a change in neural activity that gates-in modality specific information to the brain information processing centers (1, 4–6, 18, 19, 27). This activation leads to cognitive awareness of the modality. The significance of gating-in and gating-out sensory modalities is the need to attend to essential interoceptive physiological functions such as breathing. Respiratory afferent input is normally gated-out from cognitive brain systems, and individuals breathe without conscious attention to their ventilation. Changing the status of the respiratory system alters the input of this information. Respiratory information is then gated-in, eliciting a cognitive awareness of breathing. The transition between breathing that does not reach consciousness to cognitive awareness of breathing suggests neural information processing that gates-in respiratory sensory information. The neural processing of this respiratory information is reflected by the temporal and spatial brain activities that are the foundation of the peaks of the RREP. Hence, the present study is consistent with respiratory cognitive awareness of breathing as a gated process. The absence of the RREP with loads below the detection threshold and the presence of the RREP elicited by suprathreshold loads are consistent with the gating of these neural measures of respiratory mechanosensory information processing.
GRANTS
This study was supported by National Heart, Lung, and Blood Institute Grant HL-48792.
FOOTNOTES
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
REFERENCES
- 1 Arnfred SM, Eder DN, Hemmingsen RP, Glenthøj BY, Chen ACN. Gating of the vertex somatosensory and auditory evoked potential P50 and the correlation to skin conductance orienting response in healthy men. Psychiatry Res 101: 221–235, 2001.
Crossref | ISI | Google Scholar - 2 Bloch-Salisbury E, Harver A. Effects of detection and classification of resistive and elastic loads on endogenous event-related potentials. J Appl Physiol 77: 1246–1255, 1994.
Link | ISI | Google Scholar - 3 Bloch-Salisbury E, Harver A, Squires NK. Event-related potentials to inspiratory flow-resistive loads in young adults: stimulus magnitude effects. Biol Psychol 49: 165–186, 1998.
Crossref | ISI | Google Scholar - 4 Boutros NN, Zouridakis G, Overall J. Replication and extension of P50 findings in schizophrenia. Clin Encephalogr 22: 40–45, 1991.
Google Scholar - 5 Boutros NN, Belger A. Midlatency evoked potentials attenuation and augmentation reflect different aspects of sensory gating. Biol Psychiatry 45: 917–922, 1999.
Crossref | ISI | Google Scholar - 6 Boutros NN, Torello MW, Barker BA, Tueting PA, Wu SC, Nasrallah HA. The P50 evoked potential component and mismatch detection in normal volunteers: implications for the study of sensory gating. Psychiatry Res 57: 83–88, 1995.
Crossref | ISI | Google Scholar - 7 Campbell EJM, Freedman S, Smith PS, Taylor ME. The ability of man to detect added elastic loads to breathing. Clin Sci (Lond) 20: 223–231, 1961.
ISI | Google Scholar - 8 Campbell EJM, Gandevia SC, Killian KJ, Mahutte CK, Rigg JR. Changes in the perception of inspiratory resistive loads during partial curarization. J Physiol 309: 93–100, 1980.
Crossref | ISI | Google Scholar - 9 Colrain I, Webster K, Hirst G. Multi-channel EEG analysis of respiratory and auditory evoked potentials while awake and in NREM sleep. J Sleep Res 5, Suppl 1: 39, 1996.
Google Scholar - 10 Davenport PW, Friedman WA, Thompson FJ, Franzen O. Respiratory-related cortical potentials evoked by inspiratory occlusion in humans. J Appl Physiol 60: 1843–1848, 1986.
Link | ISI | Google Scholar - 11 Davenport PW, Thompson FJ, Reep R, Freed A. Projection of phrenic nerve afferents to the cat sensorimotor cortex. Brain Res 328: 150–153, 1986.
ISI | Google Scholar - 12 Davenport PW, Dalziel DJ, Webb B, Bellah JR, Vierck CJ Jr. Inspiratory resistive load detection in conscious dogs. J Appl Physiol 70: 1284–1289, 1991.
Link | ISI | Google Scholar - 13 Davenport PW, Shannon R, Mercak A, Reep RL, Lindsey BG. Cerebral cortical evoked potentials elicited by cat intercostal muscle mechanoreceptors. J Appl Physiol 74: 799–804, 1993.
Link | ISI | Google Scholar - 14 Davenport PW, Colrain IM, Hill PM. Scalp topography of the short-latency components of the respiratory-related evoked potential in children. J Appl Physiol 80: 1785–1791, 1996.
Link | ISI | Google Scholar - 15 Davenport PW, Cruz M, Stecenko AA, Kifle Y. Respiratory-related evoked potentials in children with life-threatening asthma. Am J Respir Crit Care Med 161: 1830–1835, 2000.
Crossref | ISI | Google Scholar - 16 Davenport PW, Kifle Y. Inspiratory resistive load detection in children with life-threatening asthma. Pediatr Pulmonol 32: 44–48, 2001.
ISI | Google Scholar - 17 Davenport PW, Hutchison AA. Cerebral cortical respiratory-related evoked potentials elicited by inspiratory occlusion in lambs. J Appl Physiol 93: 31–36, 2002.
Link | ISI | Google Scholar - 18 Gandevia SC, Macefield G. Projection of low-threshold afferents from human intercostal muscles to the cerebral cortex. Respir Physiol 77: 203–214, 1989.
Crossref | PubMed | Google Scholar - 19 Grunwald T, Boutros NN, Pezer N, von Oertzen J, Fernandez G, Schaller C, Elger CE. Neuronal substrates of sensory gating within the human brain. Biol Psychiatry 53: 511–519, 2003.
Crossref | ISI | Google Scholar - 20 Hammond CS, Gaeta H, Sapienza C, Davenport PW. Respiratory-related evoked potential elicited by expiratory occlusion. J Appl Physiol 87: 835–842, 1999.
Link | ISI | Google Scholar - 21 Hansen JC, Hillyard SA. Endogenous brain potentials associated with selective auditory attention. Electroencephalogr Clin Neurophysiol 49: 277–290, 1980.
Crossref | PubMed | Google Scholar - 22 Harver A, Squires NK, Bloch-Salisbury E, Katkin ES. Event-related potentials to airway occlusion in young and old subjects. Psychophysiology 32: 121–129, 1995.
Crossref | ISI | Google Scholar - 23 Knafelc M, Davenport PW. Relationship between resistive loads and P1 peak of respiratory-related evoked potential. J Appl Physiol 83: 918–926, 1997.
Link | ISI | Google Scholar - 24 Knafelc M, Davenport PW. Relationship between magnitude estimation of resistive loads, inspiratory pressures, and the RREP P1 peak. J Appl Physiol 87: 516–522, 1997.
ISI | Google Scholar - 25 Kellerman BA, Martin AD, Davenport PW. Inspiratory strengthening effect on resistive load detection and magnitude estimation. Med Sci Sports Exerc 32: 1859–1867, 2000.
Crossref | ISI | Google Scholar - 26 Killian KJ, Mahutte CK, Howell JBL, Campbell EJM. Effect of timing, flow, lung volume, and threshold pressures on resistive load detection. J Appl Physiol 49: 958–963, 1997.
Google Scholar - 27 Kisley MA, Noecker TL, Guinther PM. Comparison of sensory gating mismatch negativity and self-reported perceptual phenomena in healthy adults. Psychophysiology 41: 604–612, 2004.
Crossref | ISI | Google Scholar - 28 Logie SL, Colrain IM, Webster KE. Source localization of the early components of the respiratory-related evoked potentials. Brain Topography 11: 153–164, 1998.
Crossref | ISI | Google Scholar - 29 Menon V, Ford JM, Lim KO, Glover GH, Pfefferbaum A. Combined event-related fMRI and EEG evidence for temporal-parietal cortex activation during target detection. Neuro Report 8: 3029–3037, 1997.
Crossref | PubMed | ISI | Google Scholar - 30 Milic-Emili J, Zin WA. Breathing responses to imposed mechanical loads. In: Handbook of Physiology. The Respiratory System. Control of Breathing. Bethesda, MD: Am. Physiol. Soc., 1986, sect. 3, vol. II, pt. 2, chapt. 23, p. 751–769.
Google Scholar - 31 Näätäanen R, Gillard AW, Mantysalo S. Early selective-attention effect on evoked potential reinterpreted. Acta Psychol 42: 313–329, 1978.
Crossref | PubMed | ISI | Google Scholar - 32 Näätäanen R, Picton T. The N1 wave of the human electric and magnetic response to sound: a review and an analysis of the component structure. Psychophysiology 24: 375–425, 1987.
Crossref | PubMed | ISI | Google Scholar - 33 Näätäanen R. The role of attention in auditory information processing as revealed by event-related potentials and other brain measures of cognitive function. Behav Brain Sci 13: 201–288, 1990.
Crossref | ISI | Google Scholar - 34 Puddy A, Giesbrecht G, Sanii R, Younes M. Mechanism of detection of resistive loads in conscious humans. J Appl Physiol 72: 2267–2270, 1992.
Link | ISI | Google Scholar - 35 Revelette WR, Davenport PW. Effects of timing of inspiratory occlusion on cerebral evoked potentials in humans. J Appl Phsyiol 68: 282–288, 1990.
Link | ISI | Google Scholar - 36 Soininen K, Jarvilehto T. Somatosensory evoked potentials associated with tactile stimulation at detection threshold in man. Electroencephalogr Clin Neurophysiol 56: 494–500, 1983.
Crossref | Google Scholar - 37 Strobel RJ, Daubenspeck JA. Early and late respiratory-related cortical potentials evoked by pressure pulse stimuli in humans. J Appl Physiol 74: 1484–1491, 1993.
Link | ISI | Google Scholar - 38 Tapper DP, Duncan SR, Kraft S, Kagawa FT, Marshall S, Theodore J. Detection of inspiratory resistive loads by heart-lung transplant recipients. Am Rev Respir Dis 145: 458–460, 1992.
Crossref | ISI | Google Scholar - 39 Teder W, Alho K, Reinikainen R, Näätäanen R. Interstimulus interval and the selective-attention effect on auditory ERP's: “N1 enhancement” versus processing negativity. Psychophysiology 30: 71–81, 1993.
PubMed | ISI | Google Scholar - 40 Webster KE, Colrain IM. Multi-channel EEG analysis of respiratory evoked potential components during wake and NREM sleep. J Appl Physiol 85: 1727–1735, 1998.
Link | ISI | Google Scholar - 41 Webster KE, Colrain IM. The respiratory-related evoked potentials: Effects of attention and occlusion duration. Psychophysiology 37: 810–818, 2000.
Google Scholar - 42 Webster KE, Colrain IM. The relationship between respiratory-related evoked potentials and the perception of inspiratory resistive loads. Psychophysiology 37: 831–841, 2000.
Crossref | Google Scholar - 43 Webster KE, Adey SA, Clorain IM. The effect of stimulus probability on P3 in the respiratory-related evoked potentials. Psychophysiology 39: 9–15, 2002.
Crossref | ISI | Google Scholar - 44 Webster KE, Colrain IM. P3-specific amplitude reductions to respiratory and auditory stimuli in subjects with asthma. Am J Respir Crit Care Med 166: 47–52, 2002.
Crossref | ISI | Google Scholar - 45 Webster KE, Colrain IM, Davenport PW. P300 from inspiratory occlusion reflects orienting but not startle. Biol Psychol 66: 21–33, 2004.
Crossref | ISI | Google Scholar - 46 Wiley RL, Zechman FW. Perception of added airflow resistance in humans. Respir Physiol 2: 73–87, 1966.
Crossref | Google Scholar - 47 Zechman FW, Hall FG, Hull WE. Effects of graded resistance to tracheal airflow in man. J Appl Physiol 10: 356–362, 1957.
Link | ISI | Google Scholar - 48 Zechman FW, Davenport PW. Temporal differences in the detection of resistive and elastic loads to breathing. Respir Physiol 34: 267–277, 1978.
Crossref | Google Scholar - 49 Zechman FW, Muza SR, Davenport PW, Wiley RL. Relationship of transdiaphragmatic pressure and latencies for detecting inspiratory loads. J Appl Physiol 58: 236–243, 1985.
Link | ISI | Google Scholar - 50 Zhang W, Davenport PW. Activation of neurons in the ventroposteriolateral nucleus of the thalamus by phrenic nerve afferents in cats and rats. J Appl Physiol 94: 220–226, 2003.
Link | ISI | Google Scholar - 51 Zhao W, Martin AD, Davenport PW. Respiratory-related evoked potentials in response to inspiratory occlusions in double-lung transplant recipients. J Appl Physiol 93: 894–902, 2002.
Link | ISI | Google Scholar

